R Shiny/R Leaflet/ R : 使用R Shiny中的sliderInput动态呈现choropleth地图 2

huangapple
go评论92阅读模式
英文:

R Shiny/R Leaflet/ R : Dynamically render choropleth map with sliderInput in R shiny 2

问题

我正在尝试使用滑块功能在Leaflet地图上使用Shiny来选择并仅显示选定的多边形。以下是您代码的翻译部分：

  1. 我正在尝试使用滑块功能在Leaflet地图上使用Shiny来选择并仅显示选定的多边形。
  3. 以下是我的代码：
  5. library(shiny)
  6. library(leaflet)
  7. library(rgdal)
  8. library(sf)
  9. library(dplyr)
  10. library(RColorBrewer)
  12. # 读取波浪数据
  13. wave_data <- read_sf("~/path/Wave.shp")
  15. # 转换坐标系
  16. wave_data <- st_transform(wave_data, crs = '+proj=longlat +datum=WGS84')
  18. # 处理数据
  19. wave_data2 <- wave_data %>%
  20.   mutate_at(vars(An_mn_P_OD), ~ as.integer(round(.x)))
  22. # 创建Leaflet地图
  23. wave_data_map <- leaflet() %>%
  24.   addProviderTiles(providers$Esri.WorldTopoMap) %>%
  25.   setView(lng = -4.2026458, lat = 56.4906712, zoom = 5)
  27. # 设置深度分级
  28. bins <- c(0, 25, 50, 100, 150, 200, 300, 400, 500, 1000, 3100)
  29. pal <- colorBin("RdYlBu", domain = wave_data2$Ave_Depth, bins = bins)
  31. # 添加多边形
  32. wave_data_map <- wave_data_map %>%
  33.   addPolygons(data = wave_data2, 
  34.               weight = 1,
  35.               smoothFactor = 0.5,
  36.               color = "white", 
  37.               fillOpacity = 0.5, 
  38.               fillColor = pal(wave_data2$Ave_Depth))
  40. # 添加图例
  41. wave_data_map <- wave_data_map %>%
  42.   addLegend(pal = pal, 
  43.             values = wave_data2$Ave_Depth, 
  44.             title = "Average depth",
  45.             labFormat = labelFormat(suffix = "m"),
  46.             opacity = 0.7, 
  47.             position = "bottomright")
  49. # 定义UI
  50. ui <- bootstrapPage(
  51.   tags$style(type = "text/css", "html, body {width:100%;height:100%}"),
  52.   leafletOutput("wave_data_map", width = "100%", height = "100%"),
  53.   sliderInput("wave_data_slider", "Wave energy", 
  54.               min = 0, max = 75, 
  55.               value = c(min(wave_data2$An_mn_P_OD), max(wave_data2$An_mn_P_OD)), 
  56.               step = 5, 
  57.               round = 0.5,
  58.               dragRange = TRUE)
  59. )
  61. # 定义服务器逻辑
  62. server <- function(input, output, session) {
  63.   # 创建反应式对象
  64.   wave_energy_output <- reactive({
  65.     wave_data2 %>%
  66.       filter(An_mn_P_OD >= input$wave_data_slider[1]) %>%
  67.       filter(An_mn_P_OD <= input$wave_data_slider[2]) %>%
  68.       group_by(Lat)
  69.   })
  70.   
  71.   # 渲染Leaflet地图
  72.   output$wave_data_map <- renderLeaflet({
  73.     leaflet() %>%
  74.       addProviderTiles(providers$Esri.WorldTopoMap) %>%
  75.       setView(lng = -4.2026458, lat = 56.4906712, zoom = 5) %>%
  76.       addPolygons(data = wave_data2, 
  77.                   weight = 1,
  78.                   smoothFactor = 0.5,
  79.                   color = "white", 
  80.                   fillOpacity = 0.3, 
  81.                   fillColor = pal(wave_data2$Ave_Depth)) %>%
  82.       addLegend(pal = pal, 
  83.                 values = wave_data2$Ave_Depth, 
  84.                 title = "Average depth",
  85.                 labFormat = labelFormat(suffix = "m"),
  86.                 opacity = 0.7, 
  87.                 position = "bottomright")
  88.   })
  89.   
  90.   observeEvent({input$wave_data_slider}, {
  91.     data <- wave_energy_output()[wave_data2$An_mn_P_OD]
  92.     leafletProxy("wave_data_map", data = wave_data2) %>%
  93.       clearShapes() %>%
  94.       addPolygons(data = data,
  95.                   fillColor = pal, 
  96.                   weight = 0.0, 
  97.                   opacity = 1, 
  98.                   color = "white", 
  99.                   dashArray = 3, 
  100.                   fillOpacity = 0.7, 
  101.       ) %>%
  102.       clearControls() %>%
  103.       addLegend(pal = pal, 
  104.                 values = wave_data2$Ave_Depth, 
  105.                 title = "Average depth",
  106.                 labFormat = labelFormat(suffix = "m"),
  107.                 opacity = 0.7, 
  108.                 position = "bottomright")
  109.   })
  110. }
  112. # 运行应用
  113. shinyApp(ui = ui, server = server)

