stplanr (version 0.5.0)

crs_select_aeq: Select a custom projected CRS for the area of interest

Description

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.

Usage

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

See Also

Other geo: bbox_scale(), gclip(), geo_bb_matrix(), geo_bb(), mapshape_available(), mapshape(), quadrant(), reproject()

Examples

# NOT RUN {
data("routes_fast")
new_crs <- geo_select_aeq(routes_fast)
plot(routes_fast)
rf_projected <- sp::spTransform(routes_fast, new_crs)
plot(rf_projected)
sp::bbox(rf_projected)
line_length <- rgeos::gLength(rf_projected, byid = TRUE)
plot(line_length, rf_projected$length)
cor(line_length, rf_projected$length)
# }