- 积分
- 26283
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-6-1
- 最后登录
- 1970-1-1
|
发表于 2013-1-17 03:49:39
|
显示全部楼层
这个读取的方法是对一个文件中的多年数据吧?感觉批量的读取不太一致。- clc; %清屏
- clear; %清空
- datadir='D:\data\降水数据\CPC Unified Gauge-Based Analysis of Daily Precipitation over CONUS\'; %指定批量数据所在的文件夹
- filelist=dir([datadir,'*.nc']); %指定批量数据的类型
- a=filelist(1).name; %查看你要读取的文件的编号。filelist(1).name在window下为第一个标号数据
- b=filelist(2).name; %查看你要读取的文件的编号。filelist(2).name在window下为第二个标号数据
- k=length(filelist);
- for s=1:k
- filename=[datadir,filelist(s).name];
- ncid=netcdf.open(filename,'NC_NOWRITE');
- ncdisp('D:\data\降水数据\CPC Unified Gauge-Based Analysis of Daily Precipitation over CONUS\precip.V1.0.1948.nc'); %在命令窗中显示nc文件的变量
- %任意取其中一个来看数据中所包含的变量特征,以为下面读取数据变量做铺垫
- % ncid = netcdf.open('D:\data\降水数据\CPC Unified Gauge-Based Analysis of Daily Precipitation over CONUS\precip.V1.0.1948.nc','NOWRITE'); %打开nc文件
- % ncdisp('D:\data\降水数据\CPC Unified Gauge-Based Analysis of Daily Precipitation over CONUS\precip.V1.0.1948.nc'); %在命令窗中显示nc文件的变量
- PrecipData = ncread(filename,'precip'); %读入变量precip
- TimeData = ncread(filename,'time'); %读入变量time
- LonData = ncread(filename,'lon'); %读入变量lon
- LatData = ncread(filename,'lat'); %读入变量lat
- netcdf.close(ncid); % 关闭文件
- end;
复制代码 |
|