
Modèles Hydrologiques sur aménagements (3) dans QGIS
la chaîne complète modèle 3D → couches hydrologiques QGIS, et non le seul chargement du nuage en couche (qui est traité dans vue3d.py). Dites-moi si c’est ce dernier point qui vous intéressait.
Ce que fait la chaîne, dans l’ordre des messages précédents :
etape01_mnt.py — nuage → MNS/MNT, puis estimation de σ_Z par ajustement de plan sur vos zones de contrôle (et non par le RMS du rapport photogrammétrique), détection du dôme par ajustement d’ordre 2, CHM, fusion avec le MNT LiDAR régional pondérée par la distance au couvert, lissage dont max_diff est calé sur le σ_Z mesuré, brûlage des buses, breaching avec seuil de profondeur à 3σ.
etape02_classif.py — indices RGB, attributs géométriques (CHM, rugosité), texture GLCM, Random Forest avec kappa sur jeu de validation indépendant, repli par seuils si scikit-learn est absent.
etape03_parametres.py — rasters Manning (modulé par le CHM), Green-Ampt, SCS-CN en alternative, pertes initiales, courbes n(h). Contient une correction du double comptage stockage MNT / stockage paramétré.
etape04_routage.py — D-infinity, réseau, bassin versant, TWI, et estimation du temps de concentration avec les durées d’averse à balayer.
etape05_export.py — terrain, polygones Manning, breaklines, plus les deux fonctions de campagne : test_a_blanc() qui établit le plancher de bruit numérique, et carte_delta() qui masque en dessous.
Trois garde-fous sont codés en dur, parce qu’ils correspondent aux erreurs les plus coûteuses : la validation refuse une maille inférieure à 10 cm, le breaching ignore les dépressions sous 3σ, et parametres_utilises.json trace le jeu de paramètres à chaque exécution — sans quoi les deltas entre scénarios ne sont pas défendables.
J’ai vérifié la syntaxe de tous les modules, la validité du YAML, et testé numériquement les fonctions de calcul focal contre une implémentation naïve. Ce qui ne peut pas être testé ici, faute de QGIS installé : les appels processing.run aux fournisseurs WhiteboxTools et PDAL, dont les noms de paramètres varient légèrement selon les versions.
Presented 14 files
Le point d’entrée est config.yaml : renseignez les chemins de la section entrees, vérifiez le CRS, puis lancez depuis la console Python de QGIS. Le reste s’enchaîne.