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:

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:

r3 <- mask(r2,border)

Wyświetlamy wynik:

W efekcie naszego działania otrzymaliśmy raster dokładnie docięty do naszego poligonu.