art算法图像重构.pptVIP

图像重建概述 图像重建是图像处理中的一个重要分支，广泛地应用于物体内部结构图像的检测和观察中，它是一种无损检测技术。 关于图像处理的一些基本内容，如对图像的几何处理，图像的增强，还有复原等，均是从图像到图像，即输入的原始数据是图像，处理后输出的仍是图像。而图像重建是从数据到图像。 图像重建的三种常用检测模型：透射模型、发射模型、反射模型 ---我们从CT谈起 计算机层析成像（Computed Tomography, CT）是通过对物体进行不同角度的射线投影测量来获取物体横截面信息的成像技术。 CT的核心技术是由投影数据来重建图像的理论，其实质是由扫描所得到的的投影数据来求出成像平面上每个点的衰减系数值。 二维投影与CT值 当强度为 的x-ray通过吸收率为μ(x,y)的均匀吸收物体，由于均匀吸收，则I必是指数下降，则有 CT任意角度扫描 Radon变换 Radon变换是计算图像在某一指定角度射线方向上的投影的变换方法。二维函数f(x,y)的投影是其在确定方向上的线积分，如下图所示，二维函数f(x,y)在水平方向的线积分就是f(x,y)在y轴上的投影，二维函数f(x,y)在垂直方向的线性积分就是f(x,y)在x轴上的投影。 Radon变换（续） 由此，可以沿任意角度 计算函数的投影，计算图像f(x,y)在任意角度的Radon变换。 中心切片定理 密度函数在某一方向上的投影函数的一维傅立叶变换函数是原密度函数的二维傅立叶变换函数在平面上沿同一方向且过原点的直线上值。 滤波反投影算法的原理 1、在不同的角度下取得足够多的投影数据(Radon变换) 2、将这些投影数据做一维的Fourier变换，那么变换后的这些数据将充满整个(u,v)平面。（许多过原点成不同夹角的直线） 3、也就是说，F(u,v)的全部值都为已知，那么我们将其做一次二维的Fourier逆变换就可以得到原始的衰减系数函数f(x,y) 二维傅立叶反变换 由傅立叶变换性质可知．频域中的滤波运算可等效地在空域中用卷积运算来完成 所以要实现对投影数据实现图像重建，可以采取两步：首先将投影数据和响应脉冲滤波器进行卷积，然后由式对不同旋转角θ求和，就能实现图像重建。这就是卷积法进行图像重建的基本思路和方法。 卷积可看作一种滤波手段，卷积投影相当于对数据先滤波再将结果逆投影回来，这样可以使模糊得到校正。 反投影算法举例 基本原理是将所测得的投影值按其原路径平均的分配到每一点上，各个方向上投影值反投影后，在影像处进行叠加，从而推体出原图像。 而滤波却是要投影函数的一维Fourier加上权重因子。 滤波反投影算法 滤波反投影法采用先修正、后反投影的做法，其基本方法是：在某一投影角下取得了投影函数（一维函数）后，对此一维投影函数作滤波处理，得到一个经过修正的投影函数；然后再将此修正后的投影函数作反投影运算，得到所需的密度函数。 滤波反投影法重建图像有以下几个步骤： （1）对某一角度下的投影函数作一维傅立叶变换； （2）对（1）的变换结果乘上一维权重因子； （3）对（2）的加权结果作一维逆傅立叶变换； （4）用（3）中得出的修正过的投影函数做直接反投影； （5）改变投影角度，重复（1）~（4）的过程，直到完成全部180度的反投影。 滤波函数 滤波函数的选取是滤波反投影法的关键问题 （1）R-L滤波函数 由于在频域中用矩形函数截断了滤波函数，在相应的空域中造成振荡响应 ，重建的图像质量也不够满意 （2）S-L滤波函数 与R-L滤波函数不同的是，S-L滤波函数它的关键是把频域的陡峭截止改成缓慢截止。 用S-L滤波函数重建的图像中振荡相应较小，对含噪声的数据重建出来的图像质量也较R-L滤波函数重建的图像质量要好。但是，S-L滤波函数重建的图像在高频响应方面不如R-L滤波函数好，这是因为S-L滤波函数在高频段偏离了理想的滤波函数 一、滤波反投影matlab实现 %P=imread(lena.jpg); P = phantom(256); %P= rgb2gray(O); R = radon(P,0:179); I0 = iradon(R,0:179,linear,Ram-Lak); I1 = iradon(R,0:179,linear,Shepp-Logan); I2=iradon(R,0:179,linear,cosine); I3 = iradon(R,0:179,linear,none); subplot(2,3,1), imshow(P), title(Original) subplot(2,3,2), imshow(I0,[]), title(FBP R-L) subplot(2,3,3), imsho