Select a custom projected CRS for the area of interest

Source: R/crs-funs.R

This function takes a spatial object with a geographic (WGS84) CRS and returns a custom projected CRS focussed on the centroid of the object. This function is especially useful for using units of metres in all directions for data collected anywhere in the world.

crs_select_aeq(shp)

Arguments

shp

A spatial object with a geographic (WGS84) coordinate system

Details

The function is based on this stackexchange answer: http://gis.stackexchange.com/questions/121489

Examples

data("routes_fast")
new_crs <- crs_select_aeq(routes_fast)
plot(routes_fast)
#> Error in as.double(y): cannot coerce type 'S4' to vector of type 'double'
rf_projected <- sp::spTransform(routes_fast, new_crs)
plot(rf_projected)
#> Error in as.double(y): cannot coerce type 'S4' to vector of type 'double'
sp::bbox(rf_projected)
#>         min      max
#> x -1545.979 1086.642
#> y -1833.567 1275.756
line_length <- rgeos::gLength(rf_projected, byid = TRUE)
plot(line_length, rf_projected$length)
cor(line_length, rf_projected$length)
#> [1] NA