[R] Plotting multiple layers(maps) on same page
LCOG1
jroll at lcog.org
Tue Aug 3 18:54:52 CEST 2010
Hey guys and gals,
I searched through the forum and a bunch of R-mapping dedicated sites but
have not found what i know is quite elementary process, mapping more than
one layer on the same plot. I need to show some reference lines for the map
to make sense. I know the below wont work for anyone but to show what im
trying to do:
TazFile <- "data/input/TAZ.shp"
UGBFile<- "data/input/metugb.shp"
TazShape <- readShapeSpatial(TazFile)
UGBShape <- readShapeSpatial(UGBFile)
#Write out to pdf
pdf("Results/HhEmpForecast.pdf", width=8.2, height=11, onefile=TRUE)
#Households
plot(TazShape, col=colsHh[findInterval(TazShape$Hh, brksHh,
all.inside=TRUE)])
plot(UGBShape)
title(main="Households Forecast by Taz")
legend(4199277,860791.2,title="Total Households by
Taz",LegBrksHh,fill=colsHh,cex=.8)
#Employment
plot(TazShape, col=colsEmp[findInterval(TazShape$Emp, brksEmp,
all.inside=TRUE)])
title(main="Employment Forecast by Taz")
legend(4199277,860791.2,title="Total Employment by
Taz",LegBrksEmp,fill=colsEmp,cex=.8)
#Close file
dev.off()
I need Households and employment on separate maps but i need the UGB file to
be on each for reference. Doing it the above way obviously puts the
reference map (UGB) on a separate map. So how do i plot them on the same
page. I thought this was easy enough but im not finding a simple answer.
Thanks for the help.
In Solidarity,
JR
--
View this message in context: http://r.789695.n4.nabble.com/Plotting-multiple-layers-maps-on-same-page-tp2312223p2312223.html
Sent from the R help mailing list archive at Nabble.com.
More information about the R-help
mailing list