R/projection_proj4string.R

In jae0/aegis: Assessment of ecosystems using GIS

 #' @export
projection_proj4string = function(x) {
    # for examples, see:
    # http://www.progonos.com/furuti/MapProj/Normal/TOC/cartTOC.html
    out  = switch( x,
      lonlat_wgs84 = "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs",
      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",
      omerc_nova_scotia = "+proj=omerc +lat_0=44.0 +lonc=-63.0 +gamma=0.0 +k=1 +alpha=325 +x_0=0 +y_0=0 +ellps=WGS84 +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 )
  }