数字图像处理试验汇报(全答案)

《数字图像处理试验汇报(全答案)》由会员分享，可在线阅读，更多相关《数字图像处理试验汇报(全答案)（29页珍藏版）》请在装配图网上搜索。

1、数字图像处理试验汇报(全答案) 试验一常见BmATLAB图像处理命令一、试验目标1、熟悉并掌握mATLAB工具的使用；2、实现图像的读取、显示、代数运算和简单变换。二、试验环境以上版本、wINXP或wIN2021计算机三、常见函数l读写图像文件1imreaddimread函数用于读入多种图像文件,如：a=imread(e:)2imwriteeimwrite函数用于写入图像文件,如：imwrite(a,e:,tif)3imfinfooimfinfo函数用于读取图像文件的相关信息,如：imfinfo(e:)l图像的显示11imageeimage函数是BmATLAB提供的最原始的图像显示函数,如：a

2、=1,2,3,4;4,5,6,7;8,9,10,11,12;image(a);2imshowwimshow函数用于图像文件的显示,如：i=imread(e:);imshow(i);title（原图像）%加上图像标题3colorbarrcolorbar函数用显示图像的颜色条,如：i=imread(e:);imshow(i);colorbar;4figureefigure函数用于设定图像显示窗口,如：figure(1)；/figure(2)；5subplot把图形窗口分成多个矩形部分,每个部分能够分别用来进行显示。SSubplot（m,n,p）分成nm*n个小窗口,在第pp个窗口中创立坐标轴为目前

3、坐标轴,用于显示图形。6plot绘制二维图形框架三层框架三层,局部为二层钢构。本工程外脚手架采取落地式钢管脚手架局部为二层钢构。本工程外脚手架采取落地式钢管脚手架,外架随主体结构上升外架随主体结构上升,同时搭设同时搭设,比操作面高出一步比操作面高出一步,确保主体及外装修的正常安全施工。确保主体及外装修的正常安全施工。plot（yy）PPlot（x,y）yxy能够是向量、矩阵。l图像类型转换1rgb2gray/灰色把真彩图像转换为灰度图像i=rgb2gray（jj）2im2bw/黑白经过阈值化方法把图像转换为二值图像II=im2bw（j,level）LLlevel表示灰度阈值,取值范围01（即）

4、,表示阈值取自原图像灰度范围的n%3imresize改变图像的大小II=imresize（j,mn）将图像jj大小调整为mm行行nn列l图像运算1imadd两幅图像相加,要求一样大小,同种数据类型ZZ=imadd（x,y）表示图像x+y2imsubstract两幅图像相减,要求一样大小,同种数据类型ZZ=imsubtract（x,y）表示图像xx-yy3immultiplyZZ=immultiply（x,y）表示图像x*y4imdivideZZ=imdivide（x,y）表示图像x/y5：m=imadjust(a,;1);%图像变亮n=imadjust(a,0;);%图像变暗g=255-a;%

5、负片效果四、试验内容（请将试验程序填写在下方适宜的位置,试验图像结果拷屏粘贴）1、读入一幅RGB图像,变换为灰度图像和二值图像,并在同一个窗口内分成三个子窗口来分别显示RGB图像和灰度图像,注上文字标题。a=imread(f:)i=rgb2gray(a)I=im2bw(a,)subplot(3,1,1);imshow(a);title(原图像)subplot(3,1,2);imshow(i);title(灰度图像)subplot(3,1,3);imshow(I);title(二值图像)3/252、对两幅不一样图像实施加、减、乘、除操作,在同一个窗口内分成五个子窗口来分别显示,注上文字标题。a=

6、imread(f:)A=imresize(a,800800)b=imread(f:)B=imresize(b,800800)Z1=imadd(A,B)Z2=imsubtract(A,B)Z3=immultiply(A,B)Z4=imdivide(A,B)subplot(3,2,1);imshow(A);title(原图像A)subplot(3,2,2);imshow(B);title(原图像B)subplot(3,2,3);imshow(Z1);title(加法图像)subplot(3,2,4);imshow(Z2);title(减法图像)subplot(3,2,5);imshow(Z3);t

