技术摘要：
本发明属于信息采集和处理技术领域，以解决山地作业难度评估方式不能全面地考虑作业点距离道路的远近、坡度大小、作业点分布的问题。山地作业难度指数计算方法：获取目标地区的数字图像；统一上述图像的分辨率并合成一张多信息融合图，道路网数据转化为栅格数据，道路 全部
背景技术：
山地作业的工作难度评估因其复杂、繁琐、计算量大等原因，长久以来一直没有得 到有效的解决。实际工作中一直沿用传统的一口价方式，这种估价方式较为简单，当工作量 较小时尚可，但是当工作量很大时，它的缺点就特别突出。主要有以下几个方面：(1)没有考 虑到作业点距离道路的远近。作业地点距离道路越远，作业人员需要花费的时间和精力就 越多，特别是当携带较重的工具的时候。(2)没有考虑作业点到道路坡度大小。在相同的距 离长度上，如山地的坡度越大，那么山路就越陡峭，同样其所消耗的时间和精力也越多。(3) 没有考虑作业点在山地上的分布情况。通常来说，采用统一定价的方式会综合考虑作业地 区山地的地形情况，按照平均难度给出个估计价格。但由于作业点的分布并不是完全均衡 的，有的时候会出现作业点大多靠近道路，有的时候会出现作业点大多远离道路，这就不满 足按照平均难度给出估价的初衷。(4)主观性强，没有确定的估价标准，结果往往因人而异。 传统的估价方式的弊端日趋明显，但是，山地作业的工作难度评估因其复杂、繁琐、计算量 大等原因，有效的切实可行的评估方法一直没有出现，上述问题也一直没有得到解决。鉴于 此，提出本发明。
技术实现要素：
现有的山地作业难度评估方式较为简单，不能够全面地考虑作业点距离道路的远 近、作业点到道路坡度大小、作业点分布情况，反映在工作定价上大多采用一口价的方式， 弊端明显。针对上述问题，本发明提出一种山地作业难度指数计算方法，综合利用了地理信 息系统和计算机图形学等技术，逐块区域的评估该区域工作的难度，逐点的计算出山地上 每点的工作难度值，从而为山地作业的定价细化提供计算方法。 一种山地作业难度指数计算方法，包括以下步骤： 步骤一：获取目标地区的遥感影像数据、DEM图和道路网数据； 步骤二：将上述图像的分辨率进行统一； 步骤三：将统一分辨率后的图像合成为一张多信息融合图，其中道路网数据转化 为栅格数据，道路在图像上表征为数值1，非道路在图像上表征为数值0； 步骤四：利用道路网图像、DEM图和遥感影像图，按照下面描述的方法进行山地作 业难度计算： A.生成一张道路网图像的复制图，在该图上沿着道路收集靠近道路的边缘点，即 该点的像素值为0，但至少有一个方向的像素点值为1，将上述边缘点构成一个集合S，以第 一个边缘点为选定点开始下述计算； B.生成基本算子B，基本算子B为一个3*3的模板，在道路网图像上，由选定点和周 4 CN 111582659 A 说 明 书 2/4 页 围的8个像素值生成基本算子，基本算子B为： 其中，数值1表示为道路，数值0表示为非道路，基本算子B表示道路的位置信息； C.生成距离算子D，距离算子D为一个3*3的模板： 1.4 1 1.4 1 0 1 1.4 1 1.4 距离算子D表示每个点距离中心的距离值； D.生成高差算子H，高差算子H为一个3*3的模板，在DEM图上，计算选定点和周围8 个像素点的高度差的绝对值并加1，将此值作为该点的值，从而生成高差算子H： h1 h2 h3 h4 0 h6 h7 h8 h9 E.计算选定点的难度指数，将基本算子B、距离算子D和高差算子H进行位乘运算， 计算结果表示为： d1 d2 d3 d4 0 d6 d7 d8 d9 那么选定点的难度指数表示为： d＝min{|di|} 其中di≠0 得到d值后，用该值替换道路网图像上选定点的值； F.继续计算集合S中的边缘点，直到集合S中所有的边缘点都计算完；在道路网图 像的复制图上，将集合S中所有边缘点的值置1，并重新生成一组边缘点集S，从步骤B开始计 算； G .所有非道路点都计算完成后，生成一张难度指数图，即得到该地区的山地作业 难度指数。 进一步的，步骤一遥感影像数据的来源不局限于卫星遥感数据和无人机遥感数 据。 进一步的，为了更方便后续的计算，步骤二中图像的分辨率统一到1米，即图像上 每1个像素代表1米长度。 本发明提出的一种基于多信息的山地作业难度指数计算方法，糅合了地理信息系 统和计算机图形学等技术，设计了三个算子：基本算子、距离算子和高差算子，综合考量了 距离、坡度和道路对山地作业难度的影响，从而可以逐块区域的评估该区域工作的难度，甚 至逐点的计算出山地上每点的工作难度值，为山地作业的定价细化提供了计算方法。该方 5 CN 111582659 A 说 明 书 3/4 页 法较传统的方式能够更精确的描述山地作业的难易程度，从而为合理定价打下了基础。该 方法综合考虑了距离、坡度、道路网格等多种因素，可以根据具体情况加入权值从而对计算 结果进行调整。 附图说明 图1为整体流程图； 图2为山地作业难度指数计算流程图。
本发明属于信息采集和处理技术领域，以解决山地作业难度评估方式不能全面地考虑作业点距离道路的远近、坡度大小、作业点分布的问题。山地作业难度指数计算方法：获取目标地区的数字图像；统一上述图像的分辨率并合成一张多信息融合图，道路网数据转化为栅格数据，道路 全部
背景技术：
山地作业的工作难度评估因其复杂、繁琐、计算量大等原因，长久以来一直没有得 到有效的解决。实际工作中一直沿用传统的一口价方式，这种估价方式较为简单，当工作量 较小时尚可，但是当工作量很大时，它的缺点就特别突出。主要有以下几个方面：(1)没有考 虑到作业点距离道路的远近。作业地点距离道路越远，作业人员需要花费的时间和精力就 越多，特别是当携带较重的工具的时候。(2)没有考虑作业点到道路坡度大小。在相同的距 离长度上，如山地的坡度越大，那么山路就越陡峭，同样其所消耗的时间和精力也越多。(3) 没有考虑作业点在山地上的分布情况。通常来说，采用统一定价的方式会综合考虑作业地 区山地的地形情况，按照平均难度给出个估计价格。但由于作业点的分布并不是完全均衡 的，有的时候会出现作业点大多靠近道路，有的时候会出现作业点大多远离道路，这就不满 足按照平均难度给出估价的初衷。(4)主观性强，没有确定的估价标准，结果往往因人而异。 传统的估价方式的弊端日趋明显，但是，山地作业的工作难度评估因其复杂、繁琐、计算量 大等原因，有效的切实可行的评估方法一直没有出现，上述问题也一直没有得到解决。鉴于 此，提出本发明。
技术实现要素：
现有的山地作业难度评估方式较为简单，不能够全面地考虑作业点距离道路的远 近、作业点到道路坡度大小、作业点分布情况，反映在工作定价上大多采用一口价的方式， 弊端明显。针对上述问题，本发明提出一种山地作业难度指数计算方法，综合利用了地理信 息系统和计算机图形学等技术，逐块区域的评估该区域工作的难度，逐点的计算出山地上 每点的工作难度值，从而为山地作业的定价细化提供计算方法。 一种山地作业难度指数计算方法，包括以下步骤： 步骤一：获取目标地区的遥感影像数据、DEM图和道路网数据； 步骤二：将上述图像的分辨率进行统一； 步骤三：将统一分辨率后的图像合成为一张多信息融合图，其中道路网数据转化 为栅格数据，道路在图像上表征为数值1，非道路在图像上表征为数值0； 步骤四：利用道路网图像、DEM图和遥感影像图，按照下面描述的方法进行山地作 业难度计算： A.生成一张道路网图像的复制图，在该图上沿着道路收集靠近道路的边缘点，即 该点的像素值为0，但至少有一个方向的像素点值为1，将上述边缘点构成一个集合S，以第 一个边缘点为选定点开始下述计算； B.生成基本算子B，基本算子B为一个3*3的模板，在道路网图像上，由选定点和周 4 CN 111582659 A 说 明 书 2/4 页 围的8个像素值生成基本算子，基本算子B为： 其中，数值1表示为道路，数值0表示为非道路，基本算子B表示道路的位置信息； C.生成距离算子D，距离算子D为一个3*3的模板： 1.4 1 1.4 1 0 1 1.4 1 1.4 距离算子D表示每个点距离中心的距离值； D.生成高差算子H，高差算子H为一个3*3的模板，在DEM图上，计算选定点和周围8 个像素点的高度差的绝对值并加1，将此值作为该点的值，从而生成高差算子H： h1 h2 h3 h4 0 h6 h7 h8 h9 E.计算选定点的难度指数，将基本算子B、距离算子D和高差算子H进行位乘运算， 计算结果表示为： d1 d2 d3 d4 0 d6 d7 d8 d9 那么选定点的难度指数表示为： d＝min{|di|} 其中di≠0 得到d值后，用该值替换道路网图像上选定点的值； F.继续计算集合S中的边缘点，直到集合S中所有的边缘点都计算完；在道路网图 像的复制图上，将集合S中所有边缘点的值置1，并重新生成一组边缘点集S，从步骤B开始计 算； G .所有非道路点都计算完成后，生成一张难度指数图，即得到该地区的山地作业 难度指数。 进一步的，步骤一遥感影像数据的来源不局限于卫星遥感数据和无人机遥感数 据。 进一步的，为了更方便后续的计算，步骤二中图像的分辨率统一到1米，即图像上 每1个像素代表1米长度。 本发明提出的一种基于多信息的山地作业难度指数计算方法，糅合了地理信息系 统和计算机图形学等技术，设计了三个算子：基本算子、距离算子和高差算子，综合考量了 距离、坡度和道路对山地作业难度的影响，从而可以逐块区域的评估该区域工作的难度，甚 至逐点的计算出山地上每点的工作难度值，为山地作业的定价细化提供了计算方法。该方 5 CN 111582659 A 说 明 书 3/4 页 法较传统的方式能够更精确的描述山地作业的难易程度，从而为合理定价打下了基础。该 方法综合考虑了距离、坡度、道路网格等多种因素，可以根据具体情况加入权值从而对计算 结果进行调整。 附图说明 图1为整体流程图； 图2为山地作业难度指数计算流程图。
