1km分辨率逐月降雨量和最高温度数据集(1901-2022)--数据处理

1km分辨率逐月降雨量和最高温度数据集(1901-2022)的下载可以参考我的另外一篇博客:

这里的温度和降雨数据集都是NC格式的,需要将其处理为tif格式,我采用的处理软件是MATLAB。

本篇博客以处理温度数据为例,下面是完整的处理代码,注意看代码中的注释。

Matlab 复制代码
clc
clear all
%% 批读取NC文件的准备工作
datadir = 'E:\GSH\1Km温度数据\1985-2017_3\'; %指定批量数据所在的文件夹
filelist = dir([datadir,'*.nc']);    %列出所有满足指定类型的文件
% a = filelist(1).name               %查看要读取的文件的编号
% b = filelist(2).name
k=length(filelist);
for i = 1:11 %依次读取并处理,按照文件夹内文件的顺序读取
    
    %% 批量读取NC文件
    ncFilePath = ['E:\GSH\1Km温度数据\1985-2017_3\',filelist(i).name]; %设定NC路径
    num_1 = filelist(i).name(1:8); %读取数据编号,以便于保存时以此编号储存tif
    num_2=strcat('tmx_',num2str(str2num(num_1(5:8))+1 ));
    num_3=strcat('tmx_',num2str(str2num(num_1(5:8))+2 ));
    
    %% 读取变量值
    %根据ncdisp函数读取到的nc文件变量相应替换
    lon=ncread(ncFilePath,'lon');     %读取经度信息(范围、精度)
    lat=ncread(ncFilePath,'lat');      %读取维度信息
    time=ncread(ncFilePath,'time');         %读取时间序列
    pre=ncread(ncFilePath,'tmx');   %获取温度数据
    %一个文件包含三年数据,将这三年的数据单独导出
    time_1=time(1:12);
    time_2=time(13:24);
    time_3=time(25:36);
    pre_1=pre(:,:,1:12);
    pre_2=pre(:,:,13:24);
    pre_3=pre(:,:,25:36);
    sum_pre_1=sum(pre_1,3)/12;                %求第一年平均温度
    sum_pre_2=sum(pre_2,3)/12;                %求第二年平均温度
    sum_pre_3=sum(pre_3,3)/12;                %求第三年平均温度
    %% 存为tif格式
    %第一年
    data_1=rot90(sum_pre_1,1);   %逆时针旋转90°,不旋转的话最后的图像朝向是错的
    R = georasterref('RasterSize', size(data_1),'Latlim', [double(min(lat)) double(max(lat))], 'Lonlim', [double(min(lon)) double(max(lon))]);
    geotiffwrite(['E:\GSH\1Km温度数据\1985-2017_3\',num_1,'.tif'],data_1,R);
    disp([num_1,'done'])
    %第二年
    data_2=rot90(sum_pre_2,1);   %逆时针旋转90°,不旋转的话最后的图像朝向是错的
    R = georasterref('RasterSize', size(data_2),'Latlim', [double(min(lat)) double(max(lat))], 'Lonlim', [double(min(lon)) double(max(lon))]);
    geotiffwrite(['E:\GSH\1Km温度数据\1985-2017_3\',num_2,'.tif'],data_2,R);
    disp([num_2,'done'])
    %第三年
    data_3=rot90(sum_pre_3,1);   %逆时针旋转90°,不旋转的话最后的图像朝向是错的
    R = georasterref('RasterSize', size(data_3),'Latlim', [double(min(lat)) double(max(lat))], 'Lonlim', [double(min(lon)) double(max(lon))]);
    geotiffwrite(['E:\GSH\1Km温度数据\1985-2017_3\',num_3,'.tif'],data_3,R);
    disp([num_3,'done'])
    
end
disp('Done!')
相关推荐
晨非辰1 小时前
#C语言——刷题攻略:牛客编程入门训练(十一):攻克 循环控制(三),轻松拿捏!
c语言·开发语言·经验分享·学习·visual studio
Luna-player8 小时前
快速开发一类似个人网站空间的工具使用什么方式比较好,比如网页或者个Windows程序,并且使用什么技术开发比较好,区别优势局限性,分别说明一下
经验分享
阿豪38 小时前
2025 年职场转行突围:除实习外,这些硬核证书让你的简历脱颖而出(纯经验分享)
大数据·人工智能·经验分享·科技·信息可视化·产品经理
iGarment12 小时前
服装采购跟单系统的高效管理实践
大数据·经验分享·云计算
Luna-player13 小时前
unity开发类似个人网站空间
经验分享
汇能感知14 小时前
光谱相机在AI眼镜领域中的应用
经验分享·笔记·科技
汇能感知14 小时前
光谱相机的图像模式
经验分享·笔记·科技
赵谨言20 小时前
基于数据挖掘的单纯冠心病与冠心病合并糖尿病的证治规律对比研究
经验分享·数据挖掘·毕业设计
赵谨言20 小时前
基于大数据挖掘的药品不良反应知识整合与利用研究
经验分享·数据挖掘·毕业设计
月阳羊20 小时前
【硬件-笔试面试题-81】硬件/电子工程师,笔试面试题(知识点:详细讲讲同步时钟与异步时钟通信)
java·经验分享·单片机·嵌入式硬件·面试