Home    |    Concepts   |   API   |   Samples
Concepts > Geometry > Coordinate References > Coordinate Systems > Geographic Coordinate Systems
Geographic Transformation Methods

Moving your data between coordinate systems sometimes includes transforming between the geographic coordinate systems.

Because the geographic coordinate systems contain datums that are based on spheroids, a geographic transformation also changes the underlying spheroid. There are several methods, which have different levels of accuracy and ranges, for transforming between datums. The accuracy of a particular transformation can range from centimeters to meters depending on the method and the quality and number of control points available to define the transformation parameters.

A geographic transformation is always defined in a particular direction. When working with geographic transformations, if no mention is made of the direction, the command will handle the directionality automatically. For example, if converting data from WGS 1984 to NAD 1927, you can pick a transformation called NAD_1927_to_WGS_1984_3 and the software will apply it correctly.

A geographic transformation always converts geographic (longitude–latitude) coordinates. Some methods convert the geographic coordinates to geocentric (X,Y,Z) coordinates, transform the X,Y,Z coordinates, and convert the new values back to geographic coordinates.

These include the Geocentric Translation, Molodensky, Coordinate Frame, and Molodensky-Badekas methods.

Other methods, such as NADCON and NTv2 use a grid of differences and convert the longitude–latitude values directly.

Equation-based methods

Equation-based transformation methods can be classified into the following four method types. Usually the transformation parameters are defined as going from a local datum to WGS 1984 or another geocentric datum.

Three-parameter methods

The simplest datum transformation method is a geocentric, or three-parameter, transformation. The geocentric transformation models the differences between two datums in the X,Y,Z coordinate system. One datum is defined with its center at 0,0,0. The center of the other datum is defined at some distance (DX,DY,DZ) in meters away.

The three parameters are linear shifts and are always in meters.

Seven-parameter methods Top

A more complex and accurate datum transformation is possible by adding four more parameters to a geocentric transformation. The seven parameters are three linear shifts (DX,DY,DZ), three angular rotations around each axis (rx,ry,rz), and a scale factor.

The rotation values are given in decimal seconds, while the scale factor is in parts per million (ppm). The rotation values are defined in two different ways. It's possible to define the rotation angles as positive either clockwise or counterclockwise as you look toward the origin of the X,Y,Z systems.

The United States, Australia, New Zealand, and a few other countries define the equations such that the rotation values are positive counterclockwise. This method is called the Coordinate Frame Rotation transformation. Europe uses a different convention called the Position Vector transformation. Both methods are sometimes referred to as the Bursa–Wolf method. In the Projection Engine, the Coordinate Frame and Bursa–Wolf methods are the same. Both Coordinate Frame and Position Vector methods are supported, and it is easy to convert transformation values from one method to the other simply by changing the signs of the three rotation values. For example, the parameters to convert from the WGS 1972 datum to the WGS 1984 datum with the Coordinate Frame method are (in the order DX,DY,DZ,rx,ry,rz,s):

(0.0, 0.0, 4.5, 0.0, 0.0, -0.554, 0.227)

To use the same parameters with the Position Vector method, change the sign of the rotation so the new parameters are:

(0.0, 0.0, 4.5, 0.0, 0.0, +0.554, 0.227)

It's impossible to tell from the parameters alone which convention is being used. If you use the wrong method, your results can return inaccurate coordinates. The only way to determine how the parameters are defined is by checking a control point whose coordinates are known in the two systems.

The Molodensky–Badekas method is a variation of the seven-parameter methods. It has an additional three parameters that define the XYZ origin of rotation. Sometimes this point is known as the origin of the datum, or geographic coordinate system. Given the XYZ origin of rotation point, it is possible to calculate an equivalent Coordinate Frame transformation. The DX, DY, and DZ values will change but the rotation and scale values will remain the same.

Molodensky method Top

The Molodensky method converts directly between two geographic coordinate systems without actually converting to an X,Y,Z system. The Molodensky method requires three shifts (DX,DY,DZ) and the differences between the semimajor axes (Da) and the flattenings (Df) of the two spheroids. The Projection Engine automatically calculates the spheroid differences according to the datums involved.

Abridged Molodensky method

The Abridged Molodensky method is a simplified version of the Molodensky method.

Grid-based methods

Grid-based transformation methods include the following:

NADCON and HARN methods

The United States uses a grid-based method to convert between geographic coordinate systems. Grid-based methods allow you to model the differences between the systems and are potentially the most accurate method. The area of interest is divided into cells. The National Geodetic Survey (NGS) publishes grids to convert between NAD 1927 and other older geographic coordinate systems and NAD 1983. These transformations are grouped into the NADCON method. The main NADCON grid, CONUS, converts the contiguous 48 states. The other NADCON grids convert older geographic coordinate systems to NAD 1983 for:

  • Alaska
  • Hawaiian islands
  • Puerto Rico and Virgin Islands
  • St. George, St. Lawrence, and St. Paul Islands in Alaska

The accuracy is approximately 0.15 meters for the contiguous states, 0.50 for Alaska and its islands, 0.20 for Hawaii, and 0.05 for Puerto Rico and the Virgin Islands. Accuracies can vary depending on how good the geodetic data in the area was when the grids were computed (NADCON, 1999).

The Hawaiian islands were never on NAD 1927. They were mapped using several datums that are collectively known as the Old Hawaiian datums.

New surveying and satellite measuring techniques have allowed NGS and the states to update the geodetic control point networks. As each state is finished, the NGS publishes a grid that converts between NAD 1983 and the more accurate control point coordinates. Originally, this effort was called the High Precision Geodetic Network (HPGN). It is now called the High Accuracy Reference Network (HARN). Four territories and 46 states have published HARN grids as of January 2004. HARN transformations have an accuracy approximately 0.05 meters (NADCON, 2000).

The difference values in decimal seconds are stored in two files: one for longitude and the other for latitude. A bilinear interpolation is used to calculate the exact difference between the two geographic coordinate systems at a point. The grids are binary files, but a program, NADGRD, from the NGS, allows you to convert the grids to American Standard Code for Information Interchange (ASCII) format. Shown at the bottom of the page is the header and first row of the CSHPGN.LOA file. This is the longitude grid for Southern California. The format of the first row of numbers is, in order, the number of columns, number of rows, number of z–values (always one), minimum longitude, cell size, minimum latitude, cell size, and not used.

The next 37 values in this case are the longitude shifts from -122° to -113° at 32° N in 0.25°, or 15 minute, intervals in longitude.

NADCON EXTRACTED REGION NADGRD

37 21 1 -122.00000 .25 32.00000 .25 .00000

.007383 .004806 .002222 -.000347 -.002868
-.005296 -.007570 -.009609 -.011305 -.012517
-.013093 -.012901 -.011867 -.009986 -.007359
-.004301 -.001389 .001164 .003282 .004814
.005503 .005361 .004420 .002580 .000053
-.002869 -.006091 -.009842 -.014240 -.019217
-.025104 -.035027 -.050254 -.072636 -.087238
-.099279 -.110968
Top
Send your comments to: Site Administrator | Terms of Use | PRIVACY | Copyright © ESRI.