postgresql - 在 PostgreSQL/PostGIS 中为多个几何图形生成单个 KML

Question

我有一个 PostgreSQL 表，如下所示:

CREATE TABLE area
(
  area_code character varying(10) NOT NULL,
  shape geometry NOT NULL,
  CONSTRAINT pk_area PRIMARY KEY (area_code)
)

并且在几何列上有一个索引。

我想编写一个函数，传入一个area_code数组并返回一个表示这些区域的组合KML的值。要返回 KML，我使用 ST_AsKML 。要将所有区域合并在一起，我使用 ST_UNION 。为了确保联合返回单个几何图形而不是集合，我使用 ST_MULTI 。这一切给了我:

SELECT ST_AsKML(ST_MULTI(ST_UNION(shape))) as KML
FROM area
WHERE area_code = ANY(thearray)

在 pgAdmin III 中运行此命令时，我收到错误:

ERROR: lwgeom_to_kml2: 'GeometryCollection' geometry type not supported SQL state: XX000 Context: SQL function "st_askml" statement 1

但据我所知，这不应该发生。 ST_MULTI 应确保传递给 ST_AsKML 的值不是 GeometryCollection。我的方法是否正确但需要修复，还是我使用完全错误的技术来生成 KML？

Answer 1

如果您的 ST_Union 操作出于某种原因生成 Points 或 Linestrings，或者您在其中混合了几何类型，则可能会发生这种情况您对 ST_Union 的输入，即除多边形或多多边形之外的输入。正如 Google KML docs 中所述, KML 不支持 GeometryCollections。

您可以通过使用 ST_Dump 调用 ST_Union 两次来解决此问题将联合几何图形分开，并仅从第一个联合(转储)中选择那些生成多边形或多多边形的几何图形 - 您可以将它们组合在一起，而无需生成几何集合。您需要使用 ST_Union 两次，就像您在 Polygon 和 MultiPolygon 上使用 ST_Collect 一样，您将再次获得 GeometryCollection。

WITH polygons (geom) as 
   (SELECT (ST_Dump(ST_Union(shape))).geom 
       FROM area 
       WHERE area_code = any(thearray))
   SELECT ST_AsKML(ST_Multi(ST_Union(geom))) as KML 
       FROM polygons 
       WHERE ST_GeometryType(geom) in ('ST_Polygon','ST_MultiPolygon');

通过运行以下查询，您可以看到 ST_Multi 的问题，该查询尝试从点和线串中创建多几何图形:

SELECT 
   ST_AsText(
       ST_Multi(
         ST_Union(
           ST_MakePoint(0,0),
           ST_MakeLine(ST_MakePoint(5,5), ST_MakePoint(10,10))
         )
      )
 );

返回GEOMETRYCOLLECTION(POINT(0 0),LINESTRING(5 5,10 10))。 ST_Multi 不是一个聚合函数，因此本质上只是将 Multi 和一对额外的 () 添加到任何可以成为 Multi 的对象中，但如果您有两种类型的几何图形，它们的并集将强制进行几何集合，则不会有任何影响。

编辑:注意之间的区别，

 SELECT 
   ST_AsText(
     ST_Multi(
       ST_Collect(
         ST_Multi(ST_Expand(ST_MakePoint(0,0),2)),
         ST_Expand(ST_MakePoint(5,5), 1)
      )
   )
);

产生GEOMETRYCOLLECTION(MULTIPOLYGON(((-2 -2,-2 2,2 2,2 -2,-2 -2))),POLYGON((4 4,4 6,6 6, 6 4,4 4)))

和

SELECT 
   ST_AsText(
     ST_Multi(
       ST_Union(
         ST_Multi(ST_Expand(ST_MakePoint(0,0),2)),
         ST_Expand(ST_MakePoint(5,5), 1)
      )
   )
);

产生 MULTIPOLYGON(((-2 -2,-2 2,2 2,2 -2,-2 -2)),((4 4,4 6,6 6,6 4,4 4)))

因为多边形可以成为多多边形的一部分，所以在合并后，您可以避免第二种情况中可怕的 GeometryCollection。但是，如果您将任何线或点放入混合中，您将回到 GeomtryCollection 领域。

编辑2:根据OP的评论，Google KML docs声明 KML 不支持 GeometryCollections。由于无法在不生成几何集合的情况下将点和/或线与多边形组合，因此上述(多)多边形方法是唯一可行的方法。