[Gfoss] Proiezioni

giuliano curti giulianc a tiscali.it
Sab 14 Feb 2009 00:53:04 CET


Ciao,
sto seguendo la lista da un po' di tempo perche' vorrei entrare
nel fascinoso mondo dei gis......

sto cercando di capire qualcosa delle trasformazioni e in
particolare di come sono calcolati i seguenti parametri che ho
tratto da fiduciali.it:

lat(WGS84)  lon(WGS84) N(UTM32)   E(UTM32)
45.3141298  9.50299354 5017970.25 539426.36 #Lodi
45.36041241 9.32232977 5023039.67 525244.74 #F100-05-01
45.3582976  9.3243868  5022805.15 525406.79 #F100-05-02
45.35935875 9.32609001 5022923.72 525539.72 #F100-05-03
45.35656502 9.32853526 5022613.99 525732.5  #F100-12-02
(si tratta di una origine su Lodi e di alcuni punti fiduciali del
mio comune);

ho fatto i test con proj (inserendo in ordine LON e LAT), ma
settando +ellps=WGS84 +proj=utm ottenevo dati non conformi con
quelli del prospetto;
 
fortunatamente ho rivisto la mail di
Mando del 10 feb u.s. ed ho provato i parametri da lui riportati:
effettivamente con <proj +proj=tmerc +lat_0=0.0 +lon_0=9.0
+k=0.9996 +x_0=500000.0 +unit=m> (ho omesso +y_0=-4000000.0
+ellps=intl) ottengo questi risultati:

539426.36       5017970.25 #Lodi
525244.74       5023039.45 #F100-05-01
525406.79       5022805.15 #F100-05-02
525539.72       5022923.58 #F100-05-03
525732.50       5022613.99 #F100-12-02

che coincidono con la tabella;

l'euforia per il risultato ottenuto non nasconde il fatto che mi
sarei aspettato di dover usare +ellps=WGS84 +proj=utm: dove
sbaglio?

grazie, ciao,
giuliano






Maggiori informazioni sulla lista Gfoss