[Gfoss] Come fondere un set di raster TIFF4 e ricondurli a uno solo.

Salve,

invio questo caso di uso che potrebbe essere utile ad altri.
L'obiettivo, come detto nell'oggetto, e' il seguente:

Riuscire a fondere un set di ratser georeferenziati TIFF4 (1 bit level) ottenendo un unico raster georeferenziato sempre Tiff4.

E' interessante anche sapere perche' si deve fare questo :slight_smile:
Ma la storia sarebbe lunga, per cui mi limito a dire che serve per venire incontro a l'impiegodi un famoso software GIS desktop (commerciale) che non riesce (bonta' sua) a gestire la trasparenza su un unico livello di colore (il bianco) sui cataloghi di raster.
Onde per cui con tale stupendo software il risultato con i cataloghi e' quello riportato nell'immagine allegata (sigh) . :slight_smile:
Il buffo e' che se invece di un catalogo si caricassero i rasters uno per uno lo stupendo software commerciale allora riesce a mettere il bianco a trasparente.

A questo punto l'unica strada se non si vuole costringere l'utente a caricarsi uno per uno tutti i raster e' ricondurli a uno soltanto.

Certo questo introduce altri problemi, come la pesantezza del raster risultante.
Ma questo e' una altra storia. Diciamo che รจ uno scotto da pagare se si vuole usare certi software.

Comunque venendo al "caso d' uso".
In questo caso il software GFoss che viene in soccorso del software commerciale che tentenna, e' gdal.
Due utilissime utility di gdal che servono a questo scopo sono:
gdalwarp e gdal_translate.
Con il loro aiuto si riesce a fondere un set incommensurabile di raster (nel mio caso arrivo anche a un raster di 40 Gbytes) e poi tramite gdal_translate riportarlo a un tiff4 di 200 Mbytes.

Di seguito la sequenza:

gdalwarp.exe -s_srs EPSG:3003 -t_srs EPSG:3003 -r near -ot Byte
-srcnodata "0" -of GTiff -co TFW=YES ../raster/list/*.tif
../raster/raster_result.tif

gdal_translate -co NBITS=1 -co COMPRESS=CCITTFAX4 -co TFW=YES
../raster/raster_result.tif ../raster/raster_result_as_tiff4.tif

I due parametri
-s_srs EPSG:3003 e -t_srs EPSG:3003
impongono che il raster di partenza e quello di arrivo abbiano il medesimo sistema di riferimento e il parametro
-srcnodata "0"
impone che il bianco sia trattato come trasparenza per evitare quello che si vede nell'immagine allegata.

Saluti,

Andrea

gdalwarp_merge_sample.gif