佳木斯大学数字图像报告.doc

数字图像处理实验报告(一)学院名称信息电子技术专业名称通信年级/班级通信一班姓 名学 号 指导教师实验日期2015-03-21实验地点1C06-329实验成绩实验项目名称实验一：MATLAB编程基础实验目的了解MATLAB语言的编程特点，熟悉MATLAB软件环境，掌握一般的编程方法和技巧。熟悉矩阵和数组的操作及数组索引的常用方法，掌握常用的数据类型及转换的方法，熟练掌握M文件的编写，为利用MATLAB语言进行数字图像处理打下基础。 实验内容1-1）在命令窗口创建矩阵a、b、c，将3个变量存储为数据文件（文件名为:data），并在清除工作空间后调用它，同时查看矩阵a的维数。参考命令：save，load，clear all，size a= b= c= d= 3 4 5 6 -0.5 0.5 1 2 5 6 7 8 9 10 0.75 1.5 3 4 7 81-2）生成空阵，33的全0阵、全1阵、单位阵和随机矩阵；参考命令： ，zeros，ones，eye，rand1-3）利用for循环语句生成一个数组AA= 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 201-4）利用实验1-1）中的矩阵c和d实现下列操作：c+d、c-d、c*d、c.*d，数组c的转置、c求平方、c中所有元素之和，并将c和d合并为一个数组e；1-5）将实验1-1）中的矩阵b和c分别转换为uint8和logical类数组，并将矩阵b归一化为double类数组，另外将uint8型的数组48, 125；100, 248转化为double类数组；参考命令：uint8，double，im2uint8，mat2gray，im2double，im2bw1-6）按以下方式检索查找实验1-3）中矩阵A的相关元素：A(8)、A(1,7,14,20)、A(end-4:end)、A(3, 4)、A(1:3, 2:4)、A(1 2 3, 3 4)、A(1:4, 5)、A(:, 5)、A(end, :)、A(17:20)。1-7）建立一个m函数文件，实现一个两参数的加法和乘法运算，并在命令窗口中调用它；1-8）建立一个m脚本文件，比较向量化循环与for循环的执行速度，参考公式为：f(x, y) = Asin(ux + vy)，其中x、y范围可定为0-99，A、u和v参数值可自行设定；参考命令：tic、toc实验要求五号宋体，单倍行距实验原理五号宋体，单倍行距实验结果及分析1题。 a=3,4,5,6;7,8,9,10a = 3 4 5 6 7 8 9 10 b=-0.5,0.5;0.75,1.5b = -0.5000 0.5000 0.7500 1.5000 c=1,2;3,4c = 1 2 3 4 d=5,6;7,8d = 5 6 7 8 save data a b c d clear all load data a size(a)ans = 2 42题。 e=e = F=zeros(3)F = 0 0 0 0 0 0 0 0 0 G=ones(3)G = 1 1 1 1 1 1 1 1 1 H=eye(3)H = 1 0 0 0 1 0 0 0 1 I=rand(3)I = 0.9501 0.4860 0.4565 0.2311 0.8913 0.0185 0.6068 0.7621 0.82143题 for i=1:4for j=1:5A(i,j)=(i-1)*5+j;endend AA = 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 204题. load data a b c d c+dans = 6 8 10 12 c-dans = -4 -4 -4 -4 c*dans = 19 22 43 50 c.*dans = 5 1221 32 U=rot90(c,1)U = 2 4 1 3 J=c.2J = 1 4 9 16 K=sum (c)K = 4 6 L=c dL = 1 2 5 6 3 4 7 85题 load data a b c d m=uint8(b)m = 0 1 1 2 e=mat2gray(c)e = 0 0.3333 0.6667 1.0000 n=im2bw(e)n = 0 0 1 1 g=mat2gray(b,0,1)g = 0 0.5000 0.7500 1.0000 h=uint8(48,125;100,248)h = 48 125 100 248 l=double(h)l = 48 125 100 2486题 A(8)ans = 17 A(1,7,14,20)ans = 1 12 9 20 A(end-4:end)ans = 19 5 10 15 20 A(3, 4)ans = 14 A(1:3, 2:4)ans = 2 3 4 7 8 912 13 14A(1 2 3, 3 4)ans = 3 4 8 9 13 14 A(1:4, 5)ans = 5 10 15 20 A(:, 5)ans = 5 10 15 20 A(end, :)ans = 16 17 18 19 20 A(17:20)ans = 5 10 15 207题function y=myfunction(a,b)y=a*b+a*b; t=myfunction(1,2)t = 48题x=1:99;y=1:99;a=1,u=2,v=3;ticfor i=1:99 for j=1:99 f(i,j)=a*sin(u*i+v*j); endendt1=tocticf=a*sin(u*x+v*y);t2=toctimes=t1/t2 sosoa = 1u = 2t1 = 0.0110t2 = 0.0140times = 0.7857实验总结五号宋体，单倍行距指导教师意见 签名： 年 月 日填写说明： 表格可根据实际填写内容进行扩展调整； 填写表格字体为宋体五号字，行间距为单倍； 表头处 “数字图像处理实验报告( )”括号中，填写第几次实验，如第一次实验就填写“一”； 实验原理处，主要填写与实验内容相关的理论、方法，以及具体流程和程序代码等内容； 实验结果与分析部分：要标明是那一部分内容的结果，如第一个实验的第二个内容标为1-2，具体标注为:1-2结果如图几所示，接下来就本结果进行适当的分析； 实验总结部分；主要写心得、体会，通过实验解决了什么问题等。- 9 -数字图像处理实验报告(二)学院名称信息电子技术专业名称通信年级/班级姓 名 学 号指导教师实验日期2015-03-21实验地点1C06-329实验成绩实验项目名称实验二：图像基本操作实验目的了解MATLAB数字图像处理工具箱，熟悉数字图像基本处理操作原理，掌握相关工具箱函数调用及其参数设置方法，能够查看图像的相关信息，并分析图像的处理效果。实验内容2-1）、将目标图像加载到MATLAB目录下work文件夹中，读入并显示目标图像I，查看图像的基本信息及其在内存中的情况； 参考函数imread、imshow、imfinfo、whos 2-2）、创建名为exp2-1.m程序文件，在文件中实现以下的图像操作：对图像I进行缩放处理，将其缩小为0.5倍，并在当前空间下存储名称为”exp11”的图片；对图像I进行旋转处理，将其顺时针旋转30，并在当前空间下存储名称为”exp12”的图片；分别给原始图像和变换后图像赋予相应的标题名称（输入图片为”原图像”、输出图像为”缩小0.5倍图像”和”旋转30图像”），并在同一个窗口中，以1行3列的方式同时显示几幅图片进行比较； 参考函数imresize、imrotate、title、subplot、subimage2-3）、通过交互方式将图像中的人眼和人脸区域剪切提取出来，分别存名为”eye”和”face”的新的图像，并查看眼部图片数据； 参考函数imcrop2-4）、将原始图片分别转换成HSV色彩空间下的图像、灰度图像和二值图像，并在同一窗口中，以2行2列的方式同时显示几幅图片进行比较（四幅图像的标题名称分别为”原图像”、”HSV图像”、”灰度图像”、”二值图像”）； 参考函数rgb2hsv、rgb2gray、im2bw、title2-5）、在原始图像中，通过交互方式，任意选取四个像素点，查看像素点的灰度；给出原始图片在R、G、B三个色彩空间下的图像，并在同一窗口中，以2行2列的方式同时显示几幅图片进行比较（四幅图像的标题名称分别为”原图像”、”R图像”、”G图像”、”B图像”）；2-6）、将原始图片转换为灰度图片，并分别在水平和垂直方向添加颜色条，在同一个窗口中，以1行3列的方式同时显示几幅图片进行比较；实验要求1）、实验采用的原始图片要求是包含自己头像的照片，图片大小控制在640480之内；2）、实验中的当前工作目录采用MATLAB目录下的work文件夹。实验原理实验结果及分析1题 A=imread(qi.jpg); imshow(A)2题A=imread(qi.jpg);subplot(1,3,1),imshow(A);title(原图像);B=imresize(A,0.5);subplot(1,3,2),imshow(B);title(缩小0.5倍);C=imrotate(A,30);subplot(1,3,3),imshow(C);title(旋转30度);3题 C=imcrop(B) imwrite(C,eye.jpg) imshow(C)D=imcrop(A)imwrite(D,face.jpg)imshow(D)4题G=rgb2hsv(B);H=rgb2gray(B); I=im2bw(B);subplot(2,2,1);imshow(B),title(原图像); subplot(2,2,2);imshow(G),title(HSV 图像);subplot(2,2,3);imshow(H),title(灰度图像); subplot(2,2,4);imshow(I),title(二值图像); 5题6题I=imread(qi.jpg);subplot(1,3,1);imshow(I)subplot(1,3,2);imshow(I)colorbarsubplot(1,3,3);imshow(I)colorbar(horizontal)实验总结五号宋体，单倍行距指导教师意见 签名： 年 月 日填写说明： 表格可根据实际填写内容进行扩展调整； 填写表格字体为宋体五号字，行间距为单倍； 表头处 “数字图像处理实验报告( )”括号中，填写第几次实验，如第一次实验就填写“一”； 实验原理处，主要填写与实验内容相关的理论、方法，以及具体流程和程序代码等内容； 实验结果与分析部分：要标明是那一部分内容的结果，如第一个实验的第二个内容标为1-2，具体标注为:1-2结果如图几所示，接下来就本结果进行适当的分析； 实验总结部分；主要写心得、体会，通过实验解决了什么问题等。数字图像处理实验报告(三)学院名称信息电子技术专业名称电子信息工程年级/班级姓 名学 号指导教师实验日期2015-03-29实验地点1C06-329实验成绩实验项目名称图像灰度统计特性及其相关变换实验目的理解直方图的形成原理，掌握绘制灰度直方图的方法；熟悉图像灰度直方图的变换及直方图均衡化方法；理解图像灰度变换处理在图像增强中的作用；熟悉图像灰度分布统计与图像视觉质量之间的关系；通过工具箱函数调用和结构化编程两种方式实现图像的相关处理，在加深理解基本原理的同时，提高编程实践的技巧和能力。实验内容3-1读入RGB目标图像，并将其转换成灰度图像，分别赋给变量M和N，查看变换前后图像大小，并分别在变换前后的图像上返回选定图像像素点（116，248），（75，184），（254，310）的灰度值；参考函数imread、size、impixel 3-2根据直方图定义，通过结构化编程方式给出3-1）中灰度图像像素点统计形式的直方图，并与利用函数调用方式获得的直方图在两个不同窗口中进行比较，两窗口图像名称分别为”编程直方图”、”函数调用直方图”；参考函数imread、size、bar、imhist、image3-3根据直方图定义，通过结构化编程方式给出3-1）中灰度图像的概率形式的直方图；3-4根据直方图均衡化原理，通过结构化编程和函数调用两种方式，实现图像的均衡化处理，并比较处理的结果（两种均衡化图像名称分别为”编程均衡化直方图”、”函数调用均衡化直方图”）；参考函数 imread、size、find、unique、histeq3-5利用对比度调整方法来改善图像质量；参考函数 imadjust3-6通过编程方式对3-1）中灰度图像N实施分段灰度变换处理，相关变换公式为： 实验要求1）、实验采用的原始图片要求是真彩色的，图片大小控制在640480之内；2）、针对实验的每一个内容建立相应的M文件，文件名分别为exp31.m exp36.m。实验结果及分析实验总结指导教师意见 签名： 年 月 日- 30 -数字图像处理实验报告(四)学院名称信息电子技术专业名称年级/班级姓 名 学 号12指导教师杜云明实验日期2015-04-5实验地点1C06-329实验成绩实验项目名称图像平滑与锐化实验目的结合课上学习的相关图像增强（平滑和锐化）处理知识，利用编程和相关函数知识实现图像的处理和结果的比较分析。着重加强对模板概念的理解，掌握利用模板通过编程实现图像的增强操作。实验内容4-1）、将自己的照片（彩色），转换为灰度图像，然后分别加入椒盐噪声和高斯噪声，在名为” The simulation window 1”的窗口中同时显示彩图、灰度图和两个加入噪声的图片，分别命名为“原始图像”、 “灰度图像”、“椒盐噪声图像”、“高斯噪声图像”；参考函数 rgb2gray、imnoise、figure4-2）、利用33的BOX模板，通过编程方式，实现对椒盐噪声图像的平滑去噪处理（边界用原图像素），并在名为” The simulation window 2”的窗口中同时显示原图与处理后图像；4-3）、分别利用33、55和77的BOX模板，通过函数调用的方式实现对椒盐噪声图像的平滑去噪处理，并在名为” The simulation window 3”的窗口中同时显示原图与处理后图像，分析不同模板对图像的处理效果；参考函数 imfilter、figure4-4）、采用33窗口，通过编程和函数调用两种方式，实现对含有椒盐噪声的图像进行中值滤波处理，并在名为” The simulation window 4”的窗口中同时显示原图与处理后的图片，名称分别为“原图”、“编程中值滤波图像”和“函数中值滤波图像”；参考函数 medfilt2 4-5）、读取图片“rice.png”设为变量I，利用水平-垂直差分法，通过结构化编程方式实现图片锐化处理（即获得I的梯度幅度图），并在同一窗口显示原图和其对应的梯度图，两图分别命名为“原始图像”和“梯度图像”；参考函数 imfilter 4-6）、读取图片“rice.png”设为变量J，分别用梯度算子“roberts”、“prewitt”和“sobel”，通过函数调用编程方式实现图像的锐化处理，在同一窗口显示原图与处理后图像，比较分析三种算子的锐化差异；参考函数 fspecial、imfilter4-7）、利用8方向的拉普拉斯算子模板，通过结构化编程方式实现对“blood.bmp”图像的锐化处理，比较分析其处理结果与梯度算子处理结果的差异。参考函数 imfilter 实验要求1）、要求在给图像加噪时，噪声的强度设定为0.02；2）、要求实验4-1) 4-4)用一个名为exp41.m的文件来实现；3）、要求在同一个窗口显示时，图片数量3时，以1行的方式进行显示，图片数量3时，以2行的方式进行显示；4）、在关键语句上添加必要的注释，以辅助说明该语句的功能。实验结果及分析1题I=imread(qi.jpg);I1=rgb2gray(I);I2=imnoise(I1,salt pepper,0.02);I3=imnoise(I1,gaussian,0.02);figure(Name,The simulation window 1),subplot(2,2,1),imshow(I),title(原始图像)subplot(2,2,2),imshow(I1),title(灰度图像)subplot(2,2,3),imshow(I2),title(椒盐噪声图像)subplot(2,2,4),imshow(I3),title(高斯噪声图像)2题A=imread(qi.jpg);B=rgb2gray(A);I=imnoise(B,salt pepper,0.02);I=double(I);s=1/9*ones(3);m n=size(B);C=zeros(m,n);for i=2:m-1 for j=2:n-1 C(i,j)=sum(sum(I(i-1:i+1,j-1:j+1).*s); C(1,:)=I(1,:);C(m,:)=I(m,:);C(:,1)=I(:,1);C(:,n)=I(:,n); endendfigure(name,The simulation window 2),subplot(1,2,1);imshow(I),title(原始图像)subplot(1,2,2);imshow(uint8(C),title(处理后图像)3题I=imread(qi.jpg);I1=imnoise(I,salt pepper,0.02);s=1/9*ones(3);s1=1/25*ones(5);s2=1/49*ones(7);A=imfilter(I1,s);A1=imfilter(I1,s1);A2=imfilter(I1,s2);figure(name,The simulation window 3),subplot(2,2,1);imshow(I1),title(原图)subplot(2,2,2);imshow(A),title(3处理后图像)subplot(2,2,3);imshow(A1),title(5处理后图像)subplot(2,2,4);imshow(A2),title(7处理后图像) 实验总结指导教师意见 签名： 年 月 日数字图像处理实验报告(五)学院名称信息电子技术专业名称电子信息工程年级/班级姓 名学 号1指导教师杜云明实验日期2015-04-12实验地点1C06-329实验成绩实验项目名称图像分割与边缘检测实验目的通过结构化编程研究和验证图像分割与边缘检测的常用算法原理；掌握MATLAB图像域值分割与边缘检测函数的使用方法；了解边缘检测的算法和用途，比较Sobel、Prewitt、Canny等算子边缘检测的差异。实验内容5-1）、自行设定阈值，利用编程方式实现“rice.png”图像的二值化处理；A=imread(rice.png);newA=im2bw(A,150/255);imshow(newA);5-2）、利用阈值分割技术提取“rice.png”图像中的目标，要求在同一窗口显示原图、原图直方图、T=120分割图和系统自定阈值的分割图，分别命名为“原始图像”、“原始图像直方图”、“T=120阈值分割图像”和“系统自定阈值分割图像”；参考函数 im2bw、graythreshA=imread(rice.png);level = graythresh(A);BW = im2bw(A,level);newA=im2bw(A,120/255);subplot(2,2,1);imshow(A);title(yuan);subplot(2,2,2);imshow(newA);title(fg);subplot(2,2,3);imhist(A);title(zhifang);subplot(2,2,4);imshow(BW);title(xitongzid);5-3）、对“blood.bmp”图像进行阈值分割提取目标，在同一窗口给出原图，同时比较T=91、T=130和T=43的分割效果，实验图像分别命名为“原始图像”、 “T=91阈值分割图像”、“T=130阈值分割图像”和“T=43阈值分割图像”，分析阈值对图像分割的影响；参考函数 im2bw （注意函数中相关参数设置要求在0，1区间之内）A=imread(blood.bmp);newj=im2bw(A,91/255);newy=im2bw(A,130/255);news=im2bw(A,43/255);subplot(2,2,1);imshow(A);title(yuan);subplot(2,2,2);imshow(newj);title(91);subplot(2,2,3);imshow(newy);title(130);subplot(2,2,4);imshow(news);title(43);5-4）、使用edge函数对图像“rice.png”进行边缘检测，比较roberts，sobel，prewitt，canny，log算子的检测效果；I=imread(rice.png);bw1=edge(I,roberts);bw2=edge(I,sobel);bw3=edge(I,log);bw4=edge(I,prewitt);bw5=edge(I,canny);subplot(2,3,1);imshow(I),title(原图)subplot(2,3,2);imshow(bw1),title(roberts)subplot(2,3,3);imshow(bw2),title(sobel)subplot(2,3,4);imshow(bw3),title(log)subplot(2,3,5);imshow(bw4),title(prewitt)subplot(2,3,6);imshow(bw5),title(canny) 5-5）、利用8连通准则，检测图像I=1 1 1 0 0 0 0 0；1 1 1 0 1 0 0 0；1 1 1 0 1 1 0 0；1 1 1 0 0 0 1 0；1 1 1 1 0 0 1 0；1 1 1 1 0 0 1 0；1 1 1 0 0 1 1 0；1 1 1 0 1 0 0 0中的目标的个数；参考函数 bwlabel5-6）、利用结构化编程方式提取图像“bt.bmp”的边缘，要求在同一窗口内显示原图像和轮廓图像，实验图像分别命名为“原始图像”和“边缘提取图像”；5-7）、利用结构化编程方式跟踪图像0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5；0.5 0.5 1 0.5 1 0.5 0.5 0.5；0.5 1 0 1 0 1 1 0.5；0.5 1 0 0 0 0 1 0.5；0.5 1 0 1 1 0 1 0.5；0.5 0.5 1 0.5 0.5 1 0.5 0.5；0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5的轮廓，并给出其边界连码；（注：先将图像转换为灰度图像然后再进行处理，要求跟踪其8方向的轮廓）5-8）、用相机在同一场景内拍摄两幅图片，其中一幅作为背景，另一幅含有本人的图片为待处理图片，利用两幅图片研究图像差影法的处理效果。实验要求1）、上述实验要以给定图像为目标进行展开；2）、所有实验要用独立的M文件来实现，文件分别为exp51.mexp58.m；3）、同一个窗口显示时，图片数量3时，以1行的方式进行显示，图片数量3时，以2行的方式进行显示；4）、在关键语句和重要步骤上添加必要的注释，以辅助说明该语句和段落的功能。实验结果及分析实验总结指导教师意见 签名： 年 月 日数字图像处理实验报告(六)学院名称信息电子技术专业名称电子信息工程年级/班级姓 名学 号指导教师杜云明实验日期2015-04-19实验地点1C06-329实验成绩实验项目名称图像几何变换实验目的了解几何变换的概念及其实现方式；掌握利用相关矩阵运算实现图像的几何变换；熟悉多种方式实现图像几何变换处理的编程技巧。实验内容6-1）、以结构化编程方式，实现利用矩阵线性运算对“lina.bmp”图像的水平、垂直和对角镜像的几何变换处理，在一个窗口下显示原图像、水平镜像图、垂直镜像图和对角镜像图；I=imread(lina.bmp);p,q=size(I);g=zeros(p,q);f=zeros(p,q);l=zeros(p,q);fh=p;fw=q;for i=1:p for j=1:q a=1 0 0;0 -1 fw;0 0 1*i j 1; b=-1 0 fh;0 1 0;0 0 1*i j 1; c=-1 0 fh;0 -1 fw;0 0 1*i j 1; g(a(1),a(2)+1)=I(i,j); f(b(1)+1,b(2)=I(i,j); l(c(1)+1,c(2)+1)=I(i,j); endendsubplot(2,2,1);imshow(I),title(原图像)subplot(2,2,2);imshow(uint8(g),title(水平镜像图)subplot(2,2,3);imshow(uint8(f),title(垂直镜像图)subplot(2,2,4);imshow(uint8(l),title(对角镜像图)6-2）、分别用函数调用和结构化编程方式，实现对“lina.bmp”图像的缩放处理，要求将原图缩小为原来的四分之一，在一个窗口下显示原图像、结构化编程缩小图和函数编程缩小图；参考函数 imresizeI=imread(lina.bmp);p,q=size(I);g=zeros(p,q);fh=p;fw=q;for i=1:p/2 for j=1:q/2 a=2 0 0;0 2 0;0 0 1*i j 1; g(i,j)=I(a(1),a(2); endendI1=imresize(I,0.25);subplot(1,3,1);imshow(I),title(原图像)subplot(1,3,2);imshow(uint8(g),title(结构化编程缩小图)subplot(1,3,3);imshow(I1),title(函数编程缩小图)6-3）、用函数调用编程方式，实现对“lina.bmp”图像的旋转处理，要求将原图逆时针旋转25度，在一个窗口下显示原图、旋转图片和旋转后截取效果图；参考函数 imrotate6-4）、利用矩阵运算方式实现对“lina.bmp”图像的平移处理，要求图片在两个方向上分别平移50个像素点，在一个窗口下显示原图、平移图片和平移后截取效果图；6-5）、用函数调用编程方式，实现对图像“lina.bmp”复合变换处理，要求先对图像放大两倍，然后旋转30度，再在两个方向上平移20个像素点，在一个窗口下显示原图和复合变换效果图片；实验要求1）、上述实验要以给定图像为目标进行展开；2）、所有实验要用独立的M文件来实现，文件分别为exp61.mexp65.m；3）、同一个窗口显示时，图片数量3时，以1行的方式进行显示，图片数量3时，以2行的方式进行显示；4）、在关键语句和重要步骤上添加必要的注释，以辅助说明该语句和段落的功能。实验结果及分析实验总结指导教师意见 签名： 年 月 日数字图像处理实验报告(七)学院名称信息电子技术专业名称电子信息工程年级/班级姓 名 学 号指导教师杜云明实验日期2015-04-26实验地点1C06-329实验成绩实验项目名称频域变换实验目的了解图像频域分析的目的、意义和方法；掌握傅立叶变换和离散余弦变换在图像处理中的应用；熟悉图像中各种频率分量的分布；掌握频域处理的编程技巧。实验内容7-1）、编程实现对“rec1.bmp”图像的频域处理，在一个窗口下显示原图像、频谱图像、中心平移后的频谱图像和中心平移后的增强型频谱图像；参考函数 fft2、fftshift、abs、logA=imread(rec1.bmp);N=100;F=fft2(A,N,N);F2=fftshift(abs(F);F3=log(abs(F);subplot(2,2,1);imshow(A);title(yuan);subplot(2,2,2);imshow(F);title(FFT);subplot(2,2,3);imshow(F2);title(pinyi);subplot(2,2,4);imshow(F3);title(LOG);7-2）、编程实现对“rec1.bmp”和“rec2.bmp”图像的频域处理，在一个窗口下显示两幅原图像及其对应的中心平移后的增强型频谱图像，分析研究傅立叶变换的旋转不变性；参考函数 fft2、fftshift、abs、logN=100;A=imread(rec1.bmp);B=imread(rec2.bmp);F=fft2(A,N,N);F2=fftshift(abs(F);E=fft2(B,N,N);E2=fftshift(abs(F);subplot(2,2,1);imshow(A);title(Ayuan);subplot(2,2,2);imshow(F2);title(Abian);subplot(2,2,3);imshow(B);title(Byuan);subplot(2,2,4);imshow(E2);title(Bbian);7-3）、编程实现对“lina.bmp”图像的频域处理，然后在其中心平移后的频谱图像中截取出中心120120区域，进行反变换，在一个窗口下显示原图、中心平移频谱图、截取区域图和其反变换效果图；参考函数 fft2、fftshift、abs、log、ifft2A=imread(lina.bmp);N=256;a=double(A);Z=fft2(a);Z1=fftshift(Z);B=Z1(68:188,68:188);F=ifft2(B,N,N);F1=uint8(abs(F);subplot(2,2,1);imshow(A),title(原图)subplot(2,2,2);imshow(Z1),title(中心平移频谱图)subplot(2,2,3);imshow(B),title(截取区域图)subplot(2,2,4);imshow(F1,),title(反变换效果图)7-4）、编程实现对“lina.bmp”图像的频域处理，然后将其中心平移后的频谱图像中心120120区域设为0，进行反变换，在一个窗口下显示原图、中心平移频谱图、去低频频谱图和其反变换效果图；参考函数 fft2、fftshift、abs、log、ifft2A=imread(lina.bmp);N=256;Z=fft2(A);Z1=fftshift(Z);Z3=fftshift(Z);Z3(68:188,68:188)=0;F=ifft2(Z1,N,N);F1=uint8(abs(F);subplot(2,2,1);imshow(A);title(yuan);subplot(2,2,2);imshow(Z1);title(pinyi);subplot(2,2,3);imshow(Z3);title(qudipin);subplot(2,2,4);imshow(F1);title(fbh);7-5）、以图像“text.png”为例，编程实现字符“a”在图片中的定位，分别在不同窗口下显示原图、字符a、及定位后的效果图效果图；参考函数 fft2、ifft2、rot90、 （提示：字符a在图片中的(32:45,88:98)区域，可将其提取出作为模板，然后通过两图的卷积实现定位；显示定位效果时要设定相关阈值，即保留卷积后的峰值部分）bw=imread(text.png);a=bw(32:45,88:98);B=real(ifft2(fft2(bw).*fft2(rot90(a,2),256,256);subplot(2,2,1);imshow(bw);subplot(2,2,2);imshow(a);subplot(2,2,3);imshow(B,);thresh=60;subplot(2,2,4);imshow(Bthresh);7-6）、编程实现对“autumn.tif”图像的离散余弦变换，并将其中小于10的系数设为0，然后进行反变换，在一个窗口下显示原图、其灰度图、灰度图的余弦变换图和反变换效果图。参考函数 rgb2gray、dct2、idct2A=imread(autumn.tif);H=rgb2gray(A);C=dct2(H);C(abs(C)10)=0;D=idct2(C);subplot(2,2,1);imshow(A);title(yuan);subplot(2,2,2);imshow(H);title(hd);subplot(2,2,3);imshow(C);title(cos);subplot(2,2,4);imshow(D,0 255);title(arccos);实验要求1）、上述实验要以给定图像为目标进行展开；2）、所有程序采用函数调用方式实现，且程序用独立的M文件来保存，文件分别为exp71.mexp76.m；3）、同一个窗口显示时，图片数量3时，以1行的方式进行显示，图片数