平行束反投影重建.docVIP

平行束反投影重建 平行束反投影重建 PAGE 平行束反投影重建 平行束反投影重建算法 X射线管平移平行束 重建采用的是平移加旋转的扫描方式，如图所示，射线源在某一角度下水平移动，将物体全部照射后旋转一角度，如此重复，在这个过程中探测器相应地运动以接收X射线。 X射线管 平移 1、反投影重建算法的物理概念： 断层平面中某一点的密度值可以看作是这一平面内所有经过该点的射线的投影值之和（的均值）。 探测器整幅重建图像可以看作是所有方向下的投影累加而成。 探测器 平移射线标号示于图中，像素值（代表密度）分别,,,， 平移 图 平行束平移加旋转赋值如下: 图 平行束平移加旋转 ,,, 根据投影的定义(某条射线投影值为该条射线穿过的所有的像素值之和)，每条射线的投影()为： , , , , 根据反投影重建算法的物理意义，重建图像中各像素，得到： 图 断层像素值和射线，， 图 断层像素值和射线 ， (a) 原图像像素值 (b)反投影重建后图像 (c)求平均后图像 图 反投影示例 重建后的图像如图(b)所示，可以看出原图像中像素值不为零的点反投影重建后仍较突出，但原图中像素值为零的点，经反投影重建后不再为零，即有伪迹。有时为了使重建后图像的像素值更接近于原图的像素值，在求反投影时，把数据除以投影的数目（即射线数），如图(c)所示。 因此有： 该式可作为反投影重建算法的计算式。其中表示像素的值，表示经过像素的第条射线投影，表示图像内的射线条数。 图(a)表示空间中一个孤立点源A，密度为1。经过A点的三条射线也示于图中。射线束理论上可以很多，取三条示意。不经过A点的射线投影为零，经过A点的射线投影值均为1，。 (a) 孤立点源A及三条射线 (b)相应的反投影重建图，有星状伪迹 图 孤立点源的反投影重建及星状伪迹 经反投影重建后，得到A点的像素值为。A点以外的像素值原来为零，经反投影重建后不等于零，而是等于1/3。所以，经反投影重建后的图像除保留A点的像外，还有像素值为1/n的灰雾背景，后者称为星状伪迹。产生星状伪迹的原因在于：反投影的本质是把取自有限物体空间的投影均匀地回抹（反投影）到射线所及的无限空间的各点上，包括原先像素值为零的点。 反投影重建算法的数学描述 原像重建后图像取投影反投影重建我们把“取投影” 、“反投影重建” 、“重建后图像”这些环节看作是一个以原像为输入，重建后图像为输出的成像系统，如图所示，先来求该系统的点扩展函数PSF。 原像 重建后图像 取投影 反投影重建 图 “ 图 “反投影重建”成像系统 图平移加旋转扫描方式所用坐标系图中，为固定坐标，为旋转坐标，为极坐标，设位于坐标原点的点源为断面中唯一的像点，扫描方式为平移加旋转。即射线先平行移动，从物体的一侧移向另一侧；然后旋转一个角度，如此继续，直到累计的旋转角达到为止。为了计及这一几何要求我们设置一旋转坐标系统，它绕原点转动使投影总是沿着方向，的原点与的原点相同。二者的夹角为，不同的代表不同的投射方向。投影线的位置可由完全确定。 图平移加旋转扫描方式所用坐标系 设为离散取值，如，则相应的投影为： 其中，这可以从图的几何关系很容易得出。 若，则相应的投影为： 根据反投影重建的定义式,点的图像在所述坐标系统中表示为： 其中，为投影数，。若在有限区间内将射线增至不相重的无限条，即连续取投影，则有： 在忽略射束硬化的情况下，在区间内的投影值等于区间内的投影值。 在输入图像为点源的情况下，由及可得： 因为，(此公式推导可参见附录一)，其中，是的唯一的根。 令，则有： 因为，所以，于是有： 可见，相应于反投影重建算法的系统，它的点扩展函数不是函数，系统不是完美的。虽然在处能反映原图是点源的情况，但在处，像素值，上式定量地描述了反投影重建算法星状伪迹的本质。 若原像为，则将原像取投影后再按反投影算法得到的重建图像为：，其中表示二维卷积。 要去除反投影算法的星状伪迹，可以在输出端加一滤波器，使加了滤波器后的反投影重建等效成像系统的点扩展函数(PSF)为。设滤波器的PSF为，相应的传递函数为，要求(式的推导可参见附录一)：