首页 > 解决方案 > R传单GeoJSON着色

问题描述

我仍在研究这个 R Leaflet 自我项目来学习,我正在尝试在北卡罗来纳州罗利的维克县地区的一些多边形上着色。下面是我要着色的图像。

https://imgur.com/a/xdvNLvM

基本上我试图让每个多边形的颜色不同。我试过 addPolygons 但我想我没有正确的多边形数据。我看过颜色分级,但我似乎没有想法。下面是我的代码。我什至尝试取消嵌套 GeoJSON 数据并创建一个因子调色板,但这似乎不起作用。

library(shiny)
library(leaflet.extras)
library(geojsonio)
library(rgdal)

dataurl <- 'https://opendata.arcgis.com/datasets/f5c3b84a6fcc499d8f9ece78602258eb_0.geojson'
data <- geojson_read(dataurl, method = 'web', parse = FALSE, what = 'list')

wake <- readOGR(dataurl)
wake$zips <- factor(sample.int(39L, nrow(wake), TRUE))

#bikedata <- 'D:/bicycle-crash-data-chapel-hill-region.geojson'
#bike <- geojson_read(bikedata)

vtdata <- 'http://geodata.vermont.gov/datasets/4c206846699947429df59c8cb552ab5c_11.geojson'
vt <- geojson_read(vtdata)

factpal <- colorFactor(topo.colors(39), wake$zips)

ui <- shinyUI(
    fluidPage(
        leafletOutput("map", width = "100%", height = "900px")
    )
)

server <- function(input, output) {

    wakegeojson <- reactive({
        data
    })

    #bikegeojson <- reactive({
    #    bike
    #})

    vtgeojson <- reactive({
        vt
    })

    output$map <- renderLeaflet({
        leaflet() %>%
            addTiles() %>%
            setView(-93.65, 42.0285, zoom = 4)

    })

    observe({
        leafletProxy("map") %>%
            addWMSTiles("http://mesonet.agron.iastate.edu/cgi-bin/wms/nexrad/n0r.cgi",
            layers = "nexrad-n0r-900913",
            options = WMSTileOptions(format = "image/png", transparent = TRUE),
            attribution = "") %>%
            addGeoJSON(wakegeojson(), weight = 3, fill = factpal) %>%
            #addGeoJSON(bikegeojson()) %>%
            addGeoJSON(vtgeojson(), fill = FALSE, color = "black")
    })
}

app <- shinyApp(ui = ui, server = server)
runApp(app, launch.browser = TRUE)

我想我需要更多地探索 addPolygons 功能,但我不完全确定如何做到这一点或如何解析/取消嵌套我的 GeoJSON 数据以完成用不同颜色填充威克县邮政编码。任何帮助总是受到赞赏。谢谢你。

标签: rcolorsleafletpolygons

解决方案


我会切换到sf. 您可以直接加载 geojson 并生成一个 Multipolygon 和一个 Multilinestring 对象,它们的读取速度也比readOGR.

然后你可以把这些对象放在addPolygonsand中addPolylines

以下示例应该有效:

library(shiny)
library(leaflet.extras)
library(geojsonio)
library(rgdal)
library(sf)

dataurl <- 'https://opendata.arcgis.com/datasets/f5c3b84a6fcc499d8f9ece78602258eb_0.geojson'
wake <- st_read(dataurl)
wake$zips <- factor(sample.int(39L, nrow(wake), TRUE))

vtdata <- 'http://geodata.vermont.gov/datasets/4c206846699947429df59c8cb552ab5c_11.geojson'
vt <- st_read(vtdata)

factpal <- colorFactor(topo.colors(39), wake$zips)

ui <- shinyUI(
  fluidPage(
    leafletOutput("map", width = "100%", height = "900px")
  )
)

server <- function(input, output) {

  wakegeojson <- reactive({
    wake
  })

  vtgeojson <- reactive({
    vt
  })

  output$map <- renderLeaflet({
    leaflet() %>%
      addTiles() %>%
      addPolygons(data=wakegeojson(), color=factpal(wake$zips)) %>%
      addPolylines(data=vtgeojson(), color="red")
  })
}

app <- shinyApp(ui = ui, server = server)
runApp(app, launch.browser = TRUE)

推荐阅读