前言
众所周知，熟悉GIS桌面软件的同学一定都知道，想要对空间中的两个或者多个地理数据进行融合，可以使用融合的工具。常见的桌面软件，比如ArcGIS或者QGIS、SuperMap等均有成熟的工具。以QGIS为例，可以使用工具包中的合并工具，如下所示：
在日常的工作当中，可以采用这些工具进行数据的融合。只是这些数据操作通常需要人工进行干预，而且数据不能实时进行。比如有时候，我们已经将这些数据保存到了空间数据库PostGIS中，而用户可能会随机抽取一些空间数据进行动态的融合。
应对上面的场景，我们就可以采用SQL查询的方式来满足。比如在实际需求中，需要将某几个城市的空间数据进行融合形成新的数据以供查询分析。而在进行数据融合的时候，通常我们可以选用PostGIS中的ST_Union()和ST_Collect()这两个函数。
这两个函数乍一看几乎是一样的，其实两者还是有一些区别的，但是通过官网的介绍又不是特别清楚。因此本文采用实例的方式进行讲解，将ST_Union()和ST_Collect()这两个函数进行空间融合实战，我们通过实际的两个面的函数操作和实际场景融合，来看看两个函数的具体使用。
一、ST_Union()简介
在PostGIS中，对于st_union函数的定义有以下的几个入参，我们来看一下定义：
st_union函数的定义如下：
CREATE OR REPLACE FUNCTION "public"."st_union"("geom1" "public"."geometry", "geom2" "public"."geometry")
RETURNS "public"."geometry" AS '$libdir/postgis-3', 'ST_Union'
LANGUAGE c IMMUTABLE STRICT
COST 100001、方法说明
与Java的方法重载类似，st_union方法有5个重载。
|
序号 |
方法 |
|
1 |
st_union(geometry, geometry) |
|
2 |
st_union(geometry, geometry, float8) |
|
3 |
st_union(_geometry) |
|
4 |
st_union(geometry) |
|
5 |
st_union(geometry, float8) |
2、参数介绍
在上面的方法介绍中，geometry表示一个空间对象，float8的参数表示数据精度。注意，使用st_union函数进行运算之后，得到的结果是合并输入几何体，合并几何体以生成没有重叠的结果几何体。输出可以是原子几何图形、多重几何图形或几何图形集合。请一定注意这里合并的描述。
二、ST_Collect()简介
首先来看在PostGIS数据库中对于st_collect()函数的定义。
st_collect函数的定义SQL如下：
CREATE OR REPLACE FUNCTION "public"."st_collect"("geom1" "public"."geometry", "geom2" "public"."geometry")
RETURNS "public"."geometry" AS '$libdir/postgis-3', 'LWGEOM_collect'
LANGUAGE c IMMUTABLE
COST 501、方法说明
|
序号 |
函数 |
说明 |
|
1 |
geometry ST_Collect(geometry g1, geometry g2); |
接受两个输入几何 |
|
2 |
geometry ST_Collect(geometry[] g1_array); |
接受几何图形数组 |
|
3 |
geometry ST_Collect(geometry set g1field); |
接受几何图形行集的聚合函数 |
2、参数介绍
将几何图形收集（聚合）到几何图形集合中。结果是一个多*或一个GeometryCollection，这取决于输入几何的类型是相同的还是不同的(同构或异类)。输入几何图形在集合中保持不变。这里特别要注意的描述是聚合，与上面的函数相比，合并和聚合是两个概念。
3、两者区别
如果单纯的从文字上讲，合并和聚合是一个近义词，表示两者的概念是差不多的。实际上两者的区别比较大，尤其是在空间运算中。可以用这句话来概括：ST_Collect将几何图形聚合到一个集合中，而不会以任何方式更改它们。ST_UNION以几何形式合并重叠的几何图形，并在交点处拆分线串。融合边界时，它可能会返回单个几何图形。
简单来理解就是，使用st_collect进行运算，不会改变数据；而用st_union进行空间运算，会改变原有的数据。下面就进行实际的讲解。
三、实际案例实践
在这里将采用面数据进行讲解，这里选择的面数据又分两种。第一种情况是空间不重叠融合，第二种是空间重叠融合，通过实际的例子来看一下在不同的空间关系中，两个函数的融合结果是什么样的。
1、不重叠融合
空间不重叠融合中，我们采用保存在PostGIS数据库中的城市信息进行实例开展，对长株潭三个城市进行融合。
biz_city的表数据如下所示：
st_union融合sql：
select st_asgeojson(st_union(geom)) from biz_city where city_name in ('长沙市','株洲市','湘潭市');查询结果如下：
来看QGIS的结果：
通过上图可以看出，使用st_union之后，三个城市已经完全融合，从而形成了一个全新的数据。下面再来看看st_collect函数的运行效果。
select st_asgeojson(st_collect(geom)) from biz_city where city_name in ('长沙市','株洲市','湘潭市');
来看一下在QGIS中的展示结果：
从总体来看，虽然长株潭三个城市同样合成了，但是三个城市的边界还是存在的，并没有进行融合。而在geojson的type中明确了两个的区别。
结论：在空间不重叠的情况下，使用st_union和st_collect的合并结果虽然不大，但是对于数据是否改变来说，影响还是很大的。
2、空间重叠融合
首先我们准备两份空间上重叠的数据，我们可以采用一下的模拟数据进行展示。
A B 图作为测试图形，用来验证这些函数的使用
-- A 图
select ST_GeomFromGeoJSON('{"type":"Polygon","coordinates":[[[120.52757263183594,31.495432803134843],[120.59211730957031,31.31199502365151],[120.89012145996094,31.35950051982242],[120.74729919466666,31.483245492650792],[120.52757263183594,31.495432803134843]]]}')
-- B 图
select ST_GeomFromGeoJSON('{"type":"Polygon","coordinates":[[[120.74386596679688,31.487235582017444],[120.84548950195312,31.269747790889888],[121.01783752441405,31.378261512889125],[120.91827392578125,31.487821121636433],[120.74386596679688,31.487235582017444]]]}')将上述两个空间数据转成geojson，在QGIS中展示如下，首先是A图
下面是B图：
将两者空间叠加后，可以看到空间关系：
st_union融合语句：
-- ST_Union在几何形状重叠的地方对其进行几何合并，并在相交处拆分线串。溶解边界时，它可能会返回单个几何。
select ST_Union (
ST_GeomFromGeoJSON('{"type":"Polygon","coordinates":[[[120.52757263183594,31.495432803134843],[120.59211730957031,31.31199502365151],[120.89012145996094,31.35950051982242],[120.52757263183594,31.495432803134843]]]}'),
ST_GeomFromGeoJSON('{"type":"Polygon","coordinates":[[[120.74386596679688,31.487235582017444],[120.84548950195312,31.269747790889888],[121.01783752441405,31.378261512889125],[120.91827392578125,31.487821121636433],[120.74386596679688,31.487235582017444]]]}')
)
st_collect融合语句：
-- ST_Collect将几何形状聚合到一个集合中，而无需进行任何更改。
select ST_Collect (
ST_GeomFromGeoJSON('{"type":"Polygon","coordinates":[[[120.52757263183594,31.495432803134843],[120.59211730957031,31.31199502365151],[120.89012145996094,31.35950051982242],[120.74729919466666,31.483245492650792],[120.52757263183594,31.495432803134843]]]}'),
ST_GeomFromGeoJSON('{"type":"Polygon","coordinates":[[[120.74386596679688,31.487235582017444],[120.84548950195312,31.269747790889888],[121.01783752441405,31.378261512889125],[120.91827392578125,31.487821121636433],[120.74386596679688,31.487235582017444]]]}')
)
通过实验，同样可以得到上面的结论，使用st_collect对于原来的空间数据没有什么改变，而st_untion对数据的改变很大。
所以，到了这里一定要深深的理解合并和聚合的区别。在实际项目过程中，根据需要进行合理的选择。
总结
以上就是本文的主要内容，本文采用实例的方式进行讲解，将ST_Union()和ST_Collect()这两个函数进行空间融合实战，我们通过实际的两个面的函数操作和实际场景融合，希望通过实战的方式让你了解两者的区别，在具体使用空间函数的时候给予充分的评估和参考。
