Funktionsweise von "Transformieren" (Coverage)

Beschreibung

Mit dem Werkzeug Transformieren werden Daten aus einem Koordinatensystem in ein anderes konvertiert. Oftmals wird es verwendet, um in Digitalisiereinheiten (meistens Zoll) erstellte Coverages in reale, auf dem Kartenmanuskript darstellbare Einheiten zu konvertieren. Diese realen Orte werden anhand einer bestimmten Projektion dargestellt. So kann es sich bei der Projektion Ihrer Grundkarte beispielsweise um UTM mit in Metern gemessenen Passpunktpositionen handeln.

Um eine Transformation durchzuführen, müssen Sie in Ihrem Coverage mindestens zwei Tics für reale Orte registrieren können. Diese realen Orte können mithilfe des Grundkartenkoordinatensystems in Fuß oder in Meter gemessen werden. Sie können sich auch an Schnittstellen von Längen-/Breitengraden befinden. Wenn Ihre Orte mit Längen- und Breitengraden angegeben sind, können Sie mit dem Werkzeug Projizieren die Passpunkte in das Koordinatensystem konvertieren, das Ihrer in Fuß oder Meter gemessenen Grundkarte entspricht. Führen Sie dies vor dem Transformieren des digitalisierten Coverages aus.

Wenn Sie die Daten zum ersten Mal transformieren, führen Sie eine Transformation in dieselbe Projektion wie die Quellkarte aus. Sie können anschließend mit dem Werkzeug "Projizieren" das Coverage in eine andere Kartenprojektion konvertieren.

Abbildung "Transformieren"

Das Werkzeug "Transformieren" basiert auf einem Vergleich der Koordinaten von Tics im Eingabe-Coverage mit den Koordinaten entsprechender Tics im Ausgabe-Coverage. Die zu vergleichenden Tics werden durch die Tic-IDs angegeben. Die Mindestanzahl benötigter Tics hängt von der jeweils verwendeten Transformationsoption ab. Tics werden bei der Transformation als Passpunkte behandelt und nicht mit der Transformationsgleichung konvertiert. Die Koordinaten der Ausgabe-Tics sind nach der Transformation unverändert. Features, die mit den Eingabe-Tics übereinstimmen, stimmen jedoch in der Regel nicht mit den Ausgabe-Tics überein. Möglicherweise müssen Sie Features anpassen, damit sie nach der Transformation mit den Tics im Ausgabe-Coverage übereinstimmen.

Mit dem Werkzeug "Transformieren" wird standardmäßig eine Funktion für die affine Transformation verwendet, um Coverage-Koordinaten aus dem Eingabe-Coverage in das Ausgabe-Coverage zu transformieren. Zur Berechnung einer affinen Transformation sind mindestens drei Tics erforderlich.

Bei einer affinen Transformation können die Daten unterschiedlich skaliert, verzerrt, gedreht und übertragen werden. In der folgenden Grafik sind die vier möglichen Änderungen dargestellt.

Affine Transformation

Die Funktion für die affine Transformation lautet wie folgt:

x’ = Ax + By + C

Die Parameter A-F werden durch die Passpunkte festgelegt,

y’ = Dx + Ey + F

und dann auf alle Koordinaten im Coverage angewendet.

Dabei sind x und y Koordinaten des Eingabe-Coverages und x’ und y’ Koordinaten des Ausgabe-Coverages. A, B, C, D, E und F werden durch einen Vergleich der Positionen von Tics im Eingabe-Coverage und ihren Positionen im Ausgabe-Coverage bestimmt. Mit ihnen werden die Tic-Koordinaten skaliert, verzerrt, gedreht und übertragen.

Das Werkzeug "Transformieren" gibt die sechs Parameter und ihre geometrische Interpretation an. Es folgt ein Beispiel für einen Teil des ausgegebenen Berichts zu einer affinen Transformation.

Transformieren von Koordinaten für Coverage "in_cov_name"

Scale (X,Y) = (246.140,255.702)  Skew (degrees) = (-0.061)
Rotation (degrees) = (0.334)  Translation = (2890.267,3679.906)
RMS Error (input, output) = (0.084,20.592)

