Un cambiamento d datum introduce per forza degli errori (anche se fatto
con Verto della IGM).
Inoltre con tutti i sw che non usano il grigliato igm ci si deve
attendere dai 5 metri di errore quadratico medio in su, fino a 150 metri
se non si inserisce nessun cambio di datum, ma solo di proiezione.
Il 1/2 millimetro dichiarato da Luca mi lascia perplesso, a meno che i
parametri di molodensky di ogr2ogr non siano proprio centrati su Trento.
Mi fa ipotizzare, invece, che anche la carta (precisa ed accurata al 1/2
mm sul terreno) usata per misurare il 1/2 mm fosse stata derivata da un
cambio di SRS (sto scherzando!). Inoltre bisognerebbe sapere rispetto a
cosa si misura l'errore. Se Luca può dirci qualcosa in più....
L'unica misura ufficiale (oltre che maggiormente precisa ed accurata)
per misurare un errore (rispetto alla misura della trasformazione igm e
non rispetto alla realtà) la possono dare solamente i soliti grigliati
igm. Con gli altri metodi si ottiene l'errore rispetto ad una fonte
dotata di una certa accuratezza.
L'argomento è lungo, nè io sono un super esperto, ma questi
considerazioni ritengo siano corrette. Ad ogni modo fatemi sapere se ho
scritto qualche castroneria.
Ciao
Cristoforo
-------- Messaggio Originale --------
Oggetto: Re:[Gfoss] precisione reproject WGS84 UTM33N -> Gauss-Boaga
Da: Luca Manganelli <[email protected]>
A: [email protected]
Data: mercoledì 28 gennaio 2009 13.55.39
----- "Luca Delucchi" <[email protected]> ha scritto:
Da: "Luca Delucchi" <[email protected]>
A: "Enrico Lambertini" <[email protected]>
Cc: [email protected]
Inviato: Mercoledì, 28 gennaio 2009 13:31:28 GMT +01:00
Amsterdam/Berlino/Berna/Roma/Stoccolma/Vienna
Oggetto: Re: [Gfoss] precisione reproject WGS84 UTM33N -> Gauss-Boaga
se la ctr è stata riproiettata con arcgis ti consiglio di non perdere
troppo tempo; con software libero e open source non si riesce (almeno
noi non siamo riusciti, perdendo qualche settimana) a eliminare
l'errore; l'unica soluzione è andare da qualcuno che ha arcgis e fare
li la riproiezione
Anche noi per il Comune di Trento abbiamo lavorato alla riproiezione da GB.
In particolare, prima con ArcGIS avevamo fatto uno script in python che faceva
le seguente riproiezione:
Gauss Boaga -> WGS84 -> ETRF89
nota: abbiamo scoperto che ETRF89 è identico ad ETRS89, e lo conferma pure
Proj.4 (vedi i file .csv che si trovano in /usr/share/proj)
Adesso con ogr2ogr facciamo la riproiezione da:
Gauss Boaga -> WGS84 -> ETRS89 (non ETRF89)
però abbiamo visto che WGS84 e ETRF/S89 sono praticamente... identici, cambia
solo il file .prj creato.
L'errore massimo che abbiamo trovato è di mezzo millimetro.
_______________________________________________
Iscriviti all'associazione GFOSS.it: http://www.gfoss.it/drupal/iscrizione
[email protected]
http://www.faunalia.com/cgi-bin/mailman/listinfo/gfoss
Questa e' una lista di discussione pubblica aperta a tutti.
I messaggi di questa lista non rispecchiano necessariamente
le posizioni dell'Associazione GFOSS.it.
_______________________________________________
Iscriviti all'associazione GFOSS.it: http://www.gfoss.it/drupal/iscrizione
[email protected]
http://www.faunalia.com/cgi-bin/mailman/listinfo/gfoss
Questa e' una lista di discussione pubblica aperta a tutti.
I messaggi di questa lista non rispecchiano necessariamente
le posizioni dell'Associazione GFOSS.it.
