全球范围内12个典型区域提取(total runoff)

最新推荐文章于 2021-12-16 12:30:11 发布
最新推荐文章于 2021-12-16 12:30:11 发布
阅读量523 收藏 7
点赞数 1
分类专栏： MATLAB+NC
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/jerry003/article/details/109075975
版权

代码内容包括：

  • 空间插值griddata函数的使用
  • 用mask矩阵提取全球范围的12个典型区域提取
  • 替换xticklabels
% 【摘要】读取.grd文件；对不同分辨率数据网格插值；对不重合区域取公共交集区,
%         比较全球12个典型区域的total runoff的值（CLM45 .VS. CLM50）与UNH-GRDC；
cd 'D:\数据\UNH-GRDC\World Runoff Data\runoff_grd\'
%（1）读 .grd 文件的 unh_grdc 月平均值
clc; clear; close all
flistobs = dir('D:\数据\UNH-GRDC\World Runoff Data\runoff_grd\*.grd');
for i =2:13                  % 留意读取的文件序号
    filename = flistobs(i).name;
    delimiterIn = ' ';       % 设置定界符，从grd文件复制
    headerlinesIn = 6;       % 设置需要跳过的行数，跳过6行，从第7行开始
    datastrct = importdata(filename,delimiterIn,headerlinesIn);
    rof = datastrct.data;
    rof = flip(rof);         % 数据可能是倒置的，可以通过pcolor查看
    rof(rof == -9999) = NaN; % 替换缺省值为 Nan 
    rofcmp_mon{i-1} = rof;    % 以cell格式存储
end

%（2）读 clm 模拟的多年月平均值（CLM45 or CLM50）
load('D:\数据\UNH-GRDC\World Runoff Data\runoff_grd\simres_clm45.mat')

%（3）对 unh-grdc 数据进行空间插值计算（0.5*0.5转到1.25*0.9375）
fpath = 'D:\Models_Output\Clm45Sp\lnd\hist\';
flistsim = dir([fpath,'*.nc']);
fname = flistsim(1).name;
xx = ncread([fpath,fname],'lon');
yy = ncread([fpath,fname],'lat');
for m = 1:length(xx)
    if  xx(m) >= 180
        xx(m) = xx(m) - 360;
    else
        continue
    end
end
xx = [xx(145:end);xx(1:144)];  % 将经度后半段(>=180)拼接到前面
[lonsim,latsim] = meshgrid(xx,yy);

x = -180:0.5:179.5;
y = -55.5:0.5:82.5;
[lonobs,latobs] = meshgrid(x,y);
for k = 1:12
    ttmp = rofcmp_mon{k};
    cmprofup_mon{k} = griddata(lonobs(:),latobs(:),ttmp(:),lonsim,latsim,'nearest');
end

%（4）确定全局情况下两者共同交集网格（利用Nan 与 非Nan 的计算结果为 Nan 确定交集区域）
commask = cmprofup_mon{1} - rofsim_mon{1};   % comaree中非nan的网格就是两者共有网格
commask(~isnan(commask)) = 1;
for m = 1:12
    tmpcmp = cmprofup_mon{m};
    cmprofup_mon{m} = tmpcmp.*commask;
    tmpsim = rofsim_mon{m};
    rofsim_mon{m} = tmpsim.*commask;
end



%（5）划定各个典型区域坐标范围
% West_Siberia(50-70N，60-90E)  East_Siberia(50-70N,90-140E)  Canada(50-70N,90-170W)
% Amazon_Basin(10S-0，70-50W)   Congo(10S-5N,10-30E)          India(10-30N,70-90E)
% South_Africa(35-10S,10-40E)   Eastern_China(25-40N,110-125E)
% Central_Europ(40-50N,0-40E)   Sahara+Arabla(10-30N,0-50E)   Western_US(30-50N,130-110W)
% Australia(40-10S,110-160E)
reginal_name = ...
{['West Siberia(50-70N，60-90E)'] ['East Siberia(50-70N,90-140E)']...
 ['Canada(50-70N,90-170W)']       ['Amazon Basin(10S-0，70-50W)']...
 ['Congo(10S-5N,10-30E)']         ['India(10-30N,70-90E)']...
 ['South Africa(35-10S,10-40E)']  ['Eastern China(25-40N,110-125E)']...
 ['Central Europ(40-50N,0-40E)']  ['Sahara-Arabla(10-30N,0-50E)']...
 ['Western US(30-50N,130-110W)']  ['Australia(40-10S,110-160E)']};
dlon = 1.25;
dlat = 0.9375;

