Chapter 1 outlined several of the distinguishing properties of geographic data. One is that geographic data are necessarily generalized, and that generalization tends to vary with scale. A second distinguishing property is that the Earth's complex, nearly-spherical shape complicates efforts to specify exact positions on Earth's surface. This chapter explores implications of these properties by illuminating concepts of scale, Earth geometry, coordinate systems, the "horizontal datums" that define the relationship between coordinate systems and the Earth's shape, and the various methods for transforming coordinate data between 3D and 2D grids, and from one datum to another.

Compared to Chapter 1, Chapter 2 may seem long, technical, and abstract, particularly to those for whom these concepts are new.

Students who successfully complete Chapter 2 should be able to:

- Demonstrate your ability to specify geospatial locations using geographic coordinates;
- Convert geographic coordinates between two different formats;
- Explain the concept of a horizontal datum;
- Calculate the change in a coordinate location due to a change from one horizontal datum to another;
- Estimate the magnitude of "datum shift" associated with the adjustment from NAD 27 to NAD 83;
- Recognize the kind of transformation that is appropriate to georegister two or more data sets;
- Describe the characteristics of the UTM coordinate system, including its basis in the Transverse Mercator map projection;
- Plot UTM coordinates on a map;
- Describe the characteristics of the SPC system, including map projection on which it is based;
- Convert geographic coordinates to SPC coordinates;
- Interpret distortion diagrams to identify geometric properties of the sphere that are preserved by a particular projection; and
- Classify projected graticules by projection family.

Take a minute to complete any of the Try This activities that you encounter throughout the chapter. These are fun, thought-provoking exercises to help you better understand the ideas presented in the chapter.

This textbook is used as a resource in Penn State's Online Geospatial Education online degree and certificate programs. If this topic is interesting to you and you want to learn more about online GIS and GEOINT education at Penn State, check out our Geospatial Education Program Office.

You hear the word "scale" often when you work around people who produce or use geographic information. If you listen closely, you'll notice that the term has several different meanings, depending on the context in which it is used. You'll hear talk about the scales of geographic phenomena and about the scales at which phenomena are represented on maps and aerial imagery. You may even hear the word used as a verb, as in "scaling a map" or "downscaling." The goal of this section is to help you learn to tell these different meanings apart, and to be able to use concepts of scale to help make sense of geographic data.

Specifically, in this part of Chapter 2 you will learn to:

- calculate map scale using representative fractions;
- describe the general relationship between map scale, detail, and accuracy.

Often "scale" is used as a synonym for "scope" or "extent." For example, the title of an international research project called The Large Scale Biosphere-Atmosphere Experiment in Amazonia (1999) uses the term "large scale" to describe a comprehensive study of environmental systems operating across a large region. This usage is common not only among environmental scientists and activists, but also among economists, politicians, and the press. Those of us who specialize in geographic information usually use the word "scale" differently, however.

When people who work with maps and aerial images use the word "scale," they usually are talking about the sizes of things that appear on a map or air photo, relative to the actual sizes of those things on the ground.

**Map scale is the proportion between a distance on a map and a corresponding distance on the ground:**

(D_{m} / D_{g}).

By convention, the proportion is expressed as a "representative fraction" in which map distance (D_{m}) is reduced to 1. The proportion, or ratio, is also typically expressed in the form 1 : D_{g} rather than 1 / D_{g}.

The representative fraction 1:100,000, for example, means that a section of road that measures 1 unit in length on a map stands for a section of road on the ground that is 100,000 units long.

If we were to change the scale of the map such that the length of the section of road on the map was reduced to, say, 0.1 units in length, we would have created a **smaller-scale** map whose representative fraction is 0.1:100,000, or 1:1,000,000. When we talk about large- and small-scale maps and geographic data, then, we are talking about the relative sizes and levels of detail of the features represented in the data. In general, the larger the map scale, the more detail is shown. This tendency is illustrated below in Figure 2.5.1.

Adapted from Thompson, 1988.

One of the defining characteristics of topographic maps is that scale is consistent across each map and within each map series. This isn't true for aerial imagery, however, except for images that have been **orthorectified**. As discussed in Chapter 6, large scale maps are typically derived from aerial imagery. One of the challenges associated with using air photos as sources of map data is that the scale of an aerial image varies from place to place as a function of the elevation of the terrain shown in the scene. Assuming that the aircraft carrying the camera maintains a constant flying height (which pilots of such aircraft try very hard to do), the distance between the camera and the ground varies along each flight path. This causes air photo scale to be larger where the terrain is higher and smaller where the terrain is lower. An "orthorectified" image is one in which variations in scale caused by variations in terrain elevation (among other effects) have been removed.

You can calculate the average scale of an** unrectified** air photo by solving the equation **S _{p} = f / (H-h_{avg})**, where f is the focal length of the camera, H is the flying height of the aircraft above mean sea level, and h

This textbook is used as a resource in Penn State's Online Geospatial Education online degree and certificate programs. If this topic is interesting to you and you want to learn more about online GIS and GEOINT education at Penn State, check out our Geospatial Education Program Office.

Another way to express map scale is with a graphic (or "bar") scale. Unlike representative fractions, graphic scales remain true when maps are shrunk or magnified.

If they include a scale at all, most maps include a bar scale like the one shown above left (Figure 2.6.1). Some also express map scale as a representative fraction. Either way, the implication is that scale is uniform across the map. In fact, except for maps that show only very small areas, scale varies across every map. As you probably know, this follows from the fact that positions on the nearly-spherical Earth must be transformed to positions on two-dimensional sheets of paper. Systematic transformations of this kind are called **map projections**. As we will discuss in greater depth later in this chapter, all map projections are accompanied by deformation of features in some or all areas of the map. This deformation causes map scale to vary across the map. Representative fractions may, therefore, specify map scale along a line at which deformation is minimal (**nominal scale**). Bar scales denote only the nominal or average map scale. Variable scales, like the one illustrated above right, show how scale varies, in this case by latitude, due to deformation caused by map projection.

One of the special characteristics of geographic data is that phenomena shown on maps tend to be represented differently at different scales. Typically, as scale decreases, so too does the number of different features and the detail with which they are represented. Not only printed maps, but also digital geographic data sets that cover extensive areas, tend to be more **generalized** than datasets that cover limited areas.

Accuracy also tends to vary in proportion with map scale. The United States Geological Survey, for example, guarantees that the mapped positions of 90 percent of well-defined points shown on its topographic map series at scales smaller than 1:20,000 will be within 0.02 inches of their actual positions on the map (see the National Geospatial Program Standards and Specifications). Notice that this "National Map Accuracy Standard" is scale-dependent. The allowable error of well-defined points (such as control points, road intersections, and such) on 1:250,000 scale topographic maps is thus 1 / 250,000 = 0.02 inches / D_{g} or D_{g} = 0.02 inches x 250,000 = 5,000 inches or 416.67 feet. Neither small-scale maps nor the digital data derived from them are reliable sources of detailed geographic information.

Credit: Chem-Nuclear Systems, Inc., 1994

Sometimes the detail lost on small-scale maps causes serious problems. For example, a contractor hired to use GIS to find a suitable site for a low-level radioactive waste storage facility in Pennsylvania presented a series of 1:1,500,000 scale maps at public hearings around the state in the early 1990s. The scale was chosen so that disqualified areas of the entire state could be printed on a single 11 x 17-inch page. A report accompanying the map included the disclaimer that "it is possible that small areas of sufficient size for the LLRW disposal facility site may exist within regions that appear disqualified on the [map]. The detailed information for these small areas is retained within the GIS even though they are not visually illustrated..." (Chem-Nuclear Systems, Inc. 1993, p. 20). Unfortunately for the contractor, alert citizens recognized the shortcomings of the small-scale map, and newspapers published reports accusing the out-of-state company of providing inaccurate documents. Subsequent maps were produced at a scale large enough to discern 500-acre suitable areas.

The term "scale" is sometimes used as a verb. To scale a map is to reproduce it at a different size. For instance, if you photographically reduce a 1:100,000-scale map to 50 percent of its original width and height, the result would be one-quarter the area of the original. Obviously, the map scale of the reduction would be smaller too: 1/2 x 1/100,000 = 1/200,000.

Because of the inaccuracies inherent in all geographic data, particularly in small scale maps, scrupulous geographic information specialists **avoid enlarging source maps**. To do so is to exaggerate generalizations and errors. The original map used to illustrate areas in Pennsylvania disqualified from consideration for low-level radioactive waste storage shown on an earlier page, for instance, was printed with the statement "Because of map scale and printing considerations, it is not appropriate to enlarge or otherwise enhance the features on this map."

This textbook is used as a resource in Penn State's Online Geospatial Education online degree and certificate programs. If this topic is interesting to you and you want to learn more about online GIS and GEOINT education at Penn State, check out our Geospatial Education Program Office.

The word "scale" can also be used as a synonym for a ruler--a measurement scale. Because data consist of symbols that represent measurements of phenomena, it's important to understand the reference systems used to take the measurements in the first place. In this section, we'll consider a measurement scale known as the **geographic coordinate system** that is used to specify positions on the Earth's roughly spherical surface. In other sections, we'll encounter two-dimensional (plane) coordinate systems, as well as the measurement scales used to specify attribute data.

In this section of Chapter 2 you will:

- demonstrate your ability to specify geospatial locations using geographic coordinates;
- convert geographic coordinates between two different formats.

As you probably know, locations on the Earth's surface are measured and represented in terms of **coordinates.** A coordinate is a set of two or more numbers that specifies the position of a point, line, or other geometric figure in relation to some reference system. The simplest system of this kind is a Cartesian coordinate system (named for the 17th century mathematician and philosopher René Descartes). A Cartesian coordinate system is simply a grid formed by juxtaposing two measurement scales, one horizontal (x) and one vertical (y). The point at which both x and y equal zero is called the **origin** of the coordinate system. In Figure 2.10.1, above, the origin (0,0) is located at the center of the grid. All other positions are specified relative to the origin. The coordinate of upper right-hand corner of the grid is (6,3). The lower left-hand corner is (-6,-3). If this is not clear, please ask for clarification!

