多边形内的网格点：以前的解决方案对我不起作用答案

【问题标题】：Grid points within polygon: previous solution does not work for me多边形内的网格点：以前的解决方案对我不起作用
【发布时间】：2018-07-09 19:53:11
【问题描述】：

由于某种原因，我无法让@RichPauloo 提供的解决方案正常工作，并感谢一些帮助。

我有一个名为“spdf”的 SpatialPolygonsDataFrame（在下面的 Dropbox 链接中） https://www.dropbox.com/s/ibhp5mqbgfmmntz/spdf.Rda?dl=0

我使用下面帖子中的代码来获取边界内的网格数据。

Create Grid in R for kriging in gstat

library(sp)
grd <- makegrid(spdf, n = 10000)

colnames(grd) <- c('x','y');

outline <- spdf@polygons[[1]]@Polygons[[1]]@coords

library(splancs)

new_grd <- grd[inout(grd,outline), ]

这是我得到的：

黑点是来自 makegrid 的“grd”
蓝点以“轮廓”为界
红点是“new-grd”作为边界内的网格

如您所见，它没有捕获边界内的所有数据？我究竟做错了什么？

【问题讨论】：

标签： r geospatial inout

【解决方案1】：

试试这个：

# packages
library(sp)

# make grid
grd <- makegrid(spdf, n = 100)
colnames(grd) <- c('x','y') # assign names to columns

# check the class
class(grd)

# transform into spatial points
grd_pts <- SpatialPoints(coords = grd, 
                         proj4string=CRS(as.character(NA)))

# check the class again
class(grd_pts)

# show that points fall outside of polygon
plot(spdf)
points(grd_pts)

# subset for points within the polygon
grd_pts_in <- grd_pts[spdf, ]

# visualize
plot(spdf)
points(grd_pts_in)

# transform grd_pts_in back into a matrix and data frame
gm <- coordinates(grd_pts_in) # matrix
gdf <- as.data.frame(coordinates(grd_pts_in)) # data frame

【讨论】：