- 积分
- 3752
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2019-3-1
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
时常有同学在群里问二爷指导手册里面的文件循环读取怎么做,遇到困哪,今天分享一下循环读取文件以nc文件为例;下面是我实操一个例子。
%% 批量读取HadCM3 2-16ka数据
list=dir(['G:\数据库\HadCM3\未加冰盖模拟结果\','*.nc']);%1-8为2,4...16ka数据
len=length(list);%求nc文件的数量
%% 读取HadCM3降水数据并处理
for i=1:len
ncid=[list(i).folder,'\',list(i).name];%拼装每个nc文件的绝对路劲
pre=ncread(ncid,'precip');%由于precip是经度*纬度*高度*时间,时间为30年*12个月,高度项为1
pre=squeeze(pre);%压缩高度维度
pre=reshape(pre,[96,73,30,12]);%将precip处理成经度*纬度*30年*12个月
pre=mean(pre(:,:,:,[6 8]),4);%计算每年6-8月平均降水
pre=mean(pre,3);%计算30年6-8月平均降水
pre=pre-pre0;%计算pre1-8(30年6-8月平均降水)与PI(30年6-8月平均降水)差值
eval(['pre',num2str(i),'=','pre(:,:)',';']);%批量命名变量eval
end
% pre1-8依次为2,4...16ka30年6-8月平均降水数据
|
|