一种复杂地形下的地表变形监测方法

1.一种复杂地形下的地表变形监测方法，其特征在于，包括以下过程：基于局部高精度无人机模型，采用虚拟控制点重构视野盲区的无人机模型并进行盲区内的地表变形分析，然后利用提出的乘幂函数进行误差估计，所述的乘幂函数y＝a*x b ，其中a和b为未知系数，x为自变量，y为因变量；最后通过计算位移来估计实际位移；包括以下过程：S1、利用虚拟控制点重构无人机模型从某一期无人机的视野盲区部分提取一个虚拟控制点ACP，并将其赋予另一期无人机数据并进行模型重构，重构后的模型虽然发生扭曲，但能精确配准并进行地表变形分析；此时进行地表位移计算，其计算值CD和真实位置RD之间存在误差，该误差在模型沿极轴拉伸时为：RD j min ＝CD j -[ε(a j ′+CD j )-ε(a j ′)]      (1)同理可得，在沿极轴压缩时为：RD j max ＝CD j +[ε(a j ′+CD j )-ε(a j ′)]       (2)根据几何关系可知，旋转时的RD j 应满足以下关系：RD j min ≤RD j ≤RD j max (3)其中RD j ，CD j 分别表示点x j 的实际位移和计算位移，ε(d)表示在控制点范围外的无人机模型距离控制点为d时的无人机模型误差，a j ′表示在扭曲的无人机模型中测量到的x j 在滑坡前到控制点的距离，a j ′+CDj则表示在扭曲的无人机模型中测量到的x j 在滑坡后到控制点的距离；S2、ε(d)的函数分布ε(d)表示在控制点范围外的无人机模型距离控制点为d时的无人机模型误差，进行了12种控制点布置情况下的误差分布模拟，每种情况下均在一个局部区域布置了3个控制点，然后分析这些情况下无人机模型精度的变化；根据误差的分布进行统计，获得了ε(d)的分布函数：ε(d)＝(2.64×10 -4 ±6.65×10 -5 )*d 1.76±0.04 (4)S3、CD和RD之间的误差估计利用公式(1)-(4)，对CD j 与RD j 之间的误差进行估计。