[发明专利]多轨道InSAR震间形变速率场拼接方法、装置、设备及介质

权利要求书

1.一种多轨道InSAR震间形变速率场拼接方法，其特征在于，包括：

获取目标区域的全球导航卫星系统GNSS的三维形变速率和多轨道合成孔径雷达SAR数据；

将所述三维形变速率投影至目标向，得到第一目标向形变速率；

将所述SAR数据进行时序处理，得到每个轨道的第二目标向形变速率；

基于所述第一目标向形变速率对所述第二目标向形变速率进行参考基准统一，得到每个轨道的第三目标向形变速率；

对所述第三目标向形变速率进行入射角校正，得到每个轨道的最终目标向形变速率，在地理坐标系下将所述每个轨道的最终目标向形变速率进行拼接；

若基于时序InSAR技术获取的SAR卫星的形变速率包括水平运动和垂直运动，将所述GNSS的三维形变速率投影至目标向，得到第一目标向形变速率，包括：获取所述三维形变速率对应的垂直向形变速率、南北向形变速率、东西向形变速率和SAR卫星的入射角和方位角；基于第一计算公式对所述垂直向形变速率、所述南北向形变速率、所述东西向形变速率和所述SAR卫星的入射角和方位角进行计算，得到LOS视线向三维形变速率作为所述第一目标向形变速率；

所述对所述第三目标向形变速率进行入射角校正，得到每个轨道的最终目标向形变速率，包括：基于SAR卫星的方位角，将各轨道的所述第三目标向形变速率转换到SAR坐标系下，并获取相邻轨道所述第三目标向形变速率之间的第二目标差值；使用二阶二次多项式以所述第二目标差值作为因变量、SAR卫星视线向作为自变量进行拟合，并采用最小二乘求解第二待估参数；基于所述第二待估参数确定第二形变速率比值；基于所述第二形变速率比值和所述第三目标向形变速率进行计算，得到所述每个轨道的最终目标向形变速率；

若所述SAR卫星仅包括水平运动，则基于第二计算公式对所述南北向形变速率、所述东西向形变速率和所述SAR卫星的入射角和方位角进行计算，得到第一地距向形变速率作为所述第一目标向形变速率；

所述对所述第三目标向形变速率进行入射角校正，得到每个轨道的最终目标向形变速率，包括：将所述第二目标向形变速率投影到地距向，得到第二地距向形变速率；计算所述第一地距向形变速率和所述第二地距向形变速率之间的第三目标差值；使用二阶二次多项式以所述第三目标差值作为因变量、地理经纬度坐标作为自变量进行拟合，并采用最小二乘求解第三待估参数；基于所述第三待估参数确定第三形变速率比值；基于所述第三形变速率比值和所述第二地距向形变速率进行计算，得到所述每个轨道的最终目标向形变速率。

2.根据权利要求1所述的多轨道InSAR震间形变速率场拼接方法，其特征在于，所述基于所述第一目标向形变速率对所述第二目标向形变速率进行参考基准统一，得到每个轨道的第三目标向形变速率，包括：

对所述第一目标向形变速率和所述LOS向形变速率进行差值处理，得到第一目标差值；

使用二阶二次多项式以所述第一目标差值作为因变量、地理经纬度坐标作为自变量进行拟合，并基于最小二乘求解第一待估参数；

基于所述第一待估参数确定第一形变速率比值；

基于所述第一形变速率比值和所述LOS向形变速率进行计算，得到所述每个轨道的第三目标向形变速率。

3.根据权利要求1所述的多轨道InSAR震间形变速率场拼接方法，其特征在于，所述在地理坐标系下将所述每个轨道的最终目标向形变速率进行拼接，包括：

基于任意相邻两个轨道在所述地理坐标系下的最终目标向形变速率进行叠加，得到拼接结果。