S'abonner à un flux RSS
 

Notice de Cartino : Différence entre versions

De Wikhydro
(Je veux créer ma matrice de débit avec des formules simples)
(Pré-requis)
 
(39 révisions intermédiaires par 4 utilisateurs sont masquées)
Ligne 10 : Ligne 10 :
 
* Guide prise en main Outils DICARTO et utilisation de cet outil avant CARTINO.
 
* Guide prise en main Outils DICARTO et utilisation de cet outil avant CARTINO.
 
L'utilisateur doit avoir, pour l'installation des différents composants, les droits d'administrateurs sur son PC et doit avoir accès à internet. Ensuite pour toute la durée des calculs, les droits d'écriture doivent être donnés sur les répertoires  C:\Program Files (x86)\Quantum Gis\bin ou C:\Program Files \Quantum Gis\bin.  
 
L'utilisateur doit avoir, pour l'installation des différents composants, les droits d'administrateurs sur son PC et doit avoir accès à internet. Ensuite pour toute la durée des calculs, les droits d'écriture doivent être donnés sur les répertoires  C:\Program Files (x86)\Quantum Gis\bin ou C:\Program Files \Quantum Gis\bin.  
De même que pour l'outil DICARTO, CARTINO a été développé sous Matlab. Son utilisation sur un PC classique nécessite le programme MCRInstaller.exe (à récupérer sous le site FTP de la DI)
+
De même que pour l'outil DICARTO, CARTINO a été développé sous Matlab. Son utilisation sur un PC classique nécessite le programme MCRInstaller.exe (à récupérer par exemple ici : http://cete-aix.fr/dwld/NUNIEAU/MCRInstaller.exe)
 
* Créer un répertoire de travail « CARTINO_test1 » sur la racine du disque « C » ou « D ».
 
* Créer un répertoire de travail « CARTINO_test1 » sur la racine du disque « C » ou « D ».
 
=== Installation de Flutor ===
 
=== Installation de Flutor ===
Concernant le logiciel Flutor, il se peut qu'il faille le compiler pour qu'il fonctionne sur votre poste. Pour cela, il faut tout d'abord installer le logiciel MinGW (dossier fourni). A l'installation, cocher l'ensemble des compilateurs (C++ compiler, fortran compiler et objC compiler ).
+
Concernant le logiciel Flutor, il faut le compiler pour qu'il fonctionne sur votre poste. Pour cela, il faut tout d'abord installer le logiciel MinGW (dossier fourni). A l'installation, cocher l'ensemble des compilateurs (C++ compiler, fortran compiler et objC compiler ).
Puis, il faut accéder aux variables d'environnement du PC (panneaux de configuration/paramètres systèmes/onglet avancé/variables d'environnement) et modifier la variable Path, en rajoutant le chemin d'accès à l'exécutable MinGW (;C:\MinGW\bin;).
+
Puis, il faut accéder aux variables d'environnement du PC (sur win 98 : panneaux de configuration/paramètres systèmes/onglet avancé/variables d'environnement. Sur win7 : Panneau de configuration\Système\paramètres systeme avancés\variables d'environnement) et modifier la variable Path, en rajoutant le chemin d'accès à l'exécutable MinGW (;C:\MinGW\bin;).
 +
Redémarrer le PC.
 
Enfin, ouvrir une boite dos (Demarrer/tous les programmes/Accessoires/invite de commande) , se positionner sur le dossier Plani (commande cd « chemin du répertoire de travail») et taper la commande écrite dans le fichier gfortran.txt du répertoire plani.
 
Enfin, ouvrir une boite dos (Demarrer/tous les programmes/Accessoires/invite de commande) , se positionner sur le dossier Plani (commande cd « chemin du répertoire de travail») et taper la commande écrite dans le fichier gfortran.txt du répertoire plani.
 
Le nom de la version est : flutor_cartino_version_20130211_LMnmax.f90
 
Le nom de la version est : flutor_cartino_version_20130211_LMnmax.f90
 
La commande à taper est:
 
La commande à taper est:
gfortran -o flutor_cartino.exe flutor_cartino_version_20130211_LMnmax.f90
+
gfortran -o flutor_cartino.exe flutor_cartino_version_20130211_LMnmax.f90.
 +
Attention, si la commande gfortran n'est pas reconnue c'est que le programme MinGW n'a pas été installé correctement, donc recommencer l'opération si besoin.
 +
 
 
=== Vérification de Flutor et Mascaret ===
 
=== Vérification de Flutor et Mascaret ===
 
Pour vérifier que Flutor et Mascaret fonctionnent, vous devez copier dans le dossier VERIF_FLUTOR_MASCARET votre version de mascaret et de flutor. Ensuite, cliquez sur mascaret.exe et flutor.exe, un fichier Result1.txt (flutor) apparaît ainsi que un fichier avec extension opt pour mascaret.
 
Pour vérifier que Flutor et Mascaret fonctionnent, vous devez copier dans le dossier VERIF_FLUTOR_MASCARET votre version de mascaret et de flutor. Ensuite, cliquez sur mascaret.exe et flutor.exe, un fichier Result1.txt (flutor) apparaît ainsi que un fichier avec extension opt pour mascaret.
 +
 
== Préalable au lancement de Cartino ==
 
== Préalable au lancement de Cartino ==
 
L'utilisateur a besoin de 4 types de données d'entrée, à placer dans le répertoire de travail :  
 
L'utilisateur a besoin de 4 types de données d'entrée, à placer dans le répertoire de travail :  
Ligne 33 : Ligne 37 :
 
* CartinoV?.ctf
 
* CartinoV?.ctf
 
* CARTINO_parametres.txt
 
* CARTINO_parametres.txt
* dossier Plani  contenant les codes de calcul Mascaret et Flutor
+
* Dossier Plani  contenant les codes de calcul Mascaret, Flutor et Hec-Ras
 
L'ensemble de ces fichiers doivent être placés dans le répertoire des données du projet.
 
L'ensemble de ces fichiers doivent être placés dans le répertoire des données du projet.
  
Ligne 52 : Ligne 56 :
  
 
==== Données de type Débit ====
 
==== Données de type Débit ====
Importer les différents fichiers de débit dans Grass (fonction r.in.gdal.qgis ou r.in.gdal) et les renommer : fichiers Q0010 à Q1000 (le nom est important, il faut une lettre et 4 chiffres)
+
[[File:1_Rmapcalc.jpg|200px|right]]
 +
Importer les différents fichiers de débit dans Grass (fonction r.in.gdal.qgis ou r.in.gdal) et les renommer : fichiers Q0010 à Q1000.
 
Attention: cocher ne pas tenir compte de la projection dans les options avancées.
 
Attention: cocher ne pas tenir compte de la projection dans les options avancées.
  
Si l'utilisateur n'est pas satisfait des débits SHYREG (qui semblent parfois un peu faibles) ou si il souhaite travailler avec des débits différents, il a la possibilité de créer de nouvelles données de type débit (ou de modifier des données Shyreg). En effet, l'utilisateur a intérêt à lancer les calculs automatiques avec un débit maximaliste, ainsi cela évitera d'élargir manuellement tous les profils,  dans les cas où on s’aperçoit que le débit n'est pas assez fort.
+
Si l'utilisateur n'est pas satisfait des débits SHYREG (qui semblent parfois un peu faibles) ou s'il souhaite travailler avec des débits différents, il a la possibilité de créer de nouvelles données de type débit (ou de modifier des données Shyreg). En effet, l'utilisateur a intérêt à lancer les calculs automatiques avec un débit maximaliste, ainsi cela évitera d'élargir manuellement tous les profils,  dans les cas où on s’aperçoit que le débit n'est pas assez fort.
 
Dans le cas de Shyreg,  Il est donc conseillé de créer un débit supérieur deux ou trois fois à celui du débit de SHYREG pour Q1000.
 
Dans le cas de Shyreg,  Il est donc conseillé de créer un débit supérieur deux ou trois fois à celui du débit de SHYREG pour Q1000.
 
Procédure :
 
Procédure :
Ligne 62 : Ligne 67 :
 
* Sélectionner le fichier Q1000, puis Clic Droit et « Ajouter la carte sélectionné à la carte »
 
* Sélectionner le fichier Q1000, puis Clic Droit et « Ajouter la carte sélectionné à la carte »
 
* Cliquer sur [[File:0_Region.jpg|20px]] afin d'éditer la bonne région avec la bonne résolution.
 
* Cliquer sur [[File:0_Region.jpg|20px]] afin d'éditer la bonne région avec la bonne résolution.
[[File:1_Rmapcalc.jpg|200px|right]]
 
 
* Aller dans la « liste des modules » et chercher la fonction « r.mapcalc ». L'objectif est d'avoir, dans le cas ci-dessous, un débit trois fois supérieur à celui de Shyreg. Attention car dans le nom du rendu, il est préférable d'écrire Q1000_3 plutôt que Q1000X3.
 
* Aller dans la « liste des modules » et chercher la fonction « r.mapcalc ». L'objectif est d'avoir, dans le cas ci-dessous, un débit trois fois supérieur à celui de Shyreg. Attention car dans le nom du rendu, il est préférable d'écrire Q1000_3 plutôt que Q1000X3.
  
 
=== Renseignement du fichier de paramètres ===
 
=== Renseignement du fichier de paramètres ===
 
<big>'''MODIFICATION DE LA STRUCTURE. L'ANCIENNE FONCTIONNE MAIS VOUS POUVEZ GERER D'AUTRES PARAMETRES <br>'''</big>
 
<big>'''MODIFICATION DE LA STRUCTURE. L'ANCIENNE FONCTIONNE MAIS VOUS POUVEZ GERER D'AUTRES PARAMETRES <br>'''</big>
Le fichier CARTINO_parametres.txt doit être renseigné avant le lancement de l'application CARTINO, ce fichier texte se compose des éléments suivants (les 4 premières lignes sont équivalentes à l'application dicarto):
+
Le fichier CARTINO_parametres.txt doit être renseigné avant le lancement de l'application CARTINO, ce fichier texte se compose des éléments suivants (les 4 premières lignes sont équivalentes à l'application dicarto):<br>
Les valeurs en rouge doivent être modifiées.
+
Les valeurs en <span style="color: #FF0000;">'''rouge'''</span> doivent être modifiées.<br>
Les valeurs en vert sont des options.
+
Les valeurs en <span style="color: #33FF00;">'''vert'''</span> sont des options.<br>
 +
 
 +
 
 +
%---- Création du Listing (1er chiffre) et Suppression des fichiers temporaires GRASS (2ème chiffre) - 0 Non 1 Oui<br>
 +
<span style="color: #33FF00;">'''0 1 1'''</span> <br>
 +
%---- Dossier d'installation de Qgis-GRASS<br>
 +
<span style="color: #FF0000;">'''C:\Program Files (x86)\Quantum Gis\bin'''</span><br>
 +
%---- Secteur GRASS - Chemin du répertoire avec les données au format GRASS<br>
 +
<span style="color: #FF0000;">'''C:\GRASSDATA\CARTINO\Exemple'''</span><br>
 +
