PostGIS--线路合并方法比较

简介: 一 前言  工作中处理GIS数据时,最常规的操作是对多条线记录,他们的几何关系存在首尾相连的情况,需要合并成一个完整的线路。如路径分析中,路网都是一段一段的,显示的时候很不好看,合并成一个完整的显示比较好,如下图:路径合并前.png路径合并后.png  PostGIS中,线路合并有几个相似的方法,如ST_Union,ST_LineMerge,ST_MakeLine方法,几种方法看起来很相似,实际应用的差别挺大。

一 前言

  工作中处理GIS数据时,最常规的操作是对多条线记录,他们的几何关系存在首尾相连的情况,需要合并成一个完整的线路。如路径分析中,路网都是一段一段的,显示的时候很不好看,合并成一个完整的显示比较好,如下图:


img_c8eaf742951927dc0a20ceddc2474d7a.png
路径合并前.png

img_3262b8048a54039830da70bfe49dc862.png
路径合并后.png

  PostGIS中,线路合并有几个相似的方法,如ST_Union,ST_LineMerge,ST_MakeLine方法,几种方法看起来很相似,实际应用的差别挺大。本文主要阐述这几种方法的异同点和使用场景。

二 ST_Union

  ST_Union是图形批处理中的一个函数,用于将图形联合(不是合并)成一条图形记录的方法。

#方法简介
geometry ST_Union(geometry set g1field);
geometry ST_Union(geometry g1, geometry g2);
geometry ST_Union(geometry[] g1_array);

输入参数:一般是图形的记录集或者图形的几何数组,也可以指定两个图形建立联合。
输出参数:返回的一般是Multi的geometry,或者geometry的collection。
使用场景:一般和聚合函数搭配使用。
使用示例:将数据中拥有相同start_id,end_id的线数据聚合成一条线路:
先查看下测试数据:

dianzhou=# select start_vid,end_vid,st_astext(geom) from a;
 start_vid | end_vid |                                        st_astext
-----------+---------+-----------------------------------------------------------------------------------------
     14178 |   43866 | LINESTRING(110.3359653 20.03026554,110.3365341 20.030855085)
     14178 |   43866 | LINESTRING(110.3365341 20.030855085,110.33657406 20.03089509,110.33722278 20.031564555)
     14178 |   43866 | LINESTRING(110.33722278 20.031564555,110.33745228 20.03179437)
(3 行记录)

图形可视化如下:


img_08855b1ff6c81794e1e42ccf4ba4bdd2.png
待合并的线.png

执行ST_Union操作:

dianzhou=# select start_vid,end_vid,st_astext(ST_Union(geom)) from a group by start_vid,end_vid;
 start_vid | end_vid |                                                                                               st_astext                                                               
-----------+---------+--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
     14178 |   43866 | MULTILINESTRING((110.33722278 20.031564555,110.33745228 20.03179437),(110.3359653 20.03026554,110.3365341 20.030855085),(110.3365341 20.030855085,110.33657406 20.03089509,110.33722278 20.031564555))
(1 行记录)

  结果说明:拥有相同start_id,end_id的线数据聚合成一条线路,但是几何类型由LineString变成了MultiLineString。ST_Union函数只是将分离的图形不做太多的处理直接塞到了一起而已。

三 ST_LineMerge

  上图可知,几条单独的线都是首尾相连的,按道理他们一条条连接起来形成完整的线是可能的,所以很多用户希望不仅仅是图形联合了,数据类型最好也是LineString(Multi类型的图形比较复杂在分析时不好处理),那么ST_LineMerge方法就可以起到用武之地了。

#方法简介
geometry ST_LineMerge(geometry amultilinestring);

  该方法可以将MultiLineString中能合并的图形“缝补”到一起来,但是有的MultiLineString图形本来就是分离的,那就“缝补”不了了。

