资源描述：
《matlab 图像的傅立叶变换.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、实验三图像的傅立叶变换熟悉傅里叶变换的基本性质；掌握FFT方法及应用；通过实验了解二维频谱的分布特点；掌握利用MATLAB实现数字图像的傅立叶变换及滤波锐化的处理；了解理想、巴特沃兹、高斯等不同滤波器的结构及滤波效果。实验步骤1．启动MATLAB程序，读入一幅图像；对图像做FFT。使用’subplot’命令，同时显示原始图像其频谱图；IenaImg=imread('lena.jpg');%读入原图像文件fftI=fft2(double(IenaImg));%二维离散傅立叶变换sfftI=fftshift
2、(fftI);%直流分量移到频谱中心RR=real(sfftI);%取傅立叶变换的实部fftImg=log(RR);subplot(1,2,1)imshow(IenaImg);%显示原图像subplot(1,2,2)imshow(fftImg,[8,10]);%显示原图像的频谱1．读入一幅图像，分别为图像添加椒盐、高斯噪声，做FFT变换。使用’subplot’命令，将原始图像、原始图像频谱图、添加噪声后的图像，以及噪声图像的频谱图同时显示出来。lenaImg=imread('lena.jpg');%读入
3、原图像文件fftI=fft2(double(lenaImg));%二维离散傅立叶变换sfftI=fftshift(fftI);%直流分量移到频谱中心RR=real(sfftI);%取傅立叶变换的实部fftLenaImg=log(RR);加入椒盐躁声saltImg=imnoise(lenaImg,'salt&pepper',0.02);%加入椒盐躁声fftI=fft2(double(saltImg));%二维离散傅立叶变换sfftI=fftshift(fftI);%直流分量移到频谱中心RR=real(sf
4、ftI);%取傅立叶变换的实部fftSaltImg=log(RR);subplot(2,2,1);imshow(lenaImg);subplot(2,2,2);imshow(fftLenaImg,[8,10]);subplot(2,2,3);imshow(saltImg);subplot(2,2,4);imshow(fftSaltImg,[8,10]);加入椒盐躁声加入高斯躁声gaussianImg=imnoise(lenaImg,'gaussian',0.02);%加入高斯躁声fftI=fft2(do
5、uble(gaussianImg));%二维离散傅立叶变换sfftI=fftshift(fftI);%直流分量移到频谱中心RR=real(sfftI);%取傅立叶变换的实部fftGaussianImg=log(RR);subplot(2,2,1);imshow(lenaImg);subplot(2,2,2);imshow(fftLenaImg,[8,10]);subplot(2,2,3);imshow(gaussianImg);subplot(2,2,4);imshow(fftGaussianImg,[
6、8,10]);加入高斯躁声1．读入一幅图像，对图像分别进行高斯低通、巴特沃兹低通、高斯高通和巴特沃兹高通频域滤波，比较其锐化和平滑效果；%高斯低通滤波器I=imread('lena.jpg');imshow(I);s=fftshift(fft2(I));[M,N]=size(s);n=2;d0=30;%GLPF滤波，d0=5，15，30（程序中以d0=30为例）n1=floor(M/2);n2=floor(N/2);fori=1:Mforj=1:Nd=sqrt((i-n1)^2+(j-n2)^2);h=
7、1*exp(-1/2*(d^2/d0^2));s(i,j)=h*s(i,j);endends=ifftshift(s);gImg=uint8(real(ifft2(s)));figure;imshow(gImg);%巴特沃思低通滤波器l=imread('test.jpg');s=fftshift(fft2(l));[M,N]=size(s);%5级巴特沃思n=5;d0=30;%BLPF滤波，d0=5，15，30（程序中以d0=30为例）n1=floor(M/2);n2=floor(N/2);fori=1
8、:Mforj=1:Nd=sqrt((i-n1)^2+(j-n2)^2);h=1/(1+(d/d0)^(2*n));s(i,j)=h*s(i,j);endends=ifftshift(s);b5Img=uint8(real(ifft2(s)));%15级巴特沃思n=15;s=fftshift(fft2(l));[M,N]=size(s);d0=30;%BLPF滤波，d0=5，15，30（程序中以d0=30为例）n1=floor(M/2);n2
