精选文档图像处理实验报告

中原工学院实验报告专用纸第页实验一MATLAB数字图像处理初步实验目的与要求1．熟悉及掌握在MATLAB中能够处理哪些格式图像。2．熟练掌握在MATLAB中如何读取图像。3．掌握如何利用MATLAB来获取图像的大小、颜色、高度、宽度等等相关信息。4．掌握如何在MATLAB中按照指定要求存储一幅图像的方法。5.图像间如何转化。实验原理及要求一幅图像可以被定义为一个二维函数f(x,y),其中x和y是空间(平面)坐标，f在任何坐标处(x,y)处的振幅称为图像在该点的亮度。灰度是用来表示黑白图像亮度的一个术语，而彩色图像是由单个二维图像组合形成的。例如，在RGB彩色系统中，一幅彩色图像是由三幅独立的分量图像(红、绿、蓝)组成的。因此，许多为黑白图像处理开发的技术适用于彩色图像处理，方法是分别处理三副独立的分量图像即可。图像关于x和y坐标以及振幅连续。要将这样的一幅图像转化为数字形式，就要求数字化坐标和振幅。将坐标值数字化成为取样；将振幅数字化成为量化。采样和量化的过程如图1所示。因此，当f的x、y分量和振幅都是有限且离散的量时，称该图像为数字图像。作为MATLAB基本数据类型的数值数组本身十分适于表达图像，矩阵的元素和图像的像素之间有着十分自然的对应关系。——————————————————————————————————————————————————————班级：学号：姓名：同组姓名：指导老师：图1图像的采样和量化根据图像数据矩阵解释方法的不同，MATLAB把其处理为4类：亮度图像(Intensityimages)二值图像(Binaryimages)索引图像(Indexedimages)RGB图像(RGBimages)(1)亮度图像一幅亮度图像是一个数据矩阵，其归一化的取值表示亮度。若亮度图像的像素都是uint8类或uint16类，则它们的整数值范围分别是[0，255]和[0，65536]。若图像是double类，则像素取值就是浮点数。规定双精度型归一化亮度图像的取值范围是[0，1](2)二值图像一幅二值图像是一个取值只有0和1的逻辑数组。而一幅取值只包含0和1的uint8类数组，在MATLAB中并不认为是二值图像。使用logical函数可以把数值数组转化为二值数组或逻辑数组。创建一个逻辑图像，其语法为：——————————————————————————————————————————————————————班级：学号：姓名：同组姓名：指导老师：B=logical(A)其中，B是由0和1构成的数值数组。要测试一个数组是否为逻辑数组，可以使用函数:islogical(c)若C是逻辑数组，则该函数返回1；否则，返回0。(3)索引图像索引颜色通常也称为映射颜色，在这种模式下，颜色都是预先定义的，并且可供选用的一组颜色也很有限，索引颜色的图像最多只能显示256种颜色。一幅索引颜色图像在图像文件里定义，当打开该文件时，构成该图像具体颜色的索引值就被读入程序里，然后根据索引值找到最终的颜色。(4)RGB图像一幅RGB图像就是彩色像素的一个M×N×3数组，其中每一个彩色相似点都是在特定空间位置的彩色图像相对应的红、绿、蓝三个分量。按照惯例，形成一幅RGB彩色图像的三个图像常称为红、绿或蓝分量图像。令fR，fG和fB分别代表三种RGB分量图像。一幅RGB图像就利用cat(级联)操作将这些分量图像组合成彩色图像：rgb_image=cat(3,fR,fG,fB)在操作中，图像按顺序放置。实验程序及结果——————————————————————————————————————————————————————班级：学号：姓名：同组姓名：指导老师：Hiyanbaoigao——————————————————————————————————————————————————————班级：学号：姓名：同组姓名：指导老师：J=imread(——————————————————————————————————————————————————————班级：学号：姓名：同组姓名：指导老师：I=imresize(J,0.5);figure,imshow(I)imfinfo'P1.tif';imwrite(I,'P1Jpg.jpg','quality',50);imwrite(I,'P1Bmp.bmp');figure,imshow('P1Bmp.bmp');I_g=rgb2gray(I);gg=im2bw(I_g);figure,imshow(gg)——————————————————————————————————————————————————————班级：学号：姓名：同组姓名：指导老师：————————————————————————————————————————————————————————————————————————————————————————————————————————————班级：学号：姓名：同组姓名：指导老师：实验二图像增强—灰度变换及直方图变换一、实验目的：1．了解图像增强的目的及意义，加深对图像增强的感性认识，巩固所学理论知识。2．学会对图像直方图的分析，掌握灰度直方图的概念及其计算方法。3．掌握直接灰度变换的图像增强方法。4．熟练掌握直力图均衡化和直方图规定化的计算过程；5．掌握色彩直方图的概念和计算方法6．利用MATLAB程序进行图像增强。——————————————————————————————————————————————————————班级：学号：姓名：同组姓名：指导老师：——————————————————————————————————————————————————————班级：学号：姓名：同组姓名：指导老师：2.1灰度变换术语‘空间域’指的是图像平面本身，在空间与内处理图像的方法是直接对图像的像素进行处理。空间域处理方法分为两种：灰度级变换、空间滤波。空间域技术直接对像素进行操作其表达式为g(x,y)=T[f(x,y)]其中f(x,y)为输入图像，g(x,y)为输出图像，T是对图像f进行处理的操作符，定义在点(x,y)的指定领域内。定义点(x,y)的空间邻近区域的主要方法是，使用中心位于(x,y)的正方形或长方形区域，。此区域的中心从原点(如左上角)开始逐像素点移动，在移动的同时，该区域会包含不同的领域。T应用于每个位置(x,y)，以便在该位置得到输出图像g。在计算(x,y)处的g值时，只使用该领域的像素。灰度变换T的最简单形式是使用领域大小为1×1，此时,(x,y)处的g值仅由f在该点处的亮度决定，T也变为一个亮度或灰度级变化函数。当处理单设(灰度)图像时，这两个术语可以互换。由于亮度变换函数仅取决于亮度的值，而与(x,y)无关，所以亮度函数通常可写做如下所示的简单形式：s=T(r)其中，r表示图像f中相应点(x,y)的亮度，s表示图像g中相应点(x,y)的亮度。2.2直方图变换直方图是多种空间城处理技术的基础。直方图操作能有效地用于图像增强。除了提供有用的图像统计资料外，直方图固有的信息在其他图像处理应用中也是非常有用的，如图像压缩与分割。直方图在软件中易于计算，也适用于商用硬件设备，因此，它们成为了实时图像处理的一个流行工具。直方图是图像的最基本的统计特征，它反映的是图像的灰度值的分布情况。直方图均衡化的目的是使图像在整个灰度值动态变化范围内的分布均匀化，改善图像的亮度分布状态，增强图像——————————————————————————————————————————————————————班级：学号：姓名：同组姓名：指导老师：——————————————————————————————————————————————————————班级：学号：姓名：同组姓名：指导老师：图像的直方图事实上就是图像的亮度分布的概率密度函数，是一幅图像的所有象素集合的最基本的统计规律。直方图反映了图像的明暗分布规律，可以通过图像变换进行直方图调整，获得较好的视觉效果。直方图均衡化是通过灰度变换将一幅图像转换为另一幅具有均衡直方图，即在每个灰度级上都具有相同的象素点数的过程。实验程序及结果灰度变化f=imread('P2.jpg');i=rgb2gray(f);subplot(1,2,1);imshow(i,[]);g=imhist(i,256);subplot(1,2,2);plot(g);g1=imadjust(f,[01],[10]);figure,imshow(g1)g2=imadjust(i,[0.50.75],[01]);figure,imshow(g2)g=imread('P2.jpg');————————————————————————————————————————————————————————————————————————————————————————————————————————————班级：学号：姓名：同组姓名：指导老师：h=mat2gray(h);h=im2uint8(h);figure,imshow(h)灰度图像 8位图直方图变化A=imread('P2.jpg');P=rgb2gray(A)[m,n]=size(P);B=zeros(size(P));l=m*n;r=zeros(1,256);y=zeros(1,256);————————————————————————————————————————————————————————————————————————————————————————————————————————————班级：学号：姓名：同组姓名：指导老师：fori=1:mforj=1:nr(A(i,j)+1)=r(A(i,j)+1)+1;endendr=r./l;fork=1:256forj=1:ky(k)=y(k)+r(j);endendfork=1:256y(k)=floor(255*y(k)+0.5);endfori=1:mforj=1:nB(i,j)=y(A(i,j)+1);endendA=uint8(A);B=uint8(B);subplot(1,2,1);imshow(A);title('原图像')subplot(1,2,2);imhist(A);title('原图像的直方图')——————————————————————————————————————————————————————班级：学号：姓名：同组姓名：指导老师：直方图变换实验三图像增强—空域滤波及频域滤波实验目的1.进一步了解MatLab软件/语言，学会使用MatLab对图像作滤波处理，使学生有机会掌握滤波算法，体会滤波效果；2.掌握常用的空域滤波方法；3．掌握怎样利用傅立叶变换进行频域滤波；4．掌握频域滤波的概念及方法；5．熟练掌握频域空间的各类滤波器；——————————————————————————————————————————————————————班级：学号：姓名：同组姓名：指导老师：——————————————————————————————————————————————————————班级：学号：姓名：同组姓名：指导老师：实验原理及要求1.学生应当完成对于给定图像+噪声，使用平均滤波器、中值滤波器对不同强度的高斯噪声和椒盐噪声，进行滤波处理；能够正确地评价处理的结果；能够从理论上作出合理的解释。2.频域滤波分为低通滤波和高通滤波两类，对应的滤波器分别为低通滤波器和高通滤波器。频域低通过滤的基本思想： G(u,v)=F(u,v)H(u,v)F(u,v)是需要钝化图像的傅立叶变换形式，H(u,v)是选取的一个低通过滤器变换函数，G(u,v)是通过H(u,v)减少F(u,v)的高频部分来得到的结果，运用傅立叶逆变换得到钝化后的图像。理想地通滤波器(ILPF)具有传递函数：其中，为指定的非负数，为(u,v)到滤波器的中心的距离。的点的轨迹为一个圆。n阶巴特沃兹低通滤波器(BLPF)(在距离原点处出现截至频率)的传递函数为与理想地通滤波器不同的是，巴特沃兹率通滤波器的传递函数并不是在处突然不连续。——————————————————————————————————————————————————————班级：学号：姓名：同组姓名：指导老师：——————————————————————————————————————————————————————班级：学号：姓名：同组姓名：指导老师：其中，为标准差。相应的高通滤波器也包括：理想高通滤波器、n阶巴特沃兹高通滤波器、高斯高通滤波器。给定一个低通滤波器的传递函数，通过使用如下的简单关系，可以获得相应高通滤波器的传递函数：3.利用MATLAB软件实现空域及频域滤波的程序：实验程序及结果I=rgb2gray(imread('P3.tif'));J=imnoise(I,'gauss',0.02);J=imnoise(I,'salt&pepper',0.02);ave1=fspecial('average',3);ave2=fspecial('average',5);K=filter2(ave1,J)/255;L=filter2(ave2,J)/255;M=medfilt2(J,[33]);N=medfilt2(J,[44]);imshow(I);figure,imshow(J);figure,imshow(K);figure,imshow(L);figure,imshow(M);figure,imshow(N);实验四图像分割实验目的使用MatLab软件进行图像的分割。使学生通过实验体会一些主要的分割算子对图像处理的效果，以及各种因素对分割效果的影响。实验原理及要求要求学生能够自行评价各主要算子在无噪声条件下和噪声条件下的分割性能。能够掌握分割条件(阈值等)的选择。完成规定图像的处理并要求正确评价处理结果，能够从理论上作出合理的解释。(1)使用Roberts算子的图像分割实验调入并显示图像room.tif中图像；使用Roberts算子对图像进行边缘检测处理；Roberts算子为一对模板：相应的矩阵为：rh=[01；-10]；rv=[10；0-1]；这里的rh为水平Roberts算子，rv为垂直Roberts算子。分别显示处——————————————————————————————————————————————————————班级：学号：姓名：同组姓名：指导老师：理后的水平边界和垂直边界检测结果；用“欧几里德距离”和“街区距离”——————————————————————————————————————————————————————班级：学号：姓名：同组姓名：指导老师：提示：先做检测结果的直方图，参考直方图中灰度的分布尝试确定阈值；应反复调节阈值的大小，直至二值化的效果最为满意为止。分别显示处理后的水平边界和垂直边界检测结果；将处理结果转化为“白底黑线条”的方式；给图像加上零均值的高斯噪声；对于噪声图像重复步骤b~f。(2)使用Prewitt算子的图像分割实验使用Prewitt算子进行内容(1)中的全部步骤。(3)使用Sobel算子的图像分割实验使用Sobel算子进行内容(1)中的全部步骤。(4)使用LoG(拉普拉斯-高斯)算子的图像分割实验使用LoG(拉普拉斯-高斯)算子进行内容(1)中的全部步骤。提示1：处理后可以直接显示处理结果，无须另外计算梯度的模。提示2：注意调节噪声的强度以及LoG(拉普拉斯-高斯)算子的参数，观察处理结果。实验程序及结果f=imread('P4.png');f=rgb2gray(f)[gv,t1]=edge(f,'sobel','vertical');imshow(gv)[gb,t2]=edge(f,'sobel','horizontal');figure,imshow(gb)w45=[-2-10;-101;012];g45=imfilter(double(f),w45,'replicate');T=0.3*max(abs(g45(:)));g45=g45>=T;figure,imshow(g45);实验五图像变换实验目的1．了解图像变换的意义和手段；2．熟悉傅立叶变换、小波变换的基本性质；3．熟练掌握FFT、小波变换方法及应用；4．通过实验了解二维频谱的分布特点；5．掌握利用MATLAB编程实现数字图像的傅立叶变换、小波变换。6．评价人眼对图像幅频特性和相频特性的敏感度实验原理及要求2.1傅立叶变换1应用傅立叶变换进行图像处理傅里叶变换是线性系统分析的一个有力工具，它能够定量地分析诸如数字化系统、采样点、电子放大器、卷积滤波器、噪音和显示点——————————————————————————————————————————————————————班级：学号：姓名：同组姓名：指导老师：——————————————————————————————————————————————————————班级：学号：姓名：同组姓名：指导老师：2傅立叶（Fourier）变换的定义对于二维信号，二维Fourier变换定义为：逆变换：二维离散傅立叶变换为：逆变换：图像的傅立叶变换与一维信号的傅立叶变换变换一样，有快速算法，具体参见参考书目，有关傅立叶变换的快速算法的程序不难找到。————————————————————————————————————————————————————————————————————————————————————————————————————————————班级：学号：姓名：同组姓名：指导老师：2.2小波变换“小波”（wavelet)就是一种“尺度”很小的波动，并具有时间和频率特性。小波函数必须满足以下两个条件的函数：(1)小波必须是振荡的；(2)小波的振幅只能在一个很短的一段区间上非零，即是局部化的。2D尺度函数和2D小波函数可由1D尺度函数和小波函数张量积得到，即：图像的2D小波变换包括沿行向(水平方向)和列向(垂直方向)滤波和对角线方向采样。变换后图像将得到近似系数cA1,水平细节cH1,垂直细节cV1,对角细节cD1，最后将各个图像拼接显示在同一图中。2D小波逆变换是将图像的近似系数cA1,水平细节cH1,垂直细节——————————————————————————————————————————————————————班级：学号：姓名：同组姓名：指导老师：——————————————————————————————————————————————————————班级：学号：姓名：同组姓名：指导老师：实验程序及结果傅里叶变换I=imread('P5.gif');imshow(I);fftI=fft2(I);sfftI=fftshift(fftI);RR=real(sfftI);II=imag(sfftI);A=sqrt(RR.^2+II.^2);A=(A-min(min(A)))/(max(max(A))-min(min(A)))*225figure;imshow(A);原图 原图的频谱小波变换oimg=imread('P5.jpg');[cA,cH,cV,cD]=dwt2(oimg,'db1');X=size(oimg,1);cod_cA=wcodemat(cA,X);cod_cH=wcodemat(cH,X);cod_cV=wcodemat(cV,X);cod_cD=wcodemat(cD,X);new_img=[cod_cA,cod_cH;cod_cV,cod_cD];figure,imshow(oimg,[]);figure,imshow(new_img,[]);new_img2=idwt2(cA,cH,cV,cD,'db1');diff=double(oimg)-double(new_img2);figure,subplot(1,3,1);imshow(oimg,[]);subplot(1,3,2);imshow(new_img2,[]);subplot(1,3,3);imshow(diff,[]);——————————————————————————————————————————————————————班级：学号：姓名：同组姓名：指导老师：——————————————————————————————————————————————————————班级：学号：姓名：同组姓名：指导老师：——————————————————————————————————————————————————————班级：学号：