问题 如何从geoalchemy2的查询结果中获取lng lat值


例如,

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)?

谢谢


1570
2018-06-07 01:34


起源

不是答案,而是 ST_X() 是经度 - 而不是纬度,来自 ST_Y() - Mike T


答案:


你可以解析WKB(众所周知的二进制文件)点,甚至其他几何形状,使用 身材匀称

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

7
2018-05-13 00:39





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


4
2018-06-07 09:20



你的意思是lake.point.ST_AsText()?它在0x2bfc050返回<geoalchemy2.functions.ST_AsText; ST_AsText>,但不是值。谢谢 - Chung
你试过func.ST_AsText(lake.point)吗? - John Powell
刚试过,它返回了同样的东西。我发现session.scalar(lake.point.ST_AsText())将返回期望结果,但是,它要求DB进行转换,这是预期的行为吗? - Chung
对不起,这是什么意思,它要求DB进行转换? - John Powell
ST_AsText是一个Postgis函数,所以如果这是你的问题的意思,那么,是的,它是预期的行为。 - John Powell