if(require(splancs)){
data(parana)
gr <- expand.grid(seq(150, 800, by=20), seq(70, 500, by=20))
plot(gr, asp=1, pch="+")
polygon(parana$borders)
gr.in <- locations.inside(gr, parana$borders)
points(gr.in, col=2, pch="+")
}
Run the code above in your browser using DataLab