LUCAS 2018 dataset as a GeoJSON file for simple and easy visualization - https://esdac.jrc.ec.europa.eu/content/lucas-2018-topsoil-data

lucas2018-to-geojson.py

#!/usr/bin/env python3

import sys
import csv
import json
import argparse
from pathlib import Path


def main(args):
    print("reading csv from {}".format(args.input), file=sys.stderr)

    geo = { "type": "FeatureCollection", "features": [] }

    with args.input.open() as src:
        reader = csv.DictReader(src)

        for row in reader:
            lat, lng = row["TH_LAT"], row["TH_LONG"]

            geo["features"] += [{
                "type": "Feature",
                "properties": row,
                "geometry": {
                    "coordinates": [ lng, lat ],
                    "type": "Point",
                }
            }]

    print("writing geojson to {}".format(args.out), file=sys.stderr)

    with args.out.open("w") as dst:
        json.dump(geo, dst, indent=2)


if __name__ == "__main__":
    parser = argparse.ArgumentParser()

    parser.add_argument("input", type=Path, help="path to lucas2018 .csv file")
    parser.add_argument("-o", "--out", type=Path, required=True, help="path to .geojson output file")

    main(parser.parse_args())