Przycinanie rastra w R
Dzisiaj dowiecie się jak przyciąć raster do wektora. Skorzystamy w tym zadaniu z pakietów raster i rgdal. Załadujmy je (czytaj tutaj). Później wczytujemy raster (w przykładzie używamy zobrazowanie satelitarne Sentinel-2) do zmiennej r, który będziemy przycinać (czytaj tutaj). Wyświetlamy wczytany raster (czytaj tutaj):
Wczytujemy poligon do wycięcia do zmiennej border (czytaj tutaj) i wyświetlamy go na rastrze:
Przycinamy raster poligonem używając funkcji crop z pakietu raster:
|
1 |
r2 <- crop(r,border) |
Wyświetlamy wynik:
Raster został przycięty tylko do zasięgu poligonu. W celu wymaskowania pikseli znajdujących się poza poligonem trzeba użyć jeszcze funkcji mask:
|
1 |
r3 <- mask(r2,border) |
Wyświetlamy wynik:
W efekcie naszego działania otrzymaliśmy raster dokładnie docięty do naszego poligonu.