%（5-1）划定各个区域范围
West_Siberia_lat = [floor((50+90)/dlat):ceil((70+90)/dlat)];
West_Siberia_lon = [floor((60+180)/dlon):ceil((90+180)/dlon)];

East_Siberia_lat = [floor((50+90)/dlat):ceil((70+90)/dlat)];
East_Siberia_lon = [floor((90+180)/dlon):ceil((140+180)/dlon)];

Canada_lat = [floor((50+90)/dlat):ceil((70+90)/dlat)];
Canada_lon = [floor((180-170)/dlon):ceil((180-90)/dlon)];

Amazon_Basin_lat = [floor((90-10)/dlat):ceil((90-0)/dlat)];
Amazon_Basin_lon = [floor((180-70)/dlon):ceil((180-50)/dlon)];

Congo_lat = [floor((90-10)/dlat):ceil((90+5)/dlat)];
Congo_lon = [floor((10+180)/dlon):ceil((30+180)/dlon)];

India_lat = [floor((90+10)/dlat):ceil((90+30)/dlat)];
India_lon = [floor((70+180)/dlon):ceil((90+180)/dlon)];

South_Africa_lat = [floor((90-35)/dlat):ceil((90-10)/dlat)];
South_Africa_lon = [floor((10+180)/dlon):ceil((40+180)/dlon)];

Eastern_China_lat = [floor((90+25)/dlat):ceil((90+40)/dlat)];
Eastern_China_lon = [floor((110+180)/dlon):ceil((125+180)/dlon)];

Central_Europ_lat = [floor((40+90)/dlat):ceil((50+90)/dlat)];
Central_Europ_lon = [floor((180+0)/dlon):ceil((180+40)/dlon)];

Sahara_Arabla_lat = [floor((90+10)/dlat):ceil((90+30)/dlat)];
Sahara_Arabla_lon = [floor((180+0)/dlon):ceil((180+50)/dlon)];

Western_US_lat = [floor((90+30)/dlat):ceil((90+50)/dlat)];
Western_US_lon = [floor((180-130)/dlon):ceil((180-110)/dlon)];

Australia_lat = [floor((90-40)/dlat):ceil((90-10)/dlat)];
Australia_lon = [floor((110+180)/dlon):ceil((180+160)/dlon)];

%（5-2）提取各个区域网格
for i = 1:12  
    tmpobs = cmprofup_mon{i}; % 逐月提取观测值
    West_Siberia_obs{i} = tmpobs(West_Siberia_lat, West_Siberia_lon);
    East_Siberia_obs{i} = tmpobs(East_Siberia_lat, East_Siberia_lon);
    Canada_obs{i}       = tmpobs(Canada_lat,Canada_lon);
    Amazon_Basin_obs{i} = tmpobs(Amazon_Basin_lat,Amazon_Basin_lon);
    Congo_obs{i}        = tmpobs(Congo_lat,Congo_lon);
    India_obs{i}        = tmpobs(India_lat,India_lon);
    South_Africa_obs{i} = tmpobs(South_Africa_lat, South_Africa_lon);
    Eastern_China_obs{i}= tmpobs(Eastern_China_lat,Eastern_China_lon);
    Central_Europ_obs{i}= tmpobs(Central_Europ_lat,Central_Europ_lon);
    Sahara_Arabla_obs{i}= tmpobs(Sahara_Arabla_lat,Sahara_Arabla_lon);
    Western_US_obs{i}   = tmpobs(Western_US_lat,Western_US_lon);
    Australia_obs{i}    = tmpobs(Australia_lat,Australia_lon);
    
    tmpsim = rofsim_mon{i};   %  逐月提取模拟值
    West_Siberia_sim{i} = tmpsim(West_Siberia_lat, West_Siberia_lon);
    East_Siberia_sim{i} = tmpsim(East_Siberia_lat, East_Siberia_lon);
    Canada_sim{i}       = tmpsim(Canada_lat,Canada_lon);
    Amazon_Basin_sim{i} = tmpsim(Amazon_Basin_lat,Amazon_Basin_lon);
    Congo_sim{i}        = tmpsim(Congo_lat,Congo_lon);
    India_sim{i}        = tmpsim(India_lat,India_lon);
    South_Africa_sim{i} = tmpsim(South_Africa_lat, South_Africa_lon);
    Eastern_China_sim{i}= tmpsim(Eastern_China_lat,Eastern_China_lon);
    Central_Europ_sim{i}= tmpsim(Central_Europ_lat,Central_Europ_lon);
    Sahara_Arabla_sim{i}= tmpsim(Sahara_Arabla_lat,Sahara_Arabla_lon);
    Western_US_sim{i}   = tmpsim(Western_US_lat,Western_US_lon);
    Australia_sim{i}    = tmpsim(Australia_lat,Australia_lon);
