申请/专利权人：中国矿业大学;中国自然资源航空物探遥感中心

申请日：2021-11-30

公开（公告）日：2024-06-28

公开（公告）号：CN114066778B

主分类号：G06T5/70

分类号：G06T5/70;G06T5/50

优先权：

专利状态码：有效-授权

法律状态：2024.06.28#授权;2022.03.08#实质审查的生效;2022.02.18#公开

摘要：本发明公开了一种基于拉普拉斯金字塔的InSAR干涉图迭代自适应滤波方法，适用于InSAR干涉图图像处理。使用迭代自适应滤波构建拉普拉斯金字塔，对拉普拉斯金字塔的每层输入影像都采用迭代自适应滤波，迭代自适应滤波首先对输入图像进行局部频率估计，保留图像中主频部分，然后对图像剩余部分进行基于相干性的Goldstein滤波，自适应确定相关参数，最后获得每层输入图像的迭代自适应滤波结果；将下一层下采样图像和本层上采样图像相减获得的差异图进行迭代自适应滤波并存入拉普拉斯金字塔，最后对拉普拉斯金字塔进行重构获得最终的图像滤波结果，同时其易于实现，干涉图滤波的效果较好，能够为后续的相位解缠等步骤提供良好的滤波结果和数据支撑。

主权项：1.一种基于拉普拉斯金字塔的InSAR干涉图迭代自适应滤波方法，其特征在于：使用迭代自适应滤波构建拉普拉斯金字塔，对拉普拉斯金字塔的每层输入图像都采用迭代自适应滤波，迭代自适应滤波首先对输入图像进行局部频率估计，保留图像中主频部分，然后对图像剩余部分进行基于相干性的Goldstein滤波，自适应确定相关参数，最后获得每层输入图像的迭代自适应滤波结果；将下一层下采样图像和本层上采样图像相减获得的差异图进行迭代自适应滤波并存入拉普拉斯金字塔，最后对拉普拉斯金字塔进行重构获得最终的图像滤波结果；具体步骤如下：S1：获取一幅w行h列的InSAR干涉图I0，该图像为复数图像，取干涉图I0相位值进行结果展示；S2：设置迭代自适应滤波的窗口最大值为wmax，最小值为wmin；S3：建立拉普拉斯金字塔，并根据InSAR干涉图I0的行w列h，设置基于拉普拉斯金字塔架构的总层数为J；S4：将干涉图I0作为输入图像输入拉普拉斯金字塔的第0层进行处理，从而将处理后的图像作为拉普拉斯金字塔高1层的输入，直至到金字塔的J-1层；处理过程具体为：利用迭代自适应滤波对金字塔第n0≤n≤J-1层的输入图像In进行处理获得滤波结果In′,然后将滤波结果In′隔行隔列下采样获得In′的下采样图像Lnd，下采样图像Lnd的大小变为下采样前图像In′的14，利用双线性插值算法计算下采样图像Lnd的上采样结果Lnp，将拉普拉斯金字塔第n-1层图像In-1第0层采用原始图像I0与上采样的结果Lnp相减，即可获得第n层金字塔分解结果Ln，对Ln进行迭代自适应滤波获得滤波结果Ln′并存入金字塔第n层中；最后将第n层下采样后的图像作为金字塔n+1层的输入图像In+1，即In+1＝Lnd，并重复本段的内容，直至到金字塔顶层，即第J-1层，令金字塔顶层的LJ-1与次顶层的下采样结果相等LJ-2d，并将金字塔顶层LJ-1进行迭代自适应滤波获得拉普拉斯金字塔顶层的滤波结果LJ-1′并存入金字塔顶层；S5：自上而下地组合拉普拉斯金字塔中每层分解后滤波结果L0′-LJ-1′获得拉普拉斯金字塔，对存入金字塔内的滤波后图像L0′-LJ-1′统一进行重构获得最终的滤波结果。

全文数据：

权利要求：

百度查询： 中国矿业大学;中国自然资源航空物探遥感中心 一种基于拉普拉斯金字塔的InSAR干涉图迭代自适应滤波方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD