[发明专利]一种地基SAR大气干扰相位校正方法及装置

权利要求书

1.一种地基SAR大气干扰相位校正方法，其特征在于，包括以下步骤：

S1：建立地基SAR差分干涉相位模型，所述地基SAR差分干涉相位模型表示为：

其中，为主图像j和从图像i的干涉相位，为主从图像采集过程中产生的形变相位，为大气干扰相位，为噪声相位；

S2：对所述地基SAR干涉相位模型进行滤波处理，以滤除其中的噪声相位；

S3：从滤除噪声相位的地基SAR差分干涉相位模型中筛选出强永久散射体，其中，所述强永久散射体为形变相位为零的稳定控制点；

S4：对步骤S3中筛选出的强永久散射体进行双向曲线拟合操作，以得到完整的大气干扰相位，其中所述双向曲线拟合包括距离向曲线拟合以及方位向曲线拟合；

S5：根据步骤S1中的地基SAR差分干涉相位模型和步骤S4中的所述大气干扰相位求得形变相位，以进行大气干扰相位校正。

2.根据权利要求1所述的方法，其特征在于，所述滤波方法为加权圆周期中值滤波，其中加权系数为其中，

3.根据权利要求1所述的方法，其特征在于，步骤S3中筛选出所述强永久散射体的步骤包括：

S31：利用第一特定范围内的相关系数初步筛选出永久散射体的候选点；

S32：通过第二特定范围内的振幅离差指数，从所述永久散射体的候选点中探测出永久散射体，以作为强永久散射体的候选点；

S33：利用时序SAR复图像的平均相关系数从所述强永久散射体的候选点中筛选出强永久散射体。

4.根据权利要求3所述的方法，其特征在于，所述第一特定范围为(0.8,1)，所述第二特定范围为(0,0.25)。

5.根据权利要求3所述的方法，其特征在于，步骤S33包括：

S331：求解出以所述强永久散射体的候选点为中心，半径为第一半径的一圆形区域内永久散射体的干涉相位的空间统计平均值其中，

其中，M为强永久散射体的候选点的个数；

S332：根据S331中求解出的所述空间统计平均值，求解各所述强永久散射体的候选点的时序相关系数

S333：通过使所述强永久散射体的候选点的相关系数等于选定的随机相位概率密度和强永久散射体的候选点的相位概率密度函数的加权之和，以确定筛选出所述强永久散射体时间相关系数的阈值；

S334：根据步骤S333中所确定的阈值从所述强永久散射体的候选点中筛选出强永久散射体。

6.根据权利要求1所述的方法，其特征在于，步骤S3还进一步包括对所筛选出的强永久散射体进行空间插值处理。

7.根据权利要求6所述的方法，其特征在于，步骤S3中的采用的插值方法为反距离加权法。

8.根据权利要求6所述的方法，其特征在于，在所述步骤S4中，对步骤S3中的插值后的强永久散射体进行双向曲线拟合操作，以得到完整的大气干扰相位。

9.根据权利要求8所述的方法，其特征在于，步骤S4中的双向曲线拟合操作包括以下步骤：

S41：对步骤S3中的插值结果进行距离向曲线拟合，以确定距离向的大气干扰相位；

S42：基于步骤S41的距离向曲线拟合，进行方位向曲线拟合，以得到完整的大气干扰相位。

10.一种应用如权利要求1-9中任意一项所述的地基SAR大气干扰相位校正方法的装置，其特征在于，所述装置包括：

模型建立单元，其用于建立地基SAR差分干涉相位模型，所述地基SAR差分干涉相位模型表示为：

其中，为主图像j和从图像i的干涉相位，为主从图像采集过程中产生的形变相位，为大气干扰相位，为噪声相位；以及滤波单元，其配置为对所述模型建立单元所建立的地基SAR干涉相位模型进行滤波处理，以滤除其中的噪声相位；

筛选单元，其配置为从所述滤波单元的滤波结果中筛选出强永久散射体，其中，所述强永久散射体为形变相位为零的稳定控制点；

拟合单元，其配置为对筛选单元所筛选出的强永久散射体进行双向曲线拟合操作，以得到完整的大气干扰相位；

计算单元，其根据所述拟合单元所得到的大气干扰相位计算出形变相位。