第四章 CT成像系统C从投影重建图像的原理中心切片定理.ppt

卷积反投影法的过程 4.5 从投影重建图像的算法— 扇形束反投影重建算法 X射线以等角度分布形成的扇形束 检测器以等间距分布形成的扇形束 扇束投影重建算法的分类 重排算法 把一个视图中采得的扇形数据重新组合成平行的射线投影数据，然后用上一节的算法重建。 直接重建算法 不必数据重排，只需适当加权即可运用与上一节类似的算法重建 * * * * * * * * * * * * * * * * * * * * * 1 1 rect(ρ/2)函数的说明 Circ(ρ)函数的形状 4.4 从投影重建图像的算法 图像重建算法 图像重建算法是图像矩阵的求解方法。如有N×N的图像矩阵，有N×N个独立的线性方程组，并且求解N×N个矩阵中的体素的吸收系数μij。 N×N个方程组求解可以采用迭代法（逐次近似法）等。 现在应用比较多的是投影法、二维傅里叶重建法和滤波反投影法（filtering back projection，FBP）。 反投影法（back projection） 又称总和法或线性叠加法。它是利用所有X线的P值计算各个像素的μ值的二维分布。 基本原理：是将所测得的投影值按其原路径平均分配到每一点上，各个方向上投影值反投影后，在影像处进行叠加，推断出原图像。 反投影法 例：对四个体素( μ1=1， μ2=2， μ3=3， μ4=4 )矩阵的图像重建： 对四体素矩作0°、45°、90°、135°投影（扫描），再将投影值反投回原矩阵的对应位置上，即可将原矩阵中的四体素的特征参数值解出。 一、直接反投影 1.直接反投影原理 f(x,y) 投影线L X Y θ R 2. 数学上的定量分析与修正 以δ函数为密度的函数，直接反投影创建的图像，带有“长尾巴”的1/r，造成图像的伪像。 直接反投影法的图像重建以及修正的过程 缺点：影像边缘处不清晰。 如果在一均匀的组织密度内，存在吸收系数极不均匀的部分时，反投影图像与原图像会出现伪影（image artifact）。 反投影数量愈多，重建图像愈接近于原图像，但由于存在星形伪影，而使得重建图像的边缘部分模糊不清。 算法举例 根据反投影算法x1=p5 = 5 x6=p2+p3+p5=18 … 平均化处理，除以投影线数目 xi=xi/6 0 0 0 0 0 5 2 0 0 1 0 0 0 0 0 0 5 6 2 3 7 18 12 7 1 10 8 1 3 6 2 5 0.83 1 0.33 00.5 1.16 3 2 1.16 0.06 1.66 1.33 0.16 0.5 1 0.33 0.83 反投影重建后 原像素值 再除以投影线数，平均化 断层平面中某一点的密度值可看作这一平面内所有经过该点的射线投影之和的平均值 1 2 3 4 5 6 星状伪迹 反投影重建后，原来为0的点不再为0，形成伪迹 0 0 0 0 0 5 2 0 0 1 0 0 0 0 0 0 0.83 1 0.33 0.5 1.16 3 2 1.16 0.06 1.66 1.33 0.16 0.5 1 0.33 0.83 原像素值 再除以投影线数，平均化 星状伪迹 产生星状伪迹的原因在于：反投影重建的本质是把取自有限物体空间的射线投影均匀地回抹(反投影)到射线所及的无限空间的各点之上，包括原先像素值为零的点 重建图像的边缘模糊原因： 重建的 fb(x，y)与实际的 f (x，y)不一样。 为获得真实的吸收系数 ，必须对 fb(x ， y)进行修正。 反投影吸收系数 fb(x，y)与实际 f(x，y)之间存在一个1/r，1/r称为模糊因子。 滤波反投影算法 去除星状伪迹的两种方法： 取投影 取投影 一维滤波器 反投影重建 二维滤波器 反投影重建 输入图像（原图像） 输入图像（原图像） 输出图像(同原图像) 输出图像(同原图像) 星状尾迹的去除 （b） （a） 第一种方法中的二维滤波器较难实现。 第二种方法交换了滤波与反投影的顺序，此时投影数据是一维的，易于实现。此方法的理论基础就是中心切片定理。 滤波反投影算法 滤波反投影法 滤波反投影法为了消除模糊因子1/r的影响，并将二维傅里叶变换改为只进行一维傅里叶变换，既可校正失真，又简化计算，提高图像重建速度。 滤波反投影法也称卷积反投影法（convolution back projection，CBP）。 滤波反投影法是在反投影之前，对所有的投影数据进行滤波或卷积，使结果的图像无“星月状”伪影。 重建成像过程： （1）获取全部的投影数