Differences

This shows you the differences between two versions of the page.

--- tecnica:gps_cartografia_gis:shp2gpx [2008/02/01 00:54] – external edit 127.0.0.1
+++ tecnica:gps_cartografia_gis:shp2gpx [2009/03/20 18:51] (current) – niccolo
@@ Line 8: / Line 8: @@
 Esiste il programma **[[http://gpx2shp.sourceforge.jp/|gpx2shp]]**, ma contrariamente a quelle che si legge sulla home page non consente la trasformazione inversa (**shp2gpx**).
 ===== Da riga di comando =====
@@ Line 49: / Line 48: @@
   - Esportare il layer GRASS in formato shapefile.
   - Usare GPSBabel per convertire da shapefile a GPX.
 ==== Dettaglio ====
@@ Line 86: / Line 83: @@
 v.out.ogr input=AO2k_line type=line dsn=AO2k_line_shp olayer=AO2k_line_shp layer=1 format=ESRI_Shapefile
 </code>
+====== Conversione di coordinate ======
+Per una semplice conversione di coordinate si può utilizzare il comando **''cs2cs''** fornito dal pacchetto **proj**.
+Questo esempio converte una coppia di coordinate dal sistema **EPSG:3003** (Monte Mario, Italy zone 1 - Roma40) al sistema **EPSG:4326** (WGS84). La cosa è resa complicata perché la definizione di EPSG:3003 fornita da Proj non include i parametri di trasformazione verso WGS84, che vanno quindi aggiunti a mano:
+<code>
+echo "1674257 4846294" | cs2cs \
+    +init=epsg:3003 \
+    +towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68 \
+    +to +init=epsg:4326 -f "%.8f"
+.16398206     43.74906364 45.83517126
+</code>
+Il comando restituisce una terna di numeri (le coordinate **//x//**, **//y//** e **//z//**), cioè anche l'altezza.