Spatial Referencing in GIS
Geodesy, Map Projections and
      Coordinate Systems
• Geodesy - the shape of the earth and
  definition of earth datums
• Map Projection - the transformation of a
  curved earth to a flat map
• Coordinate systems - (x,y,z) coordinate
  systems for map data
                    Learning Objectives:
    By the end of this class you should be able to:
•    describe the role of geodesy as a basis for earth datums
•    list the basic types of map projection
•    identify the properties of common map projections
•    properly use the terminology of common coordinate systems
•    use spatial references in ArcMap so that geographic data is
     properly displayed
     – determine the spatial reference system associated with a feature class or
       data frame
     – use ArcGIS to convert between coordinate systems
• calculate distances on a spherical earth and in a projected
  coordinate system
                Readings: Introduction
http://resources.arcgis.com/en/help/getting-started/articles/026n0000000s000000.htm
             Readings: Further Detail
http://resources.arcgis.com/en/help/main/10.2/#/What_are_map_projections/003r00000001000000/
   Types of Coordinate Systems
• (1) Global Cartesian coordinates (x,y,z) for
  the whole earth
• (2) Geographic coordinates (f, l, z)
• (3) Projected coordinates (x, y, z) on a local
  area of the earth’s surface
• The z-coordinate in (1) and (3) is defined
  geometrically; in (2) the z-coordinate is
  defined gravitationally
Global Cartesian Coordinates (x,y,z)
                          Z
      Greenwich
       Meridian
                             O
                         •                  Y
         X
                                      Equator
          Next 7 slides are from Dr Irmak
Geographic Coordinates (f, l, z)
• Latitude (f) and Longitude (l) defined
  using an ellipsoid, an ellipse rotated about
  an axis
• Elevation (z) defined using geoid, a surface
  of constant gravitational potential
• Earth datums define standard values of the
  ellipsoid and geoid
       Shape of the Earth
We think of the       It is actually a spheroid,
earth as a sphere    slightly larger in radius at
                    the equator than at the poles
                           Ellipse
An ellipse is defined by:
                                           Z
Focal length = 
Distance (F1, P, F2) is
constant for all points
on ellipse                                 b
When  = 0, ellipse = circle           O           a        X
                                F1                   F2
For the earth:
Major axis, a = 6378 km
Minor axis, b = 6357 km
Flattening ratio, f = (a-b)/a    P
                    ~ 1/300
 Ellipsoid or Spheroid
Rotate an ellipse around an axis
                Z
             b
           a O a            Y
          Rotational axis
             Standard Ellipsoids
E llipsoid    M ajor       M inor       F lattening
              axis, a (m ) axis, b (m ) ratio, f
C larke       6,378,206     6,356,584     1/294.9 8
(1866)
G R S 80      6,378,137     6,356,752     1/298.5 7
Ref: Snyder, Map Projections, A working manual, USGS
Professional Paper 1395, p.12
           Geodetic Datums
• World Geodetic System (WGS) – is a
  global system for defining latitude and
  longitude on earth independently of tectonic
  movement (military)
• North American Datum (NAD) – is a
  system defined for locating fixed objects on
  the earth’s surface and includes tectonic
  movement (civilian)
     Horizontal Earth Datums
• An earth datum is defined by an ellipse and
  an axis of rotation
• NAD27 (North American Datum of 1927)
  uses the Clarke (1866) ellipsoid on a non
  geocentric axis of rotation
• NAD83 (NAD,1983) uses the GRS80
  ellipsoid on a geocentric axis of rotation
• WGS84 (World Geodetic System of 1984)
  uses GRS80, almost the same as NAD83
Adjustments of the NAD83 Datum
                 Slightly different (f, l) for benchmark
                Continuously Operating Reference System
                Canadian Spatial Reference System
                National Spatial Reference System
                High Accuracy Reference Network
Representations of the Earth
     Mean Sea Level is a surface of constant
     gravitational potential called the Geoid
Sea surface                         Ellipsoid
                                         Earth surface
 Geoid
     Cutting Plane of a Meridian
                 P
                          Prime Meridian
                          Equator
Meridian
   Definition of Longitude, l
l = the angle between a cutting plane on the prime meridian
and the cutting plane on the meridian through the point, P
                             180°E, W
                     -150°              150°
            -120°                              120°
          90°W                                   90°E
          (-90 °)                               (+90 °)
             -60°            P l               -60°
                    -30°                 30°
                              0°E, W
Latitude and Longitude on a Sphere
  Greenwich                  Z            Meridian of longitude
                         N
   meridian                                  Parallel of latitude
     l=0°
                                      P
                                     •
                                                     l - Geographic longitude
                                                      - Geographic latitude
    W                O                          E
                                             •        Y
                         l       R
                                                     R - Mean earth radius
          •
                  Equator    =0°
              •                                      O - Geocenter
   X
 Length on Meridians and Parallels
(Lat, Long) = (f, l)
Length on a Meridian:
AB = Re Df                              R
(same for all latitudes)         R Dl       D
                                     C
                            Re    Df B
                                 Re
