Dear All,

I have created a grid for the State of FL using following command:

locpred = expand.grid(X=seq(-117026, 582149, 10000), Y=seq(2783590, 3443337, 10000))
plot(locpred)

The Gulf of Mexico is also included in the above grid.

I am using following command to import a shape file of FL:

library(maptools)
vector = readShapePoly("County_BND.shp")
plot(vector)

Now, I am looking for some code which can help me in obtaining exact grid of FL using above imported shape file as a clipping surface.

It will be great if somebody can guide me on the same as I am not getting any success in finding the code. Alternate procedure is also welcome!

Cheers,
Puneet Dwivedi


	[[alternative HTML version deleted]]

