基于投影的图像重建

1 雷登变换和中心投影定理

1.1 雷登变换

以电脑断层扫描技术(CT)为代表的各种层析成像方法对被测对象产生的并不是其直接三维图像，而是从一组投影中重建而成的该物体截面的二维图像。对于CT而言，影像的性质就是线性衰减系数。相应的X射线沿某条路径的变换可由线积分来给出1

(1)

其中是路径长度微元，是描述被成像物体的线性衰减系数二维分布的函数。数学上可将一个投影用雷登变换来描述。雷登变换是一个将一二维函数向极坐标空间上的投影变换(见图1)，这一变换表达为2

(2)

其中即，是一个与投影射束正交的单位向量；是一个在投影射束方向上的向量。称为一个投影；同时为了一般性考虑，影像的性质用二维函数来表示。狄拉克函数()“挑选出”与投影平行的直线路径。显然对于平行投影的数据，我们得到

= (3) 记住基本的投影几何中变换由以下两个式子给出：

1或者

2

现在因为旋转是一个保面积变换3，所以我们在坐标系的投影为 (4)

图 1：图示为平行投影射束图

式可以等效写成为

(5)

1.2 傅里叶切片定理投影的傅里叶变换为

3或者该换的函数行列式可记作 1

(6)

其中我们可以利用等式(5)得到最后两个表达式。现在有的二维傅里叶变换，即记为

(7)

变换到傅氏域中的极坐标系

我们有

(8)

比较等式(6)和等式(8)我们有以下结果 (9)

这样，投影的傅里叶变换和原物体在与投影射束正交的射束方向上的频谱是完全相同的。如图2所示，这个重要的结果被称为傅里叶切片定理或者中心切片定理。换句话说，我们已经在与投影成正交的傅氏空间中确定了一组傅里叶系数。

图 2：图示为傅里叶切片定理图

2. 平行射束数据连续滤波反投影的等式等式(8)的傅里叶逆变换为(10)

该积分式可被分解为两个积分的和现有，所以可得到(见等式(8)),所以有

(11)

从傅里叶切片定理(等式(9))可以得到

(反投影) (12)

其中

(滤波) (13)

表示傅里叶变换的逆变换，表示卷积而则是滤波脉冲响应。这样，等式(13)就代表对斜坡滤波器产生的投影起到的滤波作用。这个经过滤波的投影将会对图像中所有沿投影的点起到相同的作用，进而可以被看作为是反投影(见等式(12))。这样的值将由在180°范围内所有这样的投影的和构成(见图3)。因此可将这样一个重建的过程命名为滤波反投影。

3 有限频带的作用

3.1 对一个点的重建 考虑以下的图像

(14) 显然有

(15)

图 3：图示为平行射束反投影的图标(见等式(12))

从等式(5)可以我们得到投影

(对于所有的而言) (16)这样对投影的傅里叶变换为

(对于所有的而言) (17)利用等式(10)我们可以得到 (18)

其中是对的重建。考虑到重建的有限带宽的可积情况，我们得到 (19)

其中是空域中的波数，否则系统是不敏感的。所有在这个平面上的任意角度的像点都可以用极坐标来表示为

这样重建图像可由下式给出

令，可以得到

(20)

等式(20)中的第二重积分是在一个完整的圆上，这样就有 (21)

更进一步来说，因为是偶函数，则上式可以写成

(22)

幸运的是，对积分的求解有一个标准的表解法(见Abromovitz andStegun中9.1.21，第360页)，由下式给出 (23)

其中是零阶第一类贝塞尔函数(例见Kreyszig中第208页)，这样得到 (24)

利用替换式上式可以由下式给出 (25)

幸运的是，这个积分的求解同样有一个标准的表解法(A&S中11.3.20，第484页)，即

(26)

其中是阶第一类贝塞尔函数。这样，对于，条件而言，我们可以得到

(27)

图4显示了对于不同的值的结果。可以注意到一个很有趣的结果，这与由圆孔夫琅禾费衍射得到的结果是一个形式。这样，有限带宽的作用即是对图像像素的模糊作用。

可以利用对等式(27)在上的积分来验证这个结论是正确的。于是就有

其中我们已经用过了标准结果，即

图 4：有限带宽对平面中心的图像象点重建的作用(见等式(27))。(a),(b), (c) (如果距离单位是，

则的单位是)

3.2 一般情况

我们考取更加一般的情况。令为图形性质的实际分布。例如在X射线CT的 其中是线性吸收系数的二维分布函数。任何图像都可以描述成如下式

(28)

利用这个重建的图像，(对于连续的情况)将和实际的图像通过下式相联系起来

从前面部分我们已经得到 (29)于是，我们得到

(30)