新错误提示：

警告：错误在[:]中：无法在结束之后对列进行子集选择。
ℹ 位置 43、57、57、…、48 和 46 不存在。
ℹ 只有 42 列。

希望这有所帮助。如果您需要进一步的帮助，请随时告诉我。

英文:

I am trying to use the slider function to select and only display selected polygons using Shiny on a Leaflet map.

The code below gets the slider to react, but I want the colour scale to represent the depth, and the slider to select and display the energy polygons (removing non-relevant polygons) but with the depth colour scale. The polygon files contain both the depth and energy data.

The following question points me in the right direction. But I cannot get this to work with my data.

https://stackoverflow.com/questions/53480924/dynamically-render-choropleth-map-with-sliderinput-in-r-shiny/53482005#53482005

Shapefile wave-shp (wave.shp) on the following page: https://www.renewables-atlas.info/downloads/.

Here is my code:

  1. library(shiny)
  2. library(leaflet)
  3. library(rgdal)
  4. library(sf)
  5. library(dplyr)
  6. library(RColorBrewer)
  7. wave_data &lt;- read_sf(&quot;~/path/Wave.shp&quot;)
  8. wave_data &lt;- st_transform(wave_data, crs = &#39;+proj=longlat 
  9. +datum=WGS84&#39;)
  10. wave_data2 &lt;- wave_data %&gt;% 
  11. mutate_at(vars(An_mn_P_OD), ~ as.integer(round(.x))) 
  12. ## Load map
  13. wave_data_map &lt;- leaflet() %&gt;% 
  14. addProviderTiles(providers$Esri.WorldTopoMap) %&gt;% 
  15. setView(lng = -4.2026458, lat = 56.4906712, zoom = 5)
  16. wave_data_map
  17. bins &lt;- c(0, 25, 50, 100, 150, 200, 300, 400, 500, 1000, 3100)
  18. pal &lt;- colorBin(&quot;RdYlBu&quot;, domain = wave_data2$Ave_Depth, bins = 
  19. bins)
  20. ## Add polygons
  21. wave_data_map &lt;- leaflet() %&gt;% 
  22. addProviderTiles(providers$Esri.WorldTopoMap) %&gt;% 
  23. setView(lng = -4.2026458, lat = 56.4906712, zoom = 5) %&gt;% 
  24. addPolygons(data = wave_data2, 
  25. weight = 1,
  26. smoothFactor = 0.5,
  27. color = &quot;white&quot;, 
  28. fillOpacity = 0.5, 
  29. fillColor = pal(wave_data2$Ave_Depth), 
  30. )
  31. wave_data_map
  32. ## Add legend
  33. wave_data_map &lt;- leaflet() %&gt;% 
  34. addProviderTiles(providers$Esri.WorldTopoMap) %&gt;% 
  35. setView(lng = -4.2026458, lat = 56.4906712, zoom = 5) %&gt;% 
  36. addPolygons(data = wave_data2, 
  37. weight = 1,
  38. smoothFactor = 0.5,
  39. color = &quot;white&quot;, 
  40. fillOpacity = 0.3, 
  41. fillColor = pal(wave_data2$Ave_Depth)) %&gt;% 
  42. addLegend(pal = pal, 
  43. values = wave_data2$Ave_Depth, 
  44. title = &quot;Average depth&quot;,
  45. labFormat = labelFormat (suffix = &quot;m&quot;),
  46. opacity = 0.7, 
  47. position = &quot;bottomright&quot;)
  48. wave_data_map
  49. # Define UI for application
  50. ui &lt;- bootstrapPage(
  51. tags$style(type = &quot;text/css&quot;, &quot;html, body 
  52. {width:100%;height:100%}&quot;),
  53. leafletOutput(&quot;wave_data_map&quot;, width = &quot;100%&quot;, height = 
  54. &quot;100%&quot;),
  55. sliderInput(&quot;wave_data_slider&quot;, &quot;Wave energy&quot;, 
  56. min = 0, max = 75, 
  57. value = c(min(wave_data2$An_mn_P_OD), 
  58. max(wave_data2$An_mn_P_OD)), 
  59. step = 5, 
  60. round = 0.5,
  61. dragRange = TRUE)
  62. ))
  63. # Define server logic
  64. server &lt;- function(input, output, session) {
  65. wave_energy_output &lt;- reactive ({
  66. wave_data2 %&gt;% 
  67. filter(An_mn_P_OD &gt;= input$wave_data_range[1]) %&gt;% 
  68. filter(An_mn_P_OD &lt;= input$wave_data_range[2]) %&gt;% 
  69. group_by(Lat)
  70. })
  71. output$wave_data_map &lt;- renderLeaflet(
  72. leaflet() %&gt;% 
  73. addProviderTiles(providers$Esri.WorldTopoMap) %&gt;% 
  74. setView(lng = -4.2026458, lat = 56.4906712, zoom = 5) %&gt;% 
  75. addPolygons(data = wave_data2, 
  76. weight = 1,
  77. smoothFactor = 0.5,
  78. color = &quot;white&quot;, 
  79. fillOpacity = 0.3, 
  80. fillColor = pal(wave_data2$Ave_Depth)) %&gt;% 
  81. addLegend(pal = pal, 
  82. values = wave_data2$Ave_Depth, 
  83. title = &quot;Average depth&quot;,
  84. labFormat = labelFormat (suffix = &quot;m&quot;),
  85. opacity = 0.7, 
  86. position = &quot;bottomright&quot;)
  87. )
  88. observeEvent({input$wave_data_slider}, {
  89. data = wave_energy_output()[wave_data2$An_mn_P_OD]
  90. leafletProxy(&quot;wave_data_map&quot;, data = wave_data2) %&gt;% 
  91. clearShapes() %&gt;% 
  92. addPolygons(data = wave_energy_output,
  93. fillColor = pal, 
  94. weight = 0.0, 
  95. opacity = 1, 
  96. color = &quot;white&quot;, 
  97. dashArray = 3, 
  98. fillOpacity = 0.7, 
  99. ) %&gt;% clearControls() %&gt;% 
  100. addLegend(pal = pal, 
  101. values = wave_data2$Ave_Depth, 
  102. title = &quot;Average depth&quot;,
  103. labFormat = labelFormat (suffix = &quot;m&quot;),
  104. opacity = 0.7, 
  105. position = &quot;bottomright&quot;)
  106. }
  107. )}
  108. # Run the application 
  109. shinyApp(ui = ui, server = server)

