[R-br] Selecionar coordenadas e valor de pixels contidos no interior de parcelas

[Éder Comunello]

Mauro,

Não sei o que pode estar ocorrendo. Estou rodando sem problemas na v3.02.

Em quinta-feira, 26 de dezembro de 2013, Mauro Sznelwar<sznelwar em uol.com.br>
escreveu:
>
> Eu não consegui baixar a bibliotéca sptsatat, diz que não está disponível
na versão R 3.0.1 do Windows.
>
> Alexandre, bom dia!
> Não tenho experiência com a {spatstat}, mas tentei trabalhar seu código e
acredito que possa lhe ser útil.
>
> ### <BEGIN>
> require(raster); require(spatstat); require(sp)
> ## Raster
> -----------------------------------------------------------
> set.seed(20130622)
> r <- raster(nc=100, nr=100)
> r <- setValues(r, round(runif(ncell(r))* 255))
> plot(r)
> r.spdf <- data.frame(coordinates(r), px=getValues(r)); names(r.spdf)
> gridded(r.spdf) <- ~x+y
> spplot(r.spdf)
> ## Amostras
> -------------------------------------------------
> parc.p <- spsample(r.spdf, 15, type="random")
> parc.c <- vector("list",length(parc.p))
> for(i in 1:length(parc.p)) {
>           parc.c[[i]] <- disc(radius=5, centre=c(parc.p$x[i],
parc.p$y[i]))}
> W    <- do.call(union.owin,parc.c)
> W.df <- as.data.frame(as.im(W))
> W.sp <- W.df; coordinates(W.sp) <- ~x+y
> plot(W.sp)
> W.df$px <- over(W.sp, r.spdf)
> head(W.df)
> #
> #          x         y value  px
> #1 -179.9856 -39.83543     1  11
> #2 -179.9856 -39.29907     1  76
> #3 -179.9856 -38.76271     1  76
> #4 -179.9856 -38.22635     1  76
> #5 -179.9856 -37.68998     1  34
> #6 -178.7263 -41.98087     1 157
> #
> ### <END>
> .

-- 
Éder Comunello <c <comunello.eder em gmail.com>omunello.eder em gmail.com>
Dourados, MS - [22 16.5'S, 54 49'W]
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20131226/c60b50c7/attachment.html>