第七章图象重建.pptVIP

第7章 图象重建 一、傅立叶变换重建 一个三维（或二维）物体，它的一维（或 一维）投影的傅立叶变换恰与此物体的傅立叶 变换的主体部分相等，而傅立叶变换重建方法 也正是以此为基础的。 通过将投影进行旋转和傅立叶可以首先构 造整个傅立叶变换的平面，然后只须在通过傅 立叶反变换就可以得到重建后的物体。 原理如下： 令f ( x ,y )代表一图象函数，则傅立叶变换为： 图象在轴上的投影为： 投影的一维傅氏变换为： 由上式可见，它与二维傅氏变换的表达式一致。即： 假设将函数投影到一条经过旋转的直线上，该直线的 旋转角度为 将函数投影的直线选为x轴。投影点通过距离t轴为s1处的一平行线进行函数积分 该投影的一维傅氏变换为： 反变换为： 上述结论可以推广到三维： 变换的核心部分为： §2 卷积法重建 先看下图的极坐标中傅立叶反变换表达式： 笛卡尔坐标和极坐标的关系如下： 由对称共轭特性： 令 这里： 则： ＊代表卷积运算，此卷积表达式可写为： 另外，求导也可推出一种很简单的图像 重建方法。假设将投影数据 都存放于一 等量矩形空间内，这种存放数据的方式称为 Layergram。 对于一恒定 值，我们可线性地滤出该 投影数据，即可在频域内用Rho滤波器乘以 得出。也可以在空间域内通过一个滤波器冲 激响应是Rho频率滤波器的反变换的投影数据 卷积得出，表达式为： 在一特定 值，通过Rho 滤波，对Layergram 以 积分变量作积分，计算出反投影。即： 这里， §3 代数重建方法 基本投影公式： 由 可以用级数来表示函数估值： 级数估值应满足投影方程: 或 我们希望仅通过一系列离散的点来得到估 计图像： 为了说明这些数字化重建中遇到的问题及 解决方案，在介绍一下傅立叶变换的方法。 假设待重建的图像是空间有限的，这样直接应用抽样理论，首先可推出f( x ,y )的周期延拓函数fp(x,y) 在变换域重用等间隔栅格上的函数插值方法 表示的傅立叶变换为： 其中： 等于在 点的傅立叶变 换抽样值， 在给出傅立叶变换在极坐标中计算出的投 影点后，我们要确定系数 表示成极坐标形式： 可以写出用等值向量形式表示的有限插入公 式： F=WG F是一组下标为IJ的矢量，G是下标为MN的矢量。W是按IJ×MN排列的加权矩阵。 则逆矩阵为： 可以看出：傅立叶变换法的数字计算可以 说是一种级数法。 另外一种数字法是基于f(x,y)是带限函数这 一假设的基础上的方法。即： F( u ,v)=0 沿直线上的积分为： 直线方称为 y=ax+b 通过定义 可以方便地在任意一个取样点重新定位原点。 对于一个给定的 和 ，直线方称为： 其中， 为对上式进一步求值，得到加权函数为： 这里， 我们可以得到一系列到数表达式： §4 重建的优化问题 图像重建中的问题也可以通过选择一合理 的准则很数来解决。此函数