我在开发环境中的 SpatiaLite 上运行 GeoDjango 时遇到地理查询问题。
from django.contrib.gis.db import models
class Test(models.Model):
poly = models.PolygonField()
point = models.PointField()
geom = models.GeometryField()
objects = models.GeoManager()
通过外壳测试：
>>> from geotest.models import Test
>>> from django.contrib.gis.geos import GEOSGeometry
>>>
>>> point = GEOSGeometry("POINT(0 0)")
>>> point
<Point object at 0x105743490>
>>> poly = GEOSGeometry("POLYGON((0 0, 0 1, 1 1, 1 0, 0 0))")
>>> poly
<Polygon object at 0x105743370>
>>>
有了这些定义，让我们尝试一些基本的地理查询。首先contains：
>>> Test.objects.filter(point__within=poly)
Assertion failed: (0), function appendGeometryTaggedText, file WKTWriter.cpp, line 228.
Abort trap
django shell 死了。并与within：
>>> Test.objects.filter(poly__contains=point)
GEOS_ERROR: Geometry must be a Point or LineString
Traceback (most recent call last):
...
GEOSException: Error encountered checking Coordinate Sequence returned from GEOS C function "GEOSGeom_getCoordSeq_r".
>>>
其他组合也会导致各种异常。我必须遗漏一些明显的东西，因为这些都是非常基本的。有任何想法吗？