<div dir="ltr"><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Não existe, ou não precisa existir. Gerar malhas pode ser feito com funções básicas. O mais específico da tarefa é verificar quais os pontos da malha que estão dentro do polígono de interesse. Os pacotes de estatística/visualização espacial geralmente trazem essas funções.<br><br><span style="font-family:courier new,monospace">## plot(xlim=c(0,100), ylim=c(0,100), x=NULL, y=NULL)<br>## pol <- locator(n=10, type="l")<br>## dput(pol)<br><br>pol <- <br>structure(list(x = c(6.01340996391505, 32.5768681613431, 71.4454577296386, <br>96.6416791080815, 96.0557204713735, 81.6020740992435, 62.0701195423111, <br>57.773089539786, 22.8108908828769, 9.91980087530154), y = c(13.011212611897, <br>9.73994454798204, 11.9207899239253, 44.6334705630748, 79.5269965781675, <br>95.2290832849592, 91.5216461458556, 36.5643426720846, 38.5271035104335, <br>24.5696931043964)), .Names = c("x", "y"))<br><br>plot(xlim=c(0,100), ylim=c(0,100), x=NULL, y=NULL)<br>polygon(pol, col="gray90")<br><br>xy <- expand.grid(x=seq(0,100,by=5), y=seq(0,100,by=5))<br><br>plot(xlim=c(0,100), ylim=c(0,100), x=NULL, y=NULL)<br>polygon(pol, col="gray90")<br>points(xy)<br><br>help(point.in.polygon, package="sp", help_type="html")<br>inside <- sp::point.in.polygon(xy$x, xy$y, pol$x, pol$y)<br><br>plot(xlim=c(0,100), ylim=c(0,100), x=NULL, y=NULL)<br>polygon(pol, col="gray90")<br>points(xy, col=inside)<br><br>xyin <- xy[as.logical(inside),]<br><br>plot(xlim=c(0,100), ylim=c(0,100), x=NULL, y=NULL)<br>polygon(pol, col="gray90")<br>points(xyin)<br><br></span>À disposição.<br>Walmes.<br><br></div>​</div>