[发明专利]一种基于粒子群优化随机森林的数字高程模型高程误差修正方法

说明书

本发明公开了一种基于粒子群优化随机森林的数字高程模型高程误差修正方法，属于遥感技术领域。本发明使用粒子群优化的随机森林方法进行修正，以进一步提高数字高程模型修正后的精度，选用SRTM作为试验所用数字高程模型，ICESat‑2强光束地表光子数据作为试验所用参考高程控制点，Globeland30作为试验所用全球地表覆盖数据，NEON发布的机载LIDAR DTM作为试验所用验证数据，通过对本发明提出的方法和基于多项式回归的修正方法进行试验，以均方根误差作为验证指标，本发明可以有效减小SRTM的高程误差，修正后的SRTM相对于修正前高程误差减小了42％‑46％，且修正精度优于基于多项式回归的修正方法。

技术领域

本发明涉及一种基于粒子群优化随机森林的数字高程模型高程误差修正方法，属于遥感技术领域。

背景技术

数字高程模型作为地表高程的数字化表达，被广泛应用于地貌、水文、测绘、灾害监测和控制等领域。但是数字高程模型受观测手段、地形条件、植被覆盖等因素的影响，其高程存在不可忽略的误差，且在不同区域高程精度往往差异较大。

对于数字高程模型所存在的高程误差，一些学者尝试利用各类精度更高的高程数据对数字高程模型高程误差进行修正，如高精度GPS测量点、机载激光雷达高程数据、高精度DEM数据。但是，以上高程数据受到分布区域有限、获取和制作难度大等因素的影响，难以对大范围、任一区域的SRTM进行高程误差修正。ICESat卫星测高数据因其覆盖全球、测高精度高等优势，被逐步应用于数字高程模型修正中，但ICESat在2009年失效，其数据也相应停止更新，无法实现更具现势性的SRTM修正。目前，ICESat-2测高卫星仍在运行，可以提供近全球范围内高精度的测高数据，Magruder等基于ICESat-2数据，结合Landsat8影像，采用多项式回归方法建立数字高程模型高程误差修正模型。但是，由于数字高程模型高程误差与其影响因素之间往往是复杂的非线性关系，简单数学表达的多项式回归方程难以充分表达这种关系，这就造成了应用此种方法修正的数字高程模型高程精度仍存在较大的局限性。

随机森林作为一种可用于解决非线性回归问题的机器学习算法，具有精度高、抗噪声能力强、不易发生过拟合的优势，但该方法精度会受到所设定超参数的影响，因此本发明通过融合粒子群算法寻找随机森林的最优超参数组合值，以粒子群优化随机森林的方法对数字高程模型进行高程误差的修正。

发明内容

为了实现更高精度数字高程模型修正结果，本发明提供了一种基于粒子群优化随机森林的数字高程模型高程误差修正方法。

本发明为实现上述目的所采用的技术方案是：

一种基于粒子群优化随机森林的数字高程模型高程误差修正方法，包括以下步骤：

a、从测高卫星数据中提取参考高程控制点，计算相对于参考高程控制点的数字高程模型高程误差及提取参考高程控制点对应的经纬度、地形参数和地表覆盖类型参数；

b、构建基于粒子群优化随机森林的数字高程模型高程误差修正模型；

c、用参考高程控制点处的经纬度、地形参数和地表覆盖参数作为修正模型输入数据，高程误差作为修正模型目标数据，建立训练模型的训练集；

d、用训练集对修正模型进行训练，将数字高程模型应用于训练所得修正模型中，进行高程误差的修正。

所述步骤a中，提取的地形参数为数字高程模型的坡度Sl、坡向As和地形起伏度Re，提取的地表覆盖类型参数来自全球地表覆盖数据Gl；对于参考高程控制点对应地形参数的计算，采用双线性插值法在参考高程控制点处进行提取；对于地表覆盖类型参数则在参考高程控制点处直接进行提取。

所述步骤b中，高程误差修正模型为：