macleish (version 0.3.8)

mass_gis: Retrieve elevation layers from MassGIS

Description

Retrieve elevation layers from MassGIS

Usage

mass_gis(layer = "contours250k")
macleish_intersect(x)

Arguments

layer

MassGIS layer name to import

x

an sf::sf() object

Details

This function will download shapefiles from MassGIS, unzip them, transform the projection to EPSG:4326, compute their intersection with the boundary of the MacLeish property, and return the resulting sf::sf() object.

Intersect a spatial layer with the MacLeish boundary layer

Examples

# NOT RUN {
# }
# NOT RUN {
# have to download the shapefiles...could take a while...
elevation <- mass_gis()
macleish_elevation <- macleish_intersect(elevation)
if (require(sf)) {
  plot(macleish_elevation)
}

dcr_trails <- mass_gis("dcrtrails")

# }