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 的输入，即除了 Polygons 或 MultiPolygons。如 Google KML docs 中所述, GeometryCollections 在 KML 中不受支持。

您可以使用 ST_Dump 两次调用 ST_Union 来解决此问题将合并的几何体分开，只从第一个联合体(转储)中选择那些产生多边形或多边形的几何体——你可以在不产生 GeometryCollection 的情况下组合它们。您需要使用两次 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)))

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

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