Как создать листовку, картографическую карту округа США - PullRequest
Новая
       33

Как создать листовку, картографическую карту округа США

0 голосов
/ 19 февраля 2020

Используя приведенный ниже код, я получаю свой фрейм данных с данными округа США 

library(raster)
library(leaflet)
library(tidyverse)

# Get USA polygon data
USA <- getData("GADM", country = "usa", level = 2)

### Get data
mydata <- read.csv("https://www.betydb.org/miscanthus_county_avg_yield.csv",
                   stringsAsFactors = FALSE)

Моя задача состоит в том, чтобы создать интерактивную картографическую карту листовки размером Avg_yield, поэтому сначала я укрепляю свои данные полигона США 

library(rgeos)
library(maptools)
library(ggplot2)
states.shp.f <- fortify(USA, region = "NAME_2")

Затем я поднабор своего набора данных и объединение его с укрепленным: 

mydata2<-mydata[,c("COUNTY_NAME","Avg_yield")]
colnames(mydata2)[1]<-"id"

## merge shape file with data
merge.shp.coef <- merge(states.shp.f, mydata2, by = "id")

, но теперь у меня есть набор данных с именем каждого округа много раз, а также некоторые округа имеют различные значения Avg_yield. Как правильно обрабатывать эти данные, чтобы использовать код листовки, например: 

leaflet() %>% 
  addProviderTiles("OpenStreetMap.Mapnik") %>%
  setView(lat = 39.8283, lng = -98.5795, zoom = 4) %>%
  addPolygons(data = USA, stroke = FALSE, smoothFactor = 0.2, fillOpacity = 0.3,
              fillColor = ~mypal(mydata$Avg_yield),
              popup = paste("Region: ", USA$NAME_2, "<br>",
                            "Avg_yield: ", mydata$Avg_yield, "<br>")) %>%
  addLegend(position = "bottomleft", pal = mypal, values = mydata$Avg_yield,
            title = "Avg_yield",
            opacity = 1)

1 Ответ

1 голос
/ 20 февраля 2020

Предлагаемый способ сделать это - преобразовать ваш объект многоугольника в объект sf с помощью st_as_sf ()

Вот вам рабочий пример: (Я использовал некоторые другие данные для многоугольника, я думал, что ваши слишком точный и требует много ресурсов, плюс я заставил его работать с блестящими) 

library(leaflet)
library(tidyverse)
library(ggplot2)
library(sf)
library(shiny)

USA <- st_read(dsn = '[your path]/cb_2018_us_county_500k.shp')

### Get data
mydata <- read.csv("https://www.betydb.org/miscanthus_county_avg_yield.csv",
                   stringsAsFactors = FALSE)

states_sf <- st_as_sf(USA)

mydata2<-mydata[,c("COUNTY_NAME","Avg_yield")]
colnames(mydata2)[1]<-"NAME"

## merge shape file with data
states_sf_coef <- left_join(states_sf, mydata2, by = "NAME")


ui <- fluidPage(

  leafletOutput("map", height = "100vh")
)

server <- function(input, output, session) {

  bins <- c(0, 5, 10, 15, 20, 25, 30, 35, 40)
  mypal <- colorBin("YlOrRd", domain = states_sf_coef$Avg_yield, bins = bins)

  #Sortie map
  output$map <- renderLeaflet({
    leaflet()%>%
      addProviderTiles("OpenStreetMap.Mapnik")%>%
      setView(lat = 39.8283, lng = -98.5795, zoom = 4) %>%
      addPolygons(
        data = states_sf_coef,
        fillColor = ~mypal(states_sf_coef$Avg_yield),
        stroke = FALSE,
        smoothFactor = 0.2,
        fillOpacity = 0.3,
        popup = paste("Region: ", states_sf_coef$NAME_2, "<br>",
                      "Avg_yield: ", states_sf_coef$Avg_yield, "<br>"))%>%
      addLegend(position = "bottomleft",
                pal = mypal,
                values = states_sf_coef$Avg_yield,
                title = "Avg_yield",
                opacity = 1)
  })

}

# Create Shiny app ----
shinyApp(ui = ui, server = server)
...