GeoAlchemy2:得到一个点,一点点

Question

考虑以下SQLAalchemy/ GeoAlchemy2ORM和几何字段:

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 shell中的项目时:

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 - 至少,它没有用Shapelyloads解析.

我如何获得该领域的lat/ ？longeom

Answer 1

获取的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)

Answer 2

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
}