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

Bisognerebbe creare un world file secondo la struttura documentata da 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:

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

Un trucco per creare il world file a partire da un GeoTiff usando i tool gdal è il seguente:

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

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:

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.

Una soluzione non ottimale è quella di convertire 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

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:

gdal_translate -a_srs EPSG:4326 world_topo_bathy.png world_topo_bathy.tif

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 creare le overview, dette anche piramidi. 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:

gdaladdo TrueMarble.8km.5400x2700.tif 2 4 8

Se invece vogliamo tenere le overview separate dal file originale:

MapServer e la libreria Gdal accedono automaticamente alle overview… ?