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

技术领域

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

背景技术

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

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

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

发明内容

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

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

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

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

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

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

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

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

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

式中，H

所述步骤d中，模型训练和修正的具体流程如图1所示：首先将原始训练集随机分为5组，其中4组用于随机森林模型的训练，1组用于模型精度验证，其中随机森林模型为：

H

验证评价指标采用均方差回归损失(MSE)：

式中，N为用于验证的数据个数，H

然后确定影响随机森林模型精度的超参数(决策树个数和节点划分可选择特征变量数量最大值)的解空间范围，通过粒子群算法设定的粒子在解空间中速度和位置的迭代更新，比较不同超参数组合下适应度函数值，搜索适应度函数值最小的超参数组合，将此超参数组合作为随机森林模型的最优超参数组合，用整个训练集训练最优参数组合下的随机森林模型。

最后由训练得到的模型，根据数字高程模型每个像元对应的[Lat,Lon,Sl,As,Re,Gl]，预测像元对应的高程误差结果H

本发明的有益效果：本发明使用粒子群优化的随机森林方法进行修正，以进一步提高数字高程模型修正后的精度，选用SRTM作为试验所用数字高程模型，ICESat-2强光束地表光子数据作为试验所用参考高程控制点，Globeland30作为试验所用全球地表覆盖数据，NEON发布的机载LIDAR DTM作为试验所用验证数据，通过对本发明提出的方法(PSO-RF)和基于多项式回归的修正方法(PR)进行试验，以均方根误差(RMSE)作为验证指标，结果如图2所示：本发明提出的方法可以有效减小SRTM的高程误差，修正后的SRTM相对于修正前高程误差减小了42％-46％，且修正精度优于基于多项式回归的修正方法。

附图说明

图1是数字高程模型高程误差修正流程图。

图2是三条测线的原始SRTM、PR修正的SRTM和PSO-RF修正的SRTM高程精度评价结果。

具体实施方式

以下结合实施例和附图进一步解释本发明的具体实施方式，但不用于限定本发明。

本发明实施的实现过程是采用计算机实现基于粒子群优化随机森林的数字高程模型高程误差修正。以某区域内的SRTM数字高程模型、ICESat-2卫星测高数据和Globeland30地表覆盖数据为例，采用本发明方法修正SRTM的高程误差，包括以下步骤：

步骤a，读入ICESat-2卫星测高数据，选取分类为地表的光子(classed_pc_flag＝1)作为参考高程控制点，得到参考高程控制点处SRTM相对于参考高程控制点的高程误差

步骤b，将SRTM进行地理处理得到坡度、坡向、地形起伏度数据，利用双线性内插提取得到参考高程控制点处的地形参数[Sl,As,Re]，根据参考高程控制点所处的Globeland30像元，直接将此像元值作为参考高程控制点处的地表覆盖类型参数Gl。

步骤c，构建基于粒子群优化随机森林的数字高程模型高程误差修正模型。

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

步骤e，原始训练集随机分为5组，其中4组用于随机森林模型的训练，1组用于模型精度验证，其中随机森林模型为：H

步骤f，由训练得到的模型，根据数字高程模型每个像元对应的[Lat,Lon,Sl,As,Re,Gl]，预测像元对应的高程误差结果H

本发明是通过实施例进行描述的，本领域技术人员知悉，在不脱离本发明的精神和范围的情况下，可以对这些特征和实施例进行各种改变或等效替换。另外，在本发明的教导下，可以对这些特征和实施例进行修改以适应具体的情况及材料而不会脱离本发明的精神和范围。因此，本发明不受此处所公开的具体实施例的限制，所有落入本申请的权利要求范围内的实施例都属于本发明的保护范围。