MATLAB数字图像处理基本操作及空域滤波

1、实验一数字图像处理基本操作及空域滤波实验目的1了解有关数字图像处理的基本概念，熟悉软件中关于数字图像处理的基本命令，掌握利用软件进行数字图像处理的简单方法。2、了解并掌握直方图统计方法以及分段线性拉伸、直方图均衡等亮度调整算法，通过观察对这些运算建立感性认识。3、掌握空域滤波中常用的平滑和锐化滤波器。实验内容观察各类图像的直方图；操作灰度对照表，进行分段线性拉伸；采用直方图均衡方法对低对比度的图像进行对比度增强。2掌握图像模板卷积运算的实质，认识各种模板的处理效果；掌握邻域平均及中值滤波降噪方法的特点，认识其功能及适用场合；掌握拉普拉斯算子、Sobel算子、Prewitt算子和Isotropi

2、c算子的特点，认识其功能及适用场合。实验原理1数字图像以一定的格式存放在计算机的存储器中（如磁盘），常见的格式有，等等，要进行数字图像处理，第一项工作就是把图像读到计算机的内存中，以便进行进一步的处理。在中，函数完成此项工作。下面一小段语句即可实现将图像“”显示在一个图像窗口的左边，将其轮廓图显示在该窗口的右边。这里，将磁盘上的图像文件读入内存变量中，生成一个可以横向放置幅图像的窗口，并设置下一显示位置在左边，显示图像，准备下一图像的显示位置，生成图像的轮廓并显示在窗口的右边。数字图像一般可分为二值图、灰度图和真彩图等几类。这里函数完成对二值图像（）的求补运算。（即原来黑的变白的，原来白的变黑

3、的。）该函数也可以用于灰度图像，此时图像的数据被减。在中，函数可用来改变图像的大小，其中表示图像，实数是放大倍数，用来选择计算数据的方法，其可取值及意义分别为：“”一最近邻插值法，“”一双线性插值法，“”一双三次插值法。下面这段程序调整图像的大小：也可以指定目标图像的大小，此时的调用格式为:其中表示将图像的大小调整为8函数对图像进行旋转操作，我们看如下的程序段：注意右边的图像是由参数“”实现的，该参数表示将旋转后的图像取与原图像相同的中心部分输出。2直方图的定义图象的灰度直方图是一个函数，表示数字图象中每一灰度级与该灰度级出现的频数（即具有这一灰度级的象素数目）间的对应关系：实验一数字图像处理

4、基本操作及空域滤波P(b)罟为一幅图象所包含的象素总数；为图象中灰度值为的象素总数。通常，以灰度值为横坐标，为纵坐标。直方图是图象中象素灰度值的一阶概率分布密度的一种近似。、对比度增强对比度增强又称为点运算，逐点改变输入图象的每一象素的灰度，而各象素的位置不改变，一般用来拓宽图象的灰度范围。灰度变换法（对照）典型的对比度拉伸灰度变换关系如图所示，其对应关系如下：of0f,ag=B（f-a）gaaf,b、丫（f-b）gbbf,L式中，、分别表示输入及输出图象，o、卩、丫为折线段的斜率，、决定低、中、高灰度级的范围。选择不同的o、卩、及数值，可得到各种各样的灰度拉伸效果。灰度变换前后的灰度变化范围

5、是不变的，对一部分灰度区域的扩张（感兴趣区）都是以其它区域的压缩为代价的。可见，输出和输入图象之间各点的灰度是按照一定的映射关系相联系的，这种映射关系在计算机中则是通过一个查照表（，即）实现的。通过对照改变了图象中不同灰度特性趋于的对比度或反差（），达到改善视觉效果的目的。直方图均衡直方图均衡（）就是通过点运算使输入图像的灰度分布较为均匀，使图像具有较好的视觉效果。设分别为原图和新图的灰度，申（P分别为原图及新图的概率密度函数，则均衡变换为原图像的累积分布函数：实验一数字图像处理基本操作及空域滤波s=T(r)=J,(r)0r对于离散图像，均衡转换公式为：sT(r)=迓Pkkj0Lk(r)=ma

