R/lookup.projection.params.r

 #' @export 
lookup.projection.params = function(x) {
    # for examples, see: 
    # http://www.progonos.com/furuti/MapProj/Normal/TOC/cartTOC.html
    out  = switch( x,
      utm19		= "+proj=utm +ellps=WGS84 +zone=19 +units=km ",
      utm20             = "+proj=utm +ellps=WGS84 +zone=20 +units=km ", 
      utm21             = "+proj=utm +ellps=WGS84 +zone=21 +units=km ", 
      utm20.substrate   = "+proj=utm +datum=NAD83 +zone=20 +units=km",
      lambert.conic.equidist    = "+proj=eqdc +ellps=WGS84 +lon_0=63W +lat_0=45N +lat_1=43N +lat_2=47N +units=km ",
      lambert.conic             = "+proj=lcc +ellps=WGS84  +lon_0=63W +lat_0=45N +lat_1=43N +lat_2=47N +units=km ",
      lambert.conic.4t          = "+proj=lcc +ellps=WGS84  +lon_0=63W +lat_0=47N +lat_1=45N +lat_2=49N +units=km ",
      lambert.conic.canada.east = "+proj=lcc +ellps=WGS84  +lon_0=62W +lat_0=45N +lat_1=43N +lat_2=47N +units=km ",
      lambert.conic.us.4x      = "+proj=lcc +ellips=WGS84  +lon_0=68W +lat_0=42N +lat_1=41N +lat_2=43N +units=km ",
      tmercator         = "+proj=tmerc +ellps=WGS84 +x_0=500000 +y=-400000 +lon_0=90w +units=km "
    )
    return ( out )
  }
jae0/emgis documentation built on May 28, 2019, 9:58 p.m.