Affine X = Ax + By + C
       Y = Dx + Ey + F
A =   246.135   B =   -1.763        C =   2890.267
D =   1.434     E =   255.696       F =   3679.906

In einem Abschnitt des Berichts sind die Parameter der Transformationsgleichungen aufgelistet. Sie können die Gleichungsparameter interpretieren, um zu verstehen, was in geometrischer Hinsicht geschieht. Der erste Abschnitt des Berichts gibt die interpretierten Werte für Skalierung, Verzerrung und Rotation an. Durch eine Änderung der Reihenfolge, in der diese geometrischen Aktionen angewendet werden, werden deren interpretierte Werte geändert, doch sind das Endergebnis und die Gleichungsparameter immer gleich. Die Interpretation wird in der Reihenfolge (1) Skalierung, (2) Verzerrung in X-Richtung und (3) Rotation von Punkten um die Achsen entgegen dem Uhrzeigersinn durchgeführt. Die Übertragungswerte werden nicht interpretiert, sondern sind einfach die Übertragungsparameter C und F aus den Gleichungen. Mit dem Werkzeug "Transformieren" werden die Parameter der affinen Gleichung wie folgt interpretiert:

A = mx · cos t
B = my · (k · cos t - sin t)
D = mx · sin t
E = my · (k · sin t + cos t)
C = translation in x direction
F = translation in y direction

where

mx = change of scale in x direction
my = change of scale in y direction
k = shear factor along the x-axis = tan (skew angle) (the skew angle is measured from the y-axis)
t = rotation angle, measured counterclockwise from the x-axis

Der Rotationswinkel bezieht sich auf die Rotation der Daten in einem festen XY-Koordinatensystem. In der vorherigen Abbildung wurden die ursprünglichen Daten um 45 Grad gedreht. Ein negativer Wert bedeutet, dass die Daten in Bezug auf die X-Achse im Uhrzeigersinn gedreht werden.

Mit der Ähnlichkeitstransformation werden die Daten skaliert, gedreht und übertragen. Die Achsen werden nicht unabhängig skaliert. Außerdem wird keine Verzerrung vorgenommen. Für die Ähnlichkeitstransformation werden mindestens zwei Passpunkte benötigt. Im folgenden Bericht wird eine Ähnlichkeitstransformation beschrieben. Es wurden dieselben Coverages verwendet, die zuvor mit einer affinen Transformation transformiert wurden.

Transformieren von Koordinaten für Coverage "in_cov_name"

Scale (X,Y) = (249.927,249.927)
Rotation (degrees) = (0.362)  Translation = (2855.407,3715.168)
RMS Error (input, output) = (0.118,29.398)

Similarity X =  Ax + By + C
           Y = -Bx + Ay + F
A =   249.922         B =   -1.578
C =   2855.407        F =   3715.168

Beachten Sie, dass die X- und Y-Maßstabswerte gleich sind. Ein Verzerrungswert ist nicht angegeben, da bei einer Ähnlichkeitstransformation eine nicht orthogonale Verschiebung der Daten unzulässig ist.

Die Funktion der Ähnlichkeitstransformation lautet

x’ =  Ax + By + C
y’ = -Bx + Ay + F

where

A = s · cos t
B = s · sin t
C = translation in x direction
F = translation in y direction

and

s = scale change (same in x and y directions)
t = rotation angle, measured counterclockwise from the x-axis

Wie bei der Erläuterung zur affinen Transformation bezieht sich der Rotationswinkel auf die Rotation der Daten in einem festen XY-Koordinatensystem. In der vorherigen Abbildung wurden die ursprünglichen Daten um 45 Grad gedreht. Ein negativer Wert bedeutet, dass die Daten, von der X-Achse aus gemessen, im Uhrzeigersinn gedreht werden.

Die projektive Transformation basiert auf einer komplexeren Formel, für die mindestens vier Tics erforderlich sind.

x’ = (Ax + By + C) / (Gx + Hy + 1)
y’ = (Dx + Ey + F) / (Gx + Hy + 1)