%---- Chemin du répertoire des données au format standard (shape, asc...),Attention, vous ne pouvez pas travailler sur C: ou D: directement!<br>
 +
<span style="color: #FF0000;">'''D:\FormationQgisGRASS'''</span><br>
 +
%---- CODE: Code Flutor et Mascaret O Oui 1 Non => CODE - Etape 4<br>
 +
<span style="color: #33FF00;">'''1 1'''</span><br>
 +
%---- TypBief: Type de Gestion du réseau (Paramètres non actif) et longueur minimale de la partie aval- Etape 4<br>
 +
1 500<br>
 +
%---- nomMNTasc: Nom du Raster MNT dans GRASS =>nomMNTasc - Etape 4<br>
 +
<span style="color: #FF0000;">'''MNT1m'''</span><br>
 +
%---- nomSD: Nom du Raster EXZECO dans GRASS (utile pour connaitre les parties amont-aval avec la surface drainée=> fichierExze  - Etape 3<br>
 +
<span style="color: #FF0000;">'''SDEXZ'''</span><br>
 +
%---- nomQasc: Nombre et Noms du (des) Raster(s) Debit(s) à mettre du + fort au + faible =>nomQasc  - Etape 4<br>
 +
<span style="color: #FF0000;">'''5'''</span><br>
 +
<span style="color: #FF0000;">'''Q1000'''</span><br>
 +
<span style="color: #FF0000;">'''Q0300'''</span><br>
 +
<span style="color: #FF0000;">'''Q0100'''</span><br>
 +
<span style="color: #FF0000;">'''Q0030'''</span><br>
 +
<span style="color: #FF0000;">'''Q0010'''</span><br>
 +
%---- nomEnvLisse: Nom de résultat pour l'enveloppe lissée et nom de l'enveloppe utile aux calculs de distance squelette-enveloppe => nomEnvLisse  - Etape 2 et 4<br>
 +
1<br>
 +
EnveloppeLisse<br>
 +
%---- PasCoupExz: [Resolution du resultat EXZECO => PasExzeco] [Distance de coupure pour le squelette => distcoupsque] - Etape 2<br>
 +
25 5<br>
 +
%---- Reso: Resolution du raster MNT et Raster DEBIT=> Reso - Etape 4<br>
 +
<span style="color: #FF0000;">'''1 50'''</span><br>
 +
%---- DistCoup: Distance de coupure du profil en travers pour le MNT et Debit=> coup - Etape 4<br>
 +
<span style="color: #33FF00;">'''2 10'''</span><br>
 +
%---- PasTri: Pas de Pré-calcul des distance Axe de calcul - Enveloppe =>PasTri - Etape 4<br>
 +
<span style="color: #33FF00;">'''75'''</span><br>
 +
%---- Strickler: Valeur des Strickler Lit mineur et majeur (le majeur n'est pas encore intégré) =>Strickler - Etape 4<br>
 +
<span style="color: #33FF00;">'''20 15'''</span><br>
 +
%---- DPNC: Distance entre profils non constantes:  [Ratio largeur de l'enveloppe - Longueur entre profil => factlarg] [Distance minimale entre 2 traces => distmin] [Distance maximale entre 2 traces => distmax]  <br>
 +
<span style="color: #33FF00;">'''1 50 500'''</span><br>
 +
<span style="color: #33FF00;">'''0.33 10 500'''</span><br>
 +
%---- ValFiltre: [lissage 1ères largeurs des profils en travers => nblisse] [lissage des angles des profils => nbfiltre] [Lissage des altitudes des profils en travers => nlisse]<br>
 +
<span style="color: #33FF00;">'''5 3 1'''</span><br>
 +
%---- Elarg: Ration, Valeur Minimum, Valeur maximum d'élargissement et largeur max du profil (en km paramètres non actif) pour l'enveloppe Initiale => PourcElargissement ValeurElargissement ValeurMaxElargissement%---- Ration, Valeur Minimum, Valeur maximum d'élargissement et largeur max du profil (en km paramètres non actif) suite aux calculs => PourcElargissement ValeurElargissement ValeurMaxElargissement<br>
 +
<span style="color: #33FF00;">'''0.33 5 50 10000'''</span><br>
 +
<span style="color: #33FF00;">'''0.2 10 150 10000'''</span><br>
 +
%---- Nombre de tour maximum pour agrandir les profils - Nombre de tours maximum pour lisser les profils à chaque agrandissement<br>
 +
<span style="color: #33FF00;">'''15 20'''</span><br>
 +
 
  
%---- Création du Listing (1er chiffre) et Suppression des fichiers temporaires GRASS (2ème chiffre) - 0 Non 1 Oui <br>
 
0 1 <br>
 
%---- Dossier d'installation de Qgis-GRASS <br>
 
C:\Program Files (x86)\Quantum Gis\bin <br>
 
%---- Secteur GRASS - Chemin du répertoire avec les données au format GRASS <br>
 
C:\GrassData\test1\test1 <br>
 
%---- Chemin du répertoire de travail contenant les données au format standard (shape, asc...),Attention, vous ne pouvez pas travailler sur C: ou D: directement! <br>
 
D:\CARTINO\test1 <br>
 
%---- CODE: Code Flutor et Mascaret O Oui 1 Non => CODE - Etape 4 <br>
 
1 1 <br>
 
(si on prend 1 et 2, il crée le fichier mascaret mais ne fait pas de calcul) <br>
 
%---- TypBief: Type de Gestion du réseau (Paramètres non actif) '''et longueur minimale de la partie aval-''' Etape 4 <br>
 
1 500 <br>
 
(si on met 2, il fait des plus petits biefs pour limiter les superpositions de calcul) <br>
 
%---- nomMNTasc: Nom du Raster MNT dans GRASS =>nomMNTasc - Etape 4 <br>
 
MNTtest1_G <br>
 
%---- nomSD: Nom du Raster EXZECO dans GRASS (utile pour connaitre les parties amont-aval avec la surface drainée=> fichierExze  - Etape 3 <br>
 
RExtest1_G <br>
 
%---- nomQasc: Nombre et Noms du (des) Raster(s) Debit(s) dans GRASS à mettre du + faible au + fort =>nomQasc  - Etape 4 <br>
 
5 <br>
 
Q0010 <br>
 
Q0030 <br>
 
Q0100 <br>
 
Q0300 <br>
 
Q1000 <br>
 
%---- nomEnvLisse: Nom de résultat pour l'enveloppe lissée et nom de l'enveloppe utile aux calculs de distance squelette-enveloppe => nomEnvLisse  - Etape 2 et 4 <br>
 
1 <br>
 
EnveloppeLisse <br>
 
%---- PasCoupExz: [Resolution du resultat EXZECO => PasExzeco] [Distance de coupure pour le squelette => distcoupsque] - Etape 2 <br>
 
25 5 <br>
 
%---- Reso: Resolution du raster MNT et Raster DEBIT=> Reso - Etape 4 <br>
 
1 50 =-> a priori inutile <br>
 
%---- DistCoup: Distance de coupure du profil en travers pour le MNT et Debit=> coup - Etape 4 <br>
 
2 10 <br>
 
%---- PasTri: Pas de Pré-calcul des distance Axe de calcul - Enveloppe =>PasTri - Etape 4 <br>
 
24 <br>
 
%---- Strickler: Valeur des Strickler Lit mineur et majeur (le majeur n'est pas encore intégré) =>Strickler - Etape 4 <br>
 
20 15 <br>
 
%---- DPNC: Distance entre profils non constantes:  [Ratio largeur de l'enveloppe - Longueur entre profil => factlarg] [Distance minimale entre 2 traces => distmin] [Distance maximale entre 2 traces => distmax]  - 2 lignes, le 1er coup pour calculer une emprise plus exacte que exzeco basé sur un premier calcul Flutor, le deuxième sur ce résultat avec la procédure antérieure<br>
 
1 50 500 <br>
 
0.33 10 1000<br>
 
%---- ValFiltre: [lissage 1ères largeurs des profils en travers => nblisse] [lissage des angles des profils => nbfiltre] [Lissage des altitudes des profils en travers => nlisse] <br>
 
5 3 1 <br>
 
%---- Elarg: Ratio, Valeur Minimum, Valeur maximum d'élargissement, Largeur max du profil (en m paramètres non actif) pour l'enveloppe Initiale => PourcElargissement ValeurElargissement  ValeurMaxElargissement%---- Ratio, Valeur Minimum, Valeur maximum d'élargissement, Largeur max du profil (en m paramètres non actif) suite aux calculs => PourcElargissement ValeurElargissement ValeurMaxElargissement <br>
 
0.05 5 10 3000 <br>
 
0.33 25 100 25000<br>
 
%---- Nombre de tour maximum pour agrandir les profils - Nombre de tours maximum pour lisser les profils à chaque agrandissement <br>
 
10 20 <br>
 
 
==== Dessin d'explications des paramètres ====
 
==== Dessin d'explications des paramètres ====
 
[[File:2_Dessin_expli_parametres.jpg|300px|right]]
 
[[File:2_Dessin_expli_parametres.jpg|300px|right]]
Ligne 137 : Ligne 142 :
 
* Un deuxième calcul qui reprend la largeur obtenu à la fin de premier et qui recalcule les positions des profils à pas non constants et à largeurs non constantes.
 
* Un deuxième calcul qui reprend la largeur obtenu à la fin de premier et qui recalcule les positions des profils à pas non constants et à largeurs non constantes.
 
Ces deux passages permettent un gain de qualité, moins de croisements et aussi un gain en temps.</big>'''''
 
Ces deux passages permettent un gain de qualité, moins de croisements et aussi un gain en temps.</big>'''''
 +
==== Explication du paramètre CODE ====
 +
Le paramètre CODE permet de faire des calculs avec Flutor et/ou Mascaret et de faire les cartographies avec Flutor et/ou Mascaret
 +
{| class="wikitable"
 +
|-
 +
! CODE Flutor-Mascaret !! Condition aval !! Calcul !! Cartographie
 +
|-
 +
| 1 1 || Caval fournie par Flutor || Flutor et Mascaret || Flutor si convergence sinon Mascaret
 +
|-
 +
| 2 1 || Caval fournie par Flutor || Flutor et Mascaret || Mascaret si convergence sinon Flutor
 +
|-
 +
| 1 2 || Caval fournie par Flutor || Flutor, pas de calcul Mascaret, juste  les fichiers d'entrée || Flutor si convergence
 +
|-
 +
| 1 0 || Caval fournie par Flutor || Flutor || Flutor si convergence
 +
|-
 +
| 1 0 || Caval fournie par un calcul Strickler peu fiable || Mascaret || Mascaret si convergence
 +
|}
 +
ATTENTION Les résultats sont parfois instables sur de très fortes pentes, il est toujours préférable de vérifier les lignes d'eau provenant des 2 calculs.
 +
 +
'''De manière générale le calcul Mascaret est à privilégier pour la cartographie : options 2 1'''
 +
 +
Si l'utilisateur veut également lancer un calcul HEC-RAS, il faut rajouter un troisième chiffre 1 sur cette ligne.
  
 
== Lancement de Cartino ==
 
== Lancement de Cartino ==
Ligne 195 : Ligne 221 :
  
 
On obtient ainsi un squelette nettoyé. Sauvegarder les changements en sortant du mode édition.
 
On obtient ainsi un squelette nettoyé. Sauvegarder les changements en sortant du mode édition.
Pour plus de sécurité, vous pouvez au départ enregistrer votre squelette sous un autre nom afin de ne pas avoir à recommencer l'étape de squeletisation en cas de fausses manipulations.
+
Pour plus de sécurité, vous pouvez au départ enregistrer votre squelette sous un autre nom afin de ne pas avoir à recommencer l'étape de squelettisation en cas de fausses manipulations.
  
 
=== Création du réseau des biefs de calculs ===
 
=== Création du réseau des biefs de calculs ===
Ligne 295 : Ligne 321 :
 
* les 3 fichiers cartino (cartino.exe, cartino.ctf, cartino_parametres.txt)
 
* les 3 fichiers cartino (cartino.exe, cartino.ctf, cartino_parametres.txt)
 
* le dossier plani
 
* le dossier plani
* les fichiers shape commençant par « BVBiefCalculs »  
+
* les fichiers shape commençant par « BVBiefCalculs » qui correspondent aux axes de modélisation
* le(s) fichier(s) vecteur(s) de profils que l'on va modifier.
+
* les fichiers vecteurs de profils que l'on va modifier.
  
 
Dans le fichier cartino_parametres.txt, il faut changer le chemin du répertoire puisqu'on a créé un nouveau répertoire « reprise »
 
Dans le fichier cartino_parametres.txt, il faut changer le chemin du répertoire puisqu'on a créé un nouveau répertoire « reprise »
Ligne 378 : Ligne 404 :
 
[[File:9_ExReprise3.jpg|400px|center]]<br>
 
[[File:9_ExReprise3.jpg|400px|center]]<br>
 
Sur notre exemple, on peut s'apercevoir que la zone inondable a bien été prolongée au niveau du profil 39. Par ailleurs nous n'avons plus de zone inondable en rive droite (au niveau du profil 41)
 
Sur notre exemple, on peut s'apercevoir que la zone inondable a bien été prolongée au niveau du profil 39. Par ailleurs nous n'avons plus de zone inondable en rive droite (au niveau du profil 41)
 +
 +
== Covadisation ==
 +
=== Maximum de raster, vectorisation, suppression des petites zones, suppression des créneaux, mise au format COVADIS DI table ISO_HT... ===
 +
En cours d'écriture
 +
Les résultats de la phase cartographie de la directive inondation doivent suivre un standard SIG fournit par la commission COVADIS et disponible sur http://www.cnig.gouv.fr/Front/index.php?RID=154.
 +
Cette page permet de créer à partir des fichiers résultats sous forme raster les formats demandés. ces éléments sont une combinaison des différents modules de DICARTO http://wikhydro.developpement-durable.gouv.fr/index.php/Notice_de_DICARTO:
 +
* Création du max des résultats
 +
* Vectorisation avec options
 +
* Suppression des petites aires
 +
* Suppression des créneaux
 +
* Création des champs COVADIS
 +
 +
Les étapes sont les suivantes:
 +
* choisir les suppressions des petites aires (de 1 à 100m² par pas de 1m par exemple (1:1:100) ou de 25 à 625 m² par pas de 25m2 (25:25:625)
 +
* donner la résolution du raster.
 +
 +
Cartino va alors créer un fichier covadisation.txt dans lequel seront automatiquement listés l'ensemble des fichiers asc présents dans le répertoire de données (un nom de fichier asc par ligne).
 +
Il convient alors de modifier ce fichier de la manière suivante sous Notepad++:
 +
*Pour chaque bassin versant, s'il existe plusieurs biefs (Bi001,Bi002...), il faut les regrouper en une seule ligne par la commande Bi*.
 +
*A chaque fin de ligne, on rajoute un ";" et un numéro qui correspond au choix du scénario (1 pour fréquent, 2 pour moyen, 4 pour extrême).
 +
*Au final on a donc une seule ligne par scénario pour chaque bassin versant.
 +
Après avoir enregistré le fichier covadisation.txt, lancer la covadisation (bouton "FAIT")
 +
 +
 +
Une fois cette étape terminée, il est possible de combiner les fichiers shape résultats sur tous les bassins versants.
 +
* Qgis Menu vecteur\Outils de Géométrie\Fusionner les shapefile en un seul
 +
 +
=== Finalisation de la table ISO_HT pour mapinfo ===
 +
Pour finaliser le travail dans mapinfo au besoin, il faut convertir les shapefile en mapinfo, le mieux est d'utiliser le traducteur universel.
 +
En utilisant la deuxième option, et en sélectionnant les fichiers tab que l'on veut modifier, il suffit de copier les lignes de commande mapbasic dans le menu Option\Fenetre mapbasic.
 +
 +
La routine Modif_rvb.mbx disponible dans le zip de CARTINO permet ensuite de mettre les couleur de manière automatique dans mapinfo.
 +
 +
Un travail est à faire au final pour combiner les données issues de cartino avec celles provenant des autres sources de données (PPRi...) afin d'avoir des identifiants uniques respectant le format COVADIS.
 +
 +
=== Création de la table INONDABLE_S ===
 +
Le troisième menu permet de créer la table N_[prefixTRI]_INONDABLE_S... à partir de fichier shape ou mapinfo.
 +
Il ne faut qu'un scenario par table.
  
 
== Améliorations prévues ==
 
== Améliorations prévues ==
Ligne 412 : Ligne 476 :
 
Valeur Nord +2 m
 
Valeur Nord +2 m
 
Valeur Sud -3 m
 
Valeur Sud -3 m
Valeur Est -2 m
+
Valeur Ouest -2 m
Valeur Ouest +3 m
+
Valeur Est +3 m
  
  
Ligne 427 : Ligne 491 :
  
 
=== J'ai un semi de points et je veux faire un raster ===
 
=== J'ai un semi de points et je veux faire un raster ===
Vous pouvez tester l'option r.in.xyz qui convertit des points en raster en ayant au préalable chois votre région (bords à 0, résolution à 5m par exemple). Ensuite faire r.fillnulls pour combler les vides.
+
Vous pouvez tester l'option r.in.xyz qui convertit des points en raster en ayant au préalable choisi votre région (bords à 0, résolution à 5m par exemple). Ensuite faire r.fillnulls pour combler les vides.
 +
 
 
=== J'ai des rasters avec des projections différentes ===
 
=== J'ai des rasters avec des projections différentes ===
 
On peut être amené à faire travailler sur un territoire avec des MNT qui ne sont pas en RGF93.
 
On peut être amené à faire travailler sur un territoire avec des MNT qui ne sont pas en RGF93.
Ligne 439 : Ligne 504 :
  
 
=== Je veux créer ma matrice de débit avec des formules simples ===
 
=== Je veux créer ma matrice de débit avec des formules simples ===
C'est possible de le faire en partant de la donnée EXZECO raster (Fichier FRD_TRI_MARSEILLE_sol.asc) qui est une surface et en appliquant des calculatrices raster (r.mapcalc) pour faire des choses comme Q=AS0^,8 avec S correspondant à EXZECO<br>
+
[[File:q_s0.8.jpg|150px|droite]]
Avant tout, il faut vérifier la région Gras et la résolution
+
C'est possible de le faire en partant de la donnée EXZECO raster (Fichier FRD_TRI_MARSEILLE_sol.asc) qui est une surface drainée. Il faut ensuite appliquer la calculatrice raster (r.mapcalc) pour faire des calculs comme Q=A.S0^,8 avec S correspondant à EXZECO<br>
Ouvrir le fichier "FRD_TRI_MARSEILLE_sol.asc" dans Q-Gis et l'importer dans Grass en le nommant "surface" par exemple.
+
Avant tout, il faut vérifier la région Grass et la résolution<br>
Dans les outils Grass, faire "rafraichir" et "surface" apparait
+
Ouvrir le fichier "FRD_TRI_MARSEILLE_sol.asc" dans Q-Gis et l'importer dans Grass en le nommant "surface" par exemple.<br>
Double-clic sur "surface" qui apparait dans la fençetre Q-Gis
+
Dans les outils Grass, faire "rafraichir" et "surface" apparait.<br>
 +
Double-clic sur "surface" qui apparait dans la fenêtre Q-Gis.<br>
 
A partir de là, il suffit d'utiliser la calculatrice des outils Grass "r.mapcalc"
 
A partir de là, il suffit d'utiliser la calculatrice des outils Grass "r.mapcalc"
  
Ligne 492 : Ligne 558 :
 
L'importer dans grass (v.in.ogr), le convertir en raster (v.to.rast) et ensuite faire un r.series par exemple avec l'option minimum entre votre MNT et le polygone converti en raster.
 
L'importer dans grass (v.in.ogr), le convertir en raster (v.to.rast) et ensuite faire un r.series par exemple avec l'option minimum entre votre MNT et le polygone converti en raster.
 
Dans Qgis, l'extension Point Simpling Tool disponible sur le dépot QGIS Official Repository 2 permet aussi de faire ce travail.
 
Dans Qgis, l'extension Point Simpling Tool disponible sur le dépot QGIS Official Repository 2 permet aussi de faire ce travail.
 +
 +
=== Méthode appliquée pour s’affranchir des deux premières étapes du modèle CARTINO PC ===
 +
 +
La procédure proposée ci-dessous permet à un utilisateur d’utiliser son réseau hydrographique et de s’affranchir totalement de la source de donnée EXZECO.
 +
 +
Prétraitement CARTINO
 +
 +
Lors de l’intégration d’un raster à surface drainée autre que celui d’EXZECO, il est possible d’être confronté à des problèmes de résolution N/S et E/W sous CARTINO. En imposant une emprise et une résolution connues au raster des surfaces drainées, on s’affranchit de cette erreur.
 +
 +
*Pour commencer, il faut construire un vecteur correspondant à l’emprise du raster SHYREG : Vecteur / Outil de recherche / Créer un polygone à partir de l’emprise de la couche. L’emprise obtenue est utilisée afin d’y extraire les nœuds du vecteur : Vecteur / Outil de géométrie / Extraire les nœuds.
 +
 +
*Parallèlement, à cette manipulation, il faut un réseau hydrographique répondant aux contraintes de CARTINO dont le faite que ses entités soient couvertes par le raster MNT et l’emprise nouvellement créée. Il faut extraire les nœuds de ce réseau : Vecteur / Outil de géométrie / Extraire les nœuds.
 +
 +
*Les deux vecteurs de points, l’un représentant l’emprise de SHYREG et l’autre la position du réseau hydrographique, sont unis : Vecteur / Outil de géotraitement / Unir.
 +
 +
*Puis, ce vecteur est modifié lors d’une phase d’édition. Un champ exutoire est ajouté, grâce à la calculatrice de champ, pour l’ensemble des entités avec la valeur « 0 ».
 +
 +
[[File:calculatrice_de_champ.png|200px|center]]
 +
 +
*Manuellement, les exutoires seront déterminés. Cela demande à l’utilisateur de connaître l’orientation de son réseau hydrographique. Il faut modifier la valeur des deux derniers points de l’exutoire de chaque bassin versant en attribuant la valeur « 1 ».
 +
 +
[[File:verification_attributaire.png|200px|center]]
 +
Ex : s’il y a 6 bassins versants, douze lignes sont modifiées correspondant aux deux derniers points de chaque branche du réseau hydrographique.
 +
 +
*Avant d’interpoler ce vecteur, il est important de noter les informations contenues dans le raster SHYREG tel que l’emprise de la couche, la taille des pixels et la dimension du raster. Ces information sont disponibles dans les propriétés de la couche au niveau de l’onglet « métadonnées ».
 +
 +
Ex : les valeurs ci-dessous permettent d’avoir un aperçu du type d’information à noter
 +
[[File:tableau_des_donnees_a_conserver.png|200px|center]]
 +
 +
*Suite à cette étape, il est possible de convertir le vecteur de points en raster à partir de l’information exutoire précédemment renseignée : Raster / Interpolation / Interpolation.
 +
 +
[[File:extension_interpolation.png|200px|center]]
 +
 +
Modèle CARTINO
 +
 +
Suite à cette note, l’utilisateur de CARTINO peut commencer à partir de l’étape 3 « Création du réseau des biefs de calculs ».
 +
 +
Au préalable, il devra charger une enveloppe lisse sous GRASS par le biais de la fonction v.in.ogr pouvant correspondre à un buffer de son réseau hydrographique. Le nom du fichier sous GRASS doit correspondre à l’information contenue dans le fichier CARTINO_paramètre.txt ; par défaut « Enveloppelisse ».
 +
 +
L’ensemble de ses manipulations aboutissent à une substitution du raster EXZECO et à l’emploie d’un réseau hydrographique extérieur au modèle.
 +
  
 
{{Auteur|NomAuteur=Mathieu Alquier}} {{Auteur|NomAuteur=Frédéric Pons}}
 
{{Auteur|NomAuteur=Mathieu Alquier}} {{Auteur|NomAuteur=Frédéric Pons}}
 +
{{Auteur|NomAuteur=Christophe Esposito}} {{Auteur|NomAuteur=Guillaume Le Bihan}}
  
 
[[Catégorie:Cartographie_des_zones_inondables]]
 
[[Catégorie:Cartographie_des_zones_inondables]]

Version actuelle en date du 27 avril 2017 à 12:33

Cette page est associée à la description de CARTINO PC et plus généralement aux outils définis dans la page Utilisation des données LIDAR pour la directive inondation.
CARTINO PC est disponible sur le site du CETE Méditerranée.

La notice Fichier:PriseenMainCARTINOv7.pdf, au format pdf n'est plus mise à jour depuis mars 2013. Elle a été traduite dans ce format wiki. Seule cette page est mise à jour aujourd'hui.
PriseenMainCARTINOv7.jpgCARTINO PC.jpgCARTINO PC Reprise.jpg

Sommaire

[modifier] Pré-requis

  • Cf. wikhydro : préalables pour l'utilisation de Qgis-Grass sur le Lidar

http://www.wikhydro.org/index.php/Préalables_pour_l'utilisation_de_Qgis-GRASS_sur_le_LIDAR

  • Guide prise en main Outils DICARTO et utilisation de cet outil avant CARTINO.

L'utilisateur doit avoir, pour l'installation des différents composants, les droits d'administrateurs sur son PC et doit avoir accès à internet. Ensuite pour toute la durée des calculs, les droits d'écriture doivent être donnés sur les répertoires C:\Program Files (x86)\Quantum Gis\bin ou C:\Program Files \Quantum Gis\bin. De même que pour l'outil DICARTO, CARTINO a été développé sous Matlab. Son utilisation sur un PC classique nécessite le programme MCRInstaller.exe (à récupérer par exemple ici : http://cete-aix.fr/dwld/NUNIEAU/MCRInstaller.exe)

  • Créer un répertoire de travail « CARTINO_test1 » sur la racine du disque « C » ou « D ».

[modifier] Installation de Flutor

Concernant le logiciel Flutor, il faut le compiler pour qu'il fonctionne sur votre poste. Pour cela, il faut tout d'abord installer le logiciel MinGW (dossier fourni). A l'installation, cocher l'ensemble des compilateurs (C++ compiler, fortran compiler et objC compiler ). Puis, il faut accéder aux variables d'environnement du PC (sur win 98 : panneaux de configuration/paramètres systèmes/onglet avancé/variables d'environnement. Sur win7 : Panneau de configuration\Système\paramètres systeme avancés\variables d'environnement) et modifier la variable Path, en rajoutant le chemin d'accès à l'exécutable MinGW (;C:\MinGW\bin;). Redémarrer le PC. Enfin, ouvrir une boite dos (Demarrer/tous les programmes/Accessoires/invite de commande) , se positionner sur le dossier Plani (commande cd « chemin du répertoire de travail») et taper la commande écrite dans le fichier gfortran.txt du répertoire plani. Le nom de la version est : flutor_cartino_version_20130211_LMnmax.f90 La commande à taper est: gfortran -o flutor_cartino.exe flutor_cartino_version_20130211_LMnmax.f90. Attention, si la commande gfortran n'est pas reconnue c'est que le programme MinGW n'a pas été installé correctement, donc recommencer l'opération si besoin.

[modifier] Vérification de Flutor et Mascaret

Pour vérifier que Flutor et Mascaret fonctionnent, vous devez copier dans le dossier VERIF_FLUTOR_MASCARET votre version de mascaret et de flutor. Ensuite, cliquez sur mascaret.exe et flutor.exe, un fichier Result1.txt (flutor) apparaît ainsi que un fichier avec extension opt pour mascaret.

[modifier] Préalable au lancement de Cartino

L'utilisateur a besoin de 4 types de données d'entrée, à placer dans le répertoire de travail :

  • fichier vecteur d'EXZECO extrait sur la zone (format *.ClMD.shp)
  • fichier raster d'EXZECO extrait sur la zone (format *sol.asc)
  • dalles raster de MNT
  • fichier raster de débits (Shyreg ou autre)

Il faut vérifier la résolution suite aux exports ou retraitement , pour cela, aller dans les outils GRASS, parcourir et regardez les résolutions. Si la valeur de la résolution ou des bornes est mauvaises, les résultats seront erronés.

Il a également besoin de 4 fichiers ou dossiers pour lancer le logiciel :

  • CartinoV?.exe (la livraison est parfois faite avec une modification de l'extension pour passer les antivirus, il faut parfois le renommer avec le même nom que le fichier ctf avec l'extension exe)
  • CartinoV?.ctf
  • CARTINO_parametres.txt
  • Dossier Plani contenant les codes de calcul Mascaret, Flutor et Hec-Ras

L'ensemble de ces fichiers doivent être placés dans le répertoire des données du projet.

Attention, ne pas mettre de caractères spéciaux ou de vides dans les noms des fichiers shapefile ou GRASS et commencer avec une lettre et pas un chiffre.

[modifier] Création du secteur Grass et importation des données

[modifier] Données EXZECO

  • Ouvrir Qgis et ouvrir le fichier exzeco vecteur (raster classé) : fichier Extest1_CIMD.shp
  • Dans un premier temps, il convient de garder les zones vecteurs d'exzeco correspondantes aux cours d'eau à cartographier pour la DI, il faut donc supprimer tous les affluents qu'on ne veut pas traiter afin d'avoir un axe de calcul uniquement dans le cours d'eau principal et moins de nettoyage manuel par la suite. Pour supprimer des éléments, basculer en mode édition dans Qgis, sélectionner les entités à supprimer et les supprimer, enfin sauvegarder la couche (pour couper un polygone, utiliser l'outil « remodeler les entités »).
  • Créer un nouveau jeu de donnée Grass sur ce secteur : C:\GrassData\test1\test1
  • Ouvrir le fichier raster exzeco dans Qgis (*sol.asc)
  • Importer le fichier raster exzeco dans Grass (fonction r.in.gdal.qgis si le fichier est visible dans la fenêtre Qgis sinon r.in.gdal) : fichier RExtest1_G

Attention: cocher ne pas tenir compte de la projection dans les options avancées.

[modifier] Données MNT

  • Dans Qgis, construire un raster virtuel (onglet raster, divers, construire 1 raster virtuel) pour toutes les tables MNT : fichier MNTtest1.vrt (Il peut être utile au préalable ou dans GRASS de réduire la résolution du MNT (g.region et r.resample), travailler avec un raster de 1m sur des dizaines de km engendre d'énormes fichiers (2Go/30km))
  • Importer ce fichier dans Grass (fonction r.in.gdal.qgis ou r.in.gdal) : fichier MNTtest1_G

Attention: cocher ne pas tenir compte de la projection dans les options avancées.

  • Éventuellement modifier le MNT pour changer sa résolution (voir chapitre 5.1).

[modifier] Données de type Débit

1 Rmapcalc.jpg

Importer les différents fichiers de débit dans Grass (fonction r.in.gdal.qgis ou r.in.gdal) et les renommer : fichiers Q0010 à Q1000. Attention: cocher ne pas tenir compte de la projection dans les options avancées.

Si l'utilisateur n'est pas satisfait des débits SHYREG (qui semblent parfois un peu faibles) ou s'il souhaite travailler avec des débits différents, il a la possibilité de créer de nouvelles données de type débit (ou de modifier des données Shyreg). En effet, l'utilisateur a intérêt à lancer les calculs automatiques avec un débit maximaliste, ainsi cela évitera d'élargir manuellement tous les profils, dans les cas où on s’aperçoit que le débit n'est pas assez fort. Dans le cas de Shyreg, Il est donc conseillé de créer un débit supérieur deux ou trois fois à celui du débit de SHYREG pour Q1000. Procédure : Ouvrir la région Grass concernée

  • Utiliser la boîte à outil Grass
  • Sélectionner le fichier Q1000, puis Clic Droit et « Ajouter la carte sélectionné à la carte »
  • Cliquer sur 0 Region.jpg afin d'éditer la bonne région avec la bonne résolution.
  • Aller dans la « liste des modules » et chercher la fonction « r.mapcalc ». L'objectif est d'avoir, dans le cas ci-dessous, un débit trois fois supérieur à celui de Shyreg. Attention car dans le nom du rendu, il est préférable d'écrire Q1000_3 plutôt que Q1000X3.

[modifier] Renseignement du fichier de paramètres

MODIFICATION DE LA STRUCTURE. L'ANCIENNE FONCTIONNE MAIS VOUS POUVEZ GERER D'AUTRES PARAMETRES
Le fichier CARTINO_parametres.txt doit être renseigné avant le lancement de l'application CARTINO, ce fichier texte se compose des éléments suivants (les 4 premières lignes sont équivalentes à l'application dicarto):
Les valeurs en rouge doivent être modifiées.
Les valeurs en vert sont des options.


%---- Création du Listing (1er chiffre) et Suppression des fichiers temporaires GRASS (2ème chiffre) - 0 Non 1 Oui
0 1 1
%---- Dossier d'installation de Qgis-GRASS
C:\Program Files (x86)\Quantum Gis\bin
%---- Secteur GRASS - Chemin du répertoire avec les données au format GRASS
C:\GRASSDATA\CARTINO\Exemple
%---- Chemin du répertoire des données au format standard (shape, asc...),Attention, vous ne pouvez pas travailler sur C: ou D: directement!
D:\FormationQgisGRASS
%---- CODE: Code Flutor et Mascaret O Oui 1 Non => CODE - Etape 4
1 1
%---- TypBief: Type de Gestion du réseau (Paramètres non actif) et longueur minimale de la partie aval- Etape 4
1 500
%---- nomMNTasc: Nom du Raster MNT dans GRASS =>nomMNTasc - Etape 4
MNT1m
%---- nomSD: Nom du Raster EXZECO dans GRASS (utile pour connaitre les parties amont-aval avec la surface drainée=> fichierExze - Etape 3
SDEXZ
%---- nomQasc: Nombre et Noms du (des) Raster(s) Debit(s) à mettre du + fort au + faible =>nomQasc - Etape 4
5
Q1000
Q0300
Q0100
Q0030
Q0010
%---- nomEnvLisse: Nom de résultat pour l'enveloppe lissée et nom de l'enveloppe utile aux calculs de distance squelette-enveloppe => nomEnvLisse - Etape 2 et 4
1
EnveloppeLisse
%---- PasCoupExz: [Resolution du resultat EXZECO => PasExzeco] [Distance de coupure pour le squelette => distcoupsque] - Etape 2
25 5
%---- Reso: Resolution du raster MNT et Raster DEBIT=> Reso - Etape 4
1 50
%---- DistCoup: Distance de coupure du profil en travers pour le MNT et Debit=> coup - Etape 4
2 10
%---- PasTri: Pas de Pré-calcul des distance Axe de calcul - Enveloppe =>PasTri - Etape 4
75
%---- Strickler: Valeur des Strickler Lit mineur et majeur (le majeur n'est pas encore intégré) =>Strickler - Etape 4
20 15
%---- DPNC: Distance entre profils non constantes: [Ratio largeur de l'enveloppe - Longueur entre profil => factlarg] [Distance minimale entre 2 traces => distmin] [Distance maximale entre 2 traces => distmax]
1 50 500
0.33 10 500
%---- ValFiltre: [lissage 1ères largeurs des profils en travers => nblisse] [lissage des angles des profils => nbfiltre] [Lissage des altitudes des profils en travers => nlisse]
5 3 1
%---- Elarg: Ration, Valeur Minimum, Valeur maximum d'élargissement et largeur max du profil (en km paramètres non actif) pour l'enveloppe Initiale => PourcElargissement ValeurElargissement ValeurMaxElargissement%---- Ration, Valeur Minimum, Valeur maximum d'élargissement et largeur max du profil (en km paramètres non actif) suite aux calculs => PourcElargissement ValeurElargissement ValeurMaxElargissement
0.33 5 50 10000
0.2 10 150 10000
%---- Nombre de tour maximum pour agrandir les profils - Nombre de tours maximum pour lisser les profils à chaque agrandissement
15 20


[modifier] Dessin d'explications des paramètres

2 Dessin expli parametres.jpg

L'image ci-dessous montre l'enveloppe lissée, le bief et le point montre la discrétisation le long du bief pour faire les premiers profils (valeur « PasTri ») du fichier de paramètre soit 24 m dans notre cas. Pour chaque point, la distance minimale entre ce point et le bord est calculée => DistBord. DistBord permet ensuite de calculer la DPNC: Distance entre profils non constantes avec la formule suivante: DPNC=min(max(10,0.33*2*DistBord),1000). Un profil en travers est positionné tous les 0,33 de la largeur de l'enveloppe initiale avec des limites basses et hautes de 10 et 1000 m

La première largeur du profil en travers fait l'objet d'un élargissement avec la formule: Larg1= min(DistBord+min(max(0,33*DistBord,25),100),25000) et chaque coté du profil en travers est égale à cette largeur: LargRG=Larg1 et LargRD=Larg1 La première largeur est agrandie de 0,33 avec une limite minimale à 25 et maximale à 100. La largeur finale étant au maximum de 25000m.

Après chaque tour, si un profil en travers n'est pas assez large, il est augmenté avec la formule suivante uniquement en RG, RD ou les deux au besoin avec la formule pour la rive gauche: LargRG= 2*min(LargRG+min(max(0,33*DistBord,10),150),25000).

Une modification importante a eu lieu à partir de la version 3 de CARTINOPCv3 disponible sur internet, à savoir:

  • Un premier calcul avec des pas non constant et des traces non constantes avec flutor pour améliorer la solution initiale mais sans chercher à éviter les croisements entre profils
  • Un deuxième calcul qui reprend la largeur obtenu à la fin de premier et qui recalcule les positions des profils à pas non constants et à largeurs non constantes.

Ces deux passages permettent un gain de qualité, moins de croisements et aussi un gain en temps.

[modifier] Explication du paramètre CODE

Le paramètre CODE permet de faire des calculs avec Flutor et/ou Mascaret et de faire les cartographies avec Flutor et/ou Mascaret

CODE Flutor-Mascaret Condition aval Calcul Cartographie
1 1 Caval fournie par Flutor Flutor et Mascaret Flutor si convergence sinon Mascaret
2 1 Caval fournie par Flutor Flutor et Mascaret Mascaret si convergence sinon Flutor
1 2 Caval fournie par Flutor Flutor, pas de calcul Mascaret, juste les fichiers d'entrée Flutor si convergence
1 0 Caval fournie par Flutor Flutor Flutor si convergence
1 0 Caval fournie par un calcul Strickler peu fiable Mascaret Mascaret si convergence

ATTENTION Les résultats sont parfois instables sur de très fortes pentes, il est toujours préférable de vérifier les lignes d'eau provenant des 2 calculs.

De manière générale le calcul Mascaret est à privilégier pour la cartographie : options 2 1

Si l'utilisateur veut également lancer un calcul HEC-RAS, il faut rajouter un troisième chiffre 1 sur cette ligne.

[modifier] Lancement de Cartino

CARTINO PC.jpg

Ouvrir QGIS, et ouvrir le secteur d'études GRASS et lancer les outils GRASS dans Qgis. Lancer la console Shell des outils GRASS de Qgis (Attention pas une console DOS basique) et aller dans le répertoire de travail pour lancer la commande cartinov?

[modifier] Création des îles pour vérification manuelle

L'objet de cette partie est d'identifier et de traiter les îles de la pré-enveloppe de zones inondables issue des résultats d'Exzeco. En effet, une fois la pré-enveloppe créée à partir du fichier exzeco, on s’aperçoit que ce polygone peut avoir des trous (qu'on peut appeler des îles). Dans ce cas, quand le squelette va être tracé (au milieu de la pré-enveloppe) il va être divisé en deux. Or le calcul 1D que nous allons réaliser n'admet pas de zones de divergences. Donc il y a deux solutions: soit on supprime l'ile (c'est à dire qu'on va la fusionner avec le reste de la pré-enveloppe) soit on conserve l'ile et dans ce cas là il faudra ultérieurement venir modifier le squelette pour conserver uniquement une branche d'écoulement.

[modifier] Principes et explication du code

voir http://www.wikhydro.org/index.php/Cartino_PC

[modifier] Démarche à suivre

lancer le menu à partir de la boite de dialogue Cartino choisir l’extension .shp sélectionner le fichier vecteur exzeco (raster classé) : Extest1_CIMD.shp

Résultat de l'opération: Cartino crée deux fichiers shape:

  • Extest1_ClMDFus.shp : ce fichier est un ensemble de polygones représentant le contour du fichier exzeco.
  • Extest1_ClMDMasq.shp : ce fichier est un ensemble de polygones représentant le négatif du polygone précédent, c'est à dire l'ensemble des îles (c'est à dire des zones non inondables).

L'utilisateur va devoir passer par une étape manuelle, pour modifier le fichier Masq (Extest1_ClMDMasq.shp), de manière à choisir les entités qu'ils souhaitent fusionner à la pré-enveloppe contenue dans le fichier Extest1_ClMDFus.shp. Les entités qui seront supprimées de Extest1_ClMDMasq.shp ne seront donc pas fusionnées à l'enveloppe et le cours d'eau contournera ces îles) :

  • Ouvrir le fichier Masq sous Qgis
3 Masque.jpg

Ce fichier Masq est un ensemble de polygones (en bleu) représentant les zones non inondables, dont on va tout d'abord supprimer le plus gros élément, ici en bleu, grâce aux outils de numérisation de QGIS. Il ne reste plus que des iles:

4 Iles.jpg

Parmi les polygones restant, nous souhaitons que la modélisation contourne l'îlot cerclé de rouge, donc nous allons supprimer ce polygone.

Dans le fichier Masq ainsi modifié, il ne reste donc que les îlots que nous allons fusionner avec le reste du fichier pré-enveloppe. Enregistrer ce fichier sans modifier son nom pour suivre l'exemple.

Pour plus de sécurité, vous pouvez au départ enregistrer votre masque sous un autre nom afin de ne pas avoir à recommencer cette étape en cas de fausses manipulations.

[modifier] Fusion des îles, …, création du squelette

L'objet de cette partie est la fusion de la pré-enveloppe avec les îles retenues, la squelétisation et enfin le nettoyage pour la création de tronçons de cours d'eau.

[modifier] Principes et explication du code

voir http://www.wikhydro.org/index.php/Cartino_PC

[modifier] Démarche à suivre

Lancer le menu à partir de la boite de dialogue Cartino Choisir l'extension *.shp Sélectionner les deux fichiers Fus et Masq créés précédemment Si le fichier Masq est vide (aucune île) ne sélectionner que le fichier Fus


Résultat de l'opération: Cartino crée l'enveloppe du modèle (fichier EnveloppeLisse.shp) et le squelette après plusieurs itérations (fichiers Squelette03.shp à Squelette09.shp dans notre cas)

  • ouvrir le fichier squelette optimisé squelette09.shp sous Qgis (le squelette nettoyé correspond toujours au fichier ayant le plus chiffre et créé en dernier)
5 Limites SHYREG.jpg

  • nettoyer ce fichier avec les outils de numérisation 6 Coup.jpget numérisation avancée (activer le mode édition, puis commande séparer les entités : clic gauche pour débuter la ligne de séparation, clic droit pour terminer. Ensuite sélectionner l'entité à supprimer, puis supprimer)
    • supprimer les affluents qu'on ne veut pas traiter
    • supprimer les sections qui dépassent l'enveloppe MNT (sinon vous aurez un message de bug "Boom, l'interpolation de l'exutoire a échoué" )
    • supprimer les sections qui dépassent l'enveloppe des débits shyreg (entouré en rouge). En effet, une fois les profils en travers tracés sur le squelette, ils devront obligatoirement intercepter une donnée de débit, le squelette et shyreg peuvent donc être décalés latéralement.
7 Gestion Iles.jpg

    • choisir une seule branche de modélisation quand on a un îlot (la branche du bas dans notre cas)

On obtient ainsi un squelette nettoyé. Sauvegarder les changements en sortant du mode édition. Pour plus de sécurité, vous pouvez au départ enregistrer votre squelette sous un autre nom afin de ne pas avoir à recommencer l'étape de squelettisation en cas de fausses manipulations.

[modifier] Création du réseau des biefs de calculs

L'objet de cette partie est la création des différents réseaux de biefs pour chaque bassin versant, à partir du squelette.

[modifier] Principes et explication du code

voir http://www.wikhydro.org/index.php/Cartino_PC

[modifier] Démarche à suivre

lancer le menu à partir de la boite de dialogue Cartino choisir l’extension *.shp choisir le fichier squelette propre squelette*.shp choisir le facteur de multiplication des tampons au niveau des confluences. En effet, à chaque confluence, on cherche à supprimer les interpolations hasardeuses d'un affluent sur le cours d'eau principal qui peut conduire à créer du « sur-débit » sur l'affluent. Un facteur 1 était mis par défaut, 2 à 3 sont plus appropriés pour des petits cours d'eau. Par contre, le débit est quand même gardé sur le profil en travers qui passe exactement sur la confluence.

Résultat de l'opération: Cartino crée le réseau de bief pour chaque bassin versant : fichier Bv01.shp Dans le cas présent nous avons 2 biefs.

[modifier] Calcul hydraulique

L'objet de cette partie est la préparation des fichiers de calcul, le lancement des calculs et la création des fichiers de post-traitement.

[modifier] Principes et explication du code

voir http://www.wikhydro.org/index.php/Cartino_PC

[modifier] Démarche à suivre

  • lancer le menu à partir de la boite de dialogue Cartino
  • donner éventuellement une valeur de surcote aval (sinon mettre la valeur -99 pour une condition limite aval automatique). Flutor gère la condition limite de telle sorte que si votre valeur est trop faible, elle sera augmentée jusqu'à être hydrauliquement acceptable.
  • choisir l’extension *.shp
  • choisir le ou les bassins versant à calculer

Résultat : une fois les calculs terminés, Cartino crée pour chaque BV un répertoire propre, dans lequel seront stockés pour chaque bief, les résultats suivants : les fichiers résultats de calcul sur les traces de profils en travers : Bv01_bief001QQP0010T12TNC58Ligne3D.shp les fichiers résultats de calcul après post-traitement sur les hauteurs d'eau si les profils en travers ont été assez agrandis: Bv01_bief001QQP0010T12TNC02CharDICARTO.asc des sous-répertoires regroupant les fichiers de calcul et résultat Mascaret/Flutor: Bv01_bief001QQP0010T12TNCok L'extension ok indique que le résultat a permis de réaliser des cartes de charges spécifiques, ce qui ne veut pas forcément dire que le calcul est finalisé S'il n'y a pas ok, cela veut dire que les calculs ont fonctionné mais il n'a pas été possible de faire le croisement hauteur d'eau et MNT car trop de profils se croisent.

NB1: Si l'opérateur a fait une erreur dans une des étapes et que Cartino plante , il est recommandé de supprimer tous les nouveaux fichiers créés dans le monde Grass.

NB2: Il est possible de vérifier les résultats de calculs en cours afin de voir s'il n'y a pas d'énormes anomalies (limite du MNT, squelette, proximité de la mer, ….), une fois qu'un calcul intermédiaire est terminé.

Les variables hydrauliques calculées dans les fichier résultats de profils sont les suivantes:

- pas de temps de calcul : Pas_de_tem
- numéro de Bassin Versant : NBasVers
- numéro du bief : NBief
- numéro du profil : NSection
- distance cumulée : Distance_m
- cote de fond : ZREF_m
- cote de l'eau : Z_m
- débit en lit mineur : QMIN_m3_s
- débit en lit majeur : QMAJ_m3_s
- coefficient de frottement du lit mineur : KMIN_m1_3_
- coefficient de frottement du lit mineur : KMAJ_m1_3_
- section mouillée lit mineur: S1_m2
- section mouillée lit majeur : S2_m2
- nombre de Froude : FR
- vitesse en lit mineur : VMIN_m_s
- hauteur d'eau : Y_m
- charge hydraulique : CHAR_m
- largeur au miroir : BMAX_m
- débit maximal au cours de la simulation : QMAX_m3_s
- débit total : Q_m3_s

[modifier] Gestion des axes réalisés

L'objet de ce menu est de permettre à l’utilisateur de visualiser les résultats directement dans Qgis, à l'aide de hotlinks:

  • lancer le menu à partir de la boite de dialogue Cartino
  • choisir l’extension *.txt
  • choisir les fichiers BV01.txt, BV02.txt...
  • Résultat de l'opération: le programme va faire un assemblage des résultats de cartino sous forme shape, raster et image et va créer l'ensemble des hotlinks dans le fichier BVBiefCalculs.shp

Sous Qgis :

  • ouvrir le fichier BVBiefCalculs.shp (rien ne s'affiche car il faut modifier les propriétés d'affichage)
  • Ouvrir les propriétés de la couche (double clic gauche)
  • dans l'onglet style effacer tout, puis classer, puis OK
  • On peut alors visualiser le bief de calcul, sur lequel on peut lancer les hotlinks pour afficher les résultats de calculs: 8 Actions Qgis.jpg
    • fichiers *.shp : fichiers de profils en travers
    • fichiers *.asc : fichiers ascii de résultats de hauteurs d'eau
    • fichiers *.png : fichiers images des résultats de calculs (profils en long des charges)

Ce fichier BVBiefCalculs.shp est la base pour les reprises de calculs dans l'étape suivante (Croisement Axes – Profil) A noter que chaque axe est allongé de 1m amont et aval pour permettre une interpolation des sections qui par des jeux d'arrondi peuvent ne pas être intersectées .


Dans un premier temps, l'utilisateur peut effectuer certains types de vérifications :

  • continuité zone inondable
  • profils qui ne se croisent pas
  • débit correct
  • comparaison avec ZI existantes
  • problème de profil qui sautent (capacité en lit mineur ok pour le 1er profil, insuffisante pour le second...)

[modifier] Reprise des calculs

La Boite de dialogue indique les différentes étapes nécessaires à la reprise de calculs.

CARTINO PC Reprise.jpg

Ce menu doit aussi pouvoir s'utiliser sans avoir réalisé les calculs automatiques cartino, il suffit de se créer des axes et des profils avec des champs NBasVers et NBief renseignés. Par contre, il nécessite d'avoir le fichier cartino_parametres.txt.

Pour faire des reprises de calcul, il est préférable de créer un nouveau dossier « reprise » dans lequel on va placer:

  • les 3 fichiers cartino (cartino.exe, cartino.ctf, cartino_parametres.txt)
  • le dossier plani
  • les fichiers shape commençant par « BVBiefCalculs » qui correspondent aux axes de modélisation
  • les fichiers vecteurs de profils que l'on va modifier.

Dans le fichier cartino_parametres.txt, il faut changer le chemin du répertoire puisqu'on a créé un nouveau répertoire « reprise »

Cette première étape a pour objet de modifier le positionnement des profils en travers. Les champs attributaires qu'il faudra nécessairement renseigner dans les fichiers d'axe et sections sont : NBasVers, NBief

Il existe deux possibilités pour les reprises de calcul:

  • soit on travaille bief par bief, cela signifie que si on modifie un profil sur une zone commune à plusieurs biefs , il faudra modifier les fichiers de profils de chaque bief. Puis on peut lancer le calcul bief par bief
  • soit on décide de fusionner les fichiers shp avant modifications des profils en supprimant les parties communes des biefs, dans ce cas là on aura plus qu'un fichier, on pourra alors modifier les profils en travers une seule fois.

Par exemple:

  • la première option est adaptée à un travail avec beaucoup de biefs qui se croisent sur un même bassin versant (Cas d'un PPRi avec un travail à partir de 20ha de BV)
  • la deuxième est adaptée dans le cadre d'un TRI avec 6 cours d'eau de 6 bassins versants qui sont donc sans chevauchement.

La largeur des profils en travers peut être dépendante des débits. Parfois, on veut limiter l'extension des profils pour Q10 pour ne pas aller inonder des chenaux sans relation avec le lit mineur. Dans ce cas, le travail sur les sections de calcul se fera sur chaque gamme de débit. Pour qu'il n'y ait pas de confusion dans les noms de fichiers, il est recommandé dans ce cas de créer un répertoire reprise par modèle (repriseQ10, repriseQ100, etc...)

En général, les profils en travers faits sur les débits les plus forts peuvent convenir pour les débits inférieurs, ce qui limite les reprises géométriques manuelles.


Avant de lancer les boites de dialogue, le travail de modification doit être réalisé sous Qgis sur les fichiers shape de profils en travers (ex: Bv01_Bi001QQ1000T06TNC28Ligne3D.shp):

  • On peut supprimer des profils, en ajouter, les modifier.
  • Attention à ne pas vouloir essayer de compléter un profil sur un bord, il est conseillé de le supprimer et de le refaire.
  • Il est possible de réaliser des lignes brisées
  • Attention de ne pas tracer de profils en travers en dehors de l'emprise MNT
  • Supprimer impérativement tout croisement de profil
  • Il faut entrer une donnée attributaire à chaque nouveau profil comprenant le n° de bassin versant NBasVers et le n° de Bief NBief
  • Entrer également une valeur dans la colonne Nsection (effectuer la commande rownum sur la colonne après avoir effectué toutes les modifications).

Lors de cette étape, l'utilisateur peut également en profiter pour retravailler l'axe principal d'écoulement (fichier BVBiefCalculs.shp) afin de réaliser une meilleure modélisation. Par contre il doit y avoir une seule polyligne pour chaque bief.

9 ExReprise1.jpg

Dans notre cas on choisit de modifier les profils 39, 40, 41 on rajoute un profil en amont du profil 41 :

9 ExReprise2.jpg

Après avoir enregistré le fichier pour prendre en compte les modifications, deux étapes de vérification:

[modifier] Vérification du sens des axes

L'objet de cette étape est de vérifier que l' axe d'écoulement soit bien orienté de l'amont vers l'aval. Dans le cas contraire, l'utilisateur a la possibilité d'inverser le sens de l'axe.

Ce menu n'est pas utile quand l'utilisateur a récupéré l'axe BVBiefCalculsAuto.shp qui est déjà correctement orienté après les calculs automatiques.

A partir de cette étape, il est possible de fusionner tous les fichiers shp en un seul (Qgis/ menu vecteur/Outil de gestion de données/fusionner les shapefiles en un seul) afin de lancer la commande qu'une seule fois.

[modifier] Vérification du non-croisement des sections

Cette deuxième étape permet de vérifier que le fichier de profils (Bv01_Bi001QQ10...Ligne3D.shp) a été correctement repris et qu'il n'y a plus de croisement de profils. Si cartino détecte des problèmes, il va créer un fichier BV1Bief1ProblemeSections.shp. Il sera alors indispensable de modifier les sections qui posent problème (problème de type croisement ou superposition) dans le fichier de profils initial (Bv01_Bi001QQ1000...Ligne3D.shp).

[modifier] Croisement Axes – Profil

Dans cette étape, Cartino va croiser le fichier de profils avec l'axe d'écoulement afin de prendre en compte les modifications effectuées sur le fichier de profils. La boite de dialogue demande de choisir :

  • le fichier shape des profils modifiés,
  • le fichier des axes ( choisir l'axe BVBiefCalculsAuto.shp).

Résultat de l'opération : Cartino crée le fichier Bv01_Bi001_SectionsReprises.shp dans lequel tous les champs attributaires ont pu être recalculés (numéro de profil, abscisse...) pour la suite du calcul. Si on avait un fichier initial avec plusieurs bassins versants et plusieurs biefs, on aurait autant de fichiers que de couples uniques.

Avant l'étape suivante , il est a priori utile de le fusionner (Qgis/ menu vecteur/Outil de gestion de données/fusionner les shapefiles en un seul) pour lancer en une seule fois l'étape suivante.

[modifier] Interpolation Débit

L'objet de cette étape est de lancer l'interpolation sur les débits des rasters du fichiers cartino_parametres (dans notre cas Q0010, Q0030, Q0100, Q0300, Q1000) à partir du fichier précédemment créé (Bv01_Bi001_SectionsReprises.shp...), elle peut être faite sur un fichier avec plusieurs bassins versants et plusieurs biefs. Un fichier SectionsDebits-dateheure-.shp est créé avec les colonnes en plus correspondantes aux débits interpolés sur les matrices demandées dans le fichier Cartino (Q0010, Q0030, Q0100, Q0300, Q1000). Vous pouvez modifier les attributs de ces débits ou même faire une autre colonne de débit avec vos valeurs sans aucune interpolation, par exemple celles provenant de vos PPRi.

[modifier] Interpolation MNT (Raster)

L'interpolation sur le MNT se fait sur le fichier SectionsDebits-dateheure-.shp. Des fichiers Bv01_Bi001MNTtest1_G3dPtZ.csv sont ainsi créés, le nombre correspond aux couples uniques nombre de bassins versant et nombre de biefs

[modifier] Interpolation Lit Mineur-Majeur et coefficient de Strickler

Vous devez disposer dans votre répertoire d'un fichier Strickler.shp avec une géométrie polygonale contenant:

  • un champ LIT
    • 1 indique le lit mineur
    • 2 indique le lit majeur
  • un champ Strickler
    • avec une valeur de Strickler

Il faut sélectionner la liste des fichiers Bv01_Bi001MNTtest1_G3dPtZ.csv précédemment créés. Associés à chaque fichier, un fichier Bv01_Bi001MNTtest1_G3dPtZ_Strickler.csv sera créé.

[modifier] Calcul Hydraulique simple

Il suffit ensuite de lancer le calcul de débit sur le fichier SectionsDebits-dateheure-.shp, il faut juste donner les colonnes de débits que l'on veut lancer et donner éventuellement une cote NGF pour la condition limite aval. De même que pour les calculs automatiques, les fichiers mascaret, les fichiers résultats .shp et .asc sont créés (ils ne sont pas bien rangés dans les répertoires pour le moment, mais sont bien présents)

9 ExReprise3.jpg

Sur notre exemple, on peut s'apercevoir que la zone inondable a bien été prolongée au niveau du profil 39. Par ailleurs nous n'avons plus de zone inondable en rive droite (au niveau du profil 41)

[modifier] Covadisation

[modifier] Maximum de raster, vectorisation, suppression des petites zones, suppression des créneaux, mise au format COVADIS DI table ISO_HT...

En cours d'écriture Les résultats de la phase cartographie de la directive inondation doivent suivre un standard SIG fournit par la commission COVADIS et disponible sur http://www.cnig.gouv.fr/Front/index.php?RID=154. Cette page permet de créer à partir des fichiers résultats sous forme raster les formats demandés. ces éléments sont une combinaison des différents modules de DICARTO http://wikhydro.developpement-durable.gouv.fr/index.php/Notice_de_DICARTO:

  • Création du max des résultats
  • Vectorisation avec options
  • Suppression des petites aires
  • Suppression des créneaux
  • Création des champs COVADIS

Les étapes sont les suivantes:

  • choisir les suppressions des petites aires (de 1 à 100m² par pas de 1m par exemple (1:1:100) ou de 25 à 625 m² par pas de 25m2 (25:25:625)
  • donner la résolution du raster.

Cartino va alors créer un fichier covadisation.txt dans lequel seront automatiquement listés l'ensemble des fichiers asc présents dans le répertoire de données (un nom de fichier asc par ligne). Il convient alors de modifier ce fichier de la manière suivante sous Notepad++:

  • Pour chaque bassin versant, s'il existe plusieurs biefs (Bi001,Bi002...), il faut les regrouper en une seule ligne par la commande Bi*.
  • A chaque fin de ligne, on rajoute un ";" et un numéro qui correspond au choix du scénario (1 pour fréquent, 2 pour moyen, 4 pour extrême).
  • Au final on a donc une seule ligne par scénario pour chaque bassin versant.

Après avoir enregistré le fichier covadisation.txt, lancer la covadisation (bouton "FAIT")


Une fois cette étape terminée, il est possible de combiner les fichiers shape résultats sur tous les bassins versants.

  • Qgis Menu vecteur\Outils de Géométrie\Fusionner les shapefile en un seul

[modifier] Finalisation de la table ISO_HT pour mapinfo

Pour finaliser le travail dans mapinfo au besoin, il faut convertir les shapefile en mapinfo, le mieux est d'utiliser le traducteur universel. En utilisant la deuxième option, et en sélectionnant les fichiers tab que l'on veut modifier, il suffit de copier les lignes de commande mapbasic dans le menu Option\Fenetre mapbasic.

La routine Modif_rvb.mbx disponible dans le zip de CARTINO permet ensuite de mettre les couleur de manière automatique dans mapinfo.

Un travail est à faire au final pour combiner les données issues de cartino avec celles provenant des autres sources de données (PPRi...) afin d'avoir des identifiants uniques respectant le format COVADIS.

[modifier] Création de la table INONDABLE_S

Le troisième menu permet de créer la table N_[prefixTRI]_INONDABLE_S... à partir de fichier shape ou mapinfo. Il ne faut qu'un scenario par table.

[modifier] Améliorations prévues

La liste ci-dessous non exhaustive a pour objectif de fournir les améliorations que l'on s'engage à réaliser à très court terme:

  • Garder les Xmil et Ymil des axes afin de permettre d'avoir les profils en travers centrés sur l'axe comme dans la partie automatique

La liste ci-dessous non exhaustive a pour objectif de fournir les améliorations à réaliser à moyen terme:

  • Intégrer le moyen de gérer Topo LIDAR et Topo-Bathy terrestre, faire la liaison entre les profils en travers calculés dans DICARTO et ceux de CARTINO en particulier pour faire un calcul automatique avec un raster dégradé (5m) et faire une seule fois sur le LIDAR 1m dans l'étape de reprise.
  • Faire une vérification que les bords ne sont pas en eau pour évaluer les endroits où la cartographie doit être reprise à la main (remontée aval sur un affluent).
  • Faire une lecture et jointure de tous les résultats Mascaret qui sont réalisés avec des noms de section sous forme d'entier.

La liste ci-dessous non exhaustive a pour objectif de fournir les améliorations à réaliser à long terme:

  • Gestion automatique de profils qui se croisent

La liste ci-dessous non exhaustive a pour objectif de fournir les améliorations « ergonomiques » à réaliser pour limiter les problèmes de lecture des notices:

  • faire une vérification que le fichier de paramètres a des données correctes
  • améliorer les exports qml
  • gérer les biefs non plus par la plus grande distance mais par le plus fort débit ou surface drainée amont.

Si vous avez des idées d' améliorations atteignables rapidement, merci d'en faire part.

[modifier] La Vraie Vie

Quelques questions avec des propositions de réponses qui doivent faire l'objet d'une réflexion, ce n'est pas la réponse parfaite mais quelques idées.

[modifier] J'ai un très grand territoire et ça prend beaucoup de temps

Il peut être utile au préalable ou dans GRASS de réduire la résolution du MNT (g.region et r.resample) , travailler avec un raster de 1m sur des dizaines de km engendre d’énormes fichiers (2Go/30km). Si on veut le résultat final au pas du mètre, on peut ensuite faire la reprise sur le pas métrique ! La procédure a suivre pour passer un MNT au pas de 1m à ce même MNT au pas de 5 m est la suivante :

  • Sélectionner le MNT_1m et fixer sa région avec l'outil 0 Region.jpg
10 ExplorateurQgisGRASS.jpg

  • éditer sa région qui pourrait être la suivante :
10 Region1.jpg

  • Changer les paramètres pour définir une nouvelle région au pas de 5m. En faisant attention aux bornes du MNT qui doivent être changées. Le résultat est :
10 Region2.jpg

Rem : Comme la fonction « r.resample » utilisée se sert des plus proches voisins, il suffit de faire les opérations suivantes : Valeur Nord +2 m Valeur Sud -3 m Valeur Ouest -2 m Valeur Est +3 m


Enfin, dans la console Shell, lancer la commande « r.resample » pour créer le MNT_5m. Le fichier d'entrée est le MNT_1m et le fichier de sortie MNT_5m. Cette méthode s'approche de l'IGN où D'autres méthodes plus détaillées existent mais sont plutôt à utiliser dans l'autre sens: la commande r.resamp possède plusieurs déclinaisons qui permettent de choisir la méthode la plus appropriée.

  • r.resamp.stat permet de diminuer la résolution par aggrégat statistique
  • r.resamp.interp est généralement utilisé pour augmenter la résolution à l'aide des méthodes
  • * plus proche voisin (1 cellule, moins précis)
  • * bilinéaire (4 cellules, par défaut)
  • * bicubique (16 cellules, plus lent)

r.resamp.rst utilise la méthode spline (régularisé avec tension) pour un plus grand lissage.

[modifier] J'ai un semi de points et je veux faire un raster

Vous pouvez tester l'option r.in.xyz qui convertit des points en raster en ayant au préalable choisi votre région (bords à 0, résolution à 5m par exemple). Ensuite faire r.fillnulls pour combler les vides.

[modifier] J'ai des rasters avec des projections différentes

On peut être amené à faire travailler sur un territoire avec des MNT qui ne sont pas en RGF93. Pour cela : - On fait l'étape jusqu'à la squelettisation dans un monde GRASS RGF93 - On convertit le squelette et l'enveloppe lissée dans le nouveau référentiel, puis on ferme toutes les tables et ce monde Grass - On crée un nouveau secteur grass dans le nouveau référentiel, on importe l'enveloppe lissée. - On projette notre raster shyreg et exzeco dans cette nouvelle projection (méthode plus proche voisin et surtout pas d'interpolation bi-linéaire!) Avec la fonction r.proj de GRASS (à partir de shell) et en lisant les avertissement sur l'utilisation de région dans le manuel. Attention : au préalable, bien définir les paramètres de la région (bornes et résolution) avant de lancer la commande ! - On importe le fichier MNT dans le monde Grass - On reprend cartino à partir des biefs dans le nouveau référentiel, dans le nouveau secteur Grass

[modifier] Je veux créer ma matrice de débit avec des formules simples

Q s0.8.jpg

C'est possible de le faire en partant de la donnée EXZECO raster (Fichier FRD_TRI_MARSEILLE_sol.asc) qui est une surface drainée. Il faut ensuite appliquer la calculatrice raster (r.mapcalc) pour faire des calculs comme Q=A.S0^,8 avec S correspondant à EXZECO
Avant tout, il faut vérifier la région Grass et la résolution
Ouvrir le fichier "FRD_TRI_MARSEILLE_sol.asc" dans Q-Gis et l'importer dans Grass en le nommant "surface" par exemple.
Dans les outils Grass, faire "rafraichir" et "surface" apparait.
Double-clic sur "surface" qui apparait dans la fenêtre Q-Gis.
A partir de là, il suffit d'utiliser la calculatrice des outils Grass "r.mapcalc"

[modifier] Je n'ai pas forcement des débits finaux mais je veux avancer

L'hydrologie est un élément clé de la modélisation. Les discussions sur les débits à prendre peuvent être longues et des discussions peuvent apparaître sur SHYREG, sur des débits locaux... Rien ne vous empêche de travailler sur cartino pour faire votre modèle sur lequel il faudra juste faire une reprise avec le (les) débit(s) final(aux) à utiliser. Pour cela SHYREG est une très bonne base de départ mais attention si SHYREG sous-estime vous allez faire des profils en travers peut-être trop étroit et votre travail ne sera pas finalisé. N'ayez pas peut de prendre un peu de sécurité en majorant SHYREG 1000 ans . Multiplier par 1,5, 2? à voir localement pour avoir l'ensemble de l'emprise extrême.

[modifier] Je n'ai pas forcément le MNT de l'IGN mais je vais le recevoir

Comme pour l'hydrologie, vous pouvez tester sur le MNT25m et ensuite il faudra changer le MNT pour l'étape de reprise, ceci peut éventuellement fournir une solution pour aller plus vite.

[modifier] Je veux choisir plus d'option dans mon modèle hydraulique

Vous avez la possibilité d'ouvrir les fichiers .cas, .geo, et .loi dans l'interface Fudaa-Mascaret pour modifier les limites lit mineur – lit majeur, les frottements, les formules de débordement, de faire votre modèle en mode non permanent...

[modifier] Je souhaite lancer plusieurs fichiers CARTINO ou DICARTO mais ça plante

Pour les utilisateurs avertis:

  • Ne pas travailler ni dans le même monde grass, ni dans le même répertoire de données classiques
  • Dans votre 1er projet, vous ajoutez 1 en 3ème chiffre dans le fichier paramètres de DICARTO ou CARTINO-PC

%---- Création du Listing (1er chiffre) et Suppression des fichiers temporaires GRASS (2ème chiffre) - 0 Non 1 Oui 0 1 1

  • Dans votre 2ème projet, vous ajoutez 1 en 3ème chiffre dans le fichier paramètres de DICARTO ou CARTINO-PC

%---- Création du Listing (1er chiffre) et Suppression des fichiers temporaires GRASS (2ème chiffre) - 0 Non 1 Oui 0 1 2

  • … jusqu'à 8 maximum.

Si vous avez un PC avec 4 processeurs, utilisez en 3, si vous avez 2, ne le faites pas!!!

[modifier] Qu'est-ce que je fais quand j'ai fini CARTINO et que je veux des résultats au format vecteur

Il « suffit » d'utiliser les routines de DICARTO pour convertir en iso-valeurs. ATTENTION, cela peut prendre du temps, il est conseillé de valider les résultats raster avec les Maîtrises d'ouvrages, collectivités, parties prenantes et lorsque vous êtes d'accord de faire ce travail qui prend du temps.

[modifier] Que faire lorsqu'on est dans des plaines en toit ou les lits perchés

La routine automatique peut ne pas être pertinente sur des secteurs de plaines en toît qui finissent en delta ou estuaires. Il est proposé de lancer l'outil automatique sur les parties amont tant que la zone d'écoulement est plutôt en « V ». Lors des reprises, il faut rajouter l'axe et mettre des sections de calcul en plus de manière « experte ». Il faut aussi avoir en-tête que la modélisation 1D, de plus mis en œuvre de manière rapide a des limites.

[modifier] Quels débits prendre en premier

La logique voudrait de prendre les débits du plus faible au plus fort, par exemple Q10, Q30, Q100, Q300, Q1000. Mais la quasi-obligation de reprendre les sections de calcul poussent plutôt à travailler en automatique que sur le plus gros débit, même supérieur à Q1000 SHYREG si on doit relancer des calculs par la suite. Il faut donc gérer les nombre de tour maximum pour agrandir l'enveloppe initiale.

[modifier] Que faire en cas de plantages

Vérifier que vous avez bien lu la notice! Sinon, les outils peuvent malheureusement « planter », pour les remontées, merci de copier ce qu'il y a écrit dans la console SHELL (Clic droit dans la console, Sélectionner tout, CTRL C, aller dans un éditeur de texte CTRL V).
Vous pouvez aussi ouvrir l'explorateur des outils Qgis-Grass, si des vecteurs n'ont pas de + associés, cela veut dire que le fichier est vérolé. Plusieurs solutions:

  • Fermer Qgis et GRASS
  • supprimer ces fichiers dans votre monde GRASS
  • arrêtez des processus comme dbf.exe
  • redémarrez votre PC, parfois plusieurs fois en cas d'erreur Win 32.
  • faire un autre monde GRASS

[modifier] Comment modifier le MNT

Il peut être nécessaire de diminuer des valeurs de MNT (passage de ponts par exemple. Ceci ne doit pas être le cas a priori sur les données IGN 1m. Pour le faire, définir des polygones avec des valeurs soit d'altitudes, soit de soustraction. L'importer dans grass (v.in.ogr), le convertir en raster (v.to.rast) et ensuite faire un r.series par exemple avec l'option minimum entre votre MNT et le polygone converti en raster. Dans Qgis, l'extension Point Simpling Tool disponible sur le dépot QGIS Official Repository 2 permet aussi de faire ce travail.

[modifier] Méthode appliquée pour s’affranchir des deux premières étapes du modèle CARTINO PC

La procédure proposée ci-dessous permet à un utilisateur d’utiliser son réseau hydrographique et de s’affranchir totalement de la source de donnée EXZECO.

Prétraitement CARTINO

Lors de l’intégration d’un raster à surface drainée autre que celui d’EXZECO, il est possible d’être confronté à des problèmes de résolution N/S et E/W sous CARTINO. En imposant une emprise et une résolution connues au raster des surfaces drainées, on s’affranchit de cette erreur.

  • Pour commencer, il faut construire un vecteur correspondant à l’emprise du raster SHYREG : Vecteur / Outil de recherche / Créer un polygone à partir de l’emprise de la couche. L’emprise obtenue est utilisée afin d’y extraire les nœuds du vecteur : Vecteur / Outil de géométrie / Extraire les nœuds.
  • Parallèlement, à cette manipulation, il faut un réseau hydrographique répondant aux contraintes de CARTINO dont le faite que ses entités soient couvertes par le raster MNT et l’emprise nouvellement créée. Il faut extraire les nœuds de ce réseau : Vecteur / Outil de géométrie / Extraire les nœuds.
  • Les deux vecteurs de points, l’un représentant l’emprise de SHYREG et l’autre la position du réseau hydrographique, sont unis : Vecteur / Outil de géotraitement / Unir.
  • Puis, ce vecteur est modifié lors d’une phase d’édition. Un champ exutoire est ajouté, grâce à la calculatrice de champ, pour l’ensemble des entités avec la valeur « 0 ».
Calculatrice de champ.png
  • Manuellement, les exutoires seront déterminés. Cela demande à l’utilisateur de connaître l’orientation de son réseau hydrographique. Il faut modifier la valeur des deux derniers points de l’exutoire de chaque bassin versant en attribuant la valeur « 1 ».
Verification attributaire.png

Ex : s’il y a 6 bassins versants, douze lignes sont modifiées correspondant aux deux derniers points de chaque branche du réseau hydrographique.

  • Avant d’interpoler ce vecteur, il est important de noter les informations contenues dans le raster SHYREG tel que l’emprise de la couche, la taille des pixels et la dimension du raster. Ces information sont disponibles dans les propriétés de la couche au niveau de l’onglet « métadonnées ».

Ex : les valeurs ci-dessous permettent d’avoir un aperçu du type d’information à noter

Tableau des donnees a conserver.png
  • Suite à cette étape, il est possible de convertir le vecteur de points en raster à partir de l’information exutoire précédemment renseignée : Raster / Interpolation / Interpolation.
Extension interpolation.png

Modèle CARTINO

Suite à cette note, l’utilisateur de CARTINO peut commencer à partir de l’étape 3 « Création du réseau des biefs de calculs ».

Au préalable, il devra charger une enveloppe lisse sous GRASS par le biais de la fonction v.in.ogr pouvant correspondre à un buffer de son réseau hydrographique. Le nom du fichier sous GRASS doit correspondre à l’information contenue dans le fichier CARTINO_paramètre.txt ; par défaut « Enveloppelisse ».

L’ensemble de ses manipulations aboutissent à une substitution du raster EXZECO et à l’emploie d’un réseau hydrographique extérieur au modèle.



Le créateur de cet article est Mathieu Alquier
Note : d'autres personnes peuvent avoir contribué au contenu de cet article, [Consultez l'historique].

  • Pour d'autres articles de cet auteur, voir ici.
  • Pour un aperçu des contributions de cet auteur, voir ici.

Le créateur de cet article est Frédéric Pons
Note : d'autres personnes peuvent avoir contribué au contenu de cet article, [Consultez l'historique].

  • Pour d'autres articles de cet auteur, voir ici.
  • Pour un aperçu des contributions de cet auteur, voir ici.


Le créateur de cet article est Christophe Esposito
Note : d'autres personnes peuvent avoir contribué au contenu de cet article, [Consultez l'historique].

  • Pour d'autres articles de cet auteur, voir ici.
  • Pour un aperçu des contributions de cet auteur, voir ici.

Le créateur de cet article est Guillaume Le Bihan
Note : d'autres personnes peuvent avoir contribué au contenu de cet article, [Consultez l'historique].

  • Pour d'autres articles de cet auteur, voir ici.
  • Pour un aperçu des contributions de cet auteur, voir ici.
Outils personnels