"/ndvi3g_geo_v1_",value,"_","0712",".nc4")print(paste0("正在读取...%s_%s",value,"0712"))ndvi3g<-rasterizeGimms(x=fn1)mvc<-monthlyComposite(ndvi3g,indices=monthlyIndices(fn1))for(variablein
# downloadGimms() 可以使用这个函数来下载#library("gimms","rgdal")#这个包还有些函数,可以对数据进行精度控制,这里暂时用不到。for(valuein1982:2015){fn1<-paste0("你的输入路径/ndvi3g_geo_v1_",value,"_","0106",".nc4")ndvi3g<-rasterizeGimms(x=fn1)sprintf("正在读取...%s_%s",value,"...
lon=ncread('F:\Global\NDVI3g\ndvi3g_geo_v1_2015_0712.nc4','lon');lat=ncread('F:\Global\NDVI3g\ndvi3g_geo_v1_2015_0712.nc4','lat');%tmn=ncread('C:\Users\lenovo\Documents\Tencent Files\1079192939\FileRecv\月最低气温数据集_2000_2002\tmn_2000_2002.nc','tmn');lonlim_left=-1...