% 【摘要】读取.grd文件;对不同分辨率数据网格插值;对不重合区域取公共交集区
%(1)读 .grd 文件
clc; clear; close all
flistobs = dir('D:\数据\UNH-GRDC\World Runoff Data\runoff_grd\*.grd');
%(2)计算月平均值
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
cmprof_mo{i-1} = rof; % 以cell格式存储
end
%(3)计算多年平均值
temp = importdata(flistobs(1).name,delimiterIn,headerlinesIn) % 注意文件在文件夹中的序号
temp.data(temp.data == -9999) = nan;
cmprof_yr = flip(temp.data);
%(4)读经纬度坐标
f