Lire un shapefile
Les shapefiles sont l'un des formats les plus courants pour partager des données spatiales et ils se lisent facilement dans R à l'aide de readOGR() du paquet rgdal. readOGR() a deux arguments importants : dsn et layer. Ce que vous passez à ces arguments dépend du type de données à importer. Vous avez vu dans la vidéo que, pour les shapefiles, dsn doit être le chemin du répertoire qui contient les fichiers constituant le shapefile, et layer est le nom de fichier du shapefile (sans extension).
Pour votre carte, vous souhaitez les limites des quartiers. Nous avons téléchargé les « Neighborhood Tabulation Areas », telles que définies par la Ville de New York, à partir de la plateforme de données ouvertes du Department of City Planning. Le téléchargement était une archive ZIP, et nous avons placé le résultat de la décompression dans votre répertoire de travail.
Vous utiliserez la fonction dir() de R de base pour examiner le contenu de votre répertoire de travail, puis vous importerez le shapefile dans R.
Cette activité fait partie du cours
Visualiser des données géospatiales avec R
Instructions de l’exercice
- Utilisez
dir()sans argument pour trouver le nom du répertoire contenant le shapefile. - Utilisez
dir()en lui passant le chemin du répertoire du shapefile pour voir les fichiers qu'il contient. - Vous connaissez maintenant le répertoire et le nom du fichier. Utilisez
readOGR()pour lire le shapefile des quartiers dans un objet appeléneighborhoods. - Vérifiez le contenu en appelant
summary()surneighborhoods. - Vérifiez le contenu en traçant
neighborhoods.
Exercice interactif pratique
Essayez cet exercice en complétant ce code d’exemple.
library(sp)
library(rgdal)
# Use dir() to find directory name
# Call dir() with directory name
# Read in shapefile with readOGR(): neighborhoods
# summary() of neighborhoods
# Plot neighborhoods