Точки сетки в пределах многоугольника: предыдущее решение не работает для меня

Question

По какой-то причине я не могу заставить работать решение, предоставленное @RichPauloo, и очень ценю некоторую помощь.

У меня есть SpatialPolygonsDataFrame, который называется "spdf" (в приведенной ниже ссылке на ссылку) https://www.dropbox.com/s/ibhp5mqbgfmmntz/spdf.Rda?dl=0

Я использовал код из поста, чтобы получить данные сетки в пределах границы.

Создать Grid в R для кригинга в 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), ]

Вот что я получаю:

• Черные точки "grd" от makegrid

• Синие точки являются "контуром" как граница

• Красные точки "New-Grd" как сетка в пределах границы

Как вы видите, он не захватывает все данные в пределах границы? Что я делаю неправильно?

Answer 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