GeoDjango - lon/lat PointFields距离使用几何而不是地理

Question

我正在尝试使用GeoDjango将lat/lon存储在PointField中，然后以公里为单位查询两PointField之间的距离。某些内容无法导致距离计算返回几何距离而不是地理距离。我在模型中使用geography=True，但它似乎没有帮助。我正在使用postgis和postgresql

模型：

from django.contrib.gis.db import models
class Customer(models.Model):
    location = models.CharField(max_length=100)
    gis_location = models.PointField(u"longitude/latitude",
                                     geography=True, 
                                     blank=True, 
                                     null=True)

测试：

from django.contrib.gis.geos import Point
# some set up
customer1.gis_location = Point(-79.3, 43.6)
print('Customer 1:', customer1.gis_location)
customer2.gis_location = Point(-89.2, 48.4)
print('Customer 2:', customer2.gis_location)
distance = customer1.gis_location.distance(customer2.gis_location)
print('Distance: ', distance)

输出：

Customer 1:  SRID=4326;POINT (-79.2999999999999972 43.6000000000000014)  
Customer 2:  SRID=4326;POINT (-89.2000000000000028 48.3999999999999986)  
Distance: 11.002272492535353

这仅返回两点之间的几何距离，而不是地理距离。有没有人对我如何让它返回椭球体上的 KMs 距离有任何建议？

谢谢！

Answer 1

通读文档，我发现GEOSGeometry类中的这种距离计算不尊重SRID，并且与PostGIS数据库直接提供的距离计算不同。

我使用了python libray geopy，它提供了计算距离的distance方法。

from geopy.distance import distance as geopy_distance
geopy_distance(customer1.gis_location, customer2.gis_location).kilometers

不幸的是，我仍然无法弄清楚如何使用 Django 的 GIS 库用两个Point进行距离计算。