[发明专利]一种基于Insarbm3d滤波算法的迭代反演滤波方法

权利要求书

1.一种基于Insarbm3d滤波算法的迭代反演滤波方法，其特征在于，包括如下步骤：

1）干涉相位图反演：

将干涉相位图φ采用基于N.Goodman信号模型的反演算法进行反演， N.Goodman信号模型中，当相干系数为1时，干涉相位图φ反演为复数配准图对u₁、u₂，u₁为复数配准图对主图，u₂为复数配准图对副图，u₁=A₁×v₁，u₂=A₂×v₁×e^-jφ，则复数干涉图I为：

I=u₁×（u₂^*）=A₁×A₂×e^jφ+(|v₁|²-1)×A₁×A₂×e^jφ （1）

其中A₁是配准图u₁的振幅，A₂是配准图u₂的振幅，|v₁|²是v₁的模的平方，A₁×A₂×e^jφ是干涉相位图φ经反演算法反演出的不包含噪声的复数干涉图，(|v₁|²-1)A₁×A₂×e^jφ为反演过程引入的误差；

2）Insarbm3d滤波：

将步骤1）得到的反演出的复数配准图对u₁、u₂作为Insarbm3d算法的输入数据，使用Insarbm3d滤波算法进行滤波，对输入的复数配准图对u₁、u₂共轭相乘取相位之后对应的干涉相位图φ₁进行降噪，得到滤波后的干涉相位图φ₂；其中干涉相位图φ₁与u₁、u₂的关系表达式为φ₁=angle（u₁×（u₂^*））；干涉相位图φ1为干涉相位图φ在反演算法引入噪声后的干涉相位图；具体是：

3）计算干涉相位图φ经反演滤波前后的性能指标：

3-1）计算干涉相位图φ进行反演滤波前的相位差和值SPD、相位标准偏差PSD和残数residual；

3-2）计算干涉相位图φ₂的相位差和值SPD、相位标准偏差PSD和残数residual；

4）根据干涉相位图φ与干涉相位图φ₂性能指标的变化，设置3个迭代终止条件，当性能指标满足其中任意一个条件，则终止迭代反演滤波；否则重复步骤1）至步骤3）；将干涉相位图φ的相位差和值SPD、相位标准偏差PSD和残数residual分别记为SPDφ、PSDφ、residualφ，干涉相位图φ₂的相位差和值SPD、相位标准偏差PSD和残数residual分别记为SPDφ₂、PSDφ₂、residualφ₂，则3个迭代终止条件表达式为：

条件1：(（SPDφ-SPDφ₂)÷SPDφ)0.02

条件2：(（PSDφ-PSDφ₂)÷PSDφ)0.02

条件3：（residualφ-residualφ₂）0 。

2.根据权利要1所述的一种基于Insarbm3d滤波算法的迭代反演滤波方法，其特征在于，步骤2）中，所述的Insarbm3d滤波包括两个过程，具体是：先对对噪声图像进行滤波，产生一个估计后；相同的噪声图像再次滤波；其中第二次滤波过程是由第一次滤波的估计结果驱动；两个滤波过程都包括：分组、协同滤波和聚合三个步骤；在分组步骤中，对于每个图像块，在整个图像上挑出许多相似的块，并在一个三维堆栈中收集；块的选择考虑噪声的统计量，基于块相似度度量；在协同滤波中，同时考虑局部即块内和非局部即块间的依赖关系；这些依赖关系在协同滤波步骤中被利用，以产生整个堆栈的去噪版本，堆栈中的每个块都会协同恢复所有其他块；在聚合步骤中，将堆栈中的所有过滤块返回到它们在图像中的原始位置；由于在每个堆栈中过滤了几个块，因此最终会生成大量重叠的过滤块；因此，每个图像像素都有多个估计，通过对多个估计值进行加权平均。