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: pontosSpatialPointsDataFrame[1] distancias23x23 double matrix -- Julianna Trindade "Exalta-te, Senhor, na Tua força! Nós cantaremos e louvaremos o Teu poder." Salmos 21; 13
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
Se meu EPSG eh 4291(coordenadas geograficas), então a unidade é... ? Em 14 de agosto de 2012 17:33, Paulo Justiniano <paulojus@leg.ufpr.br>escreveu:
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
_______________________________________________ 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 (2)
-
Julianna Trindade -
Paulo Justiniano