根据栅格像元中的shapefile计算不同土地利用的面积

Question

数据集：Data

我正在处理0.1 x 0.1度规则网格上的空气污染数据的栅格数据集。我还有一个shapefile，其中包含不同土地用途的分类数据。

我正在尝试建立土地利用和空气污染之间的关系，为此，我需要计算R中每个栅格单元中不同的土地利用面积。

有人知道我是怎么做到的吗？

我已经在上面的链接中包含了shapefile和我正在使用的netcdf。

Answer 1

请不要提供需要下载的数据的链接。这些都会消失，然后这个问题对其他人就不再有用了。这也是一种麻烦，因为它使帮助你变得更加困难。相反，当问一个R问题时，你应该创建一个最小的自包含的，可重复的例子，就像R帮助文件中那样(这个网站上也有上千个例子)。

示例数据

library(terra)
f <- system.file("ex/elev.tif", package="terra")
pol <- rast(f)
names(pol) <- "pollution"

f <- system.file("ex/lux.shp", package="terra")
v <- vect(f)
v$landuse <- c("a", "b", "c")
lu <- aggregate(v, "landuse")

解决方案

e <- extract(pol, lu, fun=mean, na.rm=TRUE)

cbind(landuse=lu$landuse, e[,2,drop=F])
#  landuse pollution
#1       a  370.1575
#2       b  341.3008
#3       c  336.7545

然而，将土地使用作为多边形是很奇怪的，你会认为它是一个栅格；如果它是一个大型数据集，效果会更好。