Croisement des données d'inondations avec le LIDAR
Sommaire |
Contexte
Cette page est un sous-thême de la page sur l'Utilisation des données LIDAR pour la directive inondation. Le lecteur est invité à s'y reporté pour comprendre le contexte général sur l'exploitation des données LIDAR
L'objectif de cette partie est de fournir des méthodes pour définir des hauteurs d'eau en zone inondable à partir de jeux de données divers ponctuels, linéiques ou surfaciques et/ou de leur croisement tel que:
- Plus hautes Eaux (PHE)
- Limites de zones inondées
- Résultats de modèles hydrauliques 1D, 1D casier et 2D
- Emprise de zones hydrogéomorphologiques
Principes
Les principes de cette partie vont être les suivants:
- Création d'un modèle numérique de surface en eau (MNSE)
- Différence de raster entre le MNSE et le MNT
- Réalisation d'iso-classes de hateurd 'eau avec la méthode de Classement en iso-valeurs des données LIDAR
Méthode directe dans Qgis
UNIQ6c5ad3a85fa2a077-anyweb-00000000-QINUUNIQ6c5ad3a85fa2a077-anyweb-00000001-QINU
Des explications sont fournies pour une prise en main très rapide sur cette page. Pour un travail plus fin en particulier avec des modèles à casiers, les deux documents ci-contre rentrent plus dans le détail d'un type de procédure.
Création de la couche raster de surface en eau (raster_PHE)
- ouvrir la couche vecteur des PHE
- lancer l'extension interpolation
- choisir la couche vecteur
- choisir l'attribut Altitude (cet attribut doit être de type Réel)
- l'ajouter et choisir le type d'information (point, ligne de structure ou ligne de rupture )
- choisir la méthode d'interpolation (TIN ou IDW)
- définir la taille des cellules en x et y (la taille doit être identique à celle du raster LIDAR où sera réalisé le croisement)
- donner le fichier de sortie
- lancer l'interpolation
Ouverture de la couche raster terrain naturel (raster_MNT)
Voir la page Préalables pour l'utilisation de Qgis-GRASS sur le LIDAR
Création de la couche raster de zone inondable
- lancer le menu raster/calculatrice raster
- taper l'expression raster_PHE-raster_MNT
- choisir une couche de sortie
- laisser le format de sortie GeoTiff (le format raster virtuel n'a pas fonctionné dans ce cas)
- lancer le calcul
Visualisation du résultat
- affecter une palette de couleur à la couche créée (Voir si besoin la page Préalables pour l'utilisation de Qgis-GRASS sur le LIDAR)
- ouvrir un scan IGN
Réalisation de la couche vecteur d'iso-hauteurs
Réalisation d'iso-classes de hateurs d 'eau avec la méthode de Classement en iso-valeurs des données LIDAR
Automatisation dans Qgis-GRASS
Elle a été réalisé dans l'outil DICARTO à partir de la fonction v.to.rast de GRASS. celle-ci permet de convertir des données vecteur (ligne, surface) en raster en affectant des valeurs d'attribut, une cote de casier ou une valeur à une section de calcul.
Pour le traitement des parties 1D, des sections intermédiaires ont été créés entre deux sections de résultats avec une progression linéaire de la forme. Elle ont ensuite été converties en raster. Cette étape est très rapide, permet d'éviter toute triangularisation et s'attache à fournir une cote d'eau uniquement entre deux sections de calcul.
Le travail s'effectue aussi avec une distance entre section afin de créer des plus petits fichiers raster qui sont gerés par des tables d'assemblage.
Liens externes
- http://geoinfo.metier.i2/qgis-r625.html, site intranet MEDDE sur l'outil Qgis avec version ministère, Qgis - Documentation
Qgis - Contributions des services
Qgis - Groupe de travail national
Qgis - Plugin
Qgis - Fiche pratique
Qgis - Faq
Qgis - Le coin des développeurs
Packages d’installation
GRASS
- https://github.com/qgis, les sources du logiciel
- http://forum.qgis.org/index.php, le forum Qgis
- http://grass.osgeo.org/community/support.php, le forum GRASS
- http://grass.osgeo.org/wiki/GRASS-QGIS_relevant_module_list qui permet de connaitre les modules grass présent dans les outils GRASS à travers Qgis et ceux non disponibles mais parfois utiles
- http://www.linux-nantes.org/wws/info/grass-fr, liste de discussion GRASS française, que nous remercions pour sa collaboration
- http://fr.wikipedia.org/wiki/Wikip%C3%A9dia:Atelier_graphique/Didacticiels_cartographiques/Didacticiel_de_cr%C3%A9ation_de_carte_avec_QGIS, aide à la création de carte sous Qgis
Les intervenants
Maîtrise d'Ouvrage: MEDDM/DGPR
- Responsable de l'étude: Sabine Baillarguet
Réalisation: Centre d'Études Techniques de l'Équipement Méditerranée
- Responsable de l'étude: Frédéric Pons
- Intervenants: Mathieu Alquier, Isabelle Roux, Christine Alpin, Yves Bonin, Bertrand Bouteille, Christophe Esposito, Céline Trmal , Christophe Laroche
- Stagiaire: Arnaud Alvarez (TFE ESGT), Adeline Agostini (Stage Master II GERINAT)
Partenariats:
- CPII DO Ouest: René Jonchere
- CETMEF: Bruno Bader
- CETE Sud-Ouest: Yves Nedelec
- DGPR/SRNH/SCHAPI: Christophe Astier, Jean-Luc Souldadié
- IGN: Alain Duperet, Sébastien Saur
Article rédigé par Appui CETE/CETMEF pour la Directive Inondation piloté par le PCI Inondations et Aléas Côtiers
Note : d'autres personnes peuvent avoir contribué au contenu de cet article, [Consultez l'historique]. |