Java·Applied·Geodesy·3D

Least-Squares Adjustment Software for Geodetic Sciences

Benutzer-Werkzeuge

Webseiten-Werkzeuge


least-squares-adjustment:observation

Dies ist eine alte Version des Dokuments!


Beobachtungen

JAG3D unterstützt in der Netzausgleichung verschiedene Beobachtungsarten, die innerhalb eines Projektes miteinander kombiniert werden können. Voraussetzung ist, dass die jeweilige Beobachtung zwischen Punkten definiert wird, die diesen Typ unterstützen. Die folgende Tabelle gibt einen Überblick über alle verfügbaren Beobachtungstypen und zwischen welchen Punkten diese in der Ausgleichung berücksichtigt werden können.

1D-Punkt 2D-Punkt 3D-Punkt
Nivellement x x
Richtung/Azimut x x
Horizontalstrecke x x
Raumstrecke x
Zenitwinkel x
GNSS 1D-Basislinie x x
GNSS 2D-Basislinie x x
GNSS 3D-Basislinie x
1D-Punkt x
2D-Punkt x
3D-Punkt x

Das funktionale und stochastische Modell der einzelnen Beobachtungen soll im Folgenden erläutert werden. Beobachtungen sind immer zwischen zwei Punkten definiert, dem Standpunkt Ps und dem Zielpunkt Pz. Deren Koordinaten seien Ps = [ys xs zs]T und Pz = [yz xz zz]T, wenn es sich um Raumpunkte handelt. Liegen die Punkte in einer niederen Dimension vor, so sind die entsprechenden Koordinatenkomponenten zu streichen. Mit hs und hz sind die Stand- und Zielpunkthöhe gegeben. Ferner setzt sich das stochastische Modell aus konstanten und entfernungsabhängigen Anteilen zusammen. Die Distanz zur Bestimmung der entfernungsabhängigen Anteile sei mit d bezeichnet und ρ = 200/π zur Umrechnung von Radiant in Neugrad.

JAG3D unterstützt eine integrierte 3D-Netzausgleichung für terrestrische Beobachtungen, Nivellements und GNSS-Beobachtungen. Hierbei berücksichtigt die Software eine mögliche nicht-Parallelität zwischen den Stehachsen der (Stand-)Punkte infolge von vorhandenen Lotabweichungen. Werden die Koordinatendifferenzen zwischen zwei Punkten mit

$$ \begin{pmatrix} \Delta y \\ \Delta x \\ \Delta z \\ \end{pmatrix} = \begin{pmatrix} y_z - y_s \\ x_z - x_s \\ z_z - z_s \\ \end{pmatrix} $$

bezeichnet, so ergibt sich die allg. Beobachtungsgleichung

$$\begin{pmatrix}v\\u\\w\end{pmatrix}=\mathbf{R_s}\begin{pmatrix}\Delta y\\\Delta x\\\Delta z\end{pmatrix}-\begin{pmatrix}0\\0\\h_s\end{pmatrix}+\mathbf{R_sR^T_z}\begin{pmatrix}0\\0\\h_z\end{pmatrix}$$

worin zur Modellierung der Lotabweichungen im Standpunkt die Rotationssequenz Rs und im Zielpunkt die Rotationssequenz Rz eingeführt werden. Die Rotationssequenz für jeden Punkt ergibt sich aus einer kombinierten Drehung um die x und y-Achse

$$\mathbf{R}=\begin{pmatrix}1&0&0\\0&\cos\zeta_x&-\sin\zeta_x\\0&\sin\zeta_x&\cos\zeta_x\\\end{pmatrix}\begin{pmatrix}\cos\zeta_y&0&\sin\zeta_y\\0&1&0\\-\sin\zeta_y&0&\cos\zeta_y\\\end{pmatrix}$$

worin ζx und ζy die Drehwinkel zur Beschreibung der Lotabweichungen darstellen.

Das im folgenden aufgeführte stochastische Modell wird herangezogen, wenn keine individuellen Genauigkeiten für die einzelnen Beobachtungen vorliegen. In diesem Fall greift der gruppenbasierte Ansatz.

