图像重建资料.ppt

X 射线经过物体时会发生衰减，不同的物质 衰减是不一样的。得到物体的图像最直接的方法 是沿 Y 轴经衰减直接在胶片上成像。这与 X 光透视 是一样的，这样会造成图像的混叠。 CT 是把物体在 Y 轴方向划分成小的薄 片，薄片的厚度是一 个重要的参数，一般 为 1 、 2 、 3 、 4 、 5 、 8 、 10mm 。每个薄片再 划分为小的单元，即 体素。 设某一物体体素对 X 射线的衰减系数为 μ ， 体素厚度为 d ， 和 为穿透物体前后的 X 射线的 辐射强度。射线遵循如下的衰减定律： 在断层扫描时，生成大量的数据，根据该数 据再计算出每个体素的衰减系数，然后把这些衰 减系数按一定的函数关系显示在屏幕上，这样， 就产生了断层图像。 d e ? ? ? ? ? 0 ? ? ? 0 1 ln d ? 0 ? ? d e 1 0 1 ? ? ? ? ? 假如一条直线上有 n 个体素， 第一个体素的衰减为 : 第二个体素衰减为： d ) ( d d d e e e e 2 1 2 1 2 0 0 1 2 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 对于第 n 个体素有： d ) ( n n e ? ? ? ? ? ? ? ? ? ? ? 2 1 0 n n ln d ? ? ? ? ? 0 2 1 1 ? ? ? ? ? 即： 一 般 情 况 探 测 器 只 能 测 到 ， 而 不 能 测 到 ，因此，不能直接记 录各个体素的衰减系数。但是，我们可以用数 学方法求解衰减系数。 0 1 1 ln n i i n d ? ? ? ? ? ? n ? 1 2 1 , , n ? ? ? ? n n ln d ? ? ? ? ? 0 2 1 1 ? ? ? ? ? 假如某断层有 2X2 个体素，相应的衰减系数 为 ， 22 21 12 11 , , , ? ? ? ? 分别从 X 和 Z 方向投影， 测得的衰减系数为 A ， B ， C ， D ，即： D C B A ? ? ? ? ? ? ? ? 22 12 21 11 22 21 12 11 ? ? ? ? ? ? ? ? 从而可以解出 的值来。我们 用一定的函数关系在屏幕上显示出来就可以得 到相应的断层图像。 如果图像的分辨率为 512X512 ，则图像有 262144 个独立阵元，需要解 262144 元的方程 组，计算出 μ 值，重建出图像。 22 21 12 11 , , , ? ? ? ? 7.3 投影和雷登变换 下面详细推导根据 X 射线计算机断层重建图像所需要 的数学问题： ? 笛卡尔坐标系中一条直线可以用斜截式来表示， y=ax+b ，也可以用法线来表示 - 极坐标 cos sin x y ? ? ? ? ? ? 平行线束的投影可以 由一组直线建模，投 影信号中任意一点由 沿着直线 的 射线和 给出。工作在 连续变量情况下，线 求和变成了积分 cos sin k k j x y ? ? ? ? ? ( , ) ( cos sin ) ( , ) k k j j k f x y x y dxd g y ? ? ? ? ? ? ? ? ?? ?? ? ? ? ? ? 由冲击函数性质，积分只有沿着直线进行才会不为零 ? 给出沿 xy 平面中任意一条线的 f(x,y) 的投影 ( 线积分 ) 的公式就是雷登变换 ， 符号 都可以用来表示 f 的雷登变化，雷登变换是投影重建 的基石。离散情况下： { ( , )}, { }, ( , ) R f x y R f g ? ? 1 1 0 0 ( , ) ( , ) ( cos sin ) M N x y g f x y x y ? ? ? ? ? ? ? ? ? ? ? ? ? ?? , , , ( , ) x y f x y ? ? ? ? ? ? 是离散变量，如果固定 令 变化，投影简化为 沿着这两个参数的特定值定义的线对 的像素求和。 通过增加覆盖图像的所要求的所有 的值产生一个投影， 改变 并重复这个过程产生另一个投影，类推 ... 例 2 ：使用雷登变换得到圆形区域的投影 圆的解析式 2 2