Java·Applied·Geodesy·3D

Least-Squares Adjustment Software for Geodetic Sciences

Benutzer-Werkzeuge

Webseiten-Werkzeuge


least-squares-adjustment:reliability

Unterschiede

Hier werden die Unterschiede zwischen zwei Versionen angezeigt.

Link zu dieser Vergleichsansicht

Beide Seiten der vorigen RevisionVorhergehende Überarbeitung
Nächste ÜberarbeitungBeide Seiten der Revision
least-squares-adjustment:reliability [2018/03/11 20:34] – [Genauigkeitsmaße der Parameter] Michael Löslerleast-squares-adjustment:reliability [2018/03/11 22:07] Michael Lösler
Zeile 10: Zeile 10:
 </figure> </figure>
  
-Die Genauigkeitsmaße der geschätzten Parameter z.B. die Koordinaten der Netzpunkte werden nach der [[:least-squares-adjustment|Ausgleichungsrechnung]] aus der Varianz-Kovarianz-Matrix ''**Cxx**'' der Parameter abgeleitet. +Die Genauigkeitsmaße der geschätzten Parameter z.B. die Koordinaten der Netzpunkte werden nach der [[:least-squares-adjustment|Ausgleichungsrechnung]] aus der Varianz-Kovarianz-Matrix $\mathbf{C_{\hat x \hat x}}$ der Parameter abgeleitet. 
  
 $$\mathbf{C_{\hat x \hat x}}=\hat{\sigma}_0^2\mathbf{Q_{\hat x \hat x}}$$ $$\mathbf{C_{\hat x \hat x}}=\hat{\sigma}_0^2\mathbf{Q_{\hat x \hat x}}$$
  
-Während auf der Hauptdiagonalen die Varianzen ''σ<sup>2</sup>'' der Parameter angeordnet sind, befinden sich auf den Nebendiagonalen die Kovarianzen ''ρ'', welche die (linearen) Abhängigkeiten zwischen den Parametern quantifiziert. Während für ''ρ > 0'' eine positive Korrelation besteht, sind die Parameter für ''ρ < 0'' negativ korreliert. Für ''ρ = 0'' sind die Parameter vollständig (linear) unabhängig.+Während auf der Hauptdiagonalen die Varianzen $\sigma^2der Parameter angeordnet sind, befinden sich auf den Nebendiagonalen die Kovarianzen $\rho$, welche die (linearen) Abhängigkeiten zwischen den Parametern quantifiziert. Während für $\rho \gt 0eine positive Korrelation besteht, sind die Parameter für $\rho \lt 0negativ korreliert. Für $\rho = 0sind die Parameter vollständig (linear) unabhängig.
  
-Die geschätzten Standardunsicherheiten ''σ'' der Parameter z.B. der ausgeglichenen Koordinaten eines Punktes können demnach direkt aus der Varianz-Kovarianz-Matrix ''**Cxx**'' abgelesen werden. Zu beachten ist, dass die Varianz-Kovarianz-Matrix ''**Cxx**'' datumsabhängig ist. Dies bedeutet, dass bspw. die Standardunsicherheiten ''σ'' nicht invariant gegenüber Netzverdrehungen sind und sich daher nur bedingt zur Beurteilung eignen. Ein punktbezogenes rotationsinvariantes Genauigkeitsmaß sind hingegen die Halbachsen der Konfidenzbereiche. +Die geschätzten Standardunsicherheiten $\sigma$ der Parameter z.B. der ausgeglichenen Koordinaten eines Punktes können demnach direkt aus der Varianz-Kovarianz-Matrix $\mathbf{C_{\hat x \hat x}}$ abgelesen werden. Zu beachten ist, dass die Varianz-Kovarianz-Matrix $\mathbf{C_{\hat x \hat x}}$ datumsabhängig ist. Dies bedeutet, dass bspw. die Standardunsicherheiten $\sigma$ nicht invariant gegenüber Netzverdrehungen sind und sich daher nur bedingt zur Genauigkeitsbeurteilung eignen. Ein punktbezogenes rotationsinvariantes Genauigkeitsmaß sind hingegen die Halbachsen der Konfidenzbereiche. 
  