end

%（6）比较逐月runoff数据
for i = 1:12
    tmpreg = West_Siberia_obs{i};
    West_Siberia_obsmon(i) = mean(tmpreg(~isnan(tmpreg)));
    tmpreg = West_Siberia_sim{i};
    West_Siberia_simmon(i) = mean(tmpreg(~isnan(tmpreg)));
    
    tmpreg = East_Siberia_obs{i};
    East_Siberia_obsmon(i) = mean(tmpreg(~isnan(tmpreg)));
    tmpreg = East_Siberia_sim{i};
    East_Siberia_simmon(i) = mean(tmpreg(~isnan(tmpreg)));
    
    tmpreg = Canada_obs{i};
    Canada_obsmon(i) = mean(tmpreg(~isnan(tmpreg)));
    tmpreg = Canada_sim{i};
    Canada_simmon(i) = mean(tmpreg(~isnan(tmpreg)));
    
    tmpreg = Amazon_Basin_obs{i};
    Amazon_Basin_obsmon(i) = mean(tmpreg(~isnan(tmpreg)));
    tmpreg = Amazon_Basin_sim{i};
    Amazon_Basin_simmon(i) = mean(tmpreg(~isnan(tmpreg)));
    
    tmpreg = Congo_obs{i};
    Congo_obsmon(i) = mean(tmpreg(~isnan(tmpreg)));
    tmpreg = Congo_sim{i};
    Congo_simmon(i) = mean(tmpreg(~isnan(tmpreg)));
    
    tmpreg = India_obs{i};
    India_obsmon(i) = mean(tmpreg(~isnan(tmpreg)));
    tmpreg = India_sim{i};
    India_simmon(i) = mean(tmpreg(~isnan(tmpreg)));
    
    tmpreg = South_Africa_obs{i};
    South_Africa_obsmon(i) = mean(tmpreg(~isnan(tmpreg)));
    tmpreg = South_Africa_sim{i};
    South_Africa_simmon(i) = mean(tmpreg(~isnan(tmpreg)));
    
    tmpreg = Eastern_China_obs{i};
    Eastern_China_obsmon(i) = mean(tmpreg(~isnan(tmpreg)));
    tmpreg = Eastern_China_sim{i};
    Eastern_China_simmon(i) = mean(tmpreg(~isnan(tmpreg)));
    
    tmpreg = Central_Europ_obs{i};
    Central_Europ_obsmon(i) = mean(tmpreg(~isnan(tmpreg)));
    tmpreg = Central_Europ_sim{i};
    Central_Europ_simmon(i) = mean(tmpreg(~isnan(tmpreg)));
    
    tmpreg = Sahara_Arabla_obs{i};
    Sahara_Arabla_obsmon(i) = mean(tmpreg(~isnan(tmpreg)));
    tmpreg = Sahara_Arabla_sim{i};
    Sahara_Arabla_simmon(i) = mean(tmpreg(~isnan(tmpreg)));
    
    tmpreg = Western_US_obs{i};
    Western_US_obsmon(i) = mean(tmpreg(~isnan(tmpreg)));
    tmpreg = Western_US_sim{i};
    Western_US_simmon(i) = mean(tmpreg(~isnan(tmpreg)));
    
    tmpreg = Australia_obs{i};
    Australia_obsmon(i) = mean(tmpreg(~isnan(tmpreg)));
    tmpreg = Australia_sim{i};
    Australia_simmon(i) = mean(tmpreg(~isnan(tmpreg)));
end

%(7) 画每个区域逐月runoff图，为方便画图，将此12个区域汇总起来
regset_45sim{1}  = West_Siberia_simmon;
regset_45sim{2}  = East_Siberia_simmon;
regset_45sim{3}  = Canada_simmon;
regset_45sim{4}  = Amazon_Basin_simmon;
regset_45sim{5}  = Congo_simmon;
regset_45sim{6}  = India_simmon;
regset_45sim{7}  = South_Africa_simmon;
regset_45sim{8}  = Eastern_China_simmon;
regset_45sim{9}  = Central_Europ_simmon;
regset_45sim{10} = Sahara_Arabla_simmon;
regset_45sim{11} = Western_US_simmon;
regset_45sim{12} = Australia_simmon;

