Photo by Ricardo Resende on Unsplash

This homework focuses on analysis of habitat classification using a chunk of landscape including the Rice River Center. The data were grabbed from the NLCD viewer website consisting of a (roughtly) 25 km\(^2\) chunk of landscape.

Here are the raw data url’s on GitHub.

library( raster )
library( landscapemetrics )

land_2011 <- "https://github.com/dyerlab/ENVS-Lectures/raw/master/data/NLCD_2011_Land_Cover_L48_20190424_qn2B1f8ganicJNKnJN0e.tiff"
land_2013 <- "https://github.com/dyerlab/ENVS-Lectures/raw/master/data/NLCD_2013_Land_Cover_L48_20190424_qn2B1f8ganicJNKnJN0e.tiff"
land_2016 <- "https://github.com/dyerlab/ENVS-Lectures/raw/master/data/NLCD_2016_Land_Cover_L48_20190424_qn2B1f8ganicJNKnJN0e.tiff"
land_legend <- "https://raw.githubusercontent.com/dyerlab/ENVS-Lectures/master/data/NLCD_landcover_legend_2018_12_17_qn2B1f8ganicJNKnJN0e.csv"

Questions

  1. Landscape structure can be examined by type, composition, and configuration. Load in the above three time periods for the same location and compare them in terms of Type, Composition, and Configuration (where configuration is measured by average patch area for forest types (Value %in% c(41, 42, 43) )).

  2. In 2011, the biggest difference in habitat areas was accounted for by Herbaceous ground cover. What happened to that landtype in 2013 and 2016?

LS0tCnRpdGxlOiAiSGFiaXRhdCBIb21ld29yayIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQoKIVtQaG90byBieSBSaWNhcmRvIFJlc2VuZGUgb24gVW5zcGxhc2hdKGh0dHBzOi8vbGl2ZS5zdGF0aWNmbGlja3IuY29tLzY1NTM1LzUxMDgzODI4NTMxX2ExNmU4ZDJmMmZfY19kLmpwZykKCgpUaGlzIGhvbWV3b3JrIGZvY3VzZXMgb24gYW5hbHlzaXMgb2YgaGFiaXRhdCBjbGFzc2lmaWNhdGlvbiB1c2luZyBhIGNodW5rIG9mIGxhbmRzY2FwZSBpbmNsdWRpbmcgdGhlIFtSaWNlIFJpdmVyIENlbnRlcl0oaHR0cHM6Ly9yaWNlcml2ZXJzLnZjdS5lZHUpLiAgVGhlIGRhdGEgd2VyZSBncmFiYmVkIGZyb20gdGhlIFtOTENEXShodHRwczovL25sY2QuZ292KSB2aWV3ZXIgd2Vic2l0ZSBjb25zaXN0aW5nIG9mIGEgKHJvdWdodGx5KSAyNSBrbSReMiQgY2h1bmsgb2YgbGFuZHNjYXBlLgoKSGVyZSBhcmUgdGhlIHJhdyBkYXRhIHVybCdzIG9uIEdpdEh1Yi4KCmBgYHtyfQpsaWJyYXJ5KCByYXN0ZXIgKQpsaWJyYXJ5KCBsYW5kc2NhcGVtZXRyaWNzICkKCmxhbmRfMjAxMSA8LSAiaHR0cHM6Ly9naXRodWIuY29tL2R5ZXJsYWIvRU5WUy1MZWN0dXJlcy9yYXcvbWFzdGVyL2RhdGEvTkxDRF8yMDExX0xhbmRfQ292ZXJfTDQ4XzIwMTkwNDI0X3FuMkIxZjhnYW5pY0pOS25KTjBlLnRpZmYiCmxhbmRfMjAxMyA8LSAiaHR0cHM6Ly9naXRodWIuY29tL2R5ZXJsYWIvRU5WUy1MZWN0dXJlcy9yYXcvbWFzdGVyL2RhdGEvTkxDRF8yMDEzX0xhbmRfQ292ZXJfTDQ4XzIwMTkwNDI0X3FuMkIxZjhnYW5pY0pOS25KTjBlLnRpZmYiCmxhbmRfMjAxNiA8LSAiaHR0cHM6Ly9naXRodWIuY29tL2R5ZXJsYWIvRU5WUy1MZWN0dXJlcy9yYXcvbWFzdGVyL2RhdGEvTkxDRF8yMDE2X0xhbmRfQ292ZXJfTDQ4XzIwMTkwNDI0X3FuMkIxZjhnYW5pY0pOS25KTjBlLnRpZmYiCmxhbmRfbGVnZW5kIDwtICJodHRwczovL3Jhdy5naXRodWJ1c2VyY29udGVudC5jb20vZHllcmxhYi9FTlZTLUxlY3R1cmVzL21hc3Rlci9kYXRhL05MQ0RfbGFuZGNvdmVyX2xlZ2VuZF8yMDE4XzEyXzE3X3FuMkIxZjhnYW5pY0pOS25KTjBlLmNzdiIKCmBgYAoKCgojIyBRdWVzdGlvbnMKCjEuIExhbmRzY2FwZSBzdHJ1Y3R1cmUgY2FuIGJlIGV4YW1pbmVkIGJ5IHR5cGUsIGNvbXBvc2l0aW9uLCBhbmQgY29uZmlndXJhdGlvbi4gIExvYWQgaW4gdGhlIGFib3ZlIHRocmVlIHRpbWUgcGVyaW9kcyBmb3IgdGhlIHNhbWUgbG9jYXRpb24gYW5kIGNvbXBhcmUgdGhlbSBpbiB0ZXJtcyBvZiBUeXBlLCBDb21wb3NpdGlvbiwgYW5kIENvbmZpZ3VyYXRpb24gKHdoZXJlIGNvbmZpZ3VyYXRpb24gaXMgbWVhc3VyZWQgYnkgYXZlcmFnZSBwYXRjaCBhcmVhIGZvciBmb3Jlc3QgdHlwZXMgKGBWYWx1ZSAlaW4lIGMoNDEsIDQyLCA0MykgKWApLgoKCjIuIEluIDIwMTEsIHRoZSBiaWdnZXN0IGRpZmZlcmVuY2UgaW4gaGFiaXRhdCBhcmVhcyB3YXMgYWNjb3VudGVkIGZvciBieSBIZXJiYWNlb3VzIGdyb3VuZCBjb3Zlci4gIFdoYXQgaGFwcGVuZWQgdG8gdGhhdCBsYW5kdHlwZSBpbiAyMDEzIGFuZCAyMDE2PwoK