閃亮：如何在同一張小冊子地圖上繪製幾個參數？

Question

我是新來閃亮。我想製作一個閃亮的應用程序，顯示不同參數的空間分佈。我使用sp包中的meuse數據集。

下面是我用

library(sp) 
library(rgdal) 
library(ggmap) 
library(leaflet) 
library(dplyr) 
library(shiny) 

ui <- fluidPage(
    mainPanel(
    titlePanel("Copper (ppm)"), 
    leafletOutput("copper"), 
    titlePanel("Lead (ppm)"), 
    leafletOutput("lead"), 
    titlePanel("Zinc (ppm)"), 
    leafletOutput("zinc") 
    ) 
) 
server <- function(input,output){ 
    output$copper <- renderLeaflet({ 

data(meuse) 
coordinates(meuse) <- 1:2 
proj4string(meuse) <- CRS("+init=epsg:28992")     
meuse <- spTransform(meuse, CRS("+proj=longlat +datum=WGS84")) 
meuse_df <- as.data.frame(meuse) 

leaflet(meuse_df) %>% 
    addProviderTiles("OpenStreetMap", group = "OpenStreetMap") %>% 
    addProviderTiles("Esri.WorldImagery", group = "Esri.WorldImagery") %>% 
    addCircleMarkers(~x, 
        ~y, 
        radius = ~copper/10, 
        stroke = FALSE, 
        fillOpacity = 0.4, 
        group = "meuse_df", 
        popup = ~as.character(copper)) %>% 
addLayersControl(position = "bottomleft", 
       baseGroups = c("OpenStreetMap", 
          "Esri.WorldImagery"), 
overlayGroups = "meuse_df") 
}) 
    output$lead <- renderLeaflet({ 

    leaflet(meuse_df) %>% 
    addProviderTiles("OpenStreetMap", group = "OpenStreetMap") %>% 
    addProviderTiles("Esri.WorldImagery", group = "Esri.WorldImagery") %>% 
    addCircleMarkers(~x, 
        ~y, 
        radius = ~lead/50, 
        stroke = FALSE, 
        fillOpacity = 0.4, 
        group = "meuse_df", 
        popup = ~as.character(lead)) %>% 
    addLayersControl(position = "bottomleft", 
        baseGroups = c("OpenStreetMap", 
            "Esri.WorldImagery"), 
        overlayGroups = "meuse_df") 
}) 

output$zinc <- renderLeaflet({ 

leaflet(meuse_df) %>% 
    addProviderTiles("OpenStreetMap", group = "OpenStreetMap") %>% 
    addProviderTiles("Esri.WorldImagery", group = "Esri.WorldImagery") %>% 
    addCircleMarkers(~x, 
        ~y, 
        radius = ~zinc/100, 
        stroke = FALSE, 
        fillOpacity = 0.4, 
        group = "meuse_df", 
        popup = ~as.character(zinc)) %>% 
    addLayersControl(position = "bottomleft", 
        baseGroups = c("OpenStreetMap", 
            "Esri.WorldImagery"), 
        overlayGroups = "meuse_df") 
}) 
} 
shinyApp(ui = ui, server = server) 

代碼，這裏是我

我不知道是否有繪製的所有參數（銅，鉛的方式和鋅的結果）在一張地圖上。任何建議，將不勝感激。

UPDATE

感謝@ Symbolix答案，並用checkBoxGroupInput的建議。相反，我用addCircleMarkers三次，所以我可以繪製在一個地圖上所有的金屬，我可以接通他們，並關閉

