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.
_______________________________________________ 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.
