我有一个应用程序需要一个地址字符串，将它发送到谷歌地图 API 并获取纬度/经度坐标，然后我想显示该点 X 米范围内的所有用户(纬度/经度存储在我的数据库)，然后我想过滤结果以仅显示拥有某些宠物的用户

所以首先，我有我的模型

class User(UserMixin, Base):
    first_name = Column(Unicode)

    address = Column(Unicode)
    location = Column(Geometry('POINT'))

    pets = relationship('Pet', secondary=user_pets, backref='pets') 

class Pet(Base):
    __tablename__ = 'pets'   
    id = Column(Integer, primary_key=True)
    name = Column(Unicode)

user_pets = Table('user_pets', Base.metadata,
    Column('user_id', Integer, ForeignKey('users.id')),
    Column('pet_id', Integer, ForeignKey('pets.id'))
)

我从 Google API 获取我的纬度/经度并将其存储在我的数据库中，因此我从地址字符串“London England”中获取

POINT (-0.1198244000000000 51.5112138999999871)

这在我的数据库中存储如下:

0101000000544843D7CFACBEBF5AE102756FC14940

现在一切正常，现在阅读 Geoalchemy2 文档，我似乎无法找到一个 exmaple 查询来解决我的问题。

我想要传递的是另一组经/纬度坐标到 Geoalchemy2，然后返回最近的 10 个用户。在查询这个的同时，我还将仅过滤拥有特定宠物的用户(这对于我的查询工作来说不是必需的，但我想展示查询的全部实际作用)。

我不太喜欢在不提供示例查询的情况下回答问题，但我真的不知道应该使用哪些函数来实现所需的结果。

我猜我需要使用“ST_DWithin”或“ST_DFullyWithin”，但我找不到这两个函数的完整示例。谢谢。

所以我知道有一个有效的查询

distance = 10
address_string = "London, England"
results = Geocoder.geocode(address_string)
# load long[1], lat[0] into shapely
center_point = Point(results.coordinates[1], results.coordinates[0])
print center_point
# 'POINT (-0.1198244000000000 51.5112138999999871)'
wkb_element = from_shape(center_point)
users = DBSession.query(User).\
    filter(func.ST_DWithin(User.location,  wkb_element, distance)).all()

生成如下SQL

2013-12-30 15:12:06,445 INFO  [sqlalchemy.engine.base.Engine][Dummy-2] SELECT users.first_name AS users_first_name, users.last_name AS users_last_name, users.phone AS users_phone, users.address AS users_address, users.about AS users_about, ST_AsBinary(users.location) AS users_location, users.profile_image_id AS users_profile_image_id, users.searchable AS users_searchable, users.user_password AS users_user_password, users.registered_date AS users_registered_date, users.id AS users_id, users.last_login_date AS users_last_login_date, users.status AS users_status, users.user_name AS users_user_name, users.email AS users_email, users.security_code AS users_security_code 
FROM users 
WHERE ST_DWithin(users.location, ST_GeomFromWKB(%(ST_GeomFromWKB_1)s, %(ST_GeomFromWKB_2)s), %(param_1)s)
2013-12-30 15:12:06,445 INFO  [sqlalchemy.engine.base.Engine][Dummy-2] {'ST_GeomFromWKB_1': <read-only buffer for 0x7f7d10258f70, size -1, offset 0 at 0x7f7d10258db0>, 'param_1': 10, 'ST_GeomFromWKB_2': -1}

现在这总是会返回我的所有用户，无论距离变量如何，所以我猜有些事情不完全正确，但我不知道为什么。

最佳答案

回答:

单位以度为半径，所以我必须将英里转换为度以获得最佳(粗略)估计。它不需要是精确的:

d = 90
distance = d * 0.014472
#1 mile = 0.014472 degrees  

r1 = -0.1198244
r2 = 51.5112139

# load long[1], lat[0] into shapely
center_point = Point(r1, r2)
# 'POINT (-0.1198244000000000 51.5112138999999871)'

wkb_element = from_shape(center_point)

users = DBSession.query(User).\
    filter(func.ST_DFullyWithin(User.location,  wkb_element, distance)).all()

关于python - Geoalchemy2查询X米范围内的所有用户，我们在Stack Overflow上找到一个类似的问题： https://stackoverflow.com/questions/20803878/