Načtení shapefile
Shapefile je jeden z nejrozšířenějších formátů pro sdílení prostorových dat a do R ho snadno načteš pomocí funkce readOGR() z balíčku rgdal. Funkce readOGR() má dva důležité argumenty: dsn a layer. To, co jim předáš, závisí na typu dat, která načítáš. Ve videu sis ukázal/a, že u shapefile souborů by dsn mělo být cesta k adresáři, který obsahuje soubory tvořící shapefile, a layer je název konkrétního shapefile (bez přípony).
Pro svou mapu budeš potřebovat hranice čtvrtí. Stáhli jsme tzv. Neighborhood Tabulation Areas, jak je definuje město New York, z Open Data Platform of the Department of City Planning. Stažený soubor byl ve formátu zip archivu a výsledek jeho rozbalení jsme umístili do tvého pracovního adresáře.
Pomocí funkce dir() ze základního R prozkoumáš obsah pracovního adresáře a pak shapefile načteš do R.
Toto cvičení je součástí kurzu
Vizualizace geoprostorových dat v R
Pokyny k cvičení
- Zavolej
dir()bez argumentů, abys zjistil/a název adresáře se shapefile. - Zavolej
dir()a předej mu cestu k adresáři shapefile, abys viděl/a jeho obsah. - Teď znáš název adresáře i souboru. Pomocí
readOGR()načti shapefile čtvrtí do objektuneighborhoods. - Obsah ověř tak, že zavoláš
summary()na objektuneighborhoods. - Obsah ověř také vykreslením objektu
neighborhoods.
Interaktivní cvičení na vyzkoušení si v praxi
Vyzkoušejte si toto cvičení dokončením tohoto ukázkového kódu.
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