matlab读取nc数据的一个问题

夏夜

一个sst的nc数据ncdisp之后看到如下信息：
analysed_sst  
           Size:       4096x2048x1
           Dimensions: lon,lat,time
           Datatype:   int16
           Attributes:
                       long_name     = 'analysed sea surface temperature'
                       standard_name = 'sea_surface_temperature'
                       type          = 'foundation'
                       units         = 'kelvin'
                       _FillValue    = -3.28e+04
                       add_offset    = 273
                       scale_factor  = 0.01
                       valid_min     = -300
                       valid_max     = 4.5e+03
                       source        = 'REMSS'


然后分别用ncload和ncread读取：

ncload filename.nc                       读取的analysed_sst（赋给变量sst_load） 维数是2048*4096
ncread('filename.nc','analysed_sst') 读取的 analysed_sst  (赋给变量sst_read) 维数是4096*2048

这是维数有区别，也没关系，转置一下。但是sst_load和sst_read的数值竟然不相同。

区别一：sst_load中的FillValue 显示为-3.28e+04，而sst_read对应的显示为NaN

区别二：非FillValue的地方数值也不相同，大概满足下面这个关系式：

            sst_load==(sst_read-add_offset)./scale_factor;  
         (这两者之差数量级1e-06，由于matlab计算时取的精度不同，可能会有这个误差)。



然后换了一个ssh的nc数据：

Grid_0001   
           Size:       721x1440
           Dimensions: NbLatitudes,NbLongitudes
           Datatype:   single
           Attributes:
                       _FillValue   = 1.84e+19
                       long_name    = 'H'
                       units        = 'cm'
                       Date_CNES_JD = 2.35e+04
                       date         = '2014-03-27 00:00:00.000000 UTC'


这个数据没有add_offset和scale_factor,用ncload和ncread读取的数据之后只有区别一，没有区别二。

请问大家遇到这样的问题吗？

是因为数据还是因为这两个命令的本身就存在的差异？

夏夜

help ncread：

The following attribute conventions are applied
     in sequence to VARDATA if the corresponding attribute exists for this
     variable:

        1. Values in VARDATA equal to the '_FillValue' attribute value are
           replaced with NaNs. If '_FillValue' attribute does not exist,
           ncread will query the library for the variable's fill value.
        2. VARDATA is multiplied by the value of 'scale_factor' attribute.
        3. The value of the 'add_offset' attribute is added to VARDATA.

kongfeng0824

ncread一般会直接读取成正确的数值。只有getVar之类的才需要转换。版本必须在2012a及以上

夏夜

kongfeng0824 发表于 2014-9-9 19:18
ncread一般会直接读取成正确的数值。只有getVar之类的才需要转换。版本必须在2012a及以上

ncread读出的数值，只需要转换一下单位就行了，比如sst的开尔文转换摄氏度（如果需要转换的话）。

ncload读出的还要用add_offset、scale_factor转换之后再转换单位。

是这样吧

夏夜

kongfeng0824 发表于 2014-9-9 19:18
ncread一般会直接读取成正确的数值。只有getVar之类的才需要转换。版本必须在2012a及以上

在版本R2013a看到的

Song0o

谢谢楼主的分享。。。。

午后阳光

学习了、、、

ziqiangbuxi

好东西，最近就在处理netcdf数据，解决了我的疑问。谢谢大家！

江南鹤

纠结了一下午这个问题，学习了

sunshine111222

{:eb513:}{:eb513:}{:eb513:}