Как получить значение lng lat из результатов запроса 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)?

Спасибо

Answer 1

Вы можете анализировать WKB ( хорошо известные двоичные) точки и даже другие геометрические фигуры, используя shapely.

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

Answer 2

http://geoalchemy-2.readthedocs.org/en/0.2.4/spatial_functions.html - это то, что вы ищете. Это вернет 'POINT (lng, lat)'. Однако ST_X должен работать, поэтому у вас может возникнуть другая проблема, если он не возвращает правильное значение.

Answer 3

Использование стройности.

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

    latitude = point.y
    longitude = point.x

Источник /questions/18090868/kak-poluchit-znachenie-lng-lat-iz-rezultatov-zaprosa-geoalchemy2/18090877#18090877

Answer 4

Я сделал это следующим образом

      
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]

Answer 5

Расширяя ответ Джона, вы можете использовать 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