Topotijdreis in R

[Lennart]

Hallo,

Ik wil topotijdreis ook gebruiken in mijn R script. Het lukt me helaas niet. Kan iemand mij misschien helpen?

Door dit stappenplan, van Richard Duivenvoorde, te volgen is het me wel gelukt om de kaarten te zien in QGis.
Ik heb ook deze post op het geoforum bekeken maar toen kwam ik er helaas ook niet uit.

Dit is mijn R-code:

library(leaflet)

leaflet() %>%
  # Set view to the Netherlands
  setView(5.41077,52.13012,zoom = 8) %>%
  # Add openstreetmap tiles
  addTiles(group = "OSM") %>%
  # Add topotijdreis wmts tiles
  addWMSTiles("https://tiles.arcgis.com/tiles/nSZVuSZjHpEZZbRo/arcgis/rest/services/Historische_tijdreis_1950/MapServer/WMTS?"
              ,layers = "Historische_tijdreis_1950"
              ,options = WMSTileOptions(format = "image/jpgpng"
                                        ,crs = "EPSG:28992"
                                        ,version = "1.0.0"
              )
              ,group = "Topotijdreis"
  ) %>%
  # Add layer control
  addLayersControl(
    baseGroups = c("OSM"),
    overlayGroups = c("Topotijdreis"), 
    options = layersControlOptions(collapsed = TRUE)
  )

Alvast super bedankt!

Groet,
Lennart Jongen

[Lennart]

Ik heb het antwoord gevonden

library(leaflet)
library(leaflet.esri)

epsg28992 <- leafletCRS(
  crsClass = "L.Proj.CRS"
  ,code = "EPSG:28992"
  ,proj4def = "+proj=sterea +lat_0=52.15616055555555 +lon_0=5.38763888888889 +k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel +units=m +towgs84=565.2369,50.0087,465.658,-0.406857330322398,0.350732676542563,-1.8703473836068,4.0812 +no_defs"
  ,origin=c(-3.05155E7,3.1112399999999993E7)
  ,resolutions = c(3251.206502413005,1625.6032512065026,812.8016256032513,406.40081280162565,203.20040640081282,101.60020320040641, 50.800101600203206,25.400050800101603,12.700025400050801,6.350012700025401,3.1750063500127004,1.5875031750063502,0.7937515875031751,0.39687579375158755,0.19843789687579377,0.09921894843789689,0.04960947421894844)
)

leaflet(options = leafletOptions(crs = epsg28992)) %>%
  setView(229684.400,513480.640,zoom = 10) %>%
  addEsriTiledMapLayer("https://tiles1.arcgis.com/tiles/nSZVuSZjHpEZZbRo/arcgis/rest/services/Historische_tijdreis_1862/MapServer/"
                       ,group = "Layer_1862" 
  )