[R-br] Classificação de distâncias a partir de um raio em uma matriz de distância

[Alisson Lucrécio]

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 <alexandresantosbr em yahoo.com.br>:

> 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 (TIM)   (+55) 65 9686-6970 (VIVO)
> e-mails:alexandresantosbr em yahoo.com.br
>         alexandre.santos em cas.ifmt.edu.br
> Lattes: http://lattes.cnpq.br/1360403201088680
> ======================================================================
>
>
> ---
> Este email foi escaneado pelo Avast antivírus.
> http://www.avast.com
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> 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
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150106/8a78fc04/attachment.html>