- 积分
- 2293
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2018-6-10
- 最后登录
- 1970-1-1
|
发表于 2019-12-22 00:15:58
|
显示全部楼层
%批量读取nc文件的方法:
clc; %清屏
clear; %清空
datadir='G:\Hotspot\'; %指定批量数据所在的文件夹
filelist=dir([datadir,'*.nc']); %指定批量数据的类型
a=filelist(1).name; %查看你要读取的文件的编号。filelist(1).name在window下为第一个标号数据
b=filelist(2).name; %查看你要读取的文件的编号。filelist(2).name在window下为第二个标号数据
k=length(filelist);
lon=ncread('G:\DHW\ct5km_dhw-max_v3.1_201904.nc','lon');
lat=ncread('G:\DHW\ct5km_dhw-max_v3.1_201904.nc','lat');
i_inx=find(lon>=104&lon<=122);
j_inx=find(lat>=2&lat<=24);
nx=length(i_inx);
ny=length(j_inx);
HS=zeros(360,440,418);
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文件的变量
HS(:,:,s) = ncread(filename,'hotspot',[i_inx(1) j_inx(1) 1],[nx ny 1]);; %读入变量lat
%读入变量precip
netcdf.close(ncid); % 关闭文件
end
这是nc批量读取数据的方法,网站上一个大大发出来的,具体在哪里我忘记了,隔空感谢大大指点
|
|