Discussion :

[vbain]

Bonjour,

pour un besoin de cartographie sur une zone assez vaste, je renonce à travailler avec les données lidarhd, et reviens vers le thème altimétrie du RGE.
je me heurte à ce qui pourrait être une erreur de formatage de mes requêtes au serveur https://data.geopf.fr/wms-r/wms lorsque j'interroge la couche ELEVATION.ELEVATIONGRIDCOVERAGE.HIGHRES.
Les deux captures d'écran jointes présentent un extrait du rgealti de 150 m de côté au pas de 5 m. Elles mettent en évidence un motif d'erreur d'échantillonnage (enfin ce que j'interprète comme tel...). Les requêtes sont les suivantes :
- capture 1 :

Code :

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

https://data.geopf.fr/wms-r/wms?LAYERS=ELEVATION.ELEVATIONGRIDCOVERAGE.HIGHRES&EXCEPTIONS=text/xml&FORMAT=image/geotiff&SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&STYLES=&CRS=EPSG:2154&BBOX=1007500,6471500,1007650,6471650&WIDTH=30&HEIGHT=30

- capture 2 :

Code :

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

https://data.geopf.fr/wms-r/wms?LAYERS=ELEVATION.ELEVATIONGRIDCOVERAGE.HIGHRES&EXCEPTIONS=text/xml&FORMAT=image/geotiff&SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&STYLES=&CRS=EPSG:2154&BBOX=1007500,6471485,1007650,6471635&WIDTH=30&HEIGHT=30

C'est-à-dire que si je décale ma fenêtre de requête (ici de 15 m en ordonnée), je vois que le motif se décale lui aussi, ce qui me fait dire que le souci est une erreur d'échantillonnage (et non la donnée source).
Cette anomalie vous suggère-t-elle une réponse ?
Merci,

Vincent

[vbain]

Un complément utile à mon message précédent :
j'ai récupéré la donnée sous la forme d'une archive départementale, ici la Savoie.
Si je la compare avec ce que je récupère via le flux wms, il apparaît des différences métriques, avec là encore un motif curieux :

[azarzelli]

Bonjour,

Ce sont des artefacts de reprojection, il faut utiliser pour du Lambert 93 la ressource suivante : RGEALTI-MNT_PYR-ZIP_FXX_LAMB93_WMS

Cordialement,

[vbain]

Merci Azarzelli de vous intéresser à mon sujet,

hélas j'avais aussi exploré cette option sans guère plus de succès (des différences de près de 3 m subsistent). J'ai aussi essayé d'interroger ces thèmes en pseudo-mercator (epsg:3857) mais sans succès non plus.
Je repasse bientôt un peu de temps sur cette question, et donnerai des nouvelles.

Edit du lendemain, à tête reposée

L'interrogation de la couche RGEALTI-MNT_PYR-ZIP_FXX_LAMB93_WMS renvoie bien exactement la même chose que les dalles départementales. Mon erreur venait d'un défaut de centrage de la grille de requête (les coordonnées doivent être justes modulo 0.5m... et je requêtais sur des entiers multiples de 5). Pour ma culture, les artefacts de reprojection dans le cas de l'interrogation de la couche ELEVATION.ELEVATIONGRIDCOVERAGE.HIGHRES correspondent-ils à une interpolation un peu rustique (du type plus proche voisin) ?
Enfin, désolé pour le bruit, quoique ma mauvaise expérience puisse servir à d'autres.

Bonne journée
V.

[azarzelli]

Concernant les artefacts de reprojection, c'est effectivement dû à l'interpolation. Il me semble qu'une interpolation linéaire ou bilinéraire n'est pas possible à cause du NODATA à -99999, mais j'avoue ne pas être un expert sur le sujet.

[vbain]

Il me semble qu'une interpolation linéaire ou bilinéraire n'est pas possible à cause du NODATA à -99999

En effet la grille d'interpolation rencontre un problème aux marges de la zone renseignée. Merci pour vos réponses,
V.