#########################################################
###############根据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 =&#