Editing Entfernungsberechnung mit Quaternionen und Geodäten (section)

=== Verfahren zur Entfernungsberechnung ===

==== Beschreibung der Verfahrensschritte ====

Im ersten Verfahrensschritt wird der Rotationsellipsoid (Gleichung für $\mathbf{P}_i (\vartheta_i,\varphi_i)$) derart rotiert, dass der Punkt $\mathbf{P}_1$ mit $(\vartheta_1, \varphi_1)$ auf der $x$-$z$-Ebene liegt. Dazu ist eine Drehung um $\alpha_1 = - \varphi_1$ um die Applikate notwendig. Es entsteht ein gedrehter Rotationsellipsoid<ref>Die hochgestellte Ziffer in Klammern links von $\mathbf{P}$ bzw. $P$ symbolisiert die jeweilige Stufe der Transformation (1: erste Transformation, 2: zweite Transformation, etc. pp.).</ref>
\begin{gather}
	\Prefix^{(1)}{P} = [0, \Prefix^{(1)}{\mathbf{P}}] = q_R(\alpha_1, \mathbf{u}_1) \cdot [0, \mathbf{P}] \cdot \overline{q_R}(\alpha_1, \mathbf{u}_1) \\
	\notag \text{mit } \alpha_1 = -\varphi_1, \ \mathbf{u}_1 = \begin{pmatrix} 0 \\ 0 \\ 1 \end{pmatrix}, \\
	\notag \Prefix^{(1)}{\mathbf{P}} \xlongequal{\text{def.}} \begin{pmatrix} \Prefix^{(1)}{P}_x(\vartheta, \varphi, \ldots) \\ \Prefix^{(1)}{P}_y(\vartheta, \varphi, \ldots) \\ \Prefix^{(1)}{P}_z(\vartheta, \varphi, \ldots) \end{pmatrix}. 
\end{gather}

Analog zum ersten Verfahrensschritt wird der Ellipsoid nun um die Ordinate derart gedreht, dass $\mathbf{P}_1$ auf die Abszisse ausgerichtet wird:
\begin{gather}
	\Prefix^{(2)}{P} = q_R(\alpha_2, \mathbf{u}_2) \cdot \Prefix^{(1)}{P} \cdot \overline{q_R}(\alpha_2, \mathbf{u}_2) \\
	\notag \text{mit } \alpha_2 = \arctan \left( \frac{\Prefix^{(1)}{P}_z(\vartheta = \vartheta_1, \varphi = \varphi_1, \ldots)}{\Prefix^{(1)}{P}_x(\vartheta = \vartheta_1, \varphi = \varphi_1, \ldots)} \right), \\
	\notag \mathbf{u}_2 = \begin{pmatrix} 0 \\ 1 \\ 0 \end{pmatrix}, \ \Prefix^{(2)}{\mathbf{P}} \xlongequal{\text{def.}} \begin{pmatrix} \Prefix^{(2)}{P}_x(\vartheta, \varphi, \ldots) \\ \Prefix^{(2)}{P}_y(\vartheta, \varphi, \ldots) \\ \Prefix^{(2)}{P}_z(\vartheta, \varphi, \ldots) \end{pmatrix}. 
\end{gather}

Im letzten Transformationsschritt wird $\Prefix^{(2)}{\mathbf{P}}$ so gedreht, dass der Punkt $\mathbf{P}_2$ auf der $x$-$y$-Ebene liegt:
\begin{gather}
	\Prefix^{(3)}{P} = q_R(\alpha_3, \mathbf{u}_3) \cdot \Prefix^{(2)}{P} \cdot \overline{q_R}(\alpha_3, \mathbf{u}_3) \\
	\notag \text{mit } \alpha_3 = - \arctan \left( \frac{\Prefix^{(2)}{P}_z(\vartheta = \vartheta_2, \varphi = \varphi_2, \ldots)}{\Prefix^{(2)}{P}_y(\vartheta = \vartheta_2, \varphi = \varphi_2, \ldots)} \right), \\
	\notag \mathbf{u}_3 = \begin{pmatrix} 1 \\ 0 \\ 0 \end{pmatrix}, \ \Prefix^{(3)}{\mathbf{P}} \xlongequal{\text{def.}} \begin{pmatrix} \Prefix^{(3)}{P}_x(\vartheta, \varphi, \ldots) \\ \Prefix^{(3)}{P}_y(\vartheta, \varphi, \ldots) \\ \Prefix^{(3)}{P}_z(\vartheta, \varphi, \ldots) \end{pmatrix}. 
\end{gather}

Mit $\Prefix^{(3)}{\mathbf{P}}$ steht nun eine Parametergleichung für die transformierte Erde zur Verfügung. Beide Punkte $\Prefix^{(3)}{\mathbf{P}}_1$ bzw. $\Prefix^{(3)}{\mathbf{P}}_2$ im transformierten Koordinatensystem liegen nun in der $x$-$y$-Ebene; die $z$-Komponente ist jeweils Null.

==== Herleitung des Lösungswegs ====

Die Bogenlänge zwischen zwei Punkten auf einer ebenen, rektifizierbaren Kurve mit dem Parameter $t \mapsto \left( x(t), y(t) \right)^\mathrm{T}$ kann allgemein mit
\begin{equation}
	\alpha \leq t \leq \beta: \ \int \limits_{\alpha}^{\beta} \sqrt{\dot{x}(t)^2 + \dot{y}(t)^2} \, \dd t
\end{equation}
ermittelt werden. Hierbei hängt die vorhandene Parametergleichung $\Prefix^{(3)}{\mathbf{P}}$ jedoch von zwei Parametern $\vartheta$ und $\varphi$ ab. Eine Einschränkung auf einen Parameter ist daher notwendig, so dass eine zwei-dimensionale Parametergleichung
\begin{equation}
	\Prefix^{(4)}{\mathbf{P}} = \begin{pmatrix} \Prefix^{(4)}{P}_x (\varphi) \\ \Prefix^{(4)}{P}_y (\varphi) \end{pmatrix} = \begin{pmatrix} \Prefix^{(3)}{P}_x \left(\vartheta(\varphi), \varphi \right) \\ \Prefix^{(3)}{P}_y \left(\vartheta(\varphi), \varphi \right) \end{pmatrix}
\end{equation}
entsteht. Da für das Verfahren lediglich der Schnitt durch $z=0$ relevant ist, erhält man durch Lösung von $\Prefix^{(3)}{P}_z(\ldots) \stackrel{!}= 0$ einen Zusammenhang beider Parameter: $\vartheta(\varphi) = \ldots$

Schließlich ergibt sich der Abstand beider Punkte durch Integration von
\begin{equation}
		\int \limits_{\varphi_1}^{\varphi_2} \sqrt{\left( \frac{\dd}{\dd \varphi} \ \Prefix^{(4)}{P}_x(\varphi) \right)^2 + \left( \frac{\dd}{\dd \varphi} \ \Prefix^{(4)}{P}_y(\varphi) \right)^2} \, \dd \varphi.
\end{equation}