plotlymapviewex

gistfile1.txt

# load packages 
library(sf)
library(plotly)
library(mapview)

# download data and save in memory
wsa <- read.csv("https://www.epa.gov/sites/production/files/2014-10/wsa_siteinfo_ts_final.csv")

# make an sf object
wsa <- st_as_sf(wsa, coords = c("LON_DD", "LAT_DD"), crs = 4269,agr = "constant")

# plotly scatterplot
plot_ly(wsa, x = ~XELEV, y = ~WSAREA) %>% 
  add_markers()

# plotly using ggplotly
p <- ggplot(wsa, aes(x = XELEV, y = WSAREA)) +
  geom_point() 
ggplotly(p)

# use mapview
mapview(wsa, zcol = 'XELEV')
mapview(wsa, zcol = 'WSAREA')
mapview(wsa, zcol = 'ECOWSA9')