[GRASS-Italia] R: Esportare File GeoTiff

[Paolo Zatelli]

-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1

Ciao Marco,

Marco Pinelli ha scritto:
> Ciao Paolo,
> innanzi tutto ti ringrazio per la tua risposta così rapida ed esaustiva. 
> Io uso GRASS 6.3, che ho installato seguendo le istruzioni del sito
> ufficiale: ho aggiunto la repository e poi ho installato GRASS con apt-get
> install ecc (PS: credo sia da pacchetto così, e ti confesso che non so cosa
> preveda l'installazione "compilata"). Cercando di capire che versione di
> GDAL avessi, per scrivertelo, ho scoperto di non avere installato GDAL. Ma
> credevo di avere GDAL dato che r.in.gdal di GRASS funzionava. A questo punto
> mi sa che la mia situazione e' più grave del previsto. Sapresti dirmi quante
> altre cose non ho installato usando la procedura standard che ti ho
> descritto e che invece servirebbe avere per usare GRASS al pieno delle sue
> potenzialita'?
> Comunque ora ho installato l'ultima versione di GDAL e interrogando il nuovo
> file tif che ho creato, con TIFFINFO, mi da praticamente cio' che mi avevi
> copiato tu nella precedente mail dal tuo file di esempio.
> Tutto a posto quindi?? Nient'affatto. Perché comunque importando il tif
> creato (bene) con GRASS in ArcMAP non mi vede nessun sistema di riferimento
> (uso ArcMap a volte per visualizzare i risultati ottenuti con GRASS perché,
> secondo me, ha una grafica un po' più precisa e veloce. Solo per quello
> giuro!!), mentre con il tif aspect2.tif del tutorial si (se guardo, da
> ArcMAP, le proprietà del tif compaiono tutte i riferimenti del sistema di
> proiezione). Non so se c'entra ma ho notato che nel Tag 34737 del file del
> tutorial compare (come e' pubblicato sul sito): 
> Tag 34737: UTM Zone 13, Northern Hemisphere|clark66|, 
> mentre nel tuo esempio della mail precedente, e anche nei miei file in
> uscita ora che ho GDAL, compare:  
> Tag 34737: Transverse Mercator|international|
> e forse così lui non sa che zona UTM stiamo considerando (quindi credo che
> anche importando in ArcMAP il file visibilita_1.tif che hai creato tu come
> esempio per me non verrebbe riconosciuto il sistema di riferimento). Pensi
> che possa essere questo il motivo? Se si', come si fa a creare un tif come
> quello del tutorial? Sarebbe fantastico poterlo fare così' quando ci si
> scambiano tif (tipo i CTR) non ci sarebbe più' bisogno del file TFW e non si
> rischierebbe più di non sapere a quale sistema di proiezione appartengono le
> coordinate nel tif.
> Scusa la lunghezza della risposta e grazie ancora per la disponibilita'.
Un po' di chiarezza: UTM e Gauss Boaga sono due *proiezioni*, Roma40 e
NAD27 i sistemi di riferimento associati negli esempi sopra.
La proiezione di Gauss Boaga (quella attualmente usata dalla cartografia
ufficiale italiana) ha due fusi: ovest ed est (che corrispondono, circa,
ai fusi UTM 32 e 33, in realta' sono di 30'' piu' larghi).
La *zona* UTM, contraddistinta da una lettera, indica fasce di
latitudini e solitamente non e' usata in ambito GIS. Nei tag sopra non
e' indicata, in inglese si usa invece "zone" per l'italiano "fuso".

> Tag 34737: UTM Zone 13, Northern Hemisphere|clark66|,

indica il fuso UTM 13, emisfero nord, con ellissoide Clark66 (e' il fuso
UTM in cui cade la contea di Sperfish negli USA, i cui dati sono usati
nel tutorial).

> Tag 34737: Transverse Mercator|international|

indica che la proiezione e' una Mercatore trasversa su ellissoide
Internazionale (o di Hayford): per definire la Gauss Boaga servono altri
parametri, che sono indicati qui:

Tag 34736: 0.000000,9.000000,0.999600,1500000.000000,0.000000

sono il parallelo e meridiano centrale, il fattore di contrazione e le
false origini Est e Nord (in questo caso Gauss Boaga fuso ovest).

La grossa differenza tra UTM e Gauss Boaga e' che mentre l'UTM e'
completamente definito scegliendo il fuso, dato che tutti i GIS
conoscono i parametri dell'UTM, per Gauss Boaga il piu' delle volte si
devono indicare esplicitamente i parametri.
In sostanza credo che il problema qui sia in ArcMap, il quale non sa
utilizzare i tag che specificano i parametri della proiezione Gauss Boaga.

Come controprova puoi:
1. creare una location UTM/WGS84 ed esportare un tiff, se la mia ipotesi
e' giusta ArcMap dovrebbe leggere correttamente i parametri;
2. reimportare con r.in.gal il tuo tiff in GRASS e vedere se la
proiezione ed il sistema di riferimento sono riconosciuti.

Non conosco ArcMap, magari qualcuno che lo usa puo' darti informazioni
piu' precise.

ciao		Paolo

> 
> Ciao Marco
> 
> PS: per Luca: grazie Luca anche per il tuo consiglio, ma era proprio per
> evitare di dover costruire e portarsi dietro, oltre al tif, il file TFW che
> sono così interessato a costruire un unico file che contenga gia' tutto.
> Credo che aiuterebbe a fare meno errori (con le proiezioni) e a mantenere
> più ordinati i database. Comunque l'opzione tfw e' sempre valida. Grazie.
> 
> 


- --
-----BEGIN PGP SIGNATURE-----
Version: GnuPG v1.4.9 (GNU/Linux)
Comment: Using GnuPG with Mozilla - http://enigmail.mozdev.org

iEYEARECAAYFAkkN178ACgkQsbjfGi+aQw7k4QCfTJf1n7Q91sO4Ha5R9k0zwBEk
2/gAoNEY90oW6GZEmuIKgqawB4xd3S13
=D8ot
-----END PGP SIGNATURE-----