根据掩膜求不同区域平均值

#########################################################

###############根据mask求不同区域年平均值################

#########################################################

#数据选取

library(sp)

library(raster)

library(rgdal)

mask_NA<- raster("E:\\Zhangbh\\ali ndvi\\na.tif")

mask<- raster("E:\\Zhangbh\\ali ndvi\\xianjie2.tif")#master

setwd("E:\\Zhangbh\\ali ndvi2")

list1 <- list.files("E:\\Zhangbh\\ali ndvi2",pattern=".tif")

#mask修订

mask1<-mask

mask1[mask1==0]<- 8

#mask_NA[mask_NA！=NA]<-1

mask_NA<-mask_NA/mask_NA

plot(mask_NA)

index<- c()

#修改mask个数

for(j in 1:length(list1)){

NDVI <- raster(list1[j])

test <- c()

for(i in 1:8){

mask2 <- mask1

mask2[mask2 != i] <- NA

mask2[mask2 =&#