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/

Github: https://github.com/Squiercg

Lattes: http://lattes.cnpq.br/7355685961127056