如何从geoalchemy2的查询结果中获取经纬度值

Question

例如，

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)）的正确方法是什么？

谢谢

- Chung

---

不是答案，但 ST_X() 是用于经度，而纬度则是使用 ST_Y()。 - Mike T

Answer 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

- Enrico Tiongson

Answer 2

http://geoalchemy-2.readthedocs.org/en/0.2.4/spatial_functions.html#geoalchemy2.functions.ST_AsText 是你要找的内容。它将返回 'POINT (lng, lat)'。ST_X应该也可以工作，如果它没有返回正确的值，那么你可能有另一个问题。

- John Powell

---

你是指lake.point.ST_AsText()吗？它返回了<geoalchemy2.functions.ST_AsText at 0x2bfc050; ST_AsText>，但没有值。谢谢。 - Chung

刚试了一下，结果还是一样的。我发现session.scalar(lake.point.ST_AsText())会返回预期的结果，但它要求DB进行转换，这是预期的行为吗？ - Chung

抱歉，您的意思是它向数据库请求转换吗？ - John Powell

ST_AsText是一个Postgis函数，所以如果您的问题是指这个，那么是的，这是预期的行为。 - John Powell

Answer 3

在扩展John的答案时，您可以在查询时使用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

- Atul Kumar

Answer 4

使用shapely。

from shapely import wkb
for lake in query:
    point = wkb.loads(lake.point.data.tobytes())

    latitude = point.y
    longitude = point.x

来源 https://dev59.com/w2Af5IYBdhLWcg3w9mts#30203761

- Garvit Jain

Answer 5

我是这样做的

class BusDepotGeoFence(models.Base):
    __tablename__ = "bus_geofence_data"

    id = Column(Integer, primary_key=True)
    area = Column(Geometry("POLYGON"))

查询部分如下所示

    query = select(ST_AsGeoJSON(BusDepotGeoFence.area).label("area"))
    results = await db.execute(query)
    results = results.fetchall()
    results = [i[0] for i in results]
    results = [json.loads(result) for result in results]

- Vishesh Mangla