stplanr (version 0.1.4)

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

Description

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

Usage

reproject(shp, crs = crs_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 crs_select_aeq).

Examples

rf_aeq = reproject(routes_fast)
rf_osgb = reproject(routes_fast, 27700)
cor(rgeos::gLength(rf_aeq, byid = TRUE), rgeos::gLength(rf_osgb, byid = TRUE))
rf_aeq_wgs84 = sp::spTransform(rf_aeq, CRS("+init=epsg:4326"))
rf_osgb_wgs84 = sp::spTransform(rf_osgb, CRS("+init=epsg:4326"))
plot(rf_aeq_wgs84)
plot(rf_osgb_wgs84, col = "red", add = TRUE)