Der Bericht zu einer projektiven Transformation beinhaltet den ungefähren Maßstab, den RMS-Fehler und die Gleichungsparameter. Aufgrund der Komplexität der Interpretation werden die Parameter der projektiven Transformation nicht interpretiert. Weitere Informationen finden Sie unter den photogrammetrischen Texten, die in den Referenzen am Ende dieses Themas aufgeführt sind. Es folgt ein Beispiel eines ausgegebenen Berichts zu einer projektiven Transformation.

Transformieren von Koordinaten für Coverage "in_cov_name"

Approximate scale = 1479.087
RMS Error (input, output) = (0.040,60.878)

2D Projective  X = (Ax + By + C) / (Gx + Hy + 1)
               Y = (Dx + Ey + F) / (Gx + Hy + 1)

A =         55.667   B =       -718.999   C =    2125052.558
D =       -199.525   E =       1385.541   F =     317759.475
G =         -0.001   H =          0.000

Principal point of input (xp,yp) = (2.000,16.946)
Exposure center of output(Xc,Yc) = (2127791.000,343183.000)

Der Hauptpunkt und der Lagemittelpunkt entsprechen dem ersten Tic in der Tic-Datei des Eingabe-Coverages. Wenn Sie einen anderen Punkt verwenden möchten (z. B. den Mittelpunkt des Datasets), stellen Sie sicher, dass es sich um den ersten Tic in der Datei handelt.

Ein RMS-Fehler (Root Mean Square) wird für jede ausgeführte Transformation berechnet. Er gibt an, wie gut die abgeleitete Transformation ist. Im folgenden Beispiel ist die relative Position von vier Ausgabe-Passpunkten (Tics) und den transformierten Eingabe-Passpunkten dargestellt:

Abbildung "RMS"

Der RMS-Fehler misst die Fehler zwischen den Tics des Ausgabe-Coverages und den transformierten Positionen der Tics des Eingabe-Coverages.

Die Transformation wird anhand der Methode der kleinsten Quadrate abgeleitet, sodass mehr Tics als nötig erstellt werden können.

Formel für RMS

Der RMS-Fehler für die einzelnen Transformationen wird sowohl in den Eingabe-Coverage-Einheiten (z. B. 0,031 Zoll) als auch in den Ausgabe-Coverage-Einheiten (z. B. 37,465 Fuß) angegeben.

RMS Error (input, output) = (0.031,37.465)

Eine perfekte Transformation ergibt einen RMS-Fehler von null. Sie werden zwar keinen Wert 0,000 erreichen, sollten den RMS jedoch so gering wie möglich halten. Es ist empfehlenswert, einen maximal zulässigen Fehler festzulegen und beizubehalten. Zulässige Werte sind je nach Genauigkeit der Ursprungsdaten und dem Maßstab der Quellkarte verschieden. Hohe RMS-Fehler zeigen an, dass die alten und die neuen Tics nicht denselben relativen Positionen entsprechen. Wenn der RMS-Fehler höher als der festgelegte Wert ist, kann dies auf Konvertierungsprobleme hinweisen, die zu einer falschen Streckung des Ausgabe-Coverages führen können.

Im letzten Teil des Berichts sind die Eingabe- und Ausgabe-Tic-Koordinaten sowie ihre X- und Y-Fehler aufgelistet (siehe unten). Die Transformationsparameter sind der bestmögliche Kompromiss zwischen den Eingabe- und Ausgabe-Tics. Wenn Sie mit den Transformationsparametern die tatsächlichen Eingabe-Tics transformieren, stimmen die transformierten Ausgabepositionen nicht mit den tatsächlichen Positionen der Ausgabe-Tics überein. Die X- und Y-Fehler sind eine Maßeinheit für die Übereinstimmung der tatsächlichen Positionen mit den transformierten Positionen der Ausgabe-Tics. Sie können manchmal einen Fehler entdecken, indem Sie sich die X- und Y-Fehler ansehen.

tic id  input x   input y   output x      output y     x error   y error
------------------------------------------------------------------------
1       2.000     16.946    2127791.000   343183.000    14.463    75.499
2       12.764    16.821    2143469.000   343326.000   -31.043   -85.363
3       2.052     1.976     2128000.000   320680.000   -36.290    -2.353
4       12.922    2.013     2143729.000   320912.000    20.245    -6.163
5       2.082     9.442     2127944.000   332015.000    22.016   -74.699
6       12.662    9.442     2143320.000   332015.000    10.609    93.079

