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

权利要求书

1.一种用于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；

g重复步骤(e)和(f)直到得到第L层Goldstein金字塔G_L，把第0-L层金字塔图像由下至上组合得到层数为L+1的Goldstein金字塔。

2.根据权利要求1所述的用于SAR干涉图处理的Goldstein金字塔构建方法，其特征在于：所述步骤c中最大层数L满足条件为：和式中log为取对数符号。

3.根据权利要求1所述的用于SAR干涉图处理的Goldstein金字塔构建方法，其特征在于：频率域Goldstein滤波器的核大小N，根据选择的傅里叶变换方法不同，要满足条件d1或d2任一：

d1当Goldstein滤波算法使用普通傅里叶变换实现时，N为奇数，3≤N≤min(w,h)且N∈N^*，其中N^*表示正整数，min(·)表示取最小值，特别地，N建议设置为N＝5；

d2当Goldstein滤波算法使用快速傅里叶变换实现时，N为偶数，2≤N≤min(w,h)，特别地，N建议设置为N＝4或N＝6。

4.根据权利要求1所述的用于SAR干涉图处理的Goldstein金字塔构建方法，其特征在于：步骤e中用到的金字塔层数i满足条件：0≤i≤L-1，对干涉图G_i滑移滤波处理使用的滑移步长为S＝1。

5.根据权利要求1所述的用于SAR干涉图处理的Goldstein金字塔构建方法，其特征在于：步骤f中对滤波后干涉图进行下采样的方式为直接移除图像偶数行和偶数列，利用的奇数行和奇数列组成G_i+1。

6.根据权利要求1所述的用于SAR干涉图处理的Goldstein金字塔构建方法，其特征在于：步骤g中Goldstein金字塔的组合方式为按照层数由小到大，由底层到顶层排列，第0层金字塔图像G₀在最底层，第L层金字塔图像G_L在最顶层。