使用 django.contrib.gis.measure.D 时出现 GeoDjango dwithin 错误

首先：RHEL 6.5 上的 Python 2.7.6、Django 1.6.5、Postgres 9.3.4、PostGIS 2.1.3、psycopg2 2.5.3

这是相关模型：

class Location(models.Model):
    name = models.CharField(max_length=255)
    geometry = models.MultiPolygonField(blank=True, default=None, null=True)
    objects = models.GeoManager()  # override the default manager with a GeoManager instance
    parent = models.ForeignKey('self', blank=True, default=None, null=True)

    def __unicode__(self):
        return self.name

这个查询应该有效根据文档:

touching_locations = Location.objects.filter(geometry__dwithin=(location.geometry, D(km=5)))
logging.debug(type(touching_locations))
logging.debug(len(touching_locations))

但事实并非如此。第一个调试调用有效，但第二个调用抛出ValueError:

<class 'django.contrib.gis.db.models.query.GeoQuerySet'>
ValueError: Only numeric values of degree units are allowed on geographic DWithin queries.

如果我通过改变做一个小小的改变D(km=5) to 5:

touching_locations = Location.objects.filter(geometry__dwithin=(location.geometry, 5))
logging.debug(type(touching_locations))
logging.debug(len(touching_locations))

突然间它就起作用了。我得到的输出是这样的：

<class 'django.contrib.gis.db.models.query.GeoQuerySet'>
54

有谁知道为什么这不能按预期工作？这可能是一个错误，还是我犯了一个我没有看到的错误？

[edit]
我认为这可能是 Django 的一个错误。我继续开票here。一旦我弄清楚正确的解决方法是什么，我将在此处添加答案。

我收到了对我提交的票证的回复（https://code.djangoproject.com/ticket/22830）。显然，我发现了一个看似没有记录（或者至少没有明确记录）的问题dwithin查询有Distance对象。一位开发者这样说：

由于您的对象位于地理坐标中（几何字段默认 到 WGS84），您必须提供以度为单位的距离。这是 例如匹配 PostGIS 定义：

boolean ST_DWithin(几何 g1, 几何 g2, 双精度 距离_of_srid);

distance_of_srid 是 WGS84 的度数。所以 5 适合你的 例如表示 5 度，而不是 5 公里！

看起来他们将澄清文档以使这一点更加清晰（太棒了！）。

由于我想要的是 5km，所以我需要将 5km 转换为度。 1度约为111.325公里。因此，1公里=1/111.325度。因此，5 公里约为 0.0449 或约 0.05 度。所以我只需要把我的电话改成这样：

touching_locations = Location.objects.filter(geometry__dwithin=(location.geometry, 0.05))