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+yspplot(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+yplot(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>
Em 26 de dezembro de 2013 08:37, ASANTOS <alexandresantosbr@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?
--
_______________________________________________ R-br mailing list R-br@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.
-- ====================================================================== 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@yahoo.com.br alexandre.santos@cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ======================================================================
_______________________________________________
R-br mailing list
R-br@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.
_______________________________________________ R-br mailing list R-br@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.
--
======================================================================
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@yahoo.com.br
alexandre.santos@cas.ifmt.edu.br
Lattes: http://lattes.cnpq.br/1360403201088680
======================================================================