sto sbattendo la testa con la riproiezione di alcune cartografie raster georeferenziate secondo wgs84 utm33N (epsg: 32633) che dovrei convertire in Gausss-Boaga fuso est (EPSG: 3004) per sovrapporle allla ctr.
Dopo alcuni infruttuosi tentativi con grass ho ottenuto un buon risultato con gdalwarp:
il problema è che mi permane un fastidiosissimo errore di spostamento di circa 70 mt. qualcuno sa dove sbaglio? o come eliminare quest’errore?
Altra domanda: dovendo scegliere un sitema di riferimento di base per i progetti, nel quale riproiettarre le cartografie derivanti da diverse fonti, qual’è il metodo di scelta / sistema di riferimeto che secondo voi conviene adottare?
sto sbattendo la testa con la riproiezione di alcune cartografie raster
georeferenziate secondo wgs84 utm33N (epsg: 32633) che dovrei convertire
in Gausss-Boaga fuso est (EPSG: 3004) per sovrapporle allla ctr.
Dopo alcuni infruttuosi tentativi con grass ho ottenuto un buon
risultato con gdalwarp:
il problema è che mi permane un fastidiosissimo errore di spostamento di
circa 70 mt. qualcuno sa dove sbaglio? o come eliminare quest'errore?
dovresti inserire i parametri +towgs84 della trasformazione:
-t_srs "+proj=tmerc +ellps=intl +lat_0=0 +lon_0=15 +k=0.999600
+x_0=2520000 +y_0=0 +units=m
+towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68"
in luogo di -t_srs EPSG:3004
Non lo elimina l'errore, ma lo dovrebbe ridurre in maniera consistente.
Altra domanda: dovendo scegliere un sitema di riferimento di base per i
progetti, nel quale riproiettarre le cartografie derivanti da diverse
fonti, qual'è il metodo di scelta / sistema di riferimeto che secondo
voi conviene adottare?
Tutte le strade portano a ...WGS84 e UTM derivati!
Il 28 gennaio 2009 13.12, Enrico Lambertini <lamba84@gmail.com> ha scritto:
Buongiorno a tutti,
sto sbattendo la testa con la riproiezione di alcune cartografie raster
georeferenziate secondo wgs84 utm33N (epsg: 32633) che dovrei convertire in
Gausss-Boaga fuso est (EPSG: 3004) per sovrapporle allla ctr.
Dopo alcuni infruttuosi tentativi con grass ho ottenuto un buon risultato
con gdalwarp:
il problema è che mi permane un fastidiosissimo errore di spostamento di
circa 70 mt. qualcuno sa dove sbaglio? o come eliminare quest'errore?
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
Altra domanda: dovendo scegliere un sitema di riferimento di base per i
progetti, nel quale riproiettarre le cartografie derivanti da diverse fonti,
qual'è il metodo di scelta / sistema di riferimeto che secondo voi conviene
adottare?
secondo me dipende cosa devi fare e quanti dati hai da
riproiettare...comunque ormai i più utilizzati credo che siano WGS84 e
UTM ED50
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
Questo è *assolutamente* un mito da sfatare! ArcGIS usa gli stessi
parametri che ho scritto poco fa. E poi anche con software libero lo
puoi fare in meno di un minuto.
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.
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.
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'errore è relativo al confronto dei shapefile riproiettati da ArcGIS e da OGR2OGR da GB a WGS84, partendo da un unico shapefile in Gauss Boaga.
Cosa è che ti lascia perplesso?
Ho misurato l'errore con lo strumento righello ed è lo stesso in ArcGIS Desktop, in QGIS e gvSig.
On Wed, Jan 28, 2009 at 1:28 PM, Antonio Falciano <afalciano@yahoo.it> wrote:
Enrico Lambertini ha scritto:
Buongiorno a tutti,
sto sbattendo la testa con la riproiezione di alcune cartografie raster
georeferenziate secondo wgs84 utm33N (epsg: 32633) che dovrei convertire
in Gausss-Boaga fuso est (EPSG: 3004) per sovrapporle allla ctr.
Dopo alcuni infruttuosi tentativi con grass ho ottenuto un buon
risultato con gdalwarp:
in luogo di -t_srs EPSG:3004
Non lo elimina l’errore, ma lo dovrebbe ridurre in maniera consistente.
era proprio quello che mi servia favoloso (grazie anche a Maurizio)
Altra domanda: dovendo scegliere un sitema di riferimento di base per i
progetti, nel quale riproiettarre le cartografie derivanti da diverse
fonti, qual’è il metodo di scelta / sistema di riferimeto che secondo
voi conviene adottare?
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.
visto che stiamo lavorando proprio con questi dati mi daresti il
codice che usate? noi con i valori passati dalla provincia abbiamo un
errore di 100 m. Per di più il mio collega ha notato grossissimi
errori (anche decine di km) sui punti fiduciali catastali che sono poi
il punto di partenza...
On Wed, Jan 28, 2009 at 3:26 PM, Luca Delucchi <lucadeluge@gmail.com> wrote:
visto che stiamo lavorando proprio con questi dati mi daresti il
codice che usate? noi con i valori passati dalla provincia abbiamo un
errore di 100 m. Per di più il mio collega ha notato grossissimi
errori (anche decine di km) sui punti fiduciali catastali che sono poi
il punto di partenza...
Nella stessa provincia non avevo questi problemi
E probabilmente solo un problema di impostazione programma che hai...
Nella stessa provincia non avevo questi problemi
E probabilmente solo un problema di impostazione programma che hai...
ho provato con tre programmi diversi (grass, gvsig, ogr2ogr), passando
i parametri direttamente oppure creando nel file epsg un codice nuovo,
ma l'errore era sempre lo stesso!
favoloso (grazie anche a Maurizio)