Lire des données spatiales

Visualiser des données géospatiales avec R

Charlotte Wickham

Assistant Professor at Oregon State University

Revenus médians dans le comté de New York

Les secteurs de recensement ont un nombre de personnes similaire
Objets spatiaux :
- Polygones de secteurs de recensement
- Polygones de quartiers plus vastes
- Polygones de plans d'eau

Carte de NYC

Procédure

Lire des fichiers de formes décrivant quartiers et voies d'eau
Faire correspondre deux systèmes de coordonnées
Fusionner des données d'un data frame dans un SpatialPolygonsDataFrame
Peaufiner une carte prête pour la publication

Carte de NYC

Lecture d'un fichier de formes

Données vectorielles : décrites par des points, lignes, polygones
Le fichier de formes est le format le plus courant

Lecture d'un fichier de formes

# rgdal::readOGR() reads in vector formats
library(rgdal) 
library(sp)

dir()

"water"

dir("water")

"water-areas.dbf" "water-areas.prj" 
"water-areas.shp" "water-areas.shx"

water <- readOGR("water", "water-areas")

Source de données OGR avec pilote : ESRI Shapefile 
Source : "water", couche : "water-areas"
avec 20 entités
Elle contient 5 champs

Vérifier le résultat

summary(water)

Objet de classe SpatialPolygonsDataFrame
Coordonnées :
        min       max
x -74.04731 -73.90866
y  40.68419  40.88207
Projeté : FAUX 
...

plot(water)

Tracés de NYC

Lecture d'un fichier raster

library(rgdal)  # rgdal::readGDAL() reads in raster formats to sp objects
library(raster) # raster::raster() reads in raster formats to raster objects

dir()

"usgrid_data_2000"   "usgrid_data_2000_1"

dir("usgrid_data_2000")

"metadata"     "usarea00.tif"  "usba00.tif" "usfb00.tif" "usgrid-2000-variables.xls" "usp2500.tif" 
"uspop300.tif" "uspov00.tif"  "uspvp00.tif"

total_pop <- raster("usgrid_data_2000/uspop300.tif")

Vérifier le résultat

total_pop

classe       : RasterLayer 
dimensions   : 3120, 7080, 22089600  (nrow, ncol, ncell)
résolution   : 0.008333333, 0.008333333  (x, y)
extension    : -125, -66, 24, 50  (xmin, xmax, ymin, ymax)
réf. coord.  : +proj=longlat +datum=NAD83 +no_defs +ellps=GRS80 +towgs84=0,0,0 
source des données : /Users/wickhamc/Documents/Projects/courses-visualizing-geospatial-
              data-in-r/data/census_grids/usgrid_data_2000/uspop300.tif 
noms         : uspop300 
valeurs      : 0, 65535  (min, max)

Passons à la pratique !

Visualiser des données géospatiales avec R