Wczytywanie pliku shapefile
Pliki shapefile to jeden z najpopularniejszych formatów udostępniania danych przestrzennych. Można je łatwo wczytać do R za pomocą funkcji readOGR() z pakietu rgdal. Funkcja readOGR() przyjmuje dwa ważne argumenty: dsn i layer. To, co przekazujesz do tych argumentów, zależy od rodzaju wczytywanych danych. Z materiału wideo wiesz już, że w przypadku plików shapefile dsn powinien wskazywać ścieżkę do katalogu zawierającego pliki składające się na shapefile, natomiast layer to nazwa konkretnego pliku shapefile (bez rozszerzenia).
Do swojej mapy potrzebujesz granic dzielnic. Pobraliśmy dane Neighborhood Tabulation Areas, zgodnie z definicją Miasta Nowy Jork, z platformy Open Data Departamentu Planowania Miejskiego. Pobrane archiwum ZIP zostało rozpakowane i umieszczone w twoim katalogu roboczym.
Użyjesz funkcji dir() z bazowego R, aby sprawdzić zawartość katalogu roboczego, a następnie wczytasz plik shapefile do R.
To ćwiczenie jest częścią kursu
Wizualizacja danych geoprzestrzennych w R
Instrukcje do ćwiczenia
- Wywołaj
dir()bez argumentów, aby sprawdzić nazwę katalogu zawierającego plik shapefile. - Wywołaj
dir(), przekazując ścieżkę do katalogu shapefile, aby zobaczyć znajdujące się w nim pliki. - Teraz, gdy znasz już nazwę katalogu i pliku, użyj
readOGR(), aby wczytać shapefile z granicami dzielnic do obiektu o nazwieneighborhoods. - Sprawdź zawartość obiektu, wywołując na nim funkcję
summary(). - Sprawdź zawartość obiektu, wyświetlając
neighborhoodsna wykresie.
Interaktywne ćwiczenie praktyczne
Spróbuj tego ćwiczenia, uzupełniając ten przykładowy kod.
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