[发明专利]一种InSAR干涉图像的滤波方法及装置

权利要求书

1.一种InSAR干涉图像的滤波方法，其特征在于，包括：

将InSAR干涉图像分割成多个大小相同的图像块，通过计算每个图像与其它图像块的距离，获取每个图像块的若干个相似图像块，将每个图像块和其所有的相似图像块组成图像块组；

对每个图像块组的各图像块进行快速傅立叶变换处理得到原始的频谱数据，对所述图像块组的频谱数据的幅值进行平滑处理，得到滤波后的频谱数据，对原始的频谱数据和滤波后的频谱数据进行处理，得到所述图像块组的各图像块的估计值；

将每个图像块的估计值和所述图像块作为相似图像块的估计值进行加权平均，得到每个图像块的基本估计值；

将每个图像块的基本估计值作为经验值执行最终估计，得到滤波后的InSAR干涉图像。

2.根据权利要求1所述的InSAR干涉图像的滤波方法，其特征在于，通过计算每个图像与其它图像块的距离，获取每个图像块的若干个相似图像块，将每个图像块和其所有的相似图像块组成图像块组；包括：

计算第n个图像块Y_n与第j个图像块Y_j的距离d(Y_n,Y_j)：

其中，1≤n≤N，N为分割后图像块的个数，1≤j≤N,j≠n；H为图像块的宽，图像块的长和宽相等；||·||₂表示l²的范数；

判断d(Y_n,Y_j)是否小于预设的阈值，若为是，则图像块Y_j为图像块Y_n的相似图像块；

将第n个图像块及其所有的相似图像块组成图像块组S_n：其中，M_n为第n个图像块的相似图像块的个数，为第n个图像块的所有的相似图像块；由此得到N个图像块组：S₁,S₂…S_N。

3.根据权利要求2所述的InSAR干涉图像的滤波方法，其特征在于，对每个图像块组的每个图像块进行快速傅立叶变换处理得到原始的频谱数据，对所述图像块组的频谱数据的幅值进行平滑处理，得到滤波后的频谱数据，基于原始的频谱数据和滤波后的频谱数据进行处理，得到所述图像块组的每个图像块的估计值；包括：

对图像块组S_n中的每个图像块分别进行快速傅立叶变换，得到原始频谱数据F^w(u,v)_m，其中，u和v表示空间频率，m为图像块组S_n中图像块的序号；

将频谱幅值|F^w(u,v)_m|与选取的核函数K做平滑处理，得到

对进行归一化,得到归一化后的频谱幅值A_m：

对原始频谱数据F^w(u,v)进行平滑处理，得到滤波后的频谱数据

其中，·表示点乘，α_m为频域滤波器的滤波参数，计算公式为：

其中，γ为局部相干系数；

将原始频谱数据和滤波后频谱数据进行相乘，对乘积结果进行快速傅立叶逆变换，计算相位主值，得到图像块组S_n的每个图像块的估计值：其中，为第n个图像块Y_n的估计值，为相似图像块Y_n,1的估计值，为相似图像块的估计值。

4.根据权利要求3所述的InSAR干涉图像的滤波方法，其特征在于，将InSAR干涉图像的每个图像块的估计值及其作为相似图像块的估计值进行加权平均，得到每个图像块的基本估计值；包括：

对于第n个图像块，获取其作为相似图像块的所有估计值，表示为：其中，A_n为第n个图像块作为相似图像块的数量；

计算第n个图像块的基本估计值

其中，表示的特征函数；ω_k表示的权重值，其表达式为：

其中，σ²表示噪声方差。