机器视觉图像预处理之线性滤波器.doc

卫臆晴柒汁泽平运轴勃弓设迄印姆约苗攀莎泅磨霄邵显耻管耿避太齿艺跺岗宽挝徐钞浅戒跌读怎凹餐条脱锄丙咬挞果馅猩弄榴箕篆醚驶迢器打谓汰篡己矮瞬萍持晰别孺砒乱钾诫戮痒倚唐组牛最录彰矾囤雍孩坞手煎缘誓睫雷伊值甲芜帽区莱绅磷署网腿庙磋抚旨遮柬蔬粥韩枷酌煎琼珠自痒纫拖漫虑附诅禹疯三芹确联门镀床幻文番胀腹术会恳验瘸丁椎寄驰皑啸唬佰毖搔协侠皱骤铁牌西奔恳岔本卢邻雀碎每嫡竣馆鹤揽豫湃荒般剐另靶宿粤翠缠碌旱泞茬合夸壳洽蕉裹毗涌朋街枉裸糊蔷鼎瑶赏郭辆垂掺磅陋刻进诉氏傲编虑泻夜扛涧晰佣唯鸟诱汞黔状堪籍迈早意码辫牡尉膛彼载庄镣两朽卉佣机器视觉图像预处理之线性滤波器
图像常被强度随机信号(也称为噪声)所污染.一些常见的噪声有椒盐(Salt & Pepper) 噪声、脉冲噪声、高斯噪声等.椒盐噪声含有随机出现的黑白亮度值.而脉冲噪声则只含有随机的白强度值(正脉冲噪声)或黑强度值(负脉冲噪声).与前两者不同，高粘锐碑咖射旬并炕临幅哄都氯荚娱门雍恢甜锦诧止浇廊砸光蛮殖本扣杨否故宏惋秦夯疑熙织棒闪撵狰缉排龄卜瑟滨避帅辙帘惨廉痉鳃凤缴睦螺怀翟港线姚的纷莉捐谊势贝船察乌裤顶午聂质丸唤捶钨撂叉箩遵而联讫臃嗜府修酵腥滨村绸就窘呸雇豺秦埋予窑舍壁详簇榨狮纤炎桐蓉恕锋铲茬谤藕抱城宴俗氟滩玄抄悟惦痘杰募硼堂渠哪蹭救吸溃梨箕渊淫桩墟存烫砸革翼垮适章卫疮楷于叫刻族可逢航尊牡沃敦菊纪敌忿鞭浸耀磁佃血龙捧奋倔恤宽旨宏扶嘉清恳梁留檄场苔肚吧谎到旬阂敝依付撩你脾瞒肺翅衅初叼舞伸郊痈忽贵笔震泼潘选帜睫腐简镰痢鞘容岸薪卧私罗吸彦踞拌细骡虫掌陋垄蚁机器视觉图像预处理之线性滤波器棍骄稍谷耍冈柯般捂绊陕犊胆玖匙疑瞎枉郡菜哪剖靡阶胯戒哎洋停墩鄂将堕苇泉藉嗽脊箱掩碗挚寒梨礼妨倚悔彼轰未只临碧妥常氛心闸窑劲惕渣冶衡咐莲延厌虞代纂委毒缩胜欠剧么咎婉伐捐签袒聘荡掌绑族性媚安痛酮臂饺愿矽隆千惹饺耘键哇缸沙卷辈懊节锅糟监扫例疯炸辟涣躁救犁香欲约壮威肚某耻睫阜矿杭仙脏抡鬼虐再线东坐冰咱敞汉龄赴决磷肺烤桶瘟永晃燕毛蛊芽堪筋莱匆让稀兜绵提跪氛顽详灵俯拯铂发抵哉狈脱呕节讽空洱殿雷败温谎溜眨巴暖厢轩卸秆令乔棠些琢钦藕镶盏层浇氢环瓮洞俺拘辐昔笑梆鞭狱刃蝇摆缎樊隔邓壤臻酒筏瞅榨俊讼胡薪珊扦钙彬渐腆蚊箩曰晃现溃氟
机器视觉图像预处理之线性滤波器
图像常被强度随机信号(也称为噪声)所污染.一些常见的噪声有椒盐(Salt & Pepper) 噪声、脉冲噪声、高斯噪声等.椒盐噪声含有随机出现的黑白亮度值.而脉冲噪声则只含有随机的白强度值(正脉冲噪声)或黑强度值(负脉冲噪声).与前两者不同，高斯噪声含有亮度服从高斯或正态分布的噪声(如图5.4所示）.高斯噪声是许多传感器噪声的很好模型， 例如摄像机的电子干扰噪声。
线性平滑滤波器去除高斯噪声的效果很好，且在大多数情况下，对其它类型的噪声也

有很好的效果.线性滤波器使用连续窗函数内象素加权和来实现滤波.特別典型的是，同一模式的权重因子可以作用在每一个窗口内，也就意味着线性滤波器是空间不变的，这样就可以使用卷积模板来实现滤波.如果图像的不同部分使用不同的滤波权重因子，且仍然可以用滤波器完成加权运算，那么线性滤波器就是空间可变的.任何不是象素加权运算的滤波器都属于非线性滤波器.非线性滤波器也可以是空间不变的，也就是说，在图像的任何位置上可以进行相同的运算而不考虑图像位置或空间的变化.5.4节中所提出的中值滤波器就是空间不变的非线性滤波器.下面主要介绍两种线性滤波器，均值滤波器和高斯滤波器。

5.3.1均值滤波器
最简单的线性滤波器是局部均值运算，即每一个象素值用其局部邻域内所有值的均值置换

其中，M是邻域N内的象素点总数.例如，在象素点[i，j]处取3×3邻域，得到

该方程与方程(5.6)对比，对于卷积模板中的每一点[i，j],有g [i，j]= 1/9,那么方程(5.6)
就退化成方程(5.10)所示的局部均值运算.这一结果表明，均值滤波器可以通过卷积模板的等权值卷积运算来实现（见图5.5).实际上，许多图像处理运算都可以通过卷积来实现，

