使用sitk.ResampleImageFilter()与复合变换重采样MR图像时出现灰色伪影的问题求助
使用sitk.ResampleImageFilter()与复合变换重采样MR图像时出现灰色伪影的问题求助
各位大佬好,我最近在处理一组交错采集的T2加权MR图像,为了消除旋转、平移这类几何畸变,分别用FSL和ANTs做了序列切片间的配准。之后我想用SimpleITK的CompositeTransform构建一个复合变换,来完成整体的图像校正。
我仔细检查过生成的复合变换,逻辑上应该是没问题的,但不管用FSL还是ANTs的配准结果,重采样后的图像在最后几个切片(或者切片的部分区域)都会出现灰色伪影。而且这个问题和我往复合变换里加的变换数量无关——不管只加1个,还是加40个,甚至我按「单位矩阵→第1、2个切片的变换矩阵→依次加入所有切片的变换矩阵」的顺序组合,伪影依然存在。
我试过用SimpleITK的SetUseNearestNeighborExtrapolator()方法,这确实能消除伪影,但带来了新问题:当移动图像的像素超出网格范围时,只会把参考图像的像素映射到配准后的图像上,这不符合我的需求。
我已经尝试了各种复合变换的组合方式,实在想不到其他办法了,有没有大佬能给我一些思路或者解决方案?
备注:内容来源于stack exchange,提问作者marlene schichl




