[R-br] Cálculo entre todos os pares de distâncias em um objeto espacial [RESOLVIDO]

[ASANTOS]

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 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.
>>
>
>
> _______________________________________________
> 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.


-- 
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

-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120604/cea85e58/attachment.html>