# Creating simple maps using nhdR

#### Jemma Stachelek

#### 2023-08-11

    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)