# Problema di Hansen

Il problema di Hansen (o doppia intersezione inversa) è un problema topografico, dal nome dell'astronomo Peter Andreas Hansen (1795 – 1874), che ha lavorato al rilevamento geodetico della Danimarca. Consiste nell'avere due punti di coordinate note $A$ e $B$ e due punti ignoti, da cui si effettua una stazione libera, $P$ e $Q$. Da $P$ e $Q$ l'osservatore misura gli angoli formati dalle linee di mira rispetto a ciascuno degli altri tre punti. Il problema consiste nel trovare le posizioni della stazione $P$ e della stazione (ausiliaria) $Q$. Guarda la figura; gli angoli misurati sono ($\alpha _{1}$, $\beta _{1}$, $\alpha _{2}$, $\beta _{2}$)$_.$
Poiché comporta osservazioni di angoli formati da punti sconosciuti, il problema è un esempio di intersezione inversa.
Questo problema può essere risolto in tre modi: il metodo grafico, il metodo della base (o della figura) fittizia o il metodo dell'angolo ausiliario.

## Metodo risolutivo dell'angolo ausiliario

Il metodo dell'angolo ausiliario consiste nell'individuare un angolo fittizio, spesso chiamato $k$ o $\lambda$, che non ha rispondenza reale con gli elementi del rilievo topografico, ma che ci dà la possibilità di individuare gli angoli ignoti del problema, altrimenti irrisolvibili a causa della carenza di elementi noti del poligono.
Vengono definiti i seguenti angoli: $\gamma =P{\widehat {A}}Q$, $\delta =P{\widehat {B}}Q$, $\varphi =Q{\widehat {A}}B$, $\psi =P{\widehat {B}}A$.
Dobbiamo risolvere $\varphi$ e $\psi$. Iniziamo determinando la somma di questi due angoli sconosciuti, uguale alla somma di $\beta _{1}$ e $\beta _{2}$, ottenendo l'equazione
    $\varphi +\psi =\alpha _{2}+\beta _{1}$
In alternativa si può osservare il triangolo $PAB$ e ottenere l'equazione:
    $\varphi +\psi =\pi -\left(\alpha _{1}+\gamma \right)$
Questa equazione verrà successivamente messa a sistema. Possiamo stabilire un legame tra il lato $AB$ e $PQ$ grazie al Teorema dei seni. Osservando i triangoli $AQB$ e $PQB$, possiamo scrivere le seguenti uguaglianze:
    ${\frac {AB}{QB}}={\frac {\sin \beta _{2}}{\sin \varphi }}$
    ${\frac {QB}{PQ}}={\frac {\sin \alpha _{2}}{\sin \delta }}$
Moltiplicando questi rapporti, otteniamo:
    ${\frac {AB}{QB}}\cdot {\frac {QB}{BQ}}={\frac {AB}{PQ}}={\frac {\sin \alpha _{2}\sin \beta _{2}}{\sin \varphi \sin \delta }}$
Conduciamo un ragionamento analogo sui triangoli $ABP$ e $APQ$:
    ${\frac {AB}{AP}}={\frac {\sin \alpha _{1}}{\sin \psi }}$
    ${\frac {AP}{PQ}}={\frac {\sin \beta _{1}}{\sin \gamma }}$
Moltiplicando questi rapporti, otteniamo:
    ${\frac {AB}{AP}}\cdot {\frac {AP}{PQ}}={\frac {AB}{PQ}}={\frac {\sin \alpha _{1}\sin \beta _{1}}{\sin \psi \sin \gamma }}$
Eguagliando i due risultati si ottiene
    ${\frac {\sin \alpha _{2}\sin \beta _{2}}{\sin \varphi \sin \delta }}={\frac {\sin \alpha _{1}\sin \beta _{1}}{\sin \psi \sin \gamma }}$
Isolando ${\frac {\sin \varphi }{\sin \psi }}$ si ottiene:
    ${\frac {\sin \varphi }{\sin \psi }}={\frac {\sin \alpha _{2}\sin \beta _{2}\sin \gamma }{\sin \alpha _{1}\sin \beta _{1}\sin \delta }}$

