Utilisation des flux WFS depuis R

R
géomatique

Dans cet ébauche de tutoriel, la démarche de base pour télécharger des données à partir d’un flux WFS depuis R est donnée.

Auteur·rice

Frédéric Grelot

Date de publication

10 décembre 2024

WFS, de quoi parle-t-on ?

Les Web Feature Service (WFS) sont un service internet d’intéropérabilité défini par l’open gesospatial consortium (OGC) pour permettre le téléchargement des données géomatiques.

Utilisation avec R

Installation nécessaire

sudo apt install libsodium-dev
sudo R -e "install.packages('ows4R')"

Exemple d’utilisation

library("ows4R")
service = "https://geolittoral.din.developpement-durable.gouv.fr/wxs"

client = ows4R::WFSClient[["new"]](service, "2.0.0")
client[["getFeatureTypes"]](pretty = TRUE)

feature = client[["getFeatureTypes"]](pretty = TRUE)
typeName = feature[48, 1]

sapply(
    client[["getCapabilities"]]()[["findFeatureTypeByName"]](typeName)[["getDescription"]](),
    function(x){x[["getName"]]()}
)

sapply(
    client[["getCapabilities"]]()[["getOperationsMetadata"]]()[["getOperations"]](),
    function(x){x[["getParameter"]]()}
)
# toto = ows4R::WFSGetFeature[["new"]](client, url = service, version = "2.0.0", typeName = typeName)

url = httr::parse_url(service)
url[["query"]] = list(
    service = "wfs",
    request = "GetFeature",
    typename = typeName,
    version = "2.0.0"
    )
request = httr::build_url(url)

toto = sf::read_sf(request)
Retour au sommet