从3D体积数据中提取像素的不同视图

Question

我有3D CT数据。我想提取三个不同的贴片，以对应于三个3个不同视图（轴向，冠状，矢状）的特定像素为中心。谁能告诉我如何在python中做到这一点？

Answer 1

如果您的音量为I且您想要提取补丁的点是x0，yo，zo且补丁大小为p而不是

I[xo,y0-p/2:y0+p/2,z0-p/2:z0+p/2] = Pach1
I[xo-p/2:x0+p/2,y0-p/2:y0+p/2,z0] = Pach2
I[xo-p/2:x0+p/2,y0,z0-p/2:z0+p/2] = Pach3