[R-sig-Geo] raster / polygon overlay

[laure velez]

Many thanks for the fast replies,



@ Roger : polygons and raster are ok (see attached file).

@  Johannes : I have already try this but the result was not satisfying, using a SpatialPointsDataFrame do not solve the problem as the points are not  inside the polygon (see attached file).

I guess the problem arises because of the raster resolution that is to large ....

Any suggestion would be appreciated,

Please find attached a plot showing the raster (image), the polygon and the centroid of the raster cells : 

image(grid, axes=TRUE)
plot(pol, add=TRUE)
points(coordinates(grid))

Regards,

Laure.
 		 	   		  
_________________________________________________________________
[[elided Hotmail spam]]

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20100105/2a958e2d/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ex.jpg
Type: image/jpeg
Size: 20164 bytes
Desc: not available
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20100105/2a958e2d/attachment.jpg>