-Durch spektrale Zerlegung der punktbezogenen Sub-Kofaktormatrix ''**Qx**<sub>i</sub>**x**<sub>i</sub>'' werden die Halbachsen und die Orientierung des ''n''-dimensionalen Hyperellipsoids gewonnen. Während sich für ''n = 3'' ein Ellipsoid ergibt, erhält man für ''n = 2'' eine Ellipse und für ''n = 1'' ein einfaches Intervall. Die Halbachsen bilden somit die Extremstellen des Konfidenzbereichs.+Durch spektrale Zerlegung der punktbezogenen Sub-Kofaktormatrix $\mathbf{Q_{\hat \hat x}}_{ii}$ werden die Halbachsen und die Orientierung des $n$-dimensionalen Hyperellipsoids gewonnen. Während sich für $n=3ein Ellipsoid ergibt, erhält man für $n = 2eine Ellipse und für $n = 1ein einfaches Intervall. Die Halbachsen bilden somit die Extremstellen des Konfidenzbereichs.
  
 $$\mathbf{Q_{\hat x \hat x}}= $$\mathbf{Q_{\hat x \hat x}}=
 \mathbf{\begin{bmatrix} \mathbf{\begin{bmatrix}
-Q_{\hat x_i \hat x_i}&Q_{\hat x_i \hat x_j}& \dots\\ +\mathbf{Q_{\hat \hat x}}_{ii} & \mathbf{Q_{\hat \hat x}}_{ij} & \dots\\ 
-Q_{\hat x_j \hat x_i}&Q_{\hat x_j \hat x_j}& \dots\\+\mathbf{Q_{\hat \hat x}}_{ji} & \mathbf{Q_{\hat \hat x}}_{jj} & \dots\\
 \vdots & \vdots & \ddots \vdots & \vdots & \ddots
 \end{bmatrix}}$$ \end{bmatrix}}$$
  
-{{ref>confidence_ellipse}} stellt für zwei Punkte die zugehörigen Konfidenzellipsen gegenüber. Während die Standardunsicherheiten ''σ<sub>x</sub> σ<sub>y</sub> = 5 mm'' bei beiden Punkten identisch sind, lässt sich aus der Form der Ellipsen erkennen, dass der Konfidenzbereich der roten Ellipse kleiner ist als der der blauen. Beide Punkte weisen somit nicht dieselbe Unsicherheit auf. Die Koordinatenachsen-bezogenen Standardunsicherheiten sind demnach nur aussagekräftig, wenn diese mit den Halbachsen des Konfidenzbereichs (näherungsweise) zusammenfallen, siehe blaue Ellipse. Die Forderung nach //homogenen// und //isotropen// Konfidenzbereichen bedeutet, dass im 2D-Fall die Ellipse zu einem Kreis und im 3D-Fall das Ellipsoid zu einer Kugel entartet.+{{ref>confidence_ellipse}} stellt für zwei Punkte die zugehörigen Konfidenzellipsen gegenüber. Während die Standardunsicherheiten $\sigma_x \sigma_y = 5~\rm{mm}$ bei beiden Punkten identisch sind, lässt sich aus der Form der Ellipsen erkennen, dass der Konfidenzbereich der roten Ellipse kleiner ist als der der blauen. Beide Punkte weisen somit nicht dieselbe Unsicherheit auf. Die Koordinatenachsen-bezogenen Standardunsicherheiten sind demnach nur aussagekräftig, wenn diese mit den Halbachsen des Konfidenzbereichs (näherungsweise) zusammenfallen, siehe blaue Ellipse. Die Forderung nach //homogenen// und //isotropen// Konfidenzbereichen bedeutet, dass im 2D-Fall die Ellipse zu einem Kreis und im 3D-Fall das Ellipsoid zu einer Kugel entartet.
  
-Im Folgenden soll explizit nur der räumliche Fall skizziert werden. Für ''3'' vereinfacht sich die Darstellung sinngemäß. Die ''//i//''-te punktbezogene Sub-Kofaktormatrix ''**Qx**<sub>i</sub>**x**<sub>i</sub>'' lautet+Im Folgenden soll explizit nur der räumliche Fall skizziert werden. Für $\lt 3vereinfacht sich die Darstellung sinngemäß. Die $i$-te punktbezogene Sub-Kofaktormatrix $\mathbf{Q_{\hat \hat x}}_{ii}$ lautet
  
-$$\hat{\sigma}_0^2\mathbf{Q_{\hat x_i \hat x_i}}=+$$\hat{\sigma}_0^2\mathbf{Q_{\hat \hat x}}_{ii}=
 \begin{pmatrix} \begin{pmatrix}
 \sigma_x^2&\rho_{xy}&\rho_{xz}\\ \sigma_x^2&\rho_{xy}&\rho_{xz}\\
