開発環境の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シェルが死にます。そしてと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".
>>>
他の組み合わせもさまざまな例外を引き起こします。これらは非常に基本的なものなので、私は明白な何かを見逃しているに違いありません。何か案は?