Как использовать WMS в R?

Мне нравится реализовать специальный WMS-слой с leaflet() в R. Но я не добился успеха в WMS-слое, который мне нравится использовать.

library("leaflet")
library("sp")

# this example from the tutorial works
leaflet() %>% addTiles() %>% setView(-93.65, 42.028, zoom = 4) %>%
  addWMSTiles(
    "http://mesonet.agron.iastate.edu/cgi-bin/wms/nexrad/n0r.cgi",
    layers = "nexrad-n0r",
    options = WMSTileOptions(format = "image/png", transparent = TRUE) 
)

# but for my implementation it don't work
leaflet() %>% addTiles() %>% setView(lat = 53.8, lng = 12.6, zoom = 8) %>% 
  addWMSTiles("http://www.umweltkarten.mv-regierung.de/script/mv_a3_wasserschutz_wms.php?", 
              layers = "t2_bib_p", 
              options = WMSTileOptions(format = "image/png")
)

Я предполагаю, что существует проблема с системой координат. Поскольку WMS-слой, который мне нравится использовать, не имеет стандартного CRS. Поэтому я пытаюсь экспериментировать с этим WMSTileOptions():

crs=CRS("+init=epsg:5650")

Идея состоит в том, чтобы использовать WMS-Layers в качестве фоновых изображений для моих векторных данных. Векторные данные поступают в основном из файлов формы и GPS-данных (текстовых файлов). Целью является создание карт, которые я могу использовать в отчете. Если вы знаете другие возможности рядом с leaflet для использования WMS в R, также приветствуются?

Ответ 1

Я замечаю, что если вы исследуете этот базовый URL, он говорит

msWMSDispatch(): ошибка сервера WMS. Неполный запрос WMS: VERSION отсутствует параметр

и, исследуя веб-сайт, я определил, что правильный формат:

https://www.umweltkarten.mv-regierung.de/script/mv_a3_wasserschutz_wms.php?SERVICE=WMS&REQUEST=GetCapabilities&VERSION=1.3.0

который дает довольно много данных WMS, если вы поместите их в свой браузер.

Мы можем использовать "найти", чтобы узнать, какие слои доступны для запроса. Есть несколько таких, как t3_ksg, t3_wsg и lungwms, но кажется, что t2_bib_p недоступен для mv_a3_wasserschutz_wms.