Ho una tabella con una decina di punti le cui coordinate sono espresse in gradi sessagesimali. Per far meglio digerire questi dati a QGIS voglio (devo) trasformarli in sessadecimali. A tal fine uso la formuletta classica: Gradi sessadecimali = g° +(m’/60) + (m’‘/3600) dove g° sono i gradi, m’ i minuti primi e m’’ i minuti secondi.
Ora, se mi capitano coordinate ben divisibili per 3600 mi viene fuori un numero finito (tipo: 30°43’12’’ = 30°.72) e va tutto bene ma, se mi capitano coordinate mal divisibili per 3600 mi vengono fuori numeri con più di venti cifre decimali (tipo: 30°43’13’’ = 30,720277777777777777777777777778.
La domanda è questa. Per avere una coordinata sessadecimale che mi esprima, come precisione, il centimetro, a quante cifre decimale posso permettermi di poter arrotondare il numero che ottengo dalla formuletta di cui sopra?
Ciao, 6 cifre decimali dovrebbero fornire una precisione centimetrica.
-----
Dott. For. Ludovico Frate, PhD
--
Sent from: http://osgeo-org.1560.x6.nabble.com/QGIS-Italian-User-f5250612.html
Il 21/05/2020 22:46, Marco Spaziani ha scritto:
Ho una tabella con una decina di punti le cui coordinate sono espresse in gradi sessagesimali. Per far meglio digerire questi dati a QGIS voglio (devo) trasformarli in sessadecimali. A tal fine uso la formuletta classica: Gradi sessadecimali = g° +(m'/60) + (m''/3600) dove g° sono i gradi, m' i minuti primi e m'' i minuti secondi.
Ora, se mi capitano coordinate ben divisibili per 3600 mi viene fuori un numero finito (tipo: 30°43'12'' = 30°.72) e va tutto bene ma, se mi capitano coordinate mal divisibili per 3600 mi vengono fuori numeri con più di venti cifre decimali (tipo: 30°43'13'' = 30,720277777777777777777777777778.
La domanda è questa. Per avere una coordinata sessadecimale che mi esprima, come precisione, il centimetro, a quante cifre decimale posso permettermi di poter arrotondare il numero che ottengo dalla formuletta di cui sopra?
- poniamo sferica la Terra con raggio di 6378 km
- l'arco sulla circonferenza è dato dal prodotto del raggio per l'angolo al centro espresso in radianti
- l'angolo al centro è dato dal rapporto tra l'arco e ed il raggio
Quindi:
- l'arco di meridiano corrispondente ad un secondo sessagesimale è: 1/3600/180*3.14*6378000= 30 metri circa.
- larco di longidudine: 30*cos lat (per esempio alle nostre latitudini viene circa 20 metri)
questo vuol dire che se il tuo dato è espresso al secondo sessadecimale stai dando delle geometrie "a rettangoli" di 30x20 metri
Detto questo, calcoliamo i decimali per i gradi sessadecinali
angolo = 30 metri / 6378000 / 3.14 * 180 = 2.69*e-4 cioè devi scrivere almeno quattro decimali
Se il tuo dato ha un'accuratezza al centimetro:
angolo = 0,01 / 6378000 / 3.14 * 180 = 8.98e-8 gradi sessadecimali, almeno 8 decimali
Ma tieni presente l'assunto originale: se la tabella dei punti di partenza è al secondo sessagesimale, il centimetro è pura invenzione matematica. stai sui 30 metri!
marco
--
Marco Guiducci - 055 4383194
SITA - Sistema informativo territoriale e ambientale
Regione Toscana - Via di Novoli 26 - 50127 Firenze
On 5/22/20, Marco Guiducci <marco.guiducci@regione.toscana.it> wrote:
Il 21/05/2020 22:46, Marco Spaziani ha scritto:
Ho una tabella con una decina di punti le cui coordinate sono espresse
in gradi sessagesimali..........
ciao Marco,
scusa la franchezza, ma credo che così il problema sia mal posto e
impossibile da risolvere; provo a darne una esemplificazione diversa;
sull'equatore (della Terra, sul sole sarebbe diverso, ma credo la
prima basti per il tuo problema 
) un grado corrisponde a
2.pi.R/360 = 111317m (meglio ricontrollare i conti), quindi per avere
precisione centimetrica devi dividere per 1E7; la tua precisione
angolare deve essere di 1E-7 gradi per discriminare il cm sulla faccia
della Terra (stiamo parlando di angoli decimali)
(Detto non per asserire, ma per contribuire a dibattere )
ciao,
giuliano
PS: ho appena visto che Marco Guiducci ha dato già la risposta, ma
aggiungo questa, forse più semplice
Il 22/05/2020 09:33, Ludovico ha scritto:
Ciao, 6 cifre decimali dovrebbero fornire una precisione centimetrica.
facciamo questa prova: fare un file txt così
x y
12 43
12.000001 43
12.000001 43.000001
12 43.000001
caricarlo in QGis con epsg 4326
poi impostare la proiezione al volo su 3003
misurare la distanza tra i punti: 8 cm in longitudine e 11 in latitudine.
Con altri due decimali verrebbe il mm, se i dati sono al centimetro forse è corretto tenere il mm. Nel caso abbassare a 7 (in questo caso si ottiene 8 mm in long e 11 mm in lat).
marco
--
Marco Guiducci - 055 4383194
SITA - Sistema informativo territoriale e ambientale
Regione Toscana - Via di Novoli 26 - 50127 Firenze
Il 22/05/2020 11:45, Giuliano Curti ha scritto:
On 5/22/20, Marco Guiducci <marco.guiducci@regione.toscana.it> wrote:
Il 21/05/2020 22:46, Marco Spaziani ha scritto:
Ho una tabella con una decina di punti le cui coordinate sono espresse
in gradi sessagesimali..........ciao Marco,
scusa la franchezza, ma credo che così il problema sia mal posto e
impossibile da risolvere; provo a darne una esemplificazione diversa;sull'equatore (della Terra, sul sole sarebbe diverso, ma credo la
prima basti per il tuo problema
2.pi.R/360 = 111317m (meglio ricontrollare i conti), quindi per avere
precisione centimetrica devi dividere per 1E7; la tua precisione
angolare deve essere di 1E-7 gradi per discriminare il cm sulla faccia
della Terra (stiamo parlando di angoli decimali)(Detto non per asserire, ma per contribuire a dibattere
ciao,
giuliano
grazie del tuo contributo. rileggendo la mia, ho notato una imprecisione (per non dire di peggio :-)): "il tuo dato è espresso al secondo sessadecimale".
Questo argomento è stato trattato più volte qui in lista.
Più volte, negli anni, ho visto cose raccapriccianti, spesso dovute alla bizzarria di un foglio di calcolo (non faccio nomi) che nel momento di esportare (per esempio in un file tipo csv) troncava i numeri in funzione della larghezza della cella, con buona pace dei nostri carissimi decimali. Ma, risottolineo, purché quei decimali siano realmente significativi, cioè rispondano alla reale precisione del dato da trattare.
marco
--
Marco Guiducci - 055 4383194
SITA - Sistema informativo territoriale e ambientale
Regione Toscana - Via di Novoli 26 - 50127 Firenze
Grazie a tutti.
Alla prossima
Il giorno ven 22 mag 2020 alle ore 12:37 Marco Guiducci <marco.guiducci@regione.toscana.it> ha scritto:
Il 22/05/2020 11:45, Giuliano Curti ha scritto:
On 5/22/20, Marco Guiducci <marco.guiducci@regione.toscana.it> wrote:
Il 21/05/2020 22:46, Marco Spaziani ha scritto:
Ho una tabella con una decina di punti le cui coordinate sono espresse
in gradi sessagesimali…
ciao Marco,scusa la franchezza, ma credo che così il problema sia mal posto e
impossibile da risolvere; provo a darne una esemplificazione diversa;sull’equatore (della Terra, sul sole sarebbe diverso, ma credo la
prima basti per il tuo problema
2.pi.R/360 = 111317m (meglio ricontrollare i conti), quindi per avere
precisione centimetrica devi dividere per 1E7; la tua precisione
angolare deve essere di 1E-7 gradi per discriminare il cm sulla faccia
della Terra (stiamo parlando di angoli decimali)(Detto non per asserire, ma per contribuire a dibattere
ciao,
giulianograzie del tuo contributo. rileggendo la mia, ho notato una imprecisione
(per non dire di peggio :-)): “il tuo dato è espresso al secondo
sessadecimale”.
Questo argomento è stato trattato più volte qui in lista.
Più volte, negli anni, ho visto cose raccapriccianti, spesso dovute alla
bizzarria di un foglio di calcolo (non faccio nomi) che nel momento di
esportare (per esempio in un file tipo csv) troncava i numeri in
funzione della larghezza della cella, con buona pace dei nostri
carissimi decimali. Ma, risottolineo, purché quei decimali siano
realmente significativi, cioè rispondano alla reale precisione del dato
da trattare.
marco–
Marco Guiducci - 055 4383194
SITA - Sistema informativo territoriale e ambientale
Regione Toscana - Via di Novoli 26 - 50127 Firenze
QGIS-it-user mailing list
QGIS-it-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/qgis-it-user