我使用R中的for循环从文件夹中读取netCDF文件,并提取给定经度、纬度列表的值。看起来很管用,除了一个问题。当循环按日期返回值时,它会在2月28之后创建1月29日至31日。这是我的R码:sb1 <- data.frame(longitude=1:10,latitude =1:10)
# Extractingzonal or sub-basin average rainfal
我试图从netcdf文件中提取光栅层名称,就像以前从光栅堆栈中写的那样。将光栅堆栈导出到ncdf可以工作得很好。("external/test.grd", package="raster"))r3<-r2*3
rstack<-stack(r1,r2,r3)
#Givea name - in this instance y