Beispiel

In diesem Beispiel wird ein Coverage aus Einheiten des Digitalisierungsgeräts in reale Koordinaten transformiert. Das ursprüngliche Coverage mit dem Namen INCHCOV weist die folgenden Passpunkte auf, die in Zoll des Digitalisierungsgeräts gemessen wurden:

IDTIC    XTIC       YTIC
1       2.000     16.946
2      12.764     16.821
3       2.052      1.976
4      12.922      2.013
5       2.082      9.442
6      12.662      9.442

Für die Grundkarte wurde das State-Plane-Koordinatensystem verwendet. Passpunkte wurden in Fuß gemessen. Jedes Tic entspricht den folgenden Positionen:

Tic-IDs	 X Coordinates	 Y Coordinates
1	       2,127,791     	343,183
2       	2,143,469     	343,326
3       	2,128,000     	320,680
4       	2,143,729     	320,912
5       	2,127,944     	332,015
6	       2,143,320     	332,015

Das Ausgabe-Coverage muss vorhanden sein, bevor das Werkzeug "Transformieren" ausgeführt werden kann. Mit dem Werkzeug Coverage erstellen können Sie das Ausgabe-Coverage erstellen und die TIC- und BND-Dateien gegebenenfalls aus dem Eingabe-Coverage kopieren. Optional können Sie mit dem Werkzeug Generate ein neues Coverage erstellen, das Tics an den gewünschten Positionen enthält.

Die TIC-Datei für das Ausgabe-Coverage muss XY-Koordinaten für jedes Tic enthalten, das im Ausgabe-Coverage beibehalten werden soll. Dies sind in diesem Fall die TIC-Positionen im State-Plane-Koordinatensystem. Tic-Koordinaten können in ArcMap oder in den Modulen TABLES oder INFO von ArcInfo Workstation angepasst werden.

Sobald die neuen Tics erstellt wurden, können die Features im digitalisierten Coverage in reale Einheiten transformiert werden. Das Werkzeug "Transformieren" wird verwendet, um ein State-Plane-Projektions-Coverage in Zoll (INCHCOV) in ein State-Plane-Projektions-Coverage in Fuß (STATECOV) zu konvertieren. Es folgt ein Beispiel für die affine Transformation.

Scale (X,Y) = (1452.317,1508.433)  Skew (degrees) = (0.416)
Rotation (degrees) = (0.218) Translation = (2124994.654,317664.385)
RMS Error (input, output) = (0.048,71.614)

Affine  X = Ax + By + C
        Y = Dx + Ey + F
A =     1452.230   B =       -5.526   C =   2124994.654
D =       15.858   E =     1508.462   F =    317664.385

tic id  input  x      input  y
        output x      output y   x error   y error
--------------------------------------------------
1          2.000        16.946
     2127791.000    343183.000    14.463    75.499
2         12.764        16.821
     2143469.000    343326.000   -31.043   -85.363
3          2.052         1.976
     2128000.000    320680.000   -36.290    -2.353
4         12.922         2.013
     2143729.000    320912.000    20.245    -6.163
5          2.082         9.442
     2127944.000    332015.000    22.016   -74.699
6         12.662         9.442
     2143320.000    332015.000    10.609    93.079

Es folgt das entsprechende Beispiel einer projektiven Transformation für dasselbe Coverage. Denken Sie daran, dass die Option PROJECTIVE normalerweise nur verwendet wird, wenn die ursprünglichen Coverages direkt aus Luftaufnahmen digitalisiert wurden.

Transformieren von Koordinaten für Coverage "INCHCOV"

Approximate scale = 1479.087
RMS Error (input, output) = (0.040,60.878)

2D Projective  X = (Ax + By + C) / (Gx + Hy + 1)
               Y = (Dx + Ey + F) / (Gx + Hy + 1)
A =         55.667   B =       -718.999   C =    2125052.558
D =       -199.525   E =       1385.541   F =     317759.475
G =         -0.001   H =          0.000

Principal point of input (xp,yp) = (2.000,16.946)
Exposure center of output(Xc,Yc) = (2127791.000,343183.000)

tic id  input  x      input  y
        output x      output y   x error   y error
--------------------------------------------------
1         2.000         16.946
    2127791.000     343183.000    -4.438    45.252
2        12.764         16.821
   2143469.000      343326.000   -11.447   -36.202
3         2.052          1.976
   2128000.000      320680.000   -17.300    46.421
4        12.922          2.013
   2143729.000      320912.000     1.704   -36.962
5         2.082          9.442
   2127944.000      332015.000    21.787   -93.410
6        12.662          9.442
   2143320.000      332015.000     9.694    74.901

Zum Schluss folgt dasselbe Coverage, das mit der Option SIMILARITY transformiert wurde.

Transformieren von Koordinaten für Coverage "inchcov"

Scale (X,Y) = (1483.794,1483.794)
Rotation (degrees) = (0.377)  Translation = (2124800.900,317942.729)
RMS Error (input,output) = (0.162,240.958)

Similarity  X =  Ax + By + C
Y = -Bx + Ay + F
A =         1483.762   B =           -9.765
C =      2124800.900   F =       317942.729

tic id  input  x      input  y
        output x      output y   x error   y error
--------------------------------------------------
1          2.000        16.946
     2127791.000    343183.000  -188.053   -76.916
2         12.764        16.821
     2143469.000    343326.000   106.378  -300.277
3          2.052         1.976
     2128000.000    320680.000  -173.717   214.680
4         12.922         2.013
     2143729.000    320912.000   225.411   143.724
5          2.082         9.442
     2127944.000    332015.000  -146.109   -42.262
6         12.662         9.442
     2143320.000    332015.000   176.089    61.051

Im folgenden Beispiel wird ein Coverage mit einer affinen Transformation umgekehrt.

Transformieren von Koordinaten für Coverage "flag"

Scale (X,Y) = (1.000,-1.000)  Skew (degrees) = (0.000)
*** Negative Y scaling indicates reflection around X axis. ***
Rotation (degrees) = (180.000)  Translation = (800.000,0.000)
RMS Error (input, output) = (0.000,0.000)

Affine  X = Ax + By + C
        Y = Dx + Ey + F
A =       -1.000   B =        0.000   C =      800.000
D =        0.000   E =        1.000   F =        0.000

tic id  input  x      input  y
        output x      output y   x error   y error
--------------------------------------------------
1        700.000      100.000
         100.000      100.000      0.000     0.000
2        700.000      800.000
         100.000      800.000      0.000     0.000
3        100.000      800.000
         700.000      800.000      0.000     0.000
4        100.000      100.000
         700.000      100.000      0.000     0.000

Beachten Sie die Meldung zur Reflexion im ausgegebenen Bericht. In der folgenden Grafik ist dargestellt, was mit dem Coverage geschieht. Das Rechteck mit der Beschriftung I ist das Eingabe-Coverage. Es wird ein Y-Maßstab von -1 auf die Koordinaten angewandt. Dies führt zum Rechteck II. Es ist kein Verzerrungsfaktor vorhanden, sodass die Punkte um 180° gedreht und dadurch an Position III gesetzt werden. Der abschließende Schritt (nicht dargestellt) ist es, das Coverage zu übertragen.

Abbildung "Transformieren"

Referenzliste

Maling, D.H. Coordinate Systems and Map Projections. George Philip., 1973.

Maling, D.H. "Coordinate systems and map projections for GIS." In: Maguire D.J., M.F. Goodchild, and D.W. Rhind (eds.), Geographical Information Systems: principles and applications. Vol. 1, S. 135-146. Longman Group UK Ltd., 1991.

Moffitt, F.H., and E.M. Mikhail, Photogrammetry. Dritte Ausgabe. Harper and Row, Inc., 1980.

Pettofrezzo, A.J, Matrices and Transformations. Dover Publications, Inc., 1966.

Slama, C.C., C. Theurer, and S.W. Henriksen (eds.), Manual of Photogrammetry. 4. Auflage. Kapitel XIV, S. 729-731. ASPRS, 1980.

9/11/2013