Calibration systems can be created in this section which allow the projected coordinates to be transformed into the coordinates of a specific local system.
Subsequently, they can be selected in the active project to work in the local system defined.
Based on the traditional coordinates' transformation methods and in view of their size, a classification has been carried out which includes said methods in addition to a combination of some of them.
3 transformation units are presented:
2D : Includes the 4-parameter Helmert transformation and a coordinates translation in x and y.
2D Translations
Calculates the x and y displacements by means of the average of the differences between origin and destination. Only a pair of points is required.
Formulas:
where:
x', y' = Transformed x,y coordinates
x,y = Original x,y coordinates
Tx = Translation x
Ty = Translation y
2D-Helmert
It is also known as the 4-parameter similarity transformation. The transformation process includes 3 steps: scaling, rotation and translations.
The first two are defined by one parameter each and the translations include:
At least two pairs of points are required.
Formulas:
where:
x', y' = Transformed x,y coordinates
x,y = Original x,y coordinates
S = Scale
= Rotation angleTx = Translation x
Ty = Translation y
3D : This group is formed by the 7-parameter Helmert transformation and a coordinates translation in x, y and z.
3D Translations
This type of transformation calculates the x, y and z displacements by means of the average of the differences between origin and destination. Only a pair of points is required.
Formulas:
where:
x', y', z' = Transformed x,y,z coordinates
x,y, z = Original x,y,z coordinates
Tx = Translation x
Ty = Translation y
Tz = Translation z
3D-Helmert
It is also known as the transformation of 7 parameters. The parameters involved are: three rotations, three translations and a scale factor. The rotation matrix is constructed by means of three consecutive turns around the alignments x, y, z. At least 3 pairs of points must be available.
Formulas:
where:
x', y', z' = Transformed x,y,z coordinates
x,y, z = Original x,y,z coordinates
S = Scale
Tx, Ty, Tz = Translations in x,y,z
m11...m33 = Rotation matrix coefficients
2D + 1D : 3 dimensional transformations which use 2D-Helmert for planimetry and translation in z or translation in z and slopes x and y for altimetry.
2D-Helmert + Vertical Displacement
This type of transformation is a combination of the transformation of 4-parameter Helmert and a displacement in z. This translation is the average of the differences of this component between the points of origin and destination.
Formulas:
where:
x', y' = Transformed x,y coordinates
x,y = Original x,y coordinates
S = Scale
= Rotation angle
Tx = Translation x
Ty = Translation y
Tz = Translation z
2D-Helmert + Vertical Displacement and Slopes
As with the previous transformation, this is a combination of the 4-parameter Helmert transformation and a translation in z calculated in line with a displacement and inclinations of the x and y components. At least two pairs of points are required with x and y coordinates and three with the z component.
Formulas:
where:
x', y' = Transformed x,y coordinates
x,y = Original x,y coordinates
S = Scale
= Rotation angle
Tx = Translation x
Ty = Translation y
Tz = Translation z
Px = Slope x
Py = Slope y
X0 = Original x coordinate of the first pair of points
Y0 = Original y coordinate of the first pair of points
To manage the local systems the following options' menu is presented:
Two modes are distinguished between to define the transformations:
By Parameters : Directly entering the values of the transformation equation parameters. Displacements must be entered in metres and the angle format must be sexagesimal in degrees, minutes and seconds (separated by space) and in an anti-clockwise direction with 0 in the East.
By Points : By entering a series of origin and destination coordinate pairs, 2D or 3D, and calculating the equations' systems of the selected method.
It allows the entering of the pairs of points from a file by means of the Origin File and Destination File buttons or manually by pressing the Insert button.
The files must be of the point (*.PUN) or base (*.BSE) type. For origin data geographic coordinates' files are also permitted with an extension *.W84 created under the option File Management > GPS Data.
If a point is entered manually the Name, Coordinates Origin and Destination must be indicated and whether it will be used for planimetry, Control H, altimetry, Control V or in both cases. It is also possible to edit the values of these parameters for any point, whatever the insertion method.
The programme endeavours to resolve the selected transformation equations' system every time any change is made to the list of pairs of points, either by inserting, editing or eliminating pairs. If a solution is found, the calculated values will be shown as well as the standard deviation of each parameter. The angle format is sexagesimal in degrees, minutes and seconds and in an anti-clockwise direction with 0 in the East. The displacements are shown in metres.
By contrast, if the MSE button is pressed, the screen is accessed which shows the mean square error and the maximum residuals of each x, y, z component. Depending on the size of the transformation, the following square errors are calculated:
2D : Mean square error H (x, y).
3D : Mean square error 3D (x,y,z), H (x, y), V(z).
2D + 1D : Mean square error H (x, y), V(z).
The possibility is also offered of creating a local system report through the option Create HTML Report from the context menu.
An example of this has been shown below:
Transformation Data
|
File
|
D:\Transform\Doc\H2D.htm
|
|
Type
|
Helmert (4-parameter similarity transformation) (2D)
|
|
Date and Time
|
07/04/2012 13:27:33
|
Parameters
|
TX
|
1050003.715 ± 0.123
|
|
TY
|
50542.131 ± 0.123
|
|
Rotation
|
176°46'54.97952'' ± 0°1'58.89222''
|
|
Scale
|
4.51962 ± 0.00058
|
Statistics
|
ECM
|
0.1142
|
|
Standard Deviation
|
0.140
|
|
Max X Residual
|
0.105 (Control Point 3)
|
|
Max Y Residual
|
0.106 (Control Point 3)
|
Control Points
|
Number of Control Points
|
3
|
|
Control Point
|
Used
|
X Origin
|
Y Origin
|
Z Origin
|
X Destination
|
Y Destination
|
Z Destination
|
X Residual
|
Y Residual
|
|
1
|
Yes
|
121.622
|
-128.066
|
50.000
|
1049422.400
|
51089.200
|
30.000
|
0.004
|
-0.029
|
|
2
|
Yes
|
141.228
|
187.718
|
40.000
|
1049413.950
|
49659.300
|
20.000
|
0.101
|
-0.077
|
|
3
|
Yes
|
175.802
|
135.728
|
30.000
|
1049244.950
|
49884.950
|
10.000
|
-0.105
|
0.106
|
Once the parameters have been calculated or entered manually, the accept button must be pressed to create the local system. It is an ASCII file with the extension NTR and with the following format (example of 2D-Helmert transformation):
[Type]
2D Helmert (4-parameter similarity transformation)
[Parameters]
1050003.71454
50542.13112
3.085426889
4.519620520
[Inverse Parameters]
232582.62262
-1876.44668
-3.085426889
0.221257507
[MSE]
3D 0.0000 H 0.1142 V 0.0000
[Inverse MSE]
3D 0.0000 H 0.0253 V 0.0000
[Standard Deviation]
0.140
[Max X Residual]
0.105 (Control Point 3)
[Max Y Residual]
0.106 (Control Point 3)
[Control Points]
3
[Point Used HControl VControl XSource YSource ZSource XTarget YTarget ZTarget XRes YRes ZRes]
1 1 1 1 121.622 -128.066 50.000 1049422.400 51089.200 30.000 0.004 -0.029 ?
2 1 1 1 141.228 187.718 40.000 1049413.950 49659.300 20.000 0.101 -0.077 ?
3 1 1 1 175.802 135.728 30.000 1049244.950 49884.950 10.000 -0.105 0.106 ?
After creating the transformation file, it will be necessary to go to the option Project > Properties and select it from the list Local System. If a different directory has been recorded from that of the current project, a copy must be made to the latter using the option Local System > Import. This option affords the possibility of activating it directly without having to do so in project properties.
