Wczytywanie pliku rastrowego
Pliki rastrowe najłatwiej wczytać do R za pomocą funkcji raster() z pakietu raster. Wystarczy podać nazwę pliku (wraz z rozszerzeniem) jako pierwszy argument, x.
Funkcja raster() korzysta z wbudowanych funkcji pakietu raster do odczytu określonych typów plików (na podstawie rozszerzenia) – w pozostałych przypadkach przekazuje odczyt do funkcji readGDAL() z pakietu rgdal. Zaletą korzystania z raster() zamiast bezpośrednio z readGDAL() jest to, że zwraca ona obiekt klasy RasterLayer.
Jednym z najpopularniejszych formatów plików rastrowych jest GeoTIFF, z rozszerzeniem .tif lub .tiff. Pobraliśmy raster mediany dochodów ze spisu ludności US Census i umieściliśmy go w twoim katalogu roboczym.
Przyjrzyjmy się mu i wczytajmy go.
To ćwiczenie jest częścią kursu
Wizualizacja danych geoprzestrzennych w R
Instrukcje do ćwiczenia
- Użyj
dir(), aby sprawdzić zawartość katalogu roboczego. - Użyj
dir()ponownie, aby zajrzeć do katalogunyc_grid_data. - Użyj
raster(), aby wczytać raster mediany dochodów do zmiennejincome_grid, podając pełną ścieżkę do pliku.tif. - Użyj
summary(), aby sprawdzić, czy raster jest zapisany jako obiekt klasyRasterLayer. - Użyj
plot(), aby zweryfikować zawartość rastra.
Interaktywne ćwiczenie praktyczne
Spróbuj tego ćwiczenia, uzupełniając ten przykładowy kod.
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