中南民族大学 数字图像处理实验报告 学院：电子信息工程学院 专业：光信息 学生 学号 实验一 图像变换 实验目的 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); 运 行 结 果： 2.对给定的一副图像satu rn2.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 运行结果： 3.在图像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(C>thresh); 运行结果： 0 5 1 0 4 .对一幅图进行离散余弦变换和反变换 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); -5 0 5 1 0 思考题 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)