7、itle(乘法图像)subplot(3,2,6);imshow(Z2);title(除法图像)3、对一幅图像进行灰度改变,实现图像变亮、变暗和负片效果,在同一个窗口内分成四个子窗口来分别显示,注上文字标题。a=imread(f:);m=imadjust(a,;1);%图像变亮n=imadjust(a,0;);%图像变暗g=255-a;%负片效果subplot(2,2,1);imshow(a);title(原图像)subplot(2,2,2);imshow(m);title(图像变亮)subplot(2,2,3);imshow(n);title(图像变暗)subplot(2,2,4);imsho

8、w(g);title(负片效果)框架三层框架三层,局部为二层钢构。本工程外脚手架采取落地式钢管脚手架局部为二层钢构。本工程外脚手架采取落地式钢管脚手架,外架随主体结构上升外架随主体结构上升,同时搭设同时搭设,比操作面高出一步比操作面高出一步,确保主体及外装修的正常安全施工。确保主体及外装修的正常安全施工。4、熟悉数字图像处理常见函数的使用,调出帮助文档查看其多种不一样使用方法。方法：选择函数（函数所在区变暗）,点右键弹出菜单,选择“HelponSelection”五、试验总结分析图像的代数运算结果,分别陈说图像的加、减、乘、除运算可能的应用领域。5/25试验四常见图像增强方法一、试验目标1、熟

9、悉并掌握mATLAB图像处理工具箱的使用；2、了解并掌握常见的图像的增强技术。二、试验环境以上版本、wINXP或wIN2021计算机三、相关知识1imnoiseimnoise函数用于对图像生成模拟噪声,如：i=imread(e:);j=imnoise(i,gaussian,0,);模拟均值为0方差为的高斯噪声,j=imnoise(i,salt&pepper,)模拟叠加密度为的椒盐噪声2fspecialfspecial函数用于产生预定义滤波器,如：h=fspecial(sobel);%sobel水平边缘增强滤波器h=fspecial(gaussian);%高斯低通滤波器h=fspecial(la

10、placian);%拉普拉斯滤波器h=fspecial(log);%高斯拉普拉斯（LoG）滤波器h=fspecial(average);%均值滤波器3基于卷积的图像滤波函数imfilter函数,filter2函数,二维卷积conv2滤波,全部可用于图像滤波,使用方法类似,如：i=imread(e:);h=1,2,1;0,0,0;-1,-2,-1;%产生Sobel算子的水平方向模板j=filter2(h,i);或：h=fspecial(prewitt)I=imread();imshow(I);H=fspecial(prewitt);%预定义滤波器m=imfilter(I,H);imshow(m)

11、或：i=imread(e:);h=1,1,1;1,1,1;1,1,1;h=h/9;j=conv2(i,h);4其它常见滤波举例框架三层框架三层,局部为二层钢构。本工程外脚手架采取落地式钢管脚手架局部为二层钢构。本工程外脚手架采取落地式钢管脚手架,外架随主体结构上升外架随主体结构上升,同时搭设同时搭设,比操作面高出一步比操作面高出一步,确保主体及外装修的正常安全施工。确保主体及外装修的正常安全施工。（1）中值滤波medfilt2函数用于图像的中值滤波,如：i=imread(e:);j=medfilt2(i,mN);对矩阵i进行二维中值滤波,领域为m*N,缺省值为3*3（2）利用拉氏算子锐化图像,

12、如:i=imread(e:);j=double(i);h=0,1,0;1,-4,0;0,1,0;%拉氏算子k=conv2(j,h,same);三、试验步骤1、采取二维中值滤波函数medfilt2对受椒盐噪声干扰的图像滤波,窗口分别采取3*3,5*5,7*7I=imread(f:);j=imnoise(I,salt&pepper,);k1=medfilt2(j,33);%对矩阵i进行二维中值滤波,领域为3*3k2=medfilt2(j,55);k3=medfilt2(j,77);subplot(2,2,1);imshow(j);title(椒盐噪声干扰图像)subplot(2,2,2);imsh

