一种多源数据高精度时序地表三维形变解算方法和系统与流程

本发明涉及地表形变三维动态监测与跟踪，特别涉及一种多源数据高精度时序地表三维形变解算方法和系统。

背景技术：

1、随着空间大地测量技术的发展与进步，全球导航卫星系统gnss（globalnavigation satellite system，gnss）与合成孔径雷达干涉测量（interferometricsynthetic aperture radar，insar）作为两种新型对地观测技术在地壳运动、地质灾害预警与防治等方面的应用越来越突出。

2、现有技术中，gnss、insar多源数据融合地表三维形变研究主要集中在两点，第一是灾后三维形变场重建，第二是解算三维形变速率场。前者是直接融合多种技术获得的insar和gnss形变观测值，获取瞬时三维形变场；后者是对各类形变观测值求平均速率后融合，获取三维形变速率场。无论上述哪种情况，当仅有单一insar视线向数据时，无法获取三维形变信息，大大牺牲了insar数据的时间分辨率；而且无法获取地表形变前后时刻的关系，不能反映地表形变的动态过程。

技术实现思路

1、为克服现有技术中存在的问题，本发明的目的是提供一种多源数据高精度时序地表三维形变解算方法和系统，具体的，是一种基于kalman滤波gnss-insar高精度时序地表三维形变解算方法。

2、本发明针对一般多源数据融合模型解算地表三维形变只能获取某一时刻的瞬时形变或者某段时间内的平均三维形变，而无法反映地表形变的动态过程，牺牲了观测值的时间分辨率的问题，提出基于kalman滤波gnss-insar高精度时序地表三维形变解算方法，该方法通过构建区域三维形变监测的观测方程和运动方程，时序地融合gps、in sar 这两种不同传感器在不同时刻、不同对地观测角度下获取的观测数据，得到每一个形变信息获取时刻的三维地表形变量。

3、本发明由下述技术方案实现：

4、本发明的第一方面涉及一种多源数据高精度时序地表三维形变解算方法，包括如下步骤：

5、步骤s100，对gnss数据进行预处理；

6、步骤s200，对insar数据进行预处理；

7、步骤s300，建立区域三维形变监测的运动方程和观测方程；

8、步骤s400，滤波融合解算。

9、进一步的，步骤s100中，对gnss监测值进行处理，得到gnss点变形监测结果，对变形监测结果进行插值，得到符合空间分辨率的gnss数据。

10、进一步的，步骤s200中，对insar数据去除空间相关误差处理，得到视线向形变时间序列，并对变形监测结果进行降采样处理，使其和gnss插值结果空间分辨率一致。

11、进一步的，步骤s300中，基于监测点待求形变位移值、瞬时速率、瞬时加速率建立离散型卡尔曼滤波的运动方程；根据insar视线向形变与其在地表三维方向上的投影的关系建立观测方程。

12、进一步的，步骤s400包括：

13、步骤s410，确定滤波初始值；

14、步骤s420，根据状态方程以及观测方程，在最小均方差估计的准则下进行滤波递推估计。

15、进一步的，步骤s420包括：

16、步骤s421，一步预测当前时刻形变状态值：利用上一时刻的状态量，借助状态转移矩阵，得到状态一步预测；

17、步骤s422，一步预测当前时刻形变状态值方差：利用上一时刻状态向量协方差和动力学模型噪声向量，求一步预测误差方差矩阵；

18、步骤s423，求滤波增益矩阵：利用观测向量和观测矩阵求滤波增益矩阵；

19、步骤s424，得到滤波更新状态向量和相应协方差阵。

20、利用不同时刻获取的insar或gnss观测值进行滤波解算，得到每一数据获取时刻的三维形变状态向量估值，对监测区域的三维形变动态监测。

21、本发明还涉及一种多源数据高精度时序地表三维形变解算系统，包括：

22、gnss数据预处理模块，用于对gnss数据进行预处理；

23、insar数据预处理模块，用于对insar数据进行预处理；

24、方程建立模块，用于建立区域三维形变监测的运动方程和观测方程；

25、滤波融合解算模块，用于滤波融合解算。

26、本发明还涉及一种电子设备，所述电子设备包括：

27、至少一个处理器；以及，

28、与所述至少一个处理器通信连接的存储器；其中，

29、所述存储器存储有可被所述至少一个处理器执行的指令，所述指令被所述至少一个处理器执行，以使所述至少一个处理器能够执行所述的方法。

30、本发明还涉及一种非暂态计算机可读存储介质，该非暂态计算机可读存储介质存储计算机指令，该计算机指令用于使该计算机执行所述的方法。

31、本发明的技术方案能实现如下有益的技术效果：

32、本发明的基于kalman滤波gnss-insar高精度时序地表三维形变解算算法，可以解算任一数据获取时刻的三维形变信息，能够实现对监测区域地表三维形变的实时动态估计，获取高精度时序三维形变场和三维形变速率场，提高地表三维形变监测的时间和空间分辨率。

技术特征：

1.一种多源数据高精度时序地表三维形变解算方法，其特征在于，包括如下步骤：

2.根据权利要求1所述的多源数据高精度时序地表三维形变解算方法，其特征在于，步骤s400包括：

3.根据权利要求1所述的多源数据高精度时序地表三维形变解算方法，其特征在于，步骤s420包括：

4.一种多源数据高精度时序地表三维形变解算系统，其特征在于，包括：

5.一种电子设备，其特征在于，所述电子设备包括：

6.一种非暂态计算机可读存储介质，该非暂态计算机可读存储介质存储计算机指令，该计算机指令用于使该计算机执行前述权利要求1至3中任一项所述的方法。

技术总结
一种多源数据高精度时序地表三维形变解算方法和系统，包括如下步骤：步骤S100，对GNSS数据进行预处理；步骤S200，对InSAR数据进行预处理；步骤S300，建立区域三维形变监测的运动方程和观测方程；步骤S400，滤波融合解算。本发明在基于多源形变监测数据解算地表三维形变时，可以利用上一时刻的状态量和当前时刻的观测量获得当前数据获取时刻的地表三维形变场，实现了动态三维形变监测；当前时刻仅有单一InSAR视线向观测数据时，亦可得到该时刻的地表三维形变信息。

技术研发人员：李子申,杨雪艳,孙志成,王亮亮,王宁波,潘军道,刘振耀,任仰森
受保护的技术使用者：齐鲁空天信息研究院
技术研发日：
技术公布日：2024/1/14