¿Cómo etiquetar valores ráster en la trama?

Question

¿Cómo podría agregar valores de píxeles al gráfico?Puedo obtener los valores usando click() pero quiero que aparezca en la trama.

library(raster)
r <- raster(nrow=3, ncol=3)
r[] <- 1:ncell(r)
plot(r)
click(r) 

Answer 1

Intente lo siguiente, que se basa en piezas empedradas de la función devuelta por getMethod("click", signature="Raster").

myClick <- function(x, n = Inf, id = FALSE, xy = FALSE, cell = FALSE, 
                    type = "n", show = TRUE, ...) {
    i <- 0
    n <- max(n, 1)
    while (i < n) {
        i <- i + 1
        loc <- locator(1, type, ...)
        xyCoords <- cbind(x = loc$x, y = loc$y)
        cells <- na.omit(cellFromXY(x, xyCoords))
        if (length(cells) == 0)
              break
        value <- extract(x, cells)
        text(xyCoords, labels = value)
    }
}

## Try it out
myClick(r, n=4)

Answer 2

Si desea mostrar todos los valores, puede utilizar el text método:

library(raster)
r <- raster(nrow=3, ncol=3, vals=1:9)
plot(r)
text(r)

Para un subconjunto, puedes hacer algo como:

z <- rasterToPoints(r, function(x) x > 6 )
plot(r)
text(z[,1], z[,2], z[,3])

Answer 3

Sé que esta pregunta ya está marcada como respondida, pero construyendo en la solución de Josh y la pregunta de seguimiento de Eddie, aquí hay un poco para el bucle que hace lo que Eddie estaba pidiendo (parcela de valores de trama sin números decimales y sin usar click):

r <- raster(nrow=3, ncol=3)
r[] <- runif(ncell(r))

plot(r)

for(i in 1:ncell(r)){
 xycoords <- xyFromCell(r, cell = i)
 value <- extract(r, xycoords)
 text(xycoords, labels = round(value))
}