MATLAB中利用ncread函數讀取nc文件

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 程序示例

%*************************************************************************%
% 程序目的：測試ncread的兩種讀取nc文件的方法
% 時間：2016年12月09日
% 作者：東方易
%輸入數據：F:\2科研\3書籍撰寫\2_CSDN\1_NCREAD\1_TestData\spei01.nc
%*************************************************************************%
clc;
clear all;
InPath = 'F:\2科研\3書籍撰寫\2_CSDN\1_NCREAD\1_TestData\';
ncdisp(strcat(InPath,'spei01.nc')) %獲取所讀取nc文件的基本信息

%----------------------vardata = ncread(source,varname)開始----------------%
% vardata = ncread(source,varname)
source1 = strcat(InPath,'spei01.nc');
varname1 = 'lon';
Lon_Data = ncread(source1,varname1);
%----------------------vardata = ncread(source,varname)結束----------------%

%--------vardata = ncread((source,varname,start,count,stride)開始---------%
source2 = source1;
varname2 = 'spei'; %spei是一個三維變量720*360*1368（lon*lat*time）
stlo = 100; stla = 100; stti = 100; 
start = [stlo,stla,stti]; %創建start變量，每一維開始的變量
locount = 100; lacount = 100; ticount = 3;
count = [locount,lacount,ticount]; %創建count變量，從每一維的start開始讀取的總數目
stride1 = [1,1,1]; %設置讀取的步長
stride2 = [2,2,2]; %設置讀取的步長
spei_stride1 = ncread(source2,varname2,start,count,stride1);
spei_stride2 = ncread(source2,varname2,start,count,stride2);
%--------vardata = ncread((source,varname,start,count,stride)結束---------%

disp('請注意比較spei_stride1和spei_stride2的不同！')
disp('finished!')