资料频域滤波课程设计.docVIP

1引言 当前，人类社会已经进入了信息社会。信息是事物运动状态和特征的反映，它和材料及能量一起构成社会的三个要素。但是，信息具有一些不同于材料和能量的特征。信息具有普遍性，无损性，时空独立性，等等。正是由于信息具有这些特征，因此，它和人类文明及社会发展的各个阶段都有密切的联系。 与其他学科一样，生物医学受到信息科学及技术的影响，进而相互渗透，融合。由此产生的对生物和医学学科发展的促进作用也是显而易见的：例如，X射线早在1895年已被发现，然而，只有在计算机技术快速发展和图像重建及处理方法得以实现的基础上，sfieack才发明了X射线计算机断层扫描成像技术（XCT）并应用于医学。 图像重建方法包括二维平行光束和扇形束成像，三位平行线，平行面以及锥形束成像。我们所用到的傅立叶变换就属于二维平行光束成像。一般地说，人眼所能看到的实物就是处于一个相对的空间参考系中，而在空间域中的图像就是对物体光反射特性的一种具体表现。因此，在空间域中定义的图像，其像素的坐标值确定了该像素的相对空间位置，图像中各像素的灰度值反映了物体在该点的光反射特性。正式由于图像的空间域定义直接反也容易让人接受。 但是，有时为了使某种处理或运算更为简便或快速，或者为了更明显的表现图像的另一些特征也可以在其他域中来描述和处理图像，广义地将其称为变换域（transformdomain .比如，刚获得的图像有很多噪音。这主要由于平时的工作和环境引起的，图像增强是减弱噪音，增强对比度。想得到比较干净清晰的图像并不是容易的事情。B fft2 I B fft2 I 返回图象I的二维fft变换矩阵，输入图象I和输出图象B大小相同。 （2）fftshift MATLAB提供的fftshift函数用于将变换后的图象频谱中心从矩阵的原点移到矩阵的中心，其语法格式为： 　　B fftshift I 对于矩阵I，B fftshift I 将I的一、三象限和二、四象限进行互换。 （3）ifft2 ifft2函数用于计算图象的二维傅立叶反变换，其语法格式为： B ifft2 I B ifft2 A 返回图象I的二维傅立叶反变换矩阵，输入图象I和输出图象B大小相同。其语法格式含义与fft2函数的语法格式相同，可以参考fft2函数的说明。 4.2读取任一幅灰度图像； 计 I imread 'pout.tif' ; figure 1 ; subplot 1,2,1 ,imshow I title '原图像' ; % 求离散傅立叶频谱 J fftshift fft2 I ; subplot 1,2,2 ,imshow log abs J ,[8,10] title '图像傅里叶变换所得频谱' ; 像 析 原图像经过傅立叶变换得到了该图像的频谱图像。 4.3实现频域理想低通滤波，（对比原图，原图频域信号，滤波后图，滤波后频域信号，给出滤波器幅度特性）； 通滤波器函数 一个二维的理想低通滤波器（ILPF）的转换函数满足（是一个分段函数） 其中：D0 为截止频率 D u,v 为距离函数 D u,v u2+v2 1/2 通滤波器幅度特性 计 I imread 'pout.tif' ; figure 1 ; subplot 2,2,1 ,imshow I title '原图像' ; % 求离散傅立叶频谱 J fftshift fft2 I ; subplot 2,2,2 ,imshow log abs J ,[8,10] title '图像傅里叶变换所得频谱' ; [a,b] size J ; a0 round a/2 ; b0 round b/2 ; d 10;%改变d值，所得结果不同； for i 1:a for j 1:b distance sqrt i-a0 ^2+ j-b0 ^2 ; if distance d h 1; else h 0; J i,j h*J i,j ; end; end; end; J1 uint8 real ifft2 ifftshift J ; subplot 2,2,3 ,imshow J1 ; title '低通滤波所得图像' ; J2 fftshift fft2 J1 ; subplot 2,2,4 ,imshow log abs J ,[8,10] title '低通滤波所得频谱' ; 像 析 由于傅立叶变换的实部R u,v 及虚部I u,v 随着频率u,v的升高而迅速下降，所以能量随着频率的升高而迅速减小，因此在频域平面上能量集中于频率很小的圆域内, 当D0