Boa tarde pessoal,
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:
require(plotrix)
##Cria o limite da área
limx<-c(2,2,20.5,20.5)
limy<-c(2,20.5,20.5,2)
plot(limx,limy)
lim=cbind(limx,limy)
lim<-as.matrix(lim)
polygon(lim)
##Aqui criei a posição das árvores e com área igual ao seu diâmetro
sq=seq(2.5,20,2.5)
for (k in sq) {
for(j in sq){
draw.circle(k,j,radius=0.3,col="red")
}}
Bom agora eu queria jogar sobre essa área o processo de pontos:
##Criando o processo de pontos
x <- runif(n=500,min=0, max=20)
y <- runif(n=500,min=0, max=20)
xy<-cbind(x,y)
points(xy)
##
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.
Obrigado,
--
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
_______________________________________________
R-br mailing list
[email protected]
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.
