一种SAR图像配准联系点粗差剔除方法技术

本发明专利技术涉及合成孔径雷达干涉测量数据处理技术领域，具体为一种SAR图像配准联系点粗差剔除方法，包括如下步骤：基于互相干系数，计算出两幅SAR影像配准联系点的偏移量；根据最小截尾二乘法（LTS）计算未知参数与单位权中误差的稳健初值；基于LTS稳健初值，选取合适的权函数计算等价权初值，在此基础上，迭代剔除联系点偏移量粗差，获取未知参数在抗粗差前提下的最优估值；利用未知参数在抗粗差前提下的最优估值，求解相对几何变换模型；基于相对几何变换模型，完成SAR图像配准。与现有技术相比，本发明专利技术具有粗差剔除稳健性好、运算效率高等优势。

【技术实现步骤摘要】
一种SAR图像配准联系点粗差剔除方法
本专利技术涉及测绘科学或合成孔径雷达干涉测量数据处理
，具体为一种SAR图像配准联系点粗差剔除方法。
技术介绍
合成孔径雷达干涉测量(InterferometricSyntheticApertureRadar,InSAR)技术是近年来发展迅速的空间对地观测新技术，因其覆盖范围广、测量精度高等优势，已被广泛应用于DEM生成和大范围地表形变监测。SAR图像配准是差分干涉测量处理的基础和重要环节之一。由于多次成像时卫星的轨道、姿态、SAR传感器工作状态并不完全相同，每幅影像的覆盖范围、入射角度、斜距及方位向分辨率等影像几何关系都在一定范围内有少量变化。利用影像间联系点的偏移量拟合最佳相对几何变换模型，将待配准影像重采样至主影像几何空间，使地面上回波点在两幅影像中处于同一位置，继而进行差分干涉处理生成有效的干涉相位图。影像联系点偏移量的计算通常采用基于滑动窗口的互相干系数法，互相干系数最大所对应窗口中心像元的偏移量即为影像联系点偏移量。由于窗口内像元后向反射性质决定互相干系数的计算精度，利用互相干系数计算的影像联系点偏移量不可避免的包含粗差，如果不本文档来自技高网...

【技术保护点】
一种SAR图像配准联系点粗差剔除方法，其特征在于，包括以下步骤：步骤一，利用互相干系数计算两幅SAR图像配准联系点偏移量；步骤二，利用最小截尾二乘法(LTS)计算未知参数与单位权中误差的稳健初值

【技术特征摘要】
1.一种SAR图像配准联系点粗差剔除方法，其特征在于，包括以下步骤：步骤一，利用互相干系数计算两幅SAR图像配准联系点偏移量；步骤二，利用最小截尾二乘法(LTS)计算未知参数与单位权中误差的稳健初值与步骤三，粗差剔除，基于LTS稳健初值，选取合适的权函数计算等价权初值在此基础上进行选权迭代，计算第k轮迭代未知参数估值如果满足：则认为偏移量粗差被全部剔除，即为未知参数在抗粗差前提下的最优估值；步骤四，利用未知参数在抗粗差前提下的最优估值，求解相对几何变换模型；步骤五，基于相对几何变换模型，完成SAR图像配准。2.根据权利要求1所述的一种SAR图像配准联系点粗差剔除方法，其特征在于，步骤一中所述基于互相干系数计算联系点偏移量的方法为：假设主影像为M，待配准影像为S，利用滑动窗口法计算M与S的互相干系数Γ，Γ最大所对应窗口中心像元的偏移量即为联系点偏移量；互相干系数的计算公式为：式中，m与s分别代表M与S窗口内的像元，E(·)代表数学期望，*为共轭计算符。3.根据权利要求1所述的一种SAR图像配准联系点粗差剔除方法，其特征在于，步骤二中所述未知参数与单位权中误差稳健估值的计算方法为：对配准联系点偏移量按Γ的大小重新排序，取前h个偏移量作为观测向量，利用最小二乘法(LS)计算未知参数与偏移量残差的稳健估值h的计算公式为：h＝int[n/2]+int[(b+1)/2]式中，b为未知参数个数，n为偏移量总数；单位权中误差稳健估值的计算...

【专利技术属性】
技术研发人员：朱邦彦，储征伟，陈昕，姚冯宇，
申请(专利权)人：南京市测绘勘察研究院股份有限公司，
类型：发明
国别省市：江苏,32