如何有效地联合多个3d多边形并将它们投射到另一张桌子上？

Question

我正在处理大量CityGML三维建筑数据（LoD2），并使用3d城市数据库（及其进口商/出口商工具）导入数据（请参见https://www.3dcitydb.org/3dcitydb/)它在数据库中创建一个新的模式，并使用java例程将XML样式的地理数据转换为PostGIS兼容的数据库输入

然后将每个建筑的信息（空间和非空间）写入不同的表，其中“建筑”表保存建筑元数据，如建筑功能或街道名称，“surface_geometry”表保存所有建筑的实体和三维表面几何。将列“id”、“parent_id”和“root_id”用作PK和FK，以确保一个表的元素可以与另一个表中的一个或多个元素匹配

由于我的大多数分析都需要二维建筑足迹（或鸟瞰多边形），所以我在可用FK上使用st_force2d（）对同一建筑的所有二维表面进行groupby之前创建了查询，以使曲面二维化

是否有任何想法如何改进以下代码以有效处理庞大的数据集？

查询是

CREATE TABLE citydb.building_geom_tmp AS
SELECT a.*, b.bldg_footprint AS geom
FROM Citydb.building AS a
LEFT JOIN (
        SELECT *, ST_Force2D(geometry) AS bldg_footprint FROM citydb.surface_geometry
        ) AS b
ON a.lod2_solid_id = b.root_id;

DROP TABLE IF EXISTS citydb.building_geom ;
CREATE TABLE citydb.building_geom AS
SELECT a.*, b.geom
FROM citydb.building AS a
LEFT JOIN (
        SELECT id, ST_Union(geom) AS geom
        FROM citydb.building_geom_tmp  AS c
        GROUP BY c.id
        ) AS b
ON a.id = b.id;

DROP TABLE IF EXISTS citydb.building_geom_tmp

我正在使用Win 10和< code>SELECT version()和< code > SELECT postgis _ full _ version() give的虚拟机

PostgreSQL 13.2，由Visual C构建1914编译，64位

和

POSTGIS = " 3 . 1 . 1 3 . 1 . 1 "[扩展]PGSQL = " 130 " GEOS = " 3 . 9 . 1-CAPI-1 . 14 . 1 " PROJ = " 7 . 1 . 1 " GDAL = " GDAL 3 . 2 . 1，发布于2020/12/29 " LIBXML = " 2 . 9 . 9 " LIBJSON = " 0.12 " LIBPROTOBUF = " 1 . 2 . 1 " WAGYU = " 0 . 5 . 0(内部)"栅格

分别是。3dcityDB在虚拟机上是4.1.0，上面显示了postgres/postgis版本。

CityGML LoD2数据样本可从德国联邦制图和大地测量局（此处）下载。

此外，我在github repo中向3DCityDB开发人员征求了一些想法。

提前感谢任何想法和最好的问候
安德烈

Answer 1

直接在Github上从开发人员那里得到了答案：https://github.com/3dcitydb/3dcitydb/issues/73