Hi, I have done the interpolation for my data and I was able to create the
contours in multipanel with the help of Pascal. Now, I want to clip the
contour with the shapefile. I want only the portion of contour to be
displayed which falls inside the boundary of the shapefile.

The data mydata.csv can be found on
https://www.dropbox.com/s/khi7nv0160hi68p/mydata.csv

The data for shapefile can be found on
https://www.dropbox.com/sh/ztvmibsslr9ocmc/YOtiwB8p9p

THe code I have used so far is as follows:

# Load Libraries
library(latticeExtra)
library(sp)
library(rgdal)
library(lattice)
library(gridExtra)

#Read Shapefile
hello <- readOGR("shape",
                 layer="Export_Output_4")
## Project the shapefile to the UTM 16 zone
proj4string(hello) <- CRS("+proj=utm +zone=16 +ellps=WGS84")

## Read Contour data
mydata <- read.csv("mydata.csv")
head(mydata )
summary(mydata)

#Create a contourplot
contourplot(Salinity ~ Eastings+Northings | Time, mydata,
cuts=30,pretty=TRUE)

Thank you so much. I would welcome any other ways to do this aside from
contourplot and lattice.

Best Regards,
Janesh

	[[alternative HTML version deleted]]

