堆栈滤波的递归实现.pdfVIP

堆栈滤波的递归实现 张薇沈允春张曙 (哈尔滨工程大学电子工程系哈尔滨15000D 摘要本文介绍了最小平均绝对误差准则下最优堆栈滤波嚣的递归实现方法．文中采取两种方案实 现递归：直接递归方案和量优递归方案，并分别对两种方案进行了详细阐述．量后从嗓声滤波特性 和边缘细节保护特性簪方面对递归和非递归滤波嚣进行了比较，并给出了实例． 关键词堆栈滤波递归滤波非线性滤波 f司彳垂她耳Z -rⅣ，¨_白 o2j6 。 O引言 与线性滤波器相比，非线性滤波器对于脉冲类噪声的滤除具有较好的效果，因此 近年来非线性滤波嚣广泛应用于图象处理中．堆栈滤波器是最近十几年来发展起来的 一种非线性滤波器．由于其具有并行结构，主要的滤波工作在二进制域进行，且实现 简单，因此具有重要的应用价值．从实现域角度来看，堆栈滤波器可以认为是在二进 制域实现某些非线性滤波器的一种方式；从广义上讲，秩排序统计滤波器、加权秩排 序统计滤波器、形态滤波器及中值滤波器等非线性滤波器都属于堆栈滤波器． MAE(MeanAbsolute Error)准则下最优堆栈滤波器在某一时刻的输出仅与该时刻 对应的输入序列有关，因此从根本上讲是非递归型滤波器。本文提出了直接递归和最 优递归两种递归方案对堆栈滤波器进行递归实现． 1堆栈滤波器结构和算法 1．1堆栈滤波器的结构 堆栈滤波器具有两个重要性质：阈值分解性质和堆栈性质．前者建立了非线性滤波 器的并行结构；后者则赋予并行结构各层之间一种有效的捧序关系，使滤波器实现简单． 由以上性质决定的堆栈滤波器结构如图1． 图1 堆栈滤波罂结构 一啦一 设堆栈滤波器的输入为M值非负信号矢量签=暖j，五，…，“】，z，仨{0,1，…，Mol}． 签能够分解成M-1个二进制矢量兰。，其中m=l，2’…，盯～l，为分解层层数．利用闽值 分解函数r。，置分解为： 彳钉㈤=￡凳孙 (1) 且有： 鲎=∑，。ZOo-Z-。 (2) 三一分别经各层的二进制滤波器“●滤波后，得到相应的输出矢量Y。，则十进制域的输 出兰为： r_-∑， (3) 公式(1)、(2)、(3)即为堆栈滤波器阈值分解性质的表述． 已知对于任意两个二进制矢量璺；【口。，口2，…，口．】和童=吼，6：，…，“】，若对于 q蔓以，其中Vii=l’2’…，H，存在a玉查，则称幽查之间存在“堆栈”关系．用二进 制澹波器，(●对翊垒进行滤波，令”=，(鲅v=，@，“、v为二迸制输出．当且仅当 。r幽s，@，a≤查时，称二进铡滤波嚣朋拥有堆栈性质．公式∞各个二进制矢量信 号的元素之间存在堆栈关系：_xu4蔓…s兰-s…sf．在堆栈滤波器各层上的二进制 滤波器，(．)必须是具有堆栈性质的正布尔函数，以保证堆栈滤波器的堆栈特性．正布 尔函数是指能用布尔表达式表示，表达式中不包含输入交量的“非”运算，且具有堆 栈性质的二进制滤波器，(．)．对任一已给函数集合^(n如(吐…，靠(兰)，我们可以定 义一个分离函数： 叫=∑彤以④一‰=芝’弛 其中 堡=暇，％，…，畋，畋+。】7 ≯(兰)=M(三)，如∞，…，如o，_1】】 满足条件：1)以(D是兰的正布尔函数：2)形≥0，f；l，…Ⅸ，置+1． 当且仅当1)正布尔函数能被表示为他)-c，(G∞)，孵(兰)用三表示时，正布尔函数 ，(●称为线性可分离正布尔函数；2)咖(D为非线性函数时，，(-)称为非线性可分离正 布尔函数． 所有的堆栈滤波器可以看作是带有非线性输入的线性可分离函数．由正布尔函数 ，(●定义的堆栈滤波器S，如下式所示：