doc/salome/tutorial/landCoverMap.rst

   1 ..
   2    Copyright (C) 2015-2016 EDF
   3
   4    This file is part of SALOME HYDRO module.
   5
   6    SALOME HYDRO module is free software: you can redistribute it and/or modify
   7    it under the terms of the GNU General Public License as published by
   8    the Free Software Foundation, either version 3 of the License, or
   9    (at your option) any later version.
  10
  11    SALOME HYDRO module is distributed in the hope that it will be useful,
  12    but WITHOUT ANY WARRANTY; without even the implied warranty of
  13    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  14    GNU General Public License for more details.
  15
  16    You should have received a copy of the GNU General Public License
  17    along with SALOME HYDRO module.  If not, see <http://www.gnu.org/licenses/>.
  18
  19 #########################################
  20 Land Cover Map
  21 #########################################
  22
  23
  24 .. |stricklerTable_1| image:: /_static/stricklerTable_1.png
  25    :align: middle
  26
  27 .. |importLandCoverMap| image:: /_static/importLandCoverMap.png
  28    :align: middle
  29
  30 .. |importLandCoverMap_2| image:: /_static/importLandCoverMap_2.png
  31    :align: middle
  32
  33 .. |importLandCoverMap_3| image:: /_static/importLandCoverMap_3.png
  34    :align: middle
  35
  36 .. |casLandCoverMap| image:: /_static/casLandCoverMap.png
  37    :align: middle
  38
  39 .. |BottomFriction| image:: /_static/BottomFriction.png
  40    :align: middle
  41
  42
  43 Import d'une Land Cover Map
  44 ===========================
  45 Nous allons reprendre l'exemple précédent "garonne_1" pour ajouter une carte des coefficients de Strickler.
  46 La carte à importer peur être téléchargée et éditée préalablement dans qgis. Les *Corine Land Cover* couvrent généralement
  47 un grand territoire et sont très détaillées.
  48
  49 **Remarque :** Pour que ces cartes ne soient pas trop lourdes à manipuler dans SALOME HYDRO,
  50 il est préférable de les découper dans qgis de façon à les restreindre autant que possible au domaine d'étude.
  51 On définira un polygone de découpe englobant le domaine d'étude, en vérifiant que l'on ne perd pas de zone lors de la découpe
  52 (Une carte découpée est fournie avec ce tutoriel). La carte comprend une base de données qui donne, pour chaque zone,
  53 le code de type de territoire associé, selon la nomenclature *Corine Land Cover*.
  54
  55 Avant d'importer la carte, il faut disposer d'une table définissant les codes  *Corine Land Cover*, la couleur de représentation
  56 des types de zones, les coefficients de Strickler associés.
  57
  58 Une table par défaut est proposée dans SALOME-HYDRO, visible dans l'arbre d'étude, dans la rubrique *STRICKLER TABLES*, sous le nom
  59 *Strickler table_1*. Il est possible de l'éditer avec le menu contextuel *Edit Strickler table*.
  60 La plupart des champs sont éditables,
  61 **notamment les coefficients de Strickler qu'il faudra de toutes façon adapter pour une étude réaliste**.
  62
  63   |stricklerTable_1|
  64
  65 Le champ *Attribute Name* doit correspondre au champ de codification des types de zones dans la base de donnée
  66 associée au *Corine Land Cover* importé : *CODE_06* correspond aux cartes de 2006, *CODE_12* aux cartes de 2012.
  67 Ce champ n'est pas éditable dans le dialogue ci-dessus.
  68
  69 Dans l'exemple fourni, la carte date de 2006, il nous faut donc une autre table de Strickler.
  70 Ces tables sont enregistrées dans des fichiers texte éditables.
  71 Nous chargeons une table 2006 fournie dans le répertoire d'installation de l'application SALOME HYDRO
  72 avec la commande *Import Strickler table* du menu contextuel de la rubrique *STRICKLER TABLES* ::
  73
  74   <appli_xxx>/share/salome/resources/hydro/def_strickler_table_06.txt
  75
  76 Il faut effacer la table précédente *Strickler table_1*, pour que la nouvelle table soit correctement utilisée
  77 lors de l'import de la Land Cover Map (menu contextuel *Delete*).
  78
  79 Nous pouvons maintenant charger la Land Cover Map.
  80 Nous utilisons pour cela la commande *Import land cover map from file(s)* du menu contextuel de la rubrique
  81 *LAND COVER MAPS* ::
  82
  83   <appli_xxx>/bin/salome/test/HYDRO/HYDRO/CLC_decoupe.shp
  84
  85 La carte comprend déjà un grand nombre de polygones noirs (plus de 600) dont on peut voir le contour en les sélectionnant.
  86 La découpe déborde très largement de la zone d'étude. On voit les zones "oubliées" lors de la découpe dans qgis,
  87 sur les bords. Elles sont hors de la zone d'étude.
  88
  89   |importLandCoverMap|
  90
  91 En pratique, on sélectionne tous les polygones (<crtl> A dans la liste), avant d'appuyer sur le bouton *Next>*.
  92 Il y a plusieurs attributs trouvés dans la base importée. Nous sélectionnons l'attribut qui nous intéresse,
  93 le type de zone : *CODE_06*, avant d'appuyer sur le bouton *Next>*.
  94 Les codes de zones sont correctement associés à leur définition fournie dans la nouvelle table de Strickler.
  95
  96   |importLandCoverMap_2|
  97
  98 Nous appuyons sur le bouton *Finish*. Le traitement prend quelques dizaines de secondes.
  99 La Land Cover Map ne s'affiche pas automatiquement : commande *show* du menu contextuel de l'objet *CLC_decoupe*
 100 dans la rubrique *LAND COVER MAPS*. On peut voir le cas de calcul sous la carte, en transparence.
 101
 102   |importLandCoverMap_3|
 103
 104 Création, édition de Land Cover Map
 105 ===================================
 106
 107 Il est possible de créer de toutes pièces une Land Cover Map, à l'aide de polygones créés dans SALOME-HYDRO
 108 ou importés.
 109
 110 Le menu contextuel de la Land Cover Map offre des possibilités d'édition des zones :
 111 ajouter, enlever, découper, regrouper, changer le type.
 112
 113 Nous ne détaillons pas ces opérations ici.
 114
 115 Utilisation d'une Land Cover Map dans le cas de calcul
 116 ======================================================
 117
 118 Pour exploiter la Land Cover Map, il faut éditer le cas de calcul.
 119 Nous éditons le cas de calcul *garonne_1* avec la commande *edit calculation case* du menu contextuel du cas de calcul.
 120 Il faut avancer jusqu'au panneau *Land cover map* avec les boutons *Next>* et sélectionner le Land Cover Map et la
 121 table de Strickler
 122
 123   |casLandCoverMap|
 124
 125 Nous poursuivons avec les boutons *Next>* et *Finish* pour valider le cas.
 126
 127
 128 Création d'un champ de Strickler aux noeuds du maillage
 129 =======================================================
 130
 131 Telemac exploite un champ des coefficients de Strickler aux noeuds du maillage. Ce champ a pour nom *BOTTOM FRICTION*.
 132 Ce champ est ajouté au fichier MED du maillage, comme le champ d'altitude aux noeuds.
 133
 134 Il faut adapter manuellement le script ci-dessous :
 135
 136 Il faut recopier le script ci-dessous et l'adapter en fonction du cas de calcul
 137 et des noms de fichiers en entrée et en sortie.
 138 Il est possible d'utiliser le même fichier MED en entrée et en sortie.
 139
 140 .. literalinclude:: interpolStrickler.py
 141     :lines: 1-
 142
 143 Nous pouvons voir le champ résultat avec le module MED, comme pour le champ d'altitude.
 144
 145   |BottomFriction|
 146
 147 .. only:: html
 148
 149    :ref:`ref_casParticuliers`