[工学]第11章 滤波器设计.pptVIP

* 数字图像处理 Digital Image Processing DIP 中国传媒大学信息工程学院 * 第11章　滤波器设计 11.1 引言 　　第9、10章是分析设计线性滤波器的基础，本章则是讨论设计实现特定目标的滤波器。 11.2.1 简单的低通滤波器 　　矩（箱）形滤波器：让信号与矩形脉冲做卷积，能实现局部平均，从而减弱高频噪声的影响。它称为滑动平均滤波器。 11.2 低通滤波器 矩形脉冲信号的傅立叶变换为内插函数，f(x)越窄则F(s)越宽，谱的过零点就会高于信号最高频率2倍。 1 0 a 2a 0 1/a s 见图，由于负旁瓣的影响，使不同频率的垂直棒状图案的黑白极性出现反转。 不同频率的棒状图（多波群） 矩形滤波器宽度（应不超过2个像素宽） 　　根据相似理论，只要矩形滤波器不超过两个像素宽度（够窄），它的传递函数的过零点就会落在被采样数据的最高频率以外，否则图像中的细微结构（高频成分）就存在极性反转危险。 11.2.2　高斯低通滤波器 　　高斯函数的傅立叶变换仍为高斯函数，因此高斯函数能构成一个在时域和频域都具有平滑性能的低通滤波器。可通过在时空域卷积或频域相乘来实现高斯滤波。 低通滤波 带通滤波 高通滤波（高频增强） 图像经过不同性能滤波后的效果 　　三角形滤波器：用三角脉冲?(x)做为低通滤波器的冲激响应，又称为加权平均滤波器。 　　三角脉冲的谱具有　　　的形式，不会变负（不反转），随频率提高，衰减速度快。 　　连续两次矩形滤波可得到三角波同样的效果。 　　高频截止：这种方式的低通滤波方法由三步组成a)计算信号的傅立叶变换；b)将幅度谱的高频部分置0（频域截取）；c)求傅立叶反变换。 　　这种做法等同于用一个矩形脉冲去乘幅度谱，也可认为是用　　　与信号卷积。但它会在尖峰或边界附近出现振铃，因此高频截止滤波的用途有限。 11.3　带通和带阻滤波器 11.3.1　理想带通滤波器 　　带通和带阻是实用中的需求，保留或去掉频谱中不同频段内的成分。 s ?s -f1 1 s0 -f2 G(s) f2 f1 由于G(s)是一对矩形脉冲，可令： 低频包络 高频载波 　　理想带通滤波器的冲激响应 　　当?(s)很小时，为“切口”滤波器（notch filter），也叫陷波器。 11.3.2 理想带阻滤波器 s ?s -f1 1 s0 -f2 G(s) f2 f1 　　理想带阻滤波器的冲激响应 用一个非负单峰函数K(s)与s0处的冲激偶作卷积。 11.3.3 通用带通滤波器 假定K(s)是高斯函数，则 其中： 4 s0 G(s) -s0 高斯带通滤波器的冲激响应 显然，利用上述技术也可以很容易地构造带阻滤波器。 11.4 高频增强滤波器 11.4.1 高斯差分（DOG）滤波器 　　传递函数在零频处为1，频率升高则H(s)增长，（由于不能无限增长），至某个频率值后开始衰减，直至0。 用两个不同宽度的高斯函数的差分来产生 高斯高频增强滤波器的传递函数 高斯高频增强滤波器的冲激响应 11.4.2 设计高通滤波器的经验规则 假定滤波器的冲激响应是窄脉冲减宽脉冲，则其G(s)具有高频增强滤波器的形状。下面要估算G(s)在s=0处的值及s=?时G(s)最大。 高窄脉冲 低宽脉冲 最大值：对g(t)做傅立叶变换并令s=0 这个值恰好相当于两个s=0时的G(s)的幅度差 函数下的面积 由于G(s)非负，最大值出现在G2(0)=0时，于是 由两脉冲的差构成高频增强滤波器的经验规则 低频响应：指滤波器对图像中较大物体和恒值灰度区域的响应。假定冲激响应g(t)的持续时间有限，而输入信号f(t)在g(t)期间恒定。于是系统输出： 由于s＝0时的傅立叶变换： 因此，如果G(0)=1，那么滤波器将不改变f(x)中的大而且恒定区域的幅度（低频响应为1，即低频无改变） 对二维场合，滤波器不改变输入图像中大而平坦区域的对比度。如果G(0)?1，则该参数成为增益因子，控制着h(t)和f(t)的较大分量的总体幅度关系。 11.5　最优线性滤波器设计 　　在数字图像处理中曾经出现非常多的滤波器设计，但很多都是次优的，并不是最优。 　　次优滤波器（易于用计算机实现）可能会给图像引入人为的现象，而且常常是不可预料的。例如在一个域为矩形的滤波器，在对应的域中则会产生减幅振荡（振铃）。 　　Wiener估计对于从加性噪声中恢复未知信号是最优的，而匹配检测器对于检测淹没在加性噪声中的已知信号是最优的。 11.5.1　随机变量 用随机噪声来描述未知的受到噪声污染的信号，即只有关于某信号的一些一般知识而缺乏具体了解。 11.5.1.1　遍历性随机变量 时间均值：对某成员在时间区间上作积分来计算随机变量。样本集均值：计算所有成员函数在某个特定时间点取值的平均值。 遍