SimpleITK：应用逆变换

Question

我最近开始使用SimpleITK进行图像注册。对于我当前的项目，我需要注册X射线图像和CT图像，然后将逆矩阵应用于在X射线图像上绘制的ROI蒙版。

我得到这行的逆矩阵：

inverse_transform = final_transform.GetInverse()

如何将转换应用于ROI蒙版？

Answer 1

欢迎使用SimpleITK！

假设您的X射线是配准中的固定图像（CT是移动中的），那么配准的结果就是从X射线到CT的转换映射点。您需要做的就是使用逆变换将ROI蒙版图像重新采样到CT上。

transformed_labels = sitk.Resample(xray_roi_mask, ct_image, inverse_transform, sitk.sitkNearestNeighbor, 0.0, #out of bounds pixel color xray_roi_mask.GetPixelID())

this Jupyter notebook的最后一个单元格将执行此操作。

另外两条评论：

1. 不确定X射线进行CT配准的含义。这是2D / 3D注册，目前SimpleITK不支持该注册，因此我假设您正在执行2D / 2D或3D / 3D注册。
2. 由于您是SimpleITK的新手，建议您查看我们的main Jupyter notebook repository或更简洁的IEEE ISBI'18 tutorial，以熟悉该工具包。