如何通过OSMNX高效获取国家下属所有二级行政区（州）？

问题背景 #

使用OSMNX的geometries_from_place方法获取巴西的州（对应admin_level=4）时，运行4小时仍未完成，但获取小区域（如里约热内卢的社区）时正常。核心原因是大区域查询时，默认参数会请求大量冗余数据，导致效率极低。

优化解决方案 #

1. 精准限定查询标签，减少冗余数据

默认的admin_level=4标签会匹配所有符合级别的要素，包括非边界类数据。加上boundary: administrative标签只限定行政边界，同时指定只加载必要属性，大幅减少数据量：

import osmnx as ox

# 精确限定行政边界标签
tags = {'admin_level': '4', 'boundary': 'administrative'}
# 仅加载需要的属性字段，避免无关数据
gdf = ox.geometries_from_place('Brasil', tags, attrs=['name', 'admin_level', 'iso_3166_2'])

2. 直接使用国家OSM关系ID查询，跳过地理范围解析

先获取目标国家的OSM关系ID（可通过OpenStreetMap官网搜索国家名称，进入对应relation页面获取），比如巴西的关系ID为155846。通过geometries_from_relation直接从国家关系下获取子级行政区，避免解析整个国家地理范围的开销：

import osmnx as ox

# 巴西的OSM关系ID
brazil_relation_id = 155846
tags = {'admin_level': '4', 'boundary': 'administrative'}
# 直接从国家关系节点下查询子级行政区
gdf = ox.geometries_from_relation(brazil_relation_id, tags)

3. 启用缓存与调整请求配置

开启OSMNX的本地缓存，避免重复请求OpenStreetMap服务器；同时可调整请求超时参数，提升稳定性：

import osmnx as ox

# 启用本地缓存，重复查询直接读取缓存文件
ox.config(use_cache=True, cache_folder='./osmnx_cache', timeout=60)
tags = {'admin_level': '4', 'boundary': 'administrative'}
gdf = ox.geometries_from_place('Brasil', tags)

4. 替代方案：使用预处理的行政区划数据集

如果OSMNX查询仍慢，可直接使用GADM、Natural Earth等预生成的行政区划数据集，这些数据已整理好，无需实时从OSM请求：

import geopandas as gpd

# 读取Natural Earth的全球州级数据（geopandas内置数据源）
world_states = gpd.read_file(gpd.datasets.get_path('naturalearth_states'))
# 筛选巴西的州
brazil_states = world_states[world_states['iso_a3'] == 'BRA']

原方法效率低的原因 #

geometries_from_place('Brasil')需要先解析整个巴西的地理范围，再查询该范围内所有admin_level=4的要素。默认参数会获取所有关联属性和几何类型，而巴西的州级数据在OSM中体量较大，导致请求和数据处理时间过长。

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