Length on a Parallel:
                                      A
CD = R Dl = Re Dl Cos f
(varies with latitude)
Reason-parallels converge
At poles
Example: What is the length of a 1º increment along
on a meridian and on a parallel at 30N, 90W?
Radius of the earth = 6370 km.
Solution:
• A 1º angle has first to be converted to radians
p radians = 180 º, so 1º = p/180 = 3.1416/180 = 0.0175 radians
• For the meridian, DL = Re Df = 6370 * 0.0175 = 111 km
• For the parallel, DL = Re Dl Cos f
                       = 6370 * 0.0175 * Cos 30
                       = 96.5 km
• Parallels converge as poles are approached
                        Curved Earth Distance
                                      (from A to B)
 Shortest distance is along a “Great
 Circle”                                                                Z
 A “Great Circle” is the intersection of
 a sphere with a plane going through                                             B
 its center.                                                  A
 1. Spherical coordinates converted
 to Cartesian coordinates.                                             
 2. Vector dot product used to                                         •                         Y
 calculate angle  from latitude and
 longitude
                                                X
 3. Great circle distance is R,
 where R=6378.137 km2
                           1
         Dist = R cos           [sin f A sin f B  cos f A cos f B cos( l A  l B )]
Ref: Meyer, T.H. (2010), Introduction to Geometrical and Physical Geodesy, ESRI Press, Redlands, p. 108
      Three systems for measuring
               elevation
Orthometric heights     Ellipsoidal heights     Tidal heights
(land surveys, geoid)       (lidar, GPS)      (Sea water level)
Conversion among these height systems has some uncertainty
   Representations of the-
       Earth Surfaces
        Mean Sea Level is a surface of constant
        gravitational potential called the Geoid
   Sea surface                             Spheroid
                                                Topographic
Geoid                                           surface
           Geoid is a surface of constant
           gravity.
                 Earth Surfaces-
                 Geoid and Ellipsoid
GPS Heights are from
Ellipsoid.                     Topographic surface,
Need to account Gravity
Anomaly,
So Heights are from Mean sea
level not spheroid.
         Ocean
                               Spheroid
                      Geoid               Gravity Anomaly
      Gravity anomaly is the elevation difference between
      a standard shape of the earth (ellipsoid) and a
      surface
      of constant gravitational potential (geoid)
       Vertical Earth Datums
• A vertical datum defines elevation, z
• NGVD29 (National Geodetic Vertical
  Datum of 1929)
• NAVD88 (North American Vertical Datum
  of 1988)
• takes into account a map of gravity
  anomalies between the ellipsoid and the
  geoid
           Importance of geodetic datums
                       NAVD88 – NGVD29 (cm)
                                                   NGVD29 higher
                                                      in East
More than 1 meter difference
            NAVD88 higher
               in West
          Orthometric datum height shifts are significant relative to BFE
               accuracy, so standardization on NAVD88 is justified
  Geodesy and Map Projections
• Geodesy - the shape of the earth and
  definition of earth datums
• Map Projection - the transformation of a
  curved earth to a flat map
• Coordinate systems - (x,y) coordinate
  systems for map data
 Earth to Globe to Map
     Map Scale:           Map Projection:
Representative Fraction       Scale Factor
  = Globe distance            Map distance
                          =
    Earth distance            Globe distance
    (e.g. 1:24,000)            (e.g. 0.9996)
Geographic and Projected Coordinates
    (f, l)                    (x, y)
             Map Projection
        Types of Projections
• Conic (Albers Equal Area, Lambert
  Conformal Conic) - good for East-West
  land areas
• Cylindrical (Transverse Mercator) - good
  for North-South land areas
• Azimuthal (Lambert Azimuthal Equal Area)
  - good for global views
        Types of Projections
• Conic (Albers Equal Area, Lambert
  Conformal Conic) - good for East-West
  land areas
• Cylindrical (Transverse Mercator) - good
  for North-South land areas
• Azimuthal (Lambert Azimuthal Equal Area)
  - good for global views
                     Map projection
• Map projection transforms a
  position on the Earth’s surface
  identified by latitude/longitude
  into a position in Cartesian
  coordinates (x, y)
• Impossible to have the same scale
  everywhere, or for the pixel size
  to be perfectly constant.
• But projections can preserve
  certain properties.
Conic Projections
   (Albers, Lambert)
      Cylindrical Projections
              (Mercator)
Transverse
                           Oblique
Azimuthal
 (Lambert)
Albers Equal Area Conic Projection
Lambert Conformal Conic Projection
Universal Transverse Mercator Projection
Lambert Azimuthal Equal Area Projection
           Web Mercator Projection
                        (used for ESRI Basemaps)
Web Mercator is one of the
most popular coordinate
systems used in web
applications because it fits the
entire globe into a square area
that can be covered by 256 by
256 pixel tiles.
The spatial reference for the
ArcGIS Online / Google
Maps / Bing Maps tiling
scheme is WGS 1984 Web
Mercator (Auxiliary
Sphere).
                                     Web Mercator
                                      Parameters
                                         (20037, 19971 km)
                                         = earth rad * Π
  Standard Parallel              (0,0)
