# 1读取2018年日最高气温数据
# 设置存储路径
setwd("D:/DATA/2017_min")#设置工作空间
library(raster)
library(terra)
flnames <- list.files(path = 'D:/DATA/2017_min/', pattern = '.tif$')
fl <- paste0("D:/DATA/2017_min/", flnames)
tmp2018 <- stack(fl)
names(tmp2018)
# 计算年平均气温
tmp2018_mean=stackApply(tmp2018,indices=1,fun= 'mean',na.rm=TRUE)
plot(tmp2018_mean)
terra::writeRaster(tmp2018_mean,paste0("D:\\DATA\\2017_min",names(2017),'.tif'), overwrite = T,filetype = "GTiff")
#View(tmp2018)
12-07
356