Quero plotar uma mapa unindo dois shapefiles diferentes: um da classe polígono (SpatialPolygonsDataFrame) e outro da classe linha (SpatialLinesDataFrame), mas não sei que função usar para uni-los.
# Ler arquivo shapefile da RDSM (websig - IDSM)
library(maptools)
# Toda a área
shp <- readShapeSpatial(file.choose())
crs.geo <- CRS("+proj=longlat +ellps=WGS84 +datum=WGS84")
proj4string(shp.t) <- crs.geo
par(bg='grey99', pty='s', oma=c(0.5,0.5,0.5,0.5))
plot(shp, col='lightsteelblue', axes=T)
title(main='Rio Tefé', font.main=1, cex.main=1.5)
shp1 <- readShapeSpatial(file.choose())
crs.geo <- CRS("+proj=longlat +ellps=WGS84 +datum=WGS84")
proj4string(shp1) <- crs.geo
par(bg='grey99', pty='s', oma=c(0.5,0.5,0.5,0.5))
plot(shp1, col='lightsteelblue', axes=T)
title(main='Lago Tefé', font.main=1, cex.main=1.5)