6357 km
      6378 km
Earth radius          Central Meridian
         Projections Preserve Some
             Earth Properties
• Area - correct earth surface area (Albers
  Equal Area) important for mass balances
• Shape - local angles are shown correctly
  (Lambert Conformal Conic)
• Direction - all directions are shown correctly
  relative to the center (Lambert Azimuthal
  Equal Area)
• Distance - preserved along particular lines
• Some projections preserve two properties
   Projection and Datum
Two datasets can differ in both the
projection and the datum, so it is
important to know both for every
dataset.
  Geodesy and Map Projections
• Geodesy - the shape of the earth and
  definition of earth datums
• Map Projection - the transformation of a
  curved earth to a flat map
• Coordinate systems - (x,y) coordinate
  systems for map data
         Coordinate Systems
• Universal Transverse Mercator (UTM) - a
  global system developed by the US Military
  Services
• State Plane Coordinate System - civilian
  system for defining legal boundaries
• Oregon uses a conic - a statewide
  coordinate system
          Coordinate System
     A planar coordinate system is defined by a pair
     of orthogonal (x,y) axes drawn through an origin
                     Origin                             X
                                       (xo,yo)
(fo,lo)
Universal Transverse
Mercator
• Uses the Transverse Mercator projection
• Each zone has a Central Meridian (lo),
  zones are 6° wide, and go from pole to pole
• 60 zones cover the earth from East to West
• Reference Latitude (fo), is the equator
• (Xshift, Yshift) = (xo,yo) = (500000, 0) in
  the Northern Hemisphere, units are meters
UTM Zone 14
          -99°
     -102°   -96°
     6°
   Origin
                            Equator
  -120°          -90 °   -60 °
  State Plane Coordinate System
• Defined for each State in the United States
• East-West States (e.g. Texas) use Lambert
  Conformal Conic,
• North-South States (e.g. California) use
  Transverse Mercator
• Oregon SPCS uses Lambert Conic Conformal -
  two zones (North, South) to give accurate
  representation
• Greatest accuracy for local measurements
State Plane Zones
  Oregon’s State Plane Coordinate
              System
• Lambert conformal
  conic
• Projection based on
  the mathematics of a
  cone intersecting a
  sphere at two
  different latitudes
• North-south limited to
  158 miles – minimize
  distortion
           Oregon’s Projection
• When developable surface is unwrapped the
  parallels appear as acres of concentric circles
  and meridians project toward center
• Along the standard parallels the distances on
  the projection area at a scale of 1:1
• Between parallels the scale is smaller on the
  map than on the ground
• Outside the parallels a projected distance is
  larger
Oregon’s Projection
Oregon State Plane Coordinate
  System: North and South
Oregon State Plane Coordinate System: North and South
Oregon State Projection
Oregon State Projection
    Terms to describe a projection
•   Latitude of Grid Origin (TM)
•   Standard Parallel (LCC)
•   Latitude of Local Origin (OM/RSO)
•   Central Meridian (TM, LCC)
•   Longitude of Local Origin (OM/RSO)
•   Angle Skew or Azimuth (OM/RSO)
•   False Northing, False Easting
•   Scale
 ArcGIS Spatial Reference Frames
• Defined for a feature
  dataset in ArcCatalog
• XY Coordinate System
   – Projected
   – Geographic
• Z Coordinate system
• Domain, resolution and
  tolerance
 Horizontal Coordinate Systems
• Geographic             • Projected coordinates
  coordinates (decimal     (length units, ft or
  degrees)                 meters)
   Vertical Coordinate Systems
• None for 2D
  data
• Necessary for
  3D data
ArcGIS .prj files
          Summary Concepts
• The spatial reference of a dataset comprises
  datum, projection and coordinate system.
• For consistent analysis the spatial reference
  of data sets should be the same.
• ArcGIS does projection on the fly so can
  display data with different spatial references
  properly if they are properly specified.
• ArcGIS terminology
  – Define projection. Specify the projection for
    some data without changing the data.
  – Project. Change the data from one projection
    to another.
    Summary Concepts (Cont.)
• Two basic locational systems: geometric or
  Cartesian (x, y, z) and geographic or
  gravitational (f, l, z)
• Mean sea level surface or geoid is
  approximated by an ellipsoid to define an
  earth datum which gives (f, l) and distance
  above geoid gives (z)
    Summary Concepts (Cont.)
• To prepare a map, the earth is first reduced
  to a globe and then projected onto a flat
  surface
• Three basic types of map projections: conic,
  cylindrical and azimuthal
• A particular projection is defined by a
  datum, a projection type and a set of
  projection parameters
    Summary Concepts (Cont.)
• Standard coordinate systems use particular
  projections over zones of the earth’s surface
• Types of standard coordinate systems:
  UTM, State Plane
• Web Mercator coordinate system (WGS84
  datum) is standard for ArcGIS basemaps