数字图像处理报告90754.doc

中南民族大学 数字图像处理实验报告 学院：电子信息工程学院 专业：光信息 学生 学号 实验一 图像变换 实验目的 1：了解二维fourier变换的原理 2：掌握二维fourier变换的性质 3：掌握离散余弦变换的原理 实验内容 1.绘制一个二值，并将其傅里叶函数可视化。 close all; clear; f=zeros(40,40); f(10:30,10:30)=1; subplot(1,3,1); imshow(f); F=fft2(f); subplot(1,3,2); imshow(F); D=log(1+abs(F)); subplot(1,3,3); imshow(D); 运行结果： 对给定的一副图像saturn2.tif进行傅里叶变换 load imdemos saturn2; subplot(1,2,1); imshow(saturn2); I=fftshift(fft2(saturn2)); subplot(1,2,2); imshow(log(abs(I)),[]),colormap(jet(64)),colorbar 运行结果： 在图像text,tif中定位字母a w=imread(text.png); a=w(33:45,88:98); figure; imshow(w); figure; imshow(a); C=real(ifft2(fft2(w).*fft2(rot90(a,2),256,256))); figure; imshow(C,[]); max(C(:)); thresh=60; figure; imshow(Cthresh); 运行结果： 对一幅图进行离散余弦变换和反变换 RGB=imread(autumn.tif); figure(1),imshow(RGB); I=rgb2gray(RGB); figure(2),imshow(I); J=dct2(I); figure(3),imshow(log(abs(J)),[]),colormap(jet(64)),colorbar; J(abs(J)10)=0; K=idct2(J)/255; figure(4),imshow(K); 思考题 1.编写M程序文件验证二维傅里叶变换的平移性，旋转性 clear all; image1=zeros(128,128); image1(30:93,52:67)=1; subplot(2,2,1); imshow(image1,[]); title(原始图像); image2=abs(fft2(image1)); subplot(2,2,2); imshow(log10(1+image2),[]); title(频谱); image3=image1; i=size(image1,1); j=size(image2,2); for k=1:i for r=1:j image3(k,r)=(-1)^(k+r).*image1(k,r); end end image4=abs(fft2(image3)); subplot(2,2,3); imshow(image3,[]) title(中心化后的图像); subplot(2,2,4); imshow(log10(1+image4),[]); title(中心化后的频谱); image5=zeros(i,j); i=round(i/2) j=round(j/2) for u=-i:i for v=-j:j m=cos(pi/4)*u-sin(pi/4)*v; n=sin(pi/4)*u+cos(pi/4)*v; m=round(m); n=round(n); if abs(m)iabs(n)j image5(u+i+1,v+j+1)=image1(m+i,n+j); end end end figure subplot(2,1,1); imshow(image5,[]); image6=abs(fft2(image5)); subplot(2,1,2); imshow(log10(1+image6),[])； 运行结果 实验二 图像的空间变换增强 实验目的 1：掌握图像的直方图均衡化的原理与应用 2：掌握灰度变换的原理与应用 3：掌握伪彩色图像处理增强的方法 实验内容 MATLAB提供的直方图修正函数 I=imread(pout.tif); subplot(221); imshow(I); subplot(222); imhist(I);