[Gfoss] ancora: QGIS UTM ED50 Emilia Romagna

[Enrico Lambertini]

grazie a tutti per l'aiuto, soprattutto ad Antonio per la descrizione
attenta per dettagliata

Il giorno 21 settembre 2008 19.48, Antonio Falciano ha scritto:

>
> Ti visualizza le cordinate "originali" poichè riconosce solo questa
> stringa:
> +proj=utm +zone=32 +ellps=intl +units=m
> che di fatto è l'UTM ED50 32N (EPSG:23032).
> La sottostringa "+toWGS84 +y0=-4000000" non è corretta, perchè +towgs84
> è seguito da 3 o da 7 parametri e il falso nord si definisce con +y0 ma
> con +y_0.
> Inoltre, anche se il falso nord fosse stato definito correttamente, le
> coordinate rimarrebbero ugualmente inalterate poichè in proj.4 alcuni
> sistemi di coordinate, tra i quali quelli che usano "+proj=utm" e quindi
> anche UTM ED50 32N, contengono implicitamente i valori delle false
> origini (+x_0 e +y_0).
> Per bypassare il problema, si tratta di utilizzare un trick, in modo da
> "convincere" proj.4 a fare quello che vogliamo. Basta sapere cos'è UTM.
> E' una proiezione universale trasversa di Mercatore e quindi occorre
> semplicemente sostituire "+proj=tmerc" a "+proj=utm" per ovviare al
> problema delle false origini. Così facendo "+y_0=-4000000" dovrebbe
> avere effetto.
> Ora però, avendo sostituito la proiezione (anche se, di fatto, è la
> stessa!), dobbiamo specificare tutti gli altri parametri (ellissoide,
> centro della proiezione, false origini e fattore di contrazione), poichè
> proj.4 non sa più che stiamo parlando di una proiezione UTM. E quindi,
> la stringa per 23032* dovrebbe essere:
>
> +proj=tmerc +ellps=intl +lat_0=0 +lon_0=9 +x_0=500000 +y_0=-4000000
> +k=0.9996 +units=m
>
> In definitiva, per visualizzare correttamente le coordinate (senza falso
> nord) dovresti definire 23032 per la vista, abilitare la proiezione al
> volo ed assegnare 23032* (UTM-RER) ai dati.
>
> provato, FUNZIONA... ma solo coi vettoriali, se assegno questa proiezione a
un layer raster me la colloca correttamente, ma non la visualizza (bug di
qgis??)

>
> perchè sottrarre 4000000 alle y? sono coordinate in WGS84, lasciale così
> come sono... se la vista è in 23032, dovresti aggiungere


a questo punto penso che mi convenga tenermi come sistema di riferimento
della vista UTM-RER e crearmi una proiezione personalizzata con cui traslare
solo il layer dei punti WGS84 (UTM 32N), che se ho capito il concetto
dovrebbe essere qualcosa del tipo
+proj=tmerc +ellps=intl +lat_0=0 +lon_0=9 +x_0=-500000 +y_0=4000000
+k=0.9996 +units=m

giusto?

grazie

e.l.
-------------- parte successiva --------------
Un allegato HTML è stato rimosso...
URL: <http://www.faunalia.com/pipermail/gfoss/attachments/20080922/ce86cb19/attachment.htm>