1、循环中值平均滤除遥测地形粗差摘要 对于地形遥测数据中引入的脉冲噪声，传统的滤波算法往往需要选择滤波门限，自适应能力不强，或者不能在滤除噪声的同时有效地保护信号数据，或者运算复杂。针对以上问题，文中提出了一种快速的脉冲噪声自适应滤除算法，该算法基于数理统计思想，对局部数据（滤波窗口）进行均值和方差估计，并根据估计结果自动选择检噪门限，进而实现噪声检测和平滑。实验结果显示，在脉冲噪声密度小于 5%时，该算法的滤波信噪比增益远高于常用滤波算法和其他同类算法；对平稳变化信号的处理效果较好。关键词 地形遥测 粗差探测 自适应滤波 数理统计1 引言地形遥测在地图制作、国土资源调查和地形辅助导航等民用和国防
2、领域都有广泛的应用。随着科技的进步，地形测绘已由传统的人工测绘发展为目前利用星载和机载传感器自动测绘数字地图。在传感器遥测地形数据的过程中，不可避免地会引入某种程度的脉冲噪声，对于高精度传感器而言，脉冲噪声是测图过程的主要误差源，必须加以探测和处理。以平面扫描型机载侧视激光雷达为例，其扫描过程是在载机飞行过程中，对下方地形从左至右逐行扫描并获取对地斜距，根据测距数据、飞行位置和取向可生成航带区域的数字高程地图。因此，机载激光雷达往往用于数字地图测绘和地形辅助导航。如果把雷达测距数据按扫描行顺序排列，则可视作距离数据图像。因此，图像处理的各种方法都可应用于这类地形区域覆盖型传感器。目前常见的图像
3、脉冲噪声去噪工具有中值滤波器（Median Filter） 、粗差探测算子（如 Robert 算子和 Laplacian 算子）1、低通滤波器2，以及各种自适应滤波器。文献3利用滤波窗口中间数据的平均来估计窗口均值，并采用软门限探测脉冲噪声，实质是中值和均值混合型滤波器。文献4采用了与文献3类似的脉冲噪声探测方法。文献5则利用模糊隶属度函数对窗口数据加权平均，通过迭代得到收敛的加权均值，进而实现一种不等权的均值滤波。文献6 提出了一种自适应模糊检噪平滑算法（AFND） ，该算法借鉴文献34 方法首先对粗差进行探测，然后对确定为噪声的像点使用改进的文献5 方法进行模糊加权均值平均以实现噪声点平滑
4、；文献6 对文献5方法的改进主要体现在两个方面：一是基于局部参数估计在滤波前进行检噪预处理，二是将正态信号的分布性质与模糊隶属度函数的参数联系起来，从而解决检噪和滤波中的参数自动设定问题。中值滤波、均值滤波等平滑类算法不能提供对信号数据的有效保护；而粗差探测算子（硬门限）和采用软门限的检测类算法的检噪效果取决于门限或门限参数的选取是否适当，低通滤波同样需要对信号频带的先验估计，实质上其自适应能力都十分有限；AFND 算法尽管在前两个方面做得较好，但运算复杂度高，实时性得不到保证。因此，对于需要较高滤波输出信噪比和自动化程度的实时地形测绘工作（如地形辅助导航） ，上述方法均不能得到令人满意的滤波
5、性能。本文将地形数据图像的局部参数估计与传统的均值/中值滤波的简单算法有机地融合起来，构造一种自适应检测和平滑脉冲噪声的快速算法，该算法可以自动估计数据图像中的脉冲噪声级别是否下降到了允许范围，在条件不满足时自动循环迭代。下文从讨论如何利用图像的局部统计特性进行粗差检测着手，提出了循环中值平均检噪平滑算法，并与文中提到的几种其他滤波器进行了性能对比实验。2 脉冲噪声的统计法检测不失一般性，假定无污染遥测图像数据服从正态分布。对大量高程数字地图和激光测距图像信号的直方图统计表明，其数据近似服从正态分布，因此这一假设是合理的。对于局部图像（滤波窗口）中的信号数据，可看作是在该局部范围内无限密集采样
6、总体的一组取样，因而窗口中的无污染信号数据是局部正态总体的一组样本，且（m,2） 。类似地，窗口内的脉冲噪声数据可视为局部噪声总体的一组取样。为讨论方便，假设全局图像被密度为 Pn 且服从（a, b）上均匀分布的脉冲噪声所污染，对于均匀分布，局部噪声总体具有与全局噪声总体同样的分布规律。这样，检噪就成为随机分布总体的两类划分问题，此类问题多应用概率分类法解决。滤波窗口内和的概率密度函数和概率分类的判决门限如图 1 所示：图 1 窗口数据的概率密度和判决门限从图 1 可见，总体和的分布往往是交迭的，无污染数据的判决门限界定为 mL，超出此范围的数据被判为脉冲噪声。对于两总体的交迭分布，无论判决门
7、限如何选取，总存在两种错误判决的可能：第一类错误：xT，却被误判为 xT ，设其概率为1；第二类错误：xT ，却被误判为 xT，设其概率为2。1 实质上表征了概率分类对信号数据的保护能力，而2 则是对概率分类检出脉冲噪声能力的度量。其定量描述为：（1）公式（1）中的 是标准正态概率分布函数，由公式可得到、和的定性关系：（）在 、n、a、b 不变的前提下，增大时，减小、增大；减小时，增大、减小；无论如何选取，都不可能同时减小和。（）判决门限的选取与总体的方差有很大关系，在保持 P1 不变的前提下，T 的方差 越小，L 值就越小。相对于整幅图像而言，局部总体由于离散程度 较小，因而在一定的P1 下
8、使得较小，可获得较小的 P2。由于图像的总错分概率是所有局部错分概率的平均，因而采用局部概率分类能获得比全局概率分类更优的检噪效果。为计算方便，固定各窗口的门限比常数，并根据统计推断原则将 R 取为23，这样可使得仅约 0.26%5%的数据点将被作为脉冲噪声进行平滑处理。这样，在脉冲噪声密度不高的情形下（高精度数据采集传感器工作在正常情况下） ，检噪算法对于未污染数据的保护性能较好。对于各滤波窗口中心点，用局部概率分类法进行噪声判断，需要估计局部均值 m 和局部方差 ，本文采用循环中值平均方法实现局部参数估计。3 基于循环中值平均进行检噪平滑对用于地形辅助导航的激光图像进行粗差滤波的要求是既能
9、去除图像以外的脉冲噪声，同时又要保持图像的细节不被改变；此外还要求算法实时性强。AFND 算法虽有较好的滤波效果，但涉及指数和迭代运算，运算速度低。中值滤波是最常用和简单的粗差滤波方法，但与以上激光像滤波要求相比也有差距。第一，中值滤波法及其改进型78虽能减少图像中的脉冲噪声，但同时也改变了有效的图像数据。对此，文献5采用了先进行噪声判断，再作针对性处理的方法。第二，中值滤波需要大量的排序操作，实时性还不够好。对此，文献3采用的方法是取窗口中除去最大和最小值之后的中间数据的平均来代替中值。从这两种方法得到启发，提出利用中值平均进行局部参数估计，并进行统计法检噪滤波的方法，算法如下：（1）以当前
10、输入检噪点 I（i , j）为中心选取 nn 的窗口（一般可取33 窗口） ；（2）找出窗内的最大值 max、最小值 min，并计算出窗内所有像点值之和 S 以及像点值平方和 S2；（3）用窗内中段数据估计窗口非噪声数据的均值：m=(S-max-min)/(nn-2)；（4）用窗内中段数据估计窗口非噪声数据的方差：=(S2-max2-min2)/(nn-2)-m2；（5）粗差检测和平滑输出：如果 |I(i,j)-m|2：则输出 O (i , j)=m否则：输出 O (i , j)=I(i , j)该算法假设在一个滤波窗口内的高值脉冲和低值脉冲数量各不超过一个，当窗内脉冲数量不符合这一假设时，均
11、值和方差估计误差增大，造成漏检。但由于每次估计都排除了最强脉冲的干扰，使估计更接近理想并检测平滑一定数量的噪声，如果反复执行这一过程，参数估计误差将趋于下降，滤波效果得到提高。因此设置了循环滤波方案，即当第一轮滤波检测出的噪声点总数与图像总像点数之比大于某一数值（循环阈值）时，进行第二轮滤波。以下为讨论方便，称这种方法为 CMAND 算法（Circular Median Average for Noise Detecting） 。理论上，33 的滤波窗口适合于粗差点的密度小于 10%的情形，考虑到噪声分布位置的不均匀并通过实验，将循环阈值设置为 0.05。 实验结果及讨论为了检验 CMAND
12、算法的性能，利用两组受到脉冲噪声不同程度污染的激光地形测距图像（见图 2、图 3 所示，以灰度高低表示距离远近）将CMAND 算法与其他算法进行滤波性能比较。由于从视觉直观上看不出各种滤波算法在处理效果上的明显差异，为了综合衡量滤波算法的数据保护和去噪效果，采用滤波信噪比增益（GSNR）作为滤波效果的客观评价性能指标，GSNR 定义为滤波器输出信噪比与输入信噪比的比值，当噪声总能量远小于信号能量时，GSNR 简化为输入噪声功率与输出噪声功率之比，即：公式中 x、xI、xo 分别代表无污染图像、脉冲噪声污染图像和滤波输出图像的像点数据。 图 2 第一组滤波输入图像图 3 第二组滤波输入图像采用各
13、种滤波算法，分别对上述两组输入图像进行处理，输出图像的GSNR 值对比见表 1、表 2 所示。表中 Median 表示中值滤波，Laplacian表示拉普拉斯检测算子，FWA 表示模糊加权均值滤波器5，AFND 表示自适应模糊检噪平滑算法6，CMAND 表示本文提出的循环中值平均检噪平滑算法。表 1 各种滤波器在不同粗差密度下的 GSNR 值（对图 2 滤波）DensityFilter1%2%4%6%8%10%12%16%20%Median3.917.2812.6415.7218.4622.9124.5424.3823.98Laplacian1.563.036.588.9612.5014.9818.2324.2729.63FWA3.576.4110.9613.1415.5418.9515.3412.6812.34AFND5.7112.5321.1623.8327.1442.0929.3425.8919.60CMAND12.1926.0031.1439.1747.9254.59