13、ow(k1);title(领域为3*3二维中值滤波)subplot(2,2,3);imshow(k2);title(领域为5*5二维中值滤波)subplot(2,2,4);imshow(k3);title(领域为7*7二维中值滤波)2、采取mATLAB中的函数filter2对受噪声干扰的图像进行均值滤波I=imread(f:);j=imnoise(i,gaussian,0,);%模拟均值为0方差为的高斯噪声,m=filter2(fspecial(average,9),j)/255;%模板尺寸为9subplot(2,1,1);imshow(j);title(噪声干扰图像)subplot(2,1,

14、2);imshow(m);title(改善后的图像)3、采取三种不一样算子对图像进行锐化处理。i=imread(f:)7/25I=rgb2gray(s)H=fspecial(sobel)%应用Sobel算子锐化图像I1=filter2(H,I)%Sobel算子滤波锐化H=fspecial(prewitt)%应用prewitt算子锐化图像I2=filter2(H,I)%prewitt算子滤波锐化H=fspecial(log)%应用log算子锐化图像I3=filter2(H,I)%log算子滤波锐化subplot(2,2,1);imshow(i);title(原图像)subplot(2,2,2);

15、imshow(I1);title(Sobel算子锐化图像)subplot(2,2,3);imshow(I2);title(prewitt算子锐化图像)subplot(2,2,4);imshow(I3);title(log算子锐化图像)四、试验总结1、比较不一样平滑滤波器的处理效果,分析其优缺点2、比较不一样锐化滤波器的处理效果,分析其优缺点框架三层框架三层,局部为二层钢构。本工程外脚手架采取落地式钢管脚手架局部为二层钢构。本工程外脚手架采取落地式钢管脚手架,外架随主体结构上升外架随主体结构上升,同时搭设同时搭设,比操作面高出一步比操作面高出一步,确保主体及外装修的正常安全施工。确保主体及外装修

