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:

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:

Wyświetlamy wynik:

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

Dodaj komentarz

Twój adres e-mail nie zostanie opublikowany.

Translate using Google Translate»
Social media & sharing icons powered by UltimatelySocial

Podoba Ci się nasza strona? Odwiedź nasz profil