[发明专利]一种基于局部信噪比的干涉图滤波方法

权利要求书

1.一种基于局部信噪比的干涉图滤波方法，所述方法通过窗口在所述干涉图中的滑动，依次对滑动窗口内的干涉图进行滤波，其中，在所述滑动窗口内对干涉图进行滤波包括：

对所述滑动窗口内的干涉图的干涉相位进行分解，得到所述干涉图在所述滑动窗口内的估计频率分量；

根据所述干涉相位和所述估计频率分量计算出所述干涉图在所述滑动窗口内的局部信噪比；

根据所述局部信噪比计算出滤波参数；以及

根据所述滤波参数对所述滑动窗口内的干涉图进行滤波。

2.根据权利要求1所述的方法，其中，所述对所述滑动窗口内的干涉图的干涉相位进行分解，得到所述干涉图在所述滑动窗口内的估计频率分量的步骤包括：

识别所述干涉相位的极大值点和极小值点；

根据所述极大值点和所述极小值点，计算极大值包络面和极小值包络面；

计算所述极大值包络面和极小值包络面的平均值；以及

根据所述滑动窗口内的干涉相位和所述平均值计算所述干涉图在所述滑动窗口内的估计频率分量。

3.根据权利要求2所述的方法，其中，所述根据所述滑动窗口内的干涉相位和所述平均值计算所述干涉图在所述滑动窗口内的估计频率分量的步骤包括：

从所述干涉相位中减去所述平均值，得到所述干涉图在所述滑动窗口内的估计频率分量；或者

判断所述频率分量是否满足预定条件；

如果不满足，则用所述频率分量代替所述干涉相位，从所述识别所述滑动窗口内的干涉相位的极大值点和极小值点的步骤开始迭代，得出满足所述预定条件的第一频率分量，将所述第一频率分量作为所述估计频率分量。

4.根据权利要求2所述的方法，其中，所述对所述滑动窗口内的干涉图的干涉相位进行分解，得到所述干涉图在所述滑动窗口内的估计频率分量的步骤还包括：

根据所述干涉相位和所述估计频率分量得到所述干涉相位的余量；

用所述余量代替所述干涉相位，从所述识别所述滑动窗口内的干涉相位的极大值点和极小值点的步骤开始预定次数的迭代，得出满足所述预定条件的一个或多个第二频率分量，将所述一个或多个第二频率分量作为所述估计频率分量。

5.根据权利要求1所述的方法，其中，所述局部信噪比是根据所述干涉相位和所述估计频率分量的关系计算得出的，所述局部信噪比是伪信噪比。

6.根据权利要求1所述的方法，其中，所述滤波参数是对所述局部信噪比进行归一化处理得到的。

7.根据权利要求1所述的方法，其中，所述根据所述滤波参数对所述滑动窗口内的干涉图进行滤波的步骤包括：

将所述滑动窗口内的干涉图从空间域变换到频域，得到所述干涉图的频谱；

用所述滤波参数对所述频谱进行频域平滑滤波；以及

将所述频域平滑滤波后的频谱变换到空间域。

8.根据权利要求2所述的方法，其中，所述极大值点和极小值点是通过像素点邻域搜索来识别的，所述极大值点和极小值点分别有一个或多个。

9.根据权利要求2所述的方法，其中，所述极大值包络面和所述极小值包络面是通过Delaunay三角剖分插值得到的。

10.根据权利要求1-9之一所述的方法，其中，所述干涉图是通过合成孔径雷达干涉测量InSAR得到的。