1.一种峡谷地貌实体的自动化识别方法，其特征在于该方法包括：(1)从水系图层中读取一水系要素，创建缓冲区，并提取缓冲区在水系要素两侧的边界线leftline、rightline；

(2)分别采用边界线leftline、rightline对DEM作堆栈剖面，从而得到边界线所经过的格网单元的高程序列LE、RE；

(3)分别提取高程序列LE、RE的波峰集合LP、RP；

(4)对于波峰集合LP、RP中的波峰进行配对，并基于配对的两个波峰所确定区域，生成一个外包络矩形，所述外包络矩形即为识别出的一个峡谷所在区域；

(5)返回循环执行步骤(1)到(4)，直至水系图层中每一水系要素都被处理；

(6)基于所生成的外包络矩形集合，生成峡谷实体专题图层。

2.根据权利要求1所述的峡谷地貌实体的自动化识别方法，其特征在于：所述步骤(1)具体包括：(1-1)从水系图层中读取一水系要素wi，其中，i∈{1,2,…，N}，N为水系图层中水系线要素的个数；

(1-2)创建一距离为d的缓冲区buffer；

(1-3)提取buffer在水系要素wi两侧的边界线，分别记为leftline、rightline。

3.根据权利要求1所述的峡谷地貌实体的自动化识别方法，其特征在于：步骤(2)具体包括：(2-1)采用边界线leftline对DEM作堆栈剖面，从而得到高程序列LE＝{li|i＝1,2,…,LN}，其中，li为边界线leftline所经过的第i个格网单元的高程值，LN为序列LE元素的个数；

(2-2)采用边界线rightline对DEM作堆栈剖面，从而得到高程序列RE＝{rj|j＝1,2,…,RN}，其中，rj为边界线rightline所经过的第j个格网单元的高程值，RN为序列RE元素的个数。

4.根据权利要求1所述的峡谷地貌实体的自动化识别方法，其特征在于：步骤(3)具体包括：(3-1)分别计算序列LE、RE的包络线序列LEV＝{levi|i＝1,2,…,LN}、REV＝{revj|j＝

1,2,…,RN}；

(3-2)根据以下公式计算LE的累计均值线序列LA＝{lai|i＝1,2,…,,LN}，以及RE的累计均值线序列RA＝{raj|j＝1,2,…,RN}：式中，li和rj分别为序列LE、RE的元素；

(3-3)获取包络线序列LEV在累计均值线序列LA上方的序列LEV’＝{lev’i|i＝1,2,…,LN}，以及包络线序列REV在累计均值线序列RA上方的序列REV’＝{rev’j|j＝1,2,…,RN}；

其中：

(3-4)分别提取出序列LEV’、REV’中波峰高度大于设定阈值PT的波峰，并分别加入波峰的三元组集合LP＝({lpi|i＝1,2,…,LM}、RP＝{rpj|j＝1,2,…,RM}中；其中，波峰lpi＝(lpi(start),lpi(peak),lpi(end))，波峰rpj＝(rpj(start),rpj(peak),rpj(end))，start为对应波峰左侧平衡位置、peak为对应波峰位置、end为对应波峰右侧的平衡位置，LM、RM分别为序列LEV’、REV’中波峰高度大于PT的波峰的个数。

5.根据权利要求4所述的峡谷地貌实体的自动化识别方法，其特征在于：步骤(4)具体包括：(4-1)若集合LP或RP均非空，执行步骤(4-2)；否则，执行步骤(5)；

(4-2)在集合RP中，查找与元素lpi配对的波峰rpnearest；

(4-3)针对lpi的三个元素lpi(start)、lpi(peak)、lpi(end)和其配对波峰rpnearest的三个元素rpnearest(start)、rpnearest(peak)、rpnearest(end)，在每个元素所对应网格单元的中心位置创建一个点；

(4-4)生成6个点的外包络矩形，所述外包络矩形即为识别出的一个峡谷所在区域；

(4-5)返回循环执行步骤(4-2)至(4-4)，直至集合LP中所有元素都被配对处理。

6.根据权利要求5所述的峡谷地貌实体的自动化识别方法，其特征在于：步骤(4)具体包括：步骤(4-2)具体包括：(4-2-1)计算lpi(peak)与rpj(peak)的差值nearj，并存入集合Near＝{nearj|j＝1,

2,…,RM}；

(4-2-2)获取集合Near中的最小值nearnearest，若nearnearest小于设定阈值NT，则判定波峰rpnearest与lpi配对，即rpnearest与lpi位于同一峡谷两侧，并从集合RP中剔除rpnearest，执行步骤(4-3)；否则，执行步骤(4-5)。