Thanks to both of you. I was keeping the solution through rasterization as a fallback because the involved raster would be very large. I'll try to solve the problem keeping the vectorial format (at least for the polygons) using grass through spgrass6 and will report back. Also, maybe an intersection using joinPolys() could be another solution. A suboptimal solution would be a stratified random sample of points within the grid.
Agus -- View this message in context: http://r-sig-geo.2731867.n2.nabble.com/overlay-spatialgrid-on-spatial-polygons-data-frame-tp5079242p5082854.html Sent from the R-sig-geo mailing list archive at Nabble.com. _______________________________________________ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo