OpenAerialMap è un bel progetto con un obiettivo interessante: raccogliere foto aeree di dominio pubblico fatte in tutto il mondo, e renderle consultabili tramite un’unica interfaccia coerente. Chiunque abbia a disposizione dati di questo tipo, li potrà caricare sul sito e rendere disponibili al mondo.
Il catalogo è ancora limitato, ma negli stati uniti la copertura è buona anche ad una certa risoluzione. Per l’Italia ancora nulla ad alta risoluzione, ma c’è del materiale utile per carte regionali. Si tratta di immagini i-Cubed Landsat, sulle quali purtroppo cade (al momento) quanto detto in premessa, in quanto non è chiaro con che tipo di licenza vengano rilasciate. Si legge infatti:
[...]
it is best to assume that this imagery can not be used outside of OpenAerialMap
[...]
Tutto il catalogo è accessibile tramite WMS, e quindi potrete anche scaricare le immagini sul vostro PC e farne l’uso che la licenza vi consente. Per fare il download dell’immagini potrete ancora una volta usare FWTools:
- scaricate FWTools: Linux, Windows
- installatelo
- scaricate il file xml che descrive il servizio WMS di OpenAerialMap (http://openaerialmap.org/static/gdal_wms.xml) e salvatelo nella cartella in cui avete installato FWTools
- aprite la shell di FWTools
- scrivete un comando con questa sintassi: gdal_translate -projwin 12.06979 38.59784 15.86623 36.0124 -outsize 5000 5000 -of JP2KAK gdal_wms.xml palermo.jp2
FWTools è onnivoro in questo contesto e può quindi gestire come source fonti WMS. Con questa sintassi scaricherò una foto che copre la Sicilia in formato JPEG2000, che potrei usare per una carta al 250.000 (qui sotto uno screenshot dell’area). Le coordinate che leggete nel comando sono quelle del vertice in alto a sinistra e di quello in basso a destra.
Il post è (quasi) una traduzione molto veloce di un post del bravissimo Christopher Schmidt.
Tag: fwtools gdal_translate geodata wms
Sono un geologo marino e nel 1999 ho iniziato con Antonio D'Argenio la vita da libero professionista nel mondo dei Sistemi Informativi Geografici. Nei primi anni - di tanto studio e di poche fatture - è stato fatale l’incontro con MapServer: il web-mapping e le applicazioni web sono diventati due pezzi importanti della mia professione.
Negli ultimi anni mi sono occupato anche di formazione e di recente sono entrato nello strano mondo del telerilevamento di prossimità.
Mi piace molto leggere di GIS, informazione geografica, e neogeography, e parlarne con gli amici che incontro spesso al bar dietro al router.
I contenuti potrebbero non essere più adeguati ai tempi!
By Tiziano on dic 14, 2007
Ciao,
ho provato il tuo script. E’ fantastica questa cosa. Tuttavia non appena ho cambiato coordinate non mi ha funzionato. Il comando che ho inserito è:
gdal_translate -projwin 24.93333 46.38333 24.21666 47.10000 -outsize 5000 5000 -of JP2KAK gdal_wms.xml riyadh.jp2
e l’errore che mi ha dato è stato:
Input file size is 1073741824, 536870912
Computed -srcwin 611237465 130091785 -2137551 -2137551 from projected window.
-srcwin 611237465 130091785 -2137551 -2137551 falls outside raster size of
1073741824×536870912
or is otherwise illegal
Scusa l’ignoranza, ma se mi puoi aiutare te ne sarei grato
Grazie
By andrea on dic 14, 2007
Ciao Tiziano,
credo tu abbia sbagliato l’ordine della coppia di coordinate.
Devi inserire prima le coordinate del vertice in alto a sinistra e poi quelle di quello in basso a destra:
gdal_translate -projwin 24.21666 47.10000 24.93333 46.38333 -outsize 5000 5000 -of JP2KAK gdal_wms.xml riyadh.jp2
Prova e fammi sapere.
Ciao,
a