如何将检索到的OSM建筑物数据可视化并转换为shapefile？

Question

我使用以下脚本从OSM检索数据：

import overpy
api = overpy.Overpass()
print('Obtaining Data...')
result = api.query("""
[out:json]
[timeout:25]
;
(
  relation
    ["building"]
    (51.909331730124,4.3151378631592,51.954898210091,4.405689239502);
  
    way
    ["building"]
    (51.909331730124,4.3151378631592,51.954898210091,4.405689239502);
);
out;
>;
out skel qt;
""")

在执行某些空间分析之前，我想要可视化数据并将其保存为shapefile。我想不出如何将其转换/保存为shapefile，有人能给我一些建议吗？

谢谢！

Answer 1

import osmnx as ox 
import ast

point = 'point coordinates'
dist = 'distance in m'
buildings = ox.geometries.geometries_from_point(point, {'building': True}, dist=dist)

并转换为geodataframe：

buildings_save = buildings.applymap(lambda x: str(x) if isinstance(x, list) else x)