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