雷达干涉测量图像的梯度自适应光滑子函数滤波法

付政庆,刘国林,陶秋香,刘伟科

1．山东科技大学理学院,山东青岛 266590;2．山东科技大学测绘学院,山东青岛 266590

雷达干涉测量图像的梯度自适应光滑子函数滤波法

付政庆1,刘国林2,陶秋香2,刘伟科2

1．山东科技大学理学院,山东青岛 266590;2．山东科技大学测绘学院,山东青岛 266590

提出基于梯度的自适应光滑子函数滤波的新方法。通过相位梯度判断出干涉图中的噪声点。在噪声点处,将包含噪声点的InSAR图像窗口拟合为二次曲面,然后用光滑子函数对拟合曲面进行平滑去噪处理。利用北京地区的真实ALOS数据进行滤波试验,并与几种常用的滤波方法进行比较和分析。通过均方根误差、峰值信噪比、边缘保持指数和标准相位偏差等量化指标,验证了基于梯度的自适应光滑子函数滤波方法的可行性与有效性。此方法可有效抑制斑点噪声,较好地保留边缘信息并且干涉图相位偏差较少。

雷达干涉测量;自适应;光滑子函数;滤波;斑点噪声

1 引 言

InSAR图像由于成像方式的原因,不可避免地会产生噪声,而大量的颗粒状的相干斑降低了图像质量,影响了目标的检测、分类和识别,因此噪声去除对图像后期处理非常重要[1-7]。目前常用的降噪方法有多视滤波、自适应均值滤波、自适应中值滤波、自适应小波滤波和自适应圆周期中值滤波等[8-9]。如多视滤波法可以抑制噪声,但是这种方法同时平滑了影像数据,是以损失图像空间分辨率为代价的,自适应均值滤波、自适应中值滤波和自适应圆周期中值滤波是传统滤波方法的改进,在去除噪声的同时可以保留更多的图像细节。还有一些滤波方法,比如粒子滤波[10-12],或基于某种量化指标,如文献[13]是基于信噪比的滤波方法,文献[14]是基于特征保持的滤波方法。

本文提出利用泛函分析中的光滑子函数进行滤波,这是因为光滑子函数具有保凸性和光滑性等很好的性质[15-16],可以用来对函数进行磨光处理。在其与曲面相结合时,可将带有角点的函数磨光而又与原来的函数相切,在图像发生突变处使图像变平滑,产生“峰削尖、谷填底”的效果。

将InSAR图像的相位考虑为空间中的曲面,而斑点噪声使曲面发生突变,从而利用光滑子函数的性质对其进行滤波。对每一个点求出其4个方向的梯度,利用这4个梯度来判断这个点是否为噪声点。如果是噪声点,将噪声点所在的窗口进行曲面拟合,然后用光滑子函数对拟合曲面进行平滑处理;如果不是噪声点,则不进行平滑处理。通过这种自适应滤波方法,可以较好地抑制斑点噪声同时,又较好地保留了边缘信息并使相位偏移量较小。本文采用真实数据进行了滤波试验,并将自适应光滑子函数滤波和自适应圆周期中值滤波、自适应均值滤波、自适应中值滤波、自适应小波滤波的效果进行了分析和比较。

2 光滑子函数滤波原理

2．1 曲面构造

式中,h为图像中心到边缘的距离。

如果没有噪声,干涉图的相位值具有很好的连续性。但是斑点噪声的存在破坏了这种连续性,滤波的目的就是恢复这种连续性。噪声的存在,使噪声点与周围的相位之间产生较大的梯度,因此可以利用梯度来判断干涉图中的突变点。

基于梯度的自适应法就是由干涉图的梯度来区分噪声点和非噪声点,进而决定是否要进行滤波处理。图像当中每一点的梯度为向量,本文中求出每一点梯度的模,并将梯度的模称之为图像在这一点的梯度。将带有噪声的干涉图进行分割,对于中心点,其梯度按照方向来分有4种,分别是:水平方向、左上右下方向、垂直方向、右上左下方向。公式如下

对于中心点是否为噪声,采用这个点的上述4个方向梯度进行判断[18]。具体过程是,设置阈值k,求出4个梯度的最小值,如果最小值大于k,说明4个方向梯度都比较大,则认为该点为噪声点,将这个小的曲面进行拟合,并与二维光滑子函数相结合,进行平滑;如果最小值不大于k,可以认为不是噪声点,不参与平滑,从而得到噪声去除后的干涉图。数据处理基本流程如图1所示。

图1 数据处理流程Fig．1 Data flow processing

3 定量评价指标

