Siddon算法参数说明及通俗解释

news2025/2/27 1:32:01

参数	含义
X1,Y1,Z1	射线的起点坐标。
X2,Y2,Z2	射线的终点坐标。
X(a), Y(a), Z(a)	射线的参数化表示，其中a是参数，范围从0到1。
Xpm[i]	第i个x方向的平面位置。
Ypm[j]	第j个y方向的平面位置。
Zpm[k]	第k个z方向的平面位置。
dx,dy,dz	x、y、z方向平面之间的间距。
ax[i], ay[j], az[k]	射线与第i个x平面、第j个y平面、第k个z平面的交点参数值。
αmin, αmax	射线进入和离开CT数组的参数范围。
i,j,k	体素的索引，表示体素在x、y、z方向的位置。
d12	射线的总长度。

定义射线
- 将射线从起点(X1,Y1,Z1)到终点(X2,Y2,Z2)参数化为X(a)=X1+a(X2−X1)，Y(a)=Y1+a(Y2−Y1)，Z(a)=Z1+a(Z2−Z1)，其中a是参数，范围从0到1。就像将一段线分成许多小段，每段用a来表示位置。
计算射线与CT数组边界的交点
- 确定射线进入和离开CT数组的参数值αmin和αmax。这就像找到射线进入和离开网格墙的位置。
计算射线与正交平面的交点
- 对于每个正交平面集（x、y、z方向），计算射线与这些平面的交点参数值。这些平面是等间距的，所以可以通过递归生成交点参数值。就像在网格墙上画了许多等间距的线，找出射线与这些线的交点。
合并参数值集
- 将三个平面集（x、y、z方向）的交点参数值合并成一个集合，并计算相邻参数值之间的差异，以确定射线在每个体素中的长度。这就像把所有交点按顺序排列，计算每两个相邻交点之间的距离，这段距离就是射线在对应体素中的长度。
计算体素长度和密度的乘积
- 对于每个体素的交点长度，获取对应的体素索引，并将长度与密度的乘积累加，以得到总的射线路径。这就像把每个小段的长度乘以该位置的密度，然后把所有小段的值加起来，得到射线路径的总效果。