[发明专利]一种基于山脊能量校正的山地中带状地下目标的探测方法

权利要求书

1.一种基于山脊能量校正的山地中带状地下目标的探测方法，其特征在于，所述方法包括如下步骤：

(1)获取地形数字高程信息，包括下述子步骤：

(1.1)确定数字高程信息覆盖范围的经纬度；

(1.2)根据上述经纬度计算数字高程信息覆盖范围内经纬度矩阵；

(1.3)根据上述经纬度矩阵计算经纬度矩阵中每个点的海拔，得到数字高程信息矩阵；

(2)对步骤(1)获得的数字高程信息矩阵进行去噪预处理；

(3)对降噪处理后的数字高程信息矩阵进行山脊线检测步骤；包括下述子步骤：

(3.1)沿水平方向检测山脊线；

(3.2)沿垂直方向检测山脊线；

(3.3)根据上述水平方向检测的山脊线和垂直方向检测的山脊线，提取连续的山脊线；

(4)山脊位置能量校正步骤；根据山脊位置能量分布特征对山脊位置能量进行校正；

(5)在数字高程信息矩阵中进行带状地下目标探测的步骤；包括下述子步骤：

(5.1)设置遍历探测带状地下目标所用参数；

(5.2)根据上述参数对数字高程信息矩阵进行遍历探测，并将探测的带状地下目标位置输出。

2.如权利要求1所述的方法，其特征在于，所述步骤(1.2)具体包括：

(1.2.1)计算地形范围的宽度、高度步骤：

根据步骤(1.1)中确定的矩形四个顶点的经纬度利用Google Earth提供的测距工具计算出矩形覆盖范围的宽度、高度标记为width、height单位为米，并计算出经纬度的范围；

(1.2.2)计算经纬度矩阵的步骤：

按照采样间隔step米,每次距离step米采样一次，分别计算出纵向、横向的采样点的数目，即：height/step，width/step，并计算出经纬度矩阵每相邻两个采样点的经度步长long_step和纬度步长lati_step，从而计算出经纬度矩阵中每一个采样点的经纬度。

3.如权利要求1或2所述的方法，其特征在于，所述步骤(1.3)具体包括：

将(1.2.2)中经纬度矩阵作为输入利用Google Earth得到每一个采样点的海拔数据，并根据该海拔数据生成地形数字高程信息矩阵输出。

4.如权利要求1至3任一项所述的方法，其特征在于，所述步骤(2)具体为：采用均值滤波方法对上述数字高程信息进行去噪预处理。

5.如权利要求1至3任一项所述的方法，其特征在于，所述步骤(3.1)具体为：

沿着x方向遍历，将每一个采样点的高程值与y方向预设范围内的采样点的高程值比较，如果该采样点在y方向是极大值，则该点为水平方向的山脊线的候选点。

6.如权利要求1至3任一项所述的方法，其特征在于，所述步骤(3.2)具体为：

沿着y方向遍历，将每一个采样点的高程值与x方向预设范围内的采样点的高程值比较，如果该采样点在x方向是极大值，则该点为垂直方向的山脊线的候选点。

7.如权利要求1至3任一项所述的方法，其特征在于，所述步骤(3.3)具体为：

对步骤(3.1)、(3.2)提取出的每一个山脊线候选点进行连续性判别，判别准则如下：如果该候选山脊线点的t*t邻域内不同的山脊线候选点，则该山脊线候选点为最终的山脊线点；否则，判断该孤立的山脊线候选点为非山脊线点，最终得到连续山脊线标记图，其中t为预设值。

8.如权利要求1至3任一项所述的方法，其特征在于，所述步骤(4)具体为：对山脊位置的能量进行校正，校正的具体方法为：山脊线位置采样点的能量用山脊线位置两侧的采样点的能量的均值代替。