PostGIS–线路合并方法比较

一 前言

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

路径合并前.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 行记录)

图形可视化如下:

待合并的线.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 行记录)
合并后的图形.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合并的时候,线路方向要求是很严格的,如下图:

同一个方向.png

ST_MakeLine结果.png

改变其中一条线的方向:

dianzhou=# update a set geom=ST_Reverse(geom) where gid=2;
UPDATE 1
改了其中一条线路方向.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,同向的效率快,但是不能做到对方向的纠正。

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

You may also like...

退出移动版