Die sytematische
Suche nach den optimalen Parametern, mit denen man die von S1
aus gemessenen Werte möglichst nahe an die Werte "herantransformieren"
kann, die von S0 bestimmt wurden, läuft nach der sog. Methode
der kleinsten Quadrate. Dabei findet man für die folgenden
Größen je einen Zahlenwert:
-
Maßstabsfaktor S
-
Drehwinkel
um die x-Achse rx
- Drehwinkel
um die y-Achse ry
- Drehwinkel
um die z-Achse rz
- Verschiebung
in x-Richtung DX
- Verschiebung
in y-Richtung DY
- Verschiebung
in z-Richtung DZ
Die Methode
verläuft über einen Näherungsprozess und man beginnt
mit möglichst plausiblen Anfangs-Schätzwerten für
die obigen Parameter. Im anstehenden Fall komme ich auf folgendem
Weg zu solchen Startwerten:
Ich kann davon
ausgehen, dass mein Tachymeter zuverlässige Distanzmessungen
liefert und außerdem ist anzunehmen, dass beim Wechseln
von Position S0 nach Position S1 die Abmessungen des Gebäudes
unverändert bleiben. Es ist also vernünftig, wenn ich
von einem Maßstabsfaktor 1.000 als Näherung für
den S-Wert ausgehe.
Wegen der
sorgfältigen Horizontierung des Instruments in beiden Positionen
nehme ich weiterhin für rx und ry die Werte 0° an. Für
rz nehme ich den schon erwähnten Schätzwert von 91.7°
an. Für DX gebe ich die geschätzten 45.21m, für
DY= -0.57m und für DZ=1.05m vor (siehe
letzte Seite).
Im Verlauf
des Näherungsprozesses konvergieren die vorgegebenen 7 Parameter
- wenn mir das Glück hold ist, oder anders ausgedrückt,
wenn es sich bei meinen Passpunkten um ein "gutmütiges
Punktesystem" handelt - hin zu Endwerten, die sich bei weiteren
Iterationsschritten kaum mehr verändern. Es erweist sich,
dass eine Transformation mit diesen Parameterendwerten zu einer
optimalen Einpassung der S1-Punkte in das S0-System führt,
in dem Sinn, dass die Summe der quadrierten Abweichungen minimal
wird.
Nachfolgend
können Sie sich ansehen, wie der Vektor mit den ursprünglichen
Schätzwerten zu einem Vektor mit optimierten Parameterwerten
hingeführt wird.
|
Parameterart
|
Startwerte
der 7 Parameter
|
Endwerte
der 7 Parameter
|
|
S
(Maßstab)
|
1
|
1.00084
|
|
rx
(Drehg. um x-Achse)
|
0°
|
0.12°
|
|
ry
(Drehg. um y-Achse)
|
0°
|
-0.19°
|
|
rz
(Drehg. um z-Achse)
|
91.7°
|
+84.61°
(nach links)
|
|
DX
(Verschiebg. in x-Richtung)
|
45.21m
|
43.6075m
|
|
DY
(Verschiebg. in y-Richtung)
|
-0.57m
|
-4.84487m
|
|
DZ
(Verschiebg. in z-Richtung)
|
1.05m
|
1.12189m
|
Wenn
man die sieben Parameter-Endwerte und die kartesischen Koordinaten
der von S1 aus gemessenen Punkte gemeinsam in den Vorgang der
Transformation einbringt, dann erhält man für die gemessenen
Punkte die zum S0-System passenden kartesischen Koordinaten.
Für
eigene Versuche und zur Nachprüfung obiger Resultate finden
Sie hier ein Applet als Rechenhilfe
Wenn
Sie mit einer ganzen Reihe von Punkten und mit einer selbst getroffenen
Auswahl von Transformationsparametern eine Helmerttransformation
berechnen lassen wollen,
dann
kann ich Ihnen dieses Applet zum Ausprobieren vorstellen.
|