Terrestrische Beobachtungen

Im folgenden werden die terrestrischen Beobachtungsgleichungen, wie sie von JAG3D unterstützt werden, aufgeführt. Neben dem funktionalen Modell wird das stochastische Modell vorgestellt und mögliche Gruppenparameter, die wahlweise als zusätzliche Unbekannte geschätzt werden können, genannt.

Hinweis: Das stochastische Modell setzt sich stets aus einem konstanten und zwei entfernungsabhängigen Anteilen zusammen. Für die entfernungsabhängigen Anteilen muss die Zielweite d gegeben sein. Sollte d nicht explizit vorliegen, berechnet JAG3D die Zielweite näherungsweise aus den Näherungskoordinaten. In Abhängigkeit der Güte der Näherungskoordinaten kann das stochastische Modell daher variieren, sodass die explizite Vorgabe der Zielweite empfohlen wird.

Nivellement

Funktionales Modell $$\delta h = \frac{1} m w$$
Stochastisches Modell $$\sigma_{\delta h} = \sqrt{\sigma_a^2 + \sigma_b^2 d + (\sigma_c d)^2}$$
Punktdimension 1D, 3D
Zusatzparameter Maßstab m
Einheit Meter [m]

(Bemerkung: d Nivellementsweg in [km])

Richtung/Azimut

Funktionales Modell $$t = \arctan2{ \left ( v, u \right )} \rho - o$$
Stochastisches Modell $$\sigma_t = \sqrt{\sigma_a^2 + \left(\frac{\sigma_b} {\sqrt d} \rho \right)^2 + \left(\frac{\sigma_c} d \rho \right)^2 } $$
Punktdimension 2D, 3D
Zusatzparameter Orientierung o
Einheit Neugrad [gon]

(Bemerkung: Bei Azimuten entfällt unter Umständen o)

Horizontale Strecke

Funktionales Modell $$s_{2D} = \frac{1} m \left(\sqrt{ v^2 + u^2}- a \right)$$
Stochastisches Modell $$\sigma_{s_{2D}} = \sqrt{\sigma_a^2 + \sigma_b^2 d + \sigma_c^2 d^2} $$
Punktdimension 2D, 3D
Zusatzparameter Maßstab m, Additionskonstante a
Einheit Meter [m]

(Bemerkung: d Distanz in [m])

Schrägstrecke

Funktionales Modell $$s_{3D} = \frac{1} m \left(\sqrt{ v^2 + u^2 + w^2}- a \right)$$
Stochastisches Modell $$\sigma_{s_{3D}} = \sqrt{\sigma_a^2 + \sigma_b^2 d + \sigma_c^2 d^2} $$
Punktdimension 3D
Zusatzparameter Maßstab m, Additionskonstante a
Einheit Meter [m]

(Bemerkung: d Distanz in [m])

Zenitwinkel

Funktionales Modell $$z = \arctan{ \frac{\sqrt{ v^2 + u^2}} {w}} - k \frac{s_{2D}} {2R} \rho$$
Stochastisches Modell $$\sigma_z = \sqrt{\sigma_a^2 + \left(\frac{\sigma_b} {\sqrt d} \rho \right)^2 + \left(\frac{\sigma_c} d \rho \right)^2} $$
Punktdimension 3D
Zusatzparameter Refraktionskoeffizient k
Einheit Neugrad [gon]

(Bemerkung: R = 6371km entspricht dem mittleren Erdradius)

GNSS-Basislinien

Im Gegensatz zu den terrestrischen Beobachtungen, bei denen zwischen den Punkten Ps und Pz ein einziger Messwert vorliegt, wird die GNSS-Basislinie in Abhängigkeit der Dimension in eine b1D = [δz]T, zwei b2D = [δy δx]T oder drei b3D = [δy δx δz]T Vektorkomponenten aufgesplittet.

1D-Basislinien

