第n篇 图像133.pptVIP

2) 、为了重建这一物体，需要无穷多个或连续的投影数据；因为只有这样才能完全确定傅立哀空间的物体分布，从而通过反变换确定实际物体。实际上用有限个投影角，每个投影以有限个样本就可以得到满意的重建。 6.5 反投影与累加图像 反投影是指投影之逆，形象的解释是取 的一维函数(投影数据)，把它沿 方向向整个空间均匀地抹一次，由此产生一个二维分布 。 用数学表示，一维投影 的反投影由下式给出： 极坐标下： 把许多反投影加以组合就会产生一个新的二维图像，称为累加图像。如果我们取相应于不同角度 的一组物体投影的一系列离散的反投影，则累加图像就是各个反投影的算数和。 离散情况下： 在连续情况下： 对θ的积分区间在0～π之间就可获得全部信息。b（r，α）就是反投影图像，也有称其为“层图” （Layergram）。 反投影的运算概念可由下图说明： 把一个一维的投影数据均匀地回抹到二维空间。该图表示的是四个投影 及其累加图，后者是四个相应的反投影相加形成的，原物体是具有不同尺寸的两个圆盘。 傅立哀变换重建 傅立哀变换重建是较简单的重建方法，它是1974年由shepp和logan开发的算法。 其基本原理如下： 令f（x，y）是一幅图像，其二维傅立哀变换如下式所示， 而f(x,y)在x轴上的投影可表示为下式： 这一投影的一维傅立哀变换为： 这一表达式恰好与f(x,y)的二维傅立哀变换 的表达式一致，即： 现在假设将函数投影到一条经过旋转了 角的直线上。 如果投影选作 轴，投影点通过距 轴距离为 处的平行线对函数积分有： 它的一维傅立哀变换为： 也可写成下式： 作一个代换：(为使展开式与2D傅立哀变换一致） 如果（u，v）点处在 角一定，而距原点的距离r可变的这样一条直线上，则 如果投影变换对所有的 都是已知的，图像的二维变换也是可以确定的，为了得到原图像，只需作一个反傅立哀变换就可以了 这就是用傅立哀变换的基本重建技术。 Radon反变换 由中心切片定理 写成极坐标形式： 这里写成两部分是必要的，因为傅立叶空间的极坐标定义区间为： 而数据采集坐标θ和 的定义在： 用极坐标表示傅立叶反变换为： 为了把物体函数用它的投影 清楚的表示出来，我们把上式分成两部分： 利用投影定理 此式利用了 则： 合并二式： 这一式子可理解为 的一维付里哀反变换。该式也可以写成下列形式：  式中： 应用卷积定理： 式中 代表柯西主值。 导出这一结果，利用了下边的关系式： 这个结果中利用了付里哀变换的对偶性质 由卷积的定义,容易证明: 因此,上式为： 这就是Radon反变换的一种形式。 物体空间、付里哀空间和Radon空间的关系如下图所示： 更详细的内容可参考： Harrison H.Barrett William Swindell “Radiological Imaging” “The Theory of Image Formation Detection and Processing “ 其中： 3) 极坐标 令f(x,y)是介质的衰减系数，如果衰减系数在介质中为恒定的，即 f(x,y)=μ。 如果x射线行进距离为x，则射线衰减规律将遵循比尔定律（Beer’s law)，即： 式中 ： 为入射光子积分通量 为透射的光子积分通量。 如果衰减系数不是恒定值则要沿吸收路径作积分 L为x射线源到检测器的直线路径。 为得到一个透射量的数组，可作如下处理： 即： 作线性化处理： 这里 就叫雷登变换（Radon Transformation） 或称为f(x,y)的投影 这里物体断层重建的任务就是求 的逆。 这里，我们可把 看成是在给定的投影角