使用 sqlalchemy 从 postgis 数据库中获取经度和纬度

Question

我正在使用 sqlalchemy 包对我的 postgis 数据库进行查询，该数据库中填充了城市的 .osm 数据。我想从planet_osm_point 表中检索经度和纬度值。 我声明的 sql 查询如下所示：

      SELECT st_y(st_asewkt(st_transform(way, 4326))) as lat,
             st_x(st_asewkt(st_transform(way, 4326))) as lon
             "addr:housenumber" AS husenumber,
             "addr:street" AS street,
             "addr:postcode" AS postcode
      FROM planet_osm_point

Sqlalchemy 向我抛出此错误：

sqlalchemy.exc.InternalError: (psycopg2.InternalError) FEHLER:  Argument to ST_Y() must be a point

唯一的问题是 ST_Y() 和 ST_X() 函数。

Answer 1

ST_X/ST_Y 返回浮点数。您可以使用浮点数或将它们转换为文本。

在这里使用ST_AsEWKT 是有问题的，因为ST_Y/ST_X 返回floats 和ST_AsEWKT 都需要geometry。

使用你得到的浮点数：

SELECT st_y(st_transform(way, 4326)) AS lat,
             st_x(st_transform(way, 4326)) AS lon
             "addr:housenumber" AS husenumber,
             "addr:street" AS street,
             "addr:postcode" AS postcode
      FROM planet_osm_point

或转换为文本：

SELECT cast(st_y(st_transform(way, 4326)) as text) AS lat,
             cast(st_x(st_transform(way, 4326)) as text) AS lon
             "addr:housenumber" AS husenumber,
             "addr:street" AS street,
             "addr:postcode" AS postcode
      FROM planet_osm_point