Zeile 38: Zeile 38:
 \end{pmatrix}_{ii}$$ \end{pmatrix}_{ii}$$
  
-Die drei Halbachsen ''A''''B'' und ''C'' ergeben sich aus den (sortierten) Eigenwerten ''**λ** = [λ<sub>max</sub>λ<sub>mid</sub>λ<sub>min</sub>]<sup>T</sup>''+Die drei Halbachsen $A$$Bund $Cergeben sich aus den (sortierten) Eigenwerten $\mathbf{\lambda} = [\lambda_{max}\lambda_{mid}\lambda_{min}]^T$
  
 $A = \sigma_0 \sqrt{ \lambda_{max} }$, $B = \sigma_0 \sqrt{ \lambda_{mid} }$, $C = \sigma_0 \sqrt{ \lambda_{min} }$ $A = \sigma_0 \sqrt{ \lambda_{max} }$, $B = \sigma_0 \sqrt{ \lambda_{mid} }$, $C = \sigma_0 \sqrt{ \lambda_{min} }$
  
-Die Orientierung lässt sich aus den drei korrespondierenden Eigenvektoren ''**m**'' ableiten. +Die Orientierung lässt sich aus den drei korrespondierenden Eigenvektoren $\mathbf{m}$ ableiten. 
  
 $$\mathbf{M}=\left [\begin{matrix} \mathbf{m}_{max}&\mathbf{m}_{mid}&\mathbf{m}_{min} \end{matrix} \right ]$$ $$\mathbf{M}=\left [\begin{matrix} \mathbf{m}_{max}&\mathbf{m}_{mid}&\mathbf{m}_{min} \end{matrix} \right ]$$
  
-worin die Modalmatrix ''**M**'' eine orthogonale Matrix mit ''det(**M**) = +1'' ist, die sich in drei Einzelrotationen mit den Eulerwinkeln ''α''''β'' und ''γ'' zerlegen lässt.+worin die Modalmatrix $\mathbf{M}$ eine orthogonale Matrix mit $\det \mathbf{M= +1ist, die sich in drei Einzelrotationen mit den Eulerwinkeln $\alpha$$\beta$ und $\gamma$ zerlegen lässt.
  
 $$\mathbf{M}=\mathbf{R_x}(\alpha)\mathbf{R_y}(\beta)\mathbf{R_z}(\gamma)$$ $$\mathbf{M}=\mathbf{R_x}(\alpha)\mathbf{R_y}(\beta)\mathbf{R_z}(\gamma)$$
  
-Die folgende Tabelle zeigt die maximale Sicherheitswahrscheinlichkeit in Abhängigkeit der Dimension ''n'' für die so gebildeten einfachen Konfidenzbereiche.+Die folgende Tabelle zeigt die maximale Sicherheitswahrscheinlichkeit in Abhängigkeit der Dimension $nfür die so gebildeten einfachen Konfidenzbereiche.
  
 ^  n  ^  Wahrscheinlichkeit  ^ ^  n  ^  Wahrscheinlichkeit  ^
Zeile 57: Zeile 57:
 |  3  |  19,87 %  | |  3  |  19,87 %  |
  
