能否通过OverpassAPI离线获取way信息？基于本地planet.osm的实现方法

你可以通过以下几种方案实现离线获取Way信息，适配不同的使用场景：

方案一：本地部署Overpass API（兼容原有查询语法） #

如果想继续使用熟悉的Overpass查询语句，本地部署Overpass API是最直接的方式，步骤如下：

1. 预处理数据：将Planet OSM数据（建议用.pbf格式，比.osm体积小）转换为Overpass支持的格式，可使用osmconvert工具完成格式转换。
2. 部署Overpass服务：
   • 可通过包管理器安装（如Debian/Ubuntu：sudo apt install overpass-api）
   • 或使用Docker容器快速部署，避免环境配置麻烦
3. 导入数据：执行Overpass的导入命令（如osm3s_query --meta --db-dir=/path/to/overpass/db --import /path/to/planet.osm.pbf），等待数据导入完成
4. 修改Python代码：将Overpass端点指向本地服务，其余逻辑和在线版本一致：

from OSMPythonTools.overpass import Overpass  
# 指向本地Overpass API端点，默认端口8080
overpass = Overpass(endpoint='http://localhost:8080/api/interpreter')
query = "way[highway](around:3,35.356309, 139.598544);(._;>;);out;"
result = overpass.query(query)

方案二：用PyOsmium直接解析OSM文件 #

如果不想部署服务，可直接用Python的pyosmium库解析本地OSM文件，自定义筛选逻辑：

1. 安装依赖：

pip install osmium shapely pyproj

2. 示例代码：实现筛选指定区域内的highway类型Way

import osmium
import shapely.geometry as sg
import pyproj

# 坐标转换：WGS84转UTM（适合计算短距离）
proj_wgs84 = pyproj.CRS('EPSG:4326')
# 东京区域对应UTM带为EPSG:32654，可根据目标点调整
proj_utm = pyproj.CRS('EPSG:32654')
transformer = pyproj.Transformer.from_crs(proj_wgs84, proj_utm, always_xy=True)

# 目标点与3米缓冲区
target_lon, target_lat = 139.598544, 35.356309
target_utm = transformer.transform(target_lon, target_lat)
buffer = sg.Point(target_utm).buffer(3)

# 第一步：收集所有节点的坐标映射
node_coords = {}
class NodeHandler(osmium.SimpleHandler):
    def node(self, n):
        node_coords[n.id] = (n.location.lon, n.location.lat)

node_handler = NodeHandler()
# 加载本地Planet PBF文件，开启locations参数获取节点坐标
node_handler.apply_file('planet.osm.pbf', locations=True)

# 第二步：筛选符合条件的Way
class HighwayWayHandler(osmium.SimpleHandler):
    def __init__(self, node_coords, buffer):
        super().__init__()
        self.node_coords = node_coords
        self.buffer = buffer
        self.matching_ways = []

    def way(self, w):
        # 过滤不含highway标签的Way
        if 'highway' not in w.tags:
            return
        # 构建Way的UTM坐标线串
        coords = []
        for n in w.nodes:
            if n.ref not in self.node_coords:
                continue
            lon, lat = self.node_coords[n.ref]
            coords.append(transformer.transform(lon, lat))
        if len(coords) < 2:
            return
        line = sg.LineString(coords)
        # 检查Way是否与3米缓冲区相交
        if line.intersects(self.buffer):
            self.matching_ways.append({
                'id': w.id,
                'tags': dict(w.tags),
                'node_ids': [n.ref for n in w.nodes]
            })

way_handler = HighwayWayHandler(node_coords, buffer)
way_handler.apply_file('planet.osm.pbf')

# 输出结果
for way in way_handler.matching_ways:
    print(f"Way ID: {way['id']}, Highway Type: {way['tags']['highway']}")

注意：Planet文件体积极大，直接解析会占用大量内存，建议使用分块的OSM数据（如按区域下载的切片），或使用osmium的索引功能优化性能。

方案三：导入PostGIS数据库进行空间查询 #

如果需要频繁执行复杂空间查询，推荐将OSM数据导入PostGIS数据库：

1. 准备环境：安装PostgreSQL+PostGIS扩展，以及osm2pgsql工具
2. 导入数据：

# 创建数据库并启用PostGIS
createdb osm_db
psql -d osm_db -c "CREATE EXTENSION postgis;"
# 导入Planet数据
osm2pgsql -d osm_db planet.osm.pbf

3. 执行SQL查询：通过空间函数筛选目标区域的Way

SELECT w.osm_id, w.tags->'highway' AS highway_type
FROM ways w
JOIN way_nodes wn ON w.osm_id = wn.way_id
JOIN nodes n ON wn.node_id = n.osm_id
WHERE w.tags->'highway' IS NOT NULL
AND ST_DWithin(
    ST_SetSRID(ST_MakePoint(n.lon, n.lat), 4326),
    ST_SetSRID(ST_MakePoint(139.598544, 35.356309), 4326),
    3  -- 距离单位：米
)
GROUP BY w.osm_id, w.tags->'highway';

4. Python连接查询：用psycopg2库连接数据库获取结果

import psycopg2

conn = psycopg2.connect("dbname=osm_db user=your_username")
cur = conn.cursor()
cur.execute("""
SELECT w.osm_id, w.tags->'highway' AS highway_type
FROM ways w
JOIN way_nodes wn ON w.osm_id = wn.way_id
JOIN nodes n ON wn.node_id = n.osm_id
WHERE w.tags->'highway' IS NOT NULL
AND ST_DWithin(
    ST_SetSRID(ST_MakePoint(n.lon, n.lat), 4326),
    ST_SetSRID(ST_MakePoint(139.598544, 35.356309), 4326),
    3
)
GROUP BY w.osm_id, w.tags->'highway';
""")
for row in cur.fetchall():
    print(f"Way ID: {row[0]}, Highway Type: {row[1]}")
cur.close()
conn.close()

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