一种基于经验模式分解的InSAR干涉相位滤波方法技术

技术编号:8980367 阅读:159 留言:0更新日期:2013-07-31 22:28
本发明专利技术提供了一种基于经验模式分解的InSAR干涉相位滤波的方法,克服传统的基于窗口的处理方法不能兼顾去噪和细节保持的缺点。本发明专利技术利用二维经验模式分解将干涉相位复信号分解为不同尺度的分量,各个层次的分量根据其包含干涉相位和噪声的特征不同程度的提高其指数,使得每一个分量的噪声得到抑制并且不破坏干涉相位的结构,从而克服了传统的基于窗口的滤波方法不能兼顾去噪和细节保持的缺点,实现高精度的干涉相位滤波,为高精度的干涉测量提供了条件。

【技术实现步骤摘要】

本专利技术属于干涉合成孔径雷达信号处理领域,尤其涉及一种基于经验模式分解的干涉合成孔径雷达干涉相位滤波的方法。
技术介绍
干涉合成孔径雷达(InterferometricSynthetic Aperture Radar, InSAR)是利用合成孔径雷达(SAR)复数据的相位信息获取地表的高程信息或变化信息的一项技术,具有全天时、全天候、高精度的特点,因此在地形测绘、冰川研究、海洋测绘以及地面沉降监测等多个领域都有广泛的应用。干涉测量的精度和可靠性在很大程度上取决于干涉相位图的质量。然而,在实际系统中,受热噪声去相干、时间去相干、基线去相干、配准误差等多种去相干因素的影响,干涉相位图不可避免的存在相位噪声。干涉相位噪声的存在直接影响相位解缠的效果及最终干涉测量的精度。因此,在相位解缠前必须对干涉相位进行滤波,从而获得较为准确的干涉相位估计值。目前干涉相位的滤波方法可以大致分为空间域滤波和频率域滤波。多视滤波是一种最常用的空间域滤波方法(R.Lanar1.Generation of digital elevation models byusing SIR-C/X-SAR multifrequency two-pass interferometry: The Etna case study.1EEE Transactions on Geoscience and Remote Sensing, 1996, 34(5):1097-1114.),它实现简单,运算速度快,但滤波视数难以确定,在条纹密集时容易破坏相位细节,降低分辨率。频率域滤波方法中,Goldstein滤波(R.M.Goldstein, C.L.Werner.RadarInterferogram filtering for Geophysical Application.Geophysical ResearchLetters.1998,25(21):4035-4038.)应用最为广泛,但该方法受分块大小和滤波参数的影响较大,在信噪比很低时,滤波效果较差。因此,为满足InSAR应用对干涉相位精度的要求,需要研究能有效去除噪声并保持相位细节的自适应滤波方法。
技术实现思路
为解决上述问题,本专利技术提供了一种基于经验模式分解的InSAR干涉相位滤波的方法,以克服传统的基于窗口的处理方法不能兼顾去噪和细节保持的缺点,从而提高干涉相位的精度。该基于经验模式分解的InSAR干涉相位滤波的方法,包括以下步骤:第一步:将干涉相位Φ变换到复数域中为ZeK=Cos Φ+jsin Φ ,分别取的实部为 Re {ej4>} =Cos Φ,虚部为 Im{eJ<i>} =sin Φ ;第二步:对的实部RelVI进行二维经验模式分解,得到各层内蕴模式分量Fn,n=l,...,N和一个剩余分量Rn,其中N为分解层数;第三步:分别对Fn,n=l,…,N和Rn进行小波分解,得到各层内蕴模式分量的小波系数和剩余分量的小波系数C;1,其中,η为内蕴模式分量的层数,j为小波尺度,k为小波系数的位置,N+1层代表剩余分量;第四步:根据各层内蕴模式分量和剩余分量所包含信号和噪声的特征,分别设置HSlder指数增量Λ a n,n=l,...,N+1,其中N+1层代表剩余分量;第五步:根据小波系数与H01der指数改变量的关系,分别对各层内蕴模式分量和剩余分量的小波系数(,进行调整,得到调整后的小波系数^,第六步:根据各层内蕴模式分量和剩余分量调整后的小波系数%,分别进行小波重构,得到调整后的各层内蕴模式分量和剩余分量第七步:将调整后的内蕴模式分量和剩余分量进行合成,从而得到滤波后的复干涉相位的实部,即本文档来自技高网
...

【技术保护点】
基于经验模式分解的InSAR干涉相位滤波的方法,其特征在于,包括以下步骤:第一步:将干涉相位φ变换到复数城中为:ejφ,cosφ+jsinφ,分别取ejφ的实部为Re{ejφ)=cosφ,虚部为Im{ejφ}=sinφ;第二步:对ejφ的实部Re{ejφ}进行二维经验模式分解,得到各层内蕴模式分量Fn,n=1,…,N和一个剩余分量RN,其中N为分解层数;第三步:分别对Fn,n=1,…,N和RN进行小波分解,得到各层内蕴模式分量的小波系数和剩余分量的小波系数其中,n为内蕴模式分量的层数,j为小波尺度,k为小波系数的位置;第四步:根据各层内蕴模式分量和剩余分量所包含信号和噪声的特征,分别设置指数增量Aan,n=1,…,N+l,其中N+l层代表剩余分量;第五步:根据小波系数与指数改变量的关系,分别对各层内蕴模式分量和剩余分量的小波系数进行调整,得到调整后的小波系数第六步:根据各层内蕴模式分量和剩余分量调整后的小波系数分别进行小波重构,得到调整后的各层内蕴模式分量和剩余分量第七步:将调整后的内蕴模式分量和剩余分量进行合成,从而得到滤波后的复干涉相位的实部,即第八步:对ejφ的虚部Im{ejφ}也重复第二步至第七步,得到滤波后复干涉相位的虚部Im;第九步:根据滤波后的复干涉相位的实部和虚部,计算滤波后的干涉相位为FDA00002964374000011.jpg,FDA00002964374000012.jpg,FDA00002964374000013.jpg,FDA00002964374000014.jpg,FDA00002964374000015.jpg,FDA00002964374000016.jpg,FDA00002964374000017.jpg,FDA00002964374000018.jpg,FDA00002964374000019.jpg,FDA000029643740000115.jpg,FDA000029643740000111.jpg...

【技术特征摘要】
1.基于经验模式分解的InSAR干涉相位滤波的方法,其特征在于,包括以下步骤: 第一步:将干涉相位Φ变换到复数城中为(308(^+]^11(^,分别取&_<*>的实部为Re {ej4>) =Cos Φ ,虚部为 Im{eJ<i>} =sin Φ ; 第二步:对^1*.的实部RelVI进行二维经验模式分解,得到各层内蕴模式分量Fn,n=l,…,N和一个剩余分量Rn,其中N为分解层数; 第三步:分别对Fn,n=l,…,N和Rn进行小波分解,得到各层内蕴模式分量的小波系数‘和剩余分量的小波系数^ 1,其中,η为内蕴模式分量的层数,j为小波尺度,k为小波系数的位置; 第四步:根据各层内蕴模式分量和剩余分量所包含信号和噪声的特征,分别设置Klde...

【专利技术属性】
技术研发人员:李芳芳胡东辉仇晓兰丁赤飚
申请(专利权)人:中国科学院电子学研究所
类型:发明
国别省市:

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1