Chciałbym narysować raster zawierający 4 różne wartości (1) z legendą tekstu jakościowego opisującą kategorie takie jak2 ale z kolorowymi pudełkami:

Próbowałem użyć legendy, takiej jak:

legend( 1,-20,legend = c("land","ocean/lake", "rivers","water bodies"))

ale nie wiem, jak skojarzyć jedną wartość z wyświetlanym kolorem. Czy istnieje sposób na odzyskanie koloru wyświetlanego z „działką” i użycie go w legendzie?