[R-sig-Geo] plot shape file using lattice xyplot
Zia Ahmed
zua3 at cornell.edu
Wed Dec 16 23:07:33 CET 2009
Hi,
I am trying to plot a shape file in R. I can do it easily using plot
or spplot function. But I want plot this map with lattice xyplot
function. Any one has any idea? I do not want use "maps" package.
Help will be appreciated...
Thanks
Zia
library(maptools)
library(lattice)
bound<-read.shape("Boundary.shp")
#OR
bound<-readShapePoly("Boundary.shp")
plot(bound)
spplot(bound)
# Or following way
bound<-Map2poly(bound)
xybd<-attr(bound,"maplim")
plot(xybd$x,xybd$y, type="n",asp=1, xlab="", ylab="")
for (i in 1:length(bound)) lines(bound[[i]]) # How do I apply this code
in following panel.function of xyplot?
windows(width=5, height=5)
xyplot(y~x, as.data.frame(tala),pch=1, cex=.8,col="black",
xlab=list("Easting (m)",cex=.9),ylab=list("Northing (m)",cex=.9),
aspect = "iso",
panel = function(...) {
panel.xyplot(...)
panel.abline(h = 0:4*5000 + 545000, v= 0:4*5000 + 2650000,
col = "light grey")
panel.points(coordinates(tala.nr),cex=.95,pch=20, col="black")
panel.text(2671000, 563000, cex=1, " (a)")
},
)
More information about the R-sig-Geo
mailing list