Cartesian and other two-dimensional (plane) coordinate systems are handy due to their simplicity. For obvious reasons, they are not perfectly suited to specifying geospatial positions, however. The **geographic coordinate system** is designed specifically to define positions on the Earth's roughly-spherical surface. Instead of the two linear measurement scales, x and y, the geographic coordinate systems juxtaposes two curved measurement scales. The east-west scale, called **longitude** (conventionally designated by the Greek symbol *lambda*), ranges from +180° to -180°. Because the Earth is round, +180° (or 180° E) and -180° (or 180° W) are the same grid line. That grid line is roughly the International Date Line, which has diversions that pass around some territories and island groups. Opposite the International Date Line is the **prime meridian**, the line of longitude defined by international treaty as 0°. The north-south scale, called **latitude** (designated by the Greek symbol *phi*), ranges from +90° (or 90° N) at the North pole to -90° (or 90° S) at the South pole. We'll take a closer look at the geographic coordinate system next.

Credit: David Dibase

**Longitude** specifies positions east and west as the angle between the **prime meridian** and a second **meridian** that intersects the point of interest. Longitude ranges from +180 (or 180° E) to -180° (or 180° W). 180° East and West longitude together form the International Date Line.

**Latitude** specifies positions north and south in terms of the angle subtended at the center of the Earth between two imaginary lines, one that intersects the **equator **and another that intersects the point of interest. Latitude ranges from +90° (or 90° N) at the North pole to -90° (or 90° S) at the South pole. A line of latitude is also known as a **parallel**.

At higher latitudes, the length of parallels decreases to zero at 90° North and South. Lines of longitude are not parallel but converge toward the poles. Thus, while a degree of longitude at the equator is equal to a distance of about 111 kilometers, that distance decreases to zero at the poles.

**Geographic Coordinate System Practice Application**

Nearly everyone learned latitude and longitude as a kid. But how well do you understand the geographic coordinate system, really? My experience is that while most everyone has heard of latitude and longitude, only about half can point to the location on a map that is specified by a pair of geographic coordinates. The Flash application linked below lets you test your knowledge. The application asks you to click locations on a globe as specified by randomly generated geographic coordinates.

You will notice that the application lets you choose between "easy problems" and "hard problems." Easy problems are those in which latitude and longitude coordinates are specified in 30° increments. Since the resolution of the graticule (the geographic coordinate system grid) used in the application is also 30°, the solution to every "easy" problem occurs at the intersection of a parallel and a meridian. The "easy" problems are good warm-ups.

"Hard" problems specify coordinates in 1° increments. You have to interpolate positions between grid lines. **You can consider yourself to have a good working knowledge of the geographic coordinate system if you can solve at least six "hard" problems consecutively and on the first click.**

