210127_artikkel/matsalu_wms.R

library(leaflet)

year <- format(Sys.Date(), "%Y")


leaflet() %>% addTiles() %>% setView(-93.65, 42.0285, zoom = 7) %>%addWMSTiles(
  "http://basemap.nationalmap.gov/arcgis/services/USGSHydroNHD/MapServer/WMSServer?",
  layers = "0",
  options = WMSTileOptions(format = "image/png", transparent = TRUE),
  attribution = "") 

i <- 9
# Ajaloliste kaartide kihid
layers <- c(
  "yheverstakaart_alt", "kaheverstakaart_alt", "kolmeverstakaart_alt",
  "HYBRID", "Halduspiir_raster", "vanaBaaskaart", "ekj_50T", "kk1940", "lehman",
  "ew_25T", "ew_50T", "ew_200T", "pk_trykk", "heereskarte"
)
url <- paste0("https://xgis.maaamet.ee/xgis2/service/c9me/ajal?SERVICE=WMS&VERSION=1.1.1&REQUEST=GetMap&FORMAT=image%2Fpng&TRANSPARENT=true&LAYERS=", layers[i], "&SRS=EPSG%3A3301&STYLES=&WIDTH=1236&HEIGHT=687&BBOX=450632.44047619053%2C6404438.244047619%2C459828.86904761905%2C6409549.851190476")
cat(url)

wms_map_url <- "https://kaart.maaamet.ee/wms/ajalooline?version=1.3.0"
wms_map_url <- "http://basemap.nationalmap.gov/arcgis/services/USGSHydroNHD/MapServer/WMSServer?"

x <- 23.4784
y <- 58.7864


## WMS: transpordivõrgud
wms_map_url <- "https://inspire.maaamet.ee/arcgis/rest/services/public/tn_com/MapServer/exts/InspireView/service"
leaflet() %>% addTiles() %>% setView(x, y, zoom = 14) %>% addWMSTiles(
  wms_map_url,
  layers = c("1","2"),
  options = WMSTileOptions(format = "image/png", transparent = TRUE),
  attribution = "Maa-amet") 


## WMS: Ortofoto
wms_map_url <- "https://inspire.maaamet.ee/arcgis/rest/services/public/orthophoto/MapServer/exts/InspireView/service"
leaflet() %>% addTiles() %>% setView(x, y, zoom = 15) %>%addWMSTiles(
  wms_map_url,
  layers = "0",
  options = WMSTileOptions(format = "image/png", transparent = TRUE),
  attribution = paste0("Aluskaart: Maa-amet ", year, " "))
#' Ortofotode allalaadimine: https://geoportaal.maaamet.ee/est/Ruumiandmed/Ortofotod/Laadi-ortofotod-alla-p610.html
#' Ruudustik 1:10000, näiteks 'matsalu ruut' on 62173.
#' 




## WMS: Ortofotod
k <- 0
wms_map_url <- "https://inspire.maaamet.ee/arcgis/rest/services/public/orthophoto/MapServer/exts/InspireView/service"
leaflet() %>% addTiles() %>% setView(x, y, zoom = 15) %>%addWMSTiles(
  wms_map_url,
  layers = k,
  options = WMSTileOptions(format = "image/png", transparent = TRUE),
  attribution = "Maa-amet")
k <- k + 1


# Tuumaelektrijaamad
leaflet() %>% addTiles() %>% setView(x, y, zoom = 4) %>%  addWMSTiles(
  "http://sedac.ciesin.columbia.edu/geoserver/wms",
  layers = "energy:energy-pop-exposure-nuclear-plants-locations_plants",
  options = WMSTileOptions(format = "image/png", transparent = TRUE),
  tileOptions(tms = TRUE),
  attribution = "")