数字图像处理实验报告:灰度变换与空间滤波(附带程序,不看后悔)
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
1.灰度变换与空间滤波
一种成熟的医学技术被用于检测电子显微镜生成的某类图像。为简化检测任务,技术决定采用数字图像处理技术。发现了如下问题:(1)明亮且孤立的点是不感兴趣的点;(2)清晰度不够,特别是边缘区域不明显;(3)一些图像的对比度不够;(4)技术人员发现某些关键的信息只在灰度值为I1-I2的范围,因此,技术人员想保留I1-I2区间范围的图像,将其余灰度值显示为黑色。(5)将处理后的I1-I2范围内的图像,线性扩展到0-255灰度,以适应于液晶显示器的显示。请结合本章的数字图像处理处理,帮助技术人员解决这些问题。
1.1问题分析及多种方法提出
(1)明亮且孤立的点是不够感兴趣的点
对于明亮且孤立的点,其应为脉冲且灰度值为255(uint8)噪声,即盐噪声,为此,首先对下载的细胞图像增加盐噪声,再选择不同滤波方式进行滤除。
均值滤波:
均值滤波是典型的线性滤波算法,它是指在图像上对目标像素给一个模板,该模板包括了其周围的临近像素(以目标像素为中心的周围8个像素,构成一个滤波模板,即去掉目标像素本身),再用模板中的全体像素的平均值来代替原来像素值。
优点:速度快,实现简单;
缺点:均值滤波本身存在着固有的缺陷,即它不能很好地保护图像细节,在图像去噪的同时也破坏了图像的细节部分,从而使图像变得模糊,不能很好地去除噪声点。
其公式如下:
使用矩阵表示该滤波器则为:
中值滤波:
滤除盐噪声首选的方法应为中值滤波,中值滤波法是一种非线性平滑技术,它将每一像素点的灰度值设置为该点某邻域窗口内的所有像素点灰度值的中值。
其过程为:
a 、存储像素1,像素2.....像素9的值;
b 、对像素值进行排序操作;
c 、像素5的值即为数组排序后的中值。
优点:由于中值滤波本身为一种利用统计排序方法进行的非线性滤波方法,故可以滤除在排列矩阵两边分布的脉冲噪声,并较好的保留图像的细节信息。
缺点:当噪声密度较大时,使用中值滤波后,仍然会有较多的噪声点出现。 自适应中值滤波:
自适应的中值滤波器也需要一个矩形的窗口xy S ,和常规中值滤波器不同的是这个窗口的大小会在滤波处理的过程中进行改变(增大)。需要注意的是,滤波器的输出是一个像素值,该值用来替换点),(y x 处的像素值,点),(y x 是滤波窗口的中心位置。
其涉及到以下几个参数:
其计算过程如下:
常规的中值滤波器,在噪声的密度不是很大的情况下(根据经验,噪声的出现的概率小于0.2),效果不错。但是当概率出现的概率较高时,常规的中值滤波处理后,仍然具有噪声点,并丢失了细节和边缘,效果不是很好。
优点:自适应滤波器不但能够滤除概率较大的椒盐噪声,而且能够更好的保护图像的细节,这是常规的中值滤波器做不到的。
缺点:实现复杂,运行时间长。
(2)图像本身清晰度不够,特别是边缘区域不够明显
对于图像清晰度不够,边缘区域不够明显的情况,应考虑使用图像锐化算法对图像边缘区域进行增强,以强调图像的边缘与细节。
为验证强调图像边缘与细节的最佳算法,我主要运用了拉普拉斯算子,sobel算子,prewitt算子,roberts算子对图像边缘进行了增强。
(3)一些图像的对比度不够
当图片自身对比对不够时,应使用对比度增强的方法对图像进行对比度增
强。
直方图均衡化:
直方图均衡化处理的“中心思想”是把原始图像的灰度直方图从比较集中的某个灰度区间变成在全部灰度范围内的均匀分布。直方图均衡化就是对图像进行非线性拉伸,重新分配图像像素值,使一定灰度范围内的像素数量大致相同。直方图均衡化就是把给定图像的直方图分布改变成“均匀”分布直方图分布。
优点:对于背景和前景都太亮或者太暗的图像非常有用,且计算量不大。缺点:变换后图像的灰度级减少,某些细节消失;某些图像,如直方图有高峰,经处理后对比度不自然的过分增强。
灰度分段线性变换:
分段线性变换函数可以增强原图各部分的反差,即增强输入图像中感兴趣的灰度区域,相对抑制那些不感兴趣的灰度区域。
优点:形式可以任意合成,增强输入图像中感兴趣的灰度区域,相对抑制那些不感兴趣的灰度区域。
缺点:是需要更多的用户输入。
(4)技术人员发现某些关键的信息只在灰度值为I1-I2的范围,因此,技术人员想保留I1-I2区间范围的图像,将其余灰度值显示为黑色。应提取目标区域的灰度值特性,将图像进行二值化,即可保留范围在I1-I2区
间范围的图像,将其余灰度值显示为黑色。
(5)将处理后的I1-I2范围内的图像,线性扩展到0-255灰度,以适应于液晶显示器的显示。
使用灰度线性拉伸的方法,将图像由I1-I2范围内线性扩展到0-255灰度。
1.2处理结果
图1 增加盐噪声的原始细胞图像
在网上下载了一张细胞图像,并对其加了盐噪声,模拟题目中所要求的明亮且孤立的点。
为了验证可滤除明亮且孤立的点的最佳方法,我选取了高斯滤波、均值滤波、中值滤波、自适应中值滤波,并使用MATLAB编写了函数代码(非直接使用MATLAB 自带函数程序),对原始图像进行了处理,处理结果如上图所示。高斯滤波几乎对噪声没有作用,故高斯滤波并不适合去除明亮且孤立的点;使用均值滤波可以使原有噪声变得不那么“明亮”(求了平均值,故白色噪声灰度值减小),但也牺牲了图片本身的清晰度,故均值滤波效果一般;使用中值滤波可以使用矩阵内中值代替灰度值极大或极小的点,大部分孤立且明亮的点均被滤除,故中值滤波处理效果较优,但中值滤波模板大小固定并不利于取得最佳结果(模板太小,不能有效滤除全部噪声;模板太大,运算时间变慢,图像细节信息减少);使用自适应中值滤波可以很好的弥补中值滤波自身的缺点,自适应中值滤波使用噪声敏感度判断所需模板大小,由处理结果可以看出,与中值滤波相比,自适应中值滤波所得结果边界更清晰,图像细节信息更丰富,故自适应中值滤波法对明亮且孤立的点处理效果最优。