ui <- fluidPage(
mainPanel(
titlePanel("All metals (ppm)"), 
leafletOutput("metals") 
) 
) 
server <- function(input,output){ 
output$metals <- renderLeaflet({ 

data(meuse) 
coordinates(meuse) <- 1:2 
proj4string(meuse) <- CRS("+init=epsg:28992")     
meuse <- spTransform(meuse, CRS("+proj=longlat +datum=WGS84")) 
meuse_df <- as.data.frame(meuse) 

leaflet(meuse_df) %>% 
    addProviderTiles("OpenStreetMap", group = "OpenStreetMap") %>% 
    addProviderTiles("Esri.WorldImagery", group = "Esri.WorldImagery") %>% 
    addCircleMarkers(~x, 
        ~y, 
        radius = ~copper/10, 
        color ="red", 
        stroke = FALSE, 
        fillOpacity = 0.4, 
        group = "copper", 
        popup = ~as.character(copper)) %>% 
    addCircleMarkers(~x, 
        ~y, 
        radius = ~lead/50, 
        color ="gren", 
        stroke = FALSE, 
        fillOpacity = 0.4, 
        group = "lead", 
        popup = ~as.character(lead)) %>% 
    addCircleMarkers(~x, 
        ~y, 
        radius = ~zinc/100, 
        color ="blue", 
        stroke = FALSE, 
        fillOpacity = 0.4, 
        group = "zinc", 
        popup = ~as.character(zinc)) %>% 
    addLayersControl(position = "bottomleft", 
        baseGroups = c("OpenStreetMap", 
            "Esri.WorldImagery"), 
        overlayGroups = c("copper", 
            "lead", 
            "zinc")) 
}) 
} 
shinyApp(ui = ui, server = server) 

Answer 1

融化你的資料，以便您的金屬在一列，你是好去。

這裏我用library(reshape2)來做熔化。

library(sp) 
library(rgdal) 
library(ggmap) 
library(leaflet) 
library(dplyr) 
library(shiny) 
library(reshape2) 

ui <- fluidPage(
    mainPanel(
     titlePanel("Metals"), 
     leafletOutput("all_metals") 
    ) 
) 
server <- function(input,output){ 

    output$all_metals <- renderLeaflet({ 

     data(meuse) 
     coordinates(meuse) <- 1:2 
     proj4string(meuse) <- CRS("+init=epsg:28992")     
     meuse <- spTransform(meuse, CRS("+proj=longlat +datum=WGS84")) 
     meuse_df <- as.data.frame(meuse) 

     ## melt df so 'metals' are in one column 
     ## using 'reshape2' library 
     meuse_melt <- melt(meuse_df, measure.vars = c("copper","lead","zinc"), variable.name = "metal") 

     ## specify factor levels for colours 
     meuse_melt$metal <- factor(sample.int(5L, nrow(meuse_melt), TRUE)) 
     factpal <- colorFactor(topo.colors(5), meuse_melt$metal) 

     ## now you just need one output 
     leaflet(meuse_melt) %>% 
      addProviderTiles("OpenStreetMap", group = "OpenStreetMap") %>% 
      addProviderTiles("Esri.WorldImagery", group = "Esri.WorldImagery") %>% 
      addCircleMarkers(~x, 
          ~y, 
          radius = ~value/100, 
          stroke = FALSE, 
          fillOpacity = 0.4, 
          group = "meuse_melt", 
          popup = ~metal, 
          color= ~factpal(metal)) %>% 
      addLayersControl(position = "bottomleft", 
          baseGroups = c("OpenStreetMap", "Esri.WorldImagery"), overlayGroups = "meuse_melt") 
    }) 

} 
shinyApp(ui = ui, server = server) 

Answer 2

爲了節省您的代碼行，你可以使用MapView的它提供了多層次的地圖開箱。如果你只想要某些屬性，你可以簡單地將他們的名字（或列號）提供給'zcol'參數。或者，您可以使用'burst = TRUE'來顯示屬性表中的所有圖層/列。

library(mapview) 
library(sp) 

data(meuse) 
coordinates(meuse) <- ~x+y 
proj4string(meuse) <- CRS("+init=epsg:28992")     
mapview(meuse, zcol = c("copper", "lead", "zinc")) 

## all layers 
mapview(meuse, burst = TRUE) 

mapview可以和renderLeaflet一起使用。