regset_obs{1}  = West_Siberia_obsmon;
regset_obs{2}  = East_Siberia_obsmon;
regset_obs{3}  = Canada_obsmon;
regset_obs{4}  = Amazon_Basin_obsmon;
regset_obs{5}  = Congo_obsmon;
regset_obs{6}  = India_obsmon;
regset_obs{7}  = South_Africa_obsmon;
regset_obs{8}  = Eastern_China_obsmon;
regset_obs{9}  = Central_Europ_obsmon;
regset_obs{10} = Sahara_Arabla_obsmon;
regset_obs{11} = Western_US_obsmon;
regset_obs{12} = Australia_obsmon;

%% 重复上述计算，导入 clm50 的数据
%（2）读 clm 模拟的多年月平均值（CLM45 or CLM50）
load('D:\数据\UNH-GRDC\World Runoff Data\runoff_grd\simres_clm50.mat')

%（3）对 unh-grdc 数据进行空间插值计算（0.5*0.5转到1.25*0.9375）
fpath = 'D:\Models_Output\Clm50Sp\lnd\hist\';
flistsim = dir([fpath,'*.nc']);
fname = flistsim(1).name;
xx = ncread([fpath,fname],'lon');
yy = ncread([fpath,fname],'lat');
for m = 1:length(xx)
    if  xx(m) >= 180
        xx(m) = xx(m) - 360;
    else
        continue
    end
end
xx = [xx(145:end);xx(1:144)];  % 将经度后半段(>=180)拼接到前面
[lonsim,latsim] = meshgrid(xx,yy);

x = -180:0.5:179.5;
y = -55.5:0.5:82.5;
[lonobs,latobs] = meshgrid(x,y);
for k = 1:12
    ttmp = rofcmp_mon{k};
    cmprofup_mon{k} = griddata(lonobs(:),latobs(:),ttmp(:),lonsim,latsim,'nearest');
end

%（4）确定全局情况下两者共同交集网格（利用Nan 与 非Nan 的计算结果为 Nan 确定交集区域）
commask = cmprofup_mon{1} - rofsim_mon{1};   % comaree中非nan的网格就是两者共有网格
commask(~isnan(commask)) = 1;
for m = 1:12
    tmpcmp = cmprofup_mon{m};
    cmprofup_mon{m} = tmpcmp.*commask;
    tmpsim = rofsim_mon{m};
    rofsim_mon{m} = tmpsim.*commask;
end

%（5）划定各个典型区域坐标范围
% West_Siberia(50-70N，60-90E)  East_Siberia(50-70N,90-140E)  Canada(50-70N,90-170W)
% Amazon_Basin(10S-0，70-50W)   Congo(10S-5N,10-30E)          India(10-30N,70-90E)
% South_Africa(35-10S,10-40E)   Eastern_China(25-40N,110-125E)
% Central_Europ(40-50N,0-40E)   Sahara+Arabla(10-30N,0-50E)   Western_US(30-50N,130-110W)
% Australia(40-10S,110-160E)
reginal_name = ...
{['West Siberia(50-70N，60-90E)'] ['East Siberia(50-70N,90-140E)']...
 ['Canada(50-70N,90-170W)']       ['Amazon Basin(10S-0，70-50W)']...
 ['Congo(10S-5N,10-30E)']         ['India(10-30N,70-90E)']...
 ['South Africa(35-10S,10-40E)']  ['Eastern China(25-40N,110-125E)']...
 ['Central Europ(40-50N,0-40E)']  ['Sahara-Arabla(10-30N,0-50E)']...
 ['Western US(30-50N,130-110W)']  ['Australia(40-10S,110-160E)']};
dlon = 1.25;
dlat = 0.9375;

%（5-1）划定各个区域范围
West_Siberia_lat = [floor((50+90)/dlat):ceil((70+90)/dlat)];
West_Siberia_lon = [floor((60+180)/dlon):ceil((90+180)/dlon)];

East_Siberia_lat = [floor((50+90)/dlat):ceil((70+90)/dlat)];
East_Siberia_lon = [floor((90+180)/dlon):ceil((140+180)/dlon)];

Canada_lat = [floor((50+90)/dlat):ceil((70+90)/dlat)];
Canada_lon = [floor((180-170)/dlon):ceil((180-90)/dlon)];

Amazon_Basin_lat = [floor((90-10)/dlat):ceil((90-0)/dlat)];
Amazon_Basin_lon = [floor((180-70)/dlon):ceil((180-50)/dlon)];

Congo_lat = [floor((90-10)/dlat):ceil((90+5)/dlat)];
Congo_lon = [floor((10+180)/dlon):ceil((30+180)/dlon)];

India_lat = [floor((90+10)/dlat):ceil((90+30)/dlat)];
India_lon = [floor((70+180)/dlon):ceil((90+180)/dlon)];

South_Africa_lat = [floor((90-35)/dlat):ceil((90-10)/dlat)];
South_Africa_lon = [floor((10+180)/dlon):ceil((40+180)/dlon)];

