沿轨迹体积获取HU值

Question

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

此外，我想获得轨迹尖端不同形状的密度曲线，例如：

我的想法是我制作一个3乘3的矩阵，表示尖端的形状，并将其与体素值进行卷积，得到与尖端相对应的HU值。我怎么能用ITK / VTK做到这一点？

如果您需要更多信息，请告诉我。 （我希望图像足够清晰）。

Answer 1

要在某个路径上获取个人资料，您可以使用resample filter。设置变换矩阵，将起点转换为0,0,0，将终点转换为x，0,0。将目标图像的大小设置为x，1,1，并与源图像中的间距相同。

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

Answer 2

如果要计算钻头尖端会遇到的密度，最好以高于图像的分辨率创建尖端切割面的蒙版。定义一个变换矩阵M，将您的钻头放入CT图像中的所需位置。

然后遍历掩码中的所有非零体素，transform索引到物理点，apply将M变换为它们，采样（evaluate）CT图像中的值在那一点上使用插值器，将它乘以掩码的不透明度（在非二进制掩码的情况下）并将值加到运行总和上。

最后，您的运行总和将代表遇到的总密度。此密度总和将取决于尖端切割表面的掩模分辨率。我不知道你将如何将它与一些物理量相关联（比如用牛顿抵抗力）。