dianzhou=# select start_vid,end_vid,st_astext(ST_LineMerge(ST_Union(geom))) from a group by start_vid,end_vid;
 start_vid | end_vid |                                                                st_astext
-----------+---------+------------------------------------------------------------------------------------------------------------------------------------------
     14178 |   43866 | LINESTRING(110.3359653 20.03026554,110.3365341 20.030855085,110.33657406 20.03089509,110.33722278 20.031564555,110.33745228 20.03179437)
(1 行记录)
img_1a128e07a11ab7d12b09de53a03c5616.png
合并后的图形.png

如果是“缝补”不了的图形,如下:

SELECT ST_AsText(ST_LineMerge(
ST_GeomFromText('MULTILINESTRING((-29 -27,-30 -29.7,-36 -31,-45 -33),(-45.2 -33.2,-46 -32)) '))
);
st_astext
----------------
MULTILINESTRING((-45.2 -33.2,-46 -32),(-29 -27,-30 -29.7,-36 -31,-45 -33))

缝补不了的话还是返回MULTILINESTRING类型。

四 ST_MakeLine

  ST_Union方法可以结合聚合把图形联合到一起,ST_LineMerge可以将联合的MultiLineString中能合并到一起的图形合并缝补起来。ST_MakeLine方法是图形构造函数中的一种,他可以直接从 点(Point),多义点(MultiPoint),线(LineString)这些图形集合里直接构造一条新的线。

#用法
geometry ST_MakeLine(geometry set geoms);
geometry ST_MakeLine(geometry geom1, geometry geom2);
geometry ST_MakeLine(geometry[] geoms_array);

细心的读者会发现和ST_Union方法的使用很相同,但是该方法只能接受点线,而ST_Union使用是没这种限制的,因为ST_Union方法并不仅仅是处理线的。

#案例说明
dianzhou=# select start_vid,end_vid,ST_AsText(ST_MakeLine(geom)) geom  from a group by start_vid,end_vid;
 start_vid | end_vid |                                                                   geom
-----------+---------+------------------------------------------------------------------------------------------------------------------------------------------
     14178 |   43866 | LINESTRING(110.3359653 20.03026554,110.3365341 20.030855085,110.33657406 20.03089509,110.33722278 20.031564555,110.33745228 20.03179437)
(1 行记录)

注意事项:
ST_MakeLine合并的时候,线路方向要求是很严格的,如下图:


img_6f376b10bb6f6721e5bedf41e13ffd81.png
同一个方向.png

img_56b474c353c22a0c0507a32118da43b9.png
ST_MakeLine结果.png

改变其中一条线的方向:

dianzhou=# update a set geom=ST_Reverse(geom) where gid=2;
UPDATE 1
img_3c1627d7b6c2b8f69a31059ebb0d8666.png
改了其中一条线路方向.png

img_5ff40e087d374ac915872f098277160c.png
ST_MakeLine结果.png

我们看到了一条乱线,通过示意图可以得知,ST_MakeLine的工作原理是先合并同方向的,然后合并反方向的。这并不是我们期望的结果。。。
我后来改成这种方法:

dianzhou=# select st_linemerge(st_union(geom)) geom from a  group by start_vid,end_vid;

这样倒是可以了。。。

五 总结

方法名 入参 出参 使用场景
ST_Union 任何图形 Multi类型,collection类型 聚合后的图形联合
ST_LineMerge MultiLineString LineString,MultiLineString 将MultiLineString中相连的线路合并成完整的LineString
ST_MakeLine Point,MultiPoint,LineString LineString 构造新的线路

实际应用中,希望根据某些条件对相连的LineString聚合成一个新的LineString时,有两种方式:
1 ST_Union先聚合成MultiString再ST_LineMerge成LineString,效率差,但是能做到自动纠正方向。
2 ST_MakeLine直接对聚合的LineString构建成一个新的LineString,同向的效率快,但是不能做到对方向的纠正。

针对线路合并这种场景,按实际情况选择方案。。。

