matlab批量将GPM IMERGM V06月降水数据的nc4格式转成tif文件

最新推荐文章于 2023-05-18 22:42:42 发布

会掉头发的芒果

最新推荐文章于 2023-05-18 22:42:42 发布

阅读量1.5k

点赞数 3

文章标签：经验分享 matlab Powered by 金山文档

本文链接：https://blog.csdn.net/weixin_46421446/article/details/129733500

版权

该文介绍如何从GPMIMERGMV06数据集中下载降水数据.nc4文件，并在Matlab环境中读取文件信息，然后通过ncread和geotiffwrite函数将这些.nc4文件批量转换为tif地理图像格式。

摘要由CSDN通过智能技术生成

数据下载链接：GPM IMERG Final Precipitation L3 1 month 0.1 degree x 0.1 degree V06 (GPM_3IMERGM)

1、这里下载的是GPM IMERGM V06数据集中的降水数据。

2、下载下来的是.nc4文件。在matlab中输入下面代码可以获取文件的变量信息。

%% 读取nc文件信息（变量信息）
ncdisp('D:\shiyanshuju\GPM\20010101.nc4')

3、获取变量信息后输入下面代码将.nc4格式转为tif格式。

clc;
clear;

%% 批读取NC文件的准备工作
datadir = 'D:\data\GPM2001-2020\';       %指定数据所在的文件夹
filelist = dir([datadir,'*.nc4']);       %列出所有满足指定类型的文件，数据为.nc4类型
% a = filelist(1).name;                  %查看要读取的文件的编号
k=length(filelist);                      %读取文件夹中.nc4文件个数

for i = 1:k  %依次读取并处理
    
    %% 批量读取NC文件
    ncFilePath = ['D:\data\GPM2001-2020\',filelist(i).name]; %设定NC路径
    name = filelist(i).name(); %读取数据编号，以便于保存时以此编号储存tif，括号里写上（1:6）可以表示读取文件1到6的名称
      
    %% 读取变量值
    lon=ncread(ncFilePath,'lon');          %读取经度信息（范围、精度）
    lat=ncread(ncFilePath,'lat');          %读取维度信息
    time=ncread(ncFilePath,'time');        %读取时间序列
    pre=ncread(ncFilePath,'precipitation');%获取降雨变量数据
   
    %% 存为tif格式
    R = georasterref('RasterSize', size(pre),'Latlim', [double(min(lat)) double(max(lat))], 'Lonlim', [double(min(lon)) double(max(lon))]);
    geotiffwrite(['D:\data\GPM2001-2020\nc转tif\',name,'.tif'],pre,R); %输出路径，名称
    disp([name,'done'])
    
end
disp('finish!')