如何从 geoalchemy2 的查询结果中获取 lng lat 值答案

【问题标题】：How to get lng lat value from query results of geoalchemy2如何从 geoalchemy2 的查询结果中获取 lng lat 值
【发布时间】：2014-07-28 08:56:15
【问题描述】：

例如，

class Lake(Base):
     __tablename__ = 'lake'
     id = Column(Integer, primary_key=True)
     name = Column(String)
     geom = Column(Geometry('POLYGON'))
     point = Column(Geometry('Point'))


lake = Lake(name='Orta', geom='POLYGON((3 0,6 0,6 3,3 3,3 0))', point="POINT(2 9)")
query = session.query(Lake).filter(Lake.geom.ST_Contains('POINT(4 1)'))
for lake in query:
     print lake.point

它返回了<WKBElement at 0x2720ed0; '010100000000000000000000400000000000002240'>

我也尝试过 lake.point.ST_X() 但它也没有给出预期的纬度

将值从 WKBElement 转换为可读且有用的格式（例如 (lng, lat)）的正确方法是什么？

谢谢

【问题讨论】：

不是答案，但ST_X() 是经度而不是纬度，来自ST_Y()

标签： python sqlalchemy postgis geoalchemy2

【解决方案1】：

您可以使用 shapely 解析 WKB (well-known binary) 点，甚至其他几何形状。

from shapely import wkb
for lake in query:
    point = wkb.loads(bytes(lake.point.data))
    print point.x, point.y

【讨论】：

【解决方案2】：

http://geoalchemy-2.readthedocs.org/en/0.2.4/spatial_functions.html#geoalchemy2.functions.ST_AsText 是您正在寻找的。这将返回“POINT (lng, lat)”。不过，ST_X 应该可以工作，所以如果它没有返回正确的值，您可能会遇到另一个问题。

【讨论】：

你的意思是lake.point.ST_AsText()吗？它返回，但不是一个值。谢谢
刚试了下，还是一样。我发现 session.scalar(lake.point.ST_AsText()) 会返回预期的结果，但是，它要求数据库进行转换，这是预期的行为吗？
对不起，它要求 DB 进行转换是什么意思？
ST_AsText 是一个 Postgis 函数，所以如果这就是您的问题的意思，那么，是的，这是预期的行为。

【解决方案3】：

扩展约翰的答案，您可以在这样查询时使用ST_AsText() -

import sqlalchemy as db
from geoalchemy2 import Geometry
from geoalchemy2.functions import ST_AsText

# connection, table, and stuff here...

query = db.select(
    [
        mytable.columns.id,
        mytable.columns.name,
        ST_AsText(mytable.columns.geolocation),
    ]
)

在此处查找有关使用函数的更多详细信息 - https://geoalchemy-2.readthedocs.io/en/0.2.6/spatial_functions.html#module-geoalchemy2.functions

【讨论】：