GeoDjango 上的剥离 Z 维度（强制 2D 几何）

Question

在我的项目中，我需要从 shapefiles 中导入一些几何图形。

其中一些是 MULTIPOLYGON Z 类型，但所有 Z 坐标都是 0 值。

当我尝试保存几何体时，出现错误：

"Geometry has Z dimension but column does not"

剥离 Z 维度的最佳方法是什么？

我的代码：

ds = DataSource(file_path, encoding='ISO-8859-1')
layers = ds[0]

#need something HERE to coerce geometry to 2D

obj=MyModel(geometry=GEOSGeometry(layers[0].geom.hex))
obj.save()

Answer 1

如果您有 SQL 访问权限，您可以通过 ST_Force2D (or ST_Force_2D for older versions) 获得此权限。

您也可以通过修改类型来更正来源table，例如：

ALTER TABLE my_table
  ALTER COLUMN geom TYPE geometry(MultiPolygon,4326) USING ST_Force2D(geom);

Answer 2

感谢您的回复，Mike T。

问题是我需要使用 GeoDjango framework 来创建它，而不是直接访问 Postgis 数据库。

其实经过一番努力，我找到了解决办法。我需要使用 OGRGeometry 对象中的 .clone() 方法。现在，我可以更改 coord_dim 属性。如果我在原始对象中更改 coord_dim，什么也不会发生。

这是我的代码：

ds = DataSource(file_path, encoding='ISO-8859-1')
layers = ds[0]

#HERE IS THE TRICK
new_layer = layers[0].geom.clone() 
new_layer.coord_dim = 2

obj=MyModel(geometry=GEOSGeometry(layers[0].geom.hex))
obj.save()

GeoDjango 上的剥离 Z 维度（强制 2D 几何）

Strip Z dimension on GeoDjango (Force 2D geometry)

gis

django

postgis

geodjango

geos