[发明专利]一种消除DSM中道路面高程异常的方法

权利要求书

1.一种消除DSM中道路面高程异常的方法，其特征在于，包括以下步骤：

步骤S1，通过提取道路高程算法从DSM数据中提取道路中心线逐点高程，先将道路面图层向内缓冲一定距离，缓冲距离为d1，紧接着利用掩膜提取的方法从DSM数据中提取道路面DSM数据，然后 将道路中心线按等距离法内插点，插点之间的距离为d2，从而获取插点之后道路中心线所有节点处垂线上DSM的最小高程值，由此得到赋高程后道路中心线；

步骤S2，通过行道树清理算法处理道路中心线高程异常，从而将行道树的高程还原为道路实际高程，调整坡度阈值和平地距离阈值，进而达到行道树清理的目的，并且行道树清理算法同样适用于道路两旁因具有建筑而导致的道路面高程异常；

步骤S3，通过算术平均法处理道路中心线高程异常，计算连续5个节点高程的平均值，并将计算结果赋值给中间点，从而达到消除高程异常的目的；

步骤S4，通过其他方法处理道路中心线高程异常；

步骤S5，道路中心线高程赋值给道路面，遍历道路面上的所有节点，提取离本节点最近的道路中心线的节点高程，并赋值给本道路上的节点；

步骤S6，将道路面数据与DSM数据进行叠加，依据道路面高程修改DSM高程，并更新DSM数据的高程。

2.根据权利要求1所述的一种消除DSM中道路面高程异常的方法，其特征在于，步骤S1中所述提取道路高程算法的操作步骤如下：

（1）道路面图层向内缓冲后形成道路面缓冲图层，并且以道路面缓冲图层作为输出图层，并且道路面图层向内缓冲的距离d1可以由用户自己设置，也可以选择默认值，默认值为DSM分辨率的一半；

（2）按照道路面缓冲图层的范围提取DSM数据，输出道路面DSM图层；

（3）通过道路面图层提取道路中心线；

（4）沿道路中心线间隔指定距离插入点，并为所有的道路中心线以及插入点进行编号；

（5）在每个插入点做本中心线的垂线，垂线的长度必须大于道路宽度，沿垂线按指定距离生成等分点，并记录等分点所在的道路中心线ID和插入点ID，将生成的等分点保存到高程提取点图层中；

（6）高程提取点图层与道路面DSM图层叠加后，将每个高程提取点与对应的道路面DSM图层上的高程值提取出来，并保存到高程提取点的属性中；

（7）将高程提取点图层中高程值为空的点删除，统计道路中心线ID和插入点ID一样的所有高程提取点高程的最小值，并将该值赋给道路中心线对应的插入点，最后得到赋高程后道路中心线。

3.根据权利要求1所述的一种消除DSM中道路面高程异常的方法，其特征在于，步骤S2中，所述行道树清理算法的流程如下：

第一步：为道路中心线从起点开始按节点顺序编号，编号顺序从起点开始逐点递增，从倒数第二个点开始逆序遍历所有道路节点，用户自定义坡度阈值与平地距离阈值；

第二步：令当前节点作为第i个节点，相对应的，位于当前节点的下一个节点为第i+1个节点；

第三步：计算第i个节点与第i+1个节点之间的坡度，判断坡度是否小于坡度阈值，若坡度小于坡度阈值，则进入第四步，若坡度不小于坡度阈值，则进入第五步；

第四步：判断是否存在平地起点记录，若不存在平地起点记录，则进入第六步，若存在平地起点记录，则进入第七步；

第五步：判断是否存在斜坡起点记录，若不存在斜坡起点记录，则进入第八步，若存在斜坡起点记录，则进入第九步；

第六步：将第i+1个节点作为平地起点，然后进入第七步；

第七步：计算当前节点到平地起点的距离d，并判断距离d是否大于平地距离阈值，若距离d大于平地距离阈值，则进入第十步，若距离d不大于平地距离阈值，则进入第十三步；

第八步：将第i+1个节点作为斜坡的起点，然后进入第九步；

第九步：清除平地起点记录，然后进入第十三步；

第十步：判断是否存在斜坡起点记录，若存在斜坡起点记录，则进入第十一步，若不存在斜坡起点记录，则进入第十三步；

第十一步：将平地的起点作为斜坡的终点，按距离内插斜坡起点和终点间的高程，或者删除斜坡起点和终点之间的节点；

第十二步：清理斜坡起点记录，清理平地起点记录；

第十三步：返回第二步，逆序继续处理下一个道路节点。