plot(test1, add = TRUE) But, that assumes the test1 object read from shapefile is using the same coordinate system as the raster, and that it shares at least some of the same extent. If they are in different coordinate systems (projection) one or the other must be transformed, if they don't overlap at all you need to set up the first plot to extend to the shared region. Use projection() to find out for raster object, and proj4string() for Spatial objects (sp and rgdal).
R-Sig-Geo is a more appropriate mailing list for questions on the raster, rgdal and sp packages, and please read the posting guide. Your code is not reproducible (so my answer is a big guess) and you did not declare your use of the raster package. Cheers, Mike. On Thu, Mar 22, 2012 at 7:49 AM, Komine <moma...@yahoo.fr> wrote: > Hi > I want to overlay my raster and its boundary which is a shapefile. > When I used thise code separately, all is ok: > # Open raster >>Image<-read.table("C:\\Users\\Documents\\Nouveau\\Frequence.txt",sep="",dec=",",header=TRUE) > >testo<-rasterFromXYZ(Image) >>plot(testo) >>testo2 <- aggregate(testo,fact=10, fun=mean) >>plot(testo2) > # open shapefiles= boundary >>library (rgdal) >>test1<-readOGR(dsn="C:\\Users\\Documents\\DISC D\\Nouveau",layer="pays") >>plot(test1) > > But How to overlay both layers (raster and shapefile) to have a same map?. > Thank you in advance > > -- > View this message in context: > http://r.789695.n4.nabble.com/To-overlay-my-raster-and-its-boundary-tp4493705p4493705.html > Sent from the R help mailing list archive at Nabble.com. > > ______________________________________________ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. -- Michael Sumner Institute for Marine and Antarctic Studies, University of Tasmania Hobart, Australia e-mail: mdsum...@gmail.com ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.