1

何らかの理由で、@RichPauloo が提供するソリューションを機能させることができず、助けていただければ幸いです。

「spdf」という名前の SpatialPolygonsDataFrame があります (下のドロップボックス リンクにあります) https://www.dropbox.com/s/ibhp5mqbgfmmntz/spdf.Rda?dl=0

以下の投稿のコードを使用して、境界内のグリッド データを取得しました。

gstat でのクリギング用に R でグリッドを作成する

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」です

  • 青い点は境界としての「輪郭」です

  • 赤い点は、境界内のグリッドとして「新しいグリッド」です

ご覧のとおり、境界内のすべてのデータをキャプチャしていませんか? 私は何を間違っていますか？

4

1 に答える 1

0

これを試して：

# 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
于 2018-07-10T00:08:29.460 に答える