Eastern_China_lat = [floor((90+25)/dlat):ceil((90+40)/dlat)];
Eastern_China_lon = [floor((110+180)/dlon):ceil((125+180)/dlon)];

Central_Europ_lat = [floor((40+90)/dlat):ceil((50+90)/dlat)];
Central_Europ_lon = [floor((180+0)/dlon):ceil((180+40)/dlon)];

Sahara_Arabla_lat = [floor((90+10)/dlat):ceil((90+30)/dlat)];
Sahara_Arabla_lon = [floor((180+0)/dlon):ceil((180+50)/dlon)];

Western_US_lat = [floor((90+30)/dlat):ceil((90+50)/dlat)];
Western_US_lon = [floor((180-130)/dlon):ceil((180-110)/dlon)];

Australia_lat = [floor((90-40)/dlat):ceil((90-10)/dlat)];
Australia_lon = [floor((110+180)/dlon):ceil((180+160)/dlon)];

%（5-2）提取各个区域网格
for i = 1:12  
    tmpobs = cmprofup_mon{i}; % 逐月提取观测值
    West_Siberia_obs{i} = tmpobs(West_Siberia_lat, West_Siberia_lon);
    East_Siberia_obs{i} = tmpobs(East_Siberia_lat, East_Siberia_lon);
    Canada_obs{i}       = tmpobs(Canada_lat,Canada_lon);
    Amazon_Basin_obs{i} = tmpobs(Amazon_Basin_lat,Amazon_Basin_lon);
    Congo_obs{i}        = tmpobs(Congo_lat,Congo_lon);
    India_obs{i}        = tmpobs(India_lat,India_lon);
    South_Africa_obs{i} = tmpobs(South_Africa_lat, South_Africa_lon);
    Eastern_China_obs{i}= tmpobs(Eastern_China_lat,Eastern_China_lon);
    Central_Europ_obs{i}= tmpobs(Central_Europ_lat,Central_Europ_lon);
    Sahara_Arabla_obs{i}= tmpobs(Sahara_Arabla_lat,Sahara_Arabla_lon);
    Western_US_obs{i}   = tmpobs(Western_US_lat,Western_US_lon);
    Australia_obs{i}    = tmpobs(Australia_lat,Australia_lon);
    
    tmpsim = rofsim_mon{i};   %  逐月提取模拟值
    West_Siberia_sim{i} = tmpsim(West_Siberia_lat, West_Siberia_lon);
    East_Siberia_sim{i} = tmpsim(East_Siberia_lat, East_Siberia_lon);
    Canada_sim{i}       = tmpsim(Canada_lat,Canada_lon);
    Amazon_Basin_sim{i} = tmpsim(Amazon_Basin_lat,Amazon_Basin_lon);
    Congo_sim{i}        = tmpsim(Congo_lat,Congo_lon);
    India_sim{i}        = tmpsim(India_lat,India_lon);
    South_Africa_sim{i} = tmpsim(South_Africa_lat, South_Africa_lon);
    Eastern_China_sim{i}= tmpsim(Eastern_China_lat,Eastern_China_lon);
    Central_Europ_sim{i}= tmpsim(Central_Europ_lat,Central_Europ_lon);
    Sahara_Arabla_sim{i}= tmpsim(Sahara_Arabla_lat,Sahara_Arabla_lon);
    Western_US_sim{i}   = tmpsim(Western_US_lat,Western_US_lon);
    Australia_sim{i}    = tmpsim(Australia_lat,Australia_lon);
end

