在R中以栅格格式从WorldClim创建生物气候变量NA值的掩码

时间:2016-12-16 15:39:35

标签: r mask raster

我从Bioclim下载了19个生物气候变量,我想为所有19个变量创建一个独特的掩码。掩码应排除19个变量中任何一个中具有NA值的像素。进一步的SDM分析需要此过程。我用堆栈做了这个并且工作了,但它不适用于Bioclim rasterBrick。 在R:

library(raster
bioclim<-raster::getData('worldclim',var='bio', res=10)
s<-sum(bioclim)
MASK[!is.na(s)] <- 1
plot(MASK)

1 个答案:

答案 0 :(得分:0)

坚持你的方法,你的代码在插入一行后实际定义了&#34; MASK&#34;对象

library(raster)
bioclim<-raster::getData('worldclim',var='bio', res=10)
s<-sum(bioclim)

# create MASK object with NA or other masking value
MASK <- setValues(s[[1]], NA)

MASK[!is.na(s)] <- 1
plot(MASK)

更简单一点,你可以使用:

MASK2 <- !is.na(s)
plot(MASK2)