mass_gis: Retrieve elevation layers from MassGIS

In macleish: Retrieve Data from MacLeish Field Station

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)

macleish documentation built on July 7, 2022, 1:05 a.m.