循环读取文件

微醉的阳光

时常有同学在群里问二爷指导手册里面的文件循环读取怎么做，遇到困哪，今天分享一下循环读取文件以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月平均降水数据

dm12

谢谢分享，学习了