边缘保留的图像噪声滤除方法[DSP数字信号处理]

导航:老古开发网首页→文章索引→文章分类→DSP数字信号处理→[边缘保留的图像噪声滤除方法]

| -文章搜索 - 最新文章 - |

边缘保留的图像噪声滤除方法

发布时间:2006年6月4日点击次数:733

来源:电子技术应用作者:北京理工大学电子工程系（１０００８１）刘先锋沈胜宏中国科学院电子学研究所（１０００８０）郑明洁

摘要：常用的滤除图像噪声的算法由于需要在图像边缘纹理保留和噪声滤除之间进行折中,因此不可避免的会带来边缘模糊的问题。本文提出一种既可保留图像边缘纹理又滤除图像噪声的方法。

关键词：噪声滤波边缘检测

噪声图像的滤波问题一直是图像处理的基本任务之一。常用的滤波技术主要有Ｌｅｅ滤波器,Ｆｒｏｓｔ滤波器,ＧａｍｍａＧＡＰ滤波器等。这些滤波算法都是基于对图像局部统计特性自适应的,因此滤除噪声的效果较好。但由于算法本身的原因,往往造成图像的边缘等细节信息模糊,降低了图像的质量１。针对这一问题,人们提出许多改进算法,如改进的Ｌｅｅ滤波器,变窗口滤波器等,这些算法虽然在一定程度上解决了边缘模糊的问题,但也带来新问题,比如改进的Ｌｅｅ滤波器,由于要进行边缘检测,所以选择的窗口就不能太大,但小窗口对消除斑噪声不利,又降低了效能。

针对以上情况,本文提出了一种新的滤波方法：把滤波过程分为两个步骤,先检测出图像的边缘,把图像的边缘从原图像中分离,得到无边缘纹理的图像;然后对其进行滤波,再将边缘纹理加入到滤波后的图像,这样就得到了一个即保留边缘纹理,又有好的滤波效果的图像。

１噪声图像边缘检测

最常用的边缘检测算法一般是Ｓｏｂｅｌ算子、Ｋｉｒｓｈ算子、Ｌａｐｌａｃｉａｎ算子以及Ｍａｒｒ算子等,这些算法都是基于梯度检测的。梯度检测器就是确定一个门限,当象素梯度超过门限时就认为它是边缘,在图像越亮的地方,梯度的方差就越大,如果门限不变,就会有许多虚假边缘被检测出来;既使门限变化,由于门限的变化范围不好确定,因此效果也不好,所以检测被噪声污染的图像的边缘不能用一般的方法。本文使用斜面拟合法和广义模糊算子法较好地解决了这个问题。

１.１Ｈａｒａｌｉｃｋ斜面拟合法[２]

边缘是一种灰度的起落,故可用曲面对数字图像进行拟合,然后再在拟合的曲面上进行边缘检测。由于拟合用的曲面是满足一定平滑性的有理曲面,因而可以使噪声图像得到平滑,有一定滤噪能力。本文用Ｈａｒａｌｉｃｋ曲面拟合方法作图像边缘检测,该方法用斜面拼接成表面,然后求出表面的交,从而确定物体的边缘。Ｈａｒａｌｉｃｋ边缘检测可以通过两个步骤完成：

１.１.１确定Ｈａｒａｌｉｃｋ斜面拟合参数

假设对于方形数字图像Ｒ×Ｃ中的每一点（ｘ,ｙ）,数字图像ｆ(ｉ,ｊ)的斜面模型都可以写成:

