HadISDH数据预处理

需求

对HadISDH数据进行预处理为自己需要的

代码

%% 处理HadISDH温度dat格式数据
% 原文件为一个的dat格式数据，所有数据在一起，每月数据之间会隔一行
% 结果：提取出温度Ta，相对湿度数据RH分别为一个元胞，一个月的数据为一个矩阵，从左向右排列为一行

% 加载数据
load HadISDH.dat
% 提取保存为元胞，一个月数据为一个矩阵
for i = 1:552  %1973.1-2018.12，共552个月
    HadISDH_Ta1973_2018{1,i} = HadISDH(((i-1)*37+2):i*37,:);
end
% 处理
for i = 1:552
    temp = HadISDH_Ta1973_2018{1,i};
    temp(temp<-1000) = nan;   %填充值设为nan
    temp2 = flipud(temp);     %进行上下翻转
    HadISDH_Ta1973_2018{1,i} = temp2;
end
% 显示一月的温度图，看是否正确（地图是否有旋转/翻转现象）
Taa = HadISDH_Ta1973_2018{1,1};
imshow(Taa,[]);
% 保存
save HadISDH_Ta1973_2018.mat HadISDH_Ta1973_2018
%% 处理HadISDH湿度dat格式数据
% 提取保存为元胞，一个月数据为一个矩阵
for i = 1:552  %1973.1-2018.12，共552个月
    HadISDH_RH1973_2018{1,i} = HadISDH(((i-1)*37+2):i*37,:);
end
% 处理
for i = 1:552
    temp = HadISDH_RH1973_2018{1,i};
    temp(temp<-1000) = nan;   %填充值设为nan
    temp2 = flipud(temp);     %进行上下翻转
    HadISDH_RH1973_2018{1,i} = temp2;
end
% 显示一月的图，看是否正确（地图是否有旋转/翻转现象）
RHH = HadISDH_RH1973_2018{1,1};
imshow(RHH,[]);
% 保存
save HadISDH_RH1973_2018.mat HadISDH_RH1973_2018