在PostgreSQL / PostGIS中为多个几何生成单个KML

时间:2014-10-29 12:53:15

标签: postgresql postgis

我有一个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中运行时,我收到错误:

  

错误:lwgeom_to_kml2:'GeometryCollection'几何类型不是   支持的SQL状态:XX000上下文:SQL函数“st_askml”语句   1

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

1 个答案:

答案 0 :(得分:2)

如果您的ST_Union操作由于某种原因产生PointsLinestrings,或者您在ST_Union的输入中混合了几何类型,则可能会发生这种情况,即多边形或多边形。如Google KML docs中所述,KML不支持GeometryCollections

你可以通过调用ST_Union两次来解决这个问题,使用ST_Dump将联合几何分开,只选择那些产生多边形或多边形的第一个联合(转储)中的几何 - 你可以组合而不产生一个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的问题,该查询尝试从Point和Linestring中创建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的一部分,所以在联合后,您可以避免在第二种情况下可怕的GeometryCollection。但是,如果您将任何行或点投入混合中,您将返回GeomtryCollection。

编辑2:根据OP的评论,Google KML docs表示KML不支持GeometryCollections。由于无法在不生成GeometryCollection的情况下将Point和/或Line与Polygon组合,因此上述(多)Polygon方法是唯一可用的方法。