Aan de slagGa gratis aan de slag

Een shapefile inlezen

Shapefiles zijn een van de meest gebruikte formaten om ruimtelijke gegevens te delen en kun je eenvoudig in R inlezen met readOGR() uit het rgdal-pakket. readOGR() heeft twee belangrijke argumenten: dsn en layer. Wat je precies aan deze argumenten meegeeft, hangt af van het type gegevens dat je inleest. In de video heb je geleerd dat voor shapefiles dsn het pad moet zijn naar de map die de bestanden bevat waaruit de shapefile bestaat, en dat layer de bestandsnaam is van de specifieke shapefile (zonder extensie).

Voor je kaart heb je wijkgrenzen nodig. We hebben de Neighborhood Tabulation Areas, zoals gedefinieerd door de City of New York, gedownload van het Open Data Platform van het Department of City Planning. De download was een zip-archief en we hebben het resultaat van het uitpakken in je werkmap gezet.

Je gebruikt de functie dir() uit base R om de inhoud van je werkmap te bekijken en leest daarna de shapefile in R in.

Deze oefening maakt deel uit van de cursus

Georuimtelijke data visualiseren in R

Cursus bekijken

Oefeninstructies

  • Gebruik dir() zonder argumenten om de naam van de map met de shapefile te vinden.
  • Gebruik dir() met het pad naar de shapefile-map om de bestanden erin te bekijken.
  • Je weet nu de map en de bestandsnaam. Gebruik readOGR() om de wijk-shapefile in te lezen in een object neighborhoods.
  • Controleer de inhoud door summary() op neighborhoods aan te roepen.
  • Controleer de inhoud door neighborhoods te plotten.

Praktische interactieve oefening

Probeer deze oefening eens door deze voorbeeldcode in te vullen.

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
Code bewerken en uitvoeren