Java·Applied·Geodesy·3D

Least-Squares Adjustment Software for Geodetic Sciences

Benutzer-Werkzeuge

Webseiten-Werkzeuge


least-squares-adjustment:deformationanalysis

Dies ist eine alte Version des Dokuments!


Deformationsanalyse

Das Ausgleichungsprogramm JAG3D erlaubt neben der reinen Netzausgleichung auch eine Kongruenzanalyse, bei der ein geometrischer Vergleich zweier Epochen durchgeführt wird. Auf der Basis von statistischen Tests soll im Zuge der Deformationsanalyse geprüft werden, ob ein Punkt oder eine Punktgruppe sich gegenüber eine Referenz verschoben hat. In der gegenwärtigen Literatur werden zwei Verfahren propagiert, die eine Analyse eines Überwachungsnetzes beschreiben (vgl. Jäger et al. 2005):

  • Koordinatenmethode und
  • Beobachtungsmethode.

Kongruenzmodell

Bei der Koordinatenmethode werden zunächst die Einzelepochen getrennt ausgewertet. Die anschließende Deformationsanalyse erfolgt auf der Basis der ausgeglichenen Koordinaten der Einzelepochenauswertung. Das zugrunde zulegende stochastische Modell basiert dabei auf den Varianz-Kovarianz-Matrizen der Einzelepochen. Interepochale Abhängigkeiten werden in diesem Modell vernachlässigt.

Auch bei der Beobachtungsmethode sollte zunächst die Einzelauswertung der Epochen vorausgehen. Im Rahmen dieser Auswertung ist das Beobachtungsmaterial der jeweiligen Epoche auf Fehlmessungen hin zu analysieren und das stochastische Modell festzulegen. Im Gegensatz zur Koordinatenmethode erfolgt in der Kongruenzanalyse nun eine gemeinsame Ausgleichung beider Epochen. Hierzu werden die Beobachtungen $\mathbf{l_1}$ und $\mathbf{l_2}$ beider Epochen in einem Ausgleichungsmodell zusammengeführt. Das stochastische Modell der gemeinsamen Ausgleichung leitet sich aus den a-priori Unsicherheiten der originären Beobachtungen der Einzelepochenauswertungen ab. Beim Zusammenführen beider Epochen auf der Basis der originären Beobachtungen entstehen interepochale Korrelationen, die bei der Auswertung berücksichtigt werden. JAG3D verarbeitet bei der Kongruenzanalyse die Beobachtungen der Epochen in einem Guß. Aus diesem Grund wird im Folgenden nur das Modell der Beobachtungsmethode näher beschrieben werden. Detaillierte Informationen zu den Modellen kann u.a. Pelzer (1971) oder Jäger und Drixler (1990) entnommen werden.

Das Ausgleichungsmodell der gemeinsamen Auswertung zweier Epochen im Gauß-Markov-Modell ergibt sich zu:

$$\begin{pmatrix}\mathbf{l_1} \\ \mathbf{l_2} \end{pmatrix} + \begin{pmatrix}\mathbf{v_1} \\ \mathbf{v_2} \end{pmatrix} = \begin{pmatrix}\mathbf{A}_{R,1} & \mathbf{A}_{O,1} & \mathbf{0} \\ \mathbf{A}_{R,2} & \mathbf{0} & \mathbf{A}_{O,2} \end{pmatrix} \begin{pmatrix}\mathbf{x}_R \\ \mathbf{x}_{O,1} \\ \mathbf{x}_{O,2} \end{pmatrix}$$

Bei der gemeinsamen freien Ausgleichung der Epochen werden die als stabil angenommenen Referenzkoordinaten $\mathbf{x}_R$ über alle Epochen als identische unbekannte Parameter geschätzt. Hierdurch wird ein einheitliches Datum zwischen den Einzelepochen erreicht. Die übrigen Objektpunkte werden hingegen mit $\mathbf{x}_{O,1}$ und $\mathbf{x}_{O,2}$ für jede Epoche getrennt ins Modell eingeführt.

Explizite Prüfung der Referenzpunkte

Bevor eine Analyse der Objektpunkt $\mathbf{x}_{O}$ erfolgen kann, sind die als stabil angenommenen Referenzpunkte $\mathbf{x}_R$ zunächst auf Invarianz zu prüfen. Für die explizite Suche nach veränderten Referenzpunkten ist das o.g. Gauß-Markov-Modell um den (möglichen) Deformationsanteil $\mathbf{\nabla}_{R,j}$ des zu prüfenden Referenzpunktes zu erweitert. Die Vorgehensweise ist hierbei analog zur Modellerweiterung bei der Prüfung auf Modellstörungen.

