地理方案2

Question

我想运行一个查询，返回矩形内的每个点，其中点和矩形基于真实世界的经度和纬度。

这是失败的查询：

results = session.query(Store.id).filter(func.ST_Within(Store.location, func.ST_GeomFromEWKT('SRID=4326;POLYGON((150 -33, 152 -33, 152 -31, 150 -31, 150 -33))')))

它毫无怨言地运行，但是当调用results.first（）时，我看到以下错误和警告：

sqlalchemy.exc.ProgrammingError： （psycopg2.ProgrammingError） 函数st_within（地理、几何）不存在 第 3 行：WHERE ST_Within（store.location， ST_GeomFromEWKT（'SRID=4326;P... ^ 提示：没有函数与给定的名称和参数类型匹配。可能需要添加显式类型强制转换。[SQL： '选择 store.id 作为store_id\n从存储 \nWHERE ST_Within（store.location， ST_GeomFromEWKT（%（ST_GeomFromEWKT_1）s ）） \n LIMIT %（param_1）s'][参数： {'ST_GeomFromEWKT_1'： 'SRID=4326;多边形（（150 -33， 152 -33， 152 -31， 150 -31， 150 -33））'， 'param_1'： 1}]

但是，我可以通过在查询中创建一个虚拟点（这会使每个商店都匹配）来实现查询：

results = session.query(Store.id).filter(func.ST_Within(func.ST_GeomFromEWKT('SRID=4326;POINT(151 -32)'), func.ST_GeomFromEWKT('SRID=4326;POLYGON((150 -33, 152 -33, 152 -31, 150 -31, 150 -33))')))

这表明问题出在我的Store.location字段上，但是我尝试过的所有方法都不起作用。

这是我对位置列的SQLAlchey定义：

location = Column(Geography(geometry_type='POINT', srid=4326))

这是我用来转经度的代码

stores = session.query(Store)
for store in stores:
    store.location = 'SRID=4326;POINT({} {})'.format(store.longitude, store.latitude)
session.commit()

Python告诉我Store.location的类型是“geoalchemy 2 . elements . WKB element”，这正是我从文档中所期望的。

请问有人对如何解决这个问题有什么建议吗？

仅供参考，我正在运行：

• PostgreSQL 9.6.1
• psycopg2 2.6.2
• SQLAlchemy 1.1.4和
• 地球化学2 0.4.0

Answer 1

如果你愿意的话

返回矩形内的每个点，其中点和矩形基于真实世界的经度和纬度

使用geography数据类型，可以使用ST_Interacts和ST_Covers，因为它们支持geography数据类型并提供更准确的结果。

PostGIS文档:“地理上的函数(面积、距离、长度、交叉点等)是使用球体上的弧来计算的。通过考虑世界的球体形状，这些函数提供了更精确的结果。”

另一方面，它们“需要更多的CPU时间来执行”[1]

有一个PostGIS地理支持函数的列表

Answer 2

感谢其他地方（迈克·拜耳和格雷格·贝克）的帮助，我可以发布答案。

问题是:

1. 我的观点是地理类型，我的多边形是几何类型，
2. 许多其他PostGIS功能，包括ST_Within，不支持地理（即它们只支持几何）。

答案是在查询中将“地理”转换为“几何”。以下查询有效：

results = session.query(Store.id).filter(func.ST_Within(cast(Store.location, Geometry), func.ST_GeomFromEWKT('SRID=4326;POLYGON((150 -33, 152 -33, 152 -31, 150 -31, 150 -33))')))