Le portail francophone de la géomatique

Bonjour,

Je viens d'avoir le même soucis en suivant le protocole de l'agence de l'eau pour la définition des PGSZH.
Voici ma solution pour télécharger et utiliser les données d'IDPR initialement fournies en WMS par le BRGM.

1. Enregistrer sous la couche raster depuis le WMS avec l'emprise du projet
2. Le résultat est une couche .vrt avec un affichage aux couleurs initiales via 4 bandes (RBG + Alpha) Chaque couleur correspond à une classe d'IDPR (0-200 ; 201-400 etc.)
3. Calculatrice raster pour créer une couche sommant les valeurs des trois bandes RGB, ce qui donne une valeur unique à chaque classe initiale.

Classes d'IDPR    R    G    B    Somme RGB    Classe finales
0-200        194    82    61    337            1
201-400        217    117    41    375            2
401-600        237    168    19    424            3
601-800        247    215    7    469            4
801-1000        198    247    0    445            5
1001-1200    53    227    0    280            6
1201-1400    14    196    65    275            7
1401-1600    30    158    132    320            8
1601-1800    22    109    138    269            9
1801-2000    12    47    122    181            10

4. Calculatrice raster depuis la couche de somme des RGB pour créer une couche qui redéfini les classes d'IDPR avec une formule conditionnelle :

1*("Somme_RGB_IDPR@1"=337)+
2*("Somme_RGB_IDPR@1"=375)+
etc...

Le multiplicateur en début de ligne donne la valeur de classe (1,2,3 ... 10) si la condition entre parenthèse est vraie et renvoie donc 1

Au final j'obtient une couche raster avec une seule bande dont la valeur va de 1 à 10, plus qu'à affiner l'affichage à partir des couleurs initiales.