Many GIS functions (e.g. finding the area)

reproject(shp, crs = geo_select_aeq(shp))

Arguments

shp

A spatial object with a geographic (WGS84) coordinate system

crs

An optional coordinate reference system (if not provided it is set automatically by geo_select_aeq()).

See also

Other geo: bbox_scale(), geo_bb_matrix(), geo_bb(), quadrant()

Examples

data(routes_fast)
rf_aeq <- reproject(routes_fast[1:3, ])
#> Transforming to CRS +proj=aeqd +lat_0=53.82763236 +lon_0=-1.53163002 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs
rf_osgb <- reproject(routes_fast[1:3, ], 27700)
#> Warning: Discarded datum OSGB_1936 in Proj4 definition
#> Transforming to CRS +proj=tmerc +lat_0=49 +lon_0=-2 +k=0.9996012717 +x_0=400000 +y_0=-100000 +ellps=airy +units=m +no_defs