數(shù)據(jù)來源GEE平臺

ERA5 DAILY

其中era5數(shù)據(jù)降水單位是M，轉(zhuǎn)換為mm可在gee平臺?1000，氣溫K轉(zhuǎn)換成攝氏度℃是?273.15

代碼隨后補上

降水的年度數(shù)據(jù)是月平均的累加

溫度的年度數(shù)據(jù)是一年12個月取平均

數(shù)據(jù)格式要求：年內(nèi)數(shù)據(jù)中要有相同的字符，比如1970年中的月數(shù)據(jù)都包含1970

library("raster")

tifdir <- paste("F:/data/",tifnames,sep="")

tifnames <- list.files(path = "F:/data", pattern = '.tif$')

tifdir#每個tif文件的路徑

tifnames#每個tif文件名

#思路，循環(huán)讀取每一年內(nèi)的柵格，柵格成塊進(jìn)行計算均值輸出

for (i in 1970:1972) {

? ? namefl=paste(i)

? ? allrt <-? lapply(tifdir[grep(pattern =namefl, tifnames)], raster)#

? ? all <-brick(allrt)#讀取成柵格塊，這里是包含i(如1970)的所有柵格

? ? m2<- calc(all,mean)#計算塊內(nèi)的均值，mean可以換成其他的比如sum或自定義

? ? csvfile=paste(i,".tif",sep="")# name

? ? flie=paste("F:/data/year/", csvfile ,sep = "")#path+name

? ? writeRaster(m2,flie)#寫出

? ? print(csvfile)

}

run完一個數(shù)據(jù)集要清除，不然run下一個數(shù)據(jù)集會有影響

注：run完的數(shù)據(jù)空間參考會清除，可以定義投影

參考來源：R語言柵格處理-月數(shù)據(jù)轉(zhuǎn)年均值 - 嗶哩嗶哩 (bilibili.com)