GeoAlchemy2: Получите широту и долготу точки

Рассмотрим следующее SQLAalchemy / GeoAlchemy2 ORM с полем геометрии:

from geoalchemy2 import Geometry, WKTElement

class Item(Base):

    __tablename__ = 'item'

    id = Column(Integer, primary_key=True)
    ...
    geom = Column(Geometry(geometry_type='POINTZ', srid=4326))

Когда я обновляю элемент в оболочке PostgreSQL:

UPDATE item SET geom = st_geomFromText('POINT(2 3 0)', 4326) WHERE id = 5;

Получение поля:

items = session.query(Item).\
    filter(Item.id == 3)

for item in items:
    print item.geom

Дает:

01e9030000000000000000004000000000000008400000000000000000

Это неправильный WKB — по крайней мере, он не анализирует Shapely loads.

Как получить lat/lon поля geom?


postgis gis wkt geoalchemy2 wkb
person Adam Matan    schedule 02.11.2015    source источник

Ответы (2)


arrow_upward
5
arrow_downward

Получение lat, lon через ST_X и ST_Y может быть не самым элегантным подходом, но он работает:

from sqlalchemy import func

items = session.query(Item, 
                      func.st_y(Item.geom), 
                      func.st_x(Item.geom)).\
    filter(Item.id == 3)

for item in items:
    print item.geom

Дает:

(<Item 3>, 3.0, 2.0)
person Adam Matan    schedule 02.11.2015
comment
Как это сделать без доступа к сессии? Я пытаюсь добавить свойства широты и долготы в свою точечную модель для быстрого доступа. - person Nick Sweeting; 25.02.2016

arrow_upward
4
arrow_downward

geoalchemy2 to_shape для преобразования: class:geoalchemy2.types.SpatialElement в геометрию Shapely.

в классе предметов:

from geoalchemy2.shape import to_shape

point = to_shape(self.geo)

return {
            'latitude': point.y,
            'longitude': point.x
        }
person mosi_kha    schedule 09.07.2019