Gaetanbrl
commented
7 months ago Informations générales
Le type de données OpenLayers concerné est le WebGL (GeoTiff) : https://openlayers.org/en/latest/apidoc/module-ol_layer_WebGLTile-WebGLTileLayer.html
Il permet d'utiliser une source de couche OpenLayers GeoTIFF via la librairie GeoTiff.js : https://geotiffjs.github.io/
En prérequis, un workshop permet de cerner la création / affichage / rendu de ce type de couche () :
https://openlayers.org/workshop/en/cog/colormap.html
La réalisation d'un style semble donc utiliser les "expressionValue" Openlayers qui est une syntaxe spécifique à OpenLayers. Cette syntaxe permet de réaliser un style via un ensemble d'expression logique.
On peut notamment faire une classification selon une liste de valeur, des croisements (opérations) entre plusieurs bands d'un raster, etc. (voir les URLs pour plus d'informations).
Comment représenter une couche ?
Si le flux est un WMS, il est possible d'utiliser un SLD avec GeoServer. Dans notre cas, c'est le WPS qui retourne un TIFF directement.
La documentation OpenLayers pour le type WebGL ne semble pas permettre l'utilisation d'un SLD.
Nous avons donc testés de reprendre la liste des valeurs du SLD, et de l'appliquer dans une expressionValue OpenLayers afin de représenter le TIFF.
Voici un exemple pour réaliser une classification à partir de valeurs d'une bande d'un TIFF :
const bandValue = ['band', 1];
const layer = new ol.layer.WebGLTile({
source: src,
style: {
color: [
'case',
// couleur si la valeur de band est == à -0
['==', bandValue, -0],
['color', 0, 0, 0, 0],
// couleur si la valeur de band est <= à -4
['<=', bandValue, -4],
['color', 62,55,144, 1],
['<=', bandValue, -3.75],
['color', 164, 252, 60, 1],
// valeurs par défaut
['color', 0, 0, 0, 1]
]
}
})Le résultat n'étant pas satisfaisant pour LETG, cette issue à été créée.
Il est donc demander de créer un code complexe qui nécessite de :
- définir la valeur minimum pour une bande d'un TIFF (valeur min des pixel)
- définir la valeur maximum pour une bande d'un TIFF (valeur max des pixel)
- définir une palette de couleur
- définir une série de valeurs entre l'ambplitude min/max selon un écart donné (ex: step de 0,25 pour une série comprise entre -4 et +6)
- affecter pour chaque valeur de la série une couleur selone la palette choisie
- Ne pas affecter de couleurs pour certaines valeurs comprises entre -0,25 et +0,25
En l'état, rien n'est prévu dans OpenLayers pour cette réalisation déjà annoncée comme complexe. Il semble qu'il serait toutefois possible d'avoir ce type de rendu avec une opération d'interpolation disponible dans les expressionValue OpenLayers. Mais ce code nécessite l'ajout d'une librairie externe type colormap.js inspirée de la librairie maplotLibre.
Des essais ont au préalable était réalisés en dehors du mviewer et avec un TIFF issue de la comparaison de 2 MNT avec SAGA :
https://codesandbox.io/p/sandbox/cog-no-data-interpolate-forked-j2n8qz?file=%2Fmain.js
new TileLayer({
source: source,
style: {
color: [
"interpolate",
["linear"],
layer,
// le nombre 40 correspond au nombre de couleurs différentes à obtenir par valeur
// c'est le nombre de valeurs qui existent entre -4 et +6 à des pas de 0,25:
// 10 / 0,25 = ~40
...getColorStops("jet", -4, 6, 40, false),
],
},
}),La valeur "jet" correspond à la palette suivante dans colormap:
Ce code peut sembler intéressant, et nécessite un test dans mviewer.
Réalisation dans mviewer
Pour être utilisable dans mviewer, il faut créer une extension colormap qui n'est pas disponible directement en tant que librairie buildée (type CDN). Il faut donc utiliser une extension type "module" et recopier la librairie (ce qui n'est pas très long).
Ensuite, un test peut être réalisé dans un customLayer avec un TIFF issue d'une comparaison de 2 MNT SAGA.
Analyse des premiers résultats :
- Le code semble fonctionnel pour Boutro (TIFF simple issue des MNT GDAL) mais on observe que les valeurs indésirables sont encore représentées alors qu'on souhaiterait les masquer
- Les performances semblent en plus beaucoup dégradées pour VOUGO (TIFF avec une plus grande résolution issue des MNT SAGA)
- Les couleurs ne sont pas encore satisfaisantes, il faut regarder comment colormap réalise sont rendu
Quentin-RUAUD
Description
La légende du résultat de comparaison devrait être calculée de la manière suivante de -0,02m à 0,02m : blanc ou transparent du min à -0,02 : gamme de couleur allant du bleu foncé au vert clair par pas de 0,25 de 0,02 à max : gamme de couleur allant du jaune clair au rouge foncé par pas de 0,25
min et max étant les valeurs minimum et maximum du MNT issu du résultat de la comparaison