photograph: Desert Tortoise The Advanced Resource Technology Group:
The Knowledge Base



Last updated
6/19/02

Back

By: Bill Huber
Published in: Directions Magazine
(Jul 18, 2000)

Reprojecting geographic features

Two things must happen when you draw a map: features in the real world must be "georeferenced" to a spheroid and the spheroid must be projected onto the paper.

The spheroid models the shape of the earth's surface. It is an idealization that does not account for local changes in topography.

Georeferencing assigns locations (in three dimensions!) to points on a spheroid.

Projecting is an operation that mathematically distorts and shrinks a portion of the spheroid onto flat paper. Projecting can be undone ("inverted"). "Unprojection" expands a feature on a map and plasters it back onto the spheroid. It, too, is a mathematical operation.

Georeferencing is done with a datum. A datum is usually given by a starting point and direction: it specifies where a clearly identifiable point on earth (the base point) should appear on the spheroid and it shows where a base direction, such as north, points on the spheroid at the base point. The base point and direction allow surveyors to determine the distance and angle of any other point on the earth. Moving in the corresponding direction on the spheroid for the same distance determines where the new point should go on the spheroid.

Spheroids have coordinates. They are latitude and longitude. (Geodetic) latitude is the angle made by a vertical line to the horizontal. It is not necessarily the same angle made by "straight up," because the latter is distorted by gravitational variation over the earth. It is not necessarily the angle made by a line to the center of the earth, because most spheroids have an elliptical cross-section, not a circular one.

Therefore, georeferencing endows points near the earth with latitude, longitude, and height coordinates.

Change of datum

When more than one datum is in play, the coordinates depend on the datum.

The figure above shows this in two dimensions, suppressing the third for clarity. The red ellipse represents one spheroid and the blue ellipse represents another. Notice they have different shapes and even different centers. (x,y) denotes a point near the earth. Theta is its geodetic latitude and h is its geodetic height with respect to the red spheroid. Theta' and h' are its coordinates with respect to the blue spheroid. Evidently, at this point the height and the latitude on the blue spheroid are each smaller than their counterparts on the red spheroid.

(The values a, b, a', and b' merely show the different radii for the spheroids.)

In practice, these differences in coordinates are large enough to be important in many applications. As the U.S. military says, "Ignoring the fine print in the margin of a map could get you killed!"

How to relate two maps

Two maps may use two different datums and two different projections. The process of relating features on one map to features on another is involved, but straightforward, as the next figure shows.

To shift the red feature from the yellow map ("source") to the blue map ("destination"), or reproject it, follow four steps in this order:

  1. Unproject the source map. This produces the red feature on the red spheroid.
  2. Determine the real-world location of the feature. This is the gray patch at the top of the figure.
  3. Reference the real feature to the destination datum. This is the blue feature on the blue spheroid.
  4. Project the newly referenced feature.
Often steps 2 and 3 are combined into one operation, a "datum shift." Frequently, approximate formulas are used. They depend only on the feature and the sizes, shapes, and positions of the spheroids.

The datum shift guarantees that the feature on the destination map correctly represents the same real-world object depicted by the feature on the source map.

The wrong way to do it

We have to mention one wrong way, because limitations in mapping software tempt people to do this. What happens when we ignore the datum shift?

An apparent shift in datum occurs anyway. This comes about because software (or, in earlier times, the draftsman) uses spheroid coordinates to compute the projection. The red feature on the red spheroid is determined by its coordinates. That same set of coordinates, when drawn on the blue spheroid, produces the blue feature shown. When projected, its map is the green feature in the destination map.

However, this process usually produces the wrong result.

We say "usually," because in some cases the net error is smaller than imprecision in the maps themselves. But in other cases the error can be hundreds of meters, even kilometers. That is a large distortion on detailed ("large scale") maps.

North America is a special case

Two datums for regions in North America are in common use: NAD 27 and NAD 83. They are based on different spheroids. Furthermore, fairly large errors accumulated in the surveyed locations using the earlier NAD 27 datum. These errors are difficult to describe mathematically. Instead, the U.S. National Geodetic Survey (NGS) offers the NADCON computer program, with Fortran source code, to convert coordinates between these two datums. NADCON effectively performs both steps 2 and 3 in the reprojection algorithm. Therefore you do not need to worry about the spheroids themselves when using NADCON. You do need to have the appropriate data files describing the changes. They are also available from the NADCON web site.


More about this author


Bill Huber
Company Quantitative Decisions
Email whuber@quantdec.com
URL http://www.quantdec.com/
Dr. Huber, principal in the GIS consulting firm Quantitative Decisions in Merion Station, Pennsylvania, has been among the most prolific contributors to ESRI's ArcView-L for the past four years. He is the founder of an independent ArcView discussion list, and is also on the faculty at Penn State University. Dr. Huber is Acting Editor for Directions Magazine.

Articles
GIS & Steganography - Part 3: Vector Steganography Apr 18, 2002
Why only 'Women in GIS'? Aug 22, 2001
Map Data Controversy in the U.S.: Analysis and Opinion Aug 09, 2001
Analysis and Forecast of USGS GeoData Demands on the Web Aug 09, 2001
Map Data Controversy in the U.S. Jul 31, 2001
ESRI UC 2001: On the Scene at ESRI UC Jul 10, 2001
Reducing the Info:Noise Ratio Apr 25, 2001
What do GIS consultants charge? Mar 13, 2001
Reprojecting Grids Jul 18, 2000
ESRI UC 2000: Closing Session Jun 30, 2000
ESRI UC 2000: ArcView's Transformation Jun 29, 2000
ESRI UC 2000: Technical Workshops Begin Jun 27, 2000
ESRI UC 2000: ESRI Launches The Geography Network Jun 27, 2000
ESRI UC 2000: On the Scene at ESRI UC Jun 26, 2000

Features
A Review of IDRISI32 Oct 25, 2000
The Future of Arcview; Part 2 Aug 08, 2000
The Future of ArcView; Part 1 Aug 01, 2000
Adding dimensions to GIS with VRML Jun 05, 2000
ArcView Class Hunt Results Apr 11, 2000
Convolution: Part 3 of 3 Nov 01, 1999
Sidebar to "Convolution" Oct 25, 1999
Convolution: Part 2 of 3 Oct 25, 1999
Convolution Oct 18, 1999

Back