On 18/09/2017 20:01, Heloíse Pavanato wrote:
não achei a função nearest.pixel() no spatstat. Seria nearest.raster.point()?
Veja esse exemplo:


library(spatstat)

### simula uma imagem qualquer
x0 <- seq(0, 7, 0.1)
y0 <- seq(0, 10, 0.1)
grid.z <- im(t(outer(x0, y0, function(x,y)
    cos(x) - sin(y-2))), x0, y0)

### visualiza
plot(grid.z, x0, y0)
contour(x0, y0, t(grid.z$v), add=TRUE)

### simula alguns pontos
x.loc <- runif(5, 0, 7)
y.loc <- runif(5, 0, 10)

### extrai valores dos pixels proximos
id.r.c <- nearest.pixel(x.loc, y.loc, grid.z)
str(id.r.c)

### extrai os valores
z <- grid.z[Reduce('cbind', id.r.c)]
str(z)

###junta
cbind(x.loc, y.loc, z)

### adiciona como texto
text(x.loc, y.loc, format(z, dig=2))