第7章 图像重建解说.pptVIP

衰减系数的单位H—(豪斯费尔德Hounsfield)， 豪斯费尔德将线性衰减系数分为2000个单位，称为 CT值，一个豪斯费尔德等于水的衰减系数的0.1%，标 度上选择H（水）=0对于空气 H=-1000 衰减最小（不衰减）骨骼 H=+1000 衰减最大X射线经过物体时会发生衰减，不同的物质衰减是不一样的。得到物体的图像最直接的方法是沿Y轴经衰减直接在胶片上成像。这与X光透视是一样的，这样会造成图像的混叠。CT是把物体在Y轴方向划分成小的薄片，薄片的厚度是一个重要的参数，一般为1、2、3、4、5、8、10mm。每个薄片再划分为小的单元，即体素。设某一物体体素对X射线的衰减系数为μ， 体素厚度为d ， 和 为穿透物体前后的X射线的 辐射强度。射线遵循如下的衰减定律：在断层扫描时，生成大量的数据，根据该数据再计算出每个体素的衰减系数，然后把这些衰减系数按一定的函数关系显示在屏幕上，这样，就产生了断层图像。假如一条直线上有n个体素， 第一个体素的衰减为: 第二个体素衰减为： 对于第n个体素有：即：一般情况探测器只能测到 ，而不能测到，因此，不能直接记录各个体素的衰减系数。但是，我们可以用数学方法求解衰减系数。假如某断层有2X2个体素，相应的衰减系数为， 分别从X和Z方向投影，测得的衰减系数为A，B，C，D，即： 从而可以解出的值来。我们用一定的函数关系在屏幕上显示出来就可以得到相应的断层图像。如果图像的分辨率为512X512，则图像有262144个独立阵元，需要解262144元的方程组，计算出μ值，重建出图像。 7.3 投影和雷登变换 下面详细推导根据X射线计算机断层重建图像所需要的数学问题： 笛卡尔坐标系中一条直线可以用斜截式来表示，y=ax+b，也可以用法线来表示-极坐标 平行线束的投影可以由一组直线建模，投影信号中任意一点由沿着直线的射线和给出。工作在连续变量情况下，线求和变成了积分 由冲击函数性质，积分只有沿着直线进行才会不为零 给出沿xy平面中任意一条线的f(x,y)的投影(线积分)的公式就是雷登变换，符号 都可以用来表示f的雷登变化，雷登变换是投影重建的基石。离散情况下：例2：使用雷登变换得到圆形区域的投影 圆的解析式 圆是对称的，只需求0度投影 不管 角是多少，投影都一样，说明 图像关于原点对称。 当雷登变换以为直线坐标显示为一幅图像 时，结果称之为正弦图(sinogram),概念上类似所产生 的傅里叶谱。正弦图包括重建f(x,y)所需要的数据。下 图右图是左图的正弦图，底部一行是矩形在水平方向 上的投影，中间一行是在垂直方向上的投影，水平方 向较窄，物体是对称的，且平行于x,y轴，图像平滑说 明物体具 有均匀灰 度。 如下图是一幅Shepp-Logan幻影图像，是广泛使用的 合成图像，设计用于模拟脑部主要区域的吸收，包括 几个小的肿瘤。解释这个正弦图是很困难的，其视觉 分析仅限于实际应用，但有时对于开发算法是有帮助 的。 CT的关键目的：得到物体的三维表示 1）反投影每一个投影 2）反投影求和产生一个切片图像 3）堆积所有切片图像产生三维物体的再现 公式推导： 单点开始：是全部投影的一个固定 旋转值的投影，反涂抹或者反投影后，线 上每一点的值是，对投影信号中所有重 复这个过程，角度不变，可得到 任意角度均成立，即 对所有反投影图像积分，得到最终图像 离散情况下这种方法形成的反投影图像有时候称为层图（laminogram），是对产生投影的图像的一个近似。从正弦图得到反投影图像的例子，严重的模糊说明直接利用上述推导的公式在实际中是不可接受的，必须在重建中做有意义的改善。 傅里叶切片（投影切片）定理：一个投影的傅里叶 变换是得到投影区域的二维傅里叶变换的一个 切片。如图：任何一投影的一维傅里叶变换可以沿 着一个角度提取一条线的F(u,v)来得到，该角度就 是产生投影所用的角度，得到所有一维投影的切片 组成二维傅里叶变化，反变换求原图f(x,y)。 7.4 傅里叶切片定理 关于 投影的一维傅里叶变换为 关于投影的一维傅里叶变换为表示的二维傅里叶变换。 傅里叶切片（投影切片）定理：一个投影的傅里叶 变换是得到投影区域的二维傅里叶变换的一个 切片。如图：任何一投影的一维傅里叶变换可以沿 着一个角度提取一条线的F(u,v)来得到，该角度就 是产生投影所用的角度，得到所有一维投影的切片 组成二维傅里叶变化，反变换求原图f(x,y)。 7.5 使用平行射线束滤波反投影的重建直接反投影会生成不可接受的模糊结果，如何 消除模糊——反投影之前对投影进行简单的滤波。 F(u,v)的二维反傅里叶变换为 极坐标表达 根据切片定理是一个一维傅里叶反变换的形式，附加 了一个项，这是一个一维滤波函数，