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.
A spatial object with a geographic (WGS84) coordinate system
The function is based on this stackexchange answer: http://gis.stackexchange.com/questions/121489
1 2 3 4 5 6 7 8 9
data("routes_fast") new_crs <- crs_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)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.