%（6）比较逐月runoff数据
for i = 1:12
    tmpreg = West_Siberia_obs{i};
    West_Siberia_obsmon(i) = mean(tmpreg(~isnan(tmpreg)));
    tmpreg = West_Siberia_sim{i};
    West_Siberia_simmon(i) = mean(tmpreg(~isnan(tmpreg)));
    
    tmpreg = East_Siberia_obs{i};
    East_Siberia_obsmon(i) = mean(tmpreg(~isnan(tmpreg)));
    tmpreg = East_Siberia_sim{i};
    East_Siberia_simmon(i) = mean(tmpreg(~isnan(tmpreg)));
    
    tmpreg = Canada_obs{i};
    Canada_obsmon(i) = mean(tmpreg(~isnan(tmpreg)));
    tmpreg = Canada_sim{i};
    Canada_simmon(i) = mean(tmpreg(~isnan(tmpreg)));
    
    tmpreg = Amazon_Basin_obs{i};
    Amazon_Basin_obsmon(i) = mean(tmpreg(~isnan(tmpreg)));
    tmpreg = Amazon_Basin_sim{i};
    Amazon_Basin_simmon(i) = mean(tmpreg(~isnan(tmpreg)));
    
    tmpreg = Congo_obs{i};
    Congo_obsmon(i) = mean(tmpreg(~isnan(tmpreg)));
    tmpreg = Congo_sim{i};
    Congo_simmon(i) = mean(tmpreg(~isnan(tmpreg)));
    
    tmpreg = India_obs{i};
    India_obsmon(i) = mean(tmpreg(~isnan(tmpreg)));
    tmpreg = India_sim{i};
    India_simmon(i) = mean(tmpreg(~isnan(tmpreg)));
    
    tmpreg = South_Africa_obs{i};
    South_Africa_obsmon(i) = mean(tmpreg(~isnan(tmpreg)));
    tmpreg = South_Africa_sim{i};
    South_Africa_simmon(i) = mean(tmpreg(~isnan(tmpreg)));
    
    tmpreg = Eastern_China_obs{i};
    Eastern_China_obsmon(i) = mean(tmpreg(~isnan(tmpreg)));
    tmpreg = Eastern_China_sim{i};
    Eastern_China_simmon(i) = mean(tmpreg(~isnan(tmpreg)));
    
    tmpreg = Central_Europ_obs{i};
    Central_Europ_obsmon(i) = mean(tmpreg(~isnan(tmpreg)));
    tmpreg = Central_Europ_sim{i};
    Central_Europ_simmon(i) = mean(tmpreg(~isnan(tmpreg)));
    
    tmpreg = Sahara_Arabla_obs{i};
    Sahara_Arabla_obsmon(i) = mean(tmpreg(~isnan(tmpreg)));
    tmpreg = Sahara_Arabla_sim{i};
    Sahara_Arabla_simmon(i) = mean(tmpreg(~isnan(tmpreg)));
    
    tmpreg = Western_US_obs{i};
    Western_US_obsmon(i) = mean(tmpreg(~isnan(tmpreg)));
    tmpreg = Western_US_sim{i};
    Western_US_simmon(i) = mean(tmpreg(~isnan(tmpreg)));
    
    tmpreg = Australia_obs{i};
    Australia_obsmon(i) = mean(tmpreg(~isnan(tmpreg)));
    tmpreg = Australia_sim{i};
    Australia_simmon(i) = mean(tmpreg(~isnan(tmpreg)));
end

%(7) 画每个区域逐月runoff图，为方便画图，将此12个区域汇总起来
regset_50sim{1}  = West_Siberia_simmon;
regset_50sim{2}  = East_Siberia_simmon;
regset_50sim{3}  = Canada_simmon;
regset_50sim{4}  = Amazon_Basin_simmon;
regset_50sim{5}  = Congo_simmon;
regset_50sim{6}  = India_simmon;
regset_50sim{7}  = South_Africa_simmon;
regset_50sim{8}  = Eastern_China_simmon;
regset_50sim{9}  = Central_Europ_simmon;
regset_50sim{10} = Sahara_Arabla_simmon;
regset_50sim{11} = Western_US_simmon;
regset_50sim{12} = Australia_simmon;

regset_obs{1}  = West_Siberia_obsmon;
regset_obs{2}  = East_Siberia_obsmon;
regset_obs{3}  = Canada_obsmon;
regset_obs{4}  = Amazon_Basin_obsmon;
regset_obs{5}  = Congo_obsmon;
regset_obs{6}  = India_obsmon;
regset_obs{7}  = South_Africa_obsmon;
regset_obs{8}  = Eastern_China_obsmon;
regset_obs{9}  = Central_Europ_obsmon;
regset_obs{10} = Sahara_Arabla_obsmon;
regset_obs{11} = Western_US_obsmon;
regset_obs{12} = Australia_obsmon;


%% 画图篇
figure(1)
set(gcf,'Color','w');
for i = 1:3
    subplot('Position',[0.1+0.25*(i-1) 0.78 0.22 0.17]);
    plot(regset_obs{i},'k-.*','LineWidth',1)
    hold on
    plot(regset_45sim{i},'b-','LineWidth',2)
    hold on
    plot(regset_50sim{i},'r-','LineWidth',2)
    set(gca,'FontSize',12,'LineWidth',1.5);
    title(reginal_name{i});
    xlim([0,13]);
    xticks(1:13);
    xticklabels({'J','F','M','A','M','J','J','A','S','O','N','D',''});
    if i == 1
       ylabel('runoff (mm/month)' )
       legend('GRDC','CLM45','CLM50')
    end
end

