here’s the equivalent using leaflet in R, displaying the iNat heatmap tiles over a Carto DarkMatter basemap:
library(jsonlite)
library(leaflet)
reqparams <- "project_id=gbs-2024-uruguay"
bounds <- fromJSON(paste0("https://api.inaturalist.org/v1/observations?per_page=0&return_bounds=true&",reqparams))
leaflet() %>%
fitBounds(lng1 = bounds$total_bounds$swlng, lat1 = bounds$total_bounds$swlat, lng2 = bounds$total_bounds$nelng, lat2 = bounds$total_bounds$nelat) %>%
addProviderTiles("CartoDB.DarkMatter", group = "Carto Dark Matter") %>%
addTiles(urlTemplate = paste0("https://api.inaturalist.org/v1/heatmap/{z}/{x}/{y}.png?",reqparams), group = "iNaturalist Observations Heatmap") %>%
addLayersControl(
baseGroups = c("Carto Dark Matter"),
overlayGroups = c("iNaturalist Observations Heatmap")
)