I wanted to use the NLCD 2016 raster data to calculate land cover is a drainage area. This script can also be used to calculate land cover in any boundary polygon based on the NLCD data. Let’s pull in data from north Jersey that I downloaded from: https://www.mrlc.gov/viewer/

I tried to get most of north Jersey for my selection on the MRLC viewer because this is the area I’m interested in for analysis of wild trout waters in NJ.

NLCD16 <- raster("NLCD2016/NLCD_2016_Land_Cover_L48_20190424_e1OPyurN7yAxizi2xLYK.tiff")
plot(NLCD16, axes = FALSE)

This is North Jersey for the most part. Now we need to determine what the land cover categories are. I’d like to limit it to the following (values associated with Land Cover Class Code Value from NLCD 2016):

This will make it easier to analyze and comparable to some other common Land Cover Types typically used during analyses.

vals <- unique(values(NLCD16))
recl<-matrix(c(vals, c(41,21,81,41,41,90,21,81,90,11,31,21,41,21,81, 0)),ncol=2)   ## Here we're reclassifying some of the land use types to simplify (i.e. land type 41 is representing all forest types (41, 42,and 43))
landuse <- reclassify(NLCD16, rcl=recl)
plot(landuse, legend = FALSE, axes = FALSE)

The map has less variety of color because we decreased the categories to make analysis a bit simpler.

We will read in one of our catchments that was suffering from stormwater impacts. We’ll also need to adjust the projection of the catchment to match that of the NLCD 2016 raster.

Catchment <- st_read("Catchments/GILLHILL.shp")
## Reading layer `GILLHILL' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\GILLHILL.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.53354 ymin: 40.77275 xmax: -74.50772 ymax: 40.79447
## geographic CRS: WGS 84
crs(Catchment) ## NJ State Plane, change this to WGS84
## CRS arguments: +proj=longlat +datum=WGS84 +no_defs
Catchment <- st_transform(Catchment, crs = crs(landuse))
ggplot() +
  geom_sf(data = Catchment) +
  ggtitle("Gillespie Hill", subtitle = "Catchment area upstream of temperature logger")

Looks like an odd shape but that’s our drainage area/catchment.

Clipping our raster

Next we need to clip our raster (using the crop function) to the extent of our catchment (or whatever boundary polygon we’d like to use), create a raster version of that area (using the rasterize function) and finally create a poly-raster using mask.

cr<- crop(landuse, Catchment)
fr<- rasterize(Catchment, cr)
lr<- mask(x=cr,mask=fr)

# use gplot from rasterVis
# geom_tile adds the raster, geom_polygon adds the region
# geom_text adds the labels at the centroids

gplot(lr)+
  geom_tile(aes(fill=factor(value, labels = c("Water", "Urban", "Barren Land", "Forest", "Agriculture", "Wetlands"))), alpha = 0.8)  +
   scale_fill_manual(values = c("blue", "red", "brown1", "forestgreen", "yellow", "steelblue1"), name = "Land Use Code") +
  coord_equal()

The land cover map above should be an accurate depiction of land cover in this area. It looks very pixelated, but keep in mind that the drainage area we’re looking at is very small, it’s headwater tributary to the Whippany River that contains Brook Trout.

Extracting Results for our Catchment

The next step is about quantifying proportion of land cover types for our catchment of interest.

The above results are for an individual catchment results. What if we had multiple catchments or drainages we were interested in comparing? Below shows a loop function to quickly run the above steps over multiple areas and create a dataframe with the results of Land Use calculations.

Loop function

