[Gfoss] Proiezioni

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

On Sat, 14 Feb 2009 00:53:04 +0100
giuliano curti <giulianc@tiscali.it> wrote:

Ciao,
....

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

che coincidono con la tabella;

anche questi danno lo stesso risultato:
<proj +proj=utm +zone=32 +ellps=WGS84 +datum=WGS84 +unit=m
+no_defs>

ciao,
giuliano

Ciao,

2009/2/14 giuliano curti <giulianc@tiscali.it>:
...

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

...

Non sto seguendo in dettaglio la discussione ma perché
non usate "cs2cs"? "proj" non fa trasformazione del datum
geodetico...

con
cs2cs -v ...
si vede anche che parametri vengono usati veramente (-v:
verbose).

ciao
Markus