冯丽源 聂祥飞 张 月 杨 鸿 梁 阳 陈天扬

(重庆三峡学院 电子与信息工程学院, 重庆 404000)

0 引言

合成孔径雷达干涉测量技术(Interferometric Synthetic Aperture Radar,InSAR)作为一种主动的、具有穿透性的微波遥感手段,能够全天时、全天候获取高精度的地球物理参数,如数字高程模型、地表形变和沉降、冰川运动等[1]。干涉相位是整个合成孔径雷达干涉处理中最重要的物理量[2],其质量将极大地影响获取地球物理参数的精度。然而,干涉相位受到各种去相关因素以及系统本身性能的影响,总是存在大量噪声,严重影响InSAR数据处理的后续步骤,因此需对其进行滤波进而得到更加精确的地表高程信息。

InSAR数据处理流程主要包含对合成孔径雷达SAR(Synthetic Aperture Radar)复图像的预滤波处理、影像的配准、干涉图的合成、去平地相位、干涉相位的滤波、干涉相位的解缠、绝对相位的获取和数字高程模型DEM(Digital Elevation Model)重建等步骤。其中,干涉相位滤波技术是对InSAR进行数据处理过程中一个非常重要的环节[3],经滤波后的干涉相位图质量将决定地面形变量的精度和数据处理的准确度。常见的雷达相位图滤波处理方式有很多,例如张斌等[2]人在2008年提出了一种改进的形态学干涉图滤波方法,降低了相位图噪声;郭交等[4]人提出了一种新的干涉相位图局部自适应滤波方法,使干涉相位图的边缘得到了很好的保持;2005年,汪鲁才[5]提出了一种基于小波变换和中值滤波的InSAR干涉图像滤波方法,不仅对相位图突变处的处理更加平滑,而且去除了相位图中的大量残差点。……