[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