6 Monitoring location map

6.1 Read site locations

sites=dataRetrieval::readNWISsite(siteNumbers=subset(site_list, SiteType %in% c("Elevation gauge","Causeway breach","Lake", "Lake, DBL"))$SiteNumber)
sites=merge(sites,site_list, by.x=c("site_no","station_nm"), by.y=c("SiteNumber","SiteName"))

6.2 Build site map

sites=within(sites, {
    LatitudeMeasure=dec_lat_va
    LongitudeMeasure=dec_long_va
    MonitoringLocationTypeName=SiteType
    MonitoringLocationIdentifier=site_no
    MonitoringLocationName=station_nm
})
site_map=wqTools::buildMap(sites=sites, plot_polys=F, search="") %>% 
        leaflet::fitBounds(lng1=min(sites$LongitudeMeasure), lng2=max(sites$LongitudeMeasure), lat1=min(sites$LatitudeMeasure), lat2=max(sites$LatitudeMeasure)) %>% 
        leaflet::removeMeasure()

Figure 6.1: Site location map.