[发明专利]融合SAR与光学偏移量技术的地表三维形变计算方法及系统

权利要求书

1.融合SAR与光学偏移量技术的地表三维形变计算方法，其特征在于，包括以下步骤：

步骤1，获取待监测区域形变前后的两景SAR数据，将两景SAR数据处理为两景单视复数影像，利用SAR偏移量跟踪技术，根据两景单视复数影像计算待监测区域地表的方位向形变D_AZI和视线向形变D_LOS；

步骤2，获取待监测区域形变前后的两景光学数据，两景光学数据与两景SAR数据覆盖区域相同，利用偏移量技术，根据两景光学数据计算地表的南北向形变D_NS和东西向形变D_EW；

步骤3，对步骤1获得的方位向形变D_AZI、视线向形变D_LOS和步骤2获得的南北向形变D_NS和东西向形变D_EW进行重采样，使其具有相同的空间分辨率；

步骤4，通过公式(1)，获得地表的真实垂直向形变D_up、真实南北向形变D_ns和真实东西向形变D_ew；

式中，θ表示卫星入射角，表示卫星飞行方位角；

在使用公式(1)计算D_up、D_ns和D_ew时，采用赫尔默特方差分量估计的方法迭代来最终确定SAR与光学两类形变观测值的权，进而求得地表的真实垂直向形变量D_up、真实南北向形变量D_ns和真实东西向形变量D_ew。

2.如权利要求1所述的融合SAR与光学偏移量技术的地表三维形变计算方法，其特征在于，所述的步骤2包括：

步骤2.1，对于获取的两景光学数据，判断两景光学数据的覆盖范围大小是否相同，若相同，进行步骤2.2；若不同，裁剪出两景数据相同的公共区域，对两景公共区域进行步骤2.2；

步骤2.2，使用基于ENVI软件的COSI-Corr光学偏移量计算工具获得地表的南北向形变量D_NS和东西向形变量D_EW。

3.如权利要求1所述的融合SAR与光学偏移量技术的地表三维形变计算方法，其特征在于，所述的步骤2还包括对南北向形变D_NS和东西向形变D_EW的去除误差过程：根据得到的南北向形变D_NS和东西向形变D_EW掩模形变区域，在远离形变区域的位置选取若干像素点，根据这些像素点的形变值及地理位置利用一次拟合多项式拟合整个监测区域的轨道误差趋势面，然后将原始形变图与此误差趋势面相减，得到去除误差后的南北向形变D_NS和东西向形变D_EW。