平行束反投影重建

一、 平行束反投影重建算法

平行束 重建采用的是平移加

旋转的扫描方式，如图所示，射线源在某一角度下水平移动，将物体全部

照射后旋转一角度，如此重复，在这个过程中探测器相应地运动以接收X 射线。

1、反投影重建算法的物理概念： 断层平面中某一点的密度值可以看作是这一平面内所有经过该点的射线的投影值之和（的均值）。

整幅重建图像可以看作是所有方向下的投影累加而成。

射线标号示于图中，像素值（代表密度）分别1x ,2x ,3x ,4x ，

赋值如下:

15x =,20x =,32x =,418x =

根据投影的定义(某条射线投影值为该条射线穿过的所有的像素值之和)，每条射线的投影i p (1,2

i =)为：

1215p x x =+=, 23420p x x =+=,3137p x x =+= 42418p x x =+=, 532p x ==, 61423p x x =+=

720p x ==

根据反投影重建算法的物理意义，重建图像中各像素，得到：

113635x p p p =++=，214723x p p p =++=， 323529x p p p =++=，424661x p p p =++=

(a) 原图像像素值 (b)反投影重建后图像 (c)求平均后图像

图 反投影示例

重建后的图像如图(b)所示，可以看出原图像中像素值不为零的点反投影重建后仍较突出，但原图中像素值为零的点，经反投影重建后不再为零，即有伪迹。有时为了使重建后图像的像素值更接近于原图的像素值，在求反投影时，把数据除以投影的数目（即射线数），

图 平行束平移加旋转

图 断层像素值和射线

如图(c)所示。 因此有：

,1

1p

n k k i

i p

x p

n ==

∑

该式可作为反投影重建算法的计算式。其中k x 表示像素k 的值，

i k p ,表示经过像素k

的第i 条射线投影，p n 表示图像内的射线条数。

图(a)表示空间中一个孤立点源A ，密度为1。经过A 点的三条射线也示于图中。射线束理论上可以很多，取三条示意。不经过A 点的射线投影为零，经过A 点的射线投影值均为1，

1231p p p ===。

(a) 孤立点源A 及三条射线 (b)相应的反投影重建图，有星状伪迹

图 孤立点源的反投影重建及星状伪迹

经反投影重建后，得到A 点的像素值为123()/31A f p p p =++=。A 点以外的像素值原来为零，经反投影重建后不等于零，而是等于1/3。所以，经反投影重建后的图像除保留A 点的像外，还有像素值为1/n 的灰雾背景，后者称为星状伪迹。产生星状伪迹的原因在于：反投影的本质是把取自有限物体空间的投影均匀地回抹（反投影）到射线所及的无限空间的各点上，包括原先像素值为零的点。

(1)(2)(3)

二、 反投影重建算法的数学描述

我们把“取投影” 、“反投影重建” 、“重建后图像”这些环节看作是一个以原像为输入，重建后图像为输出的成像系统，如图所示，先来求该系统的点扩展函数PSF 。

图中，(,)x y 为固定坐标，(,)r r x y 为旋转坐标，(,)r θ为极坐标，设位于坐标原点0,0x y ==的点源(,)x y δ为

x y -断面中唯一的像点，扫描方式为平移加旋转。即射

线先平行移动，从物体的一侧移向另一侧；然后旋转一个角度φ∆，如此继续，直到累计的旋转角达到180φ

-∆。为止。为了计及这一几何要求我们设置一旋转坐标系统

r r x y -，它绕原点转动使投影总是沿着r y 方向，r r

x y -的原点与x y -的原点相同。二者的夹角为φ，不同的φ代表不同的投射方向。投影线的位置可由(,)r x φ完全确定。

设φ为离散取值，如i φφ=，则相应的投影为：

()(,)(,)(,)()()()|(cos())

i i r r i r r r r r r r

r r r r i p x p x f x y dy x y dy x y dy x r φφφφδδδδδθφ+∞+∞

-∞

-∞

+∞=-∞

======-⎰

⎰⎰

其中cos()r x r θφ=-，这可以从图的几何关系很容易得出。 若n φφ=，则相应的投影为：

cos n n p r φδθφ=-[()]

根据反投影重建的定义式,点(,)r r x y 的图像在所述坐标系统中表示为：

[]

[]cos()

1

1

1

1

1

(,)(,)()|cos()1

cos()i r i i

i

N N r r r x r i

i i N i

i f r f x y p x p x r N N p r φ

φ

φ

φθφφφ

φφθθφθφφ

π

=-=====

=-=-∆∑∑∑

其中，N φ为投影数，/N φφπ∆=。若在有限区间内将射线增至不相重的无限条，即连续取投影，则有：

[]0

1

(,)cos()f r p r d π

φθθφφπ

=

-⎰

在忽略射束硬化的情况下，φ在(,2)ππ区间内的投影值等于(0,)π区间内的投影值。 在输入图像为点源的情况下，由及可得：

(,x y δ(,)

h x y 重建后图像

原像图 “反投影重建”成像系统

图平移加旋转扫描方式所用坐标系