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]>: > 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:[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. > -- 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.
