14组实验七概论.doc

实验七 医学图像频域滤波与图像复原 实验目的： 熟悉医学图像离散傅里叶变化的原理和方法； 掌握医学图像频域滤波的原理； 掌握使用Matlab中的函数实现医学图像进行频域滤波的方法; 掌握使用Matlab中的图像退化与复原的方法; 实验内容： 一、医学图像频域滤波方法与实现 使用imnoise给图像BMRI1_24bit.bmp添加概率为0.2的椒盐噪声，对原图像和加噪声后的图像进行傅氏变换并显示变换后的移中频谱图，然后分别使用Butterworth低、高通滤波器对噪声图像进行低通和高通滤波，显示D0为5，10，20，40时的滤波效果图，并说明存两种滤波效果中所存在的差异及原因。 答： f=imread(BMRI1_24bit.bmp); f=rgb2gray(f); g=imnoise(f,salt pepper,0.2); f=double(f); g=double(g); F1=fft2(f); F2=fft2(g); FC1=fftshift(F1); FC2=fftshift(F2); f=ifft2(F1); g=ifft2(F2); imshow(uint8(f));title(逆变换后原图像); figure,imshow(uint8(g));title(逆变换后的椒盐噪声原图像); figure,imshow(F1);title(原图像