我想从每个网格单元格中的栅格的随机坐标中提取非NA值。
栅格的一个例子
library(raster)
r <- raster(ncol = 10, nrow = 10, xmx = -80, xmn = -150, ymn = 20, ymx = 60)
values(r) <- runif(ncell(r))
网格示例
grid <- raster(extent(r))
res(grid) <- 15
proj4string(grid)<- proj4string(r)
gridpolygon <- rasterToPolygons(grid)
plot(r)
plot(gridpolygon, add = T)
如何为每个网格单元格内的每个栅格部分提取随机坐标值?
我真的很擅长这种东西,所以任何建议都会非常受欢迎。 感谢。
答案 0 :(得分:1)
你没有指定采样的所有条件,所以我在这里做了一些假设。 可以对每个网格多边形采样点并提取该值。以下是如何一次性完成并希望获得最佳效果的方法:
# pick random points per each grid cell and plot
set.seed(357)
pickpts <- sapply(gridpolygon@polygons, spsample, n = 1, type = "random")
sapply(pickpts, plot, add = TRUE)
# extract values of raster cells at specified points
sapply(pickpts, FUN = extract, x = r)
或者您可以循环并进行采样,直到获得非NA值。
N <- length(gridpolygon@polygons)
result <- rep(NA, times = N)
for (i in 1:N) {
message(sprintf("Trying polygon %d", i))
pl <- gridpolygon@polygons[[i]]
candval <- result[i] # start with NA
# sample until you get a non-NA hit
while (is.na(candval)) {
pickpoint <- spsample(pl, n = 1, type = "random")
candval <- extract(x = r, y = pickpoint)
}
result[i] <- candval
}
result
[1] 0.4235214 0.6081435 0.9126583 0.1710365 0.7788590 0.9413206 0.8589753
[8] 0.0376722 0.9662231 0.1421353 0.0804440 0.1969363 0.1519467 0.1398272
[15] 0.4783207