2

我正在尝试注册来自大脑的两个体积图像（PET 和 CT 甚至 PET 和 MR）。这些体积图像中的每一个都包含不同数量的 2D 图像（切片）。例如，CT有150切片和PET有100切片。我正在考虑使用插值方法来计算 CT 切片的数量并将其减少到 100。这是正确的方法吗？有谁知道任何可能对我有帮助的资源？就像一个伪代码，或者我应该通过注册两个体积图像的步骤。谢谢 ：）

4

1 回答 1

2

如果您知道 150 个 CT 切片和 100 个 PET 切片的间距信息，您可以查看 MATLAB 的interp1沿一个轴插值的函数，以将图像重新缩放到相同数量的像素。从这里可能可以使用 MATLABimregister来执行注册。

如果您想了解注册是如何在后台工作的（在像素和物理坐标之间转换、转换/重新采样图像等），我可以指导您访问的一个资源是ITK 软件指南 pdf。

特别是，尝试阅读关于图像表示的第 1 册第 4.1.4 节（pdf 第 41 页），以及关于变换的第 2 册第 3.9 节（pdf 第 532 页）。

一般来说，在配准中使用 3D 图像进行转换和插值的问题编写代码可能非常麻烦。您需要询问自己像素的间距和方向，如何转换和插入图像以使其网格重叠，并且在评估相似性度量时，您还需要决定如何处理网格中位于图像边界之外的像素.

虽然由您决定做您认为最好的事情，但我建议您使用现有的注册程序，如果它们能够做您想做的事情：

  • MATLAB imregister（我从未使用过，所以无法评论）
  • 用于 Python 的simpleITK
  • C++ 的 ITK有一个学习曲线，但可以完全控制注册过程
  • elastix是一个命令行程序，它使用参数的文本文件来执行注册。
  • 3D 切片器具有用于简单线性配准的图形用户界面
于 2015-08-11T19:47:22.047 回答