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

R如何为高度超过600米的每个多边形提取栅格的平均值？

1 个答案: