Wyświetlanie danych rastrowych w R

Wczytajmy zobrazowanie z satelity Sentinel zawierające 4 kanały spektralne. Zapis pliku opisany został we wpisie Wczytanie i zapis rastra w R.

r = stack("F:/B_all.tif")

Wyświetl najpierw jeden kanał używając:

plot( r[[1]] )

Raster wyświetlony zostanie w domyślnej palecie barw w zakładce Plots w RStudio.

Wyświetlenie rastra w odcieniach szarości wymaga dodania atrybutu col do funkcji plot:

plot( r[[1]], col = grey.colors(255) )

grey.colors() to paleta odcieni szarości, a liczba 255 określa ich liczbę.

Wyświetlenie rastra w kompozycji barwnej umożliwia funkcja plotRGB(x, r, g, b, stretch), gdzie x – nazwa zmiennej z rastrem wielokanałowym, r – numer pojedynczego rastra dla kanału R, g – numer pojedynczego rastra dla kanału G, B – numer pojedynczego rastra dla kanału B, stretch – w przypadku obrazów o rozdzielczości radiometrycznej powyżej 8 bitów konieczne jest zdegradowanie tej rozdzielczości definiując sposób przepróbkowania). Po więcej informacji odsyłam do opisu funkcji plotRGB.

plotRGB(b_all, 3, 2, 1, “lin”)

Po wykonaniu funkcji wyświetli nam się kompozycja w kolorach prawdziwych: