本发明公开了一种河曲自动提取方法，该方法包括：(1)读取一区间为[start,end]的矢量河流线图层数据，基于最短加权距离，提取出一曲流颈要素和相应的河曲要素，确定河曲区间的超始点位置,记为[neckStart,neckEnd]；(2)从数据区间[start,end]中剔除区间[neckStart,neckEnd]中间部分的数据，并进行子区间划分；(3)迭代执行步骤(1)至(2)，直至完成所有河曲及其曲流颈要素的提取；(4)从步骤(3)中得到的河曲及其曲流颈要素中提取符合用户指定阈值的部分，生成河曲线图层和曲流颈线图层。本发明执行效率高，准确率高。
1.一种河曲自动提取方法，其特征在于该方法包括：(1)读取一区间为[start,end]的矢量河流线图层数据，基于最短加权距离，提取出一曲流颈要素和相应的河曲要素，确定河曲区间的起始和终止点位置,记为[neckStart,neckEnd]；
(2)从数据区间[start,end]中剔除区间[neckStart,neckEnd]中间部分的数据，并进行子区间划分；
(3)迭代执行步骤(1)至(2)，直至完成所有河曲要素及其曲流颈要素的提取；
(4)从步骤(3)中得到的河曲要素及其曲流颈要素中提取符合用户指定阈值的部分，生成河曲线图层和曲流颈线图层。
2.根据权利要求1所述的河曲自动提取方法，其特征在于：步骤(1)具体包括：(1-1)读取一点数量为n，区间为[start,end]的矢量河流线图层数据，并写入集合L；
(1-3)计算pi和pj两点间加权距离dij；其中，pi和pj分别为位置i和j上的点；
(1-4)将i＝i+1,j＝j-1,返回执行(1-3)，直至i+2≥j时停止循环；
(1-5)循环结束后，获取加权距离最小的两点，其位置分别记为s、e；
(1-6)基于s、e所对应的两点，生成曲流颈要素；
(1-7)基于河曲区间公式，确定河曲区间的起始和终止点索引位置,记为[neckStart,neckEnd]；
(1-8)基于区间[neckStart,neckEnd]所圈定的河流分段，生成河曲要素。
3.根据权利要求2所述的河曲自动提取方法，其特征在于：步骤(1-3)中两点间加权距离的计算公式为：或者
式中，(xi，yi)、(xj，yj)分别为pi和pj两点的坐标，i、j分别为相应点在集合L中的位置,lenij指i,j两点间河流的长度。
4.根据权利要求2所述的河曲自动提取方法，其特征在于：步骤(1-7)中河曲区间公式为：式中，v为用户指定调整参数，其值域为[1,3]。
5.根据权利要求1所述的河曲自动提取方法，其特征在于：步骤(2)具体包括：(2-1)基于区间[neckStart,neckEnd]，根据以下公式计算子区间的划分位置，分别记为subLeft、subRight；
(2-2)从原始数据中剔除区间[neckStart,neckEnd]中间部分的数据，得到两个子区间[start,subLeft]和[subRight,end]。
6.根据权利要求1所述的河曲自动提取方法，其特征在于：步骤(4)具体包括：(4-1)针对一河曲要素，计算要素长度length和首尾点距离dis；
(4-3)循环执行步骤(4-1)至(4-2)，直至完成所有河曲要素的长度比系数计算；
(4-5)将剩余的河曲要素及其对应曲流颈要素生成河曲线图层和曲流颈线图层。
7.根据权利要求6所述的河曲自动提取方法，其特征在于：所述用户指定阈值R>1.3。
一种河曲自动提取方法
技术领域
[0001] 本发明涉及地理信息技术应用领域，尤其涉及一种基于矢量线状河流数据的河曲自动提取方法。
背景技术
[0002] 河流形态研究可为河流景观环境保护及水力资源开发提供地貌学依据，现代河流的几何形态参数对于古水文的重建与河流的演化研究至关重要。其中，河曲平面形态是河型研究的主要内容,与横向变化及纵向变化存在着相互制约的紧密关系。
[0003] 目前，自由河曲形态学及判别参数的研究已经取得了很大进展。然而，河曲形态学研究中，河曲的提取、分类以及河曲参数的测量等工作，自动化程度较低，仍以手工或半自动方法为主(参见张斌、艾南山、黄正文、易成波、覃发超.中国嘉陵江河曲的形态与成因.科学通报.2007,52(22):2671-2682)。如何基于矢量河流数据，自动、快速地划分并提取出河曲，对于河流形态研究具有重要的研究意义和应用价值。
发明内容
[0004] 发明目的：本发明针对现有技术存在的问题，提供一种基于矢量线状河流数据的河曲自动提取方法。
[0005] 技术方案：本发明所述的河曲自动提取方法包括：
[0006] (1)读取一区间为[start,end]的矢量河流线图层数据，基于最短加权距离，提取出一曲流颈要素和相应的河曲要素，确定河曲区间的起始和终止点位置,记为[neckStart,neckEnd]；start和end为提取区间的起点位置和终点位置；
[0007] (2)从数据区间[start,end]中剔除区间[neckStart,neckEnd]中间部分的数据，并进行子区间划分；
[0008] (3)迭代执行步骤(1)至(2)，直至完成所有河曲要素及其曲流颈要素的提取；
[0009] (4)从步骤(3)中得到的河曲要素及其曲流颈要素中提取符合用户指定阈值的部分，生成河曲线图层和曲流颈线图层。
[0010] 进一步的，步骤(1)具体包括：
[0011] (1-1)读取一点数量为n，区间为[start,end]的矢量河流线图层数据，并写入集合L；
[0012] (1-2)设置i＝start，j＝end；
[0013] (1-3)计算pi和pj两点间加权距离dij；其中，pi和pj分别为位置i和j上的点；
[0014] (1-4)将i＝i+1,j＝j-1,返回执行(1-3)，直至i+2≥j时停止循环；
[0015] (1-5)循环结束后，获取加权距离最小的两点，其位置分别记为s、e；
[0016] (1-6)基于s、e所对应的两点，生成曲流颈要素；
[0017] (1-7)基于河曲区间公式，确定河曲区间的起始和终止索引位置,记为[neckStart,neckEnd]；
[0018] (1-8)基于区间[neckStart,neckEnd]所圈定的河流分段，生成河曲要素。
[0019] 其中，步骤(1-3)中两点间加权距离的计算公式为：
[0020]
[0021] 或者
[0022]
[0023] 式中，(xi，yi)、(xj，yj)分别为pi和pj两点的坐标，i、j分别为相应点在集合L中的位置，lenij指i,j两点间河流的长度。
[0024] 其中，步骤(1-7)中河曲区间公式为：
[0025]
[0026] 式中，v为用户指定调整参数，其值域为[1,3]。
[0027] 进一步的，步骤(2)具体包括：
[0028] (2-1)基于区间[neckStart,neckEnd]，根据以下公式计算子区间的划分位置，分别记为subLeft、subRight；
[0029]
[0030] 式中，round函数用于四舍五入求值；
[0031] (2-2)从原始数据中剔除区间[neckStart,neckEnd]中间部分的数据,得到两个子区间[start,subLeft]和[subRight,end]；
[0032] 进一步的，步骤(4)具体包括：
[0033] (4-1)针对一河曲要素，计算要素长度length和首尾点距离dis；
[0034] (4-2)基于下式，计算其长度比系数ratio；
[0035]
[0036] (4-3)循环执行步骤(4-1)至(4-2)，直至完成所有河曲要素的长度比系数计算；
[0037] (4-4)基于用户指定阈值R，删除ratio
[0038] (4-5)将剩余的河曲要素及其对应曲流颈要素生成河曲线图层和曲流颈线图层。
[0039] 其中，所述用户指定阈值R>1.3。
[0040] 有益效果：本发明与现有技术相比，其显著优点是：本发明基于最短加权距离法和二分区间迭代处理算法，执行效率较高，提取效果好，与人工解译结果一致。
附图说明
[0041] 图1是本发明的一个实施例的流程示意图；
[0042] 图2是实施例中采用的实验数据图；
[0043] 图3是实施例提取的河曲示意图；
[0044] 图4是采用本发明方法基于阈值过滤后生成的河曲及其曲流颈示意图(虚线要素为曲流颈)。
具体实施方式
[0045] 下面对本发明技术方案进行详细说明，如图1所示，本实施例具体包括步骤：
[0046] 步骤1：读取一区间为[start,end]的矢量河流线图层数据，基于最短加权距离，提取出一曲流颈要素和相应的河曲要素，确定河曲区间的起始和终止点位置,记为[neckStart,neckEnd]；start和end为提取区间的起点位置和终点位置。
[0047] 其中，本实例的实验数据采用的是由Google Earth中遥感影像数字化生成的1:500000嘉陵江局部线图层数据(图2)。
[0048] 步骤1具体包括以下步骤：
[0049] (1-1)读取一点数量为n，区间为[start,end]的矢量河流线图层数据，并写入集合L；本实施例中，n为497；
[0050] (1-2)设置i＝start，j＝end；
[0051] (1-3)计算pi和pj两点间加权距离dij；本实施例使用计算公式为：
[0052]
[0053] 式中，(xi，yi)、(xj，yj)分别为pi和pj两点的坐标，i、j分别为相应点在集合L中的位置；
[0054] (1-4)将i＝i+1,j＝j-1,返回执行(1-3)，直至i+2≥j时停止循环；
[0055] (1-5)循环结束后，获取加权距离最小的两点，其位置分别记为s、e；本实施例中，首次得到s、e的值分别为416、472；
[0056] (1-6)基于s、e所对应的两点，生成曲流颈要素；
[0057] (1-7)基于河曲区间公式，确定河曲区间的起始和终止点索引位置,记为[neckStart,neckEnd]；
[0058] 其中，河曲区间公式为：
[0059]
[0060] 式中，v为用户指定调整参数，其值域为[1,3]。本实施例中，v＝1,首次得到区间为[415,473]；
[0061] (1-8)基于区间[neckStart,neckEnd]所圈定的河流分段，生成河曲要素。
[0062] 步骤2：从原始数据中剔除区间[neckStart,neckEnd]中间部分的数据，并进行子区间划分。具体包括：
[0063] (2-1)基于区间[neckStart,neckEnd]，根据以下公式计算子区间的划分位置，分别记为subLeft、subRight；
[0064]
[0065] 式中，round函数用于四舍五入求值；本实施例中，首次得到subLeft＝430、subRight＝458；
[0066] (2-2)从原始数据中剔除区间[neckStart,neckEnd]中间部分的数据,得到两个子区间[start,subLeft]和[subRight,end]；本实施例中，首次得到两个子区间分别为[0，430]、[458，496]。
[0067] 步骤3：迭代执行步骤1至2，直至完成所有河曲要素及其曲流颈要素的提取。本实施例中，处理结果如图3所示。
[0068] 步骤4：从步骤3中得到的河曲要素及其曲流颈要素中提取符合用户指定阈值的部分，生成河曲线图层和曲流颈线图层。具体包括：
[0069] (4-1)针对一河曲要素，计算要素长度length和首尾点距离dis；本实施例中，首个要素的length和dis,其值分别为0.22353592375751224、0.010305162601980757；
[0070] (4-2)基于下式，计算其长度比系数ratio；
[0071]
[0072] 本实施例中，首次得到ratio＝21.69164450782624；
[0073] (4-3)循环执行步骤(4-1)至(4-2)，直至完成所有河曲要素的长度比系数计算；
[0074] (4-4)基于用户指定阈值R，删除ratio
[0075] (4-5)将剩余的河曲要素及其对应曲流颈要素生成河曲线图层和曲流颈线图层。本实施例中，处理结果如图4所示。
[0076] 由实施例可知，本专利方法能有效地进行河曲和曲流颈提取，自动化程度高，且与人工提取结果基本一致。
[0077] 以上所揭露的仅为本发明一种较佳实施例而已，不能以此来限定本发明之权利范围，因此依本发明权利要求所作的等同变化，仍属本发明所涵盖的范围。
