使用 matlab 中的内置函数进行 3D 重建(radon 和 iradon)

3D reconstruction using inbuilt functions in matlab (radon and iradon )

我们能否将氡和氡用于 CT 扫描图像并生成 3D 切片。当然氡,iradon 与 phantom 配合得很好,它会与真正的 CT 投影配合使用吗?

我有从一个物体的不同角度拍摄的 CT 图像

.

我可以在 MATLAB 中使用内置的 radon-iradon 变换吗?我的目标是从 CT 图像重建对象。

取决于几何形状。

如果您的几何形状是平行光束几何形状,那么答案是肯定的,因为在这种情况下您只是一个接一个地进行大量二维扫描。如果将每个探测器投影行视为一个独立的数据,则可以重建许多 2D 切片 (as in the tutorial),然后将它们粘在一起以创建 3D 图像。

但是,如果您处理的是锥形束几何体,则情况有所不同。氡变换不直接适用,图像重建需要更高级的东西(应该知道,我的博士学位就是关于那个!)。如果是这样,我建议您使用 FEX toolbox for CBCT or the well known ASTRA toolbox。第一个没问题,但有一些内在错误(不限制工具箱的功能,你可以从中得到一些不错的图像),第二个被广泛使用。

实际上,最完整的(几何方面)是 ASTRA,因为它具有所有可能的不同几何作为选项,因此您始终可以独立于您的几何使用 ASTRA。

此外,如果您有兴趣获得更好的图像(如果您的投影很少或数据嘈杂,氡变换会很糟糕)您可能有兴趣研究迭代算法,例如 SART、CGLS 或其他算法。如果你在 3D 平行光束中,你可以使用 Per Christian Hansen 的 AIR tools 来做一些改变,或者再次使用 ASTRA 本身。

现在,小心我的小徒弟! 3D 断层扫描重建在计算上 非常昂贵 。非常喜欢。一些 3D 算法可能需要数周的时间 CPU。我希望你有一个糟糕的 NVIDIA GPU 来处理这个问题! (除 AIR 之外的所有链接工具箱都支持 GPU 计算)。

一般旁注:您将很难在 Whosebug 中找到有关此的更多信息。我真的在攻读 CBCT 图像重建的博士学位,这只是你不能在这里问,它的方式比那更复杂。