## Reading layer `ANDRSN' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\ANDRSN.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.93199 ymin: 40.76715 xmax: -74.90935 ymax: 40.782
## geographic CRS: WGS 84
## Reading layer `CHNGWTR' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\CHNGWTR.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -75.00959 ymin: 40.72482 xmax: -74.92998 ymax: 40.77318
## geographic CRS: WGS 84
## Reading layer `ELETRC1' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\ELETRC1.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.80771 ymin: 40.80074 xmax: -74.77733 ymax: 40.82959
## geographic CRS: WGS 84
## Reading layer `EMMONS1' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\EMMONS1.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.80538 ymin: 41.0313 xmax: -74.77987 ymax: 41.05705
## geographic CRS: WGS 84
## Reading layer `FLANDER' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\FLANDER.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.71228 ymin: 40.85002 xmax: -74.67377 ymax: 40.88412
## geographic CRS: WGS 84
## Reading layer `FRANK' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\FRANK.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -75.01836 ymin: 40.70537 xmax: -74.98354 ymax: 40.73302
## geographic CRS: WGS 84
## Reading layer `GILLHILL' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\GILLHILL.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.53354 ymin: 40.77275 xmax: -74.50772 ymax: 40.79447
## geographic CRS: WGS 84
## Reading layer `IRONIA' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\IRONIA.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.62604 ymin: 40.83132 xmax: -74.60894 ymax: 40.84291
## geographic CRS: WGS 84
## Reading layer `KREUG' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\KREUG.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.71932 ymin: 40.84138 xmax: -74.70333 ymax: 40.85617
## geographic CRS: WGS 84
## Reading layer `LAMINGT' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\LAMINGT.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.61687 ymin: 40.85813 xmax: -74.60612 ymax: 40.87461
## geographic CRS: WGS 84
## Reading layer `LEDGEB' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\LEDGEB.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.69007 ymin: 40.87192 xmax: -74.64757 ymax: 40.90556
## geographic CRS: WGS 84
## Reading layer `MAPLE' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\MAPLE.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.39136 ymin: 40.98809 xmax: -74.36524 ymax: 41.01145
## geographic CRS: WGS 84
## Reading layer `MC2' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\MC2.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -75.11902 ymin: 40.70973 xmax: -75.06776 ymax: 40.76786
## geographic CRS: WGS 84
## Reading layer `MTOLIVE3' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\MTOLIVE3.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.8234 ymin: 40.8215 xmax: -74.76584 ymax: 40.87803
## geographic CRS: WGS 84
## Reading layer `NDRAKE' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\NDRAKE.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.7804 ymin: 40.8511 xmax: -74.76782 ymax: 40.86063
## geographic CRS: WGS 84
## Reading layer `NORTON' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\NORTON.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.9884 ymin: 40.66492 xmax: -74.95964 ymax: 40.696
## geographic CRS: WGS 84
## Reading layer `OAKDLE' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\OAKDLE.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.69939 ymin: 40.78539 xmax: -74.66605 ymax: 40.802
## geographic CRS: WGS 84
## Reading layer `PCWG' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\PCWG.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -75.08309 ymin: 40.68219 xmax: -75.05552 ymax: 40.71526
## geographic CRS: WGS 84
## Reading layer `RAVINE' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\RAVINE.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.65609 ymin: 40.71906 xmax: -74.64366 ymax: 40.73889
## geographic CRS: WGS 84
## Reading layer `Rinehart' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\Rinehart.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 2 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.76688 ymin: 40.73924 xmax: -74.73384 ymax: 40.7624
## geographic CRS: WGS 84
## Reading layer `RINEHART2' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\RINEHART2.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.76688 ymin: 40.73924 xmax: -74.73703 ymax: 40.7624
## geographic CRS: WGS 84
## Reading layer `RINEHART3' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\RINEHART3.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.76688 ymin: 40.73924 xmax: -74.73798 ymax: 40.7624
## geographic CRS: WGS 84
## Reading layer `RINEHART5' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\RINEHART5.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.76688 ymin: 40.73924 xmax: -74.74652 ymax: 40.7624
## geographic CRS: WGS 84
## Reading layer `RINEHART6' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\RINEHART6.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.76628 ymin: 40.74853 xmax: -74.75147 ymax: 40.7624
## geographic CRS: WGS 84
## Reading layer `SASBURY' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\SASBURY.shp' using driver `ESRI Shapefile'
## Simple feature collection with 2 features and 13 fields
## geometry type:  MULTIPOLYGON
## dimension:      XY
## bbox:           xmin: -75.01491 ymin: 40.66672 xmax: -74.99136 ymax: 40.69258
## geographic CRS: WGS 84
## Reading layer `SR1' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\SR1.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -75.13207 ymin: 40.61584 xmax: -75.10192 ymax: 40.63189
## geographic CRS: WGS 84
## Reading layer `StoneH' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\StoneH.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.43079 ymin: 40.96514 xmax: -74.33266 ymax: 41.01065
## geographic CRS: WGS 84
## Reading layer `SUNV2' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\SUNV2.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.76094 ymin: 40.84715 xmax: -74.72163 ymax: 40.86449
## geographic CRS: WGS 84
## Reading layer `SWBudd' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\SWBudd.shp' using driver `ESRI Shapefile'
## Simple feature collection with 2 features and 127 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: 519010 ymin: 4522810 xmax: 519880 ymax: 4524260
## projected CRS:  NAD83 / UTM zone 18N
## Reading layer `THB1' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\THB1.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -75.0686 ymin: 40.63923 xmax: -75.03916 ymax: 40.66095
## geographic CRS: WGS 84
## Reading layer `TroutB' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\TroutB.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.75558 ymin: 40.74795 xmax: -74.72848 ymax: 40.77598
## geographic CRS: WGS 84
## Reading layer `WHITE' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\WHITE.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.80884 ymin: 41.27056 xmax: -74.77773 ymax: 41.30126
## geographic CRS: WGS 84
## Reading layer `WPC1' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\WPC1.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -75.05878 ymin: 40.63932 xmax: -75.0278 ymax: 40.65095
## geographic CRS: WGS 84
## Reading layer `WPC2' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\WPC2.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -75.05878 ymin: 40.63932 xmax: -75.01538 ymax: 40.67005
## geographic CRS: WGS 84

