专利名称：弱小目标二维图片噪声滤除中的匹配滤波方法
技术领域：
本发明涉及的是一种图像处理技术领域的方法，具体是一种弱小目标二维图 片噪声滤除中的匹配滤波方法。
技术背景在数字信号处理中，滤波器是不可缺少的。滤波器的一个作用是使基带信号 频谱成形，例如为了满足奈奎斯特第一准则，基带信号频谱通常采用升余弦滚降 形状。滤波器的另一个重要作用是限制噪声，将信号频带外的噪声滤掉，减小它 对信号正确判决的影响。设计最佳线性滤波器时可以有两种准则一种是使滤波后的信号波形与发送 信号之间的均方误差最小，由此而导出的最佳线性滤波器被称为维纳滤波器；另 一种是使滤波器输出信噪比在某一特定时刻上达到最大，这种滤波器称为匹配滤 波器。匹配滤波的实质是要设计一个和图像感兴趣的图像信息或目标对象相关的 滤波器，经过匹配滤波提高测量信息的信噪比，进而改善目标检测的性能，即空 间匹配滤波器的最大信噪比的设计准则。设输入图像为《x,力，其中包括待识别目标和背景，这里的背景是广泛意义 上的背景，包括大气、空中的云层、干扰点、传感器等的系统噪声、随机噪声等。 待识别目标空间分布特性为《X,力，其傅立叶频谱为7X",V);背景空间分布特性 为6(X，力，其傅立叶频谱为BO,V)。并假定目标和背景不相关。艮口/(x,y)=《x，_y) + 6(x,_y) (1) 其傅立叶变换为/(",v)/(w，v) = r(w，v)+^(",v;) (2)则由匹配滤波器光学相关识别的原理，匹配滤波输出图像为目标空间分布特性与输入图像空间分布特性的巻积F-1 [/(w， v), O， v)] = ,(x， J) + 6(x, j) ,(x，力 (3) 其中(8)表示相关，广("，v)为待识别目标的匹配滤波函数。经对现有技术的文献检索发现，大量的参考文献都是利用上述的方法设计 匹配滤波器，经典的二维空间匹配滤波器由于同时考虑到了点目标的空间分布和 背景噪声的特性，利用先验的点目标空间分布特性与红外图像作巻积，当先验点 目标与红外图像中待检测的点目标相匹配时，理论上能够极大地提高图像信噪 比。在文献Matched filters for multispectral point target detection (Proc. SPIE, Vol. 6302, 63020Z (2006))中作者针对不同的光谱图像，利用上诉的方法设计匹配滤 波器，得到良好的结果。但是，在复杂背景条件下，实际的红外图像中存在强烈 噪声的影响，点目标的空间分布特性受到噪声干扰会发生一些改变，利用上述经 典的空间匹配滤波器进行处理性能将迅速下降，往往不能有效地提高信噪比，严 重影响后续目标检测的性能，往往难以达到预期的效果。 发明内容本发明针对现有技术存在的上述不足，提供一种弱小目标二维图片噪声滤除 中的匹配滤波方法，在分析成像待滤噪目标的空间分布特性的基础上，充分利用 目标、背景、噪声之间形态上的差异，以及成像待滤噪目标的空间分布在各个方 向的截面轮廓线上均保持了在目标出现位置的局部邻域内突变的特征；同时，充 分考虑到噪声干扰的随机性和局部突变性，以及噪声干扰只可能是破坏目标出现 位置二维分布中的某一维方向上的分布特性，将经典二维匹配滤波器分解为多个 一维滤波器。通过优化设计的匹配滤波器的多重一维匹配滤波结果能够有效地平 衡目标分布特性在某一方向上受到的干扰。本发明是通过以下技术方案实现的，本发明包括以下步骤 第一步，选取窗口模板，并将窗口模板分解为若干个一维窗口； 所述的窗口模板是指因为本方法用于弱小目标检测，弱小目标在图像上往 往表现为一个多少数几个像素点，多个像素点又可以通过形态学滤波等方法将其 衰减为一个像素点，所以选择大小为3X3的矩形滤波窗口。所述的分解是指在3X3窗口中，分解是在垂直坐标系中以x轴正向、y 轴正向、y二x和y=_x共四个方向将窗口模板分为A、 A、 ￡>3和/)4四个定义域上 的1XN或NX1的矩阵，即一维窗口，这样分解的标准是使分解具有完备性和各 方向旋转不变性。第二步，设定待滤噪目标的先验一维分布特性，利用一维窗口在各个方向对 二维图片进行一维匹配滤波处理，所述的一维匹配滤波处理是指按照模板分解方向，用待滤噪目标的先验一维分布特性与待滤噪目标的待滤噪目标分布特性进行巻积运算,具体公式为
其中A = 0,1,L TV
噪目标的一维分布特性
先验一维分布特性，^^")为W^的离散傅立叶变换；^"为离散巻积结果。
所述的待滤噪目标的先验一维分布特性是指沿着第一步分解的一位窗口的
方向，设定目标的先验分布特性，比如若待检测目标点灰度比背景灰度打，则可
将先验一维分布特性设定为[4， 6， 4]，即中间元素值高于两边元素值。
第三步，在各个方向进行一维窗口匹配滤波后，定义优化设计的匹配滤波结
果为各个方向一维窗口匹配的峰值的均值，得到最终滤除噪声后的二维图片。所述的优化设计的匹配滤波结果为
W w't^,w ，
其中F(x，力表示坐标(x,力处匹配滤波结果，函数/(x,力表示原坐标(x,力处原始图像的分布特性，函数w^^ed()表示取值定义域中心点的值；g,(")表示目标沿z'方向分解的一维分布特性，其定义域为A; W")表示先验的待滤噪目标的一维分布特性，iV为分解的一维窗口的数目。
本发明的空间匹配滤波方法简单有效，它的关键在于充分考虑到噪声干扰的随机性和局部突变性，和噪声值只影响目标点某一维分布特性的特点。通过将二维滤波模板转化成多个一维滤波模板，结合图标先验分布特性，有效的检查出二维图片中的弱小目标点。
图l为实施例窗口分解示意其中图Ua)为3x3窗口模板；图1(b)为3x3窗口分解的一维窗口。图2为实施例一维匹配滤波特性示意其中图2(a)为先验目标分布，图2(b)为待滤噪目标分布，图2(c)为匹配结果。
=|>(力/<* _力=^S[G(")/Z(") Ww
-1， w为分解的一维窗口的数目，g")为A方向上待滤
，G(")为g("的离散傅立叶变换；W"为待滤噪目标的图3为实施例滤波前后图像空间分布特性示意其中图3(a)为滤波前图像空间分布特性，图3(b)为滤波后图像空间分布特性；
图4为实施例滤波前后图像及空间分布特性示意其中图4(a)为红外弱小目标图像，图4(b)为原始图像控件分布特性，图4(c)为经典二维空间匹配滤波器处理结果，图4 (d)为优化设计的空间匹配滤波器处理结果。
具体实施例方式
下面对本发明的实施例作详细说明，本实施例在以本发明技术方案为前提下进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限于下述的实施例。
本实例实施步骤具体如下
第一步、滤波模板选取
在此，主要考虑二维图片中弱小待滤噪目标成像的检测，弱小待滤噪目标在二维图片中仅占据一个像素点大小，通常选取3X3的窗口模板，其分解示意图如
l(b)所示。每个一维窗口是以当前像素为中心，令其位置为"，力，其灰度值为
/^，>0;令四个方向的定义域分别为Q、 D2、 A和A。
第二步、在模板分解方形进行一维滤波处理
优化设计的匹配滤波器在每一个一维窗口中，分别用先验的待滤噪目标的一维分布特性对其进行巻积，即进行一维匹配滤波处理
<formula>formula see original document page 7</formula>
式中，g")为沿着yt方向上待滤噪目标的分布特性，G(")为g("的离散傅
立叶变换；WW为先验目标的分布特性，根据本实例二维图片特征。目标点比周
围背景灰度值大，在形态上形成一个"突包"，取/<" = [464]， ^(")为W"的
离散傅立叶变换；^"为离散巻积结果。
图2(a)显示出了先验的空间待滤噪目标的轮廓截面图，目标出现位置局部区域内突变，形成一个类似于抛物线的形状，
图2(b)为实际二维图片中待滤噪目标空间分布进行分解得到的一个一维窗
口。为了消除光照强度引起的背景灰度值大小的整体影响，在进行巻积之前，截面轮廓图和分解的一维窗口分布均进行了 "归零化"处理，即将一维分布中的最小值归为零，只保留截面轮廓和分解窗口的"凸包"特性。
图2(c)为匹配结果，其中实线为先验目标与待滤噪目标的巻积结果，点划线为先验目标与背景的巻积结果，由此可见，巻积结果明显地分辨出了目标和背景的空间分布特性，准确地定位出了目标出现的位置。
第三步、在各个方向进行一维窗口匹配滤波后，定义优化设计的匹配滤波结果为各个方向一维窗口匹配的峰值的均值，得到最终滤除噪声后的二维图片，具体步骤为-
分别获取不同方向一维窗口巻积结果定义域中心点的峰值。定义优化设计的
匹配滤波结果为各个方向一维窗口匹配的峰值的均值，艮P:
<formula>formula see original document page 8</formula>
其中尸(X,力表示坐标(;C，力处匹配滤波结果，函数/(X，力表示原坐标(X,力处原始图像的分布特性，函数m^/尸ed()表示取值定义域中心点的值；g,(")表
示目标沿z'方向分解的一维分布特性，其定义域为A; Ww)表示先验的待滤噪目
标的一维分布特性，iV为分解的一维窗口的数目。
对整幅二维图片逐像素进行处理。滤波前后二维图片的空间分布特性对比
(见图3(a),图3(b))可以清楚地显示目标位置灰度值的增强。相比原始二维图片，经过优化设计的匹配滤波处理后目标信号大大增强。
为了验证本实施例的滤波性能，实验中采用法国Sofradir中波热成像仪采集到的320x240， 6733帧8位真实红外对空序列二维图片，截取其中的1000帧低信噪比二维图片验证。
选取1X1目标像素的红外二维图片，计算二维图片信噪比约为1.5，属于典型的低信噪比条件下的弱小目标检测问题。如4(a)为红外原始二维图片，其中白色方框标定的区域为目标出现位置，此时目标完全淹没在了杂波噪声中。4(b)为原始二维图片的空间分布特征图，可以看到在目标出现位置目标与背景灰度差异很小。4(c)为经过经典空间匹配滤波器处理后的二维图片的空间分布特性，可以看到因为原始二维图片遭噪声干扰的影响，经过匹配后目标点的位置并没有完全体现出优势，滤波后二维图片局部背景起伏特性加大，不是很平滑。4(d)为经过优化设计的空间匹配滤波后的二维图片的空间分布特性图，可以看到目标位置灰度值得到明显提升，目标出现位置像素值被提升到已经隐约可见。
对不同SVi 的待滤噪目标二维图片进行对比实验，定义信噪比
S潛,一w)/o" (3)
式中，s为目标信号幅度；"为二维图片背景灰度均值；o"为二维图片背景灰度标准差。实验结果如
表所示
表l匹配滤波前后二维图片信噪比对比表
原始二维图片经典二维匹配滤波后二维图优化匹配滤波后二维图片
SNR片SNRSNR
1. 52. 885. 9权利要求
1、一种弱小目标二维图片噪声滤除中的匹配滤波方法，其特征在于，包括以下步骤第一步，选取窗口模板，并将窗口模板分解为若干个一维窗口；第二步，设定待滤噪目标的先验一维分布特性，利用一维窗口在各个方向对二维图片进行一维匹配滤波处理；第三步，在各个方向进行一维窗口匹配滤波后，定义优化设计的匹配滤波结果为各个方向一维窗口匹配的峰值的均值，得到最终滤除噪声后的二维图片。
2、 根据权利要求1所述的弱小目标二维图片噪声滤除中的匹配滤波方法， 其特征是，所述的窗口模板是指大小为3X3像素 5X5像素的矩形窗口。
3、 根据权利要求1所述的弱小目标二维图片噪声滤除中的匹配滤波方法， 其特征是，所述的分解是指在垂直坐标系中以x轴正向、y轴正向、y^和y=-x 共四个方向将窗口模板分为A、 "2、 A和A四个定义域上的1XN或NX1的矩 阵，即一维窗口。
4、 根据权利要求1所述的弱小目标二维图片噪声滤除中的匹配滤波方法，其特征是，所述的一维匹配滤波处理是指按照模板分解方向，用待滤噪目标的先验一维分布特性与待滤噪目标的待滤噪目标分布特性进行巻积运算，具体公式为<formula>formula see original document page 2</formula>其中yt = 0,l，L iV-l， W为分解的一维窗口的数目，g("为A方向上待滤 噪目标的一维分布特性，G(")为g("的离散傅立叶变换；W^为待滤噪目标的先验一维分布特性，^(")为W"的离散傅立叶变换；W"为离散巻积结果。
5、 根据权利要求1所述的弱小目标二维图片噪声滤除中的匹配滤波方法， 其特征是，所述的待滤噪目标的先验一维分布特性是指沿着第一步分解的一位 窗口的方向，设定目标的先验分布特性。
6、 根据权利要求1所述的弱小目标二维图片噪声滤除中的匹配滤波方法，其特征是，所述的优化设计的匹配滤波结果为其中F(x，力表示……，函数wec/i^d()表示取值定义域中心点的值；表示z'方向分解的一维窗口，其定义域为"'；W")表示先验的待滤噪目标的一维 分布特性，W为分解的一维窗口的数目。
全文摘要
一种图像处理技术领域的弱小目标二维图片噪声滤除中的匹配滤波方法，包括选取窗口模板，并将窗口模板分解为若干个一维窗口；设定待滤噪目标的先验一维分布特性，利用一维窗口在各个方向对二维图片进行一维匹配滤波处理；在各个方向进行一维窗口匹配滤波后，定义优化设计的匹配滤波结果为各个方向一维窗口匹配的峰值的均值，得到最终滤除噪声后的二维图片。本发明充分考虑到噪声干扰的随机性和局部突变性，和噪声值只影响目标点某一维分布特性的特点。通过将二维滤波模板转化成多个一维滤波模板，结合图标先验分布特性，有效的检查出二维图片中的弱小目标点。
文档编号G06T5/00GK101630404SQ200910055610
公开日2010年1月20日 申请日期2009年7月30日 优先权日2009年7月30日
发明者坦 刘, 李元祥, 李建勋 申请人:上海交通大学