对InSAR干涉图用4种滤波方法进行滤波去噪处理,通过边缘保持指数、峰值信噪比、均方根误差、相位标准偏差和不连续点等量化指标来对图像质量进行比较和评价[19-21]。

3．1 边缘保持指数

边缘保持指数(edge preserve index,EPI)反映了对图像进行滤波后的边缘信息保持能力,其值越大,表明滤波器的边缘保持能力越强。其表

达式为

式中,ϕsi(,j)是滤波后的相位值;ϕoi(,j)是原始干涉图相位值。

3．2 峰值信噪比

峰值信噪比(peak signal-to-noise ratio, PSNR)是反映图像质量好坏的较为通用的指标。其值越大,去噪后的图像质量越高。其表达式为

3．3 均方根误差

均方根误差(root mean square error,RMS)用来衡量滤波后相位图和参考相位图产生的偏差,其值越小代表滤波器的保真性越好。其表达式为

式中,ϕsi(,j)是滤波后的相位值;ϕoi(,j)是参考干涉图相位值;N是滑动窗口像元个数。

3．4 相位标准偏差

相位标准差(phase standard deviation, PSD)反映了相位的平滑程度,其值越小,说明相位越平滑。其计算公式为式中,φ(i,j)是干涉图相位值(i,j)是所选滑动窗口干涉图相位值的平均值。

3．5 不连续点个数

噪声去除后,干涉图像上的不连续点个数也是噪声去除能力的重要指标之一。通过滤波使图像不连续点减少到可接受的水平,为下一步相位解缠准备。

4 试验结果与分析

4．1 试验结果

为了对光滑子函数滤波和现有滤波方法进行比较分析,本文采用了时间为2009-01-10和2009-02-25的北京ALOS数据,从中选取了400像素×400像素进行了试验。图2给出了原始干涉图及滤波后的图像,其中原始干涉图中,梯度的最小值为2．050×10－8,滤波时阈值k取为5．6× 10－5,表1给出了各个降噪算法滤波结果的量化指标。

图2 原始干涉图的降噪结果Fig．2 Filter results of original interferogram

表1 原始干涉图降噪结果的量化指标Tab．1 Quantitative index of filter in original interferogram

在原始干涉图和滤波后的干涉图中,选取第200行的相位值进行比较研究,其剖面图如图3所示。

4．2 试验结果分析

从图2(a)中可以看到,图像条纹比较清晰,但条纹的边缘仍有很多斑点噪声。经过光滑子函数滤波和其他4种方法滤波后,图像干涉条纹清晰,斑点噪声明显减少,说明这几种方法都可以有效地去除斑点噪声,同时又保留了图像的细节特征。

从图3的剖面图可以看出,与其他滤波方法相比,经本文提出的滤波方法滤波后,干涉图的条纹曲线变化平稳且毛刺较少。这说明光滑子函数滤波方法不但能去除大部分噪声,而且对于较大颗粒的斑点噪声效果也非常好。

图3 模拟的和滤波后的干涉图剖面(第200行)Fig．3 Cross section over the simulated and filtered interferogram(row 200)

通过表1中的量化指标可以对这4种滤波方法进行定量比较。在这5种滤波方法中,光滑子函数滤波方法的RMS最小,这说明用这种方法滤波后的相位图和参考相位图产生的偏差最小,即滤波器的保真性是这几种方法中最好的。从PSNR来看,经过自适应光滑子函数滤波后, PSNR要略低于小波滤波但是高于其他3种方法,说明这种方法在去噪后图像质量方面比小波滤波稍弱但是要好于其他3种滤波方法。EPI是指滤波器保持边缘信息的能力,从数据来看本文提出的方法在滤波后低于自适应圆周期中值滤波,高于其他3种滤波方法,这说明光滑子函数滤波在保持边缘信息方面的效果不如自适应圆周期中值滤波,但是要好于其他3种滤波方法。经过5种滤波方法去噪后,PSD均有下降,又以自适应光滑子函数滤波下降幅度最大,比原始图像减少了24．7%,自适应圆周期中值滤波降低很少,只减少了1．7%。这说明滤波后相位平滑程度均有所增加,但本文方法的相位平滑能力较强。

对干涉图去除噪声后,干涉图上的不连续点个数也是衡量滤波的重要指标值之一。不连续点个数在一定程度上反映了去噪后干涉图上的相位值高的颗粒噪声的数量多少。滤波后不连续点个数越少说明滤波方法抗相位畸变能力越强。表2给出了原始及滤波去噪后,真实数据干涉图上断点的个数和通过滤波不连续点减少比例。

