图像工程,图像处理,章毓晋IE1IP09.ppt

第9章 图象重建 第9章 图象重建 9.1 投影重建概述 9.1.1 投影重建方式 9.1.1 投影重建方式 9.1.1 投影重建方式 9.1.1 投影重建方式 9.1.1 投影重建方式 9.1.1 投影重建方式 9.1.1 投影重建方式 9.1.1 投影重建方式 9.1.2 投影重建原理 9.1.2 投影重建原理 9.1.2 投影重建原理 9.1.2 投影重建原理 9.2 傅里叶反变换重建 9.2 傅里叶反变换重建 9.2 傅里叶反变换重建 9.2 傅里叶反变换重建 9.2 傅里叶反变换重建 9.2 傅里叶反变换重建 9.2 傅里叶反变换重建 9.2 傅里叶反变换重建 9.2 傅里叶反变换重建 9.2 傅里叶反变换重建 9.2 傅里叶反变换重建 9.3 逆投影重建 9.3.1 卷积逆投影重建 9.3.1 卷积逆投影重建 9.3.1 卷积逆投影重建 9.3.1 卷积逆投影重建 9.3.1 卷积逆投影重建 9.3.1 卷积逆投影重建 9.3.2 其他逆投影重建 9.3.2 其他逆投影重建 9.4 级数展开重建 9.4 级数展开重建 9.4 级数展开重建 9.4 级数展开重建 9.5 综合重建方法 9.5 综合重建方法 9.5 综合重建方法 联 系 信 息 2、离散计算 逆投影过程的近似 卷积过程的近似 离散卷积 一次插值 3、扇束投影重建 缩短投影时间，使用1个发射器和1组接受器 {式(9.3.9)到式(9.3.12)} {式(9.3.14)到式(9.3.16)} 3、扇束投影重建 (a) 发射器和接受器均旋转 (b) 发射器旋转，接受器固定 4、傅里叶反变换和卷积逆投影重建法的比较 傅里叶反变换重建 ：直角坐标 卷积逆投影重建 ：极 坐标 主要优点： 卷积逆投影重建基本算法： 容易用软件和硬件实现 可重建出准确清晰的图象 傅里叶反变换重建： 计算量比较小 1、逆投影滤波 逆投影滤波器（filter of the back-projections） 先逆投影，再滤波（/卷积） 2、滤波逆投影 （backprojection of the filtered projections） 先对投影滤波，再逆投影 1、离散模型 代数重建法（ART），迭代算法，优化技术 第i条射线与 第j个象素相交 的长度aij，代 表第j个象素沿 第i条射线的贡 献的权值 2、无松弛的代数重建技术 先初始化1个图象矢量x(0)，然后用下式进行迭代： 思路： 每次取1条射线，改变图象中与该线相交的象素的值，从而把当前的图象矢量x(k)更新为x(k+1)。具体运算中就是将测量值与由当前算得的投影数据 的差正比于aij重新分配到各个象素上去 3、松弛的代数重建技术 将无松弛的代数重建技术推广： r是松弛系数，它控制了收敛的速度 重建结果示例 {图9.4.3} 重建误差： 图象的4个角处（较少的射线） 内部正方形的4个角处（吸收值急剧变化） 4、 级数法的一些特点 (1) 由于在空域中比较灵活，所以容易调整以适应新的应用 (2) 能重建出较高对比度的图象（对密度突变的材料） (3) 借助多次迭代可从较少投射（ 10）重建图象 (4) 比变换法更适合于ECT (5) 比变换法更适合于3-D重建问题 (6) 比变换法更适合于不完整投影情况 迭代变换法（iterative transform） L(s, ?n)是直线(s, ?n)与Q相交段的长度 重建工作： 给定函数 f (x, y) 的投影g(s, ?n)，其 中s可以取所有实 数，?n是1组N个离 散角度，要重新得到 f (x, y) 迭代变换法（iterative transform） 迭代公式 n = (i mod N ) +1，s = xcos?n + ysin?n g(i)(s, ?n )是 f (i)(x, y)沿?n的投影 图象序列 将收敛 *章毓晋 (TH-EE-IE) 第*页 第9讲 章毓晋 清华大学电子工程系 100084 北京 图象工程 投影