Reproject lat/long spatial object so that they are in units of 1m

Source: R/crs-funs.R

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(), gclip(), geo_bb_matrix(), geo_bb(), mapshape_available(), mapshape(), 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 CRS 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