Funktionales Modell $$\mathbf{b_{1D}} = \delta z = m\mathbf{R} \begin{pmatrix} y_z-y_s \\ x_z-x_s \\ z_z-z_s \end{pmatrix} $$
Stochastisches Modell $$\sigma_{\delta z} = \sqrt{\sigma_a^2 + \sigma_b^2 \delta z + \sigma_c^2 \delta z^2}$$
Punktdimension 1D, 3D
Zusatzparameter Maßstab m, 2 Drehwinkel in Matrix R
Einheit Meter [m]

(Bemerkung: Bei 1D-Punkten müssen die x und y-Koordinaten nur genähert bekannt sein.)

2D-Basislinien

Funktionales Modell $$\mathbf{b_{2D}} = \begin{pmatrix} \delta y \\ \delta x \end{pmatrix} = m\mathbf{R} \begin{pmatrix} y_z-y_s \\ x_z-x_s \end{pmatrix} $$
Stochastisches Modell $$\mathbf{C_b} = \begin{pmatrix} \sigma_{\delta y}^2 & 0 \\ 0 & \sigma_{\delta x}^2 \end{pmatrix}$$ mit $\sigma_{\delta} = \sqrt{\sigma_a^2 + \sigma_b^2 \delta + \sigma_c^2 \delta^2}$
Punktdimension 2D, 3D
Zusatzparameter Maßstab m, Drehwinkel in Matrix R
Einheit Meter [m]

3D-Basislinien

Funktionales Modell $$\mathbf{b_{3D}} = \begin{pmatrix} \delta y \\ \delta x \\ \delta z \end{pmatrix} = m\mathbf{R} \begin{pmatrix} y_z-y_s \\ x_z-x_s \\ z_z-z_s \end{pmatrix} $$
Stochastisches Modell $$\mathbf{C_b} = \begin{pmatrix} \sigma_{\delta y}^2 & 0 & 0 \\ 0 & \sigma_{\delta x}^2 & 0 \\ 0 & 0 &\sigma_{\delta z}^2 \end{pmatrix}$$ mit $\sigma_{\delta} = \sqrt{\sigma_a^2 + \sigma_b^2 \delta + \sigma_c^2 \delta^2}$
Punktdimension 3D
Zusatzparameter Maßstab m, 3 Drehwinkel in Matrix R
Einheit Meter [m]

Punktbeobachtungen

Im Rahmen einer weichen Lagerung bzw. dynamischen Ausgleichung werden den Anschlußpunkten Unsicherheitsbeiträge zugestanden. Dies erscheint insofern gerechtfertigt, da die Punkte, mit denen der Netzanschluß zu realisieren ist, aus vorherigen Messungen stammen und somit für jeden Punkt neben der Koordinate selbst eine Unsicherheit vorliegt. Die Matrix I bezeichnet im Folgenden eine Einheitsmatrix.

1D-Punkt

Funktionales Modell $$\mathbf{P_{1D}} = z = \mathbf{I} z$$
Stochastisches Modell $$\mathbf{C_P} = \sigma_z^2$$
Punktdimension 1D
Einheit Meter [m]

2D-Punkt

Funktionales Modell $$\mathbf{P_{2D}} = \begin{pmatrix} y \\ x \end{pmatrix} = \mathbf{I} \begin{pmatrix} y \\ x \end{pmatrix} $$
Stochastisches Modell $$\mathbf{C_P} = \begin{pmatrix} \sigma_y^2 & 0 \\ 0 & \sigma_x^2 \end{pmatrix}$$
Punktdimension 2D
Einheit Meter [m]

3D-Punkt

Funktionales Modell $$\mathbf{P_{3D}} = \begin{pmatrix} y \\ x \\ z \end{pmatrix} = \mathbf{I} \begin{pmatrix} y \\ x \\ z \end{pmatrix} $$
Stochastisches Modell $$\mathbf{C_P} = \begin{pmatrix} \sigma_y^2 & 0 & 0 \\ 0 & \sigma_x^2 & 0 \\ 0 & 0 & \sigma_z^2 \end{pmatrix}$$
Punktdimension 3D
Einheit Meter [m]
least-squares-adjustment/observation.1520791097.txt.gz · Zuletzt geändert: 2018/03/11 18:58 von Michael Lösler