GLDAS数据下载及处理（使用Matlab语言进行NC转TIF）_gldas地表径流

---

前言

随着定量遥感的不断发展，不同机构生产的不同地面参数产品也越来越成熟，但不同产品的数据格式存在学习门槛，而网上找的不同经验的适用性不同，在此分享下个人处理GLDAS数据的经验。 本文主要介绍：如何使用Matlab语言进行NC4数据格式转换为GeoTIFF数据格式

---

一、GLDAS简介

1. GLDAS(Global Land Data Assimilation System) 是描述全球陆地信息的模型，该模型包含了全球降雨量、水分蒸发量、地表径流、地下径流、土壤湿度、地表积雪的分布以及温度和热流分布等数据。模型的空间分辨率有 1∘×1∘ 和 0.25∘×0.25∘，时间分辨率有 3 个小时、1 天、一个月不等。

2. GLDAS 数据可从 GES DISC(Goddard Earth Sciences Data and Information Services Center)下载 ，结合谷歌浏览器“DownThemAll!”插件食用，风味更佳。

二、处理步骤

1.NC4转为GeoTIFF

Matlab代码如下：

clc;
clear;

%% 批读取NC文件的准备工作
datadir = 'E:\CS_You\GLDAS\'; %指定批量数据所在的文件夹
filelist = dir([datadir,'*.nc4']);       %列出所有满足指定类型的文件

k = length(filelist);

for i = 1:k  %依次读取并处理
    %% 批量读取NC文件
    ncFilePath = ['E:\CS_You\GLDAS\',filelist(i).name]; %设定NC路径
    outFileName = filelist(i).name(1:33); %输出文件名，(1:33)即为第1-33个字符
      
    %% 读取变量值
    lon = ncread(ncFilePath,'lon');          %读取经度信息（范围、精度）
    lat = ncread(ncFilePath,'lat');          %读取维度信息
    time = ncread(ncFilePath,'time');        %读取时间序列
    pre = ncread(ncFilePath,'Qle_tavg');     %获取“潜热通量”变量数据

    %% 存为tif格式
    b = flipud(pre');
    bb = rot90(b,2);
    data = fliplr(bb); %镜像反转，否则栅格的南北朝向颠倒
    
    R = georasterref('RasterSize', size(data),'Latlim', [double(min(lat)) double(max(lat))], 'Lonlim', [double(min(lon)) double(max(lon))]);
    geotiffwrite(['E:\CS_You\GLDAS\TIFF\',outFileName,'-Qle_tavg.tif'], data, R);   % 输出路径
    disp([outFileName,'done'])
    
end
disp('finish!')
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31

总结

1. GLDAS提供了多种生态因子数据，可供全球生态变化研究使用，不过其分辨率较粗，可在一定程度上揭示环境变化。
2. NC4格式的气象数据处理，颇费了些功夫，网上的教程好多不好用，虽然个人的这个数据处理版本也有人看不懂，但好的方面是：提供了一个可借鉴的经验。
3. NC数据中一般包括经纬度字段，在处理时可直接使用其构建栅格数据的投影信息。其他方式：如ENVI的Glt工具，略显繁琐。不过，哪个是自己最擅长的操作，那么他就是最好的（考虑效率时除外）

后记

写博客的初衷是分享我的一些经验，同时也方便自己在其他电脑上进行数据处理。帮了很多人，但评论区小伙伴也有遇到问题的，那么：知识付费，我的时间和经验正好可以解决你的问题。