表2 原始及滤波后干涉图不连续点个数Tab．2 The number of breakpoints of original and denoised interferogram

从表2可以看出通过5种方法滤波后,不连续点都明显减少,自适应光滑子函数滤波和自适应均值滤波后不连续点最少,要明显少于小波滤波、自适应圆周期中值滤波和自适应中值滤波后的不连续点。这说明自适应光滑子函数滤波方法抗相位畸变的能力和自适应均值滤波相当,又要明显强于其他3种滤波方法。图4给出原始干涉图及滤波后干涉图不连续点的分布图。

5 结 论

本文提出的基于梯度的光滑子函数进行滤波的方法是一种新的滤波方法。首先通过每个点的梯度来判断是否为噪声,如果是噪声,通过将包含噪声点的窗口进行曲面拟合,然后将拟合曲面与光滑子函数进行结合来平滑滤波,如果不是噪声则不进行平滑。在一定程度上消除了斑点噪声对图像造成的突变,从而达到了抑制斑点噪声的目的。利用真实数据进行试验,并与几种常用的滤波方法比较分析来看,本文的滤波方法在图像的保真性和相位平滑能力这几方面都是有优势的,在滤波后的图像质量方面好于自适应均值、中值及圆周中值滤波,在边缘保持能力方面要好于自适应中值、均值及小波滤波,抗相位畸变方面效果也较好。另外,光滑子函数还具有对空间目标进行平滑的作用,对以后的目标检测,分类和识别也具有一定作用。

图4 原始和滤波后干涉图不连续点分布图(白色为不连续点)Fig．4 The original and denoised interferogram breakpoint distribution

[1] WANG A L,ZHANG Y,GU Y F．Simultaneous Speckle Reduction and SAR Image Compression Using Multiwavelet Transform[J]．Journal of Electronic Science and Technology of China,2007,5(2):163-166．

[2] FERNANDEZ S A,LOPEZ C A．On the Estimation of the Coefficient of Variation for Anisotropic Diffusion Speckle Filtering[J]．IEEE Transactions on Image Processing,2006,15(9):2694-2701．

[3] GOLDSTEIN T,BRESSON X,OSHER S．Geometric Applications of the Split Bregman Method:Segmentation and Surface Reconstruction[J]．Journal of Scientific Computing,2010,10(1):272-293．

[4] PI Yiming,YANG Jianyu,FU Yusheng,et al．The Principle of Sythetic Aperture Radar Imaging[M]．Chengdu:Press of University of Electronic Science and Technology,2007．(皮亦鸣,杨建宇,付毓生,等．合成孔径雷达成像原理[M]．成都:电子科技大学出版社,2007)．

[5] MAITRE H．Traitement des Images de RadaràSynthèse dOuyerture[M]．SUN Hong,translate．Beijing: Publishing House of Electronics Industry,2005．(MAITRE H．合成孔径雷达图像处理[M]．孙洪,译．北京:电子工业出版社,2005．)

[6] LIAO Mingsheng,LIN Hui,ZHANG Zuxun,et al．Adaptive A1gorithm for Filtering Interferometric Phase Noise[J]．Journal of Remote Sensing,2003,7(2):98-105．(廖明生,林晖,张祖勋,等．InSAR干涉条纹图的复数空间自适应滤波[J]．遥感学报,2003,7(2):98-105．)

[7] HAN Chunming,GUO Huadong,WANG Changlin．The Essence of SAR Image Speckle Suppression[J]．Journal of Remote Sensing,2002,6(6):470-474．(韩春明,郭华东,王长林．SAR图像斑点噪声抑制的本质[J]．遥感学报, 2002,6(6):470-474．)

[8] LIAO M S,LIN H,ZHANG Z X,et al．Adaptive Algorithm for Filtering Interferometric Phase Noise[J]．Journal of Remote Sensing,2003,7(2):98-105．

[9] WANG Xingwang．The Filter Method of InSAR Phase[D]．Changsha:Central South University,2008．(王兴旺．合成孔径雷达干涉图相位滤波方法的研究[D]．长沙:中南大学,2008．)

[10] TIAN Shijun,CHEN Jun,PI Yiming．Particle Filtering and Its Application for Kinematic GPS Positioning with High Speed Movement[J]．Acta Geodaetica et Cartographica Sinica,2008,36(3):274-278．(田世君,陈俊,皮亦鸣．粒子滤波在高动态GPS定位中的应用[J]．测绘学报,2008,36(3):274-278．)