-Aufgrund der geringen Sicherheitswahrscheinlichkeit des einfachen Konfidenzbereichs von ''40 %'' bzw. ''20 %'' für einen Lage- bzw. Raumpunkt, werden die Halbachsen i.d.R. mit einem Quantil der ''χ<sup>2</sup><sub>n</sub>''-Verteilung skaliert, sodass sich ein vergrößerter Konfidenzbereich mit einer entsprechend höheren Sicherheitswahrscheinlichkeit ergibt. In JAG3D wird die Sicherheitswahrscheinlichkeit aus der gewählten [[:least-squares-adjustment:outlier#teststatistik|Irrtumswahrscheinlichkeit]] abgeleitet.+Aufgrund der geringen Sicherheitswahrscheinlichkeit des einfachen Konfidenzbereichs von 40 % bzw. 20 % für einen Lage- bzw. Raumpunkt, werden die Halbachsen i.d.R. unter der Annahme der Normalverteilung mit einem Quantil der $\chi^2_n$-Verteilung skaliert, sodass sich ein vergrößerter Konfidenzbereich mit einer entsprechend höheren Sicherheitswahrscheinlichkeit ergibt. In JAG3D wird die Sicherheitswahrscheinlichkeit aus der gewählten [[:least-squares-adjustment:outlier#teststatistik|Irrtumswahrscheinlichkeit]] abgeleitet.
  
  
 ===== Redundanzanteil ===== ===== Redundanzanteil =====
  
-Ein Maß für die Kontrolliertheit liefert die Redundanzmatrix ''**R**''+Ein Maß für die Kontrolliertheit liefert die Redundanzmatrix $\mathbf{R}$
  
 $$\mathbf{R} = \mathbf{Q_{vv}P}$$ $$\mathbf{R} = \mathbf{Q_{vv}P}$$
  
-worin ''**P**'' die Gewichtsmatrix des [[:least-squares-adjustment#stochastisches_modell|stochastischen Modells]] und ''**Q<sub>vv</sub>**'' die Kofaktormatrix der Beobachtungsresiduen sind. Auf der Hauptdiagonale von ''**R**'' stehen die sogenannten Redundanzanteile ''//r<sub>i</sub>//'' für jede Beobachtung. Der Redundanzanteil ''//r//'' einer Beobachtung ist ein normiertes Maß und zeigt an, wie gut die betreffende Beobachtung durch die übrigen Beobachtungen kontrolliert wird im Netz. Während ein Redundanzanteil von ''//r=0//'' bedeutet, dass die Beobachtung nicht kontrolliert ist, repräsentiert ''//r=1//'' eine vollständige Kontrolle. Gelegentlich wird ''//r<sub>i</sub>//'' in Prozent angeben und durch den Parameter ''//EV<sub>i</sub>[%] = 100·r<sub>i</sub>//'' ausgedrückt. Als Faustregel für den Redundanzanteile ''//r<sub>i</sub>//'' gilt: +worin $\mathbf{P}$ die Gewichtsmatrix des [[:least-squares-adjustment#stochastisches_modell|stochastischen Modells]] und $\mathbf{Q_{vv}}$ die Kofaktormatrix der Beobachtungsresiduen sind. Auf der Hauptdiagonalen von $\mathbf{R}$ stehen die sogenannten Redundanzanteile $r_i$ für jede Beobachtung. Der Redundanzanteil $r_i$ einer Beobachtung ist ein normiertes Maß und zeigt an, wie gut die betreffende Beobachtung durch die übrigen Beobachtungen kontrolliert wird. Während ein Redundanzanteil von $r_i = 0bedeutet, dass die Beobachtung nicht kontrolliert ist, repräsentiert $r_i = 1eine vollständige Kontrolle. Gelegentlich wird $r_i$ in Prozent angeben und durch den Parameter $EV_i [\%] = 100 \cdot r_i$ ausgedrückt. Als Faustregel für den Redundanzanteile $r_i$ gilt: 
-  * ''0.0 ≤ //r<sub>i</sub>//0.1'' → schlecht kontrolliert, +  * $0.0 \leq r_i \lt  0.1→ schlecht kontrolliert, 
-  * ''0.1 ≤ //r<sub>i</sub>//0.3'' → ausreichend kontrolliert, +  * $0.1 \leq r_i \lt  0.3→ ausreichend kontrolliert, 
-  * ''0.3 ≤ //r<sub>i</sub>//0.7'' → gut kontrolliert, +  * $0.3 \leq r_i \lt  0.7→ gut kontrolliert, 
-  * ''0.7 ≤ //r<sub>i</sub>// ≤ 1.0'' → sehr gut kontrolliert.+  * $0.7 \leq r_i \leq 1.0→ sehr gut kontrolliert.
  
-Die Summe aller ''//r<sub>i</sub>//'' bzw. die Spur der Matrix ''**R**'' liefert den Gesamtfreiheitsgrad ''//f// = tr(**R**)'' der [[:least-squares-adjustment|Ausgleichung]]. Eine hohes Maß an Kontrolle wird erreicht, wenn die Anzahl an Beobachtungen ''//n//'' erhöht wird. Dieser Forderung stehen jedoch vor allem wirtschaftliche Aspekte gegenüber. Mithilfe des Redundanzanteils ''//r<sub>i</sub>//'' lässt sich durch ''//1 - r<sub>i</sub>//'' der Grad der Überschüssigkeit angeben. Hierdurch lässt sich bemessen, welchen Anteil die Beobachtung am Ausgleichungsergebnis mit einbringt.+Die Summe aller $r_i$ bzw. die Spur der Matrix $\mathbf{R}$ liefert den Gesamtfreiheitsgrad $f = \tr \mathbf{R}$ der [[:least-squares-adjustment|Ausgleichung]]. Eine hohes Maß an Kontrolle wird erreicht, wenn die Anzahl an Beobachtungen $nerhöht wird. Dieser Forderung stehen jedoch vor allem wirtschaftliche Aspekte gegenüber. Mithilfe des Redundanzanteils $r_i$ lässt sich durch $1 - r_i$ der Grad der Überschüssigkeit angeben. Hierdurch lässt sich bemessen, welchen Anteil die Beobachtung am Ausgleichungsergebnis mit einbringt.
  
  
 ===== Modellstörung ===== ===== Modellstörung =====
  
-Unter einer Modellstörung ''**∇**'' verstehen wir grundsätzlich einen Fehler im [[:least-squares-adjustment#funktionales_modell|funktionalen Modell]] oder/und [[:least-squares-adjustment#stochastisches_modell|stochastischen Modell]] der Ausgleichung. Eine angezeigte Diskrepanz kann demnach durch das Fehlen eines Zusatzparameters wie bspw. einen Maßstab bei einer Streckengruppe hervorgerufen werden. Gleichzeitig kann aber auch eine zu optimistisch gewählte //a-priori// Unsicherheit für das Anschlagen des [[:least-squares-adjustment:outlier|statistischen Test]] verantwortlich sein. Eine Trennung zwischen beiden Ursachen ist daher nur bedingt möglich.+Unter einer Modellstörung $\mathbf{\nabla}$ verstehen wir grundsätzlich einen Fehler im [[:least-squares-adjustment#funktionales_modell|funktionalen Modell]] oder/und [[:least-squares-adjustment#stochastisches_modell|stochastischen Modell]] der Ausgleichung. Eine angezeigte Diskrepanz kann demnach durch das Fehlen eines Zusatzparameters wie bspw. einen Maßstab bei einer Streckengruppe hervorgerufen werden. Gleichzeitig kann aber auch eine zu optimistisch gewählte //a-priori// Unsicherheit für das Anschlagen des [[:least-squares-adjustment:outlier|statistischen Test]] verantwortlich sein. Eine Trennung zwischen beiden Ursachen ist daher nur bedingt möglich.
  
-Die Größe einer möglichen Modellstörung ''**∇**'' wird formal durch das Erweitern des Parametervektors ''**x**'' und einer Modifikation im funktionalen Modell durch die Designmatrix ''**B**'' ermittelt. In dem erweiterten Ausgleichungsmodell werden demnach ''//j//'' Zusatzparameter ''**∇**'' mitgeschätzt.+Die Größe einer möglichen Modellstörung $\mathbf{\nabla}$ wird formal durch das Erweitern des Parametervektors $\mathbf{x}$ und einer Modifikation im funktionalen Modell durch die Designmatrix $\mathbf{B}$ ermittelt. In dem erweiterten Ausgleichungsmodell werden demnach $jZusatzparameter $\mathbf{\nabla}$ mitgeschätzt.
  
 $$\mathbf{l} + \mathbf{v} = \mathbf{A}\mathbf{\hat{x}} + \mathbf{B}\mathbf{\nabla}$$ $$\mathbf{l} + \mathbf{v} = \mathbf{A}\mathbf{\hat{x}} + \mathbf{B}\mathbf{\nabla}$$
  
-Das Schätzen der Zusatzparameter ''**∇**'' und der zugehörigen Kofaktormatix ''**Q<sub>∇∇</sub>**'' erfordert keine Neuausgleichung und kann direkt aus den Resultaten des nicht-erweiterten Modells abgeleitet werden. Der Lösungsvektor ''**∇**'' ergibt sich aus+Das Schätzen der Zusatzparameter $\mathbf{\nabla}$ und der zugehörigen Kofaktormatix $\mathbf{Q_{\nabla\nabla}}$ erfordert keine Neuausgleichung und kann direkt aus den Resultaten des nicht-erweiterten Modells abgeleitet werden. Der Lösungsvektor $\mathbf{\nabla}$ ergibt sich aus
  
 $$\nabla = -\mathbf{Q_{\nabla\nabla}}\mathbf{B}^T\mathbf{Pv}$$ $$\nabla = -\mathbf{Q_{\nabla\nabla}}\mathbf{B}^T\mathbf{Pv}$$
  
-worin ''**P**'' die Gewichtsmatrix des [[:least-squares-adjustment#stochastisches_modell|stochastischen Modells]], ''**v**'' die Residuen der Beobachtungen und ''**Q<sub>∇∇</sub>**'' die zu ''**∇**'' korrespondierende Kofaktormatix darstellt.+worin $\mathbf{P}$ die Gewichtsmatrix des [[:least-squares-adjustment#stochastisches_modell|stochastischen Modells]], $\mathbf{v}$ die Residuen der Beobachtungen und $\mathbf{Q_{\nabla\nabla}}$ die zu $\mathbf{\nabla}$ korrespondierende Kofaktormatix darstellt.
  
 $$\mathbf{Q_{\nabla\nabla}} = (\mathbf{B}^T\mathbf{PQ_{vv}PB})^{-1} = (\mathbf{B}^T\mathbf{PRB})^{-1}$$ $$\mathbf{Q_{\nabla\nabla}} = (\mathbf{B}^T\mathbf{PQ_{vv}PB})^{-1} = (\mathbf{B}^T\mathbf{PRB})^{-1}$$
  
-Da eine rechen- und zeitintensive Neuausgleichung entfallen kann, lassen sich mit vertretbarem Aufwand verschiedene Störmodelle analysieren und die geschätzten Modellstörungen ''**∇<sub>i</sub>**'' in Relation zur Genauigkeit [[:least-squares-adjustment:outlier|statistisch bewerten]]. JAG3D führt für jede Beobachtung ein spezifisches Störmodell ein. Eine signifikante Modellstörung liegt vor, wenn die Nullhypothese ''**E{∇<sub>i</sub>}**=**0**|H<sub>0</sub>'' verworfen werden muss. +Da eine rechen- und zeitintensive Neuausgleichung entfallen kann, lassen sich mit vertretbarem Aufwand verschiedene Störmodelle analysieren und die geschätzten Modellstörungen $\mathbf{\nabla}_i$ in Relation zur Genauigkeit [[:least-squares-adjustment:outlier|statistisch bewerten]]. JAG3D führt für jede Beobachtung ein spezifisches Störmodell ein. Eine signifikante Modellstörung liegt vor, wenn die Nullhypothese $\rm{E}\{\mathbf{\nabla}_i\} =\mathbf{0}|H_0$ verworfen werden muss. 
  
-Offensichtlich  können nur Beobachtungen geprüft werden, deren Redundanzanteil ''//r<sub>i</sub>//0'' ist, da zur Bestimmung von ''**∇<sub>i</sub>**'' bzw. ''**Q<sub>∇∇,i</sub>**'' die [[#redundanzanteil|Redundanzmatrix]] ''**R**'' notwendig ist. Beobachtungsbezogene Störmodelle können demnach nur berechnet werden, wenn die Beobachtung durch die übrigen Beobachtungen im Netz kontrolliert wird.+Offensichtlich  können nur Beobachtungen geprüft werden, deren Redundanzanteil $r_i \gt 0ist, da zur Bestimmung von $\mathbf{\nabla}_i$ bzw. $\mathbf{Q_{\nabla\nabla}}_i$ die [[#redundanzanteil|Redundanzmatrix]] $\mathbf{R}$ notwendig ist. Beobachtungsbezogene Störmodelle können demnach nur berechnet werden, wenn die Beobachtung durch die übrigen Beobachtungen im Netz kontrolliert wird.
  
  
 ===== Einfluss auf die Punktlage ===== ===== Einfluss auf die Punktlage =====
  
-Um Abzuschätzen, welchen Einfluss eine geschätzte [[#modellstoerung|Modellstörung]] ''**∇<sub>i</sub>**'' auf die berührenden Punkte der korrespondierenden Beobachtung ausübt, ist der Einfluss auf die relative Punktlage ''**EP<sub>i</sub>**'' zur analysieren.+Um Abzuschätzen, welchen Einfluss eine geschätzte [[#modellstoerung|Modellstörung]] $\mathbf{\nabla}_i$ auf die berührenden Punkte der korrespondierenden Beobachtung ausübt, ist der Einfluss auf die relative Punktlage $\mathbf{EP}_i$ zur analysieren.
  
 $$\mathbf{EP}_i=\mathbf{A}_{\mathbf{x},i}\nabla\mathbf{x}_i$$ $$\mathbf{EP}_i=\mathbf{A}_{\mathbf{x},i}\nabla\mathbf{x}_i$$
  
-worin ''**A<sub>x,i</sub>**'' den koordinatenbezogenen Anteil der Designmatrix ''**A**'' des [[:least-squares-adjustment#funktionales_modell|funktionalen Modells]], ''**∇x<sub>i</sub>**'' die Auswirkung einer Modellstörung ''**∇<sub>i</sub>**'' auf die geschätzten Koordinaten beschreiben+worin $\mathbf{A}_{\mathbf{x},i}$ den koordinatenbezogenen Anteil der Designmatrix $\mathbf{A}$ des [[:least-squares-adjustment#funktionales_modell|funktionalen Modells]], $\nabla\mathbf{x}_i$ die Auswirkung einer Modellstörung $\nabla_i$ auf die geschätzten Koordinaten beschreiben
  
 $$\nabla\mathbf{x}_i=\mathbf{Q_{\hat{x}\hat{x}}}{\mathbf{\bar{A_{x}}}}^T{\mathbf{PB}}_i {\mathbf{\nabla}}_i$$ $$\nabla\mathbf{x}_i=\mathbf{Q_{\hat{x}\hat{x}}}{\mathbf{\bar{A_{x}}}}^T{\mathbf{PB}}_i {\mathbf{\nabla}}_i$$
Zeile 110: Zeile 110:
 $$\mathbf{\bar{A_{x}}}=\mathbf{A_x} - \mathbf{A_z}(\mathbf{A_z}^T\mathbf{PA_z})^{-1}\mathbf{A_z}^T\mathbf{PA_x}$$ $$\mathbf{\bar{A_{x}}}=\mathbf{A_x} - \mathbf{A_z}(\mathbf{A_z}^T\mathbf{PA_z})^{-1}\mathbf{A_z}^T\mathbf{PA_x}$$
  
-die um den Anteil der Zusatzunbekannten ''**A<sub>z</sub>**'' reduzierte Designmatrix ''**A**'' sind. Enthält das funktionale Modell für die betreffende Beobachtung keine Zusatzparameter, so ergibt sich ''**EP<sub>i</sub>**'' direkt aus+die um den Anteil der Zusatzunbekannten $\mathbf{A_z}$ reduzierte Designmatrix $\mathbf{A}$ sind. Enthält das funktionale Modell für die betreffende Beobachtung keine Zusatzparameter, so ergibt sich $\mathbf{EP}_i$ direkt aus
  
 $$\mathbf{EP}_i = (\mathbf{E} - \mathbf{R}_i)\mathbf{B}_i\nabla_i$$ $$\mathbf{EP}_i = (\mathbf{E} - \mathbf{R}_i)\mathbf{B}_i\nabla_i$$
  
-wodurch noch einmal die Wichtigkeit eines redundanten Netzes unterstrichen wird. Die Auswirkungen einer geschätzten Modellstörung auf die relative Punktlage hängt somit direkt vom Grad der Überschüssigkeit ''//1 - r<sub>i</sub>//'' ab. Bei einem Redundanzanteil von ''//r<sub>i</sub> = 0.9//'' würden demnach nur ''//10 %//'' der geschätzten Modellstörung ''**∇<sub>i</sub>**'' einen Einfluss auf die relative Punktlage der berührenden Punkte ausüben.+wodurch noch einmal die Wichtigkeit eines redundanten Netzes unterstrichen wird. Die Auswirkungen einer geschätzten Modellstörung auf die relative Punktlage hängt somit direkt vom Grad der Überschüssigkeit $1 - r_i$ ab. Bei einem Redundanzanteil von $r_i = 0.9würden demnach nur 10 % der geschätzten Modellstörung $\nabla_i$ einen Einfluss auf die relative Punktlage der berührenden Punkte ausüben
 + 
 +Der Einfluss auf die relative Punktlage $\mathbf{EP}$ sollte somit auch bei der Entscheidung, ob eine Beobachtung im Datenbestand verbleibt, berücksichtigt werden. Grenzwerte werden u.a. bei hoheitlichen Aufgaben von den Landesvermessungsbehörden festgesetzt. Für Liegenschaftsvermessungen in Baden-Württemberg(([[https://www.lgl-bw.de/lgl-internet/opencms/de/05_Geoinformation/Vorschriften/|Verwaltungsvorschrift für die Durchführung von Liegenschaftsvermessungen (LV-Vorschrift – VwVLV) vom 5. Dezember 2012]])) ist ein Grenzwert von 3 cm vorgesehen. Für Richtungsbeobachtungen oder Zenitwinkel ist $EP$ daher noch in eine Querabweichung umzurechnen.
  
-Der Einfluss auf die relative Punktlage ''**EP**'' sollte somit auch bei der Entscheidung, ob eine Beobachtung im Datenbestand verbleibt, berücksichtigt werden. Grenzwerte werden u.a. bei hoheitlichen Aufgaben von den Landesvermessungsbehörden festgesetzt. Für Liegenschaftsvermessungen in Baden-Württemberg(([[https://www.lgl-bw.de/lgl-internet/opencms/de/05_Geoinformation/Vorschriften/|Verwaltungsvorschrift für die Durchführung von Liegenschaftsvermessungen (LV-Vorschrift – VwVLV) vom 5. Dezember 2012]])) ist ein Grenzwert von ''//3 cm//'' vorgesehen. Für Richtungsbeobachtungen oder Zenitwinkel ist ''**EP**'' daher noch in eine Querabweichung umzurechnen. 
  
 ===== Einflussbemessung der Netzverzerrung ===== ===== Einflussbemessung der Netzverzerrung =====
  
-Die ungünstigste Auswirkung einer nicht-erkannten [[#modellstoerung|Modellstörung]] auf die ausgeglichenen Koordinaten kann durch die skalare Größe ''//EF//'' beschrieben werden. Diese Maßzahl wird häufig als Netzverzerrung oder auch Einflussfaktor (''//EF//'') bezeichnet. +Die ungünstigste Auswirkung einer nicht-erkannten [[#modellstoerung|Modellstörung]] auf die ausgeglichenen Koordinaten kann durch die skalare Größe $EFbeschrieben werden. Diese Maßzahl wird häufig als Netzverzerrung oder auch Einflussfaktor ($EF$) bezeichnet. 
  
 $$EF^2=\frac{1}{\sigma^2_0}\nabla\mathbf{x}^T\mathbf{Q_{\hat{x}\hat{x}}^{-1}}\nabla\mathbf{x}$$ $$EF^2=\frac{1}{\sigma^2_0}\nabla\mathbf{x}^T\mathbf{Q_{\hat{x}\hat{x}}^{-1}}\nabla\mathbf{x}$$
Zeile 128: Zeile 129:
 $$EF \cdot SP = EF \cdot \sigma_{P,max}$$ $$EF \cdot SP = EF \cdot \sigma_{P,max}$$
  
-worin ''//σ<sub>P,max</sub>//'' die maximale mittlere Unsicherheit der ''//k//'' Koordinaten im Netz ist und sich aus+worin $\sigma_{P,max}$ die maximale mittlere Unsicherheit der $kKoordinaten im Netz ist und sich aus
  
 $$\sigma^2_{P,\max} = \max(\tr \mathbf{Q}_{\mathbf{\hat x \hat x}_i}, i=1, 2 ,..., k)$$ $$\sigma^2_{P,\max} = \max(\tr \mathbf{Q}_{\mathbf{\hat x \hat x}_i}, i=1, 2 ,..., k)$$
  
-ergibt. Jäger et al. (2005) führen an, dass diese Einflussbemessung grundsätzlich überschätzt wird. Dennoch findet sich diese Kenngröße in den Verwaltungsvorschriften der Länder und ist bei der Netzbeurteilung mit zu berücksichtigen. In Baden-Württemberg(([[https://www.lgl-bw.de/lgl-internet/opencms/de/05_Geoinformation/Vorschriften/|Verwaltungsvorschrift für die Durchführung von Liegenschaftsvermessungen (LV-Vorschrift – VwVLV) vom 5. Dezember 2012]])) ist ein Grenzwert von ''//15 cm//'' für Liegenschaftsvermessungen noch zulässig.+ergibt. Jäger et al. (2005) führen an, dass diese Einflussbemessung grundsätzlich überschätzt wird. Dennoch findet sich diese Kenngröße in den Verwaltungsvorschriften der Länder und ist bei der Netzbeurteilung mit zu berücksichtigen. In Baden-Württemberg(([[https://www.lgl-bw.de/lgl-internet/opencms/de/05_Geoinformation/Vorschriften/|Verwaltungsvorschrift für die Durchführung von Liegenschaftsvermessungen (LV-Vorschrift – VwVLV) vom 5. Dezember 2012]])) ist ein Grenzwert von 15 cm für Liegenschaftsvermessungen noch zulässig.
least-squares-adjustment/reliability.txt · Zuletzt geändert: 2023/07/05 09:46 von Michael Lösler