From: Paul RASCLE Date: Fri, 7 Apr 2017 17:57:41 +0000 (+0200) Subject: doc Land Cover Map X-Git-Tag: Salome_8_3_Hydro_1_1rc1~47 X-Git-Url: http://git.salome-platform.org/gitweb/?a=commitdiff_plain;h=f8fdeb4e34afff51f3fa2f21036b1aa803419ec2;p=modules%2Fhydro.git doc Land Cover Map --- diff --git a/doc/salome/tutorial/_static/BottomFriction.png b/doc/salome/tutorial/_static/BottomFriction.png new file mode 100644 index 00000000..94481a91 Binary files /dev/null and b/doc/salome/tutorial/_static/BottomFriction.png differ diff --git a/doc/salome/tutorial/_static/casLandCoverMap.png b/doc/salome/tutorial/_static/casLandCoverMap.png new file mode 100644 index 00000000..a5ae5e08 Binary files /dev/null and b/doc/salome/tutorial/_static/casLandCoverMap.png differ diff --git a/doc/salome/tutorial/_static/importLandCoverMap.png b/doc/salome/tutorial/_static/importLandCoverMap.png new file mode 100644 index 00000000..936eb28d Binary files /dev/null and b/doc/salome/tutorial/_static/importLandCoverMap.png differ diff --git a/doc/salome/tutorial/_static/importLandCoverMap_2.png b/doc/salome/tutorial/_static/importLandCoverMap_2.png new file mode 100644 index 00000000..d5aba03c Binary files /dev/null and b/doc/salome/tutorial/_static/importLandCoverMap_2.png differ diff --git a/doc/salome/tutorial/_static/importLandCoverMap_3.png b/doc/salome/tutorial/_static/importLandCoverMap_3.png new file mode 100644 index 00000000..86a5aa6f Binary files /dev/null and b/doc/salome/tutorial/_static/importLandCoverMap_3.png differ diff --git a/doc/salome/tutorial/_static/stricklerTable_1.png b/doc/salome/tutorial/_static/stricklerTable_1.png new file mode 100644 index 00000000..6e472ae3 Binary files /dev/null and b/doc/salome/tutorial/_static/stricklerTable_1.png differ diff --git a/doc/salome/tutorial/interpolStrickler.py b/doc/salome/tutorial/interpolStrickler.py new file mode 100644 index 00000000..8bac9e61 --- /dev/null +++ b/doc/salome/tutorial/interpolStrickler.py @@ -0,0 +1,14 @@ +# -*- coding: utf-8 -*- + +# ===== Description du cas, a éditer ===== +#========================================= + +# --- nom du cas dans HYDRO +nomCas = 'garonne_1' + +# --- fichier med 2D(x,y) du cas, produit par SMESH +fichierMaillage = '/tmp/garonne_1F.med' + +from salome.hydrotools.interpolS import assignStrickler + +assignStrickler(nomCas, fichierMaillage, '/tmp/garonne_1F.med', 'BOTTOM_FRICTION') diff --git a/doc/salome/tutorial/landCoverMap.rst b/doc/salome/tutorial/landCoverMap.rst index b449d64a..44b72eee 100644 --- a/doc/salome/tutorial/landCoverMap.rst +++ b/doc/salome/tutorial/landCoverMap.rst @@ -20,7 +20,129 @@ Land Cover Map ######################################### -A compléter. + +.. |stricklerTable_1| image:: /_static/stricklerTable_1.png + :align: middle + +.. |importLandCoverMap| image:: /_static/importLandCoverMap.png + :align: middle + +.. |importLandCoverMap_2| image:: /_static/importLandCoverMap_2.png + :align: middle + +.. |importLandCoverMap_3| image:: /_static/importLandCoverMap_3.png + :align: middle + +.. |casLandCoverMap| image:: /_static/casLandCoverMap.png + :align: middle + +.. |BottomFriction| image:: /_static/BottomFriction.png + :align: middle + + +Import d'une Land Cover Map +=========================== +Nous allons reprendre l'exemple précedent "garonne_1" pour ajouter une carte des coefficients de Strickler. +La carte à importer peur être téléchargée et éditée préalablement dans qgis. Les *Corine Land Cover* couvrent généralement +un grand territoire et sont très détaillées. + +**Remarque :** Pour que ces cartes ne soient pas trop lourdes à manipuler dans SALOME HYDRO, +il est préférable de les découper dans qgis de façon à les restreindre autant que possible au domaine d'étude. +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 +(Une carte découpée est fournie avec ce tutoriel). La carte comprend une base de données qui donne, pour chaque zone, +le code de type de territoire associé, selon la nomenclature *Corine Land Cover*. + +Avant d'importer la carte, il faut disposer d'une table définissant les codes *Corine Land Cover*, la couleur de représentation +des types de zones, les coefficients de Strickler associés. + +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 +*Strickler table_1*. Il est possible de l'éditer avec le menu contextuel *Edit Strickler table*. +La plupart des champs sont éditables, +**notamment les coefficients de Strickler qu'il faudra de toutes façon adapter pour une étude réaliste**. + + |stricklerTable_1| + +Le champ *Attribute Name* doit correspondre au champ de codification des types de zones dans la base de donnée +associée au *Corine Land Cover* importé : *CODE_06* correspond aux cartes de 2006, *CODE_12* aux cartes de 2012. +Ce champ n'est pas éditable dans le dialogue ci-dessus. + +Dans l'exemple fourni, la carte date de 2006, il nous faut donc une autre table de Strickler. +Ces tables sont enregistrées dans des fichiers texte éditables. +Nous chargeons une table 2006 fournie dans le répertoire d'installation de l'application SALOME HYDRO +avec la commande *Import Strickler table* du menu contextuel de la rubrique *STRICKLER TABLES* :: + + /share/salome/resources/hydro/def_strickler_table_06.txt + +Il faut effacer la table précédente *Strickler table_1*, pour que la nouvelle table soit correctement utilisée +lors de l'import de la Land Cover Map (menu contextuel *Delete*). + +Nous pouvons maintenant charger la Land Cover Map. +Nous utilisons pour cela la commande *Import land cover map from file(s)* du menu contextuel de la rubrique +*LAND COVER MAPS* :: + + /bin/salome/test/HYDRO/HYDRO/CLC_decoupe.shp + +La carte comprend déjà un grand nombre de polygones noirs (plus de 600) dont on peut voir le contour en les selectionnant. +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, +sur les bords. Elles sont hors de la zone d'étude. + + |importLandCoverMap| + +En pratique, on sélectionne tous les polygones ( A dans la liste), avant d'appuyer sur le bouton *Next>*. +Il y a plusieurs attributs trouvés dans la base importée. Nous selectionnons l'attribut qui nous intéresse, +le type de zone : *CODE_06*, avant d'appuyer sur le bouton *Next>*. +Les codes de zones sont correctement associés à leur définition fournie dans la nouvelle table de Strickler. + + |importLandCoverMap_2| + +Nous appuyons sur le bouton *Finish*. Le traitement prend quelques dizaines de secondes. +La Land Cover Map ne s'affiche pas automatiquement : commande *show* du menu contextuel de l'objet *CLC_decoupe* +dans la rubrique *LAND COVER MAPS*. On peut voir le cas de calcul sous la carte, en transparence. + + |importLandCoverMap_3| + +Création, édition de Land Cover Map +=================================== + +Il est possible de créer de toutes pièces une Land Cover Map, à l'aide de polygones créés dans SALOME-HYDRO +ou importés. + +Le menu contextuel de la Land Cover Map offre des possibilités d'édition des zones : +ajouter, enlever, découper, regrouper, changer le type. + +Nous ne détaillons pas ces opérations ici. + +Utilisation d'une Land Cover Map dans le cas de calcul +====================================================== + +Pour exploiter la Land Cover Map, il faut éditer le cas de calcul. +Nous éditons le cas de calcul *garonne_1* avec la commande *edit calculation case* du menu contextuel du cas de calcul. +Il faut avancer jusqu'au panneau *Land cover map* avec les boutons *Next>* et sélectionner le Land Cover Map et la +table de Strickler + + |casLandCoverMap| + +Nous poursuivons avec les boutons *Next>* et *Finish* pour valider le cas. + + +Création d'un champ de STrickler aux noeuds du maillage +======================================================= + +Telemac exploite un champ des coefficients de Strickler aux noeuds du maillage. Ce champ a pour nom *BOTTOM_FRICTION*. +Ce champ est ajouté au fichier MED du maillage, comme le champ d'altitude aux noeuds. + +Il faut adapter manuellement le script ci-dessous : + +Il faut recopier le script ci-dessous et l'adapter en fonction du cas de calcul +et des noms de de fichiers en entrée et en sortie. +Il est possible d'utiliser le même fichier MED en entrée et en sortie. + +.. literalinclude:: interpolStrickler.py + :lines: 1- + +Nous pouvons voir le champ résultat avec le module MED, comme pour le champ d'altitude. + + |BottomFriction|   .. only:: html   diff --git a/doc/salome/tutorial/preliminaires.rst b/doc/salome/tutorial/preliminaires.rst index 5084a2c2..b0dec621 100644 --- a/doc/salome/tutorial/preliminaires.rst +++ b/doc/salome/tutorial/preliminaires.rst @@ -324,7 +324,18 @@ Il faut mettre ces objets dans le repère local avant des les importer. Tables de coefficients de Strickler, Land Cover maps ==================================================== -A compléter, la prise en compte des coefficients de strickler et des Land Cover Maps n'est pas finalisée dans la version 1.0 de mars 2016. +Il est possible définir une carte des coefficients de Strickler (frottements sur le fond) couvrant le domaine d'étude. + +Des cartes décrivant la nature des sols (Land cover Map) peuvent être récupérées sur différents sites. +Il est notamment possible de télécharger et d'éditer dans un outil de SIG (Systeme d'Information Géographique) +comme *qgis* les cartes "Corine Land Cover". +Ces cartes s'appuient sur une nomenclature standard des différents types de territoire. + +On définit en parallèle une table des coefficients de Strickler qui donne le coefficient pour chaque type défini dans la nomenclature. +Les coefficients de Strickler sont en général ajustés pour une étude donnée, pour recaler le modèle. + +Les Land Cover Map peuvent être importées depuis qgis ou créées "from scratch" dans SALOME-HYDRO. +Il est également possible d'éditer ces cartes dans SALOME-HYDRO : ajout, suppression, regroupement, modification de zones... Constitution d'un cas de calcul ================================