Läsa in en shapefil
Shapefiler är ett av de vanligaste formaten för att dela rumsliga data och kan enkelt läsas in i R med hjälp av readOGR() från paketet rgdal. readOGR() har två viktiga argument: dsn och layer. Vad du skickar till dessa argument beror på vilken typ av data du läser in. I videon lärde du dig att för shapefiler ska dsn vara sökvägen till den katalog som innehåller filerna som utgör shapefilen, och layer är filnamnet på den specifika shapefilen (utan filändelse).
För din karta behöver du stadsdelsgränser. Vi har laddat ned Neighborhood Tabulation Areas, definierade av New York City, från Open Data Platform of the Department of City Planning. Nedladdningen var ett zip-arkiv och vi har lagt resultatet av att packa upp filen i din arbetskatalog.
Du ska använda funktionen dir() från bas-R för att undersöka innehållet i din arbetskatalog och sedan läsa in shapefilen i R.
Den här övningen är en del av kursen
Visualisering av geospatial data i R
Övningsinstruktioner
- Använd
dir()utan argument för att ta reda på namnet på katalogen som innehåller shapefilen. - Använd
dir()och skicka in sökvägen till shapefile-katalogen för att se filerna inuti. - Nu känner du till katalogen och filnamnet. Använd
readOGR()för att läsa in stadsdels-shapefilen i ett objekt som heterneighborhoods. - Kontrollera innehållet genom att anropa
summary()påneighborhoods. - Kontrollera innehållet genom att plotta
neighborhoods.
Interaktiv övning med praktiskt arbete
Testa den här övningen genom att slutföra den här exempelkoden.
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