从netCDF中的组中提取变量

时间:2016-12-05 19:10:51

标签: python nco netcdf4

我有一个带有foll的netCDF文件。 ncdump:

netcdf test_nc {
dimensions:
    time = UNLIMITED ; // (20 currently)
    latitude = 360 ;
    longitude = 720 ;
    N = 3 ;
    strlen = 1 ;
variables:
    float data_array(time, latitude, longitude, N) ;
        data_array:_FillValue = -9999.f ;
        data_array:units = "1" ;
        data_array:long_name = "data_array" ;
    char N(N, strlen) ;
    double latitude(latitude) ;
        latitude:standard_name = "latitude" ;
        latitude:units = "degrees_north" ;
    double longitude(longitude) ;
        longitude:standard_name = "longitude" ;
        longitude:units = "degrees_east" ;
    double time(time) ;
        time:standard_name = "time" ;
        time:units = "days since 2000-01-01 00:00:00.0" ;
        time:calendar = "gregorian" ;
}

如何提取data_array组中的最后一个变量?如果我使用ncks,我可以像这样提取所有data_array

ncks -v data_array test_nc.nc output_nc.nc

但是,我只想在data_array中提取对应于N = 3的变量。任何使用python netCDF4或nco工具的解决方案都适合我(但不使用cdo),也需要在windows上工作。

1 个答案:

答案 0 :(得分:2)

如果我理解正确,你会使用古怪的术语并打电话给"组"和"变量"别人称之为"变量"分别是" hyperslab" NCO解决方案是对大多数运算符使用hyperslab选项(-d),例如ncks

ncks -d N,2 -v data_array in.nc out.nc

默认情况下,指数从0开始,因此N = 2会为您提供N中的最后一个切片。