GeoDjango MultiPolygonField 是否应该接受多边形几何？

Question

当我尝试在 MultiPolygonField 上设置 Polygon 时，会引发以下异常：

Traceback (most recent call last):
  File "<console>", line 1, in <module>
  File "/Users/mattrowbum/.virtualenvs/my_env/lib/python3.7/site-packages/django/contrib/gis/db/models/proxy.py", line 75, in __set__
    instance.__class__.__name__, gtype, type(value)))
TypeError: Cannot set Location SpatialProxy (MULTIPOLYGON) with value of type: <class 'django.contrib.gis.geos.polygon.Polygon'>

这是可以理解的，除了GeoDjango tutorial 中的注释指出：

...GeoDjango MultiPolygonField 将接受多边形几何体。

我查看了source of proxy.py，它正在检查值（Polygon）是否是相关几何类（MultiPolygon）的实例。我已经尝试手动进行此检查，以确认 Polygon 不会从 MultiPolygon 继承：

>>> from django.contrib.gis.geos import MultiPolygon, Polygon
>>> ext_coords = ((0, 0), (0, 1), (1, 1), (1, 0), (0, 0))
>>> int_coords = ((0.4, 0.4), (0.4, 0.6), (0.6, 0.6), (0.6, 0.4), (0.4, 0.4))
>>> poly = Polygon(ext_coords, int_coords)
>>> isinstance(poly, Polygon)
True
>>> isinstance(poly, MultiPolygon)
False

我在尝试简化现有存储的MultiPolygon 值时注意到了这一点。下面是我的模型：

from django.contrib.gis.db import models

class Location(models.Model):
    name = models.CharField(max_length=180)
    mpoly = models.MultiPolygonField(geography=True, blank=True, null=True)

我用来简化MultiPolygon 的过程如下。最后一行导致引发异常：

>>> from my_app.models import Location
>>> location = Location.objects.get(pk=1)
>>> geom = location.mpoly
>>> simplified_geom = geom.simplify(0.0002)
>>> location.mpoly = simplified_geom

如果我使用Polygon 创建MultiPolygon，它可以正常工作：

>>> multi = MultiPolygon([simplified_geom,])
>>> location.mpoly = multi

教程中的注释是否具有误导性或我做错了什么？

---

编辑：对几何形状的进一步测试。

原MultiPolygon直接来自模型字段：

>>> geom = Location.mpoly
>>> type(geom)
<class 'django.contrib.gis.geos.collections.MultiPolygon'>
>>> geom.geom_type
'MultiPolygon'
>>> geom.valid
True
>>> geom.srid
4326

应用simplify() 方法：

>>> simplified_geom = geom.simplify(0.0002)
>>> type(simplified_geom)
<class 'django.contrib.gis.geos.polygon.Polygon'>
>>> simplified_geom.geom_type
'Polygon'
>>> simplified_geom.valid
True
>>> simplified_geom.srid
4326

从简化几何创建MultiPolygon：

>>> multi = MultiPolygon([simplified_geom,])
>>> type(multi)
<class 'django.contrib.gis.geos.collections.MultiPolygon'>
>>> multi.geom_type
'MultiPolygon'
>>> multi.valid
True
>>> multi.srid
>>> print(multi.srs)
None

请注意，上面的MultiPolygon 没有 SRID。我认为这可能是被接受的原因。我使用srid=4326 参数创建了一个，但它也被该领域接受。

这是一个非常基本的问题示例：

>>> # This works
>>> location.mpoly = MultiPolygon()
>>> # This doesn't
>>> location.mpoly = Polygon()
Traceback (most recent call last):
  File "<console>", line 1, in <module>
  File "/Users/mattrowbum/.virtualenvs/musicteacher/lib/python3.7/site-packages/django/contrib/gis/db/models/proxy.py", line 75, in __set__
    instance.__class__.__name__, gtype, type(value)))
TypeError: Cannot set Location SpatialProxy (MULTIPOLYGON) with value of type: <class 'django.contrib.gis.geos.polygon.Polygon'>

Answer 1

编辑：

虽然MultiPolygon code 似乎允许将Polygon 存储为MultiPolygon 对象：

class MultiPolygon(GeometryCollection):
    _allowed = Polygon
    _typeid = 6

您提出的问题随着您的描述而上升。

我已经尝试了一些变通方法，唯一令我满意的方法是将您的 geom 字段重构为可以存储任何类型几何图形的通用 GeometryField。

另一个不接触模型的选择是将每个Polygon 转换为MultiPolygon，然后再将其插入字段：

p = Polygon()
location.mpoly = MultiPolygon(p)

在我看来，这值得一个教程更新请求或代码修复的问题。

---

在此处保留之前的答案状态以保持评论的连续性：

The issue is with the geography=True setting 而不是该字段，因为MultiPolygonField 确实接受Polygon 以及MultiPolygon：

geography 类型为使用地理坐标（例如 WGS84 经度/纬度）表示的空间特征提供本机支持。 与几何类型使用的平面不同，geography 类型使用其球面表示数据。在地理列上执行的距离和测量操作会自动采用大​​圆弧计算并返回线性单位。换句话说，当在两个地理上调用 ST_Distance 时，会返回一个以米为单位的值（与在 WGS84 中的几何列上调用度数相反）。

由于您将字段设置为期望地理类型对象，因此如果您尝试传递多边形的非地理表示，您将收到相关错误。