使用R从Web功能服务（WFS）读取信息

我喜欢使用Web功能服务（WFS）[1]获取土壤信息。但我很不确定如何做对，也许你可以给我一个提示，如何获得所有可用的信息

我可以查询

bmf:Musterstueck

并将其转换为简单的功能。但是如果我对

so:SoilHorizon

或

so:ObservedSoilProfile

执行相同操作，我会得到如下错误消息

Error: Cannot open "https://via.bund.de/bmf/inspire/so/wfs?service=wfs&version=2.0.0&request=GetFeature&typename=so%3ASoilHorizon"; The file doesn't seem to exist.
In addition: Warning message:
In CPL_read_ogr(dsn, layer, query, as.character(options), quiet,  :
  GDAL Error 1: HTTP error code : 403 

多谢各位 克里斯托夫

[1]

库（tidyverse）
图书馆（sf）
图书馆（ows4R）
图书馆（httr）
wfs=WFSClient$new（“https://via.bund.de/bmf/inspire/so/wfs“，“2.0.0”，logger=“INFO”）
wfs$getFeatureTypes（pretty=TRUE）
url为什么不使用QGIS这样的GIS？您是否尝试过以这种方式查询？403错误表示您没有查看该数据的权限
                  name                           title
      bmf:Musterstueck Musterstücke der Bodenschätzung
        so:SoilHorizon             INSPIRE SoilHorizon
so:ObservedSoilProfile     INSPIRE ObservedSoilProfile

library(tidyverse)
library(sf)
library(ows4R)
library(httr)

wfs = WFSClient$new("https://via.bund.de/bmf/inspire/so/wfs", "2.0.0", logger = "INFO")
wfs$getFeatureTypes(pretty = TRUE)

url <- parse_url("https://via.bund.de/bmf/inspire/so/wfs")

url$query <- list(service  = "wfs", 
                  version  = "2.0.0", 
                  request  = "GetFeature", 
                  typename = "bmf:Musterstueck")  

request <- build_url(url)
request

bmf_musterstueck <- read_sf(request)
bmf_musterstueck