In jsta/nlaR: Interface to the National Lakes Assessment

# rmarkdown::render("vignettes/comparison_2007-2012.Rmd", rmarkdown::html_document(pandoc_args = "--self-contained"))

library(nlaR)

library(sf)
library(maps)

par(mfrow = c(1, 2))

# nla_get(2007, use_rappdirs=TRUE)
dt <- nla_load(2007)

dt <- st_as_sf(dt$SampledLakeInformation, coords = c("LON_DD", "LAT_DD"), crs = 4326)

map("usa")
plot(dt['URBAN'], add = TRUE, reset = FALSE)
title("NLA 2007")

# nla_get(2012, use_rappdirs=TRUE)
dt <- nla_load(2012)
dt <- st_as_sf(dt$wide_siteinfo, coords = c("LON_DD83", "LAT_DD83"), crs = 4326)

map("usa")
plot(dt['URBAN'], add = TRUE, reset = FALSE)
title("NLA 2012")

jsta/nlaR documentation built on March 23, 2023, 11:53 p.m.