数字图像处理上机报告.docxVIP

数字图像处理报告 第二次上机报告 学院： 机电学院 班级： 071113-20 学号： 20111003341 姓名： 曾志鹏  2-1、计算图象的频谱函数 设计图象f4(x,y)为3*30*30/256*256，水平排列； 分析： 首先生成一个256*256的零矩阵，再按要求将其中某些部分赋值为一，即可得到图像。 程序如下： A=zeros(256); %生成一个256*256的零矩阵 for i=110:140 %设置行 for j=30:60 for p=100:130 for q=170:200 %设置列 A(i,j)=1; A(i,p)=1; A(i,q)=1; %将对应行列赋值为1, end end end end subplot(1,2,1); imshow(A); %显示得到图像 title(原图像); A1=fft2(A); %对图像进行傅里叶变换 A2=abs(A1); subplot(1,2,2); imshow(fftshift(log(A2))); title(傅里叶变换后图像); 运行结果如下： 2-2、根据计算证明傅立叶变换的性质 空域平移性: 设f1(x,y)为30*30/256*256, 左移得到f2(x,y), 求F2(u,v); 上移得到f3(x,y), 求F3(u,v); 证明F1(u,v)、F2(u,v)和F3(u,v)的绝对值相等。 分析：按照上题方法，设计图像并对图像进行傅里叶变换，再将设计好的图像的行与列的值分别改变，即可得到平移后的图像。将得到的平移后的图像在进行傅里叶变换，即可得到相应的频谱。设置变量a,当三个图的频谱相等时，a值为1，否则为0，观察a的值，即可判断呢所得结果。 程序如下： A=zeros(256); for i=100:130 for j=100:130 A(i,j)=1; end end %设置图像A A1=fft2(A); A2=abs(A1); subplot(2,3,1);imshow(A);title(原图像); subplot(2,3,4);imshow(fftshift(log(A2)));title(傅里叶变换后图像); B=zeros(256); for i=100:130 for j=50:80 B(i,j)=1; end end %设置图像B B1=fft2(B); B2=abs(B1); subplot(2,3,2);imshow(B);title(左移后图像); subplot(2,3,5);imshow(fftshift(log(B2)));title(左移后傅里叶变换图像); C=zeros(256); for i=50:80 for j=100:130 C(i,j)=1; end end %设置图像C C1=fft2(C); C2=abs(C1); subplot(2,3,3);imshow(C);title(右移后图像); subplot(2,3,6);imshow(fftshift(log(C2)));title(右移后傅里叶变换图像); [m n]=size(A2); for i=1:m for j=1:n if A2(i,j)==B2(i,j)A2(i,j)==C2(i,j) %判断频率绝对值是否相等 a=1; else a=0; end end end 结果如下： 由上图得知a值为1，即表明平移后的图像的频谱的绝对值相等，与题目意思符合。 2-3、图象变换比较 自行设计f(x,y)，（1）调用Matlab函数直接调用实现其离散傅立叶变换、离散余弦变换；（2）自行编程对f(x,y)实施Walsh变换和Hadamard变换，比较四种变换所得到的频谱。 分析：设计图像，将图像进行傅里叶变换和离散余弦变换。将得到的结果显示出来，由沃尔什变换和哈达玛变换，编写程序将算法正确的表达出来。 程序如下： N=16; n=log2(N); f=zeros(16); for i=4:8 for j=4:8 f(i,j)=1; end