[发明专利]一种山地作业难度指数计算方法

权利要求书

1.一种山地作业难度指数计算方法，其特征在于，包括以下步骤：

步骤一：获取目标地区的遥感影像数据、DEM图和道路网数据；

步骤二：将上述图像的分辨率进行统一；

步骤三：将统一分辨率后的图像合成为一张多信息融合图，其中道路网数据转化为栅格数据，道路在图像上表征为数值1，非道路在图像上表征为数值0；

步骤四：利用道路网图像、DEM图和遥感影像图，按照下面描述的方法进行山地作业难度计算：

A.生成一张道路网图像的复制图，在该图上沿着道路收集靠近道路的边缘点，即该点的像素值为0，但至少有一个方向的像素点值为1，将上述边缘点构成一个集合S，以第一个边缘点为选定点开始下述计算；

B.生成基本算子B，基本算子B为一个3*3的模板，在道路网图像上，由选定点和周围的8个像素值生成基本算子，基本算子B为：

其中，数值1表示为道路，数值0表示为非道路，基本算子B表示道路的位置信息；

C.生成距离算子D，距离算子D为一个3*3的模板：

距离算子D表示每个点距离中心的距离值；

D.生成高差算子H，高差算子H为一个3*3的模板，在DEM图上，计算选定点和周围8个像素点的高度差的绝对值并加1，将此值作为该点的值，从而生成高差算子H：

E.计算选定点的难度指数，将基本算子B、距离算子D和高差算子H进行位乘运算，计算结果表示为：

选定点的难度指数表示为：

d＝min{|di|} 其中di≠0

得到d值后，用该值替换道路网图像上选定点的值；

F.继续计算集合S中的边缘点，直到集合S中所有的边缘点都计算完；在道路网图像的复制图上，将集合S中所有边缘点的值置1，并重新生成一组边缘点集S，从步骤B开始计算；

G.所有非道路点都计算完成后，生成一张难度指数图，即得到该地区的山地作业难度指数。

2.根据权利要求1所述的计算方法，其特征在于，步骤二中图像的分辨率统一设置为1米。