Advertisement
Guest User

Untitled

a guest
Jun 24th, 2019
66
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.75 KB | None | 0 0
  1. library(tidyverse)
  2. library(mapview)
  3. library(sf)
  4.  
  5. lat <- rep(34,20) #sample latitude
  6. lon <- seq(-128, -126, length = 20) #sample longitude
  7.  
  8. set.seed(1)
  9. value <- rnorm(20, 0, 6) #20 values: some positive, some negative
  10.  
  11. #make data.frame
  12. df <- data.frame(lat,lon, value)
  13.  
  14. #make spatial object for mapview
  15. df <- st_as_sf(df, coords = c("lon", "lat"), crs = 4326) %>%
  16. mutate(value_abs = abs(value)) #value_abs intended for `cex` argument
  17.  
  18. pal <- mapviewPalette("mapviewSpectralColors") #from mapview doc. example
  19. m <- mapview(df["value"], #hover over value
  20. cex = "value", #circle diameter value
  21. legend = TRUE,
  22. col.regions = pal(100), #closest I came to achieving a red-blue divergent scale
  23. layer.name = "value")
  24. m
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement