Guest User

Untitled

a guest
Feb 25th, 2018
68
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.23 KB | None | 0 0
  1. County State Value_x
  2. .
  3. Adair Oklahoma 5
  4. Adair Missouri 2
  5. Adair Kentucky 10
  6. Adair lowas 1
  7. .
  8.  
  9. library(raster)
  10. library(leaflet)
  11. library(tidyverse)
  12.  
  13. # Get USA polygon data
  14. USA <- getData("GADM", country = "usa", level = 2)
  15.  
  16. ### Get data. I can't provide the date set. It has county, state and value _x as shown above
  17. mydata <- county,state,value_x
  18.  
  19. ### Check counties that exist in USA, but not in mydata
  20. ### Create a dummy data frame and bind it with mydata
  21.  
  22. mydata <- data.frame(COUNTY_NAME = setdiff(USA$NAME_2, mydata$county),
  23. value_x = NA,
  24. stringsAsFactors = FALSE) %>%
  25. bind_rows(mydata)
  26.  
  27. ### Create a color palette
  28. mypal <- colorNumeric(palette = "viridis", domain = mydata$value_x)
  29.  
  30. leaflet() %>%
  31. addProviderTiles("OpenStreetMap.Mapnik") %>%
  32. setView(lat = 39.8283, lng = -98.5795, zoom = 4) %>%
  33. addPolygons(data = USA, stroke = FALSE, smoothFactor = 0.2, fillOpacity = 0.3,
  34. fillColor = ~mypal(mydata$value_x),
  35. popup = paste("Region: ", USA$NAME_2, "<br>",
  36. "Value X: ", mydata$value_x, "<br>")) %>%
  37. addLegend(position = "bottomleft", pal = mypal, values = mydata$value_x,
  38. title = "Value X",
  39. opacity = 1)
Add Comment
Please, Sign In to add comment