Quei mattacchioni della NASA fanno scaricare il BlueMarble come PNG, non georiferito!

Un primo tentativo può essere quello di fornire un world file, come documentato da Georeference with World Files. Il problema è trovare i parametri! Per fortuna abbiamo un GeoTiff del TrueMarble alla stessa risoluzione, con gdalinfo vediamo i parametri:

gdalinfo TrueMarble.2km.21600x10800.tif
...
Origin = (-180.000000000000000,90.000000000000000)
Pixel Size = (0.016666666666667,-0.016666666666667)
...

Un trucco per far creare il world file con gdal, a partire da un GeoTiff georiferito è il seguente:

gdal_translate -co "TFW=YES" TrueMarble.2km.21600x10800.tif temp.tif
rm temp.tif

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:

0.0166666667
0.0000000000
0.0000000000
-0.0166666667
-179.9916666667
89.9916666667

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):

convert world_topo_bathy.png -compress none world_topo_bathy.tif

Se invece vogliamo un GeoTiff (che non ha bisogno del world file perché incorpora le informazioni geografiche) si può usare il comando:

gdal_translate world_topo_bathy.png world_topo_bathy.tif

Nel GeoTiff risultante mancano tuttavia le informazioni relative al sistema di riferimento.

Se il file è troppo grande conviene realizzare le piramidi,

gdalado