Me gustaría trazar un raster que contenga 4 valores diferentes (1) con una leyenda de texto categórica describiendo las categorías como 2 pero con cuadros de color.
He intentado usar la leyenda de la siguiente manera:
legend( 1,-20,legend = c("tierra","océano/lago", "ríos","cuerpos de agua"))
pero no sé cómo asociar un valor al color mostrado. ¿Hay alguna forma de recuperar el color mostrado con 'plot' y usarlo en la leyenda?
Raster inicial
Raster con leyenda