数字图像处理大作业.c

西安邮电大学数字图像处理大作业学生姓名：专业名称：班级：图像增强一．引言图像作为一种有效的信息载体，是人类获取和交换信息的主要来源。人类感知的外界信息80%以上是通过视觉得到的。因此，图像处理的应用领域必然涉及到人类生活和工作的方方面面。视觉是人类最重要的感知手段，图像又是视觉的基础。早期图像处理的目的是改善图像质晕，它以人为对象，以改善人的视觉效果为目的。图像处理中输入的是质量低的图像，输出的是改善质量后的图像。常用的图像处理方法有图像增强、复原、编码、压缩等。首次获得成功应用的是美国喷气推进实验室（JPL）。他们对航天探测器徘徊者7号在1964年发回的几千张月球照片进行图像处理，如：几何校正、灰度变换、去除噪声，并考虑了太阳位和月球环境的影响，由计算机成功地绘出月球表面地图，获得了巨大的成功。随后又对探测飞船发回的近十万张照片进行更为复杂的图像处理，获得月球的地形图、彩色图及全景镶嵌图，为人类登月创举奠定了坚实的基础，也推动了数字图像处理这门学科的诞生。在以后的宇航空间技术探测研究中，数字图像处理技术都发挥巨大的作用[1]数字图像处理是利用计算机对图像进行去除噪声、增强、复原、分割、提取特征等相关理论、方法和技术的总称。因为通常图像处理是用计算机和实时硬件实现的，因此也称之为计算机图像处理。一般而言，数字图像处理的主要内容包括图像获取、图像复原、图像增强、图像分割、图像分析、图像重建、图像压缩编码等等。20世纪20年代，图像处理首次应用于改善伦敦和纽约之间海底电缆发送的图片质量。直到20世纪50年代数字计算机发展到一定水平后，数字图像处理才真正引起人们的兴趣。二．国内外发展及现状从20世纪70年代中期开始，随着计算机技术和人工智能、思维科学研究的迅速发展，数字图像处理技术向更高、更深层次发展。人们已开始研究如何用计算机系统解释图像，类似人类视觉系统理解外部世界，这被称为图像理解或计算机视觉。很多国家，特别是发达国家投入更多的人力、物力到这项研究，取得了不少的重要的研究成果。其中代表性的成果是70年代末MIT的Marr提出的视觉计算理沦，这个理论成为计算机视觉领域其后十多年的主导思想[2]。20世纪80年代末期，人们开始将其应用于地理信息系统，研究海图的自动读入、自动生成方法。数字图像处理技术的应用领域不断拓展。数字图像处理技术的大发展是从20世纪90年代初开始的。自1986年以来，小波理论和变换方法迅速发展，它克服r傅里叶分析不能用于局部分析等方面的不足之处，被认为是调和分析半个世纪以来工作之结晶。Ma11at于1988年有效地将小波分析应用于图像分解和重构。小波分析被认为是信号，图像分析在数学方法上的重大突破。随后数字图像处理技术迅猛发展，到目前为止，图像处理在图像通讯、办公自动化系统、地理信息系统、医疗设备、卫星照片传输及分析和工业自动化领域的应用越来越多。进入21世纪，随着计算机技术的迅猛发展和相关理论的不断完善，图像增强处理的应用已经渗透到医学诊断、航空航天、军事侦察、指纹识别、无损探伤、卫星图片的处理等领域。如对x射线图片、CT影像、内窥镜图像进行增强，使医生更容易从中确定病变区域，从图像细节区域中发现问题；对不同时间拍摄的同一地区的遥感图片进行增强处理，侦查是否有敌人军事调动或军事装备及建筑出现；在煤矿工业电视系统中采用增强处理来提高工业电视图像的清晰度，克服因光线不足、灰尘等原因带来的图像模糊、偏差等现象，减少电视系统维护的工作量。图像增强技术的快速发展同它的广泛应用是分不开的，发展的动力来自稳定涌现的新的应用，我们可以预料，在未来社会中图像增强技术将会发挥更为重要的作用[3]。三.传统图像增强的基本方法1基本的灰度变换实践证明，灰度变换技术是一种简便而有效的提高图像对比度的方法。灰度变换也叫点运算，它不改变像素的位置，只改变像素的灰度。设输入图像为),(yxf，输出图像为),(yxg，则灰度变换的数学表达式可表示为：)],([),(yxfTyxg(式3.1)这里T为灰度变换的具体映射关系。经常出现这样的情况，图像的灰度范围],[ba没有充分利用显示装置所允许的最大灰度范围max][min,从而导致图像的对比度太低，使一些细节不易被观察到。比如：摄影过程中如果曝光过度，就会出现这样的缺陷。解决上述问题的最简单方法是进行灰度的线性变换，其数学表达式如下：min]),([minmax)],([),(ayxfabyxfTyxg(式3.2)对灰度做这样线性变换以后，把原始图像),(yxf的灰度范围],[ba强行扩展为显示装置所允许的最大灰度范围max][min,，从而提高了整幅图像的对比度，原来观察不到的一些图像细节可能更加突出了，图3.1给出了这种线性灰度变换关系：图3.1灰度范围的线性变换如果在图像处理过程中，需要突出图像中某些灰度范围内的图像的细节，同时又允许适当损失另外灰度范围内的图像处理细节，可以采用线性灰度变换的另一种形式，即分段线性变换。经过这种变换以后，可以使所关心的图像细节的灰度范围得以扩展，增强其对比度；同时又使的所不关心的图像细节所处的灰度范围得以压缩，降低其对比度。值得注意，这种分段线性变换，变换前后整幅图像总的灰度范围不变的。三段线性变换的数学表达式可写成：),(yxg=232312121),(),(),(fkgyxfkfkgyxfkyxfk)),(()),(()),(0(2211Mfyxfffyxfffyxf(式3.3)式中111fgk，)()(12122ffggk，)()(223ffggkMM在实际的处理过程中，如果图像上灰度范围的两端区域上有噪声，比如感光胶片上有划伤和黑色感光颗粒，则可以用这种变换把灰度范围的两端区域压缩，使人眼视觉对噪声的感受不明显，而对有用细节所占据的灰度区域给予线性扩展，提高这部分的对比度。如果图像上绝大部分的像素的灰度级集中在],[ba范围内，比较少的像素的灰度级超出此范围，则可用以下变换增强原图像上的],[ba范围的对比度：),(yxg=)),(max()),(min()),((min)),((minmaxMfyxfayxfbyxfaayxfab(式3.4)图3.2表示了这种变换关系。值得注意，扩展原图像灰度范围],[ba是以完全损失灰度小于a和灰度大于b的图像节为代价的。这种变换与分段线性变换实际上都是非线性变换。实际上，可能利用一些数学函数进行灰度变换，如平方、对数、指数等但这种变换必须满足以下条件，即：如果max),(minyxf，则需有max),(minyxg，也就是说，灰度变换前后的灰度范围必须在显示装置所允许的最大灰度范围之内。下面介绍以下常用的几种变换。（1）图像反转：对图像适用于求反是将原图灰度值反转，简单来说就是把黑的变百，白的变黑。嵌入于图像暗色区域的白色或者灰色细节，特别当黑色面积占主导地位时，进行图像反转是比较理想的。（2）对数变换：使窄带低灰度输入图像值映射为宽带输出值，可以利用这种变换来扩展被压缩的高值图像中的暗像素，相对的是反对数变换的调整值；对数函数有它重要的特征，就是它很大程度上压缩了图像像素值的动态范围；（3）幂次变换：幂次曲线中的部分值把输入窄带暗值映射到宽带输出值，相反，输入高值时也成立。(2)对比度增强：增强图像比度（Contraststretching)实际上是增强原图各部分之间的反差。对比度增强分为线性和非线性对比度增强两种，线性是指将对比度较差的图像灰度线性扩展，常能显著改善图像的质量。当用某些非线性函数如对数、指数函数作为映射函数时，可实现图像灰度的非线性变换。对比度增强是图像增强中最普遍的增强方法。当图像成像不足或过度曝光，图像记录设备范围太窄等，都会产生对比不足的问题，使图像的细节分辨不清。为此需对每一像素的灰度级进行变换，扩大图像灰度的范围，达到图像增强目的。2直方图处理设图像的灰度级范围为max][min,，r为此灰度范围内的任一灰度级，)(rp为这幅图像灰度级为r的像素出现的频数，可以看出，)(rp是r的函数，该函数的图形称为这幅图像的直方图。一般横坐标表示灰度级，纵坐标表示具有该灰度级的像素的频数。数字图像直方图显然是有一系列竖线条组成的图形，竖线条的高度代表了该灰度级在此幅图像出现的频数（或相对频数或频率），如下图3.2所示：图3.2数字图像方图为讨论方便，把灰度级r归一化，即)10(10代表白代表黑，rrr，假定对原始图像灰度级r做以下灰度变换：)(rTS为使这种灰度变换具有实际意义，应满足如下条件：（1）在10r区间内，)(rT为单值单调增加；（2）在10r，对应有1)(0rT。这里条件（1）使变换后的灰度值保持从黑到白的次序，条件（2）保证变换后的像素灰度级仍在允许的范围内。由S到r的反变换可用下式表示：)(1STr)10(S这里)(1ST对于变量S也应满足条件（1）和（2）。图3.3表示了满足条件（1）和（2）的变换)(rTS。图3.3具有实际意义的灰度变换接下来简单介绍以下常见的两种基于直方图灰度变换：(1)直方图均衡化：是一种借助于直方图变换实现灰度映射从而达到图像增强目的的方法。直方图表示数字图像中每一灰度级与其出现的频数(具有该灰度级的像素数目)间的统计关系。直方图能给出图像整体分布描述，如图像的灰度范围、灰度级的大致分布情况等。把原图像的直方图变换为各灰度值频率固定的直方图称为直方图均衡化。(2)直方图规定化：也是一种借助于直方图变换来增强图像方法，它通过将原始的直方图转换为期望得到的直方图，从而达到预先确定的增强效果。也就是使处理后的图像具有指定的直方图。3中值滤波传统图像增强算法中，空域滤波是直接在图像空间借助摸板卷积来实现的增强方法。根据不同的模板设计，空域滤波即可以实现对图像平滑功能。在空域滤波时即可线性的组合模板运算的结果，也可非线性的利用模板运算的结果。中值滤波的思想是对一个窗口内的所有像素灰度值进行排序，取排序结果的中间值作为原窗口中心点像素的灰度值。这种滤波也就是平滑操作，对干扰噪声的效果较好。中值滤波的关键在于选择合适的窗口大小和形状。4小波变换图像增强小波变换(WavelateTransform)是由法国科学家Molret在进行地震数据分析时提出的，在1986年Mallat将计算机视觉领域的多尺度分析思想引入到小波函数构造中，形成了统一的小波函数构造理论。小波变换在图像分解与重建过程中，按需要改变有关小波参数，并且它的多分辨率分析具有良好的空间域和频率域局部化特性，对高频采用逐渐精细的时域或空域步长，可以聚焦到分析对象的任意细节，由此可增强图像中感兴趣的部分，国内外已有部分学者开始对此方法进行研究，到目前为止小波变化在图像压缩领域应用效果最好。四结论总体来说，以上方法都有各自的使用范围，传统图像增强方法在原理上容易理解，计算也简便。图像增强是数字图像理中的重要研究内容，不仅在于它的广泛使用性，更在于它于我们的日常生活息息相关，一方面体现在医疗仪器上；另一方面在于我们通过传感器获取的大多数图像并不如我们所想象的那么理想例如噪声太大，图像模糊，背景黑或者太亮等等，都需要我们对图像进行增强处理。参考文献[1]冯安,王希常.MATLAB在数字图像增强中的应用.信息技术,2007,23(01):14-16.[2]夏德深,傅德胜.计算机图像处理及应用.南京:东南大学出版社,2004.[3]王斌.MATLAB实现数字图像增强处理.佳木斯大学学报,2005,23(1):31-34.实验部分用Matlab编程实现以下题目中的三个。(1)亮图像、暗图像、和低对比度图像的直方图均衡化。(2)用中值滤波滤除图像中的椒盐噪声。(3)用空域拉普拉斯模板实现图像锐化。(4)用频域高斯低通、高斯高通滤波器分别对图像进行平滑和锐化。实验一一、实验题目用中值滤波滤除图像中的椒盐噪声。二．实验原理椒盐噪声是由图像传感器，传输信