[发明专利]一种基于多尺度时频分析的快速非局部均值InSAR相位滤波方法

权利要求书

1.一种基于多尺度时频分析的快速非局部均值InSAR相位滤波方法，其特征在于，结合小波包变换的多尺度时频分析特性以及非局部均值滤波的良好滤波效果，对InSAR相位信号进行去噪以平滑噪声，具体包括如下步骤：

S1：对InSAR相位信号进行正余弦变换处理；

S2：对S1步骤中经正余弦变换处理后的InSAR相位信号进行小波包变换，得含噪声的小波包系数；

S3：对S2步骤中所得含噪声的小波包系数进行快速非局部均值滤波，得滤波后的小波包系数；

S4：对S3步骤中所得滤波后的小波包系数进行小波包逆变换，得滤波后的InSAR相位信号。

2.根据权利要求1所述的一种基于多尺度时频分析的快速非局部均值InSAR相位滤波方法，其特征在于，S1步骤中所述正余弦变换处理在复平面内进行。

3.根据权利要求2所述的一种基于多尺度时频分析的快速非局部均值InSAR相位滤波方法，其特征在于，步骤S1中所述正余弦变换的公式为：

其中，Φ_z表示真实相位信号。

4.根据权利要求1所述的一种基于多尺度时频分析的快速非局部均值InSAR相位滤波方法，其特征在于，S2步骤中经小波包变换处理后的InSAR相位信号噪声为加性噪声模型，并满足如下公式：

其中，N_c表示cos(v)的均值，v_c与v_s表示均值为零的加性噪声，v表示相位信号噪声项。

5.根据权利要求4所述的一种基于多尺度时频分析的快速非局部均值InSAR相位滤波方法，其特征在于，S2步骤中所述含噪声的小波包系数包含InSAR相位噪声，并满足如下公式：

其中，WPT_2D{·}表示含噪声的小波包系数，2ⁱ表示小波尺度，为无噪声相位信号，表示零均值噪声。

6.根据权利要求1所述的一种基于多尺度时频分析的快速非局部均值InSAR相位滤波方法，其特征在于，步骤S3中所述滤波后的小波包系数包含的InSAR相位噪声为加性噪声，并满足如下公式：

x(i,j)＝s(i,j)+n(i,j)；

其中，x(i,j)表示滤波后的小波包系数，s(i,j)表示真实信号，n(i,j)则表示噪声信号。

7.根据权利要求1所述的一种基于多尺度时频分析的快速非局部均值InSAR相位滤波方法，其特征在于，步骤S3中非局部均值滤波算法采用积分图像的方法进行加速以降低算法复杂度。

8.根据权利要求1所述的一种基于多尺度时频分析的快速非局部均值InSAR相位滤波方法，其特征在于，滤波后的InSAR相位信号值为-π～π。