问将GeoJSON FeatureCollection转换为GeoJSON几何对象

EN

for f in data['features']:
    if f.get("type") == "Feature":
        with open(somefile, "w") as out:
            json.dump(f, out)

import json
import sys

class Splitter:
    def __init__(self):
        self.count = 0

    def check(self, input_file):
        print("Input:", input_file)
        with open(input_file, encoding="utf-8") as file:
            data = json.load(file)
        data_type = data.get("type")
        if data_type == "FeatureCollection":
            for f in data['features']:
                obj_type = f.get("type")
                if obj_type == "Feature":
                    self.export_feature(data, f)
                else:
                    print("WARN: non-feature:", obj_type)
        elif data_type == "Feature":
            print("GeoJSON is already a feature object")
        else:
            print("WARN: unknown/other type:", data_type)

    def export_feature(self, parent, feature):
        self.count += 1
        name = f"{self.count}.geojson"
        print("Output:", name)
        if "crs" in parent and "crs" not in feature:
            # copy coordinate reference system (CRS) definition from parent
            # if not present in feature
            feature["crs"] = parent["crs"]
        with open(name, "w") as out:
            json.dump(feature, out)

##############################################################################
# main #
##############################################################################

if len(sys.argv) < 2:
    print("usage: geojson-split.py <file>")
    sys.exit(0)

splitter = Splitter()
for f in sys.argv[1:]:
    splitter.check(f)