CN114913263B - 一个基于多尺度空间相似度的线状地物自动化简方法 - Google Patents
一个基于多尺度空间相似度的线状地物自动化简方法 Download PDFInfo
- Publication number
- CN114913263B CN114913263B CN202110175507.0A CN202110175507A CN114913263B CN 114913263 B CN114913263 B CN 114913263B CN 202110175507 A CN202110175507 A CN 202110175507A CN 114913263 B CN114913263 B CN 114913263B
- Authority
- CN
- China
- Prior art keywords
- similarity
- epsilon
- scale
- reduction
- function
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 39
- 230000009467 reduction Effects 0.000 claims abstract description 32
- 230000008859 change Effects 0.000 claims abstract description 25
- 238000004364 calculation method Methods 0.000 claims description 9
- 238000012888 cubic function Methods 0.000 claims 1
- 238000012886 linear function Methods 0.000 claims 1
- 238000012887 quadratic function Methods 0.000 claims 1
- 238000010276 construction Methods 0.000 description 7
- 230000015572 biosynthetic process Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 238000003786 synthesis reaction Methods 0.000 description 5
- 238000005259 measurement Methods 0.000 description 2
- 238000011946 reduction process Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 238000007429 general method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/20—Drawing from basic elements, e.g. lines or circles
- G06T11/203—Drawing of straight lines or curves
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Processing Or Creating Images (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
本发明公开一个基于多尺度空间相似度的线状地物自动化简方法,包括:地图比例尺变化与曲线目标相似度的函数关系推导方法、曲线目标相似度与Douglas‑Peucker化简法阈值ε的函数关系推导方法,并最终得出ε与比例尺变化的函数关系,由此实现线状要素化简的全自动化。该方法不仅实现了Douglas‑Peucker算法的全自动化,而且其对曲线的化简结果与经验丰富的制图员的手工化简结果具有高度的相似性。
Description
技术领域
本发明涉及地图自动综合领域,更具体的说是一个基于多尺度空间相似度的高效的线状地物自动化简方法。
背景技术
地图上80%以上的要素(等高线、境界线、河流、道路、地类界等)以曲线形式表达。当地图由大比例尺变为较小比例尺时,就需要对曲线要素进行化简,以满足地图空间线状目标可视化的要求。曲线要素化简属于地图综合的操作之一。曲线要素化简的目的是在尽量保持化简后曲线与原始曲线形状相似的情况下,使化简后曲线的复杂程度与目标地图的比例尺相适应。为了实现地图上曲线的化简,学者们提出了许多化简方法,如临近点简化化简法、垂直距离化简法、曲率变化化简法、Douglas-Peucker化简法等。其中,在地图综合中应用比较广泛的是Douglas-Peucker化简法,以下简称其为DP化简法。
在地图综合中,理想的曲线化简法应该是全自动化的,即曲线的化简过程不受人工干预。但是,DP化简法并非一个全自动方法,因为其在曲线化简开始时需要人工输入一个距离阈值(称为ε,ε>0)。已有研究者陆续提出了改进的DP方法,例如适用于闭合曲线化简的DP化简法、可以处理化简后曲线自相交问题的DP化简法等,但它们大都聚焦于使用场景普适性的提高,并没有顾及到DP方法的非全自动化问题。Douglas-Peucker化简法非自动化的原因在于曲线化简前需要人工输入距离阈值ε。因此,解决此问题的自然设想是:(1)在程序执行前获得ε的值,把它嵌入到化简程序中;或者,(2)找到计算ε的公式,把它编写入化简程序。
对(1)进行分析:在面对特定的地图综合任务时,根据已知的原始地图比例尺和目标地图比例尺,询问有经验的地图专家,可以得到ε的值。但是,在地图数据的实际生产中,ε的确定需要顾及数量众多的地图比例尺和多样的制图区域地理特点,故很难全部枚举出可能的ε。因此,如果能按照(2)的思路,找到计算ε的公式,就可能实现化简的自动化。
基于上述原因,本发明将专注于地图上曲线状地物、地貌的化简方法,目标是提出一个全自动的Douglas-Peucker化简法。
发明内容
本发明提供了一个基于多尺度空间相似度的线状地物自动化简方法,它可以自动计算化简阈值,以实现DP化简法的全自动化。图1为本发明基于多尺度空间相似度的线状地物自动化简方法的总体流程,包括相似度与比例尺变化关系的构建、相似度与DP阈值关系的构建、比例尺变化与DP阈值关系的构建三部分。
设有一个地图综合任务,需要化简n条曲线要素,原始地图比例尺为S0,目标地图比例尺为S1。要实现DP化简的全自动化,就需要找到ε与地图比例尺变化S(S=S0/S1)之间的函数关系:
ε=f1(S) (1)
地图综合是一种相似变换,DP化简法对曲线的化简即是曲线在不同比例尺地图上的相似变化,化简后曲线与原始曲线的相似度(用Sim表示)和地图比例尺变化S之间显然存在单调函数依赖:曲线化简中,S越大,则Sim越小。所以,如果求得如下关系,就可以推导出f1。
Sim=f2(S) (2)
Sim=f3(ε) (3)
依据上述需求,本发明提供了计算ε与地图比例尺变化S函数关系的一般性操作方法:首先给出多尺度曲线相似度的计算方法,然后结合具体实验数据拟合比例尺变化与相似度的函数关系,随后拟合相似度与ε的函数关系,随后确定ε与比例尺变化的函数关系,并最终得出在不同比例尺变化S下阈值ε的值。
本发明公开提出的基于多尺度空间相似度的线状地物自动化简方法易于实现,简单高效,不仅实现了Douglas-Peucker化简法的全自动化,而且其对曲线的化简结果与经验丰富的制图员的手工化简结果具有高度的相似性。
附图与表格说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图与表格作简单介绍,显而易见地,下面描述中的附图仅仅是本发明的示意图与表格,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图获得其他的附图与表格。
图1为本发明提供的自动化简方法的流程图。
图2为最大Hausdorff距离计算示意图。
图3为1:10K河流数据示意图。
图4为3条河流化简结果示意图。
图5为比例尺变化与相似度的函数关系拟合示意图。
图6为推导ε所用的河流数据示意图。
图7为相似度与ε的函数关系拟合示意图。
图8为本发明实验结果及对比示意图。
表1为本发明得出的DP化简法最佳阈值。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
以下步骤为相似度与比例尺变化关系的构建部分:
Step 1:计算单一线要素的Hausdorff距离。设比例尺S0的地图上有n个线要素集化简到比例尺为S1的地图上得到的对应要素集/> 设原始曲线/>与化简曲线/>的构成点集为:A={a1,a2,...,am},B={b1,b2,...,bn},则其Hausdorff距离H(A,B)定义为∶
H(A,B)=max{h(A,B),h(B,A)} (4)
其中,
h(A,B)=max(a∈A)min(b∈B)||a-b|| (5)
h(B,A)=max(b∈B)min(a∈A)||b-a|| (6)
其中,||·||表示两点的距离范式,可以用欧氏距离、曼哈顿距离、切比雪夫距离等任何一种距离,这里它代表欧氏距离。
Step 2:计算单一线要素的最大Hausdorff距离。MaxH是与/>构成的度量空间内可能出现的最大Hausdorff距离,计算方法如下:由于DP化简法的特征是,化简过程不改变曲线首尾点,故可将Hausdorff距离的度量空间限定到原始曲线/>与化简曲线/>的最小面积外接矩形MABR(Minimum Area Boundary Rectangle)。由此容易得到为MABR的两条边中的较长者。如图2,ABC为原始曲线,A′C′为化简后曲线,ABC的MABR如虚点线所示,则MaxH(ABC,A′C′)显然为M2。
Step 3:计算单一线要素的相似度。根据上述单一线要素的Hausdorff距离与最大Hausdorff距离的计算方法,我们给出和/>的相似度计算方法:
Step 4:计算n条线要素的综合相似度。根据上述单一线要素的相似度计算方法,我们给出和/>的相似度计算方法:
其中,是/>的长度,/>是原始n条线要素的总长。
Step 5:构建多尺度线状要素数据库。为了得到比例尺变化和相似度之间较为普适的函数关系,需要构建多尺度线状要素数据库。这里给出一个示例:选取如图3所示的某地区1幅1∶10000矢量地形图上的15条线状河流来构建多尺度线状要素数据库。由经验丰富的制图员对这些河流分别进行化简,得到了它们在1∶25000、1∶50000、1∶100000、1∶250000、1∶500000、1∶1000000地图上的河流图形。图4给出了其中3个河流要素的在7种比例尺地图上的化简图形(篇幅所限,其余12条河流的化简图形从略)。
Step 6:推导比例尺变化与相似度的函数关系。对于每一个河流要素,分别计算原始河流(1∶10000)与化简后每个比例尺河流图形之间的相似度。这样,可得到15条河流在7种比例尺(1∶10000图形可以与原始河流自身进行相似度计算)下图形的105个相似度值。因为每个相似度(Sim)都与一个比例尺变化(S:原始图形比例尺与化简后图形比例尺的比值)相对应,所以用(S,Sim)组成了105个坐标对作为数据拟合点。考虑到相似度与比例尺变化之间的单调函数依赖关系,选取如下6个候选函数,用得到的105个坐标对作为数据拟合点,来拟合和寻找它们之间最合适的函数关系:
1.Sim=aS+b; 2.Sim=aS2+bS+c;
3.Sim=aS3+bS2+cS+d; 4.Sim=aln(S+b)+c;
5.Sim=aebS; 6.Sim=aSb+b。
拟合结果如图5所示。其中,对数函数的R2=0.92是最大的,表明其函数拟合的结果最好,故此选择对数函数作为比例尺变化与相似度的函数关系:
Sim=2.122-0.275ln(S+60.764) (9)
需要说明的是,本发明在这里给出的是一个普适性的相似度与比例尺变化关系的构建方式。发明的使用者,既可以直接使用公式(9),亦可以根据实际情况自主建库来推导二者的函数关系。
以下步骤为相似度与DP阈值关系的构建部分:
Step 7:构建多阈值DP化简数据库。下面以某特定地理区域1:10K地图水系要素的化简为例,推导其相似度与ε的函数关系。如图6所示,地图上有16条河流,相连的河流之间以圆点为分界。根据DP化简法原理可以计算得到把所有河流均化简为直线时的ε值,此即为ε最大值(本例为678.3m)。为了得到可信度较高的函数关系,取ε步长为0.5m,即ε={0.5m,1.0m,1.5m,...,678.5m},对图6中的河流数据运用DP化简法进行1356次化简,得到原始数据在各个ε下的化简结果。
Step 8:计算原始河流数据与每个ε值下化简后得到的河流数据的相似度Sim,由此得到1356个(ε,Sim)坐标对。对这些坐标对进行曲线拟合,结果如图6。其中,对数函数拟合的结果最好,R2=0.989,所以其对应的函数被认为是针对本实验数据的距离阈值与相似度的最佳函数关系:
Sim=-0.096ln(ε)+1.2097 (10)
公式(10)的反函数也可以得到:
ε=e12.601-10.471Sim (11)
需要说明的是,选取的实验数据不同,ε与相似度Sim的函数关系一般都会不同。
以下步骤为比例尺变化与DP阈值关系的构建部分:
Step 9:结合公式(9)、公式(11),可以得到地图比例尺变化S与阈值ε之间的函数关系:
ε=e12.601-10.471×(2.122-0.275ln(S+60.764)) (12)
对其进行化简可得:
ε=e0.275ln(S+60.764)-9.618 (13)
这样,就得到了针对图5数据的阈值ε与比例尺变化S之间的函数关系式。
Step 10:根据公式(13),可以计算在任何比例尺变化下,适宜的DP化简的阈值ε,下面给出一个示例:欲将图5中1:10K河流数据用DP化简法自动化简为1∶25K、1∶50K、1∶100K、1∶250K、1∶500K、1∶1M共6个小比例尺数据。为此,用公式(13)计算得出ε,见表1。如图7是本发明提出的全自动化简方法结果、有经验制图员手工化简结果、二者的叠置对比图。
表1 DP化简法最佳阈值
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本发明中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本发明所示的这些实施例,而是要符合与本发明所公开的原理和新颖特点相一致的最宽的范围。
Claims (3)
1.一个基于多尺度空间相似度的线状地物自动化简方法,其特征在于,包括:
S1:计算单一线要素的Hausdorff距离,记为H;
S2:计算单一线要素的最大Hausdorff距离,记为MaxH;
S3:使用公式计算单一线要素的相似度,其中/>为原始曲线,/>为化简曲线;
S4:使用计算n条线要素的综合相似度,其中/>是/>的长度,/>是原始n条线要素的总长;
S5:构建多尺度线状要素数据库;
S6:使用曲线拟合的方式,推导比例尺变化S与步骤S4中计算得到的相似度的函数关系,其中曲线拟合的方式可以是线性函数、二次函数、三次函数、幂函数、指数函数和对数函数;
S7:使用Douglas-Peucker算法计算把所有线状要素化简为直线时的化简阈值ε的值,得到原始数据在各个ε下的化简结果,并将各个化简结果的集合构建成数据库;
S8:使用计算原始河流数据与每个ε值下化简后得到的河流数据的相似度Sim,得到形如(ε,Sim)的坐标对,并使用与S6相同形式的函数进行曲线拟合;
S9:结合S6中得到的S与的函数关系、S8中得到的ε与Sim的函数关系,推导地图比例尺变化S与阈值ε之间的函数关系;
S10:使用S9中得到的S与ε的关系,计算在任何比例尺变化下,适宜的DP化简的阈值ε。
2.根据权利要求1所述的一个基于多尺度空间相似度的线状地物自动化简方法,其特征在于:在步骤S1到S6中,使用Hausdorff距离、最大Hausdorff距离,结合提出的单一线要素相似度计算公式综合线要素相似度计算公式推导比例尺变化与相似度的函数关系。
3.根据权利要求1或2所述的一个基于多尺度空间相似度的线状地物自动化简方法,其特征在于,在步骤S7到S9中,构建多阈值DP化简数据库,计算多阈值数据的相似度,结合权利要求2的过程,并得到地图比例尺变化S与阈值ε之间的函数关系。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110175507.0A CN114913263B (zh) | 2021-02-06 | 2021-02-06 | 一个基于多尺度空间相似度的线状地物自动化简方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110175507.0A CN114913263B (zh) | 2021-02-06 | 2021-02-06 | 一个基于多尺度空间相似度的线状地物自动化简方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114913263A CN114913263A (zh) | 2022-08-16 |
CN114913263B true CN114913263B (zh) | 2023-11-14 |
Family
ID=82762144
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110175507.0A Active CN114913263B (zh) | 2021-02-06 | 2021-02-06 | 一个基于多尺度空间相似度的线状地物自动化简方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114913263B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101707026A (zh) * | 2009-11-25 | 2010-05-12 | 中国人民解放军信息工程大学 | 数字地图线状要素化简的组合优化方法 |
JP2013174796A (ja) * | 2012-02-27 | 2013-09-05 | Toyota Mapmaster:Kk | 処理単位指定装置及びその方法、並びに処理単位を指定するためのコンピュータプログラム及びコンピュータプログラムを記録した記録媒体 |
CN104978763A (zh) * | 2015-05-13 | 2015-10-14 | 中国矿业大学(北京) | 一种基于三维Douglas-Peucker算法的河网要素与DEM的同步综合地图仿真方法 |
JP2016114970A (ja) * | 2014-12-11 | 2016-06-23 | 株式会社ゼンリン | 地図データ生成装置 |
CN106845108A (zh) * | 2017-01-19 | 2017-06-13 | 中国测绘科学研究院 | 一种顾及空间关系约束的线化简方法 |
CN108122202A (zh) * | 2018-01-05 | 2018-06-05 | 武汉大学 | 一种基于形状上下文的地图线状要素多尺度信息派生方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2007014375A2 (en) * | 2005-07-26 | 2007-02-01 | Decarta, Inc. | Generalization of features in a digital map |
-
2021
- 2021-02-06 CN CN202110175507.0A patent/CN114913263B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101707026A (zh) * | 2009-11-25 | 2010-05-12 | 中国人民解放军信息工程大学 | 数字地图线状要素化简的组合优化方法 |
JP2013174796A (ja) * | 2012-02-27 | 2013-09-05 | Toyota Mapmaster:Kk | 処理単位指定装置及びその方法、並びに処理単位を指定するためのコンピュータプログラム及びコンピュータプログラムを記録した記録媒体 |
JP2016114970A (ja) * | 2014-12-11 | 2016-06-23 | 株式会社ゼンリン | 地図データ生成装置 |
CN104978763A (zh) * | 2015-05-13 | 2015-10-14 | 中国矿业大学(北京) | 一种基于三维Douglas-Peucker算法的河网要素与DEM的同步综合地图仿真方法 |
CN106845108A (zh) * | 2017-01-19 | 2017-06-13 | 中国测绘科学研究院 | 一种顾及空间关系约束的线化简方法 |
CN108122202A (zh) * | 2018-01-05 | 2018-06-05 | 武汉大学 | 一种基于形状上下文的地图线状要素多尺度信息派生方法 |
Non-Patent Citations (5)
Title |
---|
Wang Xiao-li,etc.Selecting Optimal Threshold Value of Douglas-Peucker Algorithm Based on Curve Fit.《2010 First International Conference on Networking and Distributed Computing》.2010,全文. * |
基于三维DP算法的河网线矢量与DEM同步简化方法;窦世卿;《中国博士学位论文全文数据库基础科学辑》;A008-3 * |
基于交叉点的道路曲线化简算法研究;李世宝等;《测绘工程》;1-4+11 * |
多尺度地理空间线状目标形状相似性的度量;江浩;褚衍东;闫浩文;郭丽峰;;测绘科学(第05期);35-38 * |
线要素化简及参数自动设置的案例推理方法;何海威;钱海忠;段佩祥;谢丽敏;罗登瀚;;武汉大学学报(信息科学版)(第03期);31-39 * |
Also Published As
Publication number | Publication date |
---|---|
CN114913263A (zh) | 2022-08-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103985155B (zh) | 基于映射法的散乱点云Delaunay三角剖分曲面重构方法 | |
CN101777178B (zh) | 一种图像修复方法 | |
CN110262250B (zh) | 一种基于粒子群算法的b样条曲线拟合方法及*** | |
WO2011111680A1 (ja) | 表面加工データの作成方法および装置 | |
CN110599506B (zh) | 一种复杂异形曲面机器人三维测量的点云分割方法 | |
Urick et al. | Watertight Boolean operations: A framework for creating CAD-compatible gap-free editable solid models | |
CN109712143B (zh) | 一种基于超像素多特征融合的快速图像分割方法 | |
CN108230452B (zh) | 一种基于纹理合成的模型补洞方法 | |
Lin et al. | Quality guaranteed all-hex mesh generation by a constrained volume iterative fitting algorithm | |
CN104504693B (zh) | 一种基于义齿简单冠修复网格模型的颈缘线提取方法 | |
Wolf | Scale independent surface characterisation: Geography meets precision surface metrology | |
CN100489713C (zh) | 直缝焊管排辊成型机组全流程自动化参数建模的方法 | |
Wang et al. | Multi-dimensional dynamic programming in ruled surface fitting | |
Liu et al. | An edge-sensitive simplification method for scanned point clouds | |
CN114913263B (zh) | 一个基于多尺度空间相似度的线状地物自动化简方法 | |
CN108122202A (zh) | 一种基于形状上下文的地图线状要素多尺度信息派生方法 | |
CN111079326A (zh) | 二维各向异性网格单元度量张量场光滑化方法 | |
CN101578632B (zh) | 模糊边缘平滑度先验及其在阿尔法通道超分辨率上的应用 | |
CN113987610A (zh) | 基于网格映射的不同分辨率服装曲面网格的匹配方法 | |
CN102982552A (zh) | 一种基于里奇流的表面配准方法 | |
Steger | Subpixel-precise extraction of watersheds | |
CN110532670B (zh) | 一种适分析样条实体模型构建方法及*** | |
CN112102178A (zh) | 一种点云保特征去噪方法、装置、电子设备及存储介质 | |
CN110070120B (zh) | 基于判别采样策略的深度度量学习方法及*** | |
CN114119628B (zh) | 一种基于特征模板的点云精确分割方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |