CIE 1960 color space: Difference between revisions
en>Helpful Pixie Bot m ISBNs (Build KC) |
en>Monkbot |
||
Line 1: | Line 1: | ||
The '''Helmert transformation''' (named after [[Friedrich Robert Helmert]], 1843–1917; also called a '''seven-parameter transformation''') is a transformation method within a three-dimensional space. It is frequently used in [[geodesy]] to produce distortion-free transformations from one [[datum (geodesy)|datum]] to another using: | |||
: <math>X_T=C+\mu RX \, </math> | |||
where | |||
* ''X''<sub>''T''</sub> is the transformed vector | |||
* ''X'' is the initial vector | |||
The [[parameter]]s are: | |||
*<math>C</math> — translation vector. Contains the three [[Translation (geometry)|translations]] along the coordinate axes | |||
*<math>\mu</math> — [[scale factor]], which is unitless, and as it is usually expressed in [[parts per million|ppm]], it must be divided by 1,000,000. | |||
*<math>R</math> — [[rotation matrix]]. Consists of three axes (small [[rotation]]s around the coordinate axes) <math>r_x</math>, <math>r_y</math>, <math>r_z</math>. The rotation matrix is an [[orthogonal matrix]]. The rotation is given in [[radian]]s. | |||
Thus, the Helmert transformation is a [[similarity transformation (geometry)|similarity transformation]]. | |||
== Calculating the parameters == | |||
If the transformation parameters are unknown, they can be calculated with reference points (that is, points whose coordinates are known before and after the transformation. Since a total of seven parameters (three translations, one scale, three rotations) have to be determined, at least two points and one coordinate of a third point (for example, the Z-coordinate) must be known. This gives a [[system of linear equations]] with seven equations and seven unknowns, which can be solved. | |||
In practice, it is best to use more points. Through this correspondence, more accuracy is obtained, and a statistical assessment of the results becomes possible. In this case, the calculation is adjusted with the Gaussian [[least squares]] method. | |||
A numerical value for the accuracy of the transformation parameters is obtained by calculating the values at the reference points, and weighting the results relative to the [[centroid]] of the points. | |||
While the method is mathematically rigorous, it is entirely dependent on the rigour of the parameters that are used. In practice, these parameters are computed from the inclusion of at least three known points in the networks. However the accuracy of these will affect the following transformation parameters, as these points will contain observation errors. Therefore a "real-world" transformation will only be a best estimate and should contain a statistical measure of its quality. | |||
It is not always necessary to use the seven parameter transformation, sometimes it is sufficient to use the '''five parameter transformation''', composed of three translations, one rotation (about the Z-axis) and one change of scale. | |||
==Two-dimensional case== | |||
A special case is the two-dimensional Helmert transformation. Here, only four parameters are needed (two translations, one scaling, one rotation). These can be determined from two known points; if more points are available then checks can be made. | |||
==Application== | |||
The Helmert transformation is used, among other things, in [[geodesy]] to transform the coordinates of the point from one coordinate system into another. Using it, it becomes possible to convert regional [[surveying]] points into the [[WGS84]] locations used by [[Global Positioning System|GPS]]. | |||
In the process, the [[Gauss–Krüger coordinate system|Gauss–Krüger coordinate]], ''x'' and ''y'', plus the height, ''h'', are converted into 3D values in steps: | |||
# Calculation of the ellipsoidal width, length and height (''W'', ''L'', ''H'') | |||
# Calculation of ''X'', ''Y'' and ''Z'' relative to the [[reference ellipsoid]] of surveying | |||
# 7-parameter transformation (where ''X'', ''Y'' and ''Z'' change almost evenly, a few hundred metres at most, and the distances change a few mm per km). | |||
# Because of this, terrestrially measured positions can be compared with GPS data; these can then be brought into the surveying as new points — transformed in the opposite order. | |||
The third step consists of the application of a [[rotation matrix]], multiplication with the | |||
scale factor <math>\mu=1+s</math> (with a value near 1) and the addition of the three translations, <math>c_x</math>, <math>c_y</math>, <math>c_z</math>. | |||
The coordinates of a reference system B are derived from reference system A by the following formula:<ref>Datum Transformation Equations http://www.linz.govt.nz/geodetic/conversion-coordinates/geodetic-datum-conversion/datum-transformation-equations/index.aspx</ref> | |||
: <math>\begin{bmatrix} X \\ Y \\ Z \end{bmatrix}^B = \begin{bmatrix} c_x \\ c_y \\ c_z \end{bmatrix} + (1 + s\times10^{-6}) \cdot \begin{bmatrix} 1&-r_z&r_y \\ r_z&1&-r_x \\ -r_y & r_x & 1 \end{bmatrix} \cdot \begin{bmatrix} X \\ Y \\ Z \end{bmatrix}^A </math> | |||
or for each single parameter of the coordinate: | |||
: <math> | |||
\begin{matrix} | |||
X_B=c_x+(1+s\times10^{-6})\cdot (X_A-r_z\cdot Y_A+r_y\cdot Z_A)\\ | |||
Y_B=c_y+(1+s\times10^{-6})\cdot (r_z\cdot X_A+Y_A-r_x\cdot Z_A)\\ | |||
Z_B=c_z+(1+s\times10^{-6})\cdot (-r_y\cdot X_A+r_x\cdot Y_A+Z_A).\\ | |||
\end{matrix} | |||
</math> | |||
For the reverse transformation, each element is multiplied by -1. | |||
The seven parameters are determined for each region with three or more "identical points" of both systems. To bring them into agreement, the small inconsistencies (usually only a few cm) are adjusted using the method of [[least squares]] – that is, eliminated in a statistically plausible manner. | |||
== Standard parameters == | |||
'''Note that the rotation angles given in the table are in [[arcsecond|seconds]] and must be converted to [[radian]]s before use in the calculation.''' | |||
{| class="wikitable" | |||
!Region | |||
!Start datum | |||
!Target datum | |||
!<math>c_x</math> ([[Metre]]) | |||
!<math>c_y</math> (Metre) | |||
!<math>c_z</math> (Metre) | |||
!s ([[Parts per million|ppm]]) | |||
!r<sub>x</sub> ([[Arcsecond]]) | |||
!r<sub>y</sub> ([[Arcsecond]]) | |||
!r<sub>z</sub> ([[Arcsecond]]) | |||
|- | |||
|[[Slovenia ETRS89]] | |||
|D48 | |||
|[[D96 (datum)|D96]] | |||
| 409.545 | |||
| 72.164 | |||
| 486.872 | |||
| 17.919665 | |||
| −3.085957 | |||
| −5.469110 | |||
| 11.020289 | |||
|- | |||
|[[England]], [[Scotland]], [[Wales]] | |||
|[[WGS84]] | |||
|[[OSGB36]]<ref name=OSguide>A guide to coordinate systems in Great Britain v1.7 October 2007 D00659 Ordnance Survey</ref> | |||
| −446.448 | |||
| 125.157 | |||
| −542.06 | |||
| 20.4894 | |||
| −0.1502 | |||
| −0.247 | |||
| −0.8421 | |||
|- | |||
|[[Ireland]] | |||
|[[WGS84]] | |||
|Ireland 1965 | |||
| −482.53 | |||
| 130.596 | |||
| −564.557 | |||
| −8.15 | |||
| 1.042 | |||
| 0.214 | |||
| 0.631 | |||
|- | |||
|[[Germany]] | |||
|[[WGS84]] | |||
|[[Deutsches Hauptdreiecksnetz|DHDN]] | |||
| −591.28 | |||
| −81.35 | |||
| −396.39 | |||
| −9.82 | |||
| 1.4770 | |||
| −0.0736 | |||
| −1.4580 | |||
|- | |||
|[[Germany]] | |||
|[[WGS84]] | |||
|[[Bessel ellipsoid|Bessel 1841]] | |||
| −582 | |||
| −105 | |||
| −414 | |||
| −8.3 | |||
| −1.04 | |||
| −0.35 | |||
| 3.08 | |||
|- | |||
|[[Germany]] | |||
|[[WGS84]] | |||
|Krassovski 1940 | |||
| −24 | |||
| 123 | |||
| 94 | |||
| −1.1 | |||
| −0.02 | |||
| 0.26 | |||
| 0.13 | |||
|- | |||
| [[Austria]] (BEV) | |||
| [[WGS84]] | |||
| [[Datum Austria|MGI]] | |||
| −577.326 | |||
| −90.129 | |||
| −463.920 | |||
| −2.423 | |||
| 5.137 | |||
| 1.474 | |||
| 5.297 | |||
|- | |||
|[[USA]] | |||
|[[WGS84]] | |||
|Clarke 1866 | |||
| 8 | |||
| −160 | |||
| −176 | |||
| 0 | |||
| 0 | |||
| 0 | |||
| 0 | |||
|} | |||
These are standard parameter sets for the 7-parameter transformation (or data transformation) between two ellipsoids. For a transformation in the opposite direction, the signs of all the parameters must be changed. The translations ''c''<sub>''x''</sub>, ''c''<sub>''y''</sub>, ''c''<sub>''z''</sub> are sometimes described as ''t''<sub>''x''</sub>, ''t''<sub>''y''</sub>, ''t''<sub>''z''</sub>, or ''dx'', ''dy'', ''dz''. The rotations ''x'', ''y'', and ''z'' are sometimes also described as <math>\kappa</math>, <math>\phi</math> and <math>\omega</math>. In the United Kingdom the prime interest is the transformation between the OSGB36 datum used by the Ordnance survey for Grid References on its Landranger and Explorer maps to the WGS84 implementation used by GPS technology. The Gauss–Krüger coordinate system used in Germany normally refers to the [[Friedrich Bessel|Bessel]] ellipsoid. A further datum of interest was the [[ED50]] datum based on the [[Hayford ellipsoid]]. It is used by NATO for military maps and in future for most [[EU]] countries. | |||
The earth does not have a perfect ellipsoidal shape, but is described as a [[geoid]]. Instead, the geoid of the earth is described by many ellipsoids. Depending upon the actual location, the "locally best aligned ellipsoid" has been used for surveying and mapping purposes. The standard parameter set gives an accuracy of about 7 m for an OSGB36/WGS84 transformation. This is not precise enough for surveying, and the Ordnance Survey supplements these results by using a lookup table of further translations in order to reach 1 cm accuracy. | |||
==Restrictions== | |||
The Helmert transformation only uses one scale factor, so it is not suitable for: | |||
* The manipulation of measured drawings and [[Photography|photographs]] | |||
* The comparison of paper deformations while [[Image scanner|scanning]] old plans and maps. | |||
In these cases, use another [[affine transformation]]. | |||
==See also== | |||
* [[Galileo (satellite navigation)]] | |||
* [[Global Positioning System]] | |||
* [[Surveying]] | |||
==References== | |||
{{Reflist}} | |||
==External links== | |||
* http://www.w-volk.de/museum/mathex02.htm | |||
* http://www.webcitation.org/query?url=http://www.geocities.com/mapref/savpub/savpub-23.htm%23item40&date=2009-10-26+02:12:14 (Geometry for data exchange) | |||
* http://www.mapref.org/ | |||
* [http://trafostar.com TrafoStar] flexible 3D BestFit Transformations with: 3 Translations, 3 Rotations, 3 Scales, 3 Affine Parameters | |||
* [http://www.maths.dundee.ac.uk/gawatson/helmertrev.pdf Computing Helmert Transformations] | |||
{{Translation/Ref|de|Helmert-Transformation|oldid=40204821}} | |||
{{Use dmy dates|date=September 2010}} | |||
{{DEFAULTSORT:Helmert Transformation}} | |||
[[Category:Geodesy]] | |||
[[Category:Transformation (function)]] |
Revision as of 19:47, 3 February 2014
The Helmert transformation (named after Friedrich Robert Helmert, 1843–1917; also called a seven-parameter transformation) is a transformation method within a three-dimensional space. It is frequently used in geodesy to produce distortion-free transformations from one datum to another using:
where
- XT is the transformed vector
- X is the initial vector
The parameters are:
- — translation vector. Contains the three translations along the coordinate axes
- — scale factor, which is unitless, and as it is usually expressed in ppm, it must be divided by 1,000,000.
- — rotation matrix. Consists of three axes (small rotations around the coordinate axes) , , . The rotation matrix is an orthogonal matrix. The rotation is given in radians.
Thus, the Helmert transformation is a similarity transformation.
Calculating the parameters
If the transformation parameters are unknown, they can be calculated with reference points (that is, points whose coordinates are known before and after the transformation. Since a total of seven parameters (three translations, one scale, three rotations) have to be determined, at least two points and one coordinate of a third point (for example, the Z-coordinate) must be known. This gives a system of linear equations with seven equations and seven unknowns, which can be solved.
In practice, it is best to use more points. Through this correspondence, more accuracy is obtained, and a statistical assessment of the results becomes possible. In this case, the calculation is adjusted with the Gaussian least squares method.
A numerical value for the accuracy of the transformation parameters is obtained by calculating the values at the reference points, and weighting the results relative to the centroid of the points.
While the method is mathematically rigorous, it is entirely dependent on the rigour of the parameters that are used. In practice, these parameters are computed from the inclusion of at least three known points in the networks. However the accuracy of these will affect the following transformation parameters, as these points will contain observation errors. Therefore a "real-world" transformation will only be a best estimate and should contain a statistical measure of its quality.
It is not always necessary to use the seven parameter transformation, sometimes it is sufficient to use the five parameter transformation, composed of three translations, one rotation (about the Z-axis) and one change of scale.
Two-dimensional case
A special case is the two-dimensional Helmert transformation. Here, only four parameters are needed (two translations, one scaling, one rotation). These can be determined from two known points; if more points are available then checks can be made.
Application
The Helmert transformation is used, among other things, in geodesy to transform the coordinates of the point from one coordinate system into another. Using it, it becomes possible to convert regional surveying points into the WGS84 locations used by GPS.
In the process, the Gauss–Krüger coordinate, x and y, plus the height, h, are converted into 3D values in steps:
- Calculation of the ellipsoidal width, length and height (W, L, H)
- Calculation of X, Y and Z relative to the reference ellipsoid of surveying
- 7-parameter transformation (where X, Y and Z change almost evenly, a few hundred metres at most, and the distances change a few mm per km).
- Because of this, terrestrially measured positions can be compared with GPS data; these can then be brought into the surveying as new points — transformed in the opposite order.
The third step consists of the application of a rotation matrix, multiplication with the scale factor (with a value near 1) and the addition of the three translations, , , .
The coordinates of a reference system B are derived from reference system A by the following formula:[1]
or for each single parameter of the coordinate:
For the reverse transformation, each element is multiplied by -1.
The seven parameters are determined for each region with three or more "identical points" of both systems. To bring them into agreement, the small inconsistencies (usually only a few cm) are adjusted using the method of least squares – that is, eliminated in a statistically plausible manner.
Standard parameters
Note that the rotation angles given in the table are in seconds and must be converted to radians before use in the calculation.
Region | Start datum | Target datum | (Metre) | (Metre) | (Metre) | s (ppm) | rx (Arcsecond) | ry (Arcsecond) | rz (Arcsecond) |
---|---|---|---|---|---|---|---|---|---|
Slovenia ETRS89 | D48 | D96 | 409.545 | 72.164 | 486.872 | 17.919665 | −3.085957 | −5.469110 | 11.020289 |
England, Scotland, Wales | WGS84 | OSGB36[2] | −446.448 | 125.157 | −542.06 | 20.4894 | −0.1502 | −0.247 | −0.8421 |
Ireland | WGS84 | Ireland 1965 | −482.53 | 130.596 | −564.557 | −8.15 | 1.042 | 0.214 | 0.631 |
Germany | WGS84 | DHDN | −591.28 | −81.35 | −396.39 | −9.82 | 1.4770 | −0.0736 | −1.4580 |
Germany | WGS84 | Bessel 1841 | −582 | −105 | −414 | −8.3 | −1.04 | −0.35 | 3.08 |
Germany | WGS84 | Krassovski 1940 | −24 | 123 | 94 | −1.1 | −0.02 | 0.26 | 0.13 |
Austria (BEV) | WGS84 | MGI | −577.326 | −90.129 | −463.920 | −2.423 | 5.137 | 1.474 | 5.297 |
USA | WGS84 | Clarke 1866 | 8 | −160 | −176 | 0 | 0 | 0 | 0 |
These are standard parameter sets for the 7-parameter transformation (or data transformation) between two ellipsoids. For a transformation in the opposite direction, the signs of all the parameters must be changed. The translations cx, cy, cz are sometimes described as tx, ty, tz, or dx, dy, dz. The rotations x, y, and z are sometimes also described as , and . In the United Kingdom the prime interest is the transformation between the OSGB36 datum used by the Ordnance survey for Grid References on its Landranger and Explorer maps to the WGS84 implementation used by GPS technology. The Gauss–Krüger coordinate system used in Germany normally refers to the Bessel ellipsoid. A further datum of interest was the ED50 datum based on the Hayford ellipsoid. It is used by NATO for military maps and in future for most EU countries.
The earth does not have a perfect ellipsoidal shape, but is described as a geoid. Instead, the geoid of the earth is described by many ellipsoids. Depending upon the actual location, the "locally best aligned ellipsoid" has been used for surveying and mapping purposes. The standard parameter set gives an accuracy of about 7 m for an OSGB36/WGS84 transformation. This is not precise enough for surveying, and the Ordnance Survey supplements these results by using a lookup table of further translations in order to reach 1 cm accuracy.
Restrictions
The Helmert transformation only uses one scale factor, so it is not suitable for:
- The manipulation of measured drawings and photographs
- The comparison of paper deformations while scanning old plans and maps.
In these cases, use another affine transformation.
See also
References
43 year old Petroleum Engineer Harry from Deep River, usually spends time with hobbies and interests like renting movies, property developers in singapore new condominium and vehicle racing. Constantly enjoys going to destinations like Camino Real de Tierra Adentro.
External links
- http://www.w-volk.de/museum/mathex02.htm
- http://www.webcitation.org/query?url=http://www.geocities.com/mapref/savpub/savpub-23.htm%23item40&date=2009-10-26+02:12:14 (Geometry for data exchange)
- http://www.mapref.org/
- TrafoStar flexible 3D BestFit Transformations with: 3 Translations, 3 Rotations, 3 Scales, 3 Affine Parameters
- Computing Helmert Transformations
Template:Translation/Ref 30 year-old Entertainer or Range Artist Wesley from Drumheller, really loves vehicle, property developers properties for sale in singapore singapore and horse racing. Finds inspiration by traveling to Works of Antoni Gaudí.
- ↑ Datum Transformation Equations http://www.linz.govt.nz/geodetic/conversion-coordinates/geodetic-datum-conversion/datum-transformation-equations/index.aspx
- ↑ A guide to coordinate systems in Great Britain v1.7 October 2007 D00659 Ordnance Survey