首页 > 解决方案 > 在记录维度 (nco) 上连接 ncdf 文件和块

问题描述

我正在尝试将大量(~4000)个 ncdf 文件连接到一个文件中。每个输入文件都是一个空间栅格,具有 ax 和 y 维度。

我正在尝试使用 ncecat:

ncecat -4 -L 5 -D 2 --open_ram --cnk_csh=1000000000 \
   --cnk_dmn record,2000 --cnk_dmn x,10 --cnk_dmn y,10 \
   $input_files output.nc

这给了我这样的东西:

netcdf 测试 { 维度:记录 = UNLIMITED ;// (目前 6 个) y = 11250 ; x = 15000 ; 变量:float Band1(record, y, x) ; Band1:long_name = "GDAL 带号 1" ; Band1:_FillValue = -3.4e+38f ; Band1:grid_mapping = "transverse_mercator" ; Band1:_Storage = "chunked" ; Band1:_ChunkSizes = 1, 10, 10 ; Band1:_DeflateLevel = 5 Band1:_Filter = "|1,5 ; Band1:_Shuffle = "true" ; Band1:_Endianness = "little" ;

,并且记录维度实际上并没有被分块。我想我可以先运行这个命令,然后在输出文件上使用 ncks 来修复记录暗淡并再次重新分块,但是,由于 ncks 需要将所有内容读入 ram,而且也是另一个耗时的操作,我正在搜索告诉 ncecat 它还应该将 record-dim 视为分块维度的方法。我还没有找到一种方法来做到这一点。

标签: netcdfnetcdf4nco

解决方案


您的命令看起来格式正确,尽管我会发表一些评论。首先,您看到的行为可能是一个错误,因为该命令应生成大小为 2000 的记录维度块作为请求。其次,请阅读这里的分块文档。这导致添加--cnk_plc=cnk_xpl选项可能会有所帮助。第三,我建议你用ncrcatnot连接和分块文件ncecat。前者比后者占用更少的内存,如此所述。


推荐阅读