[发明专利]一种用于SAR干涉图处理的Goldstein金字塔构建方法

说明书

本发明公开了一种用于SAR干涉图处理的Goldstein金字塔构建方法，适用于图像处理领域。把SAR影像复数干涉图I0设为Goldstein金字塔的第0层图像G₀；设置对复数干涉图G₀进行Goldstein金字塔变换的最大层数L；设置频率域Goldstein低通滤波器的核大小N；利用大小为N的频率域Goldstein滤波器对金字塔第i层复数干涉图G_i做滑移滤波处理；对滤波后的干涉图进行下采样得到第i+1层Goldstein金字塔G_i+1；重复上述步骤直到得到第L层Goldstein金字塔G_L，把第0‑L层金字塔图像进行组合得到Goldstein金字塔。其可以生成SAR影像复数干涉图的Goldstein金字塔图像，实现对复数干涉图的多尺度描述，为进一步提升复数干涉图滤波、相位解缠等处理的效果提供了一种有效的多尺度变换方法，对SAR影像复数干涉图构建金字塔的效果好。

技术领域

本发明涉及一种Goldstein金字塔构建方法，尤其适用于图像处理技术领域中的用于SAR干涉图处理的Goldstein金字塔构建方法。

背景技术

金字塔构建属于图像处理领域，是一种图像多尺度表达技术，可以展示图像在不同分辨率下的信息，已被用于图像滤波、图像融合、图像超分辨率重建、运动识别、人脸识别等领域。一幅图像的金字塔是一组以金字塔形状排列的，由金字塔底层到顶层分辨率逐渐降低的图像，且来源于同一幅原始的图像。图像在金字塔中的层级越高，则图像越小，分辨率越低。

SAR干涉测量技术利用干涉相位信息精确测量地表目标的数字高程模型以及雷达视线向的微小形变。由于SAR传感器获取数据不受天气影响，因此能够实现全天时、全天候对地连续观测，这就使得SAR干涉测量技术在山体滑坡、沉降监测、地震形变等诸多领域有着广泛的应用。但受到时空去相干、热噪声去相干等影响，导致从实际影像得到的干涉图中存在大量相位噪声，需要通过滤波抑制干涉图中的噪声。此外，在从干涉图提取地表形变之前还需要进行相位解缠处理。在滤波、相位解缠等干涉图处理过程中构建金字塔是一种有效的多尺度、多分辨率分析方法。

高斯金字塔是最经典地图像金字塔构建方法，其首先使用二维高斯滤波器对原始图像进行滤波，然后采用一种图像下采样方法得到第一层金字塔图像。二维高斯滤波器是一种空间域滤波器，其对含有高斯噪声的普通的实数域图像具有很好的效果。但是，SAR干涉图是复数影像，由实数域向复数域扩展的高斯金字塔没有利用复数影像的频率特点

发明内容

针对上述技术的不足之处，提供一种理论严密、构建效果好的用于SAR干涉图处理的Goldstein金字塔构建方法。

为实现上述目的，本发明的用于SAR干涉图处理的Goldstein金字塔构建方法，其步骤如下：

a获取一幅列数为w、行数为h的复数干涉图I0＝A+B×j，其中A是复数干涉图I0的实部，B是复数干涉图I0的虚部，j表示虚数单位；

b建立Goldstein金字塔，根据复数干涉图I0的列数w和行数h，设置Goldstein金字塔变换的最大层数L；

c将原始SAR影像复数干涉图I0设为Goldstein金字塔的第0层图像G₀；

d设置频率域Goldstein低通滤波器的核大小N，利用Goldstein低通滤波器抑制复数干涉图的高频部分，保留复数干涉图的低频部分，把复数干涉图从空间域变换到频率域既可以使用普通傅里叶变换，也可以使用快速傅里叶变换；

e利用核大小为N的频率域Goldstein滤波器对金字塔第i层复数干涉图G_i做滑移滤波处理，得到滤波后干涉图

f对滤波后的干涉图进行下采样得到第i+1层Goldstein金字塔G_i+1；