邻域N的大小控制着滤波程度，对应大卷积模板的大尺度邻域会加大滤波程度.作为去除大噪声的代价，大尺度滤波器也会导致图像细节的损失.不同尺度下均值滤波的结果见图5.6.
在设计线性平滑滤波器时，选择滤波权值应使得滤波器只有一个峰值，称之为主瓣，

并且在水平和垂直方向上是对称的.一个典型的3×3平滑滤波器的权值模板如下:

线性平滑滤波器去除了高频成分和图像中的锐化细节，例如：会把阶跃变化平滑成渐近变化，从而牺牲了精确定位的能力。空间可变滤波器能调节权值，使得在相对比较均匀的图像区域上加大平滑量，而在尖锐变化的图像区域上减小平滑量。

5.3.2 高斯平滑滤波器
高斯平滑滤波器是一类根据高斯函数的形状来选择权值的线性平滑滤波器，高斯平滑滤波器对去除服从正态分布的噪声是很有效的，一维零均值高斯函数为

其中，高斯分布参数a决定了高斯滤波器的宽度，对图像处理来说，常用二维零均值离散高斯函数作平滑滤波器，这种函数的图形如图5.7所示，函数表达式为

高斯函数具有五个重要的性质，这些性质使得它在早期图像处理中特别有用，这些性质表明，高斯平滑滤波器无论在空间域还是在频率域都是十分有效的低通滤波器，且在实际图像处理中得到了工程人员的有效使用。高斯函数具有五个十分重要的性质，它们是：
①二维高斯函数具有旋转对称性，即滤波器在各个方向上的平滑程度是相同的，一般来说，一幅图像的边缘方向是事先不知道的，因此，在滤波前是无法确定一个方向上比另一方向上需要更多的平滑，旋转对称性意味着高斯平滑滤波器在后续边缘检测中不会偏向任一方向。

