Unità di misura angolari

Proprio oggi pomeriggio, il mio strumentista, dopo avere effettuato una stazione libera, mi fa presente alla radio: abbiamo un errore angolare di 2 cm.
Io mi sono chiesto: che cos'è un errore angolare di 2 cm? Gli angoli non si misurano in cm, ma in gradi . Per cui, che cosa avrà voluto dire?
Un grado sessagesimale corrisponde alla 360esima parte di un angolo giro, e un angolo giro è 2 pigreco in radianti. il valore in radianti, moltiplicato per una distanza in metri, mi fornisce lo scostamento trasversale in metri (facilmente convertibile in centimetri)
Ammettiamo che la stazione libera in non fosse delle migliori, i due caposaldi erano distanti tra loro 100 metri e circa 100 metri dalla stazione.
Quanto è, in gradi sessagesimali, un errore angolare di due centimetri?

Il conto è piuttosto semplice:

0.02/100 è il valore in radianti, pari a 0.0002

Facciamo una proporzione

x/2pigreco = y/360

dove x è il valore in radianti, mentre y è il valore in gradi sessagesimali

quindi y = (x * 360) /2pigreco = 0.011459

a   0.011459  corrispondono 0  0.011459  * 3600 secondi = 41 secondi

Non male per uno strumento al secondo..

Un valore abbastanza comune, nell'orientamento, è di 5mm su di una distanza di 200 metri.

In questo caso, l'errore è di circa 5 secondi.

In pratica, l'errore di un secondo corrisponde a 1 mm di precisione alla distanza di 200 metri.

Poligonali o meglio : Reti

Prendo in esempio Civil Design, in quanto è il software che conosco meglio.
Lo stesso identico ragionamento può essere effettuato con Sierra, Leonardo, Meridiana, TGO, TBC, LeicaGoeoffice, Carazzai, Tabula, ed anche Pregeo.

Prendendo il file di esempio, nella sottocartella esempi di civil design, c'è una poligonale per strati, che può essere calcolata in molti, moltissimi modi.

Citando il manuale:
- nessuna compensazione: i dati del rilievo verranno elaborati senza nessuna compensazione (nel caso di calcolo 3D non vengono applicate correzioni di curvature terrestre);

- compensazione ai minimi quadrati vincolata: si usa in genere la compensazione ai minimi quadrati quando le poligonali formano una rete complessa e si hanno misure eccedenti rispetto a quelle necessarie. Tutte le misure effettuate (angolari e distanze) rientrano nei calcoli con un peso inversamente proporzionale allo scarto quadratico medio. Lo sqm può essere definito a priori associando alla stazione (dialogo “Tabella Dati”, pulsante ) uno strumento definito precedentemente, in cui sono definiti i valori per ogni tipo di lettura (azimutale, zenitale e distanza). Nel caso di misura iterata (strati di lettura), il programma calcola in automatico media e sqm delle letture ripetute, sostituendolo allo sqm a priori. Il calcolo viene effettuato in modo che , dove ‘s’ sono gli scarti residui alle letture. Nel caso particolare di compensazioni ai minimi quadrati vincolate nel calcolo vengono introdotte come termini noti le coordinate dei punti noti.

- rototraslazione sui punti di maggiore attendibilità: i punti e l’attendibilità vengono definiti nel dialogo Gestione punti noti: il rilievo, dopo essere stato compensato ai minimi quadrati internamente (quindi senza utilizzare coordinate di Punti noti), viene dapprima traslato sul punto più attendibile e quindi ruotato sul secondo punto più attendibile.

- rototraslazione ai minimi quadrati senza variazione di scala: dopo aver eseguito il calcolo con compensazione ai minimi quadrati senza l’uso dei punti noti viene eseguita una rototraslazione rigida in cui non si modifica la mutua posizione dei punti rototraslati (il baricentro del rilievo rimane inalterato) e si minimizza la distanza fra punti di rilievo calcolati e noti;

- rototraslazione ai minimi quadrati con variazione di scala isotropa o conforme: la variazione di scala ingrandisce o rimpicciolisce il rilievo, con il risultato di ridurre ulteriormente gli scarti residui fra punti calcolati e residui;

- rototraslazione ai minimi quadrati con variazione di scala anisotropa o differenziata: la variazione di scala è diversa in X, Y;

- compensazione empirica: ogni vertice viene traslato, di entità proporzionale allo sviluppo dei lati precedenti, nella stessa direzione e nello stesso senso dell’errore riscontrato.

 Ipotizziamo, in modo da ridurre le possibilità,  che i punti abbiano la medesima attendibilità, escludiamo anche la compensazione empirica e la variazione di scala anisotropa, e la nessuna compensazione.


In pratica, ci sono tre casi, di cui riporto il risultato del calcolo:

compensazione ai minimi quadrati vincolata
V100 3224.531 12293.662
V101 3311.788 12452.826
V102 3380.806 12610.046
V103 3473.772 12774.243
V104 3546.166 12938.506
V106 3610.158 13058.925
V108 3682.868 13148.988
V110 3800.831 13257.140
V111 3925.470 13361.292
V113 4027.717 13486.648
V114 4120.441 13640.513
V115 4225.596 13780.785
V116 4346.671 13976.673
V118 4442.427 14114.256
V120 4529.065 14203.915
V121 4649.753 14293.583
V123 4744.799 14394.817
V125 4839.838 14553.156
V126 4882.787 14701.559
V128 4937.742 14863.534
V130 5000.001 14999.999
V132 5114.689 15175.456
V133 5253.521 15359.984


