如何使用psql命令通过osm2pgsql按Way OSM ID查询经纬度数组？

没问题，我来帮你搞定这个需求！要从osm2pgsql导入的PostGIS数据库里，通过Way OSM ID提取道路的经纬度（lat, lon）数组，你可以按照下面的步骤操作：

1. 先确认目标表 #

osm2pgsql默认会把道路这类线性要素存入planet_osm_line表；如果导入时用了--roads参数，部分道路可能会在planet_osm_roads表。你可以先跑个简单查询确认目标Way是否存在：

-- 替换123456789为你的Way ID
SELECT name, way FROM planet_osm_line WHERE osm_id = 123456789 LIMIT 1;

2. 核心查询：提取经纬度数组 #

因为osm2pgsql默认导入的是Web墨卡托（EPSG:3857）坐标系，我们需要先转成WGS84经纬度（EPSG:4326），再拆分点并聚合为数组。直接用下面的SQL即可：

SELECT 
  array_agg(ARRAY[ST_Y(point_geom), ST_X(point_geom)]::numeric[]) AS lat_lon_array
FROM (
  SELECT 
    ST_DumpPoints(ST_Transform(way, 4326)) AS point_dump
  FROM planet_osm_line
  WHERE osm_id = 123456789 -- 替换为你的Way OSM ID
) AS points
CROSS JOIN LATERAL (SELECT (point_dump).geom AS point_geom) AS point_geom;

关键部分解释 #

• ST_Transform(way, 4326)：把数据库里的Web墨卡托坐标转换成标准的WGS84经纬度，确保拿到的是真实的lat/lon值
• ST_DumpPoints：把道路的LINESTRING几何对象拆分成单个的点要素
• array_agg：将每个点的[纬度, 经度]数组合并成一个二维数组，结果就是整条道路的经纬度点序列

额外注意事项 #

• 如果你的道路在planet_osm_roads表，只需把查询里的planet_osm_line替换成对应表名即可
• 确保你的PostgreSQL数据库已经启用PostGIS扩展，可通过SELECT PostGIS_Version();验证
• 如果目标ID对应的是OSM Relation（比如复杂的组合道路），需要关联planet_osm_rels和planet_osm_members表查询，不过绝大多数普通道路都是Way类型，上面的查询足够覆盖

内容的提问来源于stack exchange，提问作者madan kandula