
19 Set
2017
19 Set
'17
00:08
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))