Un raster est un fichier géographique dont les données sont stockées sous forme d’une grille régulière. Il existe plusieurs format de stockage : GeoTiff, netcdf, asc, …
Dans R, il existe 3 packages principaux qui peuvent vous permettre de lire un fichier raster :
- l’historique : {raster}
- la version remastérisée avec du C++ : {terra}
- un package récent avec une philosophie différente : {stars}
Lecture d’un fichier raster avec une seule couche
# récupération du chemin d'un Geotiff
tif <- system.file("tif/L7_ETMs.tif", package = "stars")
# lecture avec {raster}
r_raster <- raster(tif)
# lecture avec {terra}
r_terra <- rast(tif)
# lecture avec {stars}
r_stars <- read_stars(tif)
Lecture d’un fichier raster avec plusieurs couches
# récupération du chemin d'un Geotiff
multiband <- system.file("ex/logo.tif", package="terra")
# lecture avec {raster}
r_raster <- stack(multiband)
# lecture avec {terra}
r_terra <- rast(multiband)
# lecture avec {stars}
r_stars <- read_stars(multiband)
Pour le choix du package à utiliser, tout dépend de votre utilisation du raster. Explorez les fonctions existantes et la doc pour choisir.
Si vous cherchiez comment lire un fichier shapefile plutôt, c’est par là .