Mas o que está acontecendo? Aparece um plot de mapa, sem pontos, ou você simplesmente não ve os pontos? Você tem alguma mensagem de erro? Tente dar um plot(dados_1) pra ver se está tudo certinho. Se os dois funcionam, confira no summary se os pontos estão dentro da área do mapa. Veja no exemplo abaixo que o que você faz era para funcionar, eu acho. library(maptools) xx <- readShapeSpatial(system.file("shapes/sids.shp", package="maptools")[1],IDvar="FIPSNO", proj4string=CRS("+proj=longlat +ellps=clrk66")) summary(xx) pontos<-data.frame(x=runif(10,-84.32385,-75.45698),y=runif(10,33.88199,36.58965)) summary(pontos) plot(xx) points(pontos,pch=19,col="red") Em 15 de setembro de 2014 00:14, Paulo Pimenta <paulopimenta6@hotmail.com> escreveu:
Prezados,
Estou tentando abrir dois shapefiles ao mesmo tempo, ou seja, sobrepor um sobre o outro e não consigo. Abaixo esta parte do meu código:
********************************************************************************************************
library(maptools)
#dados <- readShapeSpatial(file.choose())
dados <- readShapeSpatial("F:/qgis_1/exercicio_buffer_L1/estrada_2000.shp") dados_1 <- readShapePoints("F:/qgis_1/exercicio_buffer_L1/L1_SAD69UTMZ23.shp")
summary(dados) summary (dados_1)
plot(dados) points (dados_1)
********************************************************************************************************
Tentei usar o comando merge e juntar as duas shapes, mas ainda assim não dá certo. Por outro lado não sei como plotar os dois shapes. Agradeço desde já a ajuda!
Abraços.
_______________________________________________ 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.
-- Grato Augusto C. A. Ribas Site Pessoal: http://recologia.com.br/ <http://augustoribas.heliohost.org> Github: https://github.com/Squiercg Lattes: http://lattes.cnpq.br/7355685961127056