Python : lire et écrire un GeoJSON
Python : lire et écrire un GeoJSON
Objectif
À la fin, vous aurez reproduit cette étape de bout en bout sur un cas test.
Avant de commencer
- Python 3.10+ et environnement virtuel.
- Dépendances installées pour le script de l’article.
- Un dossier entrée/sortie clairement séparé.
Pas à pas
- Charger le GeoJSON et vérifier la structure
FeatureCollection. - Valider que chaque feature contient une géométrie.
- Modifier ou filtrer les features nécessaires.
- Écrire un nouveau fichier GeoJSON UTF-8.
- Contrôler l’ouverture dans QuPath/QGIS.
À copier-coller
import json
with open('input.geojson', 'r', encoding='utf-8') as f:
gj = json.load(f)
assert gj['type'] == 'FeatureCollection'
features = [ft for ft in gj['features'] if ft.get('geometry')]
with open('output.geojson', 'w', encoding='utf-8') as f:
json.dump({'type': 'FeatureCollection', 'features': features}, f, ensure_ascii=False)
Vérifier que ça marche
- Le script s’exécute sans exception.
- Les fichiers de sortie sont bien créés.
- Le résultat est cohérent sur un petit lot test.
En cas de problème
- Relancer dans un environnement virtuel propre.
- Vérifier chemins d’entrée/sortie et permissions.