Comment changer la couleur du géotiff «sans données» en blanc / transparent plutôt qu'en noir?

Je prépare les données à inclure dans un nœud de géo-réseau, où les cartes sont affichées via le géoserveur.

Le problème est que les cellules qui n'ont "aucune donnée" s'affichent en noir dans le géotiff et je souhaite que la couleur soit définie sur blanc ou transparent.

J'utilise ArcView et je n'ai accès ni à Spatial Analyst ni à 3D Analyst.

Je préférerais ne pas avoir à attribuer un nouveau fichier sld et j'utilise actuellement "Utiliser la palette de couleurs" lors de l'exportation à partir d'ArcGIS Desktop.

Dans le concept, vous devez déterminer quelles valeurs de pixels représentent des cellules sans données, puis modifier votre table de couleurs pour attribuer les valeurs ROUGE, VERT, BLEU (RVB) dans la table de couleurs à 255255255 pour cette valeur de cellule.

Si vous ne disposez pas de la licence appropriée pour le faire dans Arc, vous pouvez utiliser l'outil gdal_translate (l'un des outils de ligne de commande gdal / ogr)

1. exporter la table des couleurs de votre tif sous forme de fichier XML.
   gdal_translate -of VRT myImage.tif outColorTable.vrt

2. ouvrez le fichier .vrt dans un éditeur de texte et recherchez l'entrée de la valeur de vos pixels sans données. Il n'y a pas de numéro d'index pour les valeurs des cellules, vous devez donc compter à rebours jusqu'à l'entrée correcte.

3. modifier cette valeur pour être <Entry c1="255" c2="255" c3="255" c4="255"/>

4. maintenant, créez un nouveau tif basé sur la nouvelle table des couleurs
   gdal_translate outColorTable.vrt newCorrected.tif

(exemple ajouté sur la base du commentaire ci-dessous)

Lorsque j'exécute la commande de l'élément n ° 2 ci-dessus et ouvre le fichier .vrt dans un éditeur de texte, je vois le contenu XML.

Une partie du fichier comprend des lignes qui ressemblent à ceci, j'ai coupé la plupart d'entre elles.

J'ai écrit un petit script gdal python pour cela, gdalsetnull.py , un simple programme en ligne de commande pour définir la valeur de raster spécifiée NODATA, sans créer de nouveau raster. Le chemin le plus simple pour installer GDAL sur Windows est via OSGeo4W . Exemples:

python gdalsetnull.py foobar.tif 0            # pure black is transparent
python gdalsetnull.py foobar.tif 0 255 0      # pure green is
python gdalsetnull.py foobar.tif 50 23 77 100 # arbitrary value in 4 band image

Si cela ne vous dérange pas de traiter les frais généraux et de créer une nouvelle image plutôt que de les modifier sur place, il y a aussi gdal_calc qui permet de changer les valeurs des cellules plutôt que de simplement mettre à jour les métadonnées.

Définissez les valeurs de zéro et inférieures à null:

gdal_calc.py -A input.tif --outfile=result.tif --calc="A*(A>0)" --NoDataValue=0

Outre le fait que vous préfériez NE PAS créer un nouveau fichier SLD, voici une entrée de blog sur la façon de convertir la palette de couleurs interne d'un GeoTIFF en SymbologyEncoding / SLD et de définir une autre couleur à cet endroit . Dans le SLD dérivé, il est facile de définir exactement quelles valeurs rendre avec une couleur 100% transparente.