MATLAB中利用ncread函数读取nc文件 |
您所在的位置:网站首页 › MATLAB读取文件函数差异 › 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!')
|
今日新闻 |
点击排行 |
|
推荐新闻 |
图片新闻 |
|
专题文章 |
CopyRight 2018-2019 实验室设备网 版权所有 win10的实时保护怎么永久关闭 |