Question

我有两个文件：

一个名为m_stack的RasterStack（每层代表1个空气质量指数的模拟地图）;
一个名为shph的SpatialPointsDataFrame（其中每个点代表一个居住地）;

我想通过模拟提取每个居住地的平均空气质量指数。为了计算平均空气质量指数，我使用一个小缓冲区（60米）。

我正在设法为每个居住地提取平均空气质量指数。但我想通过模拟得到每个居住地的一个平均值。

接下来，我将尝试举一个简单的例子：

# create 2 raster maps (as example of 2 air quality index simulations)
map_r1 = raster(ncol = 10, nrow = 10, xmn = 0, xmx = 100, ymn = 0, ymx = 100)
values(map_r1) = seq(1:ncell(map_r1))
map_r2 = raster(ncol = 10, nrow = 10, xmn = 0, xmx = 100, ymn = 0, ymx = 100)
values(map_r2) = 1:ncell(map_r2)*2

# create RasterStack adding both raster maps
map_stack<-stack(map_r1,map_r2)

# create SpatialPoint (as example of 3 places of residence)
x <- c(20,40,60)
y <- c(20,40,60)
v1 <- c(1.0, 2.0, 3.0)
map_p<-as.data.frame(cbind(x,y,v1))
coordinates(map_p) <- ~x + y

然后，为了提取每个居住地的平均值，通过模拟我尝试了以下内容：

# extract mean value (buffer=15 as example) for each point
buff<-extract(map_stack,map_p,buffer=15)
mean<-sapply(buff,mean)

结果只给出了每个居住地的一个平均值（计算所有模拟的平均值）。我很乐意通过模拟了解如何提取每个住宅的平均空气质量指数。亲切的问候，曼努埃尔

Answer 1

这是一个小例子，向您展示这是如何工作的（摘自？摘录）。

r <- raster(ncol=36, nrow=18)
r[] <- 1:ncell(r)
s <- stack(r, r*2, r^2)

xy <- cbind(-50, seq(-80, 80, by=20))
e <- extract(s, xy[2:3,], buffer=1000000)
e

#[[1]]
#     layer.1 layer.2 layer.3
#[1,]     517    1034  267289
#[2,]     518    1036  268324
#[3,]     552    1104  304704
#[4,]     553    1106  305809
#[5,]     554    1108  306916
#[6,]     555    1110  308025
#
#[[2]]
#     layer.1 layer.2 layer.3
#[1,]     445     890  198025
#[2,]     446     892  198916
#[3,]     481     962  231361
#[4,]     482     964  232324

正如您所看到的，对于每个点，您都会得到一个矩阵，其中包含每个图层的值。您现在可以使用sapply来计算您想要的内容。如果你想要逐层的意思。你可以做到

sapply(e, colMeans)
            [,1]     [,2]
#layer.1    541.5    463.5
#layer.2   1083.0    927.0
#layer.3 293511.2 215156.5

通过模拟提取缓冲区（RasterStack）

1 个答案: