[Gfoss] Assegnare proiezione in gvSIG

Pietro Blu Giandonato p.giandonato at gmail.com
Tue Sep 11 13:05:30 CEST 2007 

Il 11/09/07, Antonio Falciano ha scritto:
>
> Pietro Blu Giandonato ha scritto:
>
> > Qui appresso i parametri delle trasformazioni tra:
> > i) CRS 2) e 3), tipo: Geocentric Translation, parametri: dx = -104.0; dy
> > = -101.0; dz = -140.0
>
> ED_1950_To_WGS_1984_1 1133 Geocentric_Translation -87 -98 -121


Evidentemente ho sbagliato a riportare i parametri, perchè è proprio la
ED_1950_To_WGS_1984_1 che uso.

> ii) CRS 1) e 3), tipo: Position Vector, parametri: dx = -104.0; dy =
> > -49.1; dz = 9.9; rx = 0.971; ry = -2.917; rz = 0.714 ; s = -11.68
>
> Monte_Mario_To_WGS_1984_4 1660 Position_Vector -104.1 -49.1 -9.9 0.971
> -2.917 0.714 -11.68


E' lei. La utilizzo appunto per riproiettare GB in  WGS84 e viceversa.

> iii) CRS 1) e 2) non è possibile definire una trasformazione.
>
> Dovresti utilizzare "in serie" le due citate in precedenza mediante
> project tool.


Evito se posso di riproiettare i dati, preferisco usare la "on the fly" in
modo da lasciarli intonsi il più possibile...

> Appare chiaro dunque che la mancata sovrapposizione tra i layer 1) e 2)
> > settando uno o l'altro dei CRS è dovuta al fatto che manca una
> > trasformazione tra i due sistemi di riferimento. Il problema invece non
> > sussiste se scegliamo il WGS84 UTM33N come sistema di rappresentazione
> > dell'intera mappa, poichè 1) e 2) possono essere trasformati nel 3).
>
> Infatti, le trasformazioni "composte" di coordinate ricorrono ad datum
> geocentrico di appoggio (WGS84).


Buono a sapersi.


> > Molto simile a quello che nelle versioni seguenti di ArcView 3.x sarebbe
> > poi stato il GB, finalmente supportato, ma in cui falso il est è
> > 2.520.000 metri e il falso nord 0 metri. Questa piccola differenza,
> > usando il GB "ufficiale", mi portava a riproiezioni di layer che
> > generavano uno shift diagonale che andava da 50 a oltre 200 metri,
> > mentre con il "mio" GB lo shift era quasi inesistente.
>
Al momento in ambito
> GFOSS, GRASS e Qgis che utilizzano le librerie PROJ consentono di
> applicare le trasformazioni aggiungendo al WKT di proj il parametro
> +towgs84, ovvero le traslazioni, le rotazioni ed un fattore di scala
> (questi ultimi sono presenti solo nelle trasformazioni a 7 parametri tipo
> Position Vector e Coordinate Frame).


 Lasciamo perdere allora... di fatto GB e ED50 UTM33N sono "incompatibili" a
meno di non effettuare riproiezioni in serie.

C'è anche da dire che le trasformazioni dell'EPSG presentano accuratezze
> fino a 4 m e quindi per grandi scale a partire dal 25k lasciano il tempo
> che trovano...


Sempre meglio che niente. Per il lavoro che faccio mi interessa avere il
massimo possibile di accuratezza dal p. di vista spaziale, l'importante è
essere sempre in grado di tenere traccia delle trasformazioni applicate ai
dati.

Comunque, tirando le fila del discorso, di fronte al fatto che applicativi
FOSS capaci di gestire riproiezioni al volo di più raster con CRS diversi
pare non ve ne siano, l'unica strada possibile è la riproiezione materiale
dei dati, con tutte le conseguenze del caso.

Con ArcGIS è invece possibile farlo, ma:
1. impostando il CRS del progetto/vista unicamente in WGS84, poichè GB e
ED50 convergono correttamente verso WGS84, con errore di posizionemento di
circa 1-2 metri con dati come ortofoto con risoluzione pixel di 50cm;
2. o usando un CRS per GB "corretto" (come quello ottenuto da me), qui
l'errore è più alto, aggirandosi attorno ai 4-5 metri sempre con dette
ortofoto. L'unico vantaggio è la possibilità di poter scegliere
indifferentemente uno qualunque dei 3 CRS per il progetto/vista, i dati si
sovrapporranno sempre tra loro, tenendo conto degli errori ovviamente.

Saluti e grazie per le info,
PB



>
> >
> >
> >
> > Il 10/09/07, *Markus Neteler*  ha scritto:
> >
> >     On Sun, Sep 09, 2007 at 06:44:57AM -0700, mando wrote:
> >      >
> >      > Salve,
> >      > ho iniziato a dare un'occhiata a gvSIG, ma non riesco a capire
> come
> >      > assegnare ad un layer la proiezione corretta.
> >      > Ho una base con raster 1:25000 e 1:5000 della provincia di rimini
> >     e con QGIS
> >      > li carico correttamente con Monte Mario (Rome) Italy zone 2.
> >      >
> >      > Ora in QGIS carico anche i confini ISTAT dei comuni italiani che
> >     sono in
> >      > ELD79/UTM zone 32 e tramite proiezione al volo vengono
> correttamente
> >      > sovrapposti ai dati sopra descritti.
> >      >
> >      > In gvSIG non riesco a fare questo. Quale proiezione devo
> >     assegnare o meglio,
> >      > esistono delle tabelle di confronto dei vari nomi e codici (EPSG,
> >     postgis,
> >      > etc) in modo da poter sapere quale proiezione usare?
> >
> >     Credo che sia il solito problema che il software non contiene
> >     i datum geodetici di GaussBoaga. Non ho verificato con gvSIG pero'.
> >
> >     Missa che GRASS e' l'unico GIS libero che contiene questi
> >     parametri :) Infatti, e' un po triste che tanti software hanno
> >     ancora questo problema.
> >
> >     Se mancano i datum italiani in gvSIG, bisogna definire la proiezione
>
> >     "a mano" includendo il parametro "towgs84...". Nell'archivio di
> questa
> >     lista dovrebbono esserci questi parametri (si trovano anche GRASS o
> >     sul sito http://crs.bkg.bund.de/crs-eu/
> >     <http://crs.bkg.bund.de/crs-eu/>   -> CRS Description -> national ->
> I).
> >
> >     ciao
> >     Markus
> >
> >     ------------------
> >     ITC -> dall'1 marzo 2007 Fondazione Bruno Kessler
> >     ITC -> since 1 March 2007 Fondazione Bruno Kessler
> >     ------------------
> >
> >     _______________________________________________
> >     Iscriviti all'associazione GFOSS.it:
> >     http://www.gfoss.it/drupal/iscrizione
> >     Gfoss a faunalia.com <mailto:Gfoss a faunalia.com>
> >     http://www.faunalia.com/cgi-bin/mailman/listinfo/gfoss
> >
> >
> >
> > ------------------------------------------------------------------------
> >
> > _______________________________________________
> > Iscriviti all'associazione GFOSS.it:
> http://www.gfoss.it/drupal/iscrizione
> > Gfoss a faunalia.com
> > http://www.faunalia.com/cgi-bin/mailman/listinfo/gfoss
> >
> >
> > ------------------------------------------------------------------------
> >
> > No virus found in this incoming message.
> > Checked by AVG Free Edition.
> > Version: 7.5.485 / Virus Database: 269.13.13/998 - Release Date:
> 10/09/2007 8.48
>
>
>
> _______________________________________________
> Iscriviti all'associazione GFOSS.it: http://www.gfoss.it/drupal/iscrizione
>
> Gfoss a faunalia.com
> http://www.faunalia.com/cgi-bin/mailman/listinfo/gfoss
>
-------------- parte successiva --------------
Un allegato HTML ? stato rimosso...
URL: http://www.faunalia.com/pipermail/gfoss/attachments/20070911/b460594f/attachment.htm

More information about the Gfoss mailing list