### Determinazione dell'angolo ausiliario

L’ultima espressione che abbiamo ottenuto rappresenta un valore finito, compreso tra $-\infty$ e $\infty$, che non trova una corrispondenza nelle misure plano-altimetriche del rilievo. Possiamo dire con certezza però, che esiste un angolo la cui tangente assume questo valore, visto che l'intervallo $(-\infty ;+\infty )$ è compreso nel codominio della funzione tangente. Impostiamo quindi:
    ${\frac {\sin \varphi }{\sin \psi }}=\tan \lambda$
Utilizzando la proprietà del comporre e scomporre delle proporzioni, possiamo riscrivere l'uguaglianza precedente come:
    ${\frac {\sin \varphi -\sin \psi }{\sin \varphi +\sin \psi }}={\frac {\tan \lambda -1}{\tan \lambda +1}}$
Riscriviamo ora il primo membro utilizzando la prima e la seconda formula di prostaferesi e semplifichiamo.
    ${\frac {2\cos {\frac {\varphi +\psi }{2}}\sin {\frac {\varphi -\psi }{2}}}{2\sin {\frac {\varphi +\psi }{2}}\cos {\frac {\varphi -\psi }{2}}}}={\frac {\tan \lambda -1}{\tan \lambda +1}}$
    $\tan {\frac {\varphi -\psi }{2}}\cot {\frac {\varphi +\psi }{2}}={\frac {\tan \lambda -1}{\tan \lambda +1}}$
    $\tan {\frac {\varphi -\psi }{2}}={\frac {\tan \lambda -1}{\tan \lambda +1}}\tan {\frac {\varphi +\psi }{2}}$
Risolviamo per $\varphi -\psi$ in modo da poter mettere a sistema somma e differenza.
    $\varphi -\psi =2\arctan \left({\frac {\tan \lambda -1}{\tan \lambda +1}}\tan {\frac {\varphi +\psi }{2}}\right)$
Volendo, si può semplificare anche ${\frac {\tan \lambda -1}{\tan \lambda +1}}$, riscrivendola come:
    ${\frac {\tan \lambda -1}{\tan \lambda +1}}={\frac {\tan \lambda -\tan {\frac {\pi }{4}}}{1+1\tan \lambda }}={\frac {\tan \lambda -\tan {\frac {\pi }{4}}}{1+\tan {\frac {\pi }{4}}\tan \lambda }}=\tan \left(\lambda -{\frac {\pi }{4}}\right)$
Ottenendo, come differenza di $\varphi$ e $\psi$:
    $\varphi -\psi =2\arctan \left[\tan \left(\lambda -{\frac {\pi }{4}}\right)\tan {\frac {\varphi +\psi }{2}}\right]$
Mettiamo quindi a sistema le due equazioni, impostando per comodità due valori $M$ e $N$.
    ${\begin{cases}M=\pi -(\alpha _{1}+\gamma )=\alpha _{2}+\beta _{1}\\N=2\arctan \left[\tan \left(\lambda -{\frac {\pi }{4}}\right)\tan {\frac {\varphi +\psi }{2}}\right]\\\end{cases}}$${\begin{cases}\varphi +\psi =M\\\varphi -\psi =N\\\end{cases}}$${\begin{cases}\varphi ={\frac {M+N}{2}}\\\psi ={\frac {M-N}{2}}\\\end{cases}}$
Il problema ora è ridotto ad una risoluzione di triangoli e si considera risolto.

#### Algoritmo risolutivo

