基于投影的图像重建.docVIP

基于投影的图像重建 1 雷登变换和中心投影定理 1.1 雷登变换 以电脑断层扫描技术(CT)为代表的各种层析方法对象直接三维图像重建截面。对于CT线性系数 (1) 其中是路径长度微元，是描述被成像物体的线性系数数学上上的投影变换(见图1)，这一变换表达为2 (2) 其中即，是一个与投影束是一个在投影束称为一个投影；同时为了一般性来表示。狄拉克)“挑选出”与投影平行的直线路径。显然对于平行投影的数据，我们得到 = (3) 记住基本的投影几何中变换由以下两个式子给出： 1或者 2 现在因为旋转是一个保面积坐标系的投影为 (4) 图 1：图示为平行投影射束图 后一个等式可以等效写成为 (5) 1.2 傅里叶切片定理 投影的傅里叶变换为 3或者该换的函数行列式可记作 1 (6) 其中我们可以利用等式(5)得到最后两个表达式。 现在有的二维傅里叶变换，即记为 (7) 变换到傅氏域中的极坐标系 我们有 (8) 比较等式(6)和等式(8)我们有以下结果 (9) 这样，投影的傅里叶变换和原物体在与投影射束正交的射束方向上的频谱是完全相同的。如图2所示，这个重要的结果被称为傅里叶切片定理或者中心切片定理。换句话说，我们已经在与投影成正交的傅氏空间中确定了一组傅里叶系数。 图 2：图示为傅里叶切片定理图 平行射束数据连续滤波反投影的等式 等式(8)的傅里叶逆变换为 (10) 该积分式可被分解为两个积分的和 现有，所以可得到(见等式(8)),所以有 (11) 从傅里叶切片定理(等式(9))可以得到 (反投影) (12) 其中 (滤波) (13) 表示傅里叶变换的逆变换，表示卷积而则是滤波脉冲响应。这样，等式(13)就代表对斜坡滤波器产生的投影起到的滤波作用。这个经过滤波的投影将会对图像中所有沿投影的点起到相同的作用，进而可以被看作为是反投影(见等式(12))。这样的值将由在180°范围内所有这样的投影的和构成(见图3)。因此可将这样一个重建的过程命名为滤波反投影。 3 有限频带的作用 3.1 对一个点的重建 考虑以下的图像 (14) 显然有 (15) 图 3：图示为平行射束反投影的图标(见等式(12)) 从等式(5)可以我们得到投影 (对于所有的而言)