本发明涉及一种等高线负向地貌智能识别方法，更具体说它是一种用于国标地形图数据转换为军标数据进行地图出版时的等高线负向地貌(洼地)全自动判断方法。
背景技术：
采用国标地形图数据转换为军标数据进行地图出版时，需要对等高线负向地貌进行判断，并添加洼地示坡线。由于(1)国标数据中等高线没有负向属性项，人工进行判断工作量大且容易遗漏；(2)目前现有等高线判断和识别方法识别率不高，且人工参与比重较大；(3)现有生产软件中尚没有等高线负向地貌全自动判断的功能。
因此，为了在军用地形图生产中充分利用现有的国标地形图上数据，提高军标地形图出版的正确性和效率，进行等高线负向地貌自动判断十分必要。
技术实现要素：
本发明的目的是提供一种等高线负向地貌智能识别方法，计算和判断方法较快，识别率较高，大大降低了等高线负向地貌识别过程中人工参与的程度，满足国标地形数据转换为军标数据进行地图生产的应用。
为了实现上述目的，本发明的技术方案为：一种等高线负向地貌智能识别方法，其特征在于：包括如下步骤：
步骤1：进行等高线封闭处理；
步骤2：顺序查找封闭等高线，设为la；
步骤3：基于线缓冲区获取la相邻等高线lb；
步骤4：根据等高线la和lb的高程值，相互包含关系，以及lb的闭合性判断la和lb的类型；判断完后对下一条封闭等高线进行循环处理(循环步骤2至步骤4)；
步骤5：对存疑标记的等高线再进行如上判断(循环步骤2至步骤4)，仍不能确定存疑标记的等高线类型时进行人工判定，直至全图所有等高线判断完毕；实现等高线负向地貌的智能识别。
在上述技术方案中，步骤1中，当同高的等高线之间距离小于设定值(设定值一般设定为0.5mm)时，将同高的等高线进行连接或闭合；
当同高的等高线之间距离大于设定值，不进行处理，标记此等高线为不封闭。
在上述技术方案中，步骤2中，对所有等高线按照数据存储顺序进行标号，顺序查找等高线，若封闭，将此封闭的等高线设为la。
在上述技术方案中，步骤3包括如下内容：
(1)设置线缓冲区大小，构建la的线缓冲区；
(2)判断线缓冲区与其它等高线的相交关系，若线缓冲区内没有相交等高线，增大线缓冲区大小，继续进行判断；若线缓冲区内有多条相交等高线，减小线缓冲区大小，继续进行判断；直至线缓冲区内有且仅有一条相交等高线，将此相交等高线设为lb。
在上述技术方案中，步骤4中，等高线类型包括正向地貌，负向地貌和存疑待定。
在上述技术方案中，步骤4中，判断la和lb的类型包括如下内容：
(1)高程值la<lb
如果la不包含lb，则为负向地貌，la类型中输入负向；
如果la包含lb，则为正向地貌，la类型中输入正向，lb类型中输入正向；
(2)高程值la>lb
如果la不包含lb，则为正向地貌，la类型中输入正向；
如果la包含lb，则为负向地貌，la类型中输入负向，lb类型中输入负向；
(3)高程值la＝lb，继续判断
(ⅰ)lb为封闭等高线：
当lb包含la时，
若lb为正向地貌则la类型中输入负向；
若lb为负向地貌则la类型中输入正向；
当lb不包含la时，
若lb为正向地貌则la类型中输入正向；
若lb为负向地貌则la类型中输入负向；
lb为封闭等高线时，若lb的类型属性不明时，继续查找另一等高线，即跳转到步骤2，获取下一条等高线，设为lc，进行如下判断：
①当lc与la相邻，令lb＝lc，循环进行步骤3，当循环次数达到5次时，la类型中输入正向，并做存疑标记；循环结束；
②当lc与la不相邻时，
若lc闭合，
当la的高程低于lc的高程，la类型中输入正向；
当la的高程高于lc的高程，la类型中输入负向；
当la的高程等于lc的高程，lc为正向时，la类型中输入负向；lc为负向时，la类型中输入正向；
lc的类型未知时，跳转到步骤2，查找下一条等高线，循环所述步骤①、步骤②；
若lc不闭合，
当la的高程低于lc的高程，la类型中输入正向；
当la的高程高于lc的高程，la类型中输入负向；
当la的高程等于lc的高程，跳转到步骤2，查找下一条等高线，循环所述步骤①、步骤②；
(ⅱ)lb为不封闭等高线：跳转到步骤2，查找下一条等高线，循环步骤3，当判断lb为不封闭等高线的次数累计达到3次时，将la类型中输入正向，并做存疑标记；循环结束。
在上述技术方案中，步骤5中，当全图等高线处理完后，再对存在疑问的闭合等高线进行步骤2至步骤4的循环判断，仍不能确定时作出标记，人工判定；其中，图边等高线读入邻图等高线接边后进行判断。
本发明具有如下优点：
(1)计算和判断方法简单，充分利用了缓冲区分析技术，且速度快，通用性强；
(2)大大降低了等高线负向地貌识别过程中人工参与的程度，值得应用推广；
(3)满足了国标地形数据转换为军标数据进行地图生产的应用需求；本发明是一种用于国标地形图数据转换为军标数据进行地图出版时的等高线负向地貌全自动判断方法；克服了现有技术国标地形图数据转换为军标数据进行地图出版时，不能自动识别等高线负向地貌，并添加洼地示坡线的缺点。
附图说明
图1为本发明总体设计流程图。
图2为基于线缓冲区的相邻等高线查找示意图。
具体实施方式
下面结合附图详细说明本发明的实施情况，但它们并不构成对本发明的限定，仅作举例而已。同时通过说明使本发明的优点更加清楚和容易理解。
参阅附图可知：一种等高线负向地貌智能识别方法，其特征在于：包括如下步骤：
步骤1：进行等高线封闭处理；
步骤2：顺序查找封闭等高线，设为la；
步骤3：基于线缓冲区获取la相邻等高线lb；
步骤4：根据等高线la和lb的高程值，相互包含关系，以及lb的闭合性判断la和lb的类型；判断完后对下一条封闭等高线进行循环处理(循环步骤2至步骤4)；；
步骤5：对存疑标记的等高线再进行如上判断(循环步骤2至至步骤4)，仍不能确定存疑标记的等高线的类型时进行人工判定，直至全图所有等高线判断完毕；实现等高线负向地貌的智能识别。
步骤1中，当同高的等高线之间距离小于设定值时(设定值一般设定为0.5mm)，将同高的等高线进行连接或闭合；
当同高的等高线之间距离大于设定值，不进行处理，标记此等高线为不封闭，在后续等高线类型判断中需要根据等高线的闭合性进行分别判断。
步骤2中，对所有等高线按照数据存储顺序进行标号，顺序查找等高线，若封闭，将此封闭的等高线设为la。
步骤3包括如下内容：
(1)设置线缓冲区大小，构建la的线缓冲区；
(2)判断线缓冲区与其它等高线的相交关系，若线缓冲区内没有相交等高线，增大线缓冲区大小，继续进行判断；若线缓冲区内有多条相交等高线，减小线缓冲区大小，继续进行判断；直至线缓冲区内有且仅有一条相交等高线，将此相交等高线设为lb；由于等高线是由许多顶点连接构成的折线段，若利用常用的计算折线段之间的最近距离获取相邻线的方式需要遍历折线段上的所有点，耗时较长，无法满足实用性的要求；而利用等高线la的线缓冲区，通过判断与线缓冲区相交的等高线条数最终确定la的相邻等高线的方式，操作对象是线和面，判断效率要优于传统的距离计算方式，能够满足实际应用的需求。
步骤4中，等高线类型包括正向地貌，负向地貌和存疑待定。
步骤4中，判断la和lb的类型包括如下内容：
(1)高程值la<lb
如果la不包含lb，则为负向地貌，la类型中输入负向；
如果la包含lb，则为正向地貌，la类型中输入正向，lb类型中输入正向；
(2)高程值la>lb
如果la不包含lb，则为正向地貌，la类型中输入正向；
如果la包含lb，则为负向地貌，la类型中输入负向，lb类型中输入负向；
(3)高程值la＝lb，继续判断
(ⅰ)lb为封闭等高线：
当lb包含la时，
若lb为正向地貌则la类型中输入负向；
若lb为负向地貌则la类型中输入正向；
当lb不包含la时，
若lb为正向地貌则la类型中输入正向；
若lb为负向地貌则la类型中输入负向；
lb为封闭等高线时，若lb的类型属性不明时，继续查找另一等高线，即跳转到步骤2，获取下一条等高线，设为lc，进行如下判断：
①当lc与la相邻，令lb＝lc，循环进行步骤3，当循环次数达到5次时，la类型中输入正向，并做存疑标记；循环结束；
②当lc与la不相邻时，
若lc闭合，
当la的高程低于lc的高程，la类型中输入正向；
当la的高程高于lc的高程，la类型中输入负向；
当la的高程等于lc的高程，lc为正向时，la类型中输入负向；lc为负向时，la类型中输入正向；
lc的类型未知时，跳转到步骤2，查找下一条等高线，循环所述步骤①、步骤②；
若lc不闭合，
当la的高程低于lc的高程，la类型中输入正向；
当la的高程高于lc的高程，la类型中输入负向；
当la的高程等于lc的高程，跳转到步骤2，查找下一条等高线，循环所述步骤①、步骤②；
(ⅱ)lb为不封闭等高线：跳转到步骤2，查找下一条等高线，循环步骤3，当判断lb为不封闭等高线的次数累计达到3次时，将la类型中输入正向，并做存疑标记；循环结束。通常情况下，当循环判断相邻等高线为不封闭等高线的次数达到3次时，需要进行人工辅助判定，且将la的类型输入正向的目的是为了不影响后续等高线类型的判断。
步骤5中，当全图等高线处理完后，再对存在疑问的闭合等高线进行步骤2至步骤4的循环判断处理，仍不能确定时作出标记，人工判定；其中，图边等高线读入邻图等高线接边后进行判断。
其它未说明的部分均属于现有技术。