Deu certo?<br>Fiquei curioso.<br><br>Aliás, não faz diferença no atual problema, mas como eu venho treinando o uso de funções da família apply para evitar loopings, resolvi mexer no código informado na primeira mensagem.<br>
<br>Troquei:<br>for (k in sq) {<br> for(j in sq){<br> draw.circle(k,j,radius=0.3,col="red")<br>}}<br><br>por dois sapply, um dentro do outro:<br>sapply(sq, function(k) sapply(sq,function(j) {draw.circle(k,j,radius=0.3,col="red")}))<br>
<br>O resultado foi o mesmo. Eu só não soube dizer se isso é mais eficiente do que como estava inicialmente, com os dois for.<br>Para bases grandes (talvez gigantes) os sapply reduziriam bastante o tempo de processamento ou de memória utilizada?<br>
<br>Abs.<br>Paulo<br><br><br><div class="gmail_quote">Em 16 de novembro de 2011 17:03, ASANTOS <span dir="ltr"><<a href="mailto:alexandresantosbr@yahoo.com.br">alexandresantosbr@yahoo.com.br</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<div bgcolor="#FFFFFF" text="#000000">
Prof. Paulo,<br>
<br>
Perfeito, as árvores tem aproximadamente o mesmo diâmetro
sim, pois são clones de eucalipto. Vou montar o algoritmo sugerido e
calcular as distâncias com a função outer().<br>
<br>
Novamente obrigado,<div class="im"><br>
<br>
<pre 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: <a href="tel:%2B55%2035%2092230304" value="+553592230304" target="_blank">+55 35 92230304</a></pre>
<br></div>
Em 16-11-2011 15:29, Paulo Justiniano escreveu:
<div><div></div><div class="h5"><blockquote type="cite">Alexandre
<br>
<br>
todas as arvores tem o mesmo " diametro" (area)
<br>
<br>
se positivo o caminho seria uma simulação dentro de um loop
"while"
<br>
até conseguir completar o numero de pontos requerido
<br>
<br>
o algoritmo seria
<br>
<br>
enquanto o numero de pontos nao é atingido:
<br>
<br>
1. sorteie 1 ponto
<br>
2. calcule a distancia deste 'as arvores presentes e tome a menor
<br>
3. verifique se ele está no raio da árvoce
<br>
4 se nao estiver guarde o ponto e aumente o valor do contador
<br>
caso contrario simule outro sem mudar o contador
<br>
<br>
<br>
<br>
<br>
On Wed, 16 Nov 2011, ASANTOS wrote:
<br>
<br>
<blockquote type="cite">Boa tarde pessoal,
<br>
<br>
Estou tentando criar um processo de pontos aleatório sobre
uma área regular representada por árvores regularmente
espaçadas, sendo que o processo de pontos só pode ocorrer onde
não existam árvores e estou quebrando a cabeça com isso, bom
primeiro criei a área com as árvores:
<br>
<br>
require(plotrix)
<br>
##Cria o limite da área
<br>
limx<-c(2,2,20.5,20.5)
<br>
limy<-c(2,20.5,20.5,2)
<br>
plot(limx,limy)
<br>
lim=cbind(limx,limy)
<br>
lim<-as.matrix(lim)
<br>
polygon(lim)
<br>
<br>
##Aqui criei a posição das árvores e com área igual ao seu
diâmetro
<br>
sq=seq(2.5,20,2.5)
<br>
for (k in sq) {
<br>
for(j in sq){
<br>
draw.circle(k,j,radius=0.3,col="red")
<br>
}}
<br>
<br>
<br>
Bom agora eu queria jogar sobre essa área o processo de pontos:
<br>
<br>
##Criando o processo de pontos
<br>
x <- runif(n=500,min=0, max=20)
<br>
y <- runif(n=500,min=0, max=20)
<br>
xy<-cbind(x,y)
<br>
points(xy)
<br>
##
<br>
<br>
Porém, os pontos também caem nas mesmas posições ocupadas pelas
árvores (obviamente), onde gostaria de saber se alguém
conheceria alguma solução para eu fazer com que não haja
sobreposição entre os raios que representam os diâmetros das
árvores e os pontos criados? Na verdade, preciso definir uma
maneira de fazer com que as coordenadas contidas nos diâmetros
das árvores não sejam validas no momento da geração do processo
de pontos, atualmente estou tentando soluções com o pacote
spatstat mais sem sucesso, pois não consigo fazer com que a área
criada seja um objeto ppp valido.
<br>
Obrigado,
<br>
<br>
-- <br>
Alexandre DOS SANTOS
<br>
Engenheiro Florestal, Msc.
<br>
Laboratório de Entomologia Florestal
<br>
Departamento de Entomologia
<br>
Universidade Federal de Lavras
<br>
Caixa Postal 3037
<br>
37200-000 - Lavras - Minas Gerais - Brasil
<br>
Tel: <a href="tel:%2B55%2035%2092230304" value="+553592230304" target="_blank">+55 35 92230304</a>
<br>
<br>
_______________________________________________
<br>
R-br mailing list
<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">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>
<br>
</blockquote>
<br>
<fieldset></fieldset>
<br>
<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>
</blockquote>
<br>
<br>
<pre 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: <a href="tel:%2B55%2035%2092230304" value="+553592230304" target="_blank">+55 35 92230304</a></pre>
</div></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>