pablo77 ha scritto:
Ciao a tutti
innanzitutto grazie per lo sbattimento di spiegare ai comuni mortali come scaricare ste razzo di cartografie senza impazzire con dei mostruosi printscr-incolla in photoshop et simili.
Cmq volevo chiedervi un favore. Qualcuno di voi saprebbe indicarmi l'indirizzo ecpw da cui accedere alla ortofoto
2003 della Lombardia? Mi serve per una mia amica che deve fare la tesi sul parco delle groane. Le ho già tirato giù la ver. 1998 seguendo il metodo FWT, ma non ho la più pallida idea di dove trovare la 2003. Sul portale della regione Lombardia c'è, ma se devo fare copia-incolla su una cartina di 20x10km divento pazzo !
grazie ciaoo
Guarda il primo messaggio di questo post che trovi gli indirizzi di tutte le regioni.
Mazzokan ha scritto:
Segnalo, a chi non lo sapesse, che senza comprare programmi commerciali con i semplici tool si possono fare molte più cose nella creazione di file ECW.
Mi riferisco in particolare a tre aspetti
1) superare il limite di 500 mega
2) mantenere la georeferenziazione (proiezione e datum).
3) modulare il livello di compressione dell'ECW (di dafault il 75%).
Per esempio con la linea seguente
gdal_translate -of ecw -co "LARGE_OK=YES" -co "PROJ=NUTM33" -co "DATUM=WGS84" -co "TARGET=85" C:\FWTools1.0.0a8\immagine.tif C:\FWTools1.0.0a8\immagine_percento.ecw
Comando di superare il limite di 500 mega (comando "LARGE_OK=YES") come descritto già da altri. Indico anche di usare la proiezione UTM fuso 33 (comando -co "PROJ=NUTM33") ed il datum WGS84 ( - co "DATUM=WGS84"). Infine dico anche di crearmi un file ECW compresso di circa l'85% (-co "TARGET=85").
Aggiungo di fare attenzione al fuso (l'Italia è divisa fra il 32 ed il 33). Per quanto riguarda la compressione se non specificata è pari ad un valore di 75 che corrisponde ad una immagine compressa apri ad 1:4 dell'originale (25%). Impostando il valore ad 85 si ottiene un file di quasi il 20% dell'originale (anche di meno se ci sono distese di acqua e zone bianche di contorno) con qualità mantenuta. Arrivando al valore di 90 si ottiene un file di circa il 15% dell'originale con modesta puntinatura solo per elevati ingrandimenti.
Non posso che confermare! Ci sono varie opzioni aggiuntive che permettono di potenziare ed ottimizzare notevolmente il funzionamento di gdal. E' comodo anche trasformare una mappa in una semplice immagine JPEG o GIF da utilizzare successivamente.
Segnalerei sicuramente
gdalinfo.exe nomefile.ext . Come si può intuire, mostra le informazioni di un file mappa. E' possibile vedere ad esempio:
- DATUM e Proiezione
- Dimensione in pixel ( che con alcuni programmi è possibile utilizzarli come coordinate per calibrare la carta tramite gli angoli)
- Le coordinate dei 4 angoli + quella centrale molto utile se si vuole calibrare la mappa manualmente, come magari l'immagine create prima.
- molti altri
Capita a volte che aprendo la carta si ottenga l'errore perché non contiene alcun dato, ma aprendola tramite importa mappa e specificando datum/proiezione questa viene letta. Come è possibile se prima non si riusciva? Perché la carta conteneva
solo le coordinate dei vari punti di calibrazione, ma non le informazioni sul datum/proiezione. Ecco perché per alcune carte è possibile aprirle semplicemente, altre no. Gdalinfo mostra se una mappa contiene tutte queste info
![Pinkbanana :yeah!: :yeah!:](/styles/uix/xenforo/smilies1/pinkbanana.gif)
. Se mancano, e sono cmq note, possiamo aggiungerle come ha mostrato prima Mazzokan in modo da includerle permanentemente nel file.
Esempio, un mio pezzo di mappa verificato con gdalinfo. Dopo il primo comando si vede chiaramente che non contiene /datum/proiezione:
Coordinate System is ''. Il secondo comando è successivo alla trasformazione con gdal (il nome finale è diverso, ma il file lo stesso). L'output è maggiormente dettagliato, si vedono le informazioni riguardo datum e proiezione, nonché le coordinate in gradi che sono gli angoli della carta. Georeferenziando un carta tramite gli angoli si utilizza un punto
assolutamente preciso perché è un pixel con una sua precisa coordinata. Basta inserire le relative coordinate prese da gdal è il gioco è fatto. Calibrando un punto a mano libera, per quanto ci si impegni c'é sempre qualche errore.
Compegps, ad esempio, permette di specificare la posizione dei punti di calibrazione non solo cliccando su un punto della mappa (ma quanto siete precisi a riportare un punto sulla carta??), ma anche specificando le coordinate in pixel della stessa. Ad esempio l'angolo in alto a sinistra, che è l'origine, ha coordinate in pixel
X/Y di
0,0. Nel primo punto specificate queste e poi inserite le coordinate vere e proprie che vi detto gdalinfo in
Upper left. Passando al punto successivo appare il punto di calibrazione.
L'angolo opposto, in
basso a destra lo si può ricavare tramite l'output di gdal
Size is x, y che è la dimensione. Nella mia immagine è 10000 per 10000 pixel. Quindi la posizione è X/Y 10000/10000.
Mazzokan ha scritto:
Queste cose le ho scoperte semplicemente leggendo l'help in linea del programma.
Consiglio a tutti di mantenere la georeferenziazione originale delle immagini ECW perchè in un secondo tempo sarà possibile riottenere delle coordinate per eventuali punti di interesse.
Spero di essere stato utile alla discussione
![Icon Biggrin :-) :-)](/styles/uix/xenforo/smilies1/icon_biggrin.gif)
.
PS: il trucchetto per superare i 500 mega è descritto direttamente nell'help del programma (ma sono pazzi!
![Smile No0004 :ueh: :ueh:](/styles/uix/xenforo/smilies1/smile_no0004.gif)
).
Se non fosse stato per l'help ( e i molti Gb di file utilizzati per le prove
![Pc :pc: :pc:](/styles/uix/xenforo/smilies1/pc.gif)
) sarei rimasto all'oscuro di queste preziose info. Come dici tu, consiglio pure io di fare attenzione alla presenza dei dati di georeferenziazione originali. Quando stavo provando ad unire 25 pezzi di ortofoto con Ermapper, non riusciva mai a caricarne 9. Prova e riprova, alla fine ho scoperto che mancavano proprio alcuni di questi dati...prontamente aggiunti con gdal.
Per aggiungere i metadati, oltre al metodo di Mazzokan si può utilizzare l'opzione
-a_srs in questa forma
-a_srs "+proj=UTM +zone=32 +datum=WGS84". Ovviamente inserire i dati corretti per la carta che state usando.
Per chi non si vuole sbattere, tramite il sw
ECW SPY http://www.digitalearth.com.au/downloads/ecw-spy/ è possibile controllare la presenza dei dati
![Richiamapapere :richiamo_papere: :richiamo_papere:](/styles/uix/xenforo/smilies1/richiamapapere.gif)
e quali sono.
In fatto di dimensione, sono arrivato ad ottenere file ecw di 1,2 Gb con gdal. La possibilità di superare il limite dei 500 mb con un sw free, non è poco.