用小波变换和中值滤波研究差分干涉图的去噪 a study on the differential interogram by wavelet transform and median filtering.pdfVIP

用小波变换和中值滤波研究 差分干涉图的去噪 余景波1’2，刘国林1’2，王建波1，葛振坦1’2 (1．山东科技大学测绘科学与工程学院．山东青岛266510； 2．海岛(礁)测绘技术国家测绘局重点实验室，山东青岛266510) 摘 要：介绍了小波变换的基本原理和图像去噪常见的滤波方法，采用几种常见滤波分 别对模拟差分干涉图和EVISAT卫星获取的矿区真实合成孔径雷达(ASAR)数据的差分干 涉图分别进行滤波去噪处理，并对其去噪效果进行分析。采用小波变换和中值滤波相结合的 方法对矿区真实ASAR数据差分干涉图进行去噪处理，并对先中值滤波再d、波变换和先小波 变换再中值滤波两种方式去噪结果分别进行了分析比较，结果表明：先小波变换再中值滤波去 噪后，图像保真效果较好。 关键词：小波变换；中值滤波；模拟干涉图；ENVISAT卫星；矿区真实ASAR数据差分干 涉图；去噪 中图分类号：TP751文献标志码：A 文章编号：1008—9268(2011)02—0029-07 方法对矿区真实ASAR数据差分干涉图进行滤波 0 引 言 处理。 先进的合成孔径雷达(ASAR)数据差分干涉 图在平地效应消除后，由于受到噪声的污染，从而 1小波变换的基本原理 存在大量残差点，给相位解缠带来了难度，所以，对 小波变换是一种时间窗和频率窗都可以改变 ASAR数据差分干涉图进行滤波处理是十分必要 的时频局部分析方法[73。设妒(￡)∈L2(R)表示 的。差分干涉图的噪声主要来源于以下几个方 能量有限的信号空间。其傅里叶变换为驴(∞)，当 面D-s]：InSAR系统本身的热噪声；相位图的相干 斑点噪声；基线失相干和时间失相干等。许多学者 其满足允许条件时q；f哗如∞，称 JR I叫l 给出了不同对差分干涉图进行滤波处理的方法。 ≯(￡)为一个基本小波。将其经伸缩和平移后，就 Lee等人[4]提出了应用加性相位噪声模型去除干 可以得到一个小波序列。 涉图相位图噪声的方法，该方法能够较好的保持相 位条纹的连续性，但是计算工作量极大。Eichel等 人[5]提出了圆周均值滤波方法，它在理论上具有最 为[83 大似然最优的滤波效果，但在保持相位条纹连续性 W，(口，6)=(，(￡)，咖，6(￡) 上较差。Lanari和Fornaro等人[6]提出的圆周中 值滤波方法，能够保持相位条纹连续性，由于没有 使用信号统计规律，其滤波效果不是最好的。 ，(t)=古．『矿，R嘉w，(口，6)9(气尘)dadb 本文采用常见的滤波对模拟干涉图和矿区 ASAR数据差分干涉图进行滤波去噪处理，比较 (2) 其滤波效果；最后，使用小波变换和中值滤波结合 式(2)称为小波的逆变换。 收稿日期：2010—12-07 基金项目：海岛(礁)测绘技术国家测绘局重点实验室资助项目(2010A01) 联系人：余景波E-maill918ky(孕sina．corn ·29· 万方数据 式(1)中，缈(￡)为雪(f)的共轭函数，若 2去噪的方法 9(￡)为