图象重建讲义.ppt

3.1 卷积逆投影重建 3、扇束投影重建 (a) 发射器和接受器均旋转 (b) 发射器旋转，接受器固定 3.1 卷积逆投影重建 4、傅里叶反变换和卷积逆投影重建法的比较 傅里叶反变换重建 ：直角坐标 卷积逆投影重建 ：极 坐标 主要优点： 卷积逆投影重建基本算法： 容易用软件和硬件实现 可重建出准确清晰的图象 傅里叶反变换重建： 计算量比较小 3.2 其他逆投影重建 1、逆投影滤波 逆投影滤波器（filter of the back-projections） 先逆投影，再滤波（/卷积） 3.2 其他逆投影重建 2、滤波逆投影 （backprojection of the filtered projections） 先对投影滤波，再逆投影 4 级数展开重建 1、离散模型 代数重建法（ART），迭代算法，优化技术 第i条射线与 第j个象素相交 的长度aij，代 表第j个象素沿 第i条射线的贡 献的权值 4 级数展开重建 2、无松弛的代数重建技术 先初始化1个图象矢量x(0)，然后用下式进行迭代： 思路： 每次取1条射线，改变图象中与该线相交的象素的值，从而把当前的图象矢量x(k)更新为x(k+1)。具体运算中就是将测量值与由当前算得的投影数据 的差正比于aij重新分配到各个象素上去 4 级数展开重建 3、松弛的代数重建技术 将无松弛的代数重建技术推广： r是松弛系数，它控制了收敛的速度 重建结果示例 {图4.3} 重建误差： 图象的4个角处（较少的射线） 内部正方形的4个角处（吸收值急剧变化） 4 级数展开重建 4、 级数法的一些特点 (1) 由于在空域中比较灵活，所以容易调整以适应新的应用 (2) 能重建出较高对比度的图象（对密度突变的材料） (3) 借助多次迭代可从较少投射（ 10）重建图象 (4) 比变换法更适合于ECT (5) 比变换法更适合于3-D重建问题 (6) 比变换法更适合于不完整投影情况 5 综合重建方法 迭代变换法（iterative transform） L(s, ?n)是直线(s, ?n)与Q相交段的长度 重建工作： 给定函数 f (x, y) 的投影g(s, ?n)，其 中s可以取所有实 数，?n是1组N个离 散角度，要重新得到 f (x, y) 5 综合重建方法 迭代变换法（iterative transform） 迭代公式 n = (i mod N ) +1，s = xcos?n + ysin?n g(i)(s, ?n )是 f (i)(x, y)沿?n的投影 图象序列 将收敛 5 综合重建方法 迭代变换法（iterative transform） “半离散”函数g( ·, ?n)的离散形式 引进插值函数q(·)从采样 g (m?s, ?n)估计 g (s, ?n) 引进基函数B(x, y)从采样 f (k?x, l?y)估计 f (x, y) 图象重建 投影重建 从一个物体的投影图重建目标图象的过程 输入投影图 ?（图象处理）? 输出重建图 一类特殊的图象恢复技术 投影：退化过程 （丢失了沿射线方向的分辨能力） 重建：复原过程 （恢复了2-D空间的分辨能力） 图象重建 1 投影重建概述 2 傅里叶反变换重建 3 逆投影重建 4 级数展开重建 5 综合重建方法 6 离散周期重建 1 投影重建概述 {P.10，表1.2.2} 1.1 投影重建方式 1.2 投影重建原理 利用投影重建方式工作的系统有许多种类，原理不完全相同但有一些共性 1.1 投影重建方式 1、透射断层成象? TCT（transmission computed tomography，CT） CAT 接受器获得的 射线强度实际上反 映了物体各部分对 射线的吸收情况 1.1 投影重建方式 2、发射断层成象 ECT（emission computed tomography） 发射源在物体内部，接受器在物体外部 检测到与生理有关的状况/信息 ① PET（positron emission tom