## ----message=FALSE------------------------------------------------------------
library(wikilake)
library(nhdR)

## ----message=FALSE, eval=FALSE,cache=FALSE------------------------------------
#  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)

## ----message=FALSE, fig.width=6, fig.height=5, fig.align='center'-------------
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)