1

给定以下（简化的）模型：

from django.contrib.gis.db import models

class City(models.Model):
    center = models.PointField(spatial_index=True, null=True)
    objects = models.GeoManager()

class Place(models.Model):
    city = models.ForeignKey(City, null=True)
    lat = models.FloatField(null=True)
    lng = models.FloatField(null=True)
    objects = models.GeoManager()

暂时忘记应该将 Place 中的 lat/lng 移动到 a PointField()，我试图浏览所有的Places并找到最近的城市。目前，我正在做：

from django.contrib.gis.geos import Point

places = Property.objects.filter(lat__isnull=False, lng__isnull=False)
for place in places:
    point = Point(place.lng, place.lat, srid=4326) # setting srid just to be safe
    closest_city = City.objects.distance(point).order_by('distance')[0]

这会导致以下错误：

DatabaseError: geometry_distance_spheroid: Operation on two GEOMETRIES with different SRIDs

假设 SRID 没有默认为 4326，我srid=4326在上面的代码中包含并验证了所有城市City.center的 SRID 为 4326：

In [6]: [c['center'].srid for c in City.objects.all().values('center')]
Out[6]: [4326, 4326, 4326, ...]

关于可能导致这种情况的任何想法？

更新：

sql 查询的创建方式似乎存在问题。报错后，查看sql显示：

In [9]: from django.db import connection
In [10]: print connection.queries[-1]['sql']
SELECT (ST_distance_sphere("model_city"."center", 
    ST_GeomFromEWKB(E'\\001\\001...\\267C@'::bytea))) AS "distance", 
    "model_city"."id", "model_city"."name", "listing_city"."center" 
    FROM "model_city" ORDER BY "model_city"."name" ASC LIMIT 21

看起来 django 正在将 的point参数distance()转换为 Extended Well-Known Binary。 如果我然后更改ST_GeomFromEWKB为ST_GeomFromText一切正常。 例子：

# SELECT (ST_distance_sphere("listing_city"."center",
          ST_GeomFromText('POINT(-118 38)',4326))) AS "distance", 
          "model_city"."name", "model_city"."center" FROM "model_city" 
          ORDER BY "listing_city"."name" ASC LIMIT 5;

     distance     |    name     |                       center                       
------------------+-------------+----------------------------------------------------
 3124059.73265751 | Akron       | 0101000020E6100000795DBF60376154C01CB62DCA6C8A4440
  3742978.5514446 | Albany      | 0101000020E6100000130CE71A667052C038876BB587534540
 1063596.35270877 | Albuquerque | 0101000020E6100000CC0D863AACA95AC036E7E099D08A4140

我在文档中找不到任何关于如何GeoQuerySet.distance()翻译成 SQL 的内容。我当然可以在查询中使用原始 SQL 来让事情正常工作，但我更愿意将所有内容都很好地保存在 Django 框架中。

4

1 回答 1

0

我认为这个错误：“对具有不同 SRID 的两个 GEOMETRIES 进行操作”

数据库上的“geometry_columns”表在要处理的表名之间设置了不同的 srid

***** 你应该自己改变它 在此处输入图像描述

于 2014-02-25T06:58:58.260 回答