Citirea unui shapefile
Shapefile-urile sunt unul dintre cele mai comune formate pentru partajarea datelor spațiale și pot fi citite ușor în R folosind funcția readOGR() din pachetul rgdal. readOGR() are două argumente importante: dsn și layer. Ce anume transmiți acestor argumente depinde de tipul de date pe care îl citești. Ai aflat din video că, în cazul shapefile-urilor, dsn trebuie să fie calea către directorul care conține fișierele ce alcătuiesc shapefile-ul, iar layer este numele fișierului shapefile specific (fără extensie).
Pentru harta ta, ai nevoie de limitele cartierelor. Am descărcat Neighborhood Tabulation Areas, conform definiției Orașului New York, de pe platforma Open Data a Departamentului de Planificare Urbană. Descărcarea a fost sub forma unei arhive zip, iar rezultatul dezarhivării a fost plasat în directorul tău de lucru.
Vei folosi funcția dir() din R de bază pentru a examina conținutul directorului de lucru, apoi vei citi shapefile-ul în R.
Acest exercițiu face parte din cursul
Vizualizarea datelor geospațiale în R
Instrucțiuni pentru exercițiu
- Folosește
dir()fără argumente pentru a afla numele directorului shapefile-ului. - Folosește
dir(), transmițând calea către directorul shapefile-ului, pentru a vedea fișierele din interior. - Acum că știi directorul și numele fișierului, folosește
readOGR()pentru a citi shapefile-ul cartierelor într-un obiect numitneighborhoods. - Verifică conținutul apelând
summary()peneighborhoods. - Verifică conținutul reprezentând grafic
neighborhoods.
Exercițiu interactiv practic
Încearcă acest exercițiu completând acest cod de exemplu.
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