leaflet.extras2

• Based on https://stackoverflow.com/questions/73527927/r-shiny-leaflet-how-to-query-a-web-map-service-layer-with-a-fake-click?rq=1
• using addWMS
• from https://github.com/trafficonese/leaflet.extras2

library(magrittr)
library(leaflet)
library(leaflet.extras2)
library(stringr)

wms_layer <- "https://catalogue.ceh.ac.uk/maps/51bcb92a-dd88-4034-ba65-a9d432dd632a?request=getCapabilities&service=WMS&cache=false&"

xyz_layer <- "https://tiles.globalforestwatch.org/wri_trees_in_mosaic_landscapes/v20220922/default/{z}/{x}/{y}.png"

leaflet() %>%
  addProviderTiles("Esri.WorldImagery",
                   group = "Esri.WorldImagery",
                   options = providerTileOptions(zIndex = 0)) %>%
  setView(lng = 0, lat = 30, zoom = 2.5) %>%
  addWMS(
    wms_layer,
    layers = "LC.10m.GB", # Or "LC.10m.NI" for northern Ireland
    group = "Great Britain Land cover map",
    options = WMSTileOptions(
      format = "image/png",
      version = "1.3.0",
      transparent = TRUE,
      opacity = 0.5,# Add some transparency so that we can still see the satellite image
      info_format = "application/vnd.ogc.gml"
    ),
    popupOptions = popupOptions(maxWidth = 300, closeOnClick = T)) %>%
  addTiles(
    urlTemplate = xyz_layer,
    group = "WRI Trees in Mosaic Landscapes",) %>%
  addLayersControl(
    overlayGroups = c("Great Britain Land cover map", "WRI Trees in Mosaic Landscapes"),
    options = layersControlOptions(collapsed = FALSE)
  )