Edina WMS

library(leaflet.extras)
library(mapview)
library(giscourse)

mapviewOptions(basemaps =c("Esri.WorldImagery","OpenStreetMap","Esri.WorldStreetMap", "CartoDB.Positron", "OpenTopoMap", "OpenStreetMap.BlackAndWhite", "OpenStreetMap.HOT","Stamen.TonerLite", "Esri.WorldShadedRelief", "Stamen.Terrain"))

#ky<-"place your own api token key here and uncomment"
ws1<-"http://openstream.edina.ac.uk/openstream/wms"
ws2<-"https://environment.data.gov.uk/spatialdata/lidar-composite-digital-surface-model-dtm-2m/wms"

layer1<-"osopendata"
layer2<- "LIDAR_Composite_DTM_2m"
m<-qmap() 
m@map %>% 
  addWMSTiles(group=layer1, ws1, layers=layer1,options = WMSTileOptions(format = "image/png", transparent = TRUE, token=ky)) %>% 
  addWMSTiles(group=layer2, ws2, layers=layer2,options = WMSTileOptions(format = "image/png", transparent = TRUE))%>% 
   addMeasure(primaryLengthUnit ='meters', secondaryLengthUnit='kilometers',primaryAreaUnit='hectares', secondaryAreaUnit='acres',position="topleft") %>% 
  mapview:::mapViewLayersControl (names = c(layer2,layer1))