NEW ERROR:

Warning: Error in [: Can't subset columns past the end.
ℹ Locations 43, 57, 57, …, 48, and 46 don't exist.
ℹ There are only 42 columns.

Any help would really be appreciated.

答案1

得分: 1

以下是您提供的代码的翻译部分：

  1. 库(shiny)
  2. 库(leaflet)
  3. 库(rgdal)
  4. 库(sf)
  5. 库(dplyr)
  6. 库(RColorBrewer)
  8. wave_data <- read_sf("~/path/Wave.shp")
  10. wave_data <- st_transform(wave_data, crs = '+proj=longlat +datum=WGS84')
  12. wave_data2 <- wave_data %>%
  13.    mutate_at(vars(An_mn_P_OD), ~ as.integer(round(.x))) 
  15. ## 加载地图
  17. wave_data_map <- leaflet() %>%
  18.    addProviderTiles(providers$Esri.WorldTopoMap) %>%
  19.    setView(lng = -4.2026458, lat = 56.4906712, zoom = 5)
  20. wave_data_map
  22. bins <- c(0, 25, 50, 100, 150, 200, 300, 400, 500, 1000, 3100)
  23. pal <- colorBin("RdYlBu", domain = wave_data2$Ave_Depth, bins = bins)
  25. # 定义应用程序的用户界面
  27. ui <- bootstrapPage(
  28.    tags$style(type = "text/css", "html, body {width:100%;height:100%}"),
  29.    leafletOutput("wave_data_map", width = "100%", height = "100%"),
  30.    absolutePanel(top = 10, right = 10, 
  31.    sliderInput("wave_data_range", "Wave energy", 
  32.            min = 0, max = 75, 
  33.            value = c(min(wave_data2$An_mn_P_OD), 
  34.                      max(wave_data2$An_mn_P_OD)), 
  35.            step = 5, 
  36.            round = 0.5,
  37.            dragRange = TRUE)
  38. ))
  40. # 定义服务器逻辑
  42. server <- function(input, output, session) {
  43.    wave_energy_output <- reactive({
  44.       data <- wave_data2 %>% 
  45.          filter(An_mn_P_OD >= input$wave_data_range[1],
  46.             An_mn_P_OD <= input$wave_data_range[2])
  47.       data %>% group_by(Lat)
  48.    })
  50.   output$wave_data_map <- renderLeaflet({
  51.    leaflet() %>% 
  52.      addProviderTiles(providers$Esri.WorldTopoMap) %>% 
  53.      setView(lng = -4.2026458, lat = 56.4906712, zoom = 5) %>% 
  54.      addPolygons(data = wave_data2, 
  55.                  weight = 1,
  56.                  smoothFactor = 0.5,
  57.                  color = "white", 
  58.                  fillOpacity = 0.3, 
  59.                  fillColor = pal(wave_data2$Ave_Depth)) %>% 
  60.      addLegend(pal = pal, 
  61.                values = wave_data2$Ave_Depth, 
  62.                title = "Average depth",
  63.                labFormat = labelFormat (suffix = "m"),
  64.                opacity = 0.7, 
  65.                position = "bottomright")
  66.  })
  68.  observeEvent(input$wave_data_range, {
  69.     data <- wave_energy_output()
  70.     leafletProxy("wave_data_map", data = data) %>% 
  71.      clearShapes() %>% 
  72.      addPolygons(data = data,
  73.                  fillColor = pal(wave_energy_output()$Ave_Depth), 
  74.                  weight = 1,
  75.                  smoothFactor = 0.5,
  76.                  color = "white", 
  77.                  fillOpacity = 0.3)
  78.    })
  80. }
  82. # 运行应用程序
  84. shinyApp(ui = ui, server = server)

希望这有所帮助。如果您有任何其他问题或需要进一步的协助，请随时告诉我。

英文:

Working code:

  1. library(shiny)
  2. library(leaflet)
  3. library(rgdal)
  4. library(sf)
  5. library(dplyr)
  6. library(RColorBrewer)
  7. wave_data &lt;- read_sf(&quot;~/path/Wave.shp&quot;)
  8. wave_data &lt;- st_transform(wave_data, crs = &#39;+proj=longlat 
  9. +datum=WGS84&#39;)
  10. wave_data2 &lt;- wave_data %&gt;% 
  11. mutate_at(vars(An_mn_P_OD), ~ as.integer(round(.x))) 
  12. ## Load map
  13. wave_data_map &lt;- leaflet() %&gt;% 
  14. addProviderTiles(providers$Esri.WorldTopoMap) %&gt;% 
  15. setView(lng = -4.2026458, lat = 56.4906712, zoom = 5)
  16. wave_data_map
  17. bins &lt;- c(0, 25, 50, 100, 150, 200, 300, 400, 500, 1000, 3100)
  18. pal &lt;- colorBin(&quot;RdYlBu&quot;, domain = wave_data2$Ave_Depth, bins = bins)
  19. # Define UI for application
  20. ui &lt;- bootstrapPage(
  21. tags$style(type = &quot;text/css&quot;, &quot;html, body 
  22. {width:100%;height:100%}&quot;),
  23. leafletOutput(&quot;wave_data_map&quot;, width = &quot;100%&quot;, height = 
  24. &quot;100%&quot;),
  25. absolutePanel(top = 10, right = 10, 
  26. sliderInput(&quot;wave_data_range&quot;, &quot;Wave energy&quot;, 
  27. min = 0, max = 75, 
  28. value = c(min(wave_data2$An_mn_P_OD), 
  29. max(wave_data2$An_mn_P_OD)), 
  30. step = 5, 
  31. round = 0.5,
  32. dragRange = TRUE)
  33. ))
  34. # Define server logic
  35. server &lt;- function(input, output, session) {
  36. wave_energy_output &lt;- reactive({
  37. data &lt;- wave_data2 %&gt;% 
  38. filter(An_mn_P_OD &gt;= input$wave_data_range[1],
  39. An_mn_P_OD &lt;= input$wave_data_range[2])
  40. data %&gt;% group_by(Lat)
  41. })
  42. output$wave_data_map &lt;- renderLeaflet({
  43. leaflet() %&gt;% 
  44. addProviderTiles(providers$Esri.WorldTopoMap) %&gt;% 
  45. setView(lng = -4.2026458, lat = 56.4906712, zoom = 5) %&gt;% 
  46. addPolygons(data = wave_data2, 
  47. weight = 1,
  48. smoothFactor = 0.5,
  49. color = &quot;white&quot;, 
  50. fillOpacity = 0.3, 
  51. fillColor = pal(wave_data2$Ave_Depth)) %&gt;% 
  52. addLegend(pal = pal, 
  53. values = wave_data2$Ave_Depth, 
  54. title = &quot;Average depth&quot;,
  55. labFormat = labelFormat (suffix = &quot;m&quot;),
  56. opacity = 0.7, 
  57. position = &quot;bottomright&quot;)
  58. })
  59. observeEvent(input$wave_data_range, {
  60. data &lt;- wave_energy_output()
  61. leafletProxy(&quot;wave_data_map&quot;, data = data) %&gt;% 
  62. clearShapes() %&gt;% 
  63. addPolygons(data = data,
  64. fillColor = pal(wave_energy_output()$Ave_Depth), 
  65. weight = 1,
  66. smoothFactor = 0.5,
  67. color = &quot;white&quot;, 
  68. fillOpacity = 0.3)
  69. })
  70. }
  71. # Run the application 
  72. shinyApp(ui = ui, server = server)

huangapple
  • 本文由 发表于 2023年2月19日 23:39:12
  • 转载请务必保留本文链接：https://go.coder-hub.com/75501282.html
  • qgis
  • r
  • r-leaflet
  • shiny
  • shinyapps
匿名

发表评论

匿名网友

:?: :razz: :sad: :evil: :!: :smile: :oops: :grin: :eek: :shock: :???: :cool: :lol: :mad: :twisted: :roll: :wink: :idea: :arrow: :neutral: :cry: :mrgreen:

确定