Differences

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

--- tecnica:gps_cartografia_gis:georeferencing [2011/11/15 10:24] – niccolo
+++ tecnica:gps_cartografia_gis:georeferencing [2011/11/15 10:30] (current) – removed niccolo
@@ Line 1: / Line 1: @@
-====== Raster Georeferencing ======
-===== 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 (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''**.