[Gfoss] conversione LL->UTM wgs84, risultati diversi tra proj.4 e verto. perche?!

Local Chapters GFOSS.it
1

Scusate. Sono impazzito o la proiezione di un LL wgs84 come UTM è solo
questione di applicare una... proiezione, appunto, che è univocamente
definita?
Ho fatto dei test su dei punti GPS, riproiettandoli con proj.4 (usando
ogr), con arcgis e con Verto 3.

LL 40.18282800 9.32401828
Verto 4462069.308 546268.668
proj.4 4448099.986 527584.223
arcgis 4448042.856 527617.210 (a partire da 9,323870 40,182596,
perché così mi sono stati salvati in shapefile a partire da ogr...)

Lo scarto è minore tra arcgis e proj.4, mentre si distanzia molto in
Verto. Ma non dovrebbe essere univoca la trasformazione???
Ho un lapsus di topografia. Troppo tempo che non faccio riproiezioni! :slight_smile:

giovanni

2

Per passare da coordinate polari a cartesiane
le formule sono univoche, le trovi su wikipedia
(devi usare quelle per l'ellissoide e non quel,le per la sfera)
Da polari (geografiche) a cartesiane (geocentriche)
  è un processo puramente analitico
e non si fanno approssimazioni di sorta
(quindi il risultato della conversione non è soggetto a errori)
nel caso inverso si applica un processo iterativo per il calcolo della latitudine.
Diverso invece il discorso, se vuoi passare da coordinate lonlat (geografiche) a utm (cartografiche)
basate sulla proiezione della carta di gauss.
In questo caso si applicano varie formule tra cui le formule di Hirvonen
sono rapidamente convergenti e forniscono una approssimazione millimetrica.

Se trovo un po di tempo le converto seguendo i due tipi di conversioni
(non riproiezioni ... quelle si usano per i passaggi di datum)

spero di esserti stato un minimo utile ... e di non aver detto sciocchezze :wink:

ciao,

Massimo Di Stefano
massimodisasha@yahoo.it

epifanio on irc.freenode.net /join gfoss

p.s. passami lo shape con i punti in coordinate lonlat

Il giorno 26/mar/09, alle ore 00:00, G. Allegri ha scritto:

Scusate. Sono impazzito o la proiezione di un LL wgs84 come UTM è solo
questione di applicare una... proiezione, appunto, che è univocamente
definita?
Ho fatto dei test su dei punti GPS, riproiettandoli con proj.4 (usando
ogr), con arcgis e con Verto 3.

LL 40.18282800 9.32401828
Verto 4462069.308 546268.668
proj.4 4448099.986 527584.223
arcgis 4448042.856 527617.210 (a partire da 9,323870 40,182596,
perché così mi sono stati salvati in shapefile a partire da ogr...)

Lo scarto è minore tra arcgis e proj.4, mentre si distanzia molto in
Verto. Ma non dovrebbe essere univoca la trasformazione???
Ho un lapsus di topografia. Troppo tempo che non faccio riproiezioni! :slight_smile:

Chiacchiera con i tuoi amici in tempo reale! http://it.yahoo.com/mail_it/foot/*http://it.messenger.yahoo.com

3

Diverso invece il discorso, se vuoi passare da coordinate lonlat
(geografiche) a utm (cartografiche)
basate sulla proiezione della carta di gauss.
In questo caso si applicano varie formule tra cui le formule di Hirvonen
sono rapidamente convergenti e forniscono una approssimazione millimetrica.

Quindi questo è il caso del passaggio da geografiche wgs84 a
utms32/whs84. Credevo fosse univoco anche questo...
Ti allego uno shape con soli 2 punti...

Grazie Massimo,
Giova

testmulti.tar.gz (400 Bytes)

4

G. Allegri ha scritto:

Scusate. Sono impazzito o la proiezione di un LL wgs84 come UTM è solo
questione di applicare una... proiezione, appunto, che è univocamente
definita?

Come afferma giustamente Massimo, si tratta di una *conversione*
(passaggio da coordinate geografiche a proiettate a parita' di datum) ed
in particolare, osservando la longitudine, si tratta di passare da
EPSG:4326 a EPSG:32632.

Ho fatto dei test su dei punti GPS, riproiettandoli con proj.4 (usando
ogr), con arcgis e con Verto 3.

LL 40.18282800 9.32401828
Verto 4462069.308 546268.668
proj.4 4448099.986 527584.223
arcgis 4448042.856 527617.210 (a partire da 9,323870 40,182596,
perché così mi sono stati salvati in shapefile a partire da ogr...)

Lo scarto è minore tra arcgis e proj.4, mentre si distanzia molto in
Verto. Ma non dovrebbe essere univoca la trasformazione???
Ho un lapsus di topografia. Troppo tempo che non faccio riproiezioni! :slight_smile:

Male! Un esperto in GIS dovrebbe farle molto spesso! :slight_smile:
Nel tuo caso si tratta di un classico errore legato al formato
*sessadecimale* (gg.dddddddd) delle tue coordinate geografiche. proj.4 se
le divora senza problemi per default, mentre altre applicazioni
potrebbero richiedere altri formati (es. *sessagesimale*, gg.ppssdddd) e
quindi prima di effettuare conversioni e/o trasformazioni e' opportuno
controllarne le impostazioni sulle unita' di misura degli angoli.
In definitiva, le coordinate UTM WGS84 33N sono 4448099.99 527584.22.
Spero di averti chiarito un pò le idee.

ciao
Antonio

5

Come afferma giustamente Massimo, si tratta di una *conversione*
(passaggio da coordinate geografiche a proiettate a parita' di datum) ed
in particolare, osservando la longitudine, si tratta di passare da
EPSG:4326 a EPSG:32632.

Sì, scusate. Ho usato il termine sbagliato, ma è chiaro che non si
tratta di una riproiezione...

Male! Un esperto in GIS dovrebbe farle molto spesso! :slight_smile:

Eheh, per fortuna lavoro su dati che non lo richiedono spesso! :))))

