<div dir="ltr">Alexandre, bom dia!<div><br></div><div>Não tenho experiência com a {spatstat}, mas tentei trabalhar seu código e acredito que possa lhe ser útil.<br><div><br></div><div><div><font face="courier new, monospace">### <BEGIN></font></div>
<div><font face="courier new, monospace">require(raster); require(spatstat); require(sp) </font></div><div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">## Raster </font></div>
<div><font face="courier new, monospace">----------------------------------------------------------- </font></div><div><font face="courier new, monospace">set.seed(20130622)</font></div><div><font face="courier new, monospace">r <- raster(nc=100, nr=100) </font></div>
<div><font face="courier new, monospace">r <- setValues(r, round(runif(ncell(r))* 255))</font></div><div><font face="courier new, monospace">plot(r)</font></div><div><font face="courier new, monospace"><br></font></div>
<div><font face="courier new, monospace">r.spdf <- data.frame(coordinates(r), px=getValues(r)); names(r.spdf)</font></div><div><font face="courier new, monospace">gridded(r.spdf) <- ~x+y</font></div><div><font face="courier new, monospace">spplot(r.spdf)</font></div>
<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">## Amostras</font></div><div><font face="courier new, monospace">------------------------------------------------- </font></div>
<div><font face="courier new, monospace">parc.p <- spsample(r.spdf, 15, type="random")</font></div><div><font face="courier new, monospace">parc.c <- vector("list",length(parc.p)) </font></div><div>
<font face="courier new, monospace">for(i in 1:length(parc.p)) {</font></div><div><font face="courier new, monospace"> parc.c[[i]] <- disc(radius=5, centre=c(parc.p$x[i], parc.p$y[i]))}</font></div><div><font face="courier new, monospace"><br>
</font></div><div><font face="courier new, monospace">W <- do.call(union.owin,parc.c)</font></div><div><font face="courier new, monospace">W.df <- as.data.frame(<a href="http://as.im">as.im</a>(W))</font></div><div>
<font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">W.sp <- W.df; coordinates(W.sp) <- ~x+y</font></div><div><font face="courier new, monospace">plot(W.sp)</font></div><div>
<font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">W.df$px <- over(W.sp, r.spdf)</font></div><div><font face="courier new, monospace">head(W.df)</font></div></div><div><font face="courier new, monospace">#</font></div>
<div><font face="courier new, monospace"><div># x y value px</div><div>#1 -179.9856 -39.83543 1 11</div><div>#2 -179.9856 -39.29907 1 76</div><div>#3 -179.9856 -38.76271 1 76</div><div>#4 -179.9856 -38.22635 1 76</div>
<div>#5 -179.9856 -37.68998 1 34</div><div>#6 -178.7263 -41.98087 1 157</div><div>#</div></font></div><div><font face="courier new, monospace">### <END></font></div></div></div><div class="gmail_extra"><br clear="all">
<div><div dir="ltr">Éder Comunello <<a href="mailto:comunello.eder@gmail.com" target="_blank">c</a><a href="mailto:omunello.eder@gmail.com" target="_blank">omunello.eder@gmail.com</a>> <br>Dourados, MS - [22 16.5'S, 54 49'W]<br>
</div></div>
<br><br><div class="gmail_quote">Em 26 de dezembro de 2013 08:37, ASANTOS <span dir="ltr"><<a href="mailto:alexandresantosbr@yahoo.com.br" target="_blank">alexandresantosbr@yahoo.com.br</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div text="#000000" bgcolor="#FFFFFF">
Bom dia Éder,<br>
<br>
Conforme sugerido, adicionei ao meu CRM:<br>
<br>
W <- do.call(union.owin,CIR)<br>
#<br>
imW<-<a href="http://as.im" target="_blank">as.im</a>(W)<br>
grids<-as(imW, "SpatialGridDataFrame")<br>
fullgrid(grids)=FALSE<br>
plot(grids)<br>
points(df.Im1)<br>
over(grids,df.Im1))##<br>
<br>
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<br>
<br>
Alexandre<br>
<br>
<br>
<div>Em 26/12/2013 05:53, Éder Comunello
escreveu:<br>
</div>
<blockquote type="cite"><div><div class="h5">Bom dia, Alexandre!<br>
<br>
Já deu uma olhada em sp::over?<br>
<br>
-- <br>
<div dir="ltr">Éder Comunello <<a href="mailto:comunello.eder@gmail.com" target="_blank">c</a><a href="mailto:omunello.eder@gmail.com" target="_blank">omunello.eder@gmail.com</a>> <br>
Dourados, MS - [22 16.5'S, 54 49'W]<br>
</div>
<br>
<br>
<fieldset></fieldset>
<br>
</div></div><div class="im"><pre>_______________________________________________
R-br mailing list
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</pre>
</div></blockquote><div class="im">
<br>
<pre cols="72">--
======================================================================
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: <a href="tel:%28%2B55%29%2065%208132-8112" value="+556581328112" target="_blank">(+55) 65 8132-8112</a> (TIM) <a href="tel:%28%2B55%29%2065%209686-6970" value="+556596866970" target="_blank">(+55) 65 9686-6970</a> (VIVO)
<a href="mailto:e-mails:alexandresantosbr@yahoo.com.br" target="_blank">e-mails:alexandresantosbr@yahoo.com.br</a>
<a href="mailto:alexandre.santos@cas.ifmt.edu.br" target="_blank">alexandre.santos@cas.ifmt.edu.br</a>
Lattes: <a href="http://lattes.cnpq.br/1360403201088680" target="_blank">http://lattes.cnpq.br/1360403201088680</a>
======================================================================
</pre>
</div></div>
<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br></div>