资源描述：
《数字图像处理实验图像频谱分析.doc》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、姓名：朱慧娟班级：电子二班学号：410109060325实验2图像频谱分析一、实验目的1、了解图像变换的意义和手段。2、熟悉及掌握图像的变换原理及性质，实现图像的傅里叶变换。二、实验内容1、分别显示图像Bridge.bmp、cameraman.tif（自带图像）、blood.tif及其频谱，观察图像频谱的特点。2、生成一幅图像，图像中背景黑色，目标为一亮条；平移亮条，观察其频谱的变化。3、对lena.bmp图像进行旋转，显示原始图像与旋转后图像，及其傅里叶频谱，分析旋转前、后傅里叶频谱的对应关系。三、实验程序及结果1.1实验程序cle
2、ar;%清除以前实验变量a=imread('e:ZHJBridge.bmp');%读入图像Bridge.bmp，并记为ab=imread('cameraman.tif');%读入图像cameraman.tif，并记为bc=imread('e:ZHJblood.tif');%读入图像blood.tif，并记为cd=fft2(a);%对图像a进行傅里叶变换，并记为de=fftshift(d);%将变换后图像频谱中心从矩阵的原点移动到矩阵的中心，并记为eA=abs(e);%对e取绝对值，及得到图像a的幅度谱，并记为AB=log(1+
3、A);%对幅度谱A取对数，并记为Bf=fft2(b);%对图像b进行傅里叶变换，并记为fg=fftshift(f);%将变换后图像频谱中心从矩阵的原点移动到矩阵的中心，并记为gC=abs(g);%对g取绝对值，及得到图像b的幅度谱，并记为CD=log(1+C);%对幅度谱C取对数，并记为Dh=fft2(c);%对图像c进行傅里叶变换，并记为hi=fftshift(h);%将变换后图像频谱中心从矩阵的原点移动到矩阵的中心，并记为iE=abs(i);%对i取绝对值，及得到图像c的幅度谱，并记为EF=log(1+E);%对幅度谱E取对数，并
4、记为Ffigure(1);%建立图表1subplot(2,1,1);%将图表1分成两部分，第一部分imshow(a);%显示图像atitle('Bridge.bmp');%给图像a加标题‘Bridge.bmp’subplot(2,1,2);%将图表1分成两部分，第二部分imshow(B,[]);%显示B即图像a的频谱图title('Bridge.bmp频谱图');%给图像B加标题‘Bridge.bmp频谱图’figure(2);%建立图表2subplot(2,1,1);%将图表2分成两部分，第一部分imshow(b);%显示图像bti
5、tle('cameraman.tif');%给图像b加标题‘cameraman.tif’subplot(2,1,2);%将图表2分成两部分，第二部分imshow(D,[]);%显示D即图像b的频谱图title('cameraman.tif频谱图');%给图像D加标题‘cameraman.tif频谱图’figure(3);%建立图表3subplot(2,1,1);%将图表3分成两部分，第一部分imshow(c);%显示图像ctitle('blood.tif');%给图像c加标题‘blood.tif’subplot(2,1,2);%将图表
6、3分成两部分，第二部分imshow(F,[]);%显示F即图像c的频谱图title('blood.tif频谱图');%给图像F加标题‘blood.tif频谱图’1.2实验结果2.1实验程序clear;%清除以前实验变量A=zeros(256,256);%建立行列都是256的0矩阵，即建立黑色图，并记为AA(10:20,:)=256;%矩阵A中第十到二十行数据改为256，即在黑色图像上加上亮条纹B=circshift(A,[50,0]);%将矩阵A行向移动50行，得到新矩阵记为Ba=fft2(A);%对矩阵A进行傅里叶变换，并记为ab=
7、fftshift(a);%将变换后图像频谱中心从矩阵的原点移动到矩阵的中心，并记为bM=abs(b);%对b取绝对值，及得到矩阵A的幅度谱，并记为MN=log(1+M);%对幅度谱M取对数，并记为Nc=fft2(B);%对矩阵B进行傅里叶变换，并记为cd=fftshift(c);%将变换后图像频谱中心从矩阵的原点移动到矩阵的中心，并记为dS=abs(d);%对d取绝对值，及得到矩阵B的幅度谱，并记为ST=log(1+S);%对幅度谱S取对数，并记为Tfigure;%建立图表subplot(2,2,1);%将图表分成四部分，第一部分im
8、show(A);%显示图像Atitle('原图像');%给所显示图像加标题‘原图像’subplot(2,2,2);%将图表分成四部分，第二部分imshow(B);%显示图像Btitle('平移后图像');%给所显示图像加
