leaflet 並排用於 R 中的 2 個光柵圖像

Question

{leaflet.extras2}實現了 leaflet並排插件並提供了這個最小示例：

library(leaflet)
library(leaflet.extras2)

leaflet(quakes) %>%
  addMapPane("left", zIndex = 0) %>%
  addMapPane("right", zIndex = 0) %>%
  addTiles(group = "base", layerId = "baseid",
           options = pathOptions(pane = "right")) %>%
  addProviderTiles(providers$CartoDB.DarkMatter, group="carto", layerId = "cartoid",
                   options = pathOptions(pane = "left")) %>%
  addCircleMarkers(data = breweries91[1:15,], color = "blue", group = "blue",
                   options = pathOptions(pane = "left")) %>%
  addCircleMarkers(data = breweries91[15:20,], color = "yellow", group = "yellow") %>%
  addCircleMarkers(data = breweries91[15:30,], color = "red", group = "red",
                   options = pathOptions(pane = "right")) %>%
  addLayersControl(overlayGroups = c("blue","red", "yellow")) %>%
  addSidebyside(layerId = "sidecontrols",
                rightId = "baseid",
                leftId = "cartoid")

但是，在適應兩個raster圖像時，我遇到了錯誤：

Error in addRasterImage(., r2, colors = pal2, options = pathOptions(pane = "right"),  : 
  unused argument (options = pathOptions(pane = "right"))

這是一個可重現的示例：

library(raster)
library(leaflet)
library(leaflet.extras2)
library(rcartocolor)

# example raster data
r1 <- r2 <- raster(system.file("external/test.grd", package = "raster"))

# color palettes
pal1 <- colorNumeric(carto_pal(name = "OrYel"), 
                     values(r1), na.color = "transparent")
pal2 <- colorNumeric(carto_pal(name = "BluYl"), 
                     values(r2), na.color = "transparent")

# side by side map
leaflet() %>% 
  addMapPane("left",  zIndex = 0) %>%
  addMapPane("right", zIndex = 0) %>%
  addTiles(group = "base", layerId = "baseid",
           options = pathOptions(pane = "right")) %>%
  addTiles(group = "carto", layerId = "cartoid",
           options = pathOptions(pane = "left")) %>%
  addRasterImage(r1, colors = pal1, 
                 options = pathOptions(pane = "right"), group = "r1") %>%
  addRasterImage(r2, colors = pal2, 
                 options = pathOptions(pane = "right"), group = "r2") %>%
  addLayersControl(overlayGroups = c("r1", "r2")) %>% 
  addSidebyside(layerId = "sidecontrols",
                rightId = "baseid",
                leftId  = "cartoid")

Answer 1

@RexParsons 在這里指出了解決方案。
您將需要addRasterImage的options = leafletOptions(pane = "")將在下一個版本中添加到{leaflet} 。
第 1 步：現在您可以安裝raster-options ：

remotes::install_github("rstudio/leaflet", ref="joe/feature/raster-options")

第 2 步： Session -> Restart R

第 3 步：添加options = leafletOptions(pane = "") ：

leaflet() |> 
  addMapPane("right", zIndex = 0) |> 
  addMapPane("left",  zIndex = 0) |>
  addTiles(group = "base", layerId = "baseid1", options = pathOptions(pane = "right")) |> 
  addTiles(group = "base", layerId = "baseid2", options = pathOptions(pane = "left")) |> 
  addRasterImage(x = r1, colors = pal1, options = leafletOptions(pane = "right"), group = "r1") |> 
  addRasterImage(x = r2, colors = pal2, options = leafletOptions(pane = "left"), group = "r2") |> 
  addLayersControl(overlayGroups = c("r1", "r2")) |>  
  addSidebyside(layerId = "sidecontrols",
                rightId = "baseid1",
                leftId  = "baseid2")

Output：

選項 2（使用addProviderTiles ）：

leaflet() |> 
  addMapPane("right", zIndex = 0) |> 
  addMapPane("left",  zIndex = 0) |>
  addTiles(group = "base", layerId = "baseid", options = pathOptions(pane = "right")) |>
  addProviderTiles(providers$CartoDB, group="carto", layerId = "cartoid", options = pathOptions(pane = "left")) |> 
  addRasterImage(x = r1, colors = pal1, options = leafletOptions(pane = "right"), group = "r1") |> 
  addRasterImage(x = r2, colors = pal2, options = leafletOptions(pane = "left"), group = "r2") |> 
  addLayersControl(overlayGroups = c("r1", "r2")) |>  
  addSidebyside(layerId = "sidecontrols",
                rightId = "baseid",
                leftId  = "cartoid")

Output：

對你起作用嗎？

<-- 下面對@oatmilkyway 問題的部分回答-->
有人需要找到一種方法來（有條件地）添加/刪除leaflet.extras2::addSidebyside的layerId作為可能的解決方案之一（如果已經存在，請告訴我）。

leaflet() |> 
  addMapPane("right", zIndex = 0) |> 
  addMapPane("left", zIndex = 0) |>
  addTiles(group = "Base map", layerId = "baseid1", options = pathOptions(pane = "right")) |> 
  addTiles(group = "Base map", layerId = "baseid2", options = pathOptions(pane = "left")) |> 
  addProviderTiles(provider = providers$Esri.WorldImagery, group = "Satellite map", layerId = "satelliteid1", options = pathOptions(pane = "right")) |>
  addProviderTiles(provider = providers$Esri.WorldImagery, group = "Satellite map", layerId = "satelliteid2", options = pathOptions(pane = "left")) |>
  addRasterImage(x = r1, colors = pal1, options = leafletOptions(pane = "right"), group = "raster image 1") |> 
  addRasterImage(x = r2, colors = pal2, options = leafletOptions(pane = "left"), group = "raster image 2") |> 
  addLayersControl(baseGroups = c("Base map", "Satellite map"),
                   overlayGroups = c("raster image 1", "raster image 2")) |>  
  addSidebyside(layerId = "sidecontrols2", rightId = "satelliteid1", leftId  = "satelliteid2") |>
  addSidebyside(layerId = "sidecontrols", rightId = "baseid1", leftId  = "baseid2")

leaflet() |> 
  addMapPane("right", zIndex = 0) |> 
  addMapPane("left", zIndex = 0) |>
  addTiles(group = "Base and Carto Positron maps", layerId = "baseid", options = pathOptions(pane = "right")) |>
  addProviderTiles(provider = providers$CartoDB.Positron, group = "Base and Carto Positron maps", layerId = "cartoid", options = pathOptions(pane = "left")) |>
  addProviderTiles(provider = providers$Esri.WorldImagery, group = "Satellite and Toner Lite maps", layerId = "satelliteid", options = pathOptions(pane = "right")) |>
  addProviderTiles(provider = providers$Stamen.TonerLite, group = "Satellite and Toner Lite maps", layerId = "tonerlightid", options = pathOptions(pane = "left")) |>
  addRasterImage(x = r1, colors = pal1, options = leafletOptions(pane = "right"), group = "raster image 1") |> 
  addRasterImage(x = r2, colors = pal2, options = leafletOptions(pane = "left"), group = "raster image 2") |>
  addLayersControl(baseGroups = c("Base and Carto Positron maps", "Satellite and Toner Lite maps"),
                   overlayGroups = c("raster image 1", "raster image 2"),
                   options = layersControlOptions(collapsed = FALSE)) |> 
  addSidebyside(layerId = "sidecontrols2", rightId = "satelliteid", leftId  = "tonerlightid") |>
  addSidebyside(layerId = "sidecontrols", rightId = "baseid", leftId  = "cartoid")

Output：

盡管這充其量是一個笨拙的答案，但我希望它可以幫助您。