<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
Obrigado Prof. Paulo,<br>
<br>
Vou converter para objetos Spatial do sp e tentar a função
overlay, estava lendo a função é bem isso que preciso. Já tinha
tentando <span class="aett-12">quadrat.count() do Rcitrus</span> e
quadratcount() do spatstat, mais o problema é que estou começando o
estudo com simulações e as janelas são quadradas, porém quando for
para campo e como as áreas possuem contornos irregulares, já estou
prevendo problemas com as inúmeras parcelas que se localizam na
borda e não vão formar quadrats de mesmo tamanho, por isso preciso
estudar abordagem prevendo isso.<br>
<br>
Redobrados agradecimentos,<br>
<br>
Alexandre<br>
<br>
<pre class="moz-signature" cols="72">--
Alexandre DOS SANTOS
Engenheiro Florestal, Msc.
Laboratório de Entomologia Florestal
Departamento de Entomologia
Universidade Federal de Lavras
Caixa Postal 3037
37200-000 - Lavras - Minas Gerais - Brasil
Tel: +55 35 92230304</pre>
<br>
Em 16-11-2011 16:34, Leonard de Assis escreveu:
<blockquote cite="mid:4EC40252.4090209@gmail.com" type="cite">
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
Prof Paulo<br>
<br>
Aproveitando o email, já que tenho uma dúvida no mesmo assunto:<br>
<br>
1. Qual seria uma bibliografia sugerida para se entender essa
abordagem "geoestatística" da análise de dados<br>
2. Algum curso (presencial ou EAD) previsto para ano que vem?<br>
<br>
Pergunto isso pois tenho muitas idéias de aplicação destes
conceitos, mas me falta referência (e até incentivo) para
pesquisar mais sobre o assunto e formar alguma opinião.<br>
<pre class="moz-signature" cols="72">[]s
Leonard de Assis
assis <dot> leonard <at> gmail <dot> com
</pre>
<br>
Em 16/11/2011 15:25, Paulo Justiniano escreveu:
<blockquote
cite="mid:alpine.DEB.2.00.1111161522120.3228@pataxo.est.ufpr.br"
type="cite">Alexandre <br>
<br>
sem fornecer codigo especifico aqui vao algums ideias de por
ojnde eu comecaria <br>
<br>
algumas alternativas caso nao queira programar isto: <br>
<br>
1. represetnar os dados no formato sp. Um como SpatialPoints <br>
e outro como SpatialGrid <br>
Depois a operaçõa de overlay (over() ) com tabloe vai te
fornecer as contagens <br>
<br>
2. <br>
este problema é o mesmo do chamado "quadrat counts" <br>
que é tratado em pacotes como spatstat (oficial) <br>
ou Rcitrus (nao oficial) <br>
<br>
aqui vai um exemplinho usando o 1o <br>
<br>
require(spatstat) <br>
## Carregando um conjunto de dados <br>
data(swedishpines) <br>
X <- swedishpines <br>
plot(X) <br>
summary(X) <br>
## Contagem de Quadrats <br>
Q <- quadratcount(X, nx = 4, ny = 3) <br>
Q <br>
plot(X) <br>
plot(Q, add = TRUE, cex = 2) <br>
<br>
<br>
<br>
On Tue, 15 Nov 2011, ASANTOS wrote: <br>
<br>
<blockquote type="cite">Boa tarde pessoal, <br>
<br>
Estou tentando calcular a frequência de pontos dentro de
parcelas selecionadas em uma área através da função table() e
estou tendo dificuldades, <br>
bom consegui calcular para todas as parcelas da área mais não
estou conseguindo utilizar o comando para selecionar algumas
parcelas de interesse fiz: <br>
<br>
require(splancs) <br>
#Geração dos pontos aleatórios <br>
x <- runif(n=500,min=0, max=20) <br>
y <- runif(n=500,min=0, max=20) <br>
<br>
plot(x=x,y=y,axes=F, xlab="", ylab="", <br>
xlim=c(-1,20),ylim=c(0,21))##Plota os pontos <br>
<br>
##Divide em parcelas <br>
segments(x0=rep(0,20),y0=0:20, <br>
x1=rep(20,20),y1=0:20, col="gray", lty=2) <br>
segments(x0=0:20,y0=rep(0,20), <br>
x1=0:20,y1=rep(20,20),col="gray", lty=2) <br>
<br>
pt0 <- c(0.5,0.5)##Define o centro da primeira parcela <br>
<br>
N_S=seq(pt0[1], by=1, length.out=20)##Define o centro das
parcelas <br>
L_O=seq(pt0[2],by=1, length.out=20) <br>
grid0 <- expand.grid(N_S=seq(pt0[1], by=1, length.out=20),
<br>
L_O=seq(pt0[2],by=1, length.out=20)) <br>
<br>
##Cria o limite da área <br>
limx<-c(2.5,2.5,18,18) <br>
limy<-c(2.5,18,18,2.5) <br>
lim=cbind(limx,limy) <br>
lim<-as.matrix(lim) <br>
polygon(lim)##Desenha o limite <br>
<br>
grid0<-as.matrix(grid0) <br>
m.int<-pip(grid0,lim,out=FALSE)##Seleciona as parcelas
dentro do limite <br>
points(m.int,col="blue", pch=10) <br>
<br>
##Frequencia em cada parcela <br>
x.parcela <- cut(x,breaks=0:20) <br>
y.parcela <- cut(y,breaks=0:20) <br>
contag <- data.frame(table(x.parcela,y.parcela)) <br>
<br>
Aqui esta minha dificuldade, se uso o cut(), só consigo
dividir a área toda, porém não estou conseguindo selecionar
as parcelas localizadas dentro do limite, <br>
alguém poderia me dar uma luz, <br>
Obrigado, <br>
<br>
-- <br>
Alexandre DOS SANTOS <br>
Engenheiro Florestal, Msc. <br>
Laboratório de Pragas Florestais <br>
Departamento de Entomologia <br>
Universidade Federal de Lavras <br>
Caixa Postal 3037 <br>
37200-000 - Lavras - Minas Gerais - Brasil <br>
Tel: +55 35 92230304 <br>
<br>
_______________________________________________ <br>
R-br mailing list <br>
<a moz-do-not-send="true" class="moz-txt-link-abbreviated"
href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a>
<br>
<a moz-do-not-send="true" class="moz-txt-link-freetext"
href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
<br>
Leia o guia de postagem (<a moz-do-not-send="true"
class="moz-txt-link-freetext"
href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>)
e forneça código mínimo reproduzível. <br>
<br>
</blockquote>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
R-br mailing list
<a moz-do-not-send="true" class="moz-txt-link-abbreviated" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a>
<a moz-do-not-send="true" class="moz-txt-link-freetext" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a moz-do-not-send="true" class="moz-txt-link-freetext" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</pre>
</blockquote>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
R-br mailing list
<a class="moz-txt-link-abbreviated" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a>
<a class="moz-txt-link-freetext" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a class="moz-txt-link-freetext" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</pre>
</blockquote>
<br>
<br>
<pre class="moz-signature" cols="72">
</pre>
</body>
</html>