ｆ(ｉ,ｊ)＝αｉ＋βｊ＋γ＋η(ｉ,ｊ)(１）

    其中η(ｉ,ｊ)是噪声项,是独立的零均值随机变量。用最小平方误差求α,β,γ,即最小化下述误差：

    将ε２分别对α,β,γ求偏微分,并令结果等于零：

    通常Ｒ,Ｃ均为奇数,设拟合窗口Ｒ×Ｃ的中间点坐标为（０,０）,由于对称性,有：

    求解上面的方程组,得：

    求出α,β,γ,后,就可以很容易地进行边缘检测。

由于１仅于Ｒ、Ｃ有关,故令

１．１．２斜面交界的确定

考虑两个相邻的斜面,它们分别由斜面参数α１,β１,γ１　和α２,β２,γ２　表征,设它们大小相等,斜面中心点连线的中点坐标为（０,０）。于是这两个中心点的坐标可设为－δｉ／２－δｊ／２　和δｉ／２δｊ／２　,则两个斜面在一个平面内的条件是：

α１＝α２,β１＝β２,(α１－α２)δｉ／２＋(β１－β２)δｊ／２＋(γ１－γ２)＝０(８)

因为噪声是正态分布的,则α,β,γ亦服从正态分布,所以：

都是服从零均值的正态分布随机变量。

构造Ｆ分布作为统计量：

Ｆ值越大,这两个拟合邻域之间存在边缘的可能性越大。

１．２广义模糊算子法[３]

传统的边缘检测方法是将边缘点理解为灰度的突变点,从而通过邻域象素之间的代数运算来求取边缘点。广义模糊算子法认为：图像灰度差异是由于光照不均而产生的,在二维图像中,边缘是高频成分,但图像在边缘处灰度值是连续的,也就是说,边缘包括了图像其它部分的灰度信息。但由于量化的影响,图像灰度值在边缘处产生突变。

定义论域Ｕ上广义模糊集合Ａ表征为：

其中υＡ(ｘ)∈[－１,１]称为Ｕ上的广义隶属函数;称υＡ(ｘ)∈[－１,０]为Ｕ上ｘ完全不属于Ａ的广义隶属函数;υＡ(ｘ)∈[－１,０]为Ｕ上ｘ完全属于Ａ的广义隶属函数;而υＡ(ｘ)＝０为Ｕ上Ａ的模糊分界点函数。

若Ｕ是由有限个元素构成的有限域,则广义模糊集合Ａ也是有限的。如果一个广义模糊集中仅有一个元素,则称为广义模糊单敦。

于是一个由灰度级为Ｌ的Ｍ×Ｎ二维灰度图像可看成是由一个广义模糊单敦构成的阵列,其每个元素的广义隶属函数的绝对值表示相对于最大亮度Ｌ－１的亮度程度,故可记为：

其中Ｐｉｊ／ｘｉｊ,(－１≤Ｐｉｊ≤１)表示图像中的元素（ｉ,ｊ）完全拥有或不拥有性质Ｐｉｊ的程度。

定义一个广义模糊算子（ＧＦＯ）,它作用在广义模糊集Ａ上可以产生另一个模糊集Ａ’,即Ａ’＝ＧＦＯ(Ａ)。给出ＧＦＯ的表达式如下：

其中β＞１,α＞０。在给定了β之后,α可以由上式中第二第三项在分界点[０．５,１]上进行耦合求出。图１为用广义模糊算子实现图像边缘检测的框图。

选择合适映射,将图１中Ｘ映射到Ｐ,本文采用正弦映射,经过ＧＦＯ变换,发现：

当象素点Ｘ（ｉ,ｊ）→０时,由ＧＦＯ操作后,Ｐ(ｉ,ｊ)→－１,Ｐ’(ｉ,ｊ)→１,而Ｘ’(ｉ,ｊ)→Ｘｍａｘ’,即低灰阶区域映射到高灰阶区域;

当象素点Ｘ（ｉ,ｊ）→１时,由ＧＦＯ操作后,Ｐ(ｉ,ｊ)→１,Ｐ’(ｉ,ｊ)→１,而Ｘ’(ｉ,ｊ)→Ｘｍａｘ’,即高灰阶区域映射到高灰阶区域;

当象素点Ｘ（ｉ,ｊ）为→１中间灰度集时,这通常是真正的边缘所在。由ＧＦＯ操作后,Ｐ(ｉ,ｊ)→０,Ｐ’(ｉ,ｊ)→０,而Ｘ’(ｉ,ｊ)→Ｘｍａｘ—Ｄ（Ｄ是由所选映射函数决定的常数）。

至此,就可用广义模糊算子（ＧＦＯ）检测出图像的边缘。

１．３用Ｈａｒａｌｉｃｋ斜面拟合法和广义模糊算子法综合得到图像边缘

用Ｈａｒａｌｉｃｋ斜面拟合法和广义模糊算子法结合检测边缘,可以先用这两种边缘监测器分别独立地对图像边缘进行检测,然后把得到的图像边缘取逻辑与,得到用两种方法相结合的最终图像的边缘。

２无纹理图像的噪声滤除

对于一般图像的噪声,由噪声的统计特性可知,其幅度Ａ是高斯分布,其亮度ｕ是指数分布。对亮度来讲,其分布函数可写成[４]：

因此实际上最终获得的噪声图像像素代表场景的反射强度,图像亮度Ｉ（ｔ）可写成

Ｉ（ｔ）＝Ｒ（ｔ）ｕ（ｔ）(１４)

其中ｔ＝（ｘ,ｙ）为图像上一点,Ｒ（ｔ）为理想图像亮度,ｕ（ｔ）为亮度,反应了目标的反射特性,ｕ（ｔ）与Ｒ（ｔ）统计独立。

这里我们采用Ｌｅｅ滤波器进行噪声滤除。Ｌｅｅ滤波器假设图像噪声是乘性噪声,并采用自适应迭代法。这种方法有利于保持边缘信息。其数学表达式为：

Ｒ＝Ｉ＋Ｋ(ＣＰ×Ｉ)

Ｋ＝(ＵＶ)／(ＶＵ２＋Ｉ２Ｍ),

Ｍ＝(ＳＤ／Ｉ)／２ (１５)

其中,Ｒ为处理后像元值;Ｉ为平滑窗口像元的平均值;ＣＰ为平滑窗口中心像元的灰度值;Ｕ是相乘噪声均值（基于假定,一般取Ｕ＝１）;Ｖ是平滑窗口像元的方差;ＳＤ是平滑窗口的噪声标准差。

３基于边缘检测的噪声滤除

对未处理的图像（如图２所示）用合适的方法进行边缘检测（使其纹理不丢失）。不是边缘的象素令其值为零,边缘象素保留其原值,得到边缘图像（如图３所示）。把边缘象素点从原始未处理的图像中减去,令边缘被减去处的值等于临近非边缘象素值的平均,这样就得到了非边缘的均匀的图像。用合适的滤波器（本文为Ｌｅｅ滤波器）对改图像滤波,滤除噪声,然后恢复边缘,即用边缘位置的象素值代替被滤波过的对应位置的象素值,得到最终图像（如图４所示）。图５为仅用Ｌｅｅ滤波的结果。从图中可以很明显的看到该算法即保留了边缘又滤除了噪声,效果要明显好于Ｌｅｅ滤波。

欢迎进入老古论坛进行讨论
[DSP数字信号处理] 相关文章:

带RF输出的全数字调制器设计

简介:
软件定义无线电(SDR)终端促进了物理层功能的可编程实现。很多研究工作已经应用DSP和FPGA实现物理层的基带功能。 SDR无线电是如此定义的电台，其天线后面的某段实现了数字化。其后电台可用柔性及可配置的功能模块来实现DSP算法。随着技术的进步，数字化可以在天线后，或非常接近于天线，如此几乎所有的无线电功能都可以使用高速可编程的DSP引擎通过软件实现。目前的无线电结合了模拟和数字构造模块。RF功能还是适合模拟电路实现，而基带功能更适合于DSP的实现。数字化前端(DFE)常用来连接RF及基带处理。DEF通常能够处理10MH左右的频率信号，常被称为数字中频(IF)。因此，重要......

多通道故障信号的采集处理及实现

视频字符叠加技术的发展及四种实现方案

基于双内存的图形采集卡的设计

用ＴＭＳ３２０Ｃ６２０１实现多路ＩＴＵ－ＴＧ．７２８语音编码标准

下一个:[新闻热点]中关村电脑节9月登场六大创新凸现全新亮点

简介:
中关村电脑节组委会于5月22日宣布，第九届中关村电脑节将于2006年9月6日～10日在北京举行。同时组委会发言人表示，由于采取了多项改革和创新措施，本届电脑节将以一个全新的面貌出现。在过去成功举办八届电脑节的基础之上，组委会经过广泛调研和认真总结论证，针对创建创新型国家和创新型城市的需求和相关产业的发展现状和趋势，对电脑节的定位做出了一定调整，并在组织方式和活动内容两个方面都进行了整合和创新。第九届电脑节不仅将突破以往区域性活动的局限，而且还将推出部分准公共产品的全球招标的创新举措，为将中关村电脑节打造成为一个全国乃至世界知名品牌活动奠定基础。本届电脑节以“创新和谐卓越--打造产......

上一个:[新闻热点]Agilent首款具有可调谐环路带宽时钟数据恢复选件