Načtení rastrového souboru
Rastrové soubory nejsnáze načteš do R pomocí funkce raster() z balíčku raster. Stačí jako první argument x předat název souboru včetně přípony.
Funkce raster() využívá pro čtení určitých typů souborů (podle přípony) vlastní funkce balíčku raster, jinak předá načítání funkci readGDAL() z balíčku rgdal. Výhoda oproti přímému použití readGDAL() spočívá v tom, že raster() vrací objekt třídy RasterLayer.
Běžným typem rastrového souboru je GeoTIFF s příponou .tif nebo .tiff. Ze sčítání lidu USA jsme stáhli rastr s mediánovými příjmy a umístili ho do tvého pracovního adresáře.
Podívejme se na něj a načtěme ho.
Toto cvičení je součástí kurzu
Vizualizace geoprostorových dat v R
Pokyny k cvičení
- Pomocí
dir()se podívej, co je v tvém pracovním adresáři. - Pomocí
dir()se znovu podívej dovnitř složkynyc_grid_data. - Pomocí
raster()načti rastr mediánových příjmů do proměnnéincome_grid– předej jako argument celou cestu k souboru.tif. - Pomocí
summary()ověř, že je rastr uložen jako objekt třídyRasterLayer. - Pomocí
plot()ověř obsah rastru.
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(raster)
# Call dir()
# Call dir() on the directory
# Use raster() with file path: income_grid
# Call summary() on income_grid
# Call plot() on income_grid