相关文章
|
存储 缓存 JSON
实战干货 | 分布式多级缓存设计方案
分布式多级缓存设计方案,解决海量数据读取的性能问题,包含多级缓存的存储设计,流程设计;利用多数据副本保证数据的可用性,同时通过不同数据源特点提供更高性能、更多场景数据差异化的支持
1820 0
实战干货 | 分布式多级缓存设计方案
|
数据可视化 前端开发 关系型数据库
基于Mybatis-Plus实现Geometry字段在PostGis空间数据库中的使用
本文讲解在mybatis-plus中操作geometry空间字段,同时实现查询和插入操作​。通过geojson,结合前端可视化组件即可完成​矢量数据的空间可视化。
3094 0
基于Mybatis-Plus实现Geometry字段在PostGis空间数据库中的使用
|
运维 负载均衡 网络协议
linux网络管理(链路聚合、桥接网络、故障排查、常用工具)
网卡的链路聚合就是将多块网卡连接起来,当一块网卡损坏,网络依旧可以正常运行,可以有效的防止因为网卡损坏带来的损失,同时也可以提高网络访问速度。
1531 0
linux网络管理(链路聚合、桥接网络、故障排查、常用工具)
|
10月前
|
小程序 测试技术 API
鸿蒙原生开发手记:03-元服务开发全流程(开发元服务,只需要看这一篇文章)
本文详细介绍元服务的开发及上架全流程,涵盖元服务的特点、创建项目、服务卡片、签名打包、开发测试及上架审核等环节,帮助开发者轻松掌握从零开始开发并发布元服务的全过程。元服务以其轻量、免安装、易于使用等特点,成为未来服务提供的重要形式。
742 13
鸿蒙原生开发手记:03-元服务开发全流程(开发元服务,只需要看这一篇文章)
|
11月前
|
机器学习/深度学习 人工智能 小程序
【一步步开发AI运动小程序】三、运动识别处理流程
随着人工智能技术的发展,阿里体育等公司推出的“乐动力”、“天天跳绳”等AI运动APP备受关注。本文将引导您从零开始开发一个AI运动小程序,使用“云智AI运动识别小程序插件”。文章介绍了视频帧、帧率FPS、抽帧和人体识别等基本概念,并详细说明了处理流程,包括抽帧、人体识别检测、骨骼图绘制和运动分析等步骤。下篇将介绍如何在小程序中实现抽帧。
|
数据采集 监控 安全
园区 AR 导航系统:技术创新与功能拓展
**维小帮园区AR导航系统**借助高精地图与AR技术,打破传统二维导航局限,提供三维立体地图和沉浸式导航体验。系统包括AR实景指引、VR全景预览、便捷预约参观和智能停车场导航等功能,大幅提升访客定位与通行效率,推动园区管理智能化。未来,该技术有望在更多领域广泛应用。
351 0
园区 AR 导航系统:技术创新与功能拓展
|
Java Apache Spring
若依配Mapper,若依修改Caused by: org.springframework.beans.factory.BeanCreationException: Error creating
若依配Mapper,若依修改Caused by: org.springframework.beans.factory.BeanCreationException: Error creating
|
人工智能 数据库 开发者
Python中的atexit模块:优雅地处理程序退出
Python中的atexit模块:优雅地处理程序退出
274 3
|
Java 测试技术 Spring
c3p0存在严重bug “APPARENT DEADLOCK“的问题
c3p0存在严重bug “APPARENT DEADLOCK“的问题
514 0
|
存储 NoSQL
TableStore: 海量结构化数据分层存储方案
### 前言 表格存储是阿里云自研分布式存储系统,可以用来存储海量结构化、半结构化的数据。表格存储支持高性能和容量型两种实例类型。高性能使用SSD的存储介质,针对读多写多的场景都有较好的访问延时。容量型使用的是SSD和SATA混合的存储介质。
9670 0

热门文章

最新文章