[R-br] Criar um spatial points e shapefile com coordenadas e seus atributos
ASANTOS
alexandresantosbr em yahoo.com.br
Quinta Dezembro 17 20:01:17 BRST 2015
Caros Listeiros,
Estou tentando criar um spatial points e posteriormente um
shapefile com quatro variáveis, longitude (long), latitude (lat),
espécie(esp) e cobertura (cob). Quando faço um cbind() apenas das
informações geográficas, tudo OK, mas quando tento inserir os atributos
esp e cob, a função SpatialPoints() não funciona. Alguém com experiência
poderia me dar uma luz? Segue CRM:
#Pacotes
require(shapefiles); require(sp)
#Pontos e variáveis esp e cob
long<-c(-54.8956314,-54.8986246,-54.8985713,
-54.8987823,-54.898895,-54.8988127)
lat<-c(-30.4011373,-30.4121832,-30.4121513,
-30.4120925,-30.4122421,-30.4123977)
esp<-c("a","a","b","a","b","b")
cob<-c(74.3,56.2,96.3,66.3,21.6,55.2)
final.pts<-cbind(long,lat,esp,cob)
#Transforma em spatial points
pts.sampling = SpatialPoints(final.pts, proj4string=CRS("+proj=longlat
+datum=WGS84"))
ANTcoor.df <- SpatialPointsDataFrame(pts.sampling,
data.frame(id=1:length(pts.sampling)))
#Salva em shapefile
writeOGR(ANTcoor.df, getwd(), 'ANTcoorF', 'ESRI Shapefile')
Obrigado,
--
======================================================================
Alexandre dos Santos
Proteção Florestal
IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso
Campus Cáceres
Caixa Postal 244
Avenida dos Ramires, s/n
Bairro: Distrito Industrial
Cáceres - MT CEP: 78.200-000
Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO)
e-mails:alexandresantosbr em yahoo.com.br
alexandre.santos em cas.ifmt.edu.br
Lattes: http://lattes.cnpq.br/1360403201088680
======================================================================
Mais detalhes sobre a lista de discussão R-br