Das GEO-Forum - 3D-Fehlerellipsoid, Achsen berechnen und zuordnen

3D-Fehlerellipsoid, Achsen berechnen und zuordnen (Antwort)

Barny.G — Mon, 24 Feb 2014 19:11:20 +0000

Entschuldigung, dass ich mich jetzt erst melde:

Das Ellipsoid liegt richtig. Die Mathematik hat(te) recht. Ich habe mehrere Beispiele gerechnet und im Ergebnis liegen alle so. Wahrscheinlich stimmt was mit meinem Bauchgefühl nicht

3D-Fehlerellipsoid, Achsen berechnen und zuordnen (Antwort)

MichaeL — Tue, 18 Feb 2014 07:57:20 +0000

Hi,

Möglicherweise ist es aber "nur" ein MatLab-Problem, wie es seb. schreibt:
4555

Validier es doch mit einer bekannten Lösung. Im Benning ist eine 3x3 Kovarianzmatrix drin, die mittels spektraler Analyse zerlegt wird. Hier könntest Du mal Deine ermittelten Werte mit denen aus dem Buch vergleichen.

Gruß Micha

3D-Fehlerellipsoid, Achsen berechnen und zuordnen (Antwort)

.seb — Mon, 17 Feb 2014 13:52:55 +0000

Hat funktioniert,

was du uns nicht verraten hast, dass zwei unterschiedliche Ellipsoide rauskommen, jenachdem ob man die Matrizen flippt oder nicht. Da kann was nicht stimmen. Der Ellipsoid ohne fliplr sieht mir noch am ehesten wie ein Fehlerellipsoid aus. Das andere... na ja... ist das überhaupt so vorgesehen.

Bevor wir weiter philosophieren, würde ich ein richtiges Matlab-Forum aufsuchen und dort dein Programm vorstellen. Denn es handelt sie allen Anschein nach um ein Matlab-Verarbeitungsproblem und kein Geodäsie Spezifisches.

Grüße
Sebastian

3D-Fehlerellipsoid, Achsen berechnen und zuordnen (Antwort)

Barny.G — Mon, 17 Feb 2014 13:08:12 +0000

OK, mal sehen ob das funktioniert - die gespeicherten Variablen sind hier:

http://www.workupload.com/file/QeUkbPSq

und das skript dazu hier:

http://www.workupload.com/file/XrWpP01v

Viel Spaß beim Testen...

Thomas.

3D-Fehlerellipsoid, Achsen berechnen und zuordnen (Antwort)

Barny.G — Mon, 17 Feb 2014 13:02:11 +0000

Habe mich gerade noch mal reingelesen, wenn du [V D] = eig() sind auf der Diagonalen von D (diag(D)) die Eigenwerte

... und die sollten mit den EV in V korrespondieren - oder nicht?

und die liegen ungekehrt vor als wenn du nur D = eig() eingibst!

Hm, bei meinem Bsp nicht, aber egal. Die Zuornung ist es ja, die wichtig ist. Denn diese bestimmt ja schließlich die Einheitenlänge der jeweiligen EV...

Bei meinem Versuch habe ich ja gegen meine innere Überzeugung die EV getauscht und die EW beibehalten, was diese o.g. Zuordnung verändert. Und das, glaube ich, ist falsch.