沿轨迹体积获取 HU 值

Get HU values along a trajectory volume

因此,我要做的是计算 CT 中沿轨迹(由目标 x、y、z 及其切线表示)的密度分布 (HU)。目前,我能够沿着穿过目标并与目标(入口)保持一定距离的直线获得剖面图。我想要做的是获取宽度为 1mm 左右的体积(在本例中为圆柱体)的密度分布。 我想我必须沿着体素进行某种插值,因为根据连续坐标之间的间距,几个坐标可以指向相同的索引。例如,这就是我所说的。

此外,我想获得不同形状的轨迹尖端的密度分布,例如:

我的想法是制作一个 3 x 3 矩阵,表示尖端的形状,并将其与体素值进行卷积以获得与尖端对应的 HU 值。如何使用 ITK/VTK?

执行此操作

如果您需要更多信息,请告诉我。 (希望图片足够清晰)

要沿某条路径获取配置文件,您可以使用 resample filter。设置一个变换矩阵,将起点变换为 0,0,0,将终点变换为 x,0,0。将目标图像的大小设置为 x,1,1,间距与源图像相同。

我不明白你的第二个问题。要在尖端获得 HU 值,您可以使用 high quality interpolator (example using linear interpolator) 对该点进行采样。我不明白为什么尖端的形状很重要。

如果你想计算钻尖会遇到的密度,可能最简单的方法是创建一个分辨率高于你的图像的尖端切割表面的掩码。定义一个变换矩阵 M,它将您的钻头放入 CT 图像中的所需位置。

然后遍历蒙版中的所有非零体素,transform indices to physical points, apply transform M to them, sample (evaluate) 使用插值器在该点的 CT 图像中的值,将其乘以蒙版的不透明度(在非二进制的情况下)掩码)并将值添加到 运行 总和。

最后,您的 运行 总和将代表遇到的总密度。该密度总和将取决于尖端切割表面掩模的分辨率。我不知道你如何将它与一些物理量联系起来(比如以牛顿为单位的阻力)。