在GeoDjango内置地带Z尺寸（强制2D几何体）

Question

在我的项目，我需要从形状文件导入一些几何。

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

当我尝试保存的几何形状，我得到的错误：

“几何具有Z尺寸，但列不”

什么是剥去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

感谢回答，麦克T.

事情是，我需要它使用GeoDjango framework，不必直接访问数据库的PostGIS。

事实上，经过艰苦的工作，我找到了解决办法。我需要从.clone()对象使用OGRGeometry方法。现在，我可以改变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()

Answer 2

如果您有SQL访问，你可以用ST_Force2D (or ST_Force_2D for older versions)得到这个。

您也可以通过修改类型，e.g纠正源表：

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