Obrigado Alisson,
Era exatamente isto que eu precisava,
Redobrados agradecimentos,
Alexandre
On 06/01/2015 21:27, Alisson Lucrécio wrote:
Caro,
Se me permite, seria algo parecido com isso.
#Pacotes
require(spatstat)
require(sp)
## Gera um processo de pontos
xypois=rpoispp(25,win=owin(c(0,1),c(0,1)))## ~25 pontos aleatórios
xd=xypois$x
yd=xypois$y
cod<-seq(1,length(yd))
da<-cbind(xd, yd, cod) ## Crio um vetor com coordenada x e y e um
identificador
##Transformando em um objeto da classe espacial
pontos<- SpatialPoints(cbind(da[,1],da[,2]))
pontos<- SpatialPoints(list(da[,1],da[,2]))
pontos<- SpatialPoints(data.frame(da[,1],da[,2]))
##Cria a matriz de distâncias
res<-spDists(pontos)
res <- data.frame(res)
colnames(res) <- rownames(res)
pontos <- data.frame(pontos)
library(cluster)
clusters <- as.hclust(agnes(res, diss = T))
pontos$group <- cutree(clusters, h=0.1)
library(lattice)
xyplot(da...2. ~ da...1., group = group, data = pontos)
2015-01-06 21:29 GMT-02:00 ASANTOS <[email protected]
<mailto:[email protected]>>:
Caros listeiros,
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:
#Pacotes
require(spatstat)
require(sp)
## Gera um processo de pontos
xypois=rpoispp(25,win=owin(c(0,1),c(0,1)))## ~25 pontos aleatórios
xd=xypois$x
yd=xypois$y
cod<-seq(1,length(yd))
da<-cbind(xd, yd, cod) ## Crio um vetor com coordenada x e y e um
identificador
##Transformando em um objeto da classe espacial
pontos<- SpatialPoints(cbind(da[,1],da[,2]))
pontos<- SpatialPoints(list(da[,1],da[,2]))
pontos<- SpatialPoints(data.frame(da[,1],da[,2]))
##Cria a matriz de distâncias
res<-spDists(pontos)
res
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?
Obrigado,
--
======================================================================
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: (+55) 65 8132-8112 <tel:%28%2B55%29%2065%208132-8112> (TIM)
(+55) 65 9686-6970 <tel:%28%2B55%29%2065%209686-6970> (VIVO)
e-mails:[email protected]
<mailto:e-mails%[email protected]>
[email protected]
<mailto:[email protected]>
Lattes: http://lattes.cnpq.br/1360403201088680
======================================================================
---
Este email foi escaneado pelo Avast antivírus.
http://www.avast.com
_______________________________________________
R-br mailing list
[email protected] <mailto:[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.
--
Alisson Lucrecio da Costa
_______________________________________________
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.
--
======================================================================
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: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO)
e-mails:[email protected]
[email protected]
Lattes: http://lattes.cnpq.br/1360403201088680
======================================================================
---
Este email foi escaneado pelo Avast antivírus.
http://www.avast.com
_______________________________________________
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.