MATLAB读取NC文件

一、 目的：

          了解NETCDF文件，学会利用MATLAB读取NETCDF文件

二、  撰写时间

          开始时间：2016年12月03日

          完成时间：2016年12月09日

三、知识储备

（一）NETCDF文件

   NetCDF(network Common Data Form)网络通用数据格式是由美国大学大气研究协会（University Corporation for Atmospheric Research，UCAR)的Unidata项目科学家针对科学数据的特点开发的，是一种面向数组型并适于网络共享的数据的描述和编码标准。目前，NetCDF广泛应用于大气科学、水文、海洋学、环境模拟、地球物理等诸多领域。用户可以借助多种方式方便地管理和操作 NetCDF 数据集（百度文库）。

（二）MATLAB中读取NETCDF格式数据的常用函数

1. ncread

1.1 功能

从NetCDF格式的数据源的变量中读取数据（Read data from variable in NetCDF data source）

1.2 语法结构

vardata = ncread(source,varname)

vardata = ncread(source,varname,start,count,stride)

1.3 描述

1.3.1 vardata = ncread(source,varname)

从数据源中读取变量名为varname的变量

1.3.2 vardata = ncread((source,varname,start,count,stride)

（1）start

varname所指定变量的每一维的开始读取的位置

（2）count

从start指定的开始位置算起，一共读取的每一维要素的数目

（3）stride

从start开始，每一维读取的数目为count时，每一维的读取的步长

1.4 程序示例

[plain] view plain copy

1. %*************************************************************************%  

2. % 程序目的：测试ncread的两种读取nc文件的方法  

3. % 时间：2016年12月09日  

4. % 作者：东方易  

5. %输入数据：F:2科研3书籍撰写2_CSDN1_NCREAD1_TestDataspei01.nc  

6. %*************************************************************************%  

7. clc;  

8. clear all;  

9. InPath = 'F:2科研3书籍撰写2_CSDN1_NCREAD1_TestData';  

10. ncdisp(strcat(InPath,'spei01.nc')) %获取所读取nc文件的基本信息  

11.  

12. %----------------------vardata = ncread(source,varname)开始----------------%  

13. % vardata = ncread(source,varname)  

14. source1 = strcat(InPath,'spei01.nc');  

15. varname1 = 'lon';  

16. Lon_Data = ncread(source1,varname1);  

17. %----------------------vardata = ncread(source,varname)结束----------------%  

18.  

19. %--------vardata = ncread((source,varname,start,count,stride)开始---------%  

20. source2 = source1;  

21. varname2 = 'spei'; %spei是一个三维变量720*360*1368（lon*lat*time）  

22. stlo = 100; stla = 100; stti = 100;  

23. start = [stlo,stla,stti]; %创建start变量，每一维开始的变量  

24. locount = 100; lacount = 100; ticount = 3;  

25. count = [locount,lacount,ticount]; %创建count变量，从每一维的start开始读取的总数目  

26. stride1 = [1,1,1]; %设置读取的步长  

27. stride2 = [2,2,2]; %设置读取的步长  

28. spei_stride1 = ncread(source2,varname2,start,count,stride1);  

29. spei_stride2 = ncread(source2,varname2,start,count,stride2);  

30. %--------vardata = ncread((source,varname,start,count,stride)结束---------%  

31.  

32. disp('请注意比较spei_stride1和spei_stride2的不同！')  

33. disp('finished!')