If you want to calculate impervious surface you’ll need to use a different raster that can be found from the same website as the NLCD layer (NLCD 2016 Impervious). Below calculates pervious vs impervious for each polygon run through the loop and the output is a dataframe.

## Reading layer `ANDRSN' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\ANDRSN.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.93199 ymin: 40.76715 xmax: -74.90935 ymax: 40.782
## geographic CRS: WGS 84
## Reading layer `CHNGWTR' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\CHNGWTR.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -75.00959 ymin: 40.72482 xmax: -74.92998 ymax: 40.77318
## geographic CRS: WGS 84
## Reading layer `ELETRC1' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\ELETRC1.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.80771 ymin: 40.80074 xmax: -74.77733 ymax: 40.82959
## geographic CRS: WGS 84
## Reading layer `EMMONS1' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\EMMONS1.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.80538 ymin: 41.0313 xmax: -74.77987 ymax: 41.05705
## geographic CRS: WGS 84
## Reading layer `FLANDER' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\FLANDER.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.71228 ymin: 40.85002 xmax: -74.67377 ymax: 40.88412
## geographic CRS: WGS 84
## Reading layer `FRANK' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\FRANK.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -75.01836 ymin: 40.70537 xmax: -74.98354 ymax: 40.73302
## geographic CRS: WGS 84
## Reading layer `GILLHILL' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\GILLHILL.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.53354 ymin: 40.77275 xmax: -74.50772 ymax: 40.79447
## geographic CRS: WGS 84
## Reading layer `IRONIA' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\IRONIA.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.62604 ymin: 40.83132 xmax: -74.60894 ymax: 40.84291
## geographic CRS: WGS 84
## Reading layer `KREUG' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\KREUG.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.71932 ymin: 40.84138 xmax: -74.70333 ymax: 40.85617
## geographic CRS: WGS 84
## Reading layer `LAMINGT' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\LAMINGT.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.61687 ymin: 40.85813 xmax: -74.60612 ymax: 40.87461
## geographic CRS: WGS 84
## Reading layer `LEDGEB' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\LEDGEB.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.69007 ymin: 40.87192 xmax: -74.64757 ymax: 40.90556
## geographic CRS: WGS 84
## Reading layer `MAPLE' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\MAPLE.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.39136 ymin: 40.98809 xmax: -74.36524 ymax: 41.01145
## geographic CRS: WGS 84
## Reading layer `MC2' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\MC2.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -75.11902 ymin: 40.70973 xmax: -75.06776 ymax: 40.76786
## geographic CRS: WGS 84
## Reading layer `MTOLIVE3' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\MTOLIVE3.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.8234 ymin: 40.8215 xmax: -74.76584 ymax: 40.87803
## geographic CRS: WGS 84
## Reading layer `NDRAKE' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\NDRAKE.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.7804 ymin: 40.8511 xmax: -74.76782 ymax: 40.86063
## geographic CRS: WGS 84
## Reading layer `NORTON' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\NORTON.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.9884 ymin: 40.66492 xmax: -74.95964 ymax: 40.696
## geographic CRS: WGS 84
## Reading layer `OAKDLE' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\OAKDLE.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.69939 ymin: 40.78539 xmax: -74.66605 ymax: 40.802
## geographic CRS: WGS 84
## Reading layer `PCWG' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\PCWG.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -75.08309 ymin: 40.68219 xmax: -75.05552 ymax: 40.71526
## geographic CRS: WGS 84
## Reading layer `RAVINE' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\RAVINE.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.65609 ymin: 40.71906 xmax: -74.64366 ymax: 40.73889
## geographic CRS: WGS 84
## Reading layer `Rinehart' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\Rinehart.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 2 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.76688 ymin: 40.73924 xmax: -74.73384 ymax: 40.7624
## geographic CRS: WGS 84
## Reading layer `RINEHART2' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\RINEHART2.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.76688 ymin: 40.73924 xmax: -74.73703 ymax: 40.7624
## geographic CRS: WGS 84
## Reading layer `RINEHART3' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\RINEHART3.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.76688 ymin: 40.73924 xmax: -74.73798 ymax: 40.7624
## geographic CRS: WGS 84
## Reading layer `RINEHART5' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\RINEHART5.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.76688 ymin: 40.73924 xmax: -74.74652 ymax: 40.7624
## geographic CRS: WGS 84
## Reading layer `RINEHART6' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\RINEHART6.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.76628 ymin: 40.74853 xmax: -74.75147 ymax: 40.7624
## geographic CRS: WGS 84
## Reading layer `SASBURY' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\SASBURY.shp' using driver `ESRI Shapefile'
## Simple feature collection with 2 features and 13 fields
## geometry type:  MULTIPOLYGON
## dimension:      XY
## bbox:           xmin: -75.01491 ymin: 40.66672 xmax: -74.99136 ymax: 40.69258
## geographic CRS: WGS 84
## Reading layer `SR1' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\SR1.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -75.13207 ymin: 40.61584 xmax: -75.10192 ymax: 40.63189
## geographic CRS: WGS 84
## Reading layer `StoneH' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\StoneH.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.43079 ymin: 40.96514 xmax: -74.33266 ymax: 41.01065
## geographic CRS: WGS 84
## Reading layer `SUNV2' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\SUNV2.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.76094 ymin: 40.84715 xmax: -74.72163 ymax: 40.86449
## geographic CRS: WGS 84
## Reading layer `SWBudd' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\SWBudd.shp' using driver `ESRI Shapefile'
## Simple feature collection with 2 features and 127 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: 519010 ymin: 4522810 xmax: 519880 ymax: 4524260
## projected CRS:  NAD83 / UTM zone 18N
## Reading layer `THB1' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\THB1.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -75.0686 ymin: 40.63923 xmax: -75.03916 ymax: 40.66095
## geographic CRS: WGS 84
## Reading layer `TroutB' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\TroutB.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.75558 ymin: 40.74795 xmax: -74.72848 ymax: 40.77598
## geographic CRS: WGS 84
## Reading layer `WHITE' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\WHITE.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -74.80884 ymin: 41.27056 xmax: -74.77773 ymax: 41.30126
## geographic CRS: WGS 84
## Reading layer `WPC1' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\WPC1.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -75.05878 ymin: 40.63932 xmax: -75.0278 ymax: 40.65095
## geographic CRS: WGS 84
## Reading layer `WPC2' from data source `C:\Scott20171018\MainFiles\R\LandUseAnalysis\Catchments\WPC2.shp' using driver `ESRI Shapefile'
## Simple feature collection with 1 feature and 11 fields
## geometry type:  POLYGON
## dimension:      XY
## bbox:           xmin: -75.05878 ymin: 40.63932 xmax: -75.01538 ymax: 40.67005
## geographic CRS: WGS 84