處理PML v2陸地蒸散發(fā)與總初級生產(chǎn)力數(shù)據(jù)集，由8天尺度合成月尺度、年尺度數(shù)據(jù)

處理代碼

% By Yang,2020/5/9
% 所有數(shù)據(jù)全部解壓放在一個文件夾下
%% 合成年尺度數(shù)據(jù)
for j = 1:18                                                               % 2002-2019，18年
         a   = num2str(2001+j);                                            % 數(shù)字轉(zhuǎn)字符
         dz  = ['F:\Regional_Dataset\PMLV2_005deg\PMLV2_005deg\*',a,'*.tif'];
         pml_dz = dir(dz);
         gpp = zeros(3000,7200);                                   % 生成uint16的空矩陣
         Ec  = zeros(3000,7200);
         Es  = zeros(3000,7200);
         Ei  = zeros(3000,7200);
         Ew  = zeros(3000,7200);
        parfor i = 1:size(pml_dz,1)                                        % 計算文件個數(shù)
              pml_dz2  = strcat(pml_dz(i).folder,'\',pml_dz(i).name);      % 數(shù)據(jù)的每個文件地址+文件名   
              [pml, R] = geotiffread(pml_dz2);                             % 讀取tif數(shù)據(jù)                        
              gpp      = gpp + double(pml(:,:,1));                                 % 總初級生產(chǎn)力累加
              Ec       = Ec + double(pml(:,:,2));                                  % 植被蒸騰累加
              Es       = Es + double(pml(:,:,3));                                  % 土壤蒸發(fā)累加
              Ei       = Ei + double(pml(:,:,4));                                  % 冠層截留蒸發(fā)累加
              Ew       = Ew + double(pml(:,:,5));                                  % 水體、積雪蒸發(fā)累加           
        end
              gpp(gpp==0)  = NaN;                                                    % 空值設(shè)為NaN
              Ec(Ec==0)      = NaN;
              Es(Es==0)      = NaN;
              Ei(Ei==0)       = NaN;
              Ew(Ew==0)      = NaN;
              PML_y{1,j}      = gpp*8*0.01;         % 真實值=像素值*scale_factor （scale_factor=0.01)
              PML_y{2,j}      = Ec*8*0.01;           % 每個數(shù)據(jù)對應(yīng)的是這個8天里GPP和ET要素的平均
              PML_y{3,j}      = Es*8*0.01;
              PML_y{4,j}      = Ei*8*0.01;
              PML_y{5,j}      = Ew*8*0.01;  
end     
save PML_y.mat PML_y -v7.3   

數(shù)據(jù)網(wǎng)站

國家青藏高原科學(xué)數(shù)據(jù)中心

補充

月尺度合成比年尺度還稍微要麻煩一點，其實應(yīng)該先計算月尺度的再合成年尺度的更為準(zhǔn)確些。