for i = 4:6
    subplot('Position',[0.1+0.25*(i-4) 0.55 0.22 0.17]);
    plot(regset_obs{i},'k-.*','LineWidth',1)
    hold on
    plot(regset_45sim{i},'b-','LineWidth',2)
    hold on
    plot(regset_50sim{i},'r-','LineWidth',2)
    set(gca,'FontSize',12,'LineWidth',1.5);
    title(reginal_name{i});
    xlim([0,13]);
    xticks(1:13);
    xticklabels({'J','F','M','A','M','J','J','A','S','O','N','D',''});
    if i == 4
       ylabel('runoff (mm/month)' )
    end
end

for i = 7:9
    subplot('Position',[0.1+0.25*(i-7) 0.32 0.22 0.17]);
    plot(regset_obs{i},'k-.*','LineWidth',1)
    hold on
    plot(regset_45sim{i},'b-','LineWidth',2)
    hold on
    plot(regset_50sim{i},'r-','LineWidth',2)
    set(gca,'FontSize',12,'LineWidth',1.5);
    title(reginal_name{i});
    xlim([0,13]);
    xticks(1:13);
    xticklabels({'J','F','M','A','M','J','J','A','S','O','N','D',''});
    if i == 7
       ylabel('runoff (mm/month)' )
    end
end

for i = 10:12
    subplot('Position',[0.1+0.25*(i-10) 0.09 0.22 0.17]);
    plot(regset_obs{i},'k-.*','LineWidth',1)
    hold on
    plot(regset_45sim{i},'b-','LineWidth',2)
    hold on
    plot(regset_50sim{i},'r-','LineWidth',2)
    set(gca,'FontSize',12,'LineWidth',1.5);
    title(reginal_name{i});
    xlim([0,13]);
    xticks(1:13);
    xticklabels({'J','F','M','A','M','J','J','A','S','O','N','D',''});
    if i == 10
       ylabel('runoff (mm/month)' )
    end
end

比较

大雨海深
关注
  • 1
    点赞
  • 踩
  • 7
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 2
    评论
