Learn R Programming

tmaptools (version 2.0)

set_projection: Set and get the map projection

Description

The function set_projection sets the projection of a shape file. It is a convenient wrapper of st_transform (or st_transform_proj, see details) and projectRaster with shortcuts for commonly used projections. The projection can also be set directly in the plot call with tm_shape. This function is also used to set the current projection information if this is missing. The function get_projection is used to get the projection information.

Usage

set_projection(shp, projection = NA, current.projection = NA,
  overwrite.current.projection = FALSE)

get_projection(shp, guess.longlat = FALSE, output = c("character", "crs", "epsg", "CRS"))

Arguments

shp

shape object, which is an object from a class defined by the sf, sp, or raster package.

projection

new projection. See get_proj4 for options. This argument is only used to transform the shp. Use current.projection to specify the current projection of shp.

current.projection

the current projection of shp. See get_proj4 for possible options. Only use this if the current projection is missing or wrong.

overwrite.current.projection

logical that determines whether the current projection is overwritten if it already has a projection that is different.

guess.longlat

if TRUE, it checks if the coordinates are within -180/180 and -90/90, and if so, it returns the WGS84 longlat projection (which is get_proj4("longlat")).

output

output format of the projection. One of "character", "crs" (from sf package), "epsg" or "CRS" (from sp/rgdal package)

Value

set_projection returns a (transformed) shape object with updated projection information. get_projection returns the PROJ.4 character string of shp.

Details

For sf objects, set_projection first tries to use sf::st_transform, which uses the GDAL API. For some projections, most notably Winkel Tripel ("wintri"), is doesn't work. In these cases, set_projection will use lwgeom::st_transform_proj, which uses the PROJ.4 API.

For raster objects, the projection method is based on the type of data. For numeric layers, the bilinear method is used, and for categorical layers the nearest neighbor. See projectRaster for details.