rototraslazione ai minimi quadrati senza variazione di scala

V100 3224.521936 12293.643
V101 3311.776651 12452.807
V102 3380.794667 12610.031
V103 3473.760911 12774.231
V104 3546.154819 12938.496
V106 3610.147717 13058.918
V108 3682.859322 13148.982
V110 3800.824552 13257.136
V111 3925.465661 13361.289
V113 4027.71387 13486.646
V114 4120.439722 13640.513
V115 4225.597546 13780.787
V116 4346.67504 13976.676
V118 4442.432716 14114.261
V120 4529.072286 14203.922
V121 4649.762646 14293.590
V123 4744.810806 14394.826
V125 4839.850689 14553.168
V126 4882.800328 14701.573
V128 4937.755459 14863.550
V130 5000.015061 15000.017
V132 5114.697621 15175.471
V133 5253.527237 15360.001

rototraslazione ai minimi quadrati con variazione di scala isotropa o conforme
V100 3224.535 12293.663
V101 3311.788 12452.825
V102 3380.805 12610.047
V103 3473.770 12774.245
V104 3546.163 12938.507
V106 3610.155 13058.928
V108 3682.866 13148.991
V110 3800.830 13257.143
V111 3925.469 13361.295
V113 4027.716 13486.650
V114 4120.441 13640.515
V115 4225.597 13780.787
V116 4346.673 13976.674
V118 4442.430 14114.257
V120 4529.068 14203.917
V121 4649.757 14293.584
V123 4744.804 14394.818
V125 4839.842 14553.158
V126 4882.791 14701.561
V128 4937.746 14863.536
V130 5000.005 15000.001
V132 5114.686 15175.453
V133 5253.513 15359.980

I minimi quadrati e la rototraslazione con variazione di scala sono praticamente identici anche se, a mio avviso, è da preferire la seconda in quanto, in caso di errore sulle coordinate dei caposaldi, non inficierei la bontà delle misurate, ma solamente delle coordinate.

Quest'ultima caratteristica è molto simile alla calibrazione del GPS, nel caso si impieghi un fattore di scala adattato alle coordinate dei punti

Qualora le coordinate dei caposaldi non  fossero di attendibilità sufficiente, nel GPS si impiega un fattore di scala fisso, dato dalla quota sull'ellissoide di riferimento: in modo analogo, su Civil Design si può imporre la quota di progetto, dal menù imposta del calcolo planimetrico. In questo caso, andrà impiegata la rototraslazione senza variazione di scala.

L'altezza di progetto è, generalmente, l'altezza media del rilievo. E' da considerare che, se si usano le quote sul livello del mare, non andrà utilizzato WGS84 come ellissoide, nel caso dell'italia, ma Hayford 1909 (anche se sarebbe più corretto utilizzare Roma40, e chiamarlo Roma 40); in ogni caso va inserita la latitudine del luogo.

Coordinate Piane

Coordinate Piane, coordinate rettilinee, coordinate gauss boaga uffff!!!
Insomma, che coordinate dobbiamo usare!??!?!?!
Per quel che mi riguarda, e per quel che riguarda anche TGO, esistono solo due tipi di coordinate: le coordinate suolo e le coordinate reticolo.
Le coordinate suolo sono le coordinate giuste, ovvero:


Cioè la distanza su di una geodetica che rappresenti, o meglio approssimi, la superficie reale.

Poi ci sono le coordinate reticolo, ovvero le coordinate su di un piano di proiezione.

Le coordinate reticolo dipendono dall'uso che ne dobbiamo fare:
le coordinate gauss boaga, per esempio, servono per rappresentare l'italia, divisa su due fusi, alla scala di 1:25000, ovvero con la precisione di 5 metri (ricordate l'errore di graficismo?)
le coordinate cassini servono, ma che ve lo dico a fare?
Le coordinate rettilinee sono coordinate piane che approssimano al meglio le coordinate suolo. Non saranno mai perfette perchè un ellissoide non si può rappresentare sul piano senza deformazioni.

Per esempio, quando diciamo " Nessuna proiezione nessun datum" e diamo l'altezza di progetto pari a 300 metri, diciamo che l'ellissoide si "gonfia" fino a + 300 metri (circa 260 slm) per meglio rappresentare la superficie fisica.
Da un punto di tangenza, medio tra quelli di calibrazione, si sviluppa un piano mercatore che approssima la superficie terrestre.

Su TGO, possiamo scegliere il nostro punto di partenza:

Da proprietà progetto, impostazioni sito locale, se premiamo, possiamo scegliere il punto di tangenza e l'altezza di progetto



Vedete la scala? è diversa, anzi maggiore di 1, e infatti siamo a 983 metri sull'ellissoide.

Per verificare se il sistema di cooridinate è corretto, io faccio così:

Prendo i due punti estremi del rilievo, vado su rilevamento (o CoGo sul controller) e verifico che la distanza suolo e quella reticolo coincidano in modo ragionevole.

In questo caso, 3 cm su di un chilometro, per due punti che hanno una differenza di quota di 150 metri, possiamo anche accontentarci. Sono 22 ppm ma, per prova, provate a vedere che schifezza viene con gauss boaga!