Komine,
Simply use "add = TRUE" in the plot command.

>plot(testo2)
>plot(test1, add = TRUE)


On Wed, Mar 21, 2012 at 7:34 PM, Komine <momadou@yahoo.fr> wrote:

> Hi,
> Thank Robert and Tom2 for your help,
> I did what I wanted with rasterFromXYZ.
> Now I open my shapefile to put the contours of my raster, but the raster
> disappears. I used this code:
>
> # 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)
>
> How to open the raster and the boundary in the same window.
> Thank you in advance
>
> --
> View this message in context:
> http://r-sig-geo.2731867.n2.nabble.com/Grid-of-raster-increase-resolution-tp7389731p7391985.html
> Sent from the R-sig-geo mailing list archive at Nabble.com.
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>



-- 
--
Adam Sparks
GIS Specialist and Plant Pathologist
IRRI
DAPO Box 7777, Metro Manila, Philippines
ph: +63(2) 580-5600

	[[alternative HTML version deleted]]

