Объединить и построить несколько изохрон

Действительно классный вопрос. Что вы хотите сделать, так это объединить фигуры по идентификатору, чтобы все области 0–30 минут были одной фигурой, все области 30–60 минут — другой и так далее. Есть способы сделать это с помощью других пространственных пакетов, но, похоже, он хорошо подходит для sf, который использует функции в стиле dplyr.

После того, как вы создадите isochrone, вы можете преобразовать его в объект sf, сделать метку расстояния того же типа, сгруппировать по идентификатору и вызвать summarise. По умолчанию, когда вы суммируете объекты sf, это просто пространственное объединение, поэтому вам не нужно указывать там функцию.

library(sf)
library(dplyr)

iso_sf <- st_as_sf(isochrone)

iso_union <- iso_sf %>%
  mutate(label = paste(min, max, sep = "-")) %>%
  group_by(id, label) %>%
  summarise()

У меня не было leaflet под рукой, поэтому вот только метод печати по умолчанию:

plot(iso_union["label"], pal = RColorBrewer::brewer.pal(4, "Blues"))

Я не уверен, что случилось с областями с резкими вертикальными краями, но они есть и на вашем графике.

Мне было трудно использовать метод map2, который вы использовали, потому что он выполняет как объединение, так и, я думаю, другую теорию множеств, например, функцию для создания определенных интервалов. Вместо этого я бы рекомендовал создать растровый слой из созданных вами слоев и применить одну непрозрачность к этому одному растру, как это делает пример ggmap. В блоге есть отличный пост о том, что я украл много кода у здесь (вместе с пользователем: camille).

Он использует другой API, для которого требуется mapbox, но он бесплатный. Еще одно ограничение заключается в том, что он не будет возвращать изокроны нужного вам размера, но я воссоздал его в другом месте, где три точки расположены ближе друг к другу, чтобы проверить метод.

Я также не стал векторизовать процесс создания веб-запроса isocrone, поэтому я оставляю это кому-то более умному.

# First be sure to get your mapbox token

library(fasterize)
library(sf)
library(mapboxapi)
library(leaflet)
#mapboxapi::mb_access_token("Go get the token and put it here",
#                           install = TRUE, overwrite = TRUE)

isos1 <- mb_isochrone(
  location = c("-149.883234, 61.185765"),
  profile = "driving",
  time = c(5,10,15),
)

isos2 <- mb_isochrone(
  location = c("-149.928200, 61.191227"),
  profile = "driving",
  time = c(5,10,15),
)
isos3 <- mb_isochrone(
  location = c("-149.939484, 61.160192"),
  profile = "driving",
  time = c(5,10,15),
)

library(sf)
library(dplyr)

isocrones <- rbind(isos1,isos2,isos3)

iso_sf <- st_as_sf(isocrones)

iso_union <- iso_sf %>%
  group_by(time) %>%
  summarise()

isos_proj <- st_transform(iso_sf, 32615)

template <- raster(isos_proj, resolution = 100)

iso_surface <- fasterize(isos_proj, template, field = "time", fun = "min")

pal <- colorNumeric("viridis", isos_proj$time, na.color = "transparent")
leaflet() %>%
  addTiles() %>%
  addRasterImage(iso_surface, colors = pal, opacity = 0.5) %>%
  addLegend(values = isos_proj$time, pal = pal,
            title = "Minutes of Travel") %>% 
  addMarkers(lat = c(61.185765, 61.191227, 61.160192), lng = c(-149.883234, -149.928200, -149.939484))