Abbiamo quattro angoli ($\alpha _{1}$, $\beta _{1}$, $\alpha _{2}$, $\beta _{2}$) e la distanza $AB$. Il calcolo procede come segue:
* Calcolare $\gamma =\pi -\left(\alpha _{1}+\alpha _{2}+\beta _{2}\right),\quad \delta =\pi -\left(\alpha _{2}+\beta _{1}+\beta _{2}\right)$;
* Calcolare $\tan \lambda ={\frac {\sin \varphi }{\sin \psi }}={\frac {\sin \alpha _{2}\sin \beta _{2}\sin \gamma }{\sin \alpha _{1}\sin \beta _{1}\sin \delta }}$;
* Impostare $M=\pi -(\alpha _{1}+\gamma )=\alpha _{2}+\beta _{1},\quad N=2\arctan \left[\tan \left(\lambda -{\frac {\pi }{4}}\right)\tan {\frac {\varphi +\psi }{2}}\right]$ ;
* Calcolare gli angoli $\varphi ={\frac {M+N}{2}},\quad \psi ={\frac {M-N}{2}}$;
* Calcolare il lato PQ
    $PQ=AB{\frac {\sin \phi \sin \delta }{\sin \alpha _{2}\sin \beta _{2}}}$
o in maniera equivalente
    $PQ=AB{\frac {\sin \psi \sin \gamma }{\sin \alpha _{1}\sin \beta _{1}}}$
Se una delle due frazioni ha il denominatore tendente a zero, usare l'altra al fine di minimizzare gli errori angolari.

### Determinazione diretta degli angoli incogniti

In alternativa, si può sfruttare direttamente il rapporto tra il seno degli gli angoli incogniti ${\frac {\sin \varphi }{\sin \psi }}$, quantità calcolabile che chiamiamo $k$ o $\lambda$, e metterlo a sistema con la somma degli angoli incogniti, anch'essa calcolabile, che chiamiamo $\omega$. Si ottiene:
    ${\begin{cases}k={\frac {\sin(\varphi )}{\sin(\psi )}}\\\varphi +\psi =\omega \end{cases}}$
Concentriamoci sulla prima equazione: essa può essere riscritta come:
    $k={\frac {\sin(\omega -\psi )}{\sin(\psi )}}$
Espandiamo il numeratore utilizzando la formula di sottrazione del seno e semplifichiamo:
    $k={\frac {\sin(\omega )\cos(\psi )-\sin(\psi )\cos(\omega )}{\sin(\psi )}}={\frac {\sin(\omega )}{\tan(\psi )}}-\cos(\omega )$
Infine, isoliamo $\tan(\psi )$ e utilizziamo la funzione arcotangente:
    $\psi =\arctan \left({\frac {\sin \omega }{k+\cos \omega }}\right)$
Possiamo ottenere l'angolo $\varphi$ per differenza, ottenendo il seguente sistema:
    ${\begin{cases}\psi =\arctan \left({\frac {\sin \omega }{k+\cos \omega }}\right)\\\varphi =\omega -\psi \end{cases}}$
Il problema ora è ridotto ad una risoluzione di triangoli e si considera risolto.

### Risoluzione con il metodo della base fittizia

Il metodo della base (o del poligono, della figura) fittizia, si basa sulla proprietà delle figure simili della geometria euclidea. In particolare, si sfrutta la congruenza fra il rapporto dei lati corrispondenti.
La risoluzione si impone di creare un poligono $A'B'Q'P'$, simile al poligono $ABQP$, e di impostare la lunghezza della lato $P'Q'$, corrispondente alla base di stazione $PQ$, a un valore arbitrario. Per semplicità di calcolo è comune scegliere valori multipli di 10, di ordine di grandezza corrispondente a quello delle misure del rilievo. La conoscenza di $PQ$ banalizza il problema ad una semplice risoluzione di triangoli, che possono essere risolti con il Teorema dei seni e il Teorema di Carnot.
Una volta ricavati tutti i lati, le diagonali e gli angoli del poligono fittizio, sarà sufficiente applicare la proprietà delle figure simili, precedentemente richiamata, ai poligoni $ABQP$ e $A'B'Q'P'$:
    ${\frac {AB}{A'B'}}={\frac {l}{l'}}$
Notando che il rapporto ${\frac {AB}{A'B'}}$ (fattore di scala) è costante, può essere chiamato $f$ o $\lambda _{l}$ e possiamo risolvere tutti i lati del poligono reale con la seguente relazione:
    $l=l'\cdot f$
Gli angoli del poligono $A'B'Q'P'$, per definizione, saranno anche gli angoli del poligono $ABQP$.