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

[Éder Comunello]

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]


Em 26 de dezembro de 2013 08:37, ASANTOS
<alexandresantosbr em yahoo.com.br>escreveu:

>  Bom dia Éder,
>
>   Conforme sugerido, adicionei  ao meu CRM:
>
> W <- do.call(union.owin,CIR)
> #
> imW<-as.im(W)
> grids<-as(imW, "SpatialGridDataFrame")
> fullgrid(grids)=FALSE
> plot(grids)
> points(df.Im1)
> over(grids,df.Im1))##
>
>     O over fez a listagem do valor dos pixels na imagem (stackIm1)
> contidos no interior das parcelas circulares CIR, mas meu objetivo é listar
> qual a coordenada e o valor de cada pixel em cada parcela circular CIR, com
> isso terei um resultado em colunas contendo: CIR, x, y e px
>
> Alexandre
>
>
> Em 26/12/2013 05:53, Éder Comunello escreveu:
>
> Bom dia, Alexandre!
>
> Já deu uma olhada em sp::over?
>
> --
> Éder Comunello <c <comunello.eder em gmail.com>omunello.eder em gmail.com>
> Dourados, MS - [22 16.5'S, 54 49'W]
>
>
>
> _______________________________________________
> R-br mailing listR-br em listas.c3sl.ufpr.brhttps://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
>
>
> --
> ======================================================================
> Alexandre dos Santos
> Proteção Florestal
> IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso
> Campus Cáceres
> Caixa Postal 244
> Avenida dos Ramires, s/n
> Bairro: Distrito Industrial
> Cáceres - MT                      CEP: 78.200-000
> Fone: (+55) 65 8132-8112 (TIM)   (+55) 65 9686-6970 (VIVO)e-mails:alexandresantosbr em yahoo.com.br
>         alexandre.santos em cas.ifmt.edu.br
> Lattes: http://lattes.cnpq.br/1360403201088680
> ======================================================================
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20131226/465c7c4f/attachment.html>