如何使用python提取OSM地图的建筑轮廓数据及道路数据

import osmnx as ox
import geopandas as gpd

# 获取地图数据
place = "Your Place"  # 地名
tags = {"building": True}
gdf = ox.geometries_from_place(place, tags)

# 筛选出建筑物
buildings = gdf[gdf["building"].notnull()]

# 保存为Shapefile
buildings.to_file("buildings.shp")

使用了osmnx.geometries_from_place()函数来获取指定地点的所有地理要素，然后使用geopandas库来筛选出所有“建筑物”类型的要素，并将其保存为Shapefile文件。

要在OSM地图上获取道路数据，并按照行政边界裁剪，可以使用Python中的OSMnx库。以下是一个示例代码，可以提取指定地点的道路数据，并按照指定的行政边界进行裁剪

import osmnx as ox
import geopandas as gpd

# 获取地图数据
place = "Your Place"  # 地名
graph = ox.graph_from_place(place, network_type="drive")

# 获取地理范围
boundary = ox.geocode_to_gdf(place)
polygon = boundary["geometry"][0]

# 按照行政边界裁剪道路数据
roads = ox.graph_to_gdfs(graph, nodes=False, edges=True)
roads = gpd.clip(roads, polygon)

# 保存为Shapefile
roads.to_file("roads.shp")

使用了osmnx.graph_from_place()函数来获取指定地点的道路网络数据，然后使用osmnx.geocode_to_gdf()函数来获取地点的行政边界。接着，我们使用osmnx.graph_to_gdfs()函数将道路数据转换为GeoDataFrame格式，并使用geopandas.clip()函数按照行政边界进行裁剪，最后将裁剪后的数据保存为Shapefile文件。