16、的正常安全施工。试验五图像恢复和图像分割一、试验目标1、熟悉并掌握mATLAB图像处理工具箱的使用；2、了解并掌握常见的图像的恢复和分割技术。二、试验环境以上版本、wINXP或wIN2021计算机三、相关知识1deconvwnr维纳滤波,使用方法：j=deconvwnr(I,PSF,NSR)用维纳滤波算法对图片I进行图像恢复后返回图像j。I是一个N维数组。PSF是点扩展函数的卷积。NSP是加性噪声的噪声对信号的功率比。如：I=im2double(imread();imshow(I);title(originalImage);%模拟运动模糊matlab汉字论坛LEN=21;THETA=11;PS

17、F=fspecial(motion,LEN,THETA);blurred=imfilter(I,PSF,conv,circular);figure,imshow(blurred)%恢复图像wnr2=deconvwnr(blurred_noisy,PSF);figure,imshow(wnr2)title(RestorationofBlurred)2edge检测灰度或二值图像的边缘,返回一个二值图像,1像素是检测到的边缘,0像素是非边缘。使用方法：Bw=edge(I,sobel,thresh,direction),I为检测对象；边缘检测算子可用sobel,roberts,prewitt,zero

18、cross,log,canny；thresh指定阈值,检测时忽略全部小于阈值的边缘,默认自动选择阈值；direction方向,在所指定的方向direction上,用算子进行边缘检测horizontal(水平方向)、vertical(垂直方向)或both(两个方向)。如：I=imread();Bw1=edge(I,prewitt);imshow(Bw1);3strel创立形态学结构元素。使用方法：SE=STREL(arbitrary,NHooD,HEIGHT)创立一个指定领域的非平面结构化元素。HEIGHT是一个矩阵,大小和NHooD相同,她指定了NHooD中任何非零元素的高度值。SE=STRE

19、L(ball,R,H,N)创立一个空间椭球状的结构元素,其X-y平面半径为R,高度为H。R必需为非负整数,H是一个实数。N必需为一个非负偶数。当N0时此球形结构元素由一系列空间线段结构元素来近似。SE=STREL(diamond,R)创立一个指定大小R平面钻石形状的结构化元素。R是从结构化元素原点到其点的距离,必需为非负整数。SE=STREL(disk,R,N)创立一个指定半径R的平面圆盘形的结构元素。这里R必需是非负整数.N须是0,4,6,8.当N大于0时,圆盘形结构元素由一组N(或N+2)个周期线结构元素来近似。当N等于0时,不使用近似,即结构元素的全部像素是由到中心像素距离小于等于R的像

20、素组成。N能够被忽略,此时缺省值是4。注:形态学操作在N0情况下要快于N=0的情形。如：se1=strel(square,11)%11乘以11的正方形4imerode9/25腐蚀图像使用方法：Im2=imerode(Im,SE)腐蚀灰度、二进制或压缩二进制图像Im,返回腐蚀图像Im2。参数SE是函数strel返回的一个结构元素体或是结构元素体阵列。如：使用一个盘状结构元素腐蚀一幅二进制图像。originalBw=imread();se=strel(disk,11);erodedBw=imerode(originalBw,se);imshow(originalBw),figure,imshow(

21、erodedBw)5imdilate膨胀图像使用方法：Im2=imdilate(Im,SE)膨胀灰度图像、二值图像、或打包的二值图像Im,返回膨胀图像m2。变量SE是一个结构元素或一个结构元素的数组,其是经过strel函数返回的。如：利用一个运行结构元素膨胀灰度图像。I=imread();se=strel(ball,5,5);I2=imdilate(I,se);imshow(I),title(original)figure,imshow(I2),title(Dilated)三、试验步骤1、产生运动模糊图像,利用维纳滤波进行图像恢复,显示结果。i=imread(f:)I=rgb2gray(s)I

22、=im2double(I);%模拟运动模糊LEN=21;THETA=11;PSF=fspecial(motion,LEN,THETA);blurred=imfilter(I,PSF,conv,circular);%恢复图像wnr2=deconvwnr(blurred,PSF);subplot(1,2,1);imshow(blurred);title(运动模糊图像)subplot(1,2,2);imshow(wnr2);title(恢复图像)2、采取三种不一样算子检测图像边缘,显示结果。i=imread(f:)Bw1=edge(I,prewitt);Bw2=edge(I,zerocross);B

23、w3=edge(I,canny);subplot(2,2,1);imshow(i);title(原图像)subplot(2,2,2);imshow(Bw1);title(prewitt边缘图)subplot(2,2,3);imshow(Bw2);title(zerocross边缘图)subplot(2,2,4);imshow(Bw3);title(canny边缘图)框架三层框架三层,局部为二层钢构。本工程外脚手架采取落地式钢管脚手架局部为二层钢构。本工程外脚手架采取落地式钢管脚手架,外架随主体结构上升外架随主体结构上升,同时搭设同时搭设,比操作面高出一步比操作面高出一步,确保主体及外装修的正常

24、安全施工。确保主体及外装修的正常安全施工。3、对二值图像分别进行方形模板3*3和5*5的膨胀和腐蚀操作,显示结果。a=imread(f:)i=rgb2gray(a)I=im2bw(a,)se3=strel(disk,3);erodedBw1=imerode(I,se3);se4=strel(disk,5);erodedBw2=imerode(I,se4);se1=strel(ball,3,3);I1=imdilate(a,se1);se2=strel(ball,5,5);I2=imdilate(a,se2);subplot(2,2,1);imshow(I1);title(3*3膨胀图像)sub

25、plot(2,2,2);imshow(I2);title(5*5膨胀图像)subplot(2,2,3);imshow(erodedBw1);title(3*3腐蚀图像)subplot(2,2,4);imshow(erodedBw2);title(5*5腐蚀图像)11/25试验六图像处理实际应用一、试验目标1、熟悉并掌握mATLAB图像处理工具箱的使用；2、了解并掌握常见的图像处理技术。二、试验环境以上版本、wINXP或wIN2021计算机三、试验内容调试运行下列程序,分析程序,对每条语句给出注释,并显示最终实施结果。总结算法思想及优缺点.I=imread();y,x,z=size(I);myI

26、=double(I);%RGBtoHSI%tic%y方向%Blue_y=zeros(y,1);fori=1:yforj=1:xif(myI(i,j,1)=51)&(myI(i,j,3)=119)%蓝色RGB的灰度范围Blue_y(i,1)=Blue_y(i,1)+1;endendendtempmaxy=max(Blue_y);%y方向车牌区域确定Py1=maxy;while(Blue_y(Py1,1)=5)&(Py11)Py1=Py1-1;endPy2=maxy;while(Blue_y(Py2,1)=5)&(Py2=51)&(myI(i,j,3)=119)Blue_x(1,j)=Blue_x

27、(1,j)+1;endendendPX1=1;while(Blue_x(1,PX1)PX1)PX2=PX2-1;endPX1=PX1-2;%对车牌区域的修正PX2=PX2+2;Plate=I(Py1:Py2,PX1-2:PX2,:);t=toc%读取计时%figure,imshow(I);figure,plot(Blue_y);gridfigure,plot(Blue_x);grid框架三层框架三层,局部为二层钢构。本工程外脚手架采取落地式钢管脚手架局部为二层钢构。本工程外脚手架采取落地式钢管脚手架,外架随主体结构上升外架随主体结构上升,同时搭设同时搭设,比操作面高出一步比操作面高出一步,确保

28、主体及外装修的正常安全施工。确保主体及外装修的正常安全施工。figure,imshow(Iy);添加注释和更正后的程序：I=imread(f:);%读取图片y,x,z=size(I);%给定图片大小myI=double(I);%返回双精度值%RGBtoHSI%tic%计时开始%y方向%Blue_y=zeros(y,1);%一列全零矩阵fori=1:y%给定i的范围forj=1:x%给定j的范围if(myI(i,j,1)=51)&(myI(i,j,3)=119)%蓝色RGB的灰度范围Blue_y(i,1)=Blue_y(i,1)+1;%y矩阵加一endendendtempmaxy=max(Blu

29、e_y);%y方向车牌区域确定Py1=maxy;%y矩阵的最大元素while(Blue_y(Py1,1)=5)&(Py11)%确定蓝色RGBBlue_y所在位置Py1=Py1-1;%对车牌区域的修正,向上调整endPy2=maxy;%y矩阵的最大元素while(Blue_y(Py2,1)=5)&(Py2=51)&(myI(i,j,3)=119)%调整I中的范围Blue_x(1,j)=Blue_x(1,j)+1;%Blue_x的矩阵加一endendendPX1=1;%当PX1等于1时while(Blue_x(1,PX1)PX1)%确定蓝色RGBBlue_x的位置PX2=PX2-1;%对车牌区域的

30、修正向下调整endPX1=PX1-2;%对车牌区域的修正PX2=PX2+2;%对车牌区域的修正Plate=I(Py1:Py2,PX1-2:PX2,:);%矩阵行列的范围t=toc%读取计时%figure,imshow(I);%创立图像窗口,显示图像Ifigure,plot(Blue_y);grid%创立图像窗口,绘制Blue_y图像,画出网格线figure,plot(Blue_x);grid%创立图像窗口,绘制Blue_x图像,画出网格线figure,imshow(Plate);%创立图像窗口,显示图像Platet=13/25框架三层框架三层,局部为二层钢构。本工程外脚手架采取落地式钢管脚手架

31、局部为二层钢构。本工程外脚手架采取落地式钢管脚手架,外架随主体结构上升外架随主体结构上升,同时搭设同时搭设,比操作面高出一步比操作面高出一步,确保主体及外装修的正常安全施工。确保主体及外装修的正常安全施工。15/25框架三层框架三层,局部为二层钢构。本工程外脚手架采取落地式钢管脚手架局部为二层钢构。本工程外脚手架采取落地式钢管脚手架,外架随主体结构上升外架随主体结构上升,同时搭设同时搭设,比操作面高出一步比操作面高出一步,确保主体及外装修的正常安全施工。确保主体及外装修的正常安全施工。17/25试验二图像基础操作一、试验目标1、熟悉及掌握图像的采样原理,实现图像的采样过程,进行图像的灰度转换。

32、2、了解直方图的概念及应用,实现图像直方图的显示,及经过直方图均衡和直方图要求化方法对图像进行修正。二、试验环境以上版本、wINXP或wIN2021计算机三、相关函数l图像的增强1、直方图imhist函数用于数字图像的直方图计算或显示,imhist(I,n)计算和显示图像I的直方图,n为指定的灰度级数目,默认为256。假如I是二值图像,那么n仅有两个值。counts,x=imhist(.)返回直方图数据向量counts,对应的色彩值向量x。如：i=imread(e:);imhist(i);2、直方图均衡化histeq函数用于数字图像的直方图均衡化,j=histeq(I,n)均衡化后的级数n,缺

33、省值为64。j=histeq(I,hgram)直方图要求化,马上原是图象I的直方图变换成用户指定的向量hgram(即指定另一幅图像的直方图数据向量)。如：i=imread(e:);j=histeq(i,N);对图像i实施均衡化,得到含有N个灰度级的灰度图像j,N缺省值为643、灰度调整imadjust函数用于数字图像的灰度或颜色调整,j=imadjust(I)将灰度图像I中的亮度值映射到j中的新值并使1的数据是在低高强度和饱和,这增加了输出图像j的对比度值。j=imadjust(I,low_in;high_in,low_out;high_out)将图像I中的亮度值映射到j中的新值,马上low_

34、in至high_in之间的值映射到low_out至high_out之间的值。low_in以下和high_in以上的值被剪切掉了,也就是说,low_in以下的值映射到low_out,high_in以上的值映射到high_out。它们全部能够使用空的矩阵,默认值是01。如：i=imread(e:);j=imadjust(i,);将图像i转换为j,使灰度值从和缺省值01相匹配l运算函数1、Zeros生成全0数组或矩阵如B=zeros(m,n)orB=zeros(mn)返回一个m*n全0矩阵2、取整函数floor最小取整函数框架三层框架三层,局部为二层钢构。本工程外脚手架采取落地式钢管脚手架局部为二层

35、钢构。本工程外脚手架采取落地式钢管脚手架,外架随主体结构上升外架随主体结构上升,同时搭设同时搭设,比操作面高出一步比操作面高出一步,确保主体及外装修的正常安全施工。确保主体及外装修的正常安全施工。round四舍五入取整函数ceil最大取整函数如a=-,-,+I=round(a)I=-203672+4i四、试验内容1、对一幅图像进行2倍、4倍、8倍和16倍减采样,显示结果。a=imread(f:/);b=rgb2gray(a);form=1:4figurewidth,height=size(b);quartimage=zeros(floor(width/(m),floor(height/(2*m

36、);k=1;n=1;fori=1:(m):widthforj=1:(2*m):heightquartimage(k,n)=b(i,j);n=n+1;endk=k+1;n=1;endimshow(uint8(quartimage);End19/252、显示一幅灰度图像a,改变图像亮度使其整体变暗得到图像b,显示两幅图像的直方图。a=imread(f:lena_)c=rgb2gray(a)b=c-46subplot(2,1,1);imhist(c);title(原图像)subplot(2,1,2);imhist(b);title(变暗后的图像)框架三层框架三层,局部为二层钢构。本工程外脚手架采取落

37、地式钢管脚手架局部为二层钢构。本工程外脚手架采取落地式钢管脚手架,外架随主体结构上升外架随主体结构上升,同时搭设同时搭设,比操作面高出一步比操作面高出一步,确保主体及外装修的正常安全施工。确保主体及外装修的正常安全施工。3、对图像b进行直方图均衡化,显示结果图像和对应直方图。b=imread(f:lena_)c=rgb2gray(b)j=histeq(c)subplot(2,2,1),imshow(c)subplot(2,2,2),imshow(j)subplot(2,2,3),imhist(c)%显示原始图像直方图subplot(2,2,4),imhist(j)%显示均衡化后图像的直方图4、

38、读入图像c,实施直方图要求化,使图像a的灰度分布和c大致相同,显示变换后图像及对应直方图。I=imread(f:);j=histeq(I,64);%均衡化成32个灰度级的直方图counts,x=imhist(j);%返回直方图图像向量countsb=imread(f:)Q=rgb2gray(b)m=histeq(Q,counts);%将原始图像Q的直方图变成指定向量countsfigure,21/25subplot(3,2,1),imshow(I);subplot(3,2,2),imshow(Q);subplot(3,2,3),imhist(I);subplot(3,2,4),imhist(Q

39、)subplot(3,2,5),imhist(j);subplot(3,2,6),imhist(m);框架三层框架三层,局部为二层钢构。本工程外脚手架采取落地式钢管脚手架局部为二层钢构。本工程外脚手架采取落地式钢管脚手架,外架随主体结构上升外架随主体结构上升,同时搭设同时搭设,比操作面高出一步比操作面高出一步,确保主体及外装修的正常安全施工。确保主体及外装修的正常安全施工。试验三图像变换一、试验目标1、熟悉及掌握图像的变换原理及性质,实现图像的傅里叶变换。二、试验环境以上版本、wINXP或wIN2021计算机三、相关函数l图像的变换1fft22fft2函数用于数字图像的二维傅立叶变换,如：i=

40、imread(e:);j=fft2(i);因为BmATLAB无法显示复数图像,所以变换后的结果还需进行求模运算,即调用sabs函数。以后经常还进行对数变换,即调用glog函数,以减缓傅里叶谱的快速衰减,更加好地显示高频信息。2ifft22ifft2函数用于数字图像的二维傅立叶反变换,如：i=imread(e:);j=fft2(i);k=ifft2(j);3fftshift用于将变换后图像频谱中心从矩阵的原点移动到矩阵的中心B=fftshift（ii）4利用2fft2计算二维卷积利用2fft2函数能够计算二维卷积,如:a=8,1,6;3,5,7;4,9,2;b=1,1,1;1,1,1;1,1,1

41、;a(8,8)=0;b(8,8)=0;c=ifft2(fft2(a).*fft2(b);c=c(1:5,1:5);利用conv2(二维卷积函数)校验,如：a=8,1,6;3,5,7;4,9,2;b=1,1,1;1,1,1;1,1,1;c=conv2(a,b);四、试验内容1、对一幅图像进行平移,显示原始图像和处理后图像,分别对其进行傅里叶变换,显示变换后结果,分析原图的傅里叶谱和平移后傅里叶频谱的对应关系。s=imread(f:);%读入原图像i=rgb2gray(s)23/25i=double(i)j=fft2(i);%傅里叶变换k=fftshift(j);%直流分量移到频谱中心l=log(

42、abs(k);%对数变换m=fftshift(j);%直流分量移到频谱中心RR=real(m);%取傅里叶变换的实部II=imag(m);%取傅里叶变换的虚部A=sqrt(RR.2+II.2);%计算频谱府幅值A=(A-min(min(A)/(max(max(A)*255;%归一化b=circshift(s,800450);%对图像矩阵im中的数据进行移位操作b=rgb2gray(b)b=double(b)c=fft2(b);%傅里叶变换e=fftshift(c);%直流分量移到频谱中心l=log(abs(e);%对数变换f=fftshift(c);%直流分量移到频谱中心ww=real(f);

43、%取傅里叶变换的实部BZZ=imag(f);%取傅里叶变换的虚部B=sqrt(ww.2+ZZ.2);%计算频谱府幅值B=(B-min(min(B)/(max(max(B)*255;%归一化subplot(2,2,1);imshow(s);title(原图像)subplot(2,2,2);imshow(uint8(b);title(平移图像)subplot(2,2,3);imshow(A);title(离散傅里叶频谱);subplot(2,2,4);imshow(B);title(平移图像离散傅里叶频谱)2、对一幅图像进行旋转,显示原始图像和处理后图像,分别对其进行傅里叶变换,显示变换后结果,分

44、析原图的傅里叶谱和旋转后傅里叶频谱的对应关系。框架三层框架三层,局部为二层钢构。本工程外脚手架采取落地式钢管脚手架局部为二层钢构。本工程外脚手架采取落地式钢管脚手架,外架随主体结构上升外架随主体结构上升,同时搭设同时搭设,比操作面高出一步比操作面高出一步,确保主体及外装修的正常安全施工。确保主体及外装修的正常安全施工。s=imread(f:);%读入原图像i=rgb2gray(s)i=double(i)j=fft2(i);%傅里叶变换k=fftshift(j);%直流分量移到频谱中心l=log(abs(k);%对数变换m=fftshift(j);%直流分量移到频谱中心RR=real(m);%取

45、傅里叶变换的实部II=imag(m);%取傅里叶变换的虚部A=sqrt(RR.2+II.2);%计算频谱府幅值A=(A-min(min(A)/(max(max(A)*255;%归一化b=imrotate(s,-90);%对图像矩阵im中的数据进行移位操作b=rgb2gray(b)b=double(b)c=fft2(b);%傅里叶变换e=fftshift(c);%直流分量移到频谱中心l=log(abs(e);%对数变换f=fftshift(c);%直流分量移到频谱中心ww=real(f);%取傅里叶变换的实部BZZ=imag(f);%取傅里叶变换的虚部B=sqrt(ww.2+ZZ.2);%计算频谱府幅值B=(B-min(min(B)/(max(max(B)*255;%归一化subplot(2,2,1);imshow(s);title(原图像)subplot(2,2,2);imshow(uint8(b);title(平移图像)subplot(2,2,3);imshow(A);title(离散傅里叶频谱);subplot(2,2,4);imshow(B);title(平移图像离散傅里叶频谱)25/25知识改变命运