Creating simple maps using nhdR

Joseph Stachelek

2021-01-28

library(wikilake)
library(nhdR)

gull_meta  <- wikilake::lake_wiki("Gull Lake (Michigan)")
gull <- nhd_plus_query(gull_meta$Lon, gull_meta$Lat,
              dsn = c("NHDWaterbody", "NHDFlowLine"), buffer_dist = 0.05)

plot(gull$sp$NHDWaterbody$geometry, col = "blue", main = "Gull Lake (Michigan)")
plot(gull$sp$NHDFlowLine$geometry, col = "cyan", add = TRUE)
plot(gull$pnt, col = "red", pch = 19, add = TRUE)
axis(1); axis(2)