Il 15/05/2018 11:05, a.furi...@lqt.it ha scritto:
On Tue, 15 May 2018 10:41:16 +0200, Giuliano Curti wrote:
una richiesta se posso: non sono riuscito a capire dove la nuova
definizione documenta i parametri:
...... +towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68 \
ciao Giuliano,
se per "nuova definizione" intendi quella WKT la risposta
e' molto semplice.
- la matrice a 7 parametri "+towgs94" e' esattamente uno
di quegli elementi "fuori standard" esclusivi della PROJ
che si intendono eliminare.
n.b.: quando effettui una riproiezione p.es. dal 3003 al
32632 la PROJ in effetti applica due trasformazioni:
- la prima dal 3003 al 4326 (WGS84)
- la seconda dal 4326 al 32632
ecco perche' la PROJ usa a man bassa il +towgs84; perche'
le e' strettamente indispensabile per riproiettare dal
piano all'ellissoide e viceversa.
ma purtroppo "doppia trasformazione" significa anche
"doppio errore" (arrotondamenti, troncamenti etc).
- l'approccio WKT e' del tutto differente.
come si capisce senza troppa fatica, vengono definiti
in modo minuzioso l'ellissoide di riferimento, il datum,
il meridiano fondamentale, l'orientazione degli assi,
eventuali "false easting" e "false northing", le unita'
di misura, fattori di scala etc
e quidi, essendo noti tutti i parametri su entrambi i
lati della trasformazione, non e' piu' richesta la
"doppia trasformazione"; puoi andare sparato da un
CRS all'altro in modo diretto, a tutto vantaggio della
precisione dei calcoli.
- ma non finisce neppure qua: leggo che tra le novita'
proposte per la nuova PROJ c'e' pure un meccanismo
abbastanza raffinato che dovrebbe verificare dove
cade esattamente l'area della geometria da trasformare
rispetto ai BBOXes dei due CRS, in modo tale da potere
compensare automaticamente quei fattori di distorsione
che si verificano inevitabilmente quando ci si
allontano in modo sensibile dal centro del fuso
di riferimento ... almeno sulla carta, decisamente
notevole ;-)
ciao Sandro
ciao Sandro, quel che che scrivi è giusto ma forse incompleto.
La definizione WKT che riporti nel post precedente definisce il datum e
basta, non definisce la relazione che quel datum ha con altri datum,
come invece fa la stringa dei 7 parametri.
Con la wkt che mostri tu si hanno solo le informazioni dell'ellissoide e
della proiezione associata: questi dati bastano per passare da
coordinate ellissoidiche e coordinate piane e viceversa (le formule di
gauss sono in letteratura). Mi sembra che nella definizione completa
della wkt sia ancora presente il parametro towgs84.
In parole povere: definendo i due datum non siamo a niente. Bisogna
definire in che relazione stanno se abbiamo bisogno di passare ad altri
datum, i quali non hanno una relazione matematica tra di loro. Il
passaggio si fa solo in modo statistico.
I 7 parametri realizzano un sistema di passaggio "medio" tra i sistemi.
Senza quelli non si va da nessuna parte. Quindi credo che in realtà un
meccanismo simile sia implementato anche nella nuova versione della
libreria.
Come sappiamo in Italia l'IGM ha realizzato, attraverso la conoscenza
delle coordinate di alcuni (molti) punti nei sistemi "nostrani" (ED50,
ETRF89, ETRF2000, ROMA40) un meccanismo di passaggio più raffinato dei 7
parametri, che ci permette di ottenere scarti minori perché calibrato
sulla relazione delle realizzazzioni delle reti trigonometriche
esistenti, dove per per realizzazione si intende materializzare sul
terreno un sistema di riferimento.
I famosi "grigliati", scritti in un formato ascii definito da IGM,
esistono anche nel formato NTv2 forse più idoneo ad essere implementato
in software "general purpose".
Fai riferimento al BBOX ed alla compensazione di distorsioni locali.
Bene, ma rimane sempre da capire che tipo di informazioni di maggior
dettaglio saranno contenute per distinte zone geografiche. Da dove
provengono? Chi le fornisce?
Sarà utile, al momento giusto, sfruttare e integrare le capacità e le
competenze di chi sa leggere il codice, da una parte, e di chi si
intende si sistemi di riferimento dall'altra per capire l'effettiva
valenza e potenza della libreria.
saluti
marco
--
Marco Guiducci - 055 4383194
SITA - Sistema informativo territoriale e ambientale
Regione Toscana - Via di Novoli 26 - 50127 Firenze
_______________________________________________
Gfoss@lists.gfoss.it
http://lists.gfoss.it/cgi-bin/mailman/listinfo/gfoss
Questa e' una lista di discussione pubblica aperta a tutti.
I messaggi di questa lista non hanno relazione diretta con le posizioni
dell'Associazione GFOSS.it.
796 iscritti al 28/12/2017