$$\mathbf{l+v} = \mathbf{Ax} + \mathbf{B}_{R,j} \mathbf{\nabla}_{R,j}$$

Die Größe der geschätzten $j$-ten Modellstörung $\mathbf{\nabla}_{R,j}$ lässt sich wiederum auf Signifikanz mit den beiden bekannten Teststatistiken prüfen:

$$ T_{prio,j} = \frac{\mathbf{\nabla}_{R,j}^T\mathbf{Q}_{\mathbf{\nabla\nabla}_{R,j}}^{-1}\mathbf{\nabla}_{R,j}} {m\sigma_0^2} \sim F_{m,\infty} | \mathrm{H}_0$$

und

$$ T_{post,j} = \frac{\mathbf{\nabla}_{R,j}^T\mathbf{Q}_{\mathbf{\nabla\nabla}_{R,j}}^{-1}\mathbf{\nabla}_{R,j}} {m{{\hat{\sigma}_j'}^2}} \sim F_{m,f-m} | \mathrm{H}_0$$

Hierin bezeichnet ${{\hat{\sigma}_j'}^2}$ den um den Einfluss der Modellstörung $\mathbf{\nabla}$ reduzierten a-posteriori Varianzfaktor und $m = \tr\left(\mathbf{Q}_{\mathbf{\nabla\nabla}}\right)$.

Analog zum Ausreißertest ist die Nullhypothese $\mathrm{H}_0$ - es liegt keine Veränderung im $j$-ten Referenzpunkt vor - zu verwerfen, wenn die Teststatistik $T_j$ größer als das zugehörige Quantil der F-Verteilung ist.

Ableitung der Verschiebungsvektoren der Objektpunkte

Liegen keine instabilen Referenzpunkte mehr im Modell vor, wird die Größe der Punktverschiebung $\mathbf{\nabla}_k$ des $k$-ten Objektpunktes mittels der Formmatrix $\mathbf{F}_k = \begin{pmatrix}\mathbf{0} & \cdots & \mathbf{0} & \mathbf{-I}_k & \mathbf{0} & \cdots & \mathbf{0} & \mathbf{I}_k & \mathbf{0} & \cdots \mathbf{0} \end{pmatrix}$

$$\mathbf{\nabla}_k=\mathbf{F}_k\begin{pmatrix} \mathbf{x}_{O,1} \\ \mathbf{x}_{O,2} \end{pmatrix}$$

zusammen mit der zugehörigen Varianz-Kovarianz-Matrix bestimmt.

$$\mathbf{Q}_{\mathbf{\nabla\nabla},k} = \mathbf{F}_k \begin{pmatrix}\mathbf{Q}_{\mathbf{x}_{O,1}\mathbf{x}_{O,1}} & \mathbf{Q}_{\mathbf{x}_{O,1}\mathbf{x}_{O,2}} \\ \mathbf{Q}_{\mathbf{x}_{O,2}\mathbf{x}_{O,1}} & \mathbf{Q}_{\mathbf{x}_{O,2}\mathbf{x}_{O,2}} \end{pmatrix} \mathbf{F}_k^T$$

Die Größe dieses Deformationsparametervektors $\mathbf{\nabla}_k$ lassen sich wiederum auf Signifikanz prüfen. Analog zum Referenzpunkttest lauten die beiden Teststatistiken:

$$ T_{prio,k} = \frac{\mathbf{\nabla}_k^T\mathbf{Q}_{\mathbf{\nabla\nabla},k}^{-1}\mathbf{\nabla}_k} {m\sigma_0^2} \sim F_{m,\infty} | \mathrm{H}_0$$

und

$$ T_{post,k} = \frac{\mathbf{\nabla}_k^T\mathbf{Q}_{\mathbf{\nabla\nabla},k}^{-1}\mathbf{\nabla}_k} {m\hat{\sigma}_0^2} \sim F_{m,f} | \mathrm{H}_0$$

Die detektierte Veränderung gilt wiederum als signifikant, wenn die Nullhypothese $\mathrm{H}_0$ zu verwerfen und die Teststatistik $T_k$ größer als das zugehörige Quantil der F-Verteilung ist.

least-squares-adjustment/deformationanalysis.1650792508.txt.gz · Zuletzt geändert: 2022/04/24 11:28 von Michael Lösler