[发明专利]基于FFT高精度估计干涉相位梯度的InSAR滤波方法及系统

权利要求书

1.一种基于FFT高精度估计干涉相位梯度的InSAR滤波方法，其特征在于，所述方法包括如下步骤：

步骤(1)：对配准后的主辅SAR复图像进行干涉处理，生成初始复干涉相位图s₁(x,y)，然后对s₁(x,y)进行3×3均值滤波得到粗滤波复干涉相位图s₂(x,y)；其中，x和y分别为初始复干涉相位图中任意像素的二维坐标；

步骤(2)：在粗滤波复干涉相位图s₂(x,y)中设置干涉相位梯度估计窗口，基于FFT估计得到初始干涉相位梯度估计图和估计质量图其中，p和q为初始干涉相位梯度估计图中任意像素的二维坐标，为横轴方向初始干涉相位梯度估计图，为纵轴方向初始干涉相位梯度估计图；

步骤(3)：根据估计质量图剔除初始干涉相位梯度估计图中的质量低于0.4的干涉相位梯度估计值，得到最终干涉相位梯度估计图和

步骤(4)：利用最终干涉相位梯度估计图和补偿滤波窗口中的干涉相位梯度，对初始复干涉相位图s₁(x,y)进行均值滤波，得到滤波后复干涉相位图s₃(x,y)。

2.根据权利要求1所述的基于FFT高精度估计干涉相位梯度的InSAR滤波方法，其特征在于：在所述步骤(1)中，对初始复干涉相位图s₁(x,y)进行3×3均值滤波的方法为：在初始复干涉相位图中设置滤波窗口，窗口尺寸为3×3，然后对窗口内像素进行复数平均。

3.根据权利要求2所述的基于FFT高精度估计干涉相位梯度的InSAR滤波方法，其特征在于：粗滤波复干涉相位图s₂(x,y)为N为滤波窗口中像素总数，-1≤l_x,l_y≤1为两个整数变量，随着l_x和l_y取值不同，(x+l_x,y+l_y)表示以(x,y)为中心滤波窗口中的不同位置，x和y分别为初始复干涉相位图中任意像素的二维坐标。