Differences

This shows you the differences between two versions of the page.

--- tecnica:gps_cartografia_gis:gdal_raster [2011/11/15 11:28] – [Raster con GDAL e MapServer] niccolo
+++ tecnica:gps_cartografia_gis:gdal_raster [2013/10/01 23:08] (current) – [Overview Images (piramidi)] niccolo
@@ Line 84: / Line 84: @@
 </code>
+===== Georiferire un file PNG =====
+Quei mattacchioni della NASA fanno [[http://visibleearth.nasa.gov/view_set.php?categoryID=2355|scaricare il BlueMarble]] in formato PNG, non georiferito!
+Bisognerebbe creare un //world file// secondo la struttura documentata da [[http://mapserver.org/input/raster.html#georeference-with-world-files|Georeference with World Files]]. Il problema è sapere i parametri! Per fortuna abbiamo un GeoTiff del TrueMarble alla stessa risoluzione e nello stesso sistema di riferimento, con **''gdalinfo''** vediamo i parametri del GeoTiff:
+<code>
+gdalinfo TrueMarble.2km.21600x10800.tif
+...
+Origin = (-180.000000000000000,90.000000000000000)
+Pixel Size = (0.016666666666667,-0.016666666666667)
+...
+</code>
+Un trucco per creare il world file a partire da un GeoTiff usando i tool gdal è il seguente:
+<code>
+gdal_translate -co "TFW=YES" TrueMarble.2km.21600x10800.tif temp.tif
+mv temp.tfw TrueMarble.2km.21600x10800.tfw
+rm temp.tif
+</code>
+Creiamo il world file con lo stesso nome del file ''.png'', ma estensione **''.wld''** oppure **''.tfw''** (un world file con tale estensione viene trovato sia da MapServer che dalla libreria GDAL). Le coordinate del centro del pixel in alto a sinistra sono le cordinate dell'origine riportate da ''gdalinfo'', meno la metà di un pixel:
+<file>
+.0166666667
+.0000000000
+.0000000000
+-0.0166666667
+-179.9916666667
+.9916666667
+</file>
+Peccato che la velocità di accesso non sia accettabile, il manuale di GDAL recita: //PNG files are linearly compressed, so random reading of large PNG files can be very inefficient//.
+Una **soluzione non ottimale** è quella di convertire in Tiff non compresso con ImageMagick (il world file da usare rimane esattamente lo stesso):
+<code>
+convert world_topo_bathy.png -compress none world_topo_bathy.tif
+</code>
+La **soluzione migliore** è creare un **GeoTiff completo**, che non ha bisogno del world file perché incorpora tutte le informazioni geografiche necessarie. A tale scopo si può usare il comando:
+<code>
+gdal_translate -a_srs EPSG:4326 world_topo_bathy.png world_topo_bathy.tif
+</code>
+Nel GeoTiff risultante vengono incorporate sie le informazioni relative al sistema di riferimento che quelle contenute nel world file.
+===== Overview Images (piramidi) =====
+Se il file è troppo grande conviene creare le overview, dette anche //piramidi//, altrimenti uno zoom alla massima estensione obbliga a leggere e ridimensionare al volo tutto il raster.
+In genere si generano delle copie dell'immagine con risoluzione 1/2, 1/4, 1/8, ... dell'originale. Le overview possono essere contenute all'interno dello stesso GeoTiff:
+<code>
+gdaladdo TrueMarble.8km.5400x2700.tif 2 4 8
+</code>
+Se invece vogliamo tenere le overview in un file separato (che avrà estensione **''.tif.ovr''**) ed utilizzare un algoritmo di ricampionamento migliore (//gauss//, ma anche //average// è migliore del //nearest// predefinito):
+<code>
+gdaladdo -ro -r gauss TrueMarble.8km.5400x2700.tif 2 4 8
+</code>
+Con le versioni di GDAL < 1.6.0 si deve usare una sintassi diversa, il file esterno con le overview avrà formato //Erdas Image// ed estensione **''.aux''**:
+<code>
+gdaladdo --config USE_RRD YES TrueMarble.8km.5400x2700.tif 2 4 8
+</code>
+MapServer e la libreria Gdal accedono automaticamente alle overview anche se si trovano nel file separato **''.aux''**.
+===== Creare una overview di una copertura tiled =====
+Anche se nelle tile sono state aggiunte le overview, ai fattori di scala maggiori (zoom minore) è necessario accedere a tutti i file (tile) della copertura per generare la mappa. Questo fa decadere di molto le prestazioni, è quindi opportuno creare una **versione ridotta della copertura**, contenuta in un **singolo file**.
+Tale file conviene che sia TIFF memorizzato a blocchi ed eventualmente deve comprendere le overview opportune. Ecco come usare **''gdal_merge''** per ottenere il risultato:
+<code>
+gdal_merge.py -o output.tif -of GTiff -co TILED=YES -co TFW=YES -ps 0.9765625 0.9765625 2001/*.jpg
+</code>
+La dimensione del pixel (in metri) viene calcolata dividento la larghezza della copertura per il numero di pixel desiderati. Ad esempio se si ha una superficie larga 16 km e si vuole creare un'immagine larga 16384 pixel:
+<code>
+(670510 - 654510) / 16384 = 0.9765625
+</code>