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

[aperi2007]

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 :)
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) . :)
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
-------------- parte successiva --------------
Un allegato non testuale ? stato rimosso....
Nome:        gdalwarp_merge_sample.gif
Tipo:        image/gif
Dimensione:  15112 bytes
Descrizione: non disponibile
URL:         <http://lists.gfoss.it/pipermail/gfoss/attachments/20110906/d0d53330/attachment.gif>