6、xnjMjj=0实验一数字图像处理基本操作及空域滤波实验一数字图像处理基本操作及空域滤波)5。其中，L指图像中的最大灰度值(对于灰度图像就是maxF面给出直方图均衡化增强图像对比度的MATLAB程序:上imread(原图像名.gif);%读入原图像J=histeq(I);%对原图像进行直方图均衡化处理Imshow(I);%显示原图像Title(原图像)；给原图像加标题名Figure；imshow(J);%对原图像进行屏幕控制；显示直方图均衡化后的图像Title(直方图均衡化后的图像);%给直方图均衡化后的图像加标题名Figure;subplot(1,2,1);%对直方图均衡化后的图像进行屏幕控

7、制；作一幅子图作为并排两幅图的第1幅图Imhist(I,64);%将原图像直方图显示为64级灰度Title(原图像直方图);%给原图像直方图加标题名Subplot(1,2,2);%作第2幅子图Imhist(J,64);%将均衡化后图像的直方图显示为64级灰度Title(均衡变换后的直方图);%给均衡化后图像直方图加标题名处理后的图像直方图分布更均匀了，图像在每个灰度级上都有像素点。从处理前后的图像可以看出，许多在原始图像中看不清楚的细节在直方图均衡化处理后所得到的图像中都变得十分清晰。,用离散卷积表示4、图像平滑和锐化的空间域方法都可表示为模板卷积的形式。如可将邻域平均处理看作一个作用于图象的

8、低通滤波器，该滤波器的脉冲响应为如下f(m,n)=f(mr,ns)H(r,s)r=ks=lm,n=0,1,2,N1k由邻域大小决定，一般取(即窗口)；为加权函数，称为掩模()或模板。这种离散卷积运算等价于用模板中心点逐一对准每一象素，然后将模板元素和它所“压上”的图象元素对应相乘，再求和，其结果就是该中心点象素平滑后的输出f（x,y）。图像锐化也可以用类似的方法表示。常用的平滑模板算子有：111111121H11111,H21121,H12421923111111121A111A01/40h41017h5-1/411/4485211101/40用于锐化或检测边缘的常用模板算子有：0-10拉普拉

9、斯算子：-14-10-1010-1-1-2-1Sobel算子：20-200010-112110-1-1-1-1Prewitt算子：10-100010-111110-1-1-2-1Isotropic算子：20-200010-1121实验步骤读入图像“moon.tif”将其显示在一个可以显示2幅图像的窗口中的左边，求出“moon.tif”的轮廓，并将其显示在上述窗口的右边。读入图像“cameraman.tif”将该图像显示在一个可以显示4幅图像（两行两列）的左上角显示该图像；求出该图像的补图像，将结果显示在右上角；将原图分别顺时针旋转45度和逆时针旋转45度，前者保持输出图像大小不变，后者保持图像

10、完全显示，将结果分别显示在左下角和右下角。3.读入图像“cameraman.tif”,将其分别放大到1.5倍（用“bilinear”方法）和2.6倍（用“bicubic”方法），显示在不同的图像窗口中。2502、4.读入图像“moon.tif”将其放大2倍，并将结果用BMP格式以“moonl.bmp”为文件名，写入磁盘。（用imwrite命令，该命令的用法请用helpimwrite查看）读入图像“pout.tif”对获取的数字图像进行分段线性拉伸，参数自己任意给出。显示“pout.tif”的原图、经过均衡化处理过的图像及相应直方图。对图像“cameraman.tif”采用实验原理中给出的算子进

11、行各种增强和平滑处理并显示结果。实验思考题i直方图是什么概念？它反映了图像的什么信息？2直方图均衡化是什么意思？它的主要用途是什么？3邻域平均和中值滤波各适用于那些图像平滑处理？4拉普拉斯算子、Sobel算子、Prewitt算子和Isotropic算子进行图像增强处理时各自的特点是什么？实验报告要求给出图像处理的源程序和结果。对实验结果进行简要分析。例程：1、f=imread(cameraman.tif);subplot(1,2,1);imshow(f);subplot(1,2,2);imcontour(f);5010015020050100150200250f=imread(camerama

12、n.tif);subplot(2,2,1);imshow(f);subplot(2,2,2);imshow(imcomplement(f);subplot(2,2,3);imrotate(f,45,bilinear,crop);subplot(2,2,4);imrotate(f,-45,bilinear,loose);3、f=imread(cameraman.tif);imshow(imresize(f,1.5,bilinear);figure;imshow(imresize(f,2.6,bicubic);4、f=imread(moon.tif);imwrite(imresize(f,2),m

13、oon1.bmp);5、6、f=imread(pout.tif);f1=histeq(f);imshow(f);figure,imhist(f,64);figure,imshow(f1);figure,imhist(f1,64);m,n=size(f);f2=double(f);fori=1:mforj=1:niff2(i,j)=70&f2(i,j)=160f3(i,j)=0.5*f2(i,j)+215;endendendendendfigure,imshow(f3,);2502、2502、2502、7、f=double(imread(cameraman.tif);h1=111;111;111

14、./9;h2=111;121;111./10;h3=121;242;121./16;h4=111;101;111./8;h5=00250;0.251025;00.250./2;lap=0)-10;-141;0-10;sobel1=10-1;20-2;10-1;sobel2=sobel1;prewitt1=10-1;10-1;10-1;prewitt2=prewitt1;isot1=10-1;1.4140-1.414;10-1;isot2=isot1;2502、subplot(3,4,1);imshow(imfilter(f,h1),);subplot(3,4,2);imshow(imfilter(f,h2),);subplot(3,4,3);imshow(imfilter(f,h3),);subplot(3,4,4);imshow(imfilter(f,h4),);subplot(3,4,5);imshow(imfilter(f,h5),);subplot(3,4,6);imshow(imfilter(f,lap),);subplot(3,4,7);imshow(imfilter(f