<div dir="ltr"><div>Caro,</div><div><br></div><div>Se me permite, seria algo parecido com isso.</div><div><br></div><div>#Pacotes<br></div><div>require(spatstat)</div><div>require(sp)</div><div><br></div><div>## Gera um processo de pontos</div><div>xypois=rpoispp(25,win=owin(c(0,1),c(0,1)))## ~25 pontos aleatórios</div><div>xd=xypois$x</div><div>yd=xypois$y</div><div>cod<-seq(1,length(yd))</div><div>da<-cbind(xd, yd, cod) ## Crio um vetor com coordenada x e y e um identificador</div><div><br></div><div>##Transformando em um objeto da classe espacial</div><div>pontos<- SpatialPoints(cbind(da[,1],da[,2]))</div><div>pontos<- SpatialPoints(list(da[,1],da[,2]))</div><div>pontos<- SpatialPoints(data.frame(da[,1],da[,2]))</div><div><br></div><div>##Cria a matriz de distâncias</div><div>res<-spDists(pontos)</div><div>res <- data.frame(res)</div><div>colnames(res) <- rownames(res)</div><div><br></div><div>pontos <- data.frame(pontos)</div><div><br></div><div>library(cluster)</div><div>clusters <- as.hclust(agnes(res, diss = T))</div><div>pontos$group <- cutree(clusters, h=0.1)</div><div><br></div><div>library(lattice)</div><div><br></div><div>xyplot(da...2. ~ da...1., group = group, data = pontos)</div><div class="gmail_extra"><br><div class="gmail_quote">2015-01-06 21:29 GMT-02:00 ASANTOS <span dir="ltr"><<a href="mailto:alexandresantosbr@yahoo.com.br" target="_blank">alexandresantosbr@yahoo.com.br</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Caros listeiros,<br>
<br>
Tenho um objeto com coordenadas e um identificador (da) que representa um conjunto de pontos. Calculei a matriz de distância entre todos os pontos e gostaria de a partir de um raio arbitrário, criar uma nova coluna de classificação dos pontos (em da) que estão contidos nesta distância, ou seja, quais pontos estão próximos uns dos outros em uma distância escolhida, segue CRM:<br>
<br>
#Pacotes<br>
require(spatstat)<br>
require(sp)<br>
<br>
## Gera um processo de pontos<br>
xypois=rpoispp(25,win=owin(c(<u></u>0,1),c(0,1)))## ~25 pontos aleatórios<br>
xd=xypois$x<br>
yd=xypois$y<br>
cod<-seq(1,length(yd))<br>
da<-cbind(xd, yd, cod) ## Crio um vetor com coordenada x e y e um identificador<br>
<br>
##Transformando em um objeto da classe espacial<br>
pontos<- SpatialPoints(cbind(da[,1],da[<u></u>,2]))<br>
pontos<- SpatialPoints(list(da[,1],da[,<u></u>2]))<br>
pontos<- SpatialPoints(data.frame(da[,<u></u>1],da[,2]))<br>
<br>
##Cria a matriz de distâncias<br>
res<-spDists(pontos)<br>
res<br>
<br>
<br>
Aqui eu calculei a matriz de distância, mas agora preciso de alguma função que à partir de um raio, como R=0.1, faça a classificação que eu preciso, alguém teria alguma sugestão ou conhecimento de alguma função que faça isto?<br>
<br>
Obrigado,<br>
<br>
<br>
-- <br>
==============================<u></u>==============================<u></u>==========<br>
Alexandre dos Santos<br>
Proteção Florestal<br>
IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso<br>
Campus Cáceres<br>
Caixa Postal 244<br>
Avenida dos Ramires, s/n<br>
Bairro: Distrito Industrial<br>
Cáceres - MT CEP: 78.200-000<br>
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)<br>
<a href="mailto:e-mails%3Aalexandresantosbr@yahoo.com.br" target="_blank">e-mails:alexandresantosbr@<u></u>yahoo.com.br</a><br>
<a href="mailto:alexandre.santos@cas.ifmt.edu.br" target="_blank">alexandre.santos@cas.ifmt.edu.<u></u>br</a><br>
Lattes: <a href="http://lattes.cnpq.br/1360403201088680" target="_blank">http://lattes.cnpq.br/<u></u>1360403201088680</a><br>
==============================<u></u>==============================<u></u>==========<br>
<br>
<br>
---<br>
Este email foi escaneado pelo Avast antivírus.<br>
<a href="http://www.avast.com" target="_blank">http://www.avast.com</a><br>
<br>
______________________________<u></u>_________________<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/<u></u>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-<u></u>guia</a>) e forneça código mínimo reproduzível.<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr"><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;display:inline!important;float:none;background-color:rgb(255,255,255)">Alisson Lucrecio da Costa</span></div></div>
</div></div>