【技术实现步骤摘要】
本专利技术涉及电子信息技术雷达
，尤其涉及一种结合局部频率估计的小波域InSAR干涉相位滤波方法。
技术介绍
干涉合成孔径雷达(Interferometric Synthetic Aperture Radar，InSAR)是利用合成孔径雷达(SAR)两个通道的干涉相位信息提取地表的高程信息或变化信息，将SAR的测量拓展到三维空间，具有全天时、全天候、高精度的特点，因此在地形测绘、冰川研究、海洋测绘以及地面沉降监测等多个领域都有广泛的应用。干涉测量的精度和可靠性在很大程度上取决于干涉相位图的质量。然而，在实际系统中，受热噪声去相干、时间去相干、基线去相干、配准误差等多种去相干因素的影响，干涉相位图不可避免的存在相位噪声。低质量的干涉相位将会影响后续的干涉相位解缠及高程反演的准确性。因此，在相位解缠前必须对干涉相位进行滤波，从而获取高质量的干涉相位图。目前干涉相位的滤波方法可以大致分为空间域滤波和变换域滤波两类。圆周期均值滤波或中值滤波是一种...
【技术保护点】
一种结合局部频率估计的小波域InSAR干涉相位滤波方法，其特征在于，包括：步骤A：将InSAR的干涉相位φ变换到复数域ejφ，分别取复数域干涉相位ejφ的实部和虚部；步骤B：对干涉相位φ进行局部频率估计，得到干涉相位所在的频率范围；步骤C：对复数域干涉相位的实部和虚部，分别进行尺度为s的小波分解，得到不同子带的小波系数的频率范围，其中，m，n为小波系数的位置，i为分解的尺度，其范围为1～s；步骤D：对复数域干涉相位的实部和虚部，分别根据干涉相位所在的频率范围和不同子带小波系数的频率范围，分别确定有用信息和噪声所在的子带；步骤E：对复数域干涉相位的实部和虚部，分别对噪声所在子...
【技术特征摘要】
1.一种结合局部频率估计的小波域InSAR干涉相位滤波方法，其特
征在于，包括：
步骤A：将InSAR的干涉相位φ变换到复数域ejφ，分别取复数域干涉
相位ejφ的实部和虚部；
步骤B：对干涉相位φ进行局部频率估计，得到干涉相位所在的频率范
围；
步骤C：对复数域干涉相位的实部和虚部，分别进行尺度为s的小波分
解，得到不同子带的小波系数的频率范围，其中，m，n为小波系数的位
置，i为分解的尺度，其范围为1～s；
步骤D：对复数域干涉相位的实部和虚部，分别根据干涉相位所在的
频率范围和不同子带小波系数的频率范围，分别确定有用信息和噪声所在
的子带；
步骤E：对复数域干涉相位的实部和虚部，分别对噪声所在子带的小
波系数分别进行通用阈值收缩处理；
步骤F：对复数域干涉相位的实部和虚部，分别对有用信息所在子带
的小波系数分别进行邻域阈值收缩处理；
步骤G：对复数域干涉相位的实部和虚部，分别将通用阈值收缩处理
处理后的噪声所在子带的小波系数和邻域阈值收缩处理后的有用信息所
在子带的小波系数共同进行小波重构，得到滤波后的复数域干涉相位的实
部和虚部；以及
步骤H：由滤波后的复数域干涉相位的实部和虚部得到InSAR滤波后
的干涉相位。
2.根据权利要求1所述的小波域InSAR干涉相位滤波方法，其特征
在于，所述步骤B包括：
子步骤B1：对干涉相位φ的每一像素，对其进行局部频率估计，得到
其方位向干涉相位频率和距离向干涉相位频率子步骤B2：对所有像素估计得到的方位向干涉相位频率分别取最大
值和最小值为和距离向干涉相位频率分别取最大值和最小值为
和得到整幅干涉相位图所在的频率范围为3.根据权利要求2所述的小波域InSAR干涉相位滤波方法，其特征
在于，所述子步骤B1中，对干涉相位φ的每一像素得到其方位向干涉相位
频率和距离向干涉相位频率包括：
子分步骤B1a：以该像素为中心取(2M+1)×(2N+1)的估计窗口，该估
计窗口的相位模型为：
其中，为估计窗口中心点的相位，k、l为窗口内像素相对窗口中心
的位移，fa，fr分别为该窗口沿方位向和距离向的干涉相位频率；
子分步骤B1b：通过最大化估计窗口的相位模型的代价函数来估计该
像素沿方位向的干涉相位频率和距离向的干涉相位频率4.根据权利要求1所述的小波域InSAR干涉相位滤波方法，其特征
在于，所述步骤C中，不同子带小波系数的频率范围为：
asLL∈[0,2-sπ)×[0,2-sπ)]]>diL...
【专利技术属性】
技术研发人员：李芳芳，丁赤飚，雷斌，林雪，胡东辉，仇晓兰，
申请(专利权)人：中国科学院电子学研究所，
类型：发明
国别省市：北京;11
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。