Click here to download and launch the Geographic Coordinate System practice application (5.7 Mb). (If the globe doesn't appear after the Flash application has loaded, right-click and select "Play" from the pop-up menu.)

Note: You will need to have the Adobe Flash player installed in order to complete this exercise. If you do not already have it, you can download the Flash Player for free at the Adobe website.

This textbook is used as a resource in Penn State's Online Geospatial Education online degree and certificate programs. If this topic is interesting to you and you want to learn more about online GIS and GEOINT education at Penn State, check out ourGeospatial Education Program Office.

Geographic coordinates may be expressed in decimal degrees, or in degrees, minutes, and seconds. Sometimes, you need to convert from one form to another. Steve Kiouttis (personal communication, Spring 2002), manager of the Pennsylvania Urban Search and Rescue Program, described one such situation on the course Bulletin Board: "I happened to be in the state Emergency Operations Center in Harrisburg on Wednesday evening when a call came in from the Air Force Rescue Coordination Center in Dover, DE. They had an emergency locator transmitter (ELT) activation and requested the PA Civil Air Patrol to investigate. The coordinates given to the watch officer were 39 52.5 n and -75 15.5 w. This was plotted incorrectly (treated as if the coordinates were in decimal degrees 39.525n and -75.155 w) and the location appeared to be near Vineland, New Jersey. I realized that it should have been interpreted as 39 degrees 52 minutes and 5 seconds n and -75 degrees and 15 minutes and 5 seconds w) and made the conversion (as we were taught in Chapter 2) and came up with a location on the grounds of Philadelphia International Airport, which is where the locator was found, in a parked airliner."

Here's how it works:

To convert -89.40062 from decimal degrees to degrees, minutes, seconds:

- Subtract the number of whole degrees (89°) from the total (89.40062°). (The minus sign is used in the decimal degree format only to indicate that the value is a west longitude or a south latitude.)
- Multiply the remainder by 60 minutes (.40062 x 60 = 24.0372).
- Subtract the number of whole minutes (24') from the product.
- Multiply the remainder by 60 seconds (.0372 x 60 = 2.232).
- The result, expressed in the correct number of significant figures, is 89° 24' 2.2" W or S.

To convert 43° 4' 31" from degrees, minutes, seconds to decimal degrees:

DD = Degrees + (Minutes/60) + (Seconds/3600)

- Divide the number of seconds by 60 (31 ÷ 60 = 0.5166).
- Add the quotient of step (1) to the whole number of minutes (4 + 0.5166).
- Divide the result of step (2) by 60 (4.5166 ÷ 60 = 0.0753).
- Add the quotient of step (3) to the number of whole number degrees (43 + 0.0753).
- The result, expressed in the correct number of significant figures, is 43.075°

Geographic data represent the locations and attributes of things on the Earth's surface. Locations are measured and encoded in terms of geographic coordinates (i.e., latitude and longitude) or plane coordinates (e.g., UTM). To measure and specify coordinates accurately, one first must define the geometry of the surface itself. To see what I mean, imagine a soccer ball. If you or your kids play soccer you can probably conjure up a vision of a round mosaic of 20 hexagonal (six sided) and 12 pentagonal (five sided) panels (soccer balls come in many different designs, but the 32-panel ball is used in most professional matches. Visit soccerballworld.com for more than you ever wanted to know about soccer balls). Now focus on one point at an intersection of three panels. You could use spherical (e.g., geographic) coordinates to specify the position of that point. But if you deflate the ball, the position of the point in space changes, and so must its coordinates. **The absolute** (though not the relative) **position of a point on a surface, then, depends upon the shape of the surface**.

Every position is determined in relation to at least one other position. Coordinates, for example, are defined relative to the origin of the coordinate system grid. A land surveyor measures the "corners" of a property boundary relative to a previously-surveyed control point. Surveyors and engineers measure elevations at construction sites and elsewhere. Elevations are expressed in relation to a **vertical datum**, a reference surface such as mean sea level. As you probably know, there is also such a thing as a **horizontal datum**, although this is harder to explain and to visualize than the vertical case. **Horizontal datums define the geometric relationship between a coordinate system grid and the Earth's surface.** Because the Earth's shape is complex, the relationship is too. The goal of this section is to explain the relationship.

Specifically, in this section of Chapter 2 you will learn to:

- explain the concept of a horizontal datum;
- calculate the change in a coordinate location due to a change from one horizontal datum to another;
- estimate the magnitude of "datum shift" associated with the adjustment from NAD 27 to NAD 83.

Credit: Adapted from Smith, 1988

The accuracy of coordinates that specify geographic locations depends upon how the coordinate system grid is aligned with the Earth's surface. Unfortunately for those who need accurate geographic data, defining the shape of the Earth's surface is a non-trivial problem. So complex is the problem that an entire profession, called **geodesy**, has arisen to deal with it.

Geodesists define the Earth's surface as **a surface that closely approximates global mean sea level, but across which gravity is everywhere equal.** They refer to this shape as the **geoid**. Geoids are lumpy because gravity varies from place to place in response to local differences in terrain and variations in the density of materials in the Earth's interior. Geoids are also a little squat. Sea level gravity at the poles is greater than sea level gravity at the equator, a consequence of Earth's "oblate" shape as well as the centrifugal force associated with its rotation.

Geodesists at the U.S. National Geodetic Survey describe the geoid as an "equipotential surface" because the potential energy associated with the Earth's gravitational pull is equivalent everywhere on the surface. Like fitting a trend line through a cluster of data points, the geoid is a three-dimensional statistical surface that fits as closely as possible gravity measurements taken at millions of locations around the world. As additional and more accurate gravity measurements become available, geodesists revise the geoid periodically. Some geoid models are solved only for limited areas; GEOID03, for instance, is calculated only for the continental U.S.

Recall that horizontal datums define how coordinate system grids align with the Earth's surface. Long before geodesists calculated geoids, surveyors used much simpler surrogates called **ellipsoids** to model the shape of the Earth.

Credit: Adapted from Smith, 1988

An ellipsoid is a three-dimensional geometric figure that resembles a sphere, but whose equatorial axis (*a* in Figure 2.15.1, above) is slightly longer than its polar axis (*b*). The equatorial axis of the World Geodetic System of 1984, for instance, is approximately 22 kilometers longer than the polar axis, a proportion that closely resembles the oblate spheroid that is planet Earth. **Ellipsoids are commonly used as surrogates for geoids so as to simplify the mathematics involved in relating a coordinate system grid with a model of the Earth's shape.** Ellipsoids are good, but not perfect, approximations of geoids. The map in Figure 2.15.2, below shows differences in elevation between a geoid model called GEOID96 and the WGS84 ellipsoid. The surface of GEOID96 rises up to 75 meters above the WGS84 ellipsoid over New Guinea (where the map is colored red). In the Indian Ocean (where the map is colored purple), the surface of GEOID96 falls about 104 meters below the ellipsoid surface.

Credit: National Geodetic Survey, 1997

Many ellipsoids are in use around the world. (Peter Dana has published a list at colorado.edu) Local ellipsoids minimize differences between the geoid and the ellipsoid for individual countries or continents. The Clarke 1866 ellipsoid, for example, minimizes deviations in North America. The **North American Datum of 1927** (NAD 27) associates the geographic coordinate grid with the Clarke 1866 ellipsoid. NAD 27 involved an adjustment of the latitude and longitude coordinates of some 25,000 geodetic control point locations across the U.S. The nationwide adjustment commenced from an initial control point at Meades Ranch, Kansas, and was meant to reconcile discrepancies among the many local and regional control surveys that preceded it.

**The North American Datum of 1983** (NAD 83) involved another nationwide adjustment, necessitated in part by the adoption of a new ellipsoid, called **GRS 80**. Unlike Clarke 1866, GRS 80 is a global ellipsoid centered upon the Earth's center of mass. GRS 80 is essentially equivalent to WGS 84, the global ellipsoid upon which the Global Positioning System is based. **NAD 27 and NAD 83 both align coordinate system grids with ellipsoids. They differ simply in that they refer to different ellipsoids.** Because Clarke 1866 and GRS 80 differ slightly in shape as well as in the positions of their center points, the adjustment from NAD 27 to NAD 83 involved a shift in the geographic coordinate grid. Because a variety of datums remain in use, geospatial professionals need to understand this shift, as well as how to transform data between horizontal datums.

Credit: National Geodetic Survey, 2004

Geoids, ellipsoids, and even coordinate systems are all abstractions. The fact that "horizontal datum" refers to a relationship between an ellipsoid and a coordinate system, two abstractions, may explain why the concept is so frequently misunderstood. Datums do have physical manifestations, however.

Shown above (Figure 2.16.1) is one of the approximately two million horizontal and vertical control points that have been established in the U.S. Although control point markers are fixed, the coordinates that specify their locations are liable to change. The U.S. National Geodetic Survey maintains a database of the coordinate specifications of these control points, including historical locations as well as more recent adjustments. One occasion for adjusting control point coordinates is when new horizontal datums are adopted.** Since every coordinate system grid is aligned with an ellipsoid that approximates the Earth's shape, coordinate grids necessarily shift when one ellipsoid is replaced by another.** When coordinate system grids shift, the coordinates associated with fixed control points need to be adjusted. How we account for the Earth's shape makes a difference in how we specify locations.

Here's a chance to calculate how much the coordinates of a control point change in response to an adjustment from North American Datum of 1927 (based on the Clarke 1866 ellipsoid) to the North American Datum of 1983 (based upon the GRS 80 ellipsoid).

**Find the geographic coordinates of a populated place**- Start at the USGS' Geographic Names Information System at the U.S. Board on Geographic Names.
- Follow the links labeled
**Domestic Names,**then**Search**to search place names included in the Geographic Names Information System. - At the Query Form, enter the name of your home town (or other named geographic feature) in the
**Feature Name**field, as well as your home**State**. Choose**Populated Place**(or other, as appropriate) for the**Feature Class**.- If your home is somewhere other than the U.S., enter a place name of interest or fantasy destination (e.g., "Las Vegas" ;-) ).

- Click
**Send Query**. - The result should include latitude and longitude coordinates of a
**centroid**that represents where the name your town (or other feature) would appear on a map. You'll need those coordinates for the next step.

**Find the geographic coordinates of a nearby horizontal control point**- Visit the U.S. National Geodetic Survey home page.
- Follow the link labeled
**Survey Mark D****atasheets**. - At the NGS Datasheet page, follow the link labeled
**Datasheets**.- You may wish to begin with the "Info Link" labeled "Tell me more about datasheets."

- At the NGS Datasheet Retrieval page, follow the link labeled
**Radial Search**. (You're welcome to experiment with another retrieval method if you wish.) - At the NGS Datasheet Point Radius form:
- Enter the latitude and longitude coordinates you looked up in step #1. Pay attention to the input format.
- Specify a
**Search Radius**. - Select
**Any Horz. and/or Vert. Control**from the Data Type Desired scrolling field. - Select
**Any Stability**from the Stability Desired scrolling field. - Click
**Submit**.

- The result should be a
**Station List Results**form that looks like the contents of the window pictured below. These are the results of my search on the centroid coordinates for State College PA. Note that I have highlighted the station that is both nearest to the coordinates I entered and a first-order control point (see the "1" under the column labeled "H"?)Credit: noaa.gov - Select the station nearest to the coordinates you specified that is also the highest-order horizontal control point.
- Click
**Get Datasheets**. The system should respond with a station datasheet like this example. - In the example linked above, the CURRENT SURVEY CONTROL of the station point is listed as
**NAD 83(1992) 40 48 13.83840(N) 077 51 44.25410(W) ADJUSTED.**These are the geographic coordinates of the control point relative to the NAD 83 horizontal datum. In the next step, we'll see how much the control point "moved" as a result of the adjustment of those coordinates from the earlier NAD 27 datum. (The geographic coordinates of the control point are specified to 100,000th of a second precision, or approximately 0.3 mm of longitude. Keep in mind, however, the difference between precision and accuracy; the trailing 0 suggests that the accuracy is an order of magnitude less than the precision.)

**Calculate the datum shift associated with a conversion from one horizontal datum to another**- Return to the U.S. National Geodetic Survey home page.
- Follow the link labeled
**geodetic tool kit**. - At the NGS Geodetic Tool Kit page, follow the link labeled
**NADCON**(you'll be taken to an explanatory page, where you'll need to click**NADCON**again to proceed to the utility). - At the North American Datum Conversion Utility page, read the introductory paragraphs.
- At the NADCON computations form, under the heading
**compute a datum shift for a specific location**:- Select direction of conversion:
**NAD 83 to NAD27**. - Enter the NAD 83 latitude and longitude coordinates of your control station. Pay attention to format.
- Click
**Compute Datum Shift for a Single Location**.

- Select direction of conversion:
- The result should be a
**NADCON Output**report like this example. In the State College example, the adjustment from NAD 83 to NAD 27 (associated with the replacement of the old Clarke 1866 ellipsoid by the Earth-centered GRS 80 ellipsoid, caused the geographic coordinate system grid to shift nearly 7 meters South and over 23 meters West. That grid shift is reflected in the adjustment of the coordinates that specify the control point's location. Note that the point didn't move, rather, the grid shifted. How much shift occurred at your location?

GIS specialists often need to transform data from one coordinate system and/or datum to another. For example, digital data produced by tracing paper maps over a digitizing tablet need to be transformed from the tablet's non-georeferenced plane coordinate system into a georeferenced plane or spherical coordinate system that can be **georegistered** with other digital data "layers." Raw image data produced by scanning the Earth's surface from space tend to be skewed geometrically as a result of satellite orbits and other factors; to be useful these too need to be transformed into georeferenced coordinate systems. Even the point data produced by GPS receivers, which are measured as latitude and longitude coordinates based upon the WGS84 datum, often need to be transformed to other coordinate systems or datums to match project specifications. **This section describes three categories of coordinate transformations: (1) plane coordinate transformations; (2) datum transformations; and (3) map projections.**

Students who successfully complete this section of Chapter 2 should be able to:*recognize the kind of transformation that is appropriate to georegister two or more data sets.*

Some coordinate transformations are simple. For example, the transformation from non-georeferenced plane coordinates to non-georeferenced polar coordinates shown in Figure 2.18.1, below, involves nothing more than the replacement of one kind of coordinates with another.

Credit: Adapted from Iliffe, 2000

Unfortunately, most plane coordinate transformation problems are not so simple. The geometries of non-georeferenced plane coordinate systems and georeferenced plane coordinate systems tend to be quite different, mainly because georeferenced plane coordinate systems are often projected. As you know, the act of projecting a nearly-spherical surface onto a two-dimensional plane necessarily distorts the geometry of the original spherical surface. Specifically, the scale of a projected map (or an unrectified aerial photograph, for that matter) varies from place to place. So long as the geographic area of interest is not too large, however, formulae like the ones described here can be effective in transforming a non-georeferenced plane coordinate system grid to match a georeferenced plane coordinate system grid with reasonable, and measurable, accuracy. We won't go into the math of the transformations here, since the formulae are implemented within GIS software. Instead, this section aims to familiarize you with how some common transformations work and how they may be used.

In the hypothetical illustration below (Figure 2.18.2), the spatial arrangement of six control points digitized from a paper map ("before") are shown to differ from the spatial arrangement of the same points that appear in a georeferenced aerial photograph that is referenced to a different plane coordinate system grid ("after"). If, as shown, the arrangement of the two sets of points differs only in scale, rotation, and offset, a relatively simple four-parameter **similarity transformation** may do the trick. Your GIS software should derive the parameters for you by comparing the relative positions of the common points. Note that while only six control points are illustrated, ten to twenty control points are recommended (Chrisman 2002).

Click the graphic above to view a Flash animation (transform_sim.swf) in a separate browser window.

**Note:** You will need to have the Adobe Flash player installed to complete this activity. If you have not already installed the Flash player, you can download it for free from Adobe.

Sometimes a similarity transformation doesn't do the trick. For example, because paper maps expand and contract more along the paper grain than across the grain in response to changes in humidity, the scale of a paper map is likely to be slightly greater along one axis than the other. In such cases a six-parameter **affine transformation** may be used to accommodate differences in scale, rotation, and offset along each of the two dimensions of the source and target coordinate systems. This characteristic is particularly useful for transforming image data scanned from polar-orbiting satellites whose orbits trace S-shaped paths over the rotating Earth.

Click the graphic above to view a Flash animation (transform_aff.swf) in a separate browser window.

When neither similarity nor affine transformations yield acceptable results, you may have to resort to a twelve-parameter **Second-order polynomial transformation**. Their advantage is the potential to correct data sets that are distorted in several ways at once. A disadvantage is that the stability of the results depend very much upon the quantity and arrangement of control points and the degree of dissimilarity of the source and target geometries (Iliffe 2000).

Click the graphic above to view a Flash animation (transform_poly.swf) in a separate browser window.

Even more elaborate plane transformation methods, known collectively as **rubber sheeting**, optimize the fit of a source data set to the geometry of a target data set as if the source data were mapped onto a stretchable sheet.

GIS software provides a statistical measure of how well a set of transformed control points match the positions of the same points in a target data set. Put simply, Root Mean Square (RMS) Error is the average of the distances (also known as *residuals*) between each pair of control points. What constitutes an acceptably low RMS Error depends on the nature of the project and the scale of analysis.

Point locations are specified in terms of (a) their positions relative to some coordinate system grid and (b) their heights above or below some reference surface. Obviously, the elevation of a stationary point depends upon the size and shape of the reference surface (e.g., mean sea level) upon which the elevation measurement is based. In the same way, a point's position in a coordinate system grid depends on the size and shape of the surface upon which the grid is draped.** The relationship between a grid and a model of the Earth's surface is called a horizontal datum**. GIS specialists who are called upon to merge data sets produced at different times and in different parts of the world need to be knowledgeable about datum transformations.

In the U.S., the two most frequently encountered horizontal datums are the North American Datum of 1927 (NAD 27) and the North American Datum of 1983 (NAD 83). The advent of the Global Positioning System necessitated an update of NAD 27 that included (a) adoption of a geocentric ellipsoid, GRS 80, in place of the Clarke 1866 ellipsoid; and (b) correction of many distortions that had accrued in the older datum. Bearing in mind that the realization of a datum is a network of fixed control point locations that have been specified in relation to the same reference surface,** the 1983 adjustment of the North American Datum caused the coordinate values of every control point managed by the National Geodetic Survey (NGS) to change**. Obviously, the points themselves did not shift on account of the datum transformation (although they did move a centimeter or more a year due to plate tectonics). Rather, t**he coordinate system grids based upon the datum shifted in relation to the new ellipsoid**. And because local distortions were adjusted at the same time, the magnitude of grid shift varies from place to place. The illustrations below compare the magnitude of the grid shifts associated with the NAD 83 adjustment at one location and nationwide.

Credit: Dewhurst 1990

Given the irregularity of the shift, NGS could not suggest a simple transformation algorithm that surveyors and mappers could use to adjust local data based upon the older datum. Instead, NGS created a software program called NADCON (Dewhurst 1990, Mulcare 2004) that calculates adjusted coordinates from user-specified input coordinates by interpolation from a pair of 15° correction grids generated by NGS from hundreds of thousands of previously-adjusted control points.

Try out the National Geodetic Survey's NADCON tool.

The U.S. Department of Defense created the Global Positioning System (GPS) over a period of 16 years at a startup cost of about $10 billion. GPS receivers calculate their positions in terms of latitude, longitude, and height above or below the World Geodetic System of 1984 ellipsoid (WGS 84). Developed specifically for the Global Position System, WGS 84 is an Earth-centered ellipsoid which, unlike the many regional, national, and local ellipsoids still in use, minimizes deviations from the geoid worldwide. Depending on where a GIS specialist may be working, or what data he or she may need to work with, the need to transform GPS data from WGS 84 to some other datum is likely to arise. Datum transformation algorithms are implemented within GIS software as well as in the post-processing software provided by GPS vendors for use with their receivers. Some transformation algorithms yield more accurate results than others. The method you choose will depend on what choices are available to you and how much accuracy your application requires.

Unlike the plane transformations described earlier, datum transformations involve ellipsoids and are therefore three-dimensional. The simplest is the three-parameter** Molodenski transformation**. In addition to knowledge of the size and shape of the source and target ellipsoids (specified in terms of **semimajor axis**, the distance from the ellipsoid's equator to its center, and** flattening ratio**, the degree to which the ellipsoid is flattened to approximate the Earth's oblate shape), the offset between the two ellipsoids needs to be specified along X, Y, and Z axes. The window shown below (Figure 2.19.3) illustrates ellipsoidal and offset parameters for several horizontal datums, all expressed in relation to WGS 84.

Credit: Waypoint+ software utility

For larger study areas, more accurate results may be obtained using a **seven-parameter transformation** that accounts for rotation as well as scaling and offset.

Finally,** surface-fitting transformations** like the NADCON grid interpolation described above yield the best results over the largest areas.

For routine mapping applications covering relatively small geographic areas (i.e., larger than 1:25,000), the plane transformations described earlier may yield adequate results when datum specifications are unknown and when a sufficient number of appropriately distributed control points can be identified.

Latitude and longitude coordinates specify positions in a more-or-less spherical grid called the **graticule***.* Plane coordinates like the eastings and northings in the Universal Transverse Mercator (UTM) and State Plane Coordinates (SPC) systems denote positions in flattened grids. This is why georeferenced plane coordinates are referred to as *projected*, and geographic coordinates are called **unprojected***.* The mathematical equations used to transform latitude and longitude coordinates to plane coordinates are called **map projections**. Inverse projection formulae transform plane coordinates to geographic. The simplest kind of projection, illustrated in Figure 2.20.1, below, transforms the graticule into a rectangular grid in which all grid lines are straight, intersect at right angles, and are equally spaced. More complex projections yield grids in which the lengths, shapes, and spacing of the grid lines vary.

If you are a GIS practitioner you have probably faced the need to superimpose unprojected latitude and longitude data onto projected data, and vice versa. For instance, you might have needed to merge geographic coordinates measured with a GPS receiver with digital data published by the USGS that are encoded as UTM coordinates. Modern GIS software provides sophisticated tools for projecting and unprojecting data. To use such tools most effectively, you need to understand the projection characteristics of the data sets you intend to merge. We'll examine map projections in some detail elsewhere in this chapter. Here, let's simply review the characteristics that are included in the "Spatial Reference Information" section of the metadata documents that (ideally!) accompany the data sets you might wish to incorporate in your GIS. These include:

**Projection Name**Most common in the GIS realm is the Transverse Mercator, which serves as the basis of the global UTM plane coordinate system, the U.K. and proposed U.S. National Grids and many zones in the U.S. State Plane Coordinate system (SPC). Other SPC zones are based upon the Lambert Conic Conformal projection, which like many projections is named for its inventor as well as its projection category (conic) and the geometric properties it preserves (conformal). Much map data, particularly in the form of printed paper maps, are based upon "legacy" projections (like the Polyconic in the U.S.) that are no longer widely used. A much greater variety of projection types tend to be used in small scale thematic mapping than in large scale reference mapping.**Central Meridian**Although no land masses are shown, let's assume that the graticule and projected grid shown above are centered on the intersection of the equator (0 latitude) and prime meridian (0° longitude). Most map projection formulae include a parameter that allows you to center the projected map upon any longitude.**Latitude of Projection Origin**Under certain conditions, most map projection formulae allows you to specify different*aspects*of the grid. Instead of the*equatorial aspect*illustrated above, you might specify a*polar aspect*or*oblique aspect*by varying the latitude of projection origin such that one of the poles, or any latitude between the pole and the equator, is centered in the projected map. As you might imagine, the appearance of the grid changes a lot when viewed at different aspects.**Scale Factor at Central Meridian**This is the ratio of map scale along the central meridian and the scale at a*standard*meridian, where scale distortion is zero. The scale factor at the central meridian is .9996 in each of the 60 UTM coordinate system zones since each contains two standard lines 180 kilometers west and east of the central meridian. Scale distortion increases with distance from standard lines in all projected coordinate systems.**Standard Lines**Some projections, including the Lambert Conic Conformal, include parameters by which you can specify one or two*standard lines*along which there is no scale distortion caused by the act of transforming the spherical grid into a flat grid. By the same reasoning that two standard lines are placed in each UTM zone to minimize distortion throughout the zone to a maximum of one part in 1000, two standard parallels are placed in each SPC zone that is based on a Lambert projection such that scale distortion is no worse than one part in 10,000 anywhere in the zone.

Shown below in Figure 2.21.1 is the southwest corner of a 1:24,000-scale topographic map published by the United States Geological Survey (USGS). Note that the geographic coordinates (40 45' N latitude, 77° 52' 30" W longitude) of the corner are specified. Also shown, however, are ticks and labels representing two plane coordinate systems, the Universal Transverse Mercator (UTM) system and the State Plane Coordinates (SPC) system. The tick labeled "4515" represents a UTM grid line (called a "northing") that runs parallel to and 4,515,000 meters north of, the equator. Ticks labeled "258" and "259" represent grid lines that run perpendicular to the equator and 258,000 meters and 259,000 meters east, respectively, of the origin of the UTM Zone 18 North grid. Unlike longitude lines, UTM "eastings" are straight and do not converge upon the Earth's poles. All of this begs the question, **Why are multiple coordinate system grids shown on the map? Why aren't geographic coordinates sufficient?**

Credit: USGS. "State College quadrangle, Pennsylvania"

You can think of a plane coordinate system as the juxtaposition of two measurement scales. In other words, if you were to place two rulers at right angles, such that the "0" marks of the rulers aligned, you'd define a plane coordinate system. The rulers are called "axes." The absolute location of any point in the space in the plane coordinate system is defined in terms of distance measurements along the x (east-west) and y (north-south) axes. A position defined by the coordinates (1,1) is located one unit to the right, and one unit up from the origin (0,0). The UTM grid is a widely-used type of geospatial plane coordinate system in which positions are specified as eastings (distances, in meters, east of an origin) and northings (distances north of the origin).

By contrast, the geographic coordinate system grid of latitudes and longitudes consists of two *curved* measurement scales to fit the nearly-spherical shape of the Earth. As you know, geographic coordinates are specified in degrees, minutes, and seconds of arc. **Curved grids are inconvenient to use for plotting positions on flat maps. Furthermore, calculating distances, directions and areas with spherical coordinates are cumbersome in comparison with plane coordinates.** For these reasons, cartographers and military officials in Europe and the U.S. developed the UTM coordinate system. UTM grids are now standard not only on printed topographic maps but also for the geospatial referencing of the digital data that comprise the emerging U.S. "National Map."

In this section of Chapter 2 you will learn to:

- Describe the characteristics of the UTM coordinate system, including its basis in the Transverse Mercator map projection; and
- Plot UTM coordinates on a map.

The Universal Transverse Mercator system is not really universal, but it does cover nearly the entire Earth surface. Only polar areas--latitudes higher than 84° North and 80° South--are excluded. (Polar coordinate systems are used to specify positions beyond these latitudes.) The UTM system divides the remainder of the Earth's surface into 60 zones, each spanning 6° of longitude. These are numbered west to east from 1 to 60, starting at 180° West longitude (roughly coincident with the International Date Line).

The illustration above (Figure 2.22.1) depicts UTM zones as if they were uniformly "wide" from the Equator to their northern and southern limits. In fact, since meridians converge toward the poles on the globe, every UTM zone tapers from 666,000 meters in "width" at the Equator (where 1° of longitude is about 111 kilometers in length) to only about 70,000 meters at 84° North and about 116,000 meters at 80° South.

"Transverse Mercator" refers to the manner in which geographic coordinates are transformed into plane coordinates. Such transformations are called **map projections**. The illustration below (Figure 2.22.2) shows the 60 UTM zones as they appear when projected using a Transverse Mercator map projection formula that is optimized for the UTM zone highlighted in yellow, Zone 30, which spans 6° West to 0° East longitude (the prime meridian).

As you can imagine, you can't flatten a globe without breaking or tearing it somehow. Similarly, the act of mathematically transforming geographic coordinates to plane coordinates necessarily displaces most (but not all) of the transformed coordinates to some extent. Because of this, map scale varies within projected (plane) UTM coordinate system grids.

The **distortion ellipses** plotted in red help us visualize the pattern of scale distortion associated with a particular projection. Had no distortion occurred in the process of projecting the map shown in Figure 2.22.2, below, all of the ellipses would be the same size, and circular in shape. As you can see, the ellipses centered within the highlighted UTM zone are all the same size and shape. Away from the highlighted zone the ellipses steadily increase in size, although their shapes remain uniformly circular. This pattern indicates that scale distortion is minimal within Zone 30, and that map scale increases away from that zone. Furthermore, the ellipses reveal that the character of distortion associated with this projection is that shapes of features as they appear on a globe are preserved while their relative sizes are distorted. Map projections that preserve shape by sacrificing the fidelity of sizes are called **conformal** projections. The plane coordinate systems used most widely in the U.S., UTM and SPC (the State Plane Coordinates system) are both based upon conformal projections.

The Transverse Mercator projection illustrated above (Figure 2.22.2) minimizes distortion within UTM zone 30. Fifty-nine variations on this projection are used to minimize distortion in the other 59 UTM zones. In every case, distortion is no greater than 1 part in 1,000. This means that a 1,000 meter distance measured anywhere within a UTM zone will be no worse than + or - 1 meter off.

The animation linked to the illustration in Figure 2.22.3, below, shows a series of 60 Transverse Mercator projections that form the 60 zones of the UTM system. Each zone is based upon a unique Transverse Mercator map projection that minimizes distortion within that zone. Zones are numbered 1 to 60 eastward from the international date line. The animation begins with Zone 1.

Click the graphic above in Figure 2.22.3 to download and view the animation file (utm.avi, 0.5 Mb) in a separate Microsoft Media Player window.

To download and view the same animation in QuickTime format (utm.mov, 2.9 Mb), click here. Requires QuickTime, which is available free at apple.com.

**Map projections** are mathematical formulae used to transform geographic coordinates into plane coordinates. (Inverse projection formulae transform plane coordinates back into latitudes and longitudes.) "Transverse Mercator" is one of a hypothetically infinite number of such projection formulae. A visual analog to the Transverse Mercator projection appears below in Figure 2.22.4. Conceptually, the Transverse Mercator projection transfers positions on the globe to corresponding positions on a cylindrical surface, which is subsequently cut from end to end and flattened. In the illustration, the cylinder is tangent to the globe along one line, called the standard line. As shown in the little world map beside the globe and cylinder, scale distortion is minimal along the standard line and increases with distance from it. The animation linked above (Figure 2.22.3) was produced by rotating the cylinder 59 times at an increment of 6°.

In the illustration above in Figure 2.22.4, there is one standard meridian. Some projection formulae, including the Transverse Mercator projection, allow two standard lines. **Each of the 60 variations on the Transverse Mercator projection used as the foundations of the 60 UTM zones employ not one, but two, standard lines.** These two standard lines are parallel to, and 180,000 meters east and west of, each central meridian. This scheme ensures that the maximum error associated with the projection due to scale distortion will be 1 part in 1,000 (at the outer edge of the zone at the equator). The error due to scale distortion at the central meridian is 1 part in 2,500. Distortion is zero, of course, along the standard lines.

So, what does the term "transverse" mean? This simply refers to the fact that the cylinder shown above in Figure 2.22.4 has been rotated 90° from the equatorial aspect of the standard Mercator projection, in which a single standard line coincides with 0° latitude.

Credit: U.S. Geological Survey, 2004

One disadvantage of the UTM system is that multiple coordinate systems must be used to account for large entities. The lower 48 United States, for instance, spread across ten UTM zones. The fact that there are many narrow UTM zones can lead to confusion. For example, the city of Philadelphia, Pennsylvania is east of the city of Pittsburgh. If you compare the Eastings of centroids representing the two cities, however, Philadelphia's Easting (about 486,000 meters) is less than Pittsburgh's (about 586,000 meters). Why? Because although the cities are both located in the U.S. state of Pennsylvania, they are situated in two different UTM zones. As it happens, Philadelphia is closer to the origin of its Zone 18 than Pittsburgh is to the origin of its Zone 17. If you were to plot the points representing the two cities on a map, ignoring the fact that the two zones are two distinct coordinate systems, Philadelphia would appear to the west of Pittsburgh. Inexperienced GIS users make this mistake all the time. Fortunately, GIS software is getting sophisticated enough to recognize and merge different coordinate systems automatically.

The illustration in Figure 2.23.1, below, depicts the area covered by a single UTM coordinate system grid zone. Each UTM zone spans 6° of longitude, from 84° North and 80° South. Zones taper from 666,000 meters in "width" at the Equator (where 1° of longitude is about 111 kilometers in length) to only about 70,000 meters at 84° North and about 116,000 meters at 80° South. Polar areas are covered by polar coordinate systems. Each UTM zone is subdivided along the equator into two halves, north and south.

The illustration below in Figure 2.23.2 shows how UTM coordinate grids relate to the area of coverage illustrated above in Figure 2.23.1. The north and south halves are shown side by side for comparison. Each half is assigned its own **origin**. The north south zone origins are positioned to south and west of the zone. North zone origins are positioned on the Equator, 500,000 meters west of the central meridian. Origins are positioned so that every coordinate value within every zone is a positive number. This minimizes the chance of errors in distance and area calculations. By definition, both origins are located 500,000 meters west of the central meridian of the zone (in other words, the easting of the central meridian is always 500,000 meters E). These are considered "false" origins since they are located outside the zones to which they refer. UTM eastings range from 167,000 meters to 833,000 meters at the equator. These ranges narrow toward the poles. Northings range from 0 meters to nearly 9,400,000 in North zones and from just over 1,000,000 meters to 10,000,000 meters in South zones. Note that positions at latitudes higher than 84° North and 80° South are defined in Polar Stereographic coordinate systems that supplement the UTM system.

Are you ready to try your hand at positioning within a UTM coordinate system grid? The Flash application linked below lets you test your knowledge. The application asks you to click locations within a grid zone as specified by randomly generated UTM coordinates.

You will notice that the application lets you choose between easy problems (in which the locations of possible solutions are symbolized with dots) or harder problems that require you to **interpolate** solutions. **You can consider yourself to have a good working knowledge of the geographic coordinate system if you can solve at least six "hard" problems consecutively and on the first click.**

Click here to launch the UTM Coordinate System practice application. If the UTM zone fails to appear after the Flash application loads, right-click the application and select "Play" from the pop-up menu.

Note: You will need to have the Adobe Flash player installed in order to complete this exercise. If you do not already have it, you can download the Flash Player for free at adobe.com.

See the Bibliography (last page of the chapter) for further readings about the UTM grid system.

The Transverse Mercator projection provides a basis for existing and proposed national grid systems in the United Kingdom and the United States.

In the U.K., topographic maps published by the Ordnance Survey refer to a national grid of 100 km squares, each of which is identified by a two-letter code. Positions within each grid square are specified in terms of eastings and northings between 0 and 100,000 meters. The U.K. national grid is a plane coordinate system that is based upon a Transverse Mercator projection whose central meridian is 2 West longitude, with standard meridians 180 km west and east of the central meridian. The grid is typically related to the Airy 1830 ellipsoid, a relationship known as the National Grid (OSGB36®) datum. The corresponding UTM zones are 29 (central meridian 9° West) and 30 (central meridian 3° West). One of the advantages of the U.K. national grid over the global UTM coordinate system is that it eliminates the boundary between the two UTM zones.

A similar system has been proposed for the U.S. by the Federal Geographic Data Committee. The proposed "U.S. National Grid" is the same as the Military Grid Reference System (MGRS), a worldwide grid that is very similar to the UTM system. As Phil and Julianna Muehrcke (1998, p.p. 229-230) write in the 4th edition of *Map Use*, "the military [specifically, the U.S. Department of Defense] aimed to minimize confusion when using long numerical [UTM] coordinates" by specifying UTM zones and sub-zones with letters instead of numbers. Like the UTM system, the MGRS consists of 60 zones, each spanning 6° longitude. Each UTM zone is subdivided into 19 MGRS quadrangles of 8° latitude and one (quadrangle from 72° to 84° North) of 12° latitude. The letters C through X are used to designate the grid cell rows from south to north. I and O are omitted to avoid confusion with numbers. Wikipedia offers a good entry on the MGRS here.

**Fun Demo of U.K. National Grid**

A kids-friendly information sheet about the U.K. National Grid is published by the U.K. Ordnance Survey. You can find it in the National Grid for Schools link on their website.

A less-kids-friendly video can be seen below:

Credit: Ordnance Survey

This textbook is used as a resource in Penn State's Online Geospatial Education online degree and certificate programs. If this topic is interesting to you and you want to learn more about online GIS and GEOINT education at Penn State, check out ourGeospatial Education Program Office.

Shown below in Figure 2.25.1 is the southwest corner of a 1:24,000-scale topographic map published by the United States Geological Survey (USGS). Note that the geographic coordinates (40 45' N latitude, 77° 52' 30" W longitude) of the corner are specified. Also shown, however, are ticks and labels representing two plane coordinate systems, the Universal Transverse Mercator (UTM) system and the State Plane Coordinate (SPC) system. The tick labeled "1 970 000 FEET" represents a SPC grid line that runs perpendicular to the equator and 1,970,000 feet east of the origin of the Pennsylvania North zone. The origin lies far to the west of this map sheet. Other SPC grid lines, called "northings" (not shown in the illustration), run parallel to the equator and perpendicular to SPC eastings at increments of 10,000 feet. Unlike longitude lines, SPC eastings and northings are straight and do not converge upon the Earth's poles.

Credit: USGS. "State College quadrangle, Pennsylvania"

The SPC grid is a widely-used type of geospatial plane coordinate system in which positions are specified as eastings (distances east of an origin) and northings (distances north of an origin). You can tell that the SPC grid referred to in the map illustrated above is the older 1927 version of the SPC grid system because (a) eastings and northings are specified in feet and (b) grids are based upon the North American Datum of 1927 (NAD27). The 124 zones that make up the State Plane Coordinates system of 1983 are based upon NAD 83, and generally use the metric system to specify eastings and northings.

State Plane Coordinates are frequently used to georeference large scale (small area) surveying and mapping projects because plane coordinates are easier to use than latitudes and longitudes for calculating distances and areas. And because SPC zones extend over relatively smaller areas, less error accrues to positions, distances, and areas calculated with State Plane Coordinates than with UTM coordinates.

In this section you will learn to:

- describe the characteristics of the SPC system, including map projection on which it is based; and
- convert geographic coordinates to SPC coordinates.

Plane coordinate systems pretend the world is flat. Obviously, if you flatten the entire globe to a plane surface, the sizes and shapes of the land masses will be distorted, as will distances and directions between most points. If your area of interest is small enough, however, and if you flatten it cleverly, you can get away with a minimum of distortion. The basic design problem that confronted the geodesists who designed the State Plane Coordinate System, then, was to establish coordinate system zones that were small enough to minimize distortion to an acceptable level, but large enough to be useful.

The State Plane Coordinate System of 1983 (SPC) is made up of 124 zones that cover the 50 U.S. states. As shown below in Figure 2.26.1, some states are covered with a single zone while others are divided into multiple zones. Each zone is based upon a unique map projection that minimizes distortion in that zone to 1 part in 10,000 or better. In other words, a distance measurement of 10,000 meters will be at worst one meter off (not including instrument error, human error, etc.). The error rate varies across each zone, from zero along the projection's standard lines to the maximum at points farthest from the standard lines. Errors will accrue at a rate much lower than the maximum at most locations within a given SPC zone. **SPC zones achieve better accuracy than UTM zones because they cover smaller areas, and so are less susceptible to projection-related distortion. **

Most SPC zones are based on either a **Transverse Mercator** or **Lambert Conic Conformal** map projection whose parameters (such as standard line(s) and central meridians) are optimized for each particular zone. "Tall" zones like those in New York state, Illinois, and Idaho are based upon unique Transverse Mercator projections that minimize distortion by running two standard lines north-south on either side of the central meridian of each zone. "Wide" zones like those in Pennsylvania, Kansas, and California are based on unique Lambert Conformal Conic projections that run two standard parallels west-east through each zone. (One of Alaska's zones is based upon an "oblique" variant of the Mercator projection. That means that instead of standard lines parallel to a central meridian, as in the transverse case, the Oblique Mercator runs two standard lines that are tilted so as to minimize distortion along the Alaskan panhandle.)

The two types of map projections share the property of **conformality**, which means that angles plotted in the coordinate system are equal to angles measured on the surface of the Earth. As you can imagine, conformality is a useful property for land surveyors, who make their livings measuring angles. (Surveyors measure distances too, but unfortunately there is no map projection that can preserve true distances everywhere within a plane coordinate system.) Let's consider these two types of map projections briefly.

Like most map projections, the Transverse Mercator projection is actually a mathematical transformation. The illustration below in Figure 2.26.2 may help you understand how the math works. Conceptually, the Transverse Mercator projection transfers positions on the globe to corresponding positions on a cylindrical surface, which is subsequently cut from end to end and flattened. In the illustration, the cylinder is tangent to (touches) the globe along one line, the standard line (specifically, the standard meridian). As shown in the little world map beside the globe and cylinder, scale distortion is minimal along the standard line and increases with distance from it.

The **distortion ellipses** plotted in red help us visualize the pattern of scale distortion associated with a generic Transverse Mercator projection. Had no distortion occurred in the process of projecting the map shown below, all of the ellipses would be the same size, and circular in shape. As you can see, the ellipses plotted along the central meridian are all the same size and circular shape. Away from the central meridian, the ellipses steadily increase in size, although their shapes remain uniformly circular. This pattern reflects the fact that scale distortion increases with distance from the standard line. Furthermore, the ellipses reveal that the character of distortion associated with this projection is that shapes of features as they appear on a globe are preserved while their relative sizes are distorted. By preserving true angles, conformal projections like the Mercator (including its transverse and oblique variants) also preserve shapes.

SPC zones that trend west to east (including Pennsylvania's) are based on unique Lambert Conformal Conic projections. Instead of the cylindrical projection surface used by projections like the Mercator, the Lambert Conformal Conic and map projections like it employ conical projection surfaces like the one shown below in Figure 2.26.3. Notice the two lines at which the globe and the cone intersect. Both of these are standard lines; specifically, standard parallels. The latitudes of the standard parallels selected for each SPC zones minimize scale distortion throughout that zone.

In consultation with various state agencies, the National Geodetic Survey (NGS) devised the State Plane Coordinate System with several design objectives in mind. Chief among these were:

**Plane coordinates**for ease of use in calculations of distances and areas;**All positive values**to minimize calculation errors; and- A maximum error rate of
**1 part in 10,000.**

Plane coordinates specify positions in flat grids. Map projections are needed to transform latitude and longitude coordinates to plane coordinates. The designers did two things to minimize the inevitable distortion associated with map projections. First, they divided the U.S. into 124 relatively small zones. Second, they used slightly different map projection formulae for each zone. The curved, dashed red lines in the illustration below for Figure 2.27.1 represent the two standard parallels that pass through each zone. The latitudes of the standard lines are one of the parameters of the Lambert Conic Conformal projection that can be customized to minimize distortion within the zone.

Positions in any coordinate system are specified relative to an origin. **SPC zone origins are defined so as to ensure that every easting and northing in every zone are positive numbers.** As shown in the illustration below, SPC origins are positioned south of the counties included in each zone. The origins coincide with the central meridian of the map projection upon which each zone is based. The easting and northing values at the origins are not 0, 0. Instead, eastings are defined as positive values sufficiently large to ensure that every easting in the zone is also a positive number. The false origin of the Pennsylvania North zone, for instance, is defined as 600,000 meters East, 0 meters North. Origin eastings vary from zone to zone from 200,000 to 8,000,000 meters East.

**Investigating your local State Plane Coordinate System grid zone**

In this activity you will:

- read part of an authoritative manual on State Plane Coordinates;
- look up the designation of your local SPC zone (or a would-be zone if you are from a country other than the U.S.);
- investigate the parameters of the map projection upon which your local SPC zone is based; and
- use a web-based tool provided by the U.S. National Geodetic Survey to convert geographic coordinates to SPC.

Follow this link to download the manual in Portable Document Format (PDF). Read pages 1-13 (pages 11-23 of the PDF document). Also see Appendix A, beginning on p. 62 (73), which lists map projections and other parameters of each zone.

Credit: National Geodetic Survey

- Visit the National Geodetic Survey's NGS Geodetic Toolkit. Note the various programs that NGS supports for the surveying and mapping community.
- Use the dropdown menu to navigate to the State Plane Coordinates tool.
- At the State Plane Coordinates page, follow the Interactive Conversions link labeled "Find Zone."
- Look up your local SPC zone (or adopted zone) by county or position.
- You might check the result using Rick King's list or the Stem (1990) manual you downloaded earlier.

- Refer to the Stem (1990) manual you downloaded earlier. In particular, see Appendix A, pp. 63-72 (73-83). Upon which map projection is your local (or adopted) zone based?
- Note that the appendix reports a
**central meridian**and**scale factor**for each zone that is based upon a Transverse Mercator projection.**Standard parallels**are listed for zones based upon the Lambert Conformal Conic projection. In the Stem (1990) manual, "scale factor" is expressed in terms of maximum measurement error associated with each zone. - Appendix A also lists the SPC coordinates of each zone origin as well as its corresponding geographic coordinates.

Credit: National Geodetic Survey

- You'll need to know the geographic coordinates of a place of interest to complete this part of the activity. To look up the latitude and longitude associated with a U.S. place name, visit the USGS Geographic Names Information System. If you do not reside in the U.S., use the Getty Thesaurus of Geographic Names. Pay attention; be sure to choose the correct instance of the place name. Cities and towns are listed as "Populated Places" in the Geographic Names Information System and as "Inhabited Places" in the Getty Thesaurus.
- Return to the National Geodetic Survey's NGS Geodetic Toolkit.
- Use the dropdown menu to navigate to the State Plane Coordinates tool.
- At the State Plane Coordinates page, follow the Interactive Conversions link labeled "Latitude/Longitude -> SPC"
- Specify geographic (a.k.a. "geodetic") coordinates in degrees-minutes-seconds format. You do not need to show five places to the right of the decimal, as in the example illustrated above, but you do need to add a decimal point at the end of the DMS input values; the tool allows you to input a DMS value that is accurate to fractions of a second and so is programmed to expect the decimal point.
- A correct result will include a SPC Northing, Easting, zone, convergence (a correction factor for distance calculations that compensates for the tendency of meridians to converge toward the poles), and the scale factor

Latitude and longitude coordinates specify point locations within a coordinate system grid that is fitted to sphere or ellipsoid that approximates the Earth's shape and size. To display extensive geographic areas on a page or computer screen, as well as to calculate distances, areas, and other quantities most efficiently, it is necessary to flatten the Earth.

Georeferenced plane coordinate systems like the Universal Transverse Mercator and State Plane Coordinates systems (examined elsewhere in this chapter) are created by first flattening the graticule, then superimposing a rectangular grid over the flattened graticule. The first step, transforming the geographic coordinate system grid from a more-or-less spherical shape to a flat surface, involves systems of equations called **map projections**.

Many different map projection methods exist. Although only a few are widely used in large-scale mapping, the projection parameters used vary greatly. Geographic information systems professionals are expected to be knowledgeable enough to select a map projection that is suitable for a particular mapping objective. Such professionals are expected to be able to recognize the type, amount, and distribution of geometric distortion associated with different map projections. Perhaps most important, they need to know about the parameters of map projections that must be matched in order to merge geographic data from different sources. The pages that follow introduce the key concepts. The topic is far too involved to master in one section of a single chapter, however. Indeed, Penn State offers an entire online course in "Map Projections for Geospatial Professionals" (GEOG 897K, soon to be 861). If you are or plan to become, a GIS professional, you should own at least one good book on map projections. Several recommendations follow in the bibliography at the end of this chapter.

Students who successfully complete this section should be able to:

- Interpret distortion diagrams to identify geometric properties of the sphere that are preserved by a particular projection;
- Classify projected graticules by projection family.

Many types of map projections have been devised to suit particular purposes. No projection allows us to flatten the globe without distorting it, however. **Distortion ellipses** help us to visualize what type of distortion a map projection has caused, how much distortion has occurred, and where it has occurred. The ellipses show how imaginary circles on the globe are deformed as a result of a particular projection. If no distortion had occurred in the process of projecting the map shown below in Figure 2.29.1, all of the ellipses would be the same size, and circular in shape.

When positions on the graticule are transformed to positions on a projected grid, four types of distortion can occur: distortion of sizes, angles, distances, and directions. Map projections that avoid one or more of these types of distortion are said to preserve certain properties of the globe.

So-called **equal-area** projections maintain correct proportions in the sizes of areas on the globe and corresponding areas on the projected grid (allowing for differences in scale, of course). Notice that the shapes of the ellipses in the Cylindrical Equal Area projection above (Figure 2.29.1) are distorted, but the areas each one occupies are equivalent. Equal-area projections are preferred for small-scale thematic mapping, especially when map viewers are expected to compare sizes of area features like countries and continents.

The distortion ellipses plotted on the **conformal** projection shown above in Figure 2.29.2 vary substantially in size, but are all the same circular shape. The consistent shapes indicate that conformal projections (like this Mercator projection of the world) preserve the fidelity of angle measurements from the globe to the plane. In other words, an angle measured by a land surveyor anywhere on the Earth's surface can be plotted on at its corresponding location on a conformal projection without distortion. This useful property accounts for the fact that conformal projections are almost always used as the basis for large scale surveying and mapping. Among the most widely used conformal projections are the Transverse Mercator, Lambert Conformal Conic, and Polar Stereographic.

Conformality and equivalence are mutually exclusive properties. Whereas equal-area projections distort shapes while preserving fidelity of sizes, conformal projections distort sizes in the process of preserving shapes.

**Equidistant** map projections allow distances to be measured accurately along straight lines radiating from one or two points only. Notice that ellipses plotted on the Cylindrical Equidistant (Plate Carrée) projection shown above (Figure 2.29.3) vary in both shape and size. The north-south axis of every ellipse is the same length, however. This shows that distances are true-to-scale along every meridian; in other words, the property of equidistance is preserved from the two poles. See chapters 11 and 12 of the online publication *Matching the Map Projection to the Need* to see how projections can be customized to facilitate distance measurements and to effectively depict ranges and rings of activity.

**Azimuthal** projections preserve directions (azimuths) from one or two points to all other points on the map. See how the ellipses plotted on the gnomonic projection, shown above in Figure 2.29.4, vary in size and shape, but are all oriented toward the center of the projection? In this example, that's the one point at which directions measured on the globe are not distorted on the projected graticule.

Some map projections preserve none of the properties described above, but instead seek a compromise that minimizes distortion of all kinds. The example shown above in Figure 2.29.5 is the Polyconic projection, which was used by the U.S. Geological Survey for many years as the basis of its topographic quadrangle map series until it was superceded by the conformal Transverse Mercator. Another example is the Robinson projection, which is often used for small-scale thematic maps of the world.

The term "projection" implies that the ball-shaped net of parallels and meridians is transformed by casting its shadow upon some flat, or flattenable, surface. In fact, almost all map projection methods are mathematical equations. The analogy of an optical projection onto a flattenable surface is useful, however, as a means to classify the bewildering variety of projection equations devised over the past two thousand years or more.

Imagine a model globe that is translucent, and contains a bright light bulb. Imagine the light literally casting shadows of the graticule, and of the shapes of the continents, onto another surface that touches the globe. The National Geographic Society has prepared a set of animations that may help you to visualize the analogy.

As you might imagine, the appearance of the projected grid will change quite a lot depending on the type of surface it is projected onto, and how that surface is aligned with the globe. The three surfaces shown above in Figure 2.30.1--the disk-shaped plane, the cone, and the cylinder--represent categories that account for the majority of projection equations that are encoded in GIS software. All three are shown in their normal *aspects.* The plane often is centered upon a pole. The cone is typically aligned with the globe such that its line of contact (tangency) coincides with a parallel in the mid-latitudes. And the cylinder is frequently positioned tangent to the equator (unless it is rotated 90°, as it is in the Transverse Mercator projection). The following illustrations in Figure 2.30.2 show some of the projected graticules produced by projection equations in each category.

**Cylindric** projection equations yield projected graticules with straight meridians and parallels that intersect at right angles. The example shown above at top left in Figure 2.30.2 is a Cylindrical Equidistant (also called Plate Carrée or geographic) in its normal equatorial aspect.

**Pseudocylindric** projections are variants on cylindrics in which meridians are curved. The result of a Sinusoidal projection is shown above at top right of Figure 2.30.2.

**Conic** projections yield straight meridians that converge toward a single point at the poles, parallels that form concentric arcs. The example shown above, at bottom left in Figure 2.30.2, is the result of an Albers Conic Equal Area, which is frequently used for thematic mapping of mid latitude regions.

**Planar** projections also yield meridians that are straight and convergent, but parallels form concentric circles rather than arcs. Planar projections are also called **azimuthal** because every planar projection preserves the property of azimuthality. The projected graticule shown above at bottom right of Figure 2.30.2 is the result of an Azimuthal Equidistant projection in its normal polar aspect.

Appearances can be deceiving. It's important to remember that the look of a projected graticule depends on several projection parameters, including latitude of projection origin, central meridian, standard line(s), and others. Customized map projections may look entirely different from the archetypes described above.

The Interactive Album of Map Projections 2.0 is an application developed by the Penn State Online Geospatial Education Programs and is an update of an earlier site that was inspired by the USGS Professional Paper 1453, An Album of Map Projections, by John P. Snyder and Philip M. Voxland.

Flex Projector is a free, open source software program developed in Java that supports many more projections and variable parameters than the Interactive Album. Bernhard Jenny of the Institute of Cartography at ETH Zurich created the program with assistance from Tom Patterson of the US National Park Service. You can download Flex Projector from flexprojector.com

Those who wish to explore map projections in greater depth than is possible in this text might wish to visit an informative page published by the International Institute for Geo-Information Science and Earth Observation (Netherlands), which is known by the legacy acronym ITC.

In this chapter, we've explored several connotations of the term scale. Scale is synonymous with scope when it is used to describe the extent of a phenomenon. In this sense, "large scale" means "large area." Specialists in geographic information often use the term differently, however. **Map scale** refers to the relative sizes of features on a map and of corresponding objects on the ground. In this context, "large scale" implies "small area." Large scale also implies greater detail and greater accuracy, an important point to keep in mind when using maps as sources for GIS databases. Map scale is defined mathematically as the proportion of map distance to ground distance. I hope you are now prepared to use scale equations to calculate map scale.

Scale can also be thought of as a reference system for measurement. Locations on the globe are specified with reference to the **geographic coordinate system** of latitudes and longitudes. **Plane coordinates** are often preferred over geographic coordinates because they ease calculations of distance, area, and other quantities. Georeferenced plane coordinate systems like **UTM and SPC** are established by first flattening the graticule, then superimposing a plane coordinate grid. The mathematical equations used to transform geographic coordinates into plane coordinates are called **map projections**. Both plane and geographic coordinate system grids are related to approximations of the Earth's size and shape called **ellipsoids**. Relations between grids and ellipsoids are called **horizontal datums**.

**Horizontal datum** is an elusive concept for many GIS practitioners. It is relatively easy to visualize a horizontal datum in the context of unprojected geographic coordinates. Simply drape the latitude and longitude grid over an ellipsoid and there's your horizontal datum. It is harder to think about datum in the context of a projected coordinate grid like UTM and SPC, however. Think of it this way: First drape the latitude and longitude grid on an ellipsoid. Then project that grid to a 2-D plane surface. Then, superimpose a rectangular grid of eastings and northings over the projection, using control points to georegister the grids. There you have it--a projected coordinate grid based upon a horizontal datum.

Numerous coordinate systems, datums, and map projections are in use around the world. Because we often need to combine georeferenced data from various sources, GIS professionals need to be able to **georegister** two or more data sets that are based upon different coordinate systems, datums, and/or projections. **Transformations**, including coordinate transformations, datum transformations, and map projections, are the mathematical procedures used to bring diverse data into alignment. Characteristics of the coordinate systems, datums, and projections considered in this text are outlined in the following tables.

(many other national and local systems are in use)

Coordinate System | Units | Extent | Projection Basis |
---|---|---|---|

Geographic | Angles (expressed as degrees, minutes, seconds or decimal degrees). | Global | None |

UTM | Distances (meters) | Near-global (8430' N, 80° 30' S) | Unique Transverse Mercator projection for each of 60 zones |

State Plane Coordinates | Distances (meters in SPCS 83, feet in SPCS 27) | U.S. | Unique Transverse Mercator or Lambert Conformal Conic projection for each of 123 zones (plus Oblique Mercator for Alaska panhandle) |

(many other national and local systems are in use)

Datum |
Horiztonal or vertical |
Optimized for |
Reference surface |
---|---|---|---|

NAD 27 | Horizontal | North America | Clarke 1866 ellipsoid |

NAD 83 | Horizontal | North America | GRS 80 ellipsoid |

WGS 84 | Horizontal | World | WGS 84 ellipsoid |

NAVD 88 | Vertical | North America | Sea level measured at coastal tidal stations |

(many other national and local systems are in use)

Projection name |
Properties preserved |
Class |
Distortion |
---|---|---|---|

Mercator | Conformal | Cylindrical | Area distortion increases with distance from standard parallel (typically equator). |

Transverse Mercator | Conformal | Cylindrical | Area distortion increases with distance from standard meridian. |

Lambert Conformal Conic | Conformal | Conic | Area distortion increases with distance from one or two standard parallels. |

Plate Carrée (sometimes called "Geographic" projection) | Equidistant | Cylindrical | Area and shape distortion increases with distance from standard parallel (typically equator). |

Albers Equal-Area Conic | Equivalent | Conic | Shape distortion increases with distance from one or two standard parallels. |

*Compiled from Snyder, 1997*

3-D Software (2005). *Map projections pages*. Retrieved January 8, 2005, from http://www.3dsoftware.com/Cartography/

American Congress on Surveying and Mapping (n. d.). *The North American Datum of 1983. A collection of papers describing the planning and implementation of the readjustment of the North American horizontal network. *Monograph No. 2.

Burkard, R. K. et al. (1959-2002).* Geodesy for the layman*. Retrieved October 29, 2003, from the National Imagery and Mapping Agency Web site http://www.ngs.noaa.gov/PUBS_LIB/Geodesy4Layman/toc.htm

Chem-Nuclear Systems, Inc. (1993). *Site screening interim report: Stage two -- regional disqualification.* Harrisburg PA.

Chrisman, N. (2002). *Exploring geographic information systems (2nd ed.)*. New York: John Wiley & Sons.

Clarke, K. (1995). *Analytical and computer cartography (2nd ed.)*. Upper Saddle River, NJ: Prentice Hall.

Dana, P. H. (1998). *Coordinate systems overview. The Geographer's Craft Project*. Retrieved June 25, 2004, from The University of Colorado at Boulder, Department of Geography Web site: http://www.colorado.edu/geography/gcraft/notes/coordsys/coordsys.html

Dana, P. H. (1999). *Geodetic datums overview. The Geographer's Craft Project*. Retrieved June 25, 2004, from The University of Colorado at Bolder, Department of Geography Web site: http://www.colorado.edu/geography/gcraft/notes/datum/datum.html

Dewhurst, W. T. (1990). *NADCON: The application of minimum-curvature-derived surfaces in the transformation of positional data from the North American datum of 1927 to the North American datum of 1983*. NOAA Technical Memorandum NOS NGS 50. Retrieved January 1, 2005, from http://www.ngs.noaa.gov/PUBS_LIB/NGS50.pdf

Doyle, D. (2004, February). NGS geodetic toolkit, Part 7: Computing state plane coordinates. *Professional Surveyor Magazine*, 24:, 34-36.

Dutch, S. (2003). The Universal Transverse Mercator System. Retrieved January 9, 2008, from http://www.uwgb.edu/DutchS/FieldMethods/UTMSystem.htm

Federal Geographic Data Committee. (December 2001). *United States National Grid*. Retrieved May 8, 2006, from http://www.fgdc.gov/standards/projects/FGDC-standards-projects/usng/fgdc_std_011_2001_usng.pdf

Hildebrand, B. (1997). *Waypoint+*. Retrieved January 1, 2005, from http://www.tapr.org/~kh2z/Waypoint/ (since retired).

Iliffe, J.C. (2000). *Datums and map projections for remote sensing, GIS and surveying*. Caithness, Scotland: Whittles Publishing. Distributed in U.S. by CRC Press.

Larrimore, C. (2002). *NGS Geodetic Toolkit.* Retrieved October 26, 2004, from http://www.ngs.noaa.gov/TOOLS

Muehrcke, P. C. & Muehrcke, J. O. (1992). *Map use (3rd ed.)*. Madison WI: JP Publications.

Muehrcke, P. C. & Muehrcke, J. O. (1998). *Map use (4th ed.).* Madison WI: JP Publications.

Mulcare, D. M. (2004). The National Geodetic Survey NADCON Tool. *Professional Surveyor Magazine*, February, pp. 28-33.

National Geodetic Survey. (n.d.).* North American datum conversion utility*. Retrieved April 2004, from http://www.ngs.noaa.gov/TOOLS/Nadcon/Nadcon.html

National Geodetic Survey. (1997). *Image generated from 15'x15' geoid undulations covering the planet Earth. * Retrieved 1999, from http://www.ngs.noaa.gov/GEOID/geo-index.html (since retired).

National Geodetic Survey. (2004). *Coast and geodetic survey historical image collection*. Retrieved June 25, 2004, from http://www.photolib.noaa.gov/cgs/index.html

National Geographic Society (1999). *Round earth, flat maps.* Retrieved April 18, 2006, from http://www.nationalgeographic.com/features/2000/exploration/projections/index.html

Ordnance Survey (2000). *National GPS network information. 7: Transverse mercator map projections*. Retrieved August 27, 2004, from http://www.gps.gov.uk/guide7.asp

Robinson, A. et al. (1995). *Elements of cartography (5th ed.)*. New York: John Wiley & Sons.

Robinson, A. H. & Snyder, J. P. (1997). *Matching the map projection to the need.* Retrieved January 8, 2005, from the Cartography and Geographic Information Society and the Pennsylvania State University web site: https://courseware.e-education.psu.edu/projection/

Slocum, T. A., McMaster, R. B., Kessler, F, C., & Howard, H. H. (2005). *Thematic cartography and visualization (2nd ed.)*. Upper Saddle River, NJ: Prentice Hall.

Smith, J.R. (1988). *Basic geodesy*. Rancho Cordova CA: Landmark Enterprises.

Snyder, J. P. (1987). *Map projections: A working manua*l (U.S. Geological Survey Professional Paper No. 1395). Washington DC: United States Government Printing Office.

Snyder, J. P. (1987). *Map projections: A working manual.* (USGS Professional Paper No. 1395). Washington DC: U.S. Geological Survey

(Electronic versions available at http://pubs.er.usgs.gov/djvu/PP/PP_1395.pdf)

(Electronic versions available at http://pubs.er.usgs.gov/djvu/PP/PP_1395.pdf)

Snyder, J. P. & Voxland P. M. (1989). *An album of map projections* (U.S. Geological Survey Professional Paper No. 1453). Washington DC: United States Government Printing Office.

Snyder, J. P. & Voxland, P. M. (1994). An album of map projections. (USGS Professional Paper No. 1453). Washington DC: U.S. Geological Survey. (ordering information published at http://erg.usgs.gov/isb/pubs/factsheets/fs08799.html)

Stem, J. E. (1990). *State Plane Coordinate System of 1983 (NOAA Manual NOS NGS 5)*. Rockville, MD: National Geodetic Information Center.

United States Geological Survey (2001). *The universal transverse mercator grid. Fact sheet 077-01*. Retrieved June 30, 2004, from http://mac.usgs.gov/mac/isb/pubs/factsheets/fs07701.html (since retired).

United States Geological Survey (2003). *National mapping program standards*. Retrieved October 29, 2005, from http://rockyweb.cr.usgs.gov/nmpstds/nmas647.html (since retired).

USGS. "State College Quadrangle" [map]. 7.5 minute series. Washington, D.C.: USGS, 1962.

Van Sickle, J. (2004). *Basic GIS coordinates*. Boca Raton FL: CRC Press.

Wikipedia. The free encyclopedia. (2006). *World geodetic system*. Retrieved May 8, 2006, from http://en.wikipedia.org/wiki/WGS84

Wolf, P. R. & Brinker, R. C. (1994) *Elementary Surveying (9th ed.).* New York NY: HarperCollins.