resposta parcial: a unudade de distancia depende da unidade de suas coordenadas. Por exemplo se usou UTM a distancia é em metros quanto as valores sao na mespa ordem que seus objetos adicione com rownames() e colnames() os nomes nas dimensoes usando a variavel que contem esta informacao algo como rownames(distancia) <- colnames(distancia) <- "NOME DA VARIAVEL" (MAS NOTE QUE VAI ESTAR USANDO ARMAZENAMENTO SEM NECESSÁRIAMENTE PRECISAR) On Tue, 14 Aug 2012, Julianna Trindade wrote:
Pessoal, tenho alguns pontos num shapefile e preciso saber as distâncias entre todos os pontos. Consegui fazer isso utilizando spDists mas a matriz que resulta não identifica os pontos pela coluna de atributo que eu preciso, apresenta apenas uma matriz quadrada utilizando 1,1; 1,2...2,1;2,2...e etc pra definir as linhas e as colunas e eu preciso que ao invés de 1,1 apareçam os atributos de uma coluna específica. Isso seria possível?? Além disso, qual unidade é utilizada nas distancias?Km?metros?
O codigo eh o seguinte:
rm(list=ls(all=TRUE)) library(base) library(maps) library(sp) library(maptools) library(spdep) library(classInt) library(RColorBrewer) library(shape) library(SDMTools) library(latticeExtra) pontos = readShapePoints("c:/users/Julianna/Documents/SecSaude/hospitais_produtividade.shp") proj4string(pontos) <- CRS("+init=epsg:4291") distancias = spDists(pontos) Resultando:
pontos SpatialPointsDataFrame[1] distancias 23x23 double matrix
-- Julianna Trindade
"Exalta-te, Senhor, na Tua força! Nós cantaremos e louvaremos o Teu poder."
Salmos 21; 13