[R-br] ajustar escala e interpolação

[Elias T. Krainski]

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))

-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20170918/272f1140/attachment.html>