将多维NETCDF转换为r中的TIF

Question

我的.NC文件大约为651 MB，其中有几个数据集（每日）（RR_MRG_19810101_20181231_ENACT.NC）数据。我需要将（RR_MRG_198101011_20181231_ENACT.NC）转换为多个geotiff（每个时间切片，每月一个.tif）。同样，我想阅读时间序列。但是我发现

Error in .local(x, time, ...) : 
time must has the same length as the number of layers in RasterBrick][1]

这是我所做的，

library(raster)
library(zoo)
library(rts)
       
TRF = brick("rr_mrg_19810101_20181231_ENACT.nc")
crs(TRF) <- "+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 
+towgs84=0,0,0"
Awash_Extent<-c(37,44, 7,12)
Awash_E_resize<- crop(TRF,Awash_Extent)

Awash_Month<-seq(as.Date('1981-01-01'),as.Date('2018-12-31'),by = 
"month")

rt <- rts(Awash_E_resize, Awash_Month) 

write.rts(Awash_E_resize, filename='Awash_TRF_1981_2018_mon.tif', 
overwrite=TRUE)

您能帮我解决这个问题吗？

Answer 1

这样简单的事情可以

library(raster)
b <- brick("rr_mrg_19810101_20181231_ENACT.nc")
writeRaster(b, "timeslice.tif", bylayer=TRUE)

根据您的扩​​展问题来完成此操作：

使用 Terra （不需要其他软件包）

library(terra)    
TRF = rast("rr_mrg_19810101_20181231_ENACT.nc")
Awash_Extent<-c(37,44, 7,12)
Awash_E_resize<- crop(TRF,Awash_Extent)

我在这里创建类似的内容：

A <- rast(ext(37,44, 7,12), nlyr=365)
values(A) <- runif(size(A))
terra::time(A) <- seq(as.Date('1981-01-01'), as.Date('1981-12-31'), 1)

现在按月汇总：

m <- months(time(A))
f <- factor(m, levels=unique(m))

B <- tapp(A, m, fun=sum)
B
#class       : SpatRaster 
#dimensions  : 10, 10, 12  (nrow, ncol, nlyr)
#resolution  : 0.7, 0.5  (x, y)
#extent      : 37, 44, 7, 12  (xmin, xmax, ymin, ymax)
#coord. ref. : lon/lat WGS 84 
#source      : memory 
#names       :  January, February,    March,    April,      May,     June, ... 
#min values  : 12.49764, 10.47718, 11.80974, 11.29624, 11.01105, 10.83298, ... 
#max values  : 18.90536, 16.95743, 20.57114, 18.12099, 18.46543, 18.98500, ... 
 

您可以添加 filename = tapp 的参数，但是如果要将这些层保存为单独的文件，则可以使用 writeraster 。但是有了Terra，您需要自己提供文件名。

 fnames <- paste0("rain_", 1:nlyr(B), ".tif")
 writeRaster(B, fnames, overwrite=T)

（有一个关于file_ext的警告，您可以忽略）