Select a custom projected CRS for the area of interest
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.
geo_select_aeq(shp)
shp |
A spatial object with a geographic (WGS84) coordinate system |
The function is based on this stackexchange answer: https://gis.stackexchange.com/questions/121489
sp::bbox(routes_fast) new_crs <- geo_select_aeq(routes_fast) rf_projected <- sp::spTransform(routes_fast, new_crs) sp::bbox(rf_projected) line_length <- rgeos::gLength(rf_projected, byid = TRUE) plot(line_length, rf_projected$length) shp <- zones_sf geo_select_aeq(shp)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.