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.
|
1 |
r = stack("F:/B_all.tif") |
Wyświetl najpierw jeden kanał używając:
|
1 |
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:
|
1 |
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.
|
1 |
plotRGB(b_all, 3, 2, 1, “lin”) |
Po wykonaniu funkcji wyświetli nam się kompozycja w kolorach prawdziwych:
