CN103824510A - 一种基于Voronoi图的筛选电子地图点状要素的方法 - Google Patents
一种基于Voronoi图的筛选电子地图点状要素的方法 Download PDFInfo
- Publication number
- CN103824510A CN103824510A CN201410049939.7A CN201410049939A CN103824510A CN 103824510 A CN103824510 A CN 103824510A CN 201410049939 A CN201410049939 A CN 201410049939A CN 103824510 A CN103824510 A CN 103824510A
- Authority
- CN
- China
- Prior art keywords
- grid
- point
- poi
- voronoi
- display level
- 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.)
- Pending
Links
Images
Landscapes
- Instructional Devices (AREA)
Abstract
本发明提供一种基于Voronoi图的筛选电子地图点状要素的方法,对POI点集进行网格划分;若网格中只有一个POI点,则直接转到步骤30;若网格中不存在POI点,则不对该网格进行任何操作;若网格中含有一个以上的POI点,则在该网格范围内,求取Delaunay三角网,根据所求得的Delaunay三角网求取该网格中所包含每个POI点所对应的Voronoi图;根据地图具体应用需求,确定每个POI点的最终显示等级;根据所要确定保留的POI点与周边网格中已确定的POI点的距离,最终确定每个网格所保留的POI点,本发明方法筛选保留的POI能够更好地体现其在空间分布的权重,并且提高了地图数据编译效率。
Description
技术领域
本发明涉及一种基于Voronoi图的筛选电子地图点状要素的方法。
背景技术
电子地图在不同的比例尺下,显示的POI稀疏程度不同,因此对地图厂家提供的POI数据源,在显示前有一个筛选处理过程。目前对POI的筛选多是根据地图用途,依据POI的类型划分等级,然后按等级进行筛选,没有考虑到POI的空间分布特性、相邻POI在显示时可能存在的标注相互冲突问题。Voronoi图,又叫泰森多边形或狄利克雷Dirichlet图,它是由一组由连接两邻点直线的垂直平分线组成的连续多边形组成。POI为导航地图信息,每个POI包含名称、类别、经度、纬度、附近的酒店、饭店,商铺等信息,导航地图数据是整个导航产业的基石。
发明内容
本发明要解决的技术问题,在于提供一种基于Voronoi图的筛选电子地图点状要素的方法,提高了地图数据编译效率,解决POI对应的文字标注在地图上显示的相互冲突问题。
本发明具体包括如下步骤:
步骤10、对POI点集进行网格划分;
步骤20、若网格中只有一个POI点,则直接转到步骤30;若网格中不存在POI点,则不对该网格进行任何操作;若网格中含有一个以上的POI点,则在该网格范围内,求取Delaunay三角网,根据所求得的Delaunay三角网求取该网格中所包含每个POI点所对应的Voronoi图,所求的POI点所对应的Voronoi图用于确定POI点的显示等级;
步骤30、根据地图具体应用需求,确定每个POI点的最终显示等级;
步骤40、根据所要确定保留的POI点与周边网格中已确定的POI点的距离,最终确定每个网格所保留的POI点。
进一步地,所述步骤10进一步具体为:
将地图数据看作是由一系列的方形网格组成,沿地图的经度、纬度方向按设定的间距进行划分,网格编号从小到大按顺序存放。
进一步地,所述步骤20进一步具体为:
若网格中只有一个POI点,则直接转到步骤30;若网格中不存在POI点,则不对该网格进行任何操作;若网格中含有一个以上的POI点,则在该网格范围内,采用逐点***法先计算Delaunay三角网,再根据Delaunay三角网和Voronoi图的对偶性质,求取该网格中所包含POI点集中每个POI点对应的Voronoi图。
进一步地,所述步骤30进一步具体包括如下步骤:
步骤31、根据地图具体应用需求,由POI 的类型号划分每个POI的显示等级;
步骤32、计算每个POI点对应的Voronoi多边形面积;(注:若某些点的Voronoi图和网格边界线相交,此时要把网格边界线作为该点Voronoi多边形的一部分,从而形成一个封闭的凸多边形。)
步骤33、对于显示等级相同的POI点,根据POI点对应的Voronoi多边形面积大小进行确认每个POI点的最终显示等级。
进一步地,所述步骤40进一步具体为:
把网格中显示等级最大中的Voronoi多边形面积最大的POI点作为预保留点,并判断该点与周边网格的预保留点的距离是否满足距离限制要求,若满足要求,则该点作为该网格确定的保留点;若不满足则选取显示等级次之的点作为预保留点,重新进行限制距离判定,直至找到满足距离限制要求的保留点。
本发明具有如下优点:
对POI数据作网格化划分预处理,以网格为单位计算各点的Voronoi图,大大减少了计算Voronoi图的时间,提高了地图数据编译效率,在筛选POI时加入了距离限制条件,避免了相邻POI名称标注间的相互冲突,在地图数据编译阶段就加入此限制条件筛选POI,在导航终端显示时不需要再采用复杂的避让算法来解决POI名称标注冲突问题,直接显示保留POI即可,较大提升了导航终端地图显示的性能。
附图说明
下面参照附图结合实施例对本发明作进一步的说明。
图1为本发明方法执行流程图。
图2为本发明方法具体实施例的网格划分示意图。
图3为本发明方法具体实施例的单个网格生成Voronoi示意图。
图4为本发明方法具体实施例的示意图。
具体实施方式
如图1所示,一种基于Voronoi图的筛选电子地图点状要素的方法,包括如下步骤:
步骤10、对POI点集进行网格划分;
步骤20、若网格中只有一个POI点,则直接转到步骤30;若网格中不存在POI点,则不对该网格进行任何操作;若网格中含有一个以上的POI点,则在该网格范围内,求取Delaunay三角网,根据所求得的Delaunay三角网求取该网格中所包含每个POI点所对应的Voronoi图,所求的POI点所对应的Voronoi图用于确定POI点的显示等级;
步骤30、根据地图具体应用需求,确定每个POI点的最终显示等级;
步骤40、根据所要确定保留的POI点与周边网格中已确定的POI点的距离,最终确定每个网格所保留的POI点。
所述步骤10进一步具体为:
将地图数据看作是由一系列的方形网格组成,沿地图的经度、纬度方向按设定的间距进行划分,网格编号从小到大按顺序存放,且不同比例下的地图网格大小不同。
所述步骤20进一步具体为:
若网格中只有一个POI点,则直接转到步骤30;若网格中不存在POI点,则不对该网格进行任何操作;若网格中含有一个以上的POI点,则在该网格范围内,采用逐点***法先计算Delaunay三角网,再根据Delaunay三角网和Voronoi图的对偶性质,求取该网格中所包含POI点集中每个POI点对应的Voronoi图。
所述步骤30进一步具体包括如下步骤:
步骤31、根据地图具体应用需求,由POI 的类型号划分每个POI的显示等级;
步骤32、计算每个POI点对应的Voronoi多边形面积;
步骤33、对于显示等级相同的POI点,根据POI点对应的Voronoi多边形面积大小进行确认每个POI点的最终显示等级。
所述步骤40进一步具体为:
把网格中显示等级最大中的Voronoi多边形面积最大的POI点作为预保留点,并判断该点与周边网格的预保留点的距离是否满足距离限制要求,若满足要求,则该点作为该网格确定的保留点;若不满足则选取显示等级次之的点作为预保留点,重新进行限制距离判定,直至找到满足距离限制要求的保留点。
具体实施例如下:
步骤一、POI点集网格化划分,地图在不同比例尺下的网格尺寸大小不同,加快计算Voronoi图的速度,解决计算海量POI点数据Voronoi耗时长问题。把全国地图数据沿经度、纬度方向按一定间距进行划分,把全国范围内的地图数据看作是由一系列的方形网格组成,网格编号从小到大按顺序存放,具体划分如图2所示;
步骤二、在各个网格范围内,若一个网格中只有一个POI点,则直接转到步骤30;否则,在该网格范围内,采用逐点***法先计算Delaunay三角网,再根据Delaunay三角网和Voronoi图的对偶性质,求取该网格中所包含POI点集中每个POI点对应的Voronoi图,完成各个比例下POI数据的网格化划分后,以网格为单位,计算每个网格中所包含POI点集中每个点的Voronoi图,如图3所示为单个网格计算后生成Voronoi图后的实例;
步骤三、根据地图具体应用需求,由POI的类型号,划分各个POI的第1级显示等级,计算各个POI点对应的Voronoi多边形面积,某些点的Voronoi多边形和网格边界线相交,此时要把网格边界线作为该点Voronoi多边形的一部分,从而形成一个封闭的凸多边形,其面积采用下式计算:
其中为Voronoi多边形上相邻点的x、y坐标,对第1级显示等级相同的POI点对应的Voronoi多边形面积由小到大排序,从而确定每个网格中各个POI点的最终显示等级,设某个网格中有n个POI点第1级显示等级相同,设表示i点对应的Voronoi多边形面积,表示i点对应的最终显示等级,则当时,有;
步骤四、把网格中显示等级最大中的Voronoi多边形面积最大的POI点作为预保留点,并判断该点与周边网格的预保留点的距离是否满足距离限制要求,若满足要求,则该点作为该网格确定的保留点;若不满足则选取最终显示等级次之的点作为预保留点,重新进行限制距离判定,直至找到满足距离限制要求的保留点;
如图4示例所示:对编号为491,909,751和编号为491,936,123这两个相邻网格,根据步骤二和步骤三计算出各点对应的Voronoi多边形面积可以确定,“御香苑”和“天歧祥市场”分别是这两个网格中显示等级最高的点,计算这两个点距离,如符合距离限制条件,则把这两个点作为这两个网格的最终保留点;如这两个点不满足距离限制条件,则选取网491,936,123中最终显示等级次之的“停车场”作为预保留点,判断“御香苑”和“停车场”是否满足距离限制条件,若满足则把“停车场”作为编号为491,936,123的网格的最终保留点。
本发明具有如下优点:
对POI数据作网格化划分预处理,以网格为单位计算各点的Voronoi图,大大减少了计算Voronoi图的时间,提高了地图数据编译效率,在筛选POI时加入了距离限制条件,避免了相邻POI名称标注间的相互冲突,在地图数据编译阶段就加入此限制条件筛选POI,在导航终端显示时不需要再采用复杂的避让算法来解决POI名称标注冲突问题,直接显示保留POI即可,较大提升了导航终端地图显示的性能。
虽然以上描述了本发明的具体实施方式,但是熟悉本技术领域的技术人员应当理解,我们所描述的具体的实施例只是说明性的,而不是用于对本发明的范围的限定,熟悉本领域的技术人员在依照本发明的精神所作的等效的修饰以及变化,都应当涵盖在本发明的权利要求所保护的范围内。
Claims (5)
1.一种基于Voronoi图的筛选电子地图点状要素的方法,其特征在于:包括如下步骤:
步骤10、对POI点集进行网格划分;
步骤20、若网格中只有一个POI点,则直接转到步骤30;若网格中不存在POI点,则不对该网格进行任何操作;若网格中含有一个以上的POI点,则在该网格范围内,求取Delaunay三角网,根据所求得的Delaunay三角网求取该网格中所包含每个POI点所对应的Voronoi图;
步骤30、根据地图具体应用需求,确定每个POI点的最终显示等级;
步骤40、根据所要确定保留的POI点与周边网格中已确定的POI点的距离,最终确定每个网格所保留的POI点。
2.根据权利要求1所述的一种基于Voronoi图的筛选电子地图点状要素的方法,其特征在于:所述步骤10进一步具体为:
将地图数据看作是由一系列的方形网格组成,沿地图的经度、纬度方向按设定的间距进行划分,网格编号从小到大按顺序存放。
3.根据权利要求1所述的一种基于Voronoi图的筛选电子地图点状要素的方法,其特征在于:所述步骤20进一步具体为:
若网格中只有一个POI点,则直接转到步骤30;若网格中不存在POI点,则不对该网格进行任何操作;若网格中含有一个以上的POI点,则在该网格范围内,采用逐点***法先计算Delaunay三角网,再根据Delaunay三角网和Voronoi图的对偶性质,求取该网格中所包含POI点集中每个POI点对应的Voronoi图。
4.根据权利要求1所述的一种基于Voronoi图的筛选电子地图点状要素的方法,其特征在于:所述步骤30进一步具体包括如下步骤:
步骤31、根据地图具体应用需求,由POI 的类型号划分每个POI的显示等级;
步骤32、计算每个POI点对应的Voronoi多边形面积;(注:若某些点的Voronoi图和网格边界线相交,此时要把网格边界线作为该点Voronoi多边形的一部分,从而形成一个封闭的凸多边形);
步骤33、对于显示等级相同的POI点,根据POI点对应的Voronoi多边形面积大小进行确认每个POI点的最终显示等级。
5.根据权利要求1所述的一种基于Voronoi图的筛选电子地图点状要素的方法,其特征在于:所述步骤40进一步具体为:
把网格中显示等级最大中的Voronoi多边形面积最大的POI点作为预保留点,并判断该点与周边网格的预保留点的距离是否满足距离限制要求,若满足要求,则该点作为该网格确定的保留点;若不满足则选取显示等级次之的点作为预保留点,重新进行限制距离判定,直至找到满足距离限制要求的保留点。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410049939.7A CN103824510A (zh) | 2013-12-31 | 2014-02-13 | 一种基于Voronoi图的筛选电子地图点状要素的方法 |
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310754876 | 2013-12-31 | ||
CN201310754876.0 | 2013-12-31 | ||
CN201410049939.7A CN103824510A (zh) | 2013-12-31 | 2014-02-13 | 一种基于Voronoi图的筛选电子地图点状要素的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103824510A true CN103824510A (zh) | 2014-05-28 |
Family
ID=50759542
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410049939.7A Pending CN103824510A (zh) | 2013-12-31 | 2014-02-13 | 一种基于Voronoi图的筛选电子地图点状要素的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103824510A (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104268210A (zh) * | 2014-09-12 | 2015-01-07 | 东北大学 | 基于局部超集的cpir-v最近邻隐私保护查询方法 |
CN104933437A (zh) * | 2015-05-19 | 2015-09-23 | 广州杰赛科技股份有限公司 | 动态网格的生成方法、生成装置、聚合方法和聚合装置 |
CN105046267A (zh) * | 2015-05-19 | 2015-11-11 | 广州杰赛科技股份有限公司 | 动态网格的生成方法、生成装置、聚合方法和聚合装置 |
CN105701204A (zh) * | 2016-01-12 | 2016-06-22 | 中国测绘科学研究院 | 基于道路网的电子地图兴趣点的提取方法及显示方法 |
CN106156472A (zh) * | 2015-04-20 | 2016-11-23 | 中国科学院遥感与数字地球研究所 | 一种遥感卫星对地覆盖能力分析方法 |
CN106441269A (zh) * | 2016-08-31 | 2017-02-22 | 武汉拓普伟域网络有限公司 | 一种用于景区电子语音导游的区间定位方法 |
CN107480175A (zh) * | 2017-06-30 | 2017-12-15 | 广州数鹏通科技有限公司 | 一种基于网格化的地理信息数据筛选方法、电子设备、存储介质及*** |
CN113157836A (zh) * | 2021-03-12 | 2021-07-23 | 广西壮族自治区中国科学院广西植物研究所 | 一种空间点数据集周围数据提取方法、装置及存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101655375A (zh) * | 2009-07-28 | 2010-02-24 | 北京四维图新科技股份有限公司 | 导航电子地图文字注记显示等级的自动划分方法及装置 |
CN102073981A (zh) * | 2011-02-23 | 2011-05-25 | 中国测绘科学研究院 | 一种关联要素限制下的点群地理实体选取方法 |
CN102737542A (zh) * | 2012-02-10 | 2012-10-17 | 武汉大学 | 一种顾及多重约束条件的水深注记自动选取方法 |
CN102831645A (zh) * | 2012-07-18 | 2012-12-19 | 哈尔滨工程大学 | 一种应用于海底地形的数字高程模型的建立方法 |
-
2014
- 2014-02-13 CN CN201410049939.7A patent/CN103824510A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101655375A (zh) * | 2009-07-28 | 2010-02-24 | 北京四维图新科技股份有限公司 | 导航电子地图文字注记显示等级的自动划分方法及装置 |
CN102073981A (zh) * | 2011-02-23 | 2011-05-25 | 中国测绘科学研究院 | 一种关联要素限制下的点群地理实体选取方法 |
CN102737542A (zh) * | 2012-02-10 | 2012-10-17 | 武汉大学 | 一种顾及多重约束条件的水深注记自动选取方法 |
CN102831645A (zh) * | 2012-07-18 | 2012-12-19 | 哈尔滨工程大学 | 一种应用于海底地形的数字高程模型的建立方法 |
Non-Patent Citations (2)
Title |
---|
HAOWEN YAN,ROBERT WEIBEL: "《An algorithm for point cluster generalization based on the Voronoi diagram》", 《COMPUTERS & GEOSCIENCES》 * |
于艳平等: "《点群选取与化简算法时间复杂度研究》", 《南京师大学报(自然科学版)》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104268210A (zh) * | 2014-09-12 | 2015-01-07 | 东北大学 | 基于局部超集的cpir-v最近邻隐私保护查询方法 |
CN104268210B (zh) * | 2014-09-12 | 2017-09-26 | 东北大学 | 基于局部超集的cpir‑v最近邻隐私保护查询方法 |
CN106156472A (zh) * | 2015-04-20 | 2016-11-23 | 中国科学院遥感与数字地球研究所 | 一种遥感卫星对地覆盖能力分析方法 |
CN105046267B (zh) * | 2015-05-19 | 2018-08-28 | 广州杰赛科技股份有限公司 | 动态网格的生成方法、生成装置、聚合方法和聚合装置 |
CN104933437A (zh) * | 2015-05-19 | 2015-09-23 | 广州杰赛科技股份有限公司 | 动态网格的生成方法、生成装置、聚合方法和聚合装置 |
CN105046267A (zh) * | 2015-05-19 | 2015-11-11 | 广州杰赛科技股份有限公司 | 动态网格的生成方法、生成装置、聚合方法和聚合装置 |
CN104933437B (zh) * | 2015-05-19 | 2019-10-25 | 广州杰赛科技股份有限公司 | 动态网格的生成方法、生成装置、聚合方法和聚合装置 |
CN105701204A (zh) * | 2016-01-12 | 2016-06-22 | 中国测绘科学研究院 | 基于道路网的电子地图兴趣点的提取方法及显示方法 |
CN105701204B (zh) * | 2016-01-12 | 2018-11-13 | 中国测绘科学研究院 | 基于道路网的电子地图兴趣点的提取方法及显示方法 |
CN106441269A (zh) * | 2016-08-31 | 2017-02-22 | 武汉拓普伟域网络有限公司 | 一种用于景区电子语音导游的区间定位方法 |
CN106441269B (zh) * | 2016-08-31 | 2019-10-29 | 武汉拓普伟域网络有限公司 | 一种用于景区电子语音导游的区间定位方法 |
CN107480175A (zh) * | 2017-06-30 | 2017-12-15 | 广州数鹏通科技有限公司 | 一种基于网格化的地理信息数据筛选方法、电子设备、存储介质及*** |
CN113157836A (zh) * | 2021-03-12 | 2021-07-23 | 广西壮族自治区中国科学院广西植物研究所 | 一种空间点数据集周围数据提取方法、装置及存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103824510A (zh) | 一种基于Voronoi图的筛选电子地图点状要素的方法 | |
CN102142023B (zh) | 地图数据、地图数据生成方法、存储介质以及导航设备 | |
CN102509354B (zh) | 一种与影像同步变化的投影数字高程模型制作方法 | |
CN103562681A (zh) | 产生用于导航装置的数据库的方法、输出三维地图的方法和导航装置 | |
CN103996089A (zh) | 基于gis的输电线路最优路径生成方法 | |
CN103529486B (zh) | 一种地球化学异常圈定方法 | |
CN103995281A (zh) | 一种三维观测***快速滚动布设的方法 | |
CN105844224A (zh) | 一种车载LiDAR 道路点云快速有序化方法 | |
CN113989451A (zh) | 高精地图构建方法、装置及电子设备 | |
CN108507582A (zh) | 一种导航路线的推送方法及装置 | |
CN112233205B (zh) | 基于离散数据进行分区和裁剪的电子地图制图方法及*** | |
CN105403913A (zh) | 叠前深度偏移方法和装置 | |
CN106023311B (zh) | 提高三维地形生成精度的方法 | |
CN105373538B (zh) | 一种多边形对象的显示层级确定方法及装置 | |
US20110119311A1 (en) | Map data creating device and map drawing device | |
CN105426387A (zh) | 一种基于K-means算法的地图聚合方法 | |
CN104866682A (zh) | 一种基于地面形变进行页岩气勘探区构造应力场反演的方法 | |
CN108121853A (zh) | 一种基于AutoCAD计算挖填方工程量的***与方法 | |
CN103021257A (zh) | 电子地图的生成方法及装置 | |
CN114640619B (zh) | 基于平均通信距离的天基网络拓扑设计方法及相关设备 | |
CN115830255A (zh) | 一种仿真场景生成方法、装置、电子设备和存储介质 | |
CN102496330A (zh) | 六角网格同构模型及其建模方法和应用 | |
US20190137283A1 (en) | Updating missing attributes in navigational map data via polyline geometry matching | |
JP7239053B2 (ja) | 生成装置、生成方法及び生成プログラム | |
CN104931015A (zh) | 基于dtm数据模型的横断面自动成图***与方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20140528 |
|
WD01 | Invention patent application deemed withdrawn after publication |