[R-br] união 2 shapefiles para mapa
Heloíse Pavanato
helopavanato em gmail.com
Quarta Junho 18 15:08:37 BRT 2014
Boa tarde,
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.
O CMR com os mapas plotados em separado está abaixo:
# 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)
Os arquivos Rio_Tefe.shp e Lago_tefe que estou usando estão aqui:
http://www.datafilehost.com/d/4099211a
A quem puder me ajudar, obrigada!
Abraço,
Heloise
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140618/a624ab37/attachment.html>
Mais detalhes sobre a lista de discussão R-br