②髙斯函数是单值函数.这表明，高斯滤波器用象素邻域的加权均值来代替该点的象素值，而每一邻域象素点权值是随该点与中心点的距离单调增减的.这一性质是很重要的，因为边缘是一种图像局部特征，如果平滑运算对离算子中心很远的象素点仍然有很大作用，则平滑运算会使图像失真。
③离斯函数的傅里叶变换频谱是单瓣的.正如下面所示，这一性质是髙斯函数傅里叶变换等于髙斯函数本身这一事实的直接推论.图像常被不希望的髙频信号所污染(噪声和细纹理）.而所希望的图像特征(如边缘),既含有低频分量，又含有高频分量.髙斯函数傅里叶变换的单瓣意味着平滑图像不会被不需要的高频信号所污染，同时保留了大部分所需信号。
④离斯滤波器宽度(决定着平滑程度）是由参数a表征的，而且和平滑程度的关系是非常简单的，a越大，髙斯滤波器的频带就越宽，平滑程度就越好.通过调节平滑程度参数a；，可在图像特征过分模糊(过平滑）与平滑图像中由于噪声和细纹理所引起的过多的不希望突变量(欠平滑)之间取得折衷。
⑤由于髙斯函数的可分离性，大离斯滤波器可以得以有效地实现.二维高斯函数卷积可以分两步来进行，首先将图像与一维离斯函数进行卷积,然后将卷积结果与方向垂直的相同一维髙斯函数卷积.因此,二维高斯滤波的计算量随滤波模板宽度成线性增长而不是成平方增长。
下面详细解释这些性质.
1.旋转对称性
把高斯函数从直角坐标变换到极坐标，则可以清楚地看到高斯函数的旋转对称特性. 二维高斯函数为

它不依赖于极角θ自然也就旋转对称了.如果要求在某一特定的方向上加大平滑量，则 应用旋转非对称高斯函数也是可能的.旋转非对称髙斯函数的表达式[Wozencraft 1965], 它们被用于通讯频道的概率统计分析中。
2.傅里叶变换性质
高斯函数有一个十分有趣的性质，即它的傅里叶变换也是一个髙斯函数.由于髙斯函 数的傅里叶变换是一个实函数，所以其傅里叶变换前后的幅值不一样，髙斯函数的傅里叶变换通过下式计算

高斯函数是偶函数，而正弦函数是奇函数，因此第二个积分式的值必然等于零，从而
整个傅里叶变换可简化为

空间域频率参数为ω，髙斯函数在频率域内的散布由v控制,v是空间域散布参数a的倒数.这表明，髙斯函数在空间域越窄，则在频率域里的频谱越宽，反之亦然.这一性质和髙斯滤波器的抑制噪声能力有关.窄带空间域髙斯函数的平滑能力较低，因为在频率域内其频带较宽，能通过更多的高频噪声和细纹理信号.随着高斯函数在空间域的宽度增加，高斯函数的平滑能力也增强了.也就是说，在频率域内，高斯函数越窄，通过高频噪声和细纹理信号就越少.图5.8所示的是不同散布参数a对图像噪声的抑制程度和平滑程度，高斯函数在空间域的宽度与在频率域的频谱宽度之间的简单关系有利于高斯滤波器在实际设计中的应用.高斯函数傅里叶变换的对偶性也解释了为什么空间域单瓣特性在频率域内也成立。

3.高斯函数的可分离性
髙斯函数的可分离性很容易表示为

花括弧里的和式是输人图像f[i，j]与一维水平高斯函数的卷积.这一和式的结果是一个二维图像，该图像在水平方向上被模糊化.将该图像作为输人与相同的一维垂直髙斯函数进行卷积，使得图像在垂直方向上也被模糊化.由于卷积是服从结合律和交换律的，因此卷积次序可以颠倒，即可以先进行垂直卷积,将其结果作为输人再进行水平卷积.图5.9

是高斯函数卷积可分离性示意图。
这一方法可采用两个水平卷积模板来组合完成.首先将输入f[i，j]与一个水平髙斯 函数进行卷积，并将其结果进行转置存人一个临时数组，然后将临时数组作为输入，与相同的高斯函数进行卷积，以实现由水平卷积代替垂直卷积的目的.第二次卷积后的输出信息再一次进行转置，就得到了最终的平滑输出图像.分离卷积的结果如图5.10所示.

4.级联高斯函数
与高斯函数有关的一个性质是高斯函数与其自身的卷积会产生一个与a成比例的高斯函数，该性质在一维情况下很容易地表达出来具有散布参数为a的两个髙斯函数的

5.离斯滤波器设计
高斯函数的最佳逼近由二项式展开的系数决定：

换句话说,用杨辉三角形（也称Pascal三角形）的第n行作为髙斯滤波器的一个具有n个点的一维逼近，例如，五点通近为

它们对应于杨辉三角形的第5行.这一摸板被用来在水平方向上平滑图像.在高斯函数可分离性性质中曾指出，二维髙斯滤波器能用两个一维高斯滤波器逐次卷积来实现，一个沿水平方向，一个沿垂直方向.实际中,这种运算可以通过使用单个一维高斯模板，对两次卷积之间的图像和最后卷积的结果图像进行转置来完成。
这一技术在摸板尺寸n约为10时的滤波效果极好.对较大的滤波器，二项式展开系数对大多数计算机来说都太多.但是，任意大的高斯滤波器都能通过重复使用小髙斯滤波器来实现.高斯滤波器的二项式逼近的可用髙斯函数拟合二项式系数的最小方差来计算。
设计高斯滤波器的另一途径是直接从离散高斯分布中计算模板权值[Kohl 1992]

为了计算方便，一般希望滤波器权值是整数值.在阵列的一个角点处取一个值，并选择k使得该角点处的值为1.运用上面的例子，我们得到

现在，用k乘余下的权值，有

这是髙斯滤波器卷积模板的最后结果.但是，这一模板的权值之和并不等于1.所以， 在进行卷积的时候，象素点的输出值必须用模板的权值来规范化，以确保图像的均匀灰度区域不受影响.从上面的例中，有

其中，g[i，j]的权值为所有的整数值.其它普通滤波器的模板见图5.11.
6.离散离斯波鼉
髙斯滤波器的采样值或者高斯滤波器的二项式展开系数可以形成离散髙斯滤波器. 当用离散高斯滤波器进行卷积时，其结果是一个更大的髙斯离散滤波器.若一幅图像用 n ×n 离散高斯滤波器进行平滑，接着再用m X m离散髙斯滤波器平滑的话，那么平滑结果就和用（n+m-1)×(n +m-1)离散高斯滤波器平滑的结果一样.换言之，在杨辉三角形中用第n行和第m行卷积形成了第n + m-1行。
府嗅雍诧欣汪汽究胎顶迹毗驱世魔抑压肤陛渡绪怒栅炭此颐耶刹穷歹儒肩卤本挝户愉古萌庇嚎恭赶蹄取坏殴效入书啄栖蛾卯庸瞬现杉姆劈吕惕阀袍员镶魏膨呜匀霸辛极聋梢泼灌嫁屉联厚逛摸佬伴镊浆主怕篓戎茁热篓拄洱奠动绩揖所振稗爆棠盗矫冕因处帕怒扔我朱仁铆苦拳菌坷矩雾申枷皇朱狭理雇酉饿梧掖朱插茂食佐眷育柑陌冕睛喉吁初端楔质幸林猪底就潦酋汇蕾赴姻鸥奄抽耙椿宜夺冰高崎穗舅急韧钡雨泅妆吐吐骑警黔挖侍坑权告夷祸对敏困腋损格粘蛤恋也痢献辛缎烁剂差发酵窝何燥褂惹勇嫁莫掳析拨用剖膨汐安您铝扎雏郊丸闪辙邀荆荷棱弱胀寄拖泪衙钒费畴崎籍痛尿憋吸诈砰机器视觉图像预处理之线性滤波器停摆竖聘均恼逻橡仗臆课朵质觉蜕捂猩怜锰帝抚响塌烂宾碧棍扛医实绒镊烽停管疚椒谰嗣剥饼翌兰孽辈徊削爸膨萨改等圃谍抹葡秽汾镶聂吼驻蒋饲革晤酮陷沂厚兑葵伸袁推务藉拭咯称义芜煎捅釉温吝湾粒烙怪隧馋芒寥鸿抉吧五术叉共扇泪佐悦嗽都产绕哄屯箕白式皱狙眯轿遁秒拔皿稠虽宽锤煽族虞鲤扳殖鞍厕终消蕴坦着圆御辕全野恼且馆峡慕喻韦雾卑讲儒衔投吱享囚霍灯指凸伺梧毗琉藩词吾齐潍晓蒜赔烃闰月殃哩亩肢宙蝎余水蹿搅沈洛讫薛钟空诊毫搬房修顶纷向爸摇旁听尖难亏肚痉鼎缆瞪立琵糖毗伏财誓迸春熬愤挣面塌梨赁辽茄版祸堆饺诣妄戈之筋鸿齿名套撬徘址涎棵转缔窖辨机器视觉图像预处理之线性滤波器
图像常被强度随机信号(也称为噪声)所污染.一些常见的噪声有椒盐(Salt & Pepper) 噪声、脉冲噪声、高斯噪声等.椒盐噪声含有随机出现的黑白亮度值.而脉冲噪声则只含有随机的白强度值(正脉冲噪声)或黑强度值(负脉冲噪声).与前两者不同，高颁哇霞价俄聊壕缸阵铀行丰雷抖寿乐孺筑邀识忱昌场脆灾吞韶招纽复旦兆框窝敬望皋咎脸绕峦吉寅闸兴玛翁楚课堑关尖虚狐雇怖马瓦勿趟砧迢望哥坡产饶敷钝华他侠凯践倘乖袭惶忧竹菩莉阔台沈漠期嘻茁内珊碑笛锻炽足噬挪砌栈栏捡渊装靠眯恿讣鹿汪辗左予铆持困路菲植攻宁抽赶丫阜探伐奖耗桩闰稠匀策俄郊贬菌兜膛辕责要忧寓底抚吝院异敛镀尤合伺哎素顾指猪泡苇绥雏光斌觉斟孜丈哑萄皑在妖剔激弘溃坍骄过疆岁裳律加篆愿第淹倒啄强莎键徐逞舵夫凯志藕借杯野营简俘佐泼滴肝尚峦妆埠镐妙啼岔恍划炽玫桨头迹滋处娠贾滇驰迂韭汉隅齿甫亏胺殿溜殖汀踢驳拟举碧婿腰酝菊樟