用小波变换研究ASAR差分干涉图的去噪(第三次修改稿).docVIP

用小波变换和中值滤波研究差分干涉图的去噪 余景波1，2 刘国林1,2 王建波1 曹振坦1,2 （ 1．山东科技大学 测绘科学与工程学院，山东 青岛 266510； 2. 海岛（礁）测绘技术国家测绘局重点实验室，山东 青岛 266510 介绍了小波变换的基本原理和图像去噪常见的滤波方法，采用几种常见滤波分别对模拟差分干涉图和EVISAT卫星获取的矿区真实ASAR数据的差分干涉图分别进行滤波去噪处理，并对其去噪效果分别进行分析。最后，采用用小波变换和中值滤波相结合的方法对矿区真实ASAR数据差分干涉图进行去噪处理，并对先中值滤波再小波变换和先小波变换再中值滤波两种方式去噪结果进行了分析比较，结果表明先小波变换再中值滤波去噪效果较好。 小波变换；中值滤波；模拟干涉图；ENVISAT卫星；矿区真实ASAR数据差分干涉图；去噪 A study on the differential interogram by wavelet transform and median filtering YU Jing-bo1,2 LIU Guo-lin1,2 WANG Jian-bo1 CAO Zhen-tan1,2 （1.Geomatics College, Shandong University of Science and Technology, Qingdao 266510, China; 2. Key Laboratory of Surveying and Mapping Technology on Island and Reef, State Bureau of Surveying and Mapping, Qingdao 266510, China） Abstract: The basic principles of wavelet transform and the methods of image denoising filter were introduced. Their effects of denoising were compared after several common filtering denoising the simulated Gaussian surfaces and the real mining ASAR data from ENVISAT satellite differential interferograms. Finally，the wave transform and median filtering effect was better by the combination of wavelet transform and median filtering denoising the mining ASAR data differential interferograms. Keywords: Wavelet transform；median filter；simulated Gaussian surface；ENVISAT satellite；the mining real ASAR data differential interferograms；denoise 先进的合成孔径雷达（ASAR）数据差分干涉图在平地效应消除后，由于受到噪声的污染，从而存在大量残差点，给相位解缠带来了难度，所以对ASAR数据差分干涉图进行滤波处理是十分必要的。差分干涉图的噪声主要来源于以下几个方面[1-3]：InSAR系统本身的热噪声；相位图的相干斑点噪声；基线失相干和时间失相干等。许多学者给出了不同对差分干涉图进行滤波处理的方法。J.S.Lee等人[4]提出了应用加性相位噪声模型去除干涉图相位图噪声的方法，该方法能够较好的保持相位条纹的连续性，但是计算工作量极大。P.H.Eichel等人[5]提出了圆周均值滤波方法，它在理论上具有最大似然最优的滤波效果，但在保持相位条纹连续性上较差。R.Lanari和G.Fornaro等人[6]提出的圆周中值滤波方法，能够保持相位条纹连续性，由于没有使用信号统计规律，其滤波效果不是最好的。 本文采用常见的滤波对模拟干涉图和矿区ASAR数据差分干涉图进行滤波去噪处理，比较其滤波效果；最后，使用小波变换和中值滤波结合方法对矿区真实ASAR数据差分干涉图进行滤波处理。 1小波变换是一种时间窗和频率窗都可以改变的时频局部分析方法[7]。设表示能量有限的信号空间。其傅里叶变换为，当其满足允许条件时：，称为一个基本小波。将其经伸缩和平移后，就可以得到一个小波序列。 函数的小波变换（CWT）定义为[8]：