Differences

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

--- tecnica:gps_cartografia_gis:georeferencing [2009/09/28 14:04] – 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]] come PNG, non georiferito!
-Un primo tentativo può essere quello di fornire un //world file//, come
-documentato da [[http://mapserver.org/input/raster.html#georeference-with-world-files|Georeference with World Files]]. Il problema è trovare i parametri! Per fortuna abbiamo un GeoTiff del TrueMarble alla stessa risoluzione, con **''gdalinfo''** vediamo i parametri:
-<code>
-gdalinfo TrueMarble.2km.21600x10800.tif
-...
-Origin = (-180.000000000000000,90.000000000000000)
-Pixel Size = (0.016666666666667,-0.016666666666667)
-...
-</code>
-Un trucco per far creare il world file con gdal, a partire da un GeoTiff georiferito è il seguente:
-<code>
-gdal_translate -co "TFW=YES" TrueMarble.2km.21600x10800.tif temp.tif
-rm temp.tif
-</code>
-Scriviamo 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//.
-Convertiamo quindi 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>
-Se invece vogliamo un **GeoTiff completo**, che non ha bisogno del world file perché incorpora tutte le informazioni geografiche necessarie, 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.
-Se il file è troppo grande conviene realizzare //le piramidi//,
-gdalado