Discussion :

[SgtLame]

Bonjour à tous,

La requête suivante fonctionne, elle me renvoie bien une heightmap au format tif:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

https://wxs.ign.fr/altimetrie/geoportail/r/wms?LAYERS=ELEVATION.ELEVATIONGRIDCOVERAGE.HIGHRES&EXCEPTIONS=text/xml&FORMAT=image/tiff&SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&STYLES=&CRS=EPSG:2154&BBOX=1208316,6115176,1210364,6117224&WIDTH=2049&HEIGHT=2049

Mais si je change le format d'image souhaité en image/jpeg ou image/png, une image entièrement blanche m'est retournée. Pourtant, si j'ai bien compris https://wxs.ign.fr/altimetrie/geopor...etCapabilities, les formats jpeg et png devraient être supportés, non?

[Viglino]

Chez moi ça marche : https://wxs.ign.fr/altimetrie/geopor...4314.081675609

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

https://wxs.ign.fr/altimetrie/geoportail/r/wms?SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&FORMAT=image%2Fjpeg&TRANSPARENT=true&LAYERS=ELEVATION.ELEVATIONGRIDCOVERAGE.HIGHRES&WIDTH=256&HEIGHT=256&CRS=EPSG%3A3857&STYLES=&BBOX=-78271.51696402207%2C6105178.323193599%2C-39135.75848201183%2C6144314.081675609

[xavier_dmz]

Bonjour,

Pour comparaison entre vos deux requêtes :

La requête de SgtLame qui renvoit une image blanche (je reproduis chez moi sur Firefox et Chromium) :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
* URL de base https://wxs.ign.fr/altimetrie/geoportail/r/wms
* paramètres :
    * LAYERS=ELEVATION.ELEVATIONGRIDCOVERAGE.HIGHRES
    * EXCEPTIONS=text/xml
    * FORMAT=image/jpeg
    * SERVICE=WMS
    * VERSION=1.3.0
    * REQUEST=GetMap
    * STYLES=
    * CRS=EPSG:2154
    * BBOX=1208316,6115176,1210364,6117224
    * WIDTH=2049
    * HEIGHT=2049

La version tiff de l'image contient bien de la donnée nettement visible. L'image blanche n'est donc pas expliquée par les données représentées.

La requête de Viglino qui renvoit une image en niveau de gris correspondant très nettement à une heightmap (je reproduis chez moi sur Firefox et Chromium) :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
* URL de base https://wxs.ign.fr/altimetrie/geoportail/r/wms
* paramètres :
    * SERVICE=WMS
    * VERSION=1.3.0
    * REQUEST=GetMap
    * FORMAT=image%2Fjpeg
    * TRANSPARENT=true
    * LAYERS=ELEVATION.ELEVATIONGRIDCOVERAGE.HIGHRES
    * WIDTH=256
    * HEIGHT=256
    * CRS=EPSG%3A3857
    * STYLES=
    * BBOX=-78271.51696402207%2C6105178.323193599%2C-39135.75848201183%2C6144314.081675609

J'ai gardé l'ordre original des paramètres au cas-où.

Les différences que j'observe entre les deux requêtes sont :

• L'échappement des caractères
• Les dimensions de l'image
• L'emprise géographique
• La projection
• L'ordre des paramètres
• Le paramètres "EXCEPTIONS=text/xml" dans la première requête, contre le paramètre "TRANSPARENT=true" dans al deuxième

[SgtLame]

J'ai repris la syntaxe exacte de ta requête, en y mettant mes valeurs de bbox, de CRS et de taille, peut-être qu'il y a avait un paramètre en plus qui me manquait, ou quelque chose clochait dans l'ordre des paramètres.
Et malheureusement, j'ai toujours une image blanche.

[Viglino]

contre le paramètre "TRANSPARENT=true" dans al deuxième

Au temps pour moi, TRANSPARENT=true vient du fait que la requête initiale demandait un png (et pas jpeg)...
donc pas nécessaire en jpeg

[xavier_dmz]

On m'a soufflé une cause probable dans l'oreillette.

• Dans les données originales les valeurs des pixel sont exprimées sous forme de flottants 32 bits. (signé ou pas, je ne sais pas) C'est aussi le cas dans l'image au format tiff.
• Dans les images retournées en jpeg ou png, les nombres sont exprimés en entiers 8 bits. (non signé, je crois)
• La conversiont de float32 vers int8 se fait par troncature.

Par conséquent toutes les altitudes supérieures à 255m sont corrompues lors de la conversion.

L'image blanche pourrait donc s'expliquer ainsi, mais je n'ai pas vérifié les valeurs dans les images.
Et l'image correcte obtenue par Viglino pourait s'expliquer par des altitudes comprises entre 0 et 255m, mais là encore je n'ai pas vérifié le contenu des images.

EDIT après vérification des histogrammes dans qgis :
L'image TIFF pour la première requête contient uniquement des valeurs supérieures à 886, donc hors de la plage [0; 255], qui sont bien corrompues par la conversion.
L'image TIFF pour la deuxième requête contient majoritairement des valeurs entre 50 et 250, donc qui seront conservées lors de la conversion.

[xavier_dmz]

Du coup si votre but est de récupérer la heightmap depuis ce service WMS, restez en TIFF.
Si votre but est de visualiser les altitudes dans des images plus classiques, je vois des solutions en changeant de couche ou de style, mais elles risquent de ne pas vous satisfaire :

• "LAYERS=ELEVATION.ELEVATIONGRIDCOVERAGE.HIGHRES&FORMAT=image/png&STYLES=hypso" :
  
  • même couche haute résolution, affichage en teintes hypsométriques.
  • On voit un peut de variation de relief, mais l'image me semble beaucoup moins exploitable que le TIFF en style "normal" (la valeur par défaut du paramètre "STYLES" pour cette couche.)
  • https://wxs.ign.fr/altimetrie/geopor...49&HEIGHT=2049
• "LAYERS=ELEVATION.ELEVATIONGRIDCOVERAGE.SHADOW&FORMAT=image/png&STYLES=estompage_grayscale"
  
  • couche basse résolution, affichage en estompage en niveaux de gris.
  • Exploitable, mais moins net à cause de la résolution
  • https://wxs.ign.fr/altimetrie/geopor...49&HEIGHT=2049

[SgtLame]

Merci beaucoup pour vos investigations, c'est très clair!
ça m'aurait bien arrangé de pouvoir travailler sur du jpeg ou du png, mais je vais donc rester sur du tiff.