mass_gis {macleish}

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

Source

https://www.mass.gov/info-details/massgis-data-layers

Examples

## 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")


## End(Not run)