In definitiva, le coordinate UTM WGS84 33N sono 4448099.99 527584.22.

Ah, fa piacere. Verto richiede che li metta esattamente come li ho
riportati... e torno con quel risultato. Possibile?!

6

Ah, fa piacere. Verto richiede che li metta esattamente come li ho
riportati... e torno con quel risultato. Possibile?!

La stessa frase, in italiano, sarebbe: verto richiede che i dati
vengano forniti nel formato che ho indicato nella prima email, ma il
risultato (anch'esso riportato nella mia mail) è ben diverso da quello
che tu mi indiche, e che tornerebbe anche a me. Possibile?!

7

G. Allegri ha scritto:

In definitiva, le coordinate UTM WGS84 33N sono 4448099.99 527584.22.

scusa, qui è 32N (la forza dell'abitudine...)

Ah, fa piacere. Verto richiede che li metta esattamente come li ho
riportati... e torno con quel risultato. Possibile?!

Sei sicuro? Prova pure con Cartlab e Traspunto oppure fai la
trasformazione al volo con un desktop gis e vedrai.

ciao
Antonio

8

G. Allegri ha scritto:

Ah, fa piacere. Verto richiede che li metta esattamente come li ho
riportati... e torno con quel risultato. Possibile?!

La stessa frase, in italiano, sarebbe: verto richiede che i dati
vengano forniti nel formato che ho indicato nella prima email, ma il
risultato (anch'esso riportato nella mia mail) è ben diverso da quello
che tu mi indiche, e che tornerebbe anche a me. Possibile?!

Certo che no!

Antonio

9

Confermo che Cartlab rende lo stesso risultato di proj.4
Mi sa che chi mi ha mandato i risultati abbia impostato male qualcosa
di errato in Verto... Verificherò.
Grazie!

Il 26 marzo 2009 11.13, Antonio Falciano <afalciano@yahoo.it> ha scritto:

G. Allegri ha scritto:

In definitiva, le coordinate UTM WGS84 33N sono 4448099.99 527584.22.

scusa, qui č 32N (la forza dell'abitudine...)

Ah, fa piacere. Verto richiede che li metta esattamente come li ho
riportati... e torno con quel risultato. Possibile?!

Sei sicuro? Prova pure con Cartlab e Traspunto oppure fai la
trasformazione al volo con un desktop gis e vedrai.

ciao
Antonio

10

On Thu, Mar 26, 2009 at 12:00:10AM +0100, G. Allegri wrote:

Scusate. Sono impazzito o la proiezione di un LL wgs84 come UTM è solo
questione di applicare una... proiezione, appunto, che è univocamente
definita?

Come utm con che ellissoide e datum? Occhio che magari stai usando ellissoidi
o datum differenti nei vari programmi e il risultato ovviamente cambia.

Ho fatto dei test su dei punti GPS, riproiettandoli con proj.4 (usando
ogr), con arcgis e con Verto 3.

LL 40.18282800 9.32401828
Verto 4462069.308 546268.668
proj.4 4448099.986 527584.223
arcgis 4448042.856 527617.210 (a partire da 9,323870 40,182596,
perché così mi sono stati salvati in shapefile a partire da ogr...)

--
Francesco P. Lovergine

11

Come utm con che ellissoide e datum? Occhio che magari stai usando ellissoidi
o datum differenti nei vari programmi e il risultato ovviamente cambia.

No no, si sta parlando sempre di UTM/WGS84. Mi viene però il dubbio
che ci mi ha mandato i dati di Verto abbia fatto questa
considerazione...

12

G. Allegri ha scritto:

Come utm con che ellissoide e datum? Occhio che magari stai usando ellissoidi
o datum differenti nei vari programmi e il risultato ovviamente cambia.

No no, si sta parlando sempre di UTM/WGS84. Mi viene però il dubbio
che ci mi ha mandato i dati di Verto abbia fatto questa
considerazione...

Per fugare ogni dubbio (forse stamattina sono stato poco chiaro):
le coordinate piane in Verto sono state ottenute in maniera errata
considerando le coordinate geografiche nel formato sessagesimale e non
in quello sessadecimale, cosi' come dovrebbero essere. E' possibile
verificarlo anche in Cartlab, impostando il formato delle coordinate
geografiche come sessagesimale, oppure in proj.4 formattando
adeguatamente le coordinate in ingresso:

cs2cs +proj=latlong +datum=WGS84 +to +proj=utm +zone=32 +datum=WGS84

# formato sessagesimale (gg.ppssdddd) --> coordinate errate in Verto
9d32'40.1828" 40d18'28.2800"
546268.67 4462069.31 0.00

# formato sessadecimale (gg.dddddddd) --> coordinate esatte
9.32401828 40.18282800
527584.22 4448099.99 0.00

ciao
Antonio

13

Grazie Antonio per il tuo test in Verto. Questo mi conferma che chi mi
ha passato i dati ha commesso qualche errore. L'unica cosa di cui sono
certo ¨¨ che abbiamo lavorato tutti in gradi decimali. Come puoi vedere
dalle precedenti email non ho mai parlato di gradi-minuti-secondi. Pi¨´
priobabile errore di Fuso, o di settaggio del datum... Bho, domani
capir¨°.

Il 26 marzo 2009 20.04, Antonio Falciano <afalciano@yahoo.it> ha scritto:

G. Allegri ha scritto:

Come utm con che ellissoide e datum? Occhio che magari stai usando ellissoidi
o datum differenti nei vari programmi e il risultato ovviamente cambia.

No no, si sta parlando sempre di UTM/WGS84. Mi viene per¨¾ il dubbio
che ci mi ha mandato i dati di Verto abbia fatto questa
considerazione...

Per fugare ogni dubbio (forse stamattina sono stato poco chiaro):
le coordinate piane in Verto sono state ottenute in maniera errata
considerando le coordinate geografiche nel formato sessagesimale e non
in quello sessadecimale, cosi' come dovrebbero essere. E' possibile
verificarlo anche in Cartlab, impostando il formato delle coordinate
geografiche come sessagesimale, oppure in proj.4 formattando
adeguatamente le coordinate in ingresso:

cs2cs +proj=latlong +datum=WGS84 +to +proj=utm +zone=32 +datum=WGS84

# formato sessagesimale (gg.ppssdddd) --> coordinate errate in Verto
9d32'40.1828" 40d18'28.2800"
546268.67 4462069.31 0.00

# formato sessadecimale (gg.dddddddd) --> coordinate esatte
9.32401828 40.18282800
527584.22 4448099.99 0.00

ciao
Antonio