微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

如何维护 SimpleITK 图像中的方向到 numpy 数组转换?

如何解决如何维护 SimpleITK 图像中的方向到 numpy 数组转换?

我有一个物体的三个不同的各向同性 MRI DICOM 体积,每个体积都有不同的方向(同一物体的正交矢状、冠状和横向采集)。

我想将它们转换为 numpy 数组并以它们的索引匹配的方式绘制它们。假设我有从 sik 图像发出的三个 numpy 数组:

sag_array = sitk.GetArrayFromImage( sag_sitk )
dors_array = sitk.GetArrayFromImage( dors_sitk )
trans_array = sitk.GetArrayFromImage( trans_sitk )

我希望能够使用相同的索引来绘制它们,以便切片

sag_array[:,:,index]
dors_array[:,index]
trans_array[:,index]

对应相同的视图,没有轴的翻转或反转。

我猜这个信息包含在 SimpleITK 图像的方向中,有没有办法在转换后将其传输到 numpy 数组?

Direction 属性一般对 numpy 转换有什么影响,还是丢失了?

解决方法

我通过使用 sitk.Resample() 函数将所有图像预处理为共同的原点和方向来解决它。这样,当转换为 numpy 数组时,由于它们占用相同的物理空间,因此它们之间被连贯地切片。

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。