Cálculo entre todos os pares de distâncias em um objeto espacial
Boa noite pessoal, Gostaria de calcular todas as distâncias de 6 coordenadas em UTM, à partir de um objeto da classe espacial, para isto fiz: ### Cálculo da distancia entre todos os pontos em UTM ninhos<-c("N01","N02","N03","N04","N05","N06") lat<-c(0503252,0503319,0503716,0503682,0503688,0503707) long<-c(7651950,7651694,7652672,7652670,7652627,7652662) ##Cálculo que compara a primeira coordenada com a segunda e assim sucessivamente .... lat.diff<-diff(lat) long.diff<-diff(long) d=(long.diff^2+lat.diff^2)^.5 d ##Porém não era bem isto que eu queria ainda ##Transformando em um objeto da classe espacial require(sp) pontos<- SpatialPoints(cbind(lat,long)) pontos<- SpatialPoints(list(lat,long)) pontos<- SpatialPoints(data.frame(lat,long)) Agora, gostaria de baseado no objeto espacial pontos, calcular todas as combinações de distâncias entre as coordenadas fornecidas. Alguém saberia me indicar uma função em algum pacote que faz isso, Obrigado, -- Alexandre dos Santos Engenheiro Florestal, Dr. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 9223-0304
Via library(help=sp) encontrei spDists() e spDistsN1()... E se isso nao fosse suficiente, tentaria: help.search("distance") b 2012/6/4 ASANTOS <alexandresantosbr@yahoo.com.br>:
Boa noite pessoal,
Gostaria de calcular todas as distâncias de 6 coordenadas em UTM, à partir de um objeto da classe espacial, para isto fiz:
### Cálculo da distancia entre todos os pontos em UTM ninhos<-c("N01","N02","N03","N04","N05","N06") lat<-c(0503252,0503319,0503716,0503682,0503688,0503707) long<-c(7651950,7651694,7652672,7652670,7652627,7652662)
##Cálculo que compara a primeira coordenada com a segunda e assim sucessivamente ....
lat.diff<-diff(lat) long.diff<-diff(long) d=(long.diff^2+lat.diff^2)^.5 d ##Porém não era bem isto que eu queria ainda
##Transformando em um objeto da classe espacial require(sp) pontos<- SpatialPoints(cbind(lat,long)) pontos<- SpatialPoints(list(lat,long)) pontos<- SpatialPoints(data.frame(lat,long))
Agora, gostaria de baseado no objeto espacial pontos, calcular todas as combinações de distâncias entre as coordenadas fornecidas. Alguém saberia me indicar uma função em algum pacote que faz isso,
Obrigado,
-- Alexandre dos Santos Engenheiro Florestal, Dr. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 9223-0304
_______________________________________________ R-br mailing list R-br@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.
voce deu uma olhada em dist() ? me parece que tem algo já apra classes espaciais no sp On Mon, 4 Jun 2012, ASANTOS wrote:
Boa noite pessoal,
Gostaria de calcular todas as distâncias de 6 coordenadas em UTM, à partir de um objeto da classe espacial, para isto fiz:
### Cálculo da distancia entre todos os pontos em UTM ninhos<-c("N01","N02","N03","N04","N05","N06") lat<-c(0503252,0503319,0503716,0503682,0503688,0503707) long<-c(7651950,7651694,7652672,7652670,7652627,7652662)
##Cálculo que compara a primeira coordenada com a segunda e assim sucessivamente ....
lat.diff<-diff(lat) long.diff<-diff(long) d=(long.diff^2+lat.diff^2)^.5 d ##Porém não era bem isto que eu queria ainda
##Transformando em um objeto da classe espacial require(sp) pontos<- SpatialPoints(cbind(lat,long)) pontos<- SpatialPoints(list(lat,long)) pontos<- SpatialPoints(data.frame(lat,long))
Agora, gostaria de baseado no objeto espacial pontos, calcular todas as combinações de distâncias entre as coordenadas fornecidas. Alguém saberia me indicar uma função em algum pacote que faz isso,
Obrigado,
-- Alexandre dos Santos Engenheiro Florestal, Dr. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 9223-0304
_______________________________________________ R-br mailing list R-br@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.
Obrigado Prof. Paulo e Benilton pelas sugestões, problema resolvido com função spDists(). Segue solução:
##Transformando em um objeto da classe espacial require(sp) Carregando pacotes exigidos: sp pontos<- SpatialPoints(cbind(lat,long)) pontos<- SpatialPoints(list(lat,long)) pontos<- SpatialPoints(data.frame(lat,long)) spDists(pontos) [,1] [,2] [,3] [,4] [,5] [,6] [1,] 0.0000 264.6224 858.24239 838.62983 805.24841 844.96686 [2,] 264.6224 0.0000 1055.50604 1041.31888 1003.31949 1042.86528 [3,] 858.2424 1055.5060 0.00000 34.05877 53.00000 13.45362 [4,] 838.6298 1041.3189 34.05877 0.00000 43.41659 26.24881 [5,] 805.2484 1003.3195 53.00000 43.41659 0.00000 39.82462 [6,] 844.9669 1042.8653 13.45362 26.24881 39.82462 0.00000
Em 04/06/2012 18:51, Paulo Justiniano escreveu:
voce deu uma olhada em dist() ?
me parece que tem algo já apra classes espaciais no sp
On Mon, 4 Jun 2012, ASANTOS wrote:
Boa noite pessoal,
Gostaria de calcular todas as distâncias de 6 coordenadas em UTM, à partir de um objeto da classe espacial, para isto fiz:
### Cálculo da distancia entre todos os pontos em UTM ninhos<-c("N01","N02","N03","N04","N05","N06") lat<-c(0503252,0503319,0503716,0503682,0503688,0503707) long<-c(7651950,7651694,7652672,7652670,7652627,7652662)
##Cálculo que compara a primeira coordenada com a segunda e assim sucessivamente ....
lat.diff<-diff(lat) long.diff<-diff(long) d=(long.diff^2+lat.diff^2)^.5 d ##Porém não era bem isto que eu queria ainda
##Transformando em um objeto da classe espacial require(sp) pontos<- SpatialPoints(cbind(lat,long)) pontos<- SpatialPoints(list(lat,long)) pontos<- SpatialPoints(data.frame(lat,long))
Agora, gostaria de baseado no objeto espacial pontos, calcular todas as combinações de distâncias entre as coordenadas fornecidas. Alguém saberia me indicar uma função em algum pacote que faz isso,
Obrigado,
-- Alexandre dos Santos Engenheiro Florestal, Dr. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 9223-0304
_______________________________________________ R-br mailing list R-br@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.
_______________________________________________ R-br mailing list R-br@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.
-- Alexandre dos Santos Engenheiro Florestal, Dr. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 9223-0304
Olá, como eu seto um nome de coluna para ser a referencia e aparecer o nome do objeto ao invés de 1,2,3,4,5 e 6 ???? Em 4 de junho de 2012 20:26, ASANTOS <alexandresantosbr@yahoo.com.br>escreveu:
Obrigado Prof. Paulo e Benilton pelas sugestões, problema resolvido com função spDists().
Segue solução:
##Transformando em um objeto da classe espacial require(sp) Carregando pacotes exigidos: sp pontos<- SpatialPoints(cbind(lat,long)) pontos<- SpatialPoints(list(lat,long)) pontos<- SpatialPoints(data.frame(lat,long)) spDists(pontos) [,1] [,2] [,3] [,4] [,5] [,6] [1,] 0.0000 264.6224 858.24239 838.62983 805.24841 844.96686 [2,] 264.6224 0.0000 1055.50604 1041.31888 1003.31949 1042.86528 [3,] 858.2424 1055.5060 0.00000 34.05877 53.00000 13.45362 [4,] 838.6298 1041.3189 34.05877 0.00000 43.41659 26.24881 [5,] 805.2484 1003.3195 53.00000 43.41659 0.00000 39.82462 [6,] 844.9669 1042.8653 13.45362 26.24881 39.82462 0.00000
Em 04/06/2012 18:51, Paulo Justiniano escreveu:
voce deu uma olhada em dist() ?
me parece que tem algo já apra classes espaciais no sp
On Mon, 4 Jun 2012, ASANTOS wrote:
Boa noite pessoal,
Gostaria de calcular todas as distâncias de 6 coordenadas em UTM, à partir de um objeto da classe espacial, para isto fiz:
### Cálculo da distancia entre todos os pontos em UTM ninhos<-c("N01","N02","N03","N04","N05","N06") lat<-c(0503252,0503319,0503716,0503682,0503688,0503707) long<-c(7651950,7651694,7652672,7652670,7652627,7652662)
##Cálculo que compara a primeira coordenada com a segunda e assim sucessivamente ....
lat.diff<-diff(lat) long.diff<-diff(long) d=(long.diff^2+lat.diff^2)^.5 d ##Porém não era bem isto que eu queria ainda
##Transformando em um objeto da classe espacial require(sp) pontos<- SpatialPoints(cbind(lat,long)) pontos<- SpatialPoints(list(lat,long)) pontos<- SpatialPoints(data.frame(lat,long))
Agora, gostaria de baseado no objeto espacial pontos, calcular todas as combinações de distâncias entre as coordenadas fornecidas. Alguém saberia me indicar uma função em algum pacote que faz isso,
Obrigado,
-- Alexandre dos Santos Engenheiro Florestal, Dr. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 9223-0304
_______________________________________________ R-br mailing list R-br@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.
_______________________________________________ R-br mailing listR-br@listas.c3sl.ufpr.brhttps://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 Engenheiro Florestal, Dr. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 9223-0304
_______________________________________________ R-br mailing list R-br@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.
-- Julianna Trindade "Exalta-te, Senhor, na Tua força! Nós cantaremos e louvaremos o Teu poder." Salmos 21; 13
participantes (4)
-
ASANTOS -
Benilton Carvalho -
Julianna Trindade -
Paulo Justiniano