使用OSMnx的element_type标签提取多边形内way对象失败问题
问题分析与解决方案
核心问题
element_type不能作为查询标签:osmnx.geometries_from_polygon的tags参数仅支持过滤OpenStreetMap原生标签(比如building、highway这类),而element_type是osmnx为返回数据添加的分类属性(用于区分node/way/relation),并非OSM本身存在的标签,所以直接传{"element_type":"way"}会返回空GeoDataFrame。KeyError源于索引结构:osmnx返回的GeoDataFrame中,
element_type是MultiIndex的组成部分,不是普通数据列,因此直接用new_gdf["element_type"]会触发键不存在的错误。
正确实现方式
步骤1:获取目标地理要素
先通过OSM原生标签(比如building=True)获取多边形内的要素:
polygon = disolved_ellipses_gpd.iloc[0]["geometry"] tags = {"building": True} new_gdf = ox.geometries_from_polygon(polygon, tags)
步骤2:筛选way类型条目
利用MultiIndex的get_level_values方法过滤目标类型:
# 提取element_type为way的行 way_gdf = new_gdf.loc[new_gdf.index.get_level_values('element_type') == 'way'] # 清理全空列 way_gdf_cleaned = way_gdf.dropna(axis=1, how="any")
拓展:获取所有way类型要素
如果不需要限定具体OSM标签,可传空字典获取所有要素后再过滤:
new_gdf = ox.geometries_from_polygon(polygon, tags={}) way_gdf = new_gdf.loc[new_gdf.index.get_level_values('element_type') == 'way']
内容的提问来源于stack exchange,提问作者GrumpyGrizzle
相关产品推荐
相关产品推荐

