R统计每个区域使用栅格::分区

Question

我想计算区域内的分区统计数据，以符合等宽。

得到区段列表后，想法是分配块号，因为栅格::区段函数需要一个带有表示区段代码的栅格层。

当我尝试用分区编号填充范围时，填充的区域与范围不对应（请参见图）。为什么会这样？

library(raster)
library(foreach)

filename <- system.file("external/test.grd", package="raster")
r=raster(filename)

xmin = seq(178000, 181000, 1000)
ymin = seq(329000, 333000, 1000)

e=foreach(j=ymin,.combine=c) %:% 
  foreach(i=xmin) %do% {
    e=extent(i,i+1000,j,j+1000)
  }


# this should be going into the foreach loop
n=length(r[e[[1]]])
r[e[[1]]]=rep(1,n)

plot(r)
plot(e[[1]],add=T)

Answer 1

这就是你想要的吗？如果你不想修改你的渔网，我已经修改了光栅的范围。

library(raster)
library(foreach)

filename <- system.file("external/test.grd", package="raster")
r=raster(filename)
extent(r) <- c(178000,182000,329000,334000 )

xmin = seq(178000, 181000, 1000)
ymin = seq(329000, 333000, 1000)

e=foreach(j=ymin,.combine=c) %:% 
  foreach(i=xmin) %do% {
    e=extent(i,i+1000,j,j+1000)

    n=length(r[e])
    r[e]=rep((j-i),n)#zonal stat instead of j-i?
  }
plot(r);plot(raster(filename), add=T, alpha=0.5, legend=F)

Answer 2

您可以使用QGIS生成一个形状文件，其中包含作为常规栅格的栅格单元。之后，您可以将shapefile和rasterfile与光栅：：提取一起使用(https://www.rdocumentation.org/packages/raster/versions/3.0-7/topics/extract)或者使用exactextractr：：exact\u extract更快(https://cran.r-project.org/web/packages/exactextractr/index.html)