====== 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 [[:least-squares-adjustment#stochastisches_modell|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 [[:least-squares-adjustment#funktionales_modell|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 [[:least-squares-adjustment:configuration#freie_netzausgleichung|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 [[:least-squares-adjustment:outlier|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 [[:least-squares-adjustment:outlier|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}_{R,j}$ reduzierten //a-posteriori// Varianzfaktor und $m = \tr\left(\mathbf{Q}_{\mathbf{\nabla\nabla}}\right)$. Analog zum [[:least-squares-adjustment:outlier|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.