专栏目录
GRDC径流数据处理
hengcall的博客
09-27 4792
GRDC是国际径流共享网站。下载下来是txt格式如下图： 提供一个python小脚本，输入文件，处理成excel年-月-径流的格式。 import pandas as pd import numpy as np pth = 'F:/STA-RUNOFF/AMU/NUSUR.txt' # orignal file save_pth = 'F:/STA-RUNOFF/AMU/' ...
language-roff：Atom对Unix手册页（GroffTroff）和历史RUNOFF语言的支持
02-01
罗夫语言支持 支持Atom中的Unix手册页（Troff / Groff），以及整个Roff排版语言。...覆盖范围 每个宏包： ， ， ， ， ， ， 每个预处理： ， ， ， ， ， ， ， GNU Troff扩展 古老的语言，支持
Matlab 读取.grd文件
jerry003的博客
10-10 9752
clc; clear; close all flistobs = dir('D:\数据\UNH-GRDC\World Runoff Data\runoff_grd\*.grd'); % compsite monthly for i =2:13 % 文件序号 filename = flistobs(i).name; delimiterIn = ' '; headerlinesIn = 6; datastrct = importdata(filename,delimite
全球12个典型区域的surface runoff & subsrf runoff比较
jerry003的博客
10-14 479
% 【摘要】比较全球12个典型区域的surface runoff & subsrf runoff的值（CLM45 .VS. CLM50） clc; clear; close all %（1）读 clm 模拟的多年月平均值（CLM45 or CLM50） load('D:\数据\UNH-GRDC\World Runoff Data\runoff_grd\clm45_srf_subsrf.mat') srfrof_clm45_mon = srfrofsim_mon; subsrfrof_clm45_mo
5.3.5 区域提取
kateyabc的博客
06-10 622
注：此文知识学习笔记，仅记录完整程序和实现结果，具体原理参见：https://blog.csdn.net/www_doling_net/article/details/8541534https://blog.csdn.net/shenziheng1/article/category/6114053/4vtkExtractVOI：感兴趣区域（Volum of Interest，VOI）是指图像内部的...
水文概念之：runoff 和streamflow
水滴石穿，绳锯木断
06-18 7688
runoff是坡面产流、河道汇流两个部分，而streamflow是指的河道汇流部分。 通俗理解：streamflow一般是汇流之后的，径流，单位一般是m3/s。runoff一般理解为净流深，单位是mm，可以没有汇流。GLDAS给的径流结果一般用runoff，不用streamflow，因为没有汇流模块。 streamflow是出口断面的径流量，runoff是流域整个汇流面积平均的径流深。一般...
MISDc_code.rar_rainfall runoff_径流 预测_径流预测_洪水_降雨
07-15
MISDc是可以用于洪水预测的降雨径流模型，采用Matlab编写，带有注释。
Surface runoff in urban catchments morphological identification
05-13
an example of hydrology modeling
A combined deep CNN-RNN network for rainfall-runoff modelling
最新发布
07-11
本研究选取了六个不同的人工智能模型对Bardha流域进行月度降雨径流模拟，包括： - 人工神经网络（Artificial Neural Network, ANN） - K近邻回归模型（K-Nearest Neighbors Regression, KNN） - 极限梯度提升回归...
RunOff:进行径流选举的程序
03-10
在等级选择系统中，选民可以投票多个​​候选人。 他们不仅可以为他们的首选投票，还可以按照优先顺序对候选人进行排名。 优先选择投票系统是即时径流系统。 在即时径流选举中，选民可以根据自己的意愿对候选人进行...
LTE学习---PLMN选网（23.122协议R9）
heyJoie的博客
09-02 9632
PLMN选择与漫游 MS正常情况下都接入到HPLMN或者EHPLMN，但也可能选择VPLMN。比如在丢失覆盖后，（自动选网模式下）可用的最高优先级PLMN将会被选择，（手动模式下）由用户从MS指示的可用的PLMNs中选择获取服务的PLMN。 为了防止MS不断地在not allowed area尝试漫游服务，当MS得知该区域禁止漫游服务后，将对应的LA或者TA添加到list of "forbidden LAs for roaming" or "forbidden TAs for roaming"中, 该列表在
深入理解查准率与查全率
热门推荐
zhaojianting的博客
03-22 2万+
欢迎关注鄙人公众号，技术干货随时看！ 查准率与查全率（召回率）是在信息检索与机器学习领域常用的衡量指标，书籍或网络上有非常多的定义，敝人在重温周志华先生《机器学习》一书时，发现书中的定义是非常科学的，而且把这两个指标根据混淆矩阵进行了数学化！ 在信息检索领域的定义：检索出的n个文档中正确结果（m个结果是正确的）的比率即是查准率，即p=n/m；假设在所有的数据中正确结果有n个，一次检索...
由ERA5逐小时数据获取逐日数据——三种方法
SAREbenjamin的博客
08-13 1万+
关于欧洲中心的ERA5数据，最令人想吐槽的就是官方未给出方便快捷的逐日数据下载，需要自己下载逐小时数据再进行处理，本篇文章共总结了3种方法，通过ERA5逐小时数据获取逐日数据： 一、官方提供的在线数据计算 网址：https://cds.climate.copernicus.eu/cdsapp#!/software/app-c3s-daily-era5-statistics?tab=app 该方法的优点是指向性强，但是如果要下载多年日平均数据，可能需要每一年每一个月手动下...
PCA数据降维、数据重建计算步骤
jngwe的博客
12-11 3483
PCA基本原理：
Matlab 批量修改nc文件中某一个变量的数值
jerry003的博客
12-15 8277
场景：发现nc文件中某个变量的数值存在问题，修改在源文件基础上修改。 注意事项： 注意nc文件中，变量和属性的区别。可以单独修改变量，也可单独修改某一个变量的属性； 修改之后的变量的形状（维度）不能变，否则需要修改维度信息； 修改之后的变量的属性仍是原来的信息，若也要修改，则需针对属性进行修改； 本次以修改nc文件中的 “time” 变量为例说明 % read original nc file DirPath = '\Wind'; inFileList = dir([DirPath, '\*.nc
关于Matlab插值的问题，这些应该够用了吧
jerry003的博客
08-25 7319
Matlab插值方法
是时候放弃pcolor了
jerry003的博客
12-16 3402
pcolor .vs. imagesc
制作nc文件流程——以站点降水数据为例
jerry003的博客
03-04 3363
这里写目录标题netcdf文件制作流程1. define dimensions2. define variables3. put variables into nc4. put attribute into nc5. define global attribution6.附上代码全文 netcdf文件制作流程 【前提】 提前设计好需要存储变量的结构及维度，之类需要注意的是创建 nc 文件的 类型\color{#FF0000}{类型}类型，以下供参考： Value of cmode Descript
Array subscript for RUNOFF_DAY must have same size as source expression.
04-21
这是一个编程中的错误提示信息，通常表示您正在尝试使用一个数组的索引，但是该数组的大小与源表达式的大小不相同。可能是在您的代码中使用了错误的数组索引或者数组大小被错误地定义。请仔细检查您的代码并修复此错误。
评论 2
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

大雨海深

感谢您的支持和鼓励

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值