[11] NIE Jianliang,YANG Yuanxi,WU Fumei．An Algorithm of Dynamic Precise Point Positioning Based on Modified Particle Filtering[J]．Acta Geodaetica et Cartographica Sinica,2010,39(4):338-343．(聂建亮,杨元喜,吴富梅．一种基于改进粒子滤波的动态精密单点定位算法[J]．测绘学报,2010,39(4):338-343．)

[12] GONG Yisong,GUI Qingming,LI Baoli,et al．Design of Particle Filtering Algorithm Based on Mean Shift and Its Application in Navigation Data Processing[J]．Acta Geodaetica et Cartographica Sinica,2011,40 (supplement):120-126．(宫轶松,归庆明,李保利,等．基于均值漂移的粒子滤波算法设计及其在导航数据处理中的应用[J]．测绘学报,2011,40(增刊):120-126．)

[13] SUN Qian,ZHU Jianjun,LI Zhiwei,et al．A New Adaptive InSAR Interferogram Filter Based on SNR[J]．Acta Geodaetica et Cartographica Sinica,2009,38(5):437-443．(孙倩,朱建军,李志伟,等．基于信噪比的InSAR干涉图自适应滤波[J]．测绘学报,2009,38(5):437-443．)

[14] YANG Shenbin,LI Bingbai,SHEN Shuanghe,et al．Structure Retaining Linear Multi Channel SAR Image Speckle Filter [J]．Acta Geodaetica et Cartographica Sinica,2006,35 (4):364-370．(杨沈斌,李秉柏,申双和,等．基于特征保持的线性多通道最优求和SAR图像滤波算法[J]．测绘学报,2006,35(4):364-370．)[15] YAO Zongyuan,QIU Chunhui,ZHONG Chunping．A New Leray Formula for Smooth Functions on Bounded Domains in Cn[J]．Science in China:Series A,2002,45 (5):610-619．

[16] MAZROUI A,SBIBIH D,TIJINI A．A Simple Method for Smoothing Functions and Compressing Hermite Data [J]．Advances in Computational Mathematics,2005,10 (3):279-297．

[17] ZHANG Mengjun,SHU Hong,LIU Yan,et al．An Adaptive Thresholding Approach Based on Spatial Curved Surface Fitting[J]．Geomatics and Information Science of Wuhan University,2006,31(5):395-398．(张孟军,舒红,刘艳,等．基于空间曲面拟合的自适应阈值选取方法[J]．武汉大学学报:信息科学版,2006,31(5):395-398．)

[18] BOVIK A C．On Detecting Edges in Speckle Imagers[J]．IEEE Transactions on Signal Processing,1988,36(10): 1618-1627．

[19] GHIGLIA D C,MASTIN G A,ROMERO L A．Cellularautom Atemethod for Phase Unwrapping[J]．Journal of the Optical Society of America,1987,4(1):267-280．

[20] LIU Guoxiang．Monitoring of Ground Deformations with Radar Interferometry[M]．Beijing:Surveying and Mapping Press,2006．

[21] GOLDSTEIN R M,WERNER C L．Radar Interferogram Filtering for Geophysical Applications[J]．Geophysical Research Letters,1998,25(21):4035-4038．

(责任编辑:丛树平)

Adaptive Smoothing Function Filtering Based on Gradient in InSAR Image

FU Zhengqing1,LIU Guolin2,TAO Qiuxiang2,LIU Weike2
1．College of Science,Shandong University of Science and Technology,Qingdao 266590,China;2．Geomatics College,Shandong University of Science and Technology,Qingdao 266590,China

A new adaptive filtering algorithm using smoothing function based on phase gradient is proposed．The noise points in interferogram can be judged out through gradient of phases．At the points InSAR image window is fitted to quadric surface,which can be processed by smoothing function．This method is compared with other common filtering methods using the genuine ALOS data of Beijing area and the feasibility and effectiveness of which are verified by quantitative indicators such as RMS,PSNR,EPI and PSD,which suggest that it reduces the speckle noises effectively and keep the edge information well．The phase deviation of interference in figure is less．

InSAR;adaptive;smoothing function;filter;speckle noise

FU Zhengqing(1978—),male,PhD, lecturer,majors in surveying data processing．

LIU Guolin

P237

A

1001-1595(2014)03-0263-05

国家自然科学基金(41274007);高等学校博士学科点专项科研基金(博导类) (20123718110001);山东省自然科学基金(ZR2012DM001);山东科技大学创新团队基金(2011KYTD103);山东省优秀中青年科学家科研奖励基金(BS2013SF013)

2012-11-26

付政庆(1978—),男,博士,讲师,研究方向为测量数据处理。

E-mail:fzhqing＠163．com

刘国林

修回日期:2013-06-07