matlab 序列时间尺度变换程序,使用Matlab合成月尺度、年尺度数据

处理PML v2陆地蒸散发与总初级生产力数据集,由8天尺度合成月尺度、年尺度数据

处理代码

% By Yang,2020/5/9

% 所有数据全部解压放在一个文件夹下

%% 合成年尺度数据

for j = 1:18 % 2002-2019,18年

a = num2str(2001+j); % 数字转字符

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) % 计算文件个数

pml_dz2 = strcat(pml_dz(i).folder,'\',pml_dz(i).name); % 数据的每个文件地址+文件名

[pml, R] = geotiffread(pml_dz2); % 读取tif数据

gpp = gpp + double(pml(:,:,1)); % 总初级生产力累加

Ec = Ec + double(pml(:,:,2)); % 植被蒸腾累加

Es = Es + double(pml(:,:,3)); % 土壤蒸发累加

Ei = Ei + double(

  • 0
    点赞
  • 21
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值