如何在R语言Leaflet地图中添加日期滑块选择RasterLayer

我有一个ERA5降水数据的RasterStack，目前只能提取单个RasterLayer绘制，代码如下：

precipitation_data <- raster::stack(era5_precip_filepath)
precipitation_data <- raster::projectRaster(precipitation_data, crs = "+proj=utm +zone=32")

selected_month <- precipitation_data$X2021.10.01

leaflet::leaflet() %>% 
leaflet::addTiles() %>%
leaflet::addRasterImage(x = selected_month, opacity = 1)

我希望通过日期滑块选择RasterStack中的图层进行绘制，所有图层名称都包含日期信息：

> names(era5_precip)
 [1] "X2021.10.01" "X2021.10.02" "X2021.10.03" "X2021.10.04" "X2021.10.05" "X2021.10.06" "X2021.10.07" "X2021.10.08" "X2021.10.09" "X2021.10.10" "X2021.10.11" "X2021.10.12" "X2021.10.13" "X2021.10.14" "X2021.10.15" "X2021.10.16"
[17] "X2021.10.17" "X2021.10.18" "X2021.10.19" "X2021.10.20" "X2021.10.21" "X2021.10.22" "X2021.10.23" "X2021.10.24" "X2021.10.25" "X2021.10.26" "X2021.10.27" "X2021.10.28" "X2021.10.29" "X2021.10.30" "X2021.10.31"

想知道有没有简便的实现方法？

最简便的方式是用shiny框架结合leaflet实现交互式滑块，直接绑定图层日期选择，代码示例如下：

library(shiny)
library(raster)
library(leaflet)

# 加载并投影数据
precipitation_data <- stack(era5_precip_filepath)
precipitation_data <- projectRaster(precipitation_data, crs = "+proj=utm +zone=32")

# 处理图层名称，提取日期（去掉开头的X，替换.为-）
date_labels <- gsub("^X", "", names(precipitation_data))
date_labels <- gsub("\\.", "-", date_labels)

# 构建Shiny应用
ui <- fluidPage(
  titlePanel("ERA5降水数据日期滑块选择"),
  sidebarLayout(
    sidebarPanel(
      sliderInput("date_slider", "选择日期：",
                  min = 1, max = nlayers(precipitation_data),
                  value = 1, step = 1,
                  # 用处理后的日期作为滑块标签
                  animate = TRUE,
                  labels = date_labels)
    ),
    mainPanel(
      leafletOutput("precip_map")
    )
  )
)

server <- function(input, output) {
  output$precip_map <- renderLeaflet({
    # 根据滑块选择获取对应图层
    selected_raster <- precipitation_data[[input$date_slider]]
    
    leaflet() %>%
      addTiles() %>%
      addRasterImage(x = selected_raster, opacity = 1)
  })
}

shinyApp(ui, server)

关键说明： #

• 先处理图层名称，把X2021.10.01转换成2021-10-01，让滑块标签更直观
• 开启sliderInput的animate参数，可一键播放时序降水数据
• 服务器端根据滑块索引直接提取RasterStack对应图层，实时渲染到地图上

如果不需要完整的Shiny应用，也可以用leaflet.extras的图层控件做切换，但滑块交互的实现还是Shiny最直接高效。

内容的提问来源于stack exchange，提问作者oxidising