R如何提取海拔高于600m的每个多边形栅格的平均值？

Question

首先，我设法提取平均光栅温度值为每个多边形，而下面的程序：

您可以下载这个链接上的GIS图层：

https://depots.univ-perp.fr/get?k=iTzEDSUkdyZVw2st78G

## load packages 

library(raster); library(rgdal) 

## Read rasters 

ras_temp<-raster("ras_temp.tif") 
plot(ras_temp) 
ras_alti<-raster("ras_alti.tif") 

## read polygon 

polygon <- readOGR(dsn = getwd(), layer = "polygon") 
plot(polygon,add=TRUE) 

## extract mean value for each polygon 

v1 <- extract(ras_temp, polygon, fun=mean, na.rm=TRUE) 
nom <- sapply(polygon@polygons, slot, "ID") 
v1 <- data.frame(ID = nom, Value = v1) 
View(v1) 

然后，我想提取每个多边形的温度平均值，但仅限于超过600米高度的表面。

不幸的是，我不能这样做，我的问题如何在我的函数“提取”中整合高度条件？

在此先感谢

Answer 1

您可以轻松地做到这一点是这样的：

# first resample the altitude raster to the temperature one so that they are 
# "aligned" 
ras_alti.new = resample(ras_alti, ras_temp, "bilinear") 

# set to NA all data in ras_temp corresponding to cells in ras_alti.new below 600 
# metre 
ras_temp.new = ras_temp 
ras_temp.new[ras_alti.new <= 600] = NA 

# extract the data 
v2 <- extract(ras_temp.new, polygon, fun=mean, na.rm=TRUE, sp = T) 
v2@data 

ID ras_temp 
0 417 64.11342 
1 433 68.53541