read_GLDAS.m.zip_5UV_GLDAS土壤单位_matlab gldas_土壤_读取gldas数据

%% 查看.nc数据文件，定义相关路径 clc,clear % Address='H:\GLDAS_Data\GLDAS1_NC\'; AddressOut='H:\GLDAS_Data\Nc_Out\'; AddressOut2='H:\GLDAS_Data\Nc_Out2\'; gldas=ncinfo('H:\GLDAS_Data\GLDAS1_NC\GLDAS_NOAH10_M.A200505.001.nc'); %地址+文件分开读取需要加[] %%读Noah时一定要用它的数据文件 ncdisp('H:\GLDAS_Data\GLDAS1_NC\GLDAS_NOAH10_M.A200505.001.nc') %查看变量的详细信息 GFA=dir(fullfile(Address,'*.nc')); k=length(GFA); %文件个数 all_data=[]; %% 第一列 lon1=(-179.5:179.5)'; lat1=(-59.5:89.5)'; l=ones(360*150,1); nn=0; ii=0; l2=zeros(360*150,1); ll2=0; for i=1:length(lat1) for j=1:length(lon1) ll2=ll2+1; l2(ll2)=lon1(j); end end %% 重复数组for第二列 lat2=lat1'; B = repmat(lat2,360,1); BB=reshape(B,150*360,1); %% 读取NC文件 for i=1:k i; filepath=Address; file=[filepath,GFA(i).name]; %GFA(I).name 即为第i个.nc文件 % 新建输出文件 AddressOut3=strcat(AddressOut,GFA(i).name(1:27),'txt'); fid=fopen(AddressOut3,'wt'); AddressOut4=strcat(AddressOut2,GFA(i).name(1:27),'txt'); fid2=fopen(AddressOut4,'wt'); % 新建完成 %%%%%% 读取数据 xx=ncinfo(file); %ncdisp(file); lon=ncread(file,xx.Variables(1).Name);%行是经度 lat=ncread(file,xx.Variables(2).Name);%列是纬度 Soil_Moist_data=ncread(file,xx.Variables(18).Name); % 全球土壤水数据(360*150*4) data_Soil=Soil_Moist_data(:,:,1)+Soil_Moist_data(:,:,2)+Soil_Moist_data(:,:,3)+Soil_Moist_data(:,:,4); data_Soil1=reshape(data_Soil,150*360,1); % data3=reshape([-179.5:1:179.5 -59.5:1:89.5],150*360,2); % all_data=[all_data data2]; for j=1:360*150 if isnan(data_Soil1(j)) data_Soil1(j)=0; else data_Soil1(j)=data_Soil1(j); end end data_Soil=reshape(data_Soil1,360,150); Snow_Water=ncread(file,xx.Variables(21).Name); % 全球雪水量数据(360*150*1) Snow_Water1=reshape(Snow_Water,360*150,1); for j=1:360*150 if isnan(Snow_Water1(j)) Snow_Water1(j)=0; else Snow_Water1(j)=Snow_Water1(j); end end Snow_Water=reshape(Snow_Water1,360,150); Canopy_Water=ncread(file,xx.Variables(20).Name); % 全球冠层水数据(360*150*1) Canopy_Water1=reshape(Canopy_Water,360*150,1); for j=1:360*150 if isnan(Canopy_Water1(j)) Canopy_Water1(j)=0; else Canopy_Water1(j)=Canopy_Water1(j); end end Canopy_Water=reshape(Canopy_Water1,360,150); Total_Water1=data_Soil1+Snow_Water1+Canopy_Water1; %全球陆地水总量 Total_Water=data_Soil+Snow_Water+Canopy_Water; %全球陆地水总量 all_data=[all_data Total_Water1]; for j=1:360*150 fprintf(fid,'%6.1f %6.1f %11.5f\n',l2(j),BB(j),Total_Water1(j)); end for p=1:360 for q=1:150 fprintf(fid2,' %6.1f',Total_Water(p,q)); end fprintf(fid2,'\n'); end fclose(fid); fclose(fid2); end % for i=length(all_data(:,1)) % for j=1:length(all_data(1,:)) % if isnan(all_data(i,j)) % all_data(i,j)=0; % else % all_data(i,j)=all_data(i,j); % end % end % end % save(all_data.txt,all_data);