使用Python将OSM XML转换为GEOJSON

pip install osmium geojson

import osmium
import geojson

class OSMHandler(osmium.SimpleHandler):
    def __init__(self):
        osmium.SimpleHandler.__init__(self)
        self.features = []

    def node(self, n):
        if 'highway' in n.tags:
            point = geojson.Point((n.location.lon, n.location.lat))
            feature = geojson.Feature(geometry=point, properties={"highway": n.tags['highway']})
            self.features.append(feature)

    def way(self, w):
        if 'highway' in w.tags:
            line = geojson.LineString([(n.location.lon, n.location.lat) for n in w.nodes])
            feature = geojson.Feature(geometry=line, properties={"highway": w.tags['highway']})
            self.features.append(feature)

    def area(self, a):
        if 'highway' in a.tags:
            polygon = geojson.Polygon([[(n.location.lon, n.location.lat) for n in a.outer_ring]])
            feature = geojson.Feature(geometry=polygon, properties={"highway": a.tags['highway']})
            self.features.append(feature)

# 输入的OSM XML文件路径
osm_file = "path/to/osm.xml"
# 输出的GEOJSON文件路径
geojson_file = "path/to/output.geojson"

handler = OSMHandler()
handler.apply_file(osm_file)

feature_collection = geojson.FeatureCollection(handler.features)

with open(geojson_file, 'w') as f:
    geojson.dump(feature_collection, f)