第六讲-CT成像技术.docVIP

第六讲-CT成像技术（3） X-CT反投影重建成像原理 X-CT图像的质量评价 一、 反投影重建图像的原理— 中心切片定理 CT成像本质上是人体组织的衰减系数μ成像。 成像物理原理为通过CT扫描机构扫描获取求解μ的方程组； 解方程组获得人体某一体层面各个体素的μ值；再将μ值转换为CT值； 最后将CT值变换成能视觉识别的灰度图像。 线性吸收系数 CT成像中物体对X线的吸收起主要作用，忽略对X线的散射作用。 线性吸收系数 在均匀物体中，X线衰减服从指数规律 由于人体器官或组织由多种物质成分和不同密度构成的，X线穿透人体时各点对X线的吸收系数μ不同。 入射第一个体素的X线强度为I0时，透过第一个体素的X线强度Il： 设第二个体素的吸收系数为μ2，X线经第二个体素透射出的强度I2： 第n个体素透射出的X线强度In： CT图像重建 为建立CT图像，必须求出每个体素的吸收系数μ1、μ2、μ3…… μn。 求n个吸收系数，需要建立n个或n个以上的独立方程。 CT装置从不同方向上进行扫描（scanning），来获取足够的数据建立求解μ的方程。 μ与X线能量之间有着依赖关系，即μ随X线能量的增大而减小。 一般X线束以单一频率、固定能量线束穿透物体，可检测到比较稳定的μ。 本质是吸收系数重建。 投照受检体后出射的X线束强度 I 称投影(projection，P) ，投影值的分布称为投影函数。 每个小单元体按照扫描过程中的顺序进行排列和编号，形成一个有序的数组；这些有序的数组在图像平面上形成图像矩阵。 CT图像重建按照这些有序数组计算和重建图像。N×N矩阵中的元素用μij 表示，代表组织的吸收系数或CT值。 头部CT采用256×256或320×320矩阵； 全身CT图像选256×256或512×512矩阵； 显示脊椎骨等结构的细节采用512×512或640×640矩阵。 图像重建的数学基础 1．吸收系数 上式写成对沿着X线束路径随S连续变化的吸收系数连续函数f (s)变化的积分形式： 设断层平面在X-Y中， 断层平面上每一点的吸收系数是坐标(x，y)的函数 f (x，y)。 X线束在平移和旋转扫描中，X线的投影P总是与X线来路径 l 有关，用极坐标(R，θ）来描述X线束路径 l 的位置。 P是随X线束扫描方向和路径的不同而变化，经过坐标变换后，X线束穿过吸收系数 f (x，y)的物体，在R-θ坐标平面上的投影的是函数P(R，θ)。 当在某一θ角度时，将 式表示为平面坐标上的投影P(R，θ)： 分析： 数据采集得到X线束在各个方向上的投影Pθ(R，θ) ； CT图像重建就是要从积分方程式中解出吸收系数 f（x，y）。 根据投影Pθ(R，θ)求解出断面上线性衰减系数f (x，y)分布 ，就是CT图像重建的数学方法问题。 2. δ-函数 又称单位脉冲函数。 性质：若f (t)为连续函数： δ函数推广到二维， δ(x，y)有： 意义：CT图像重建中可用δ-函数来校正反投影法重建中产生的图像模糊。 傅里叶变换重建图像的过程 图像重建算法 图像重建算法是图像矩阵的求解方法。如有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/r，造成图像的伪像。 直接反投影法的图像重建以及修正的过程 直接反投影缺点：影像边缘处不清晰。 如果在一均匀的组织密度内，存在吸收系数极不均匀的部分时，反投影图像与原图像会出现伪影（image artifact）。 反投影