免费的气象数据（GLDAS、ERA5）下载步骤及数据处理的详细过程

谁践踏了优雅 2022-09-10 00:21 1593阅读 0赞

**需要气象数据的小伙伴们看过来呦！首先声明，能够免费下载的气象数据一般分辨率比较低哟！**  
**一、GLDAS数据**  
**1.下载地址：**（[GLADS下载地址][GLADS]）GLADS数据分有3h、一天、一个月的产品，大家可以根据自己的需求进行下载。

**2.下载步骤**：选择完数据产品、研究区范围、日期，将其生成下载的链接，网址会直接将链接放在文本里，用户可以进行下载。  
（1）首先，比较容易地下载步骤是将文本导入火狐浏览器中下载管理器中，具体可以参考——[参考文章][Link 1]，此方法的缺点是一般不能链接到网址，下载总是出错。  
（2）采用wegt进行下载，下载的步骤特别推荐参考文章——[wegt下载GLDAS数据][wegt_GLDAS]，但是此篇文章的最后输入命令那点，可能说的有些含糊，所以我在此将其命令说一下：  
`wget --load-cookies D\.urs_cookies --save-cookies D:\.urs_cookies --auth-no-challenge=on --keep-session-cookies --user=<*****> --ask-password --content-disposition <你的文本名称>`  
**大家要注意，直接将用户名和密码输入进去即可，不需要如参考文章中的命令步骤。还有最重要的是不能将用户名的<>也作为命令输入进去（不要犯低级错误）。还有，要将所有生成的文件都放在一个盘下面，不需要单独建立文件夹了。**

**3.数据处理**：由于下载的数据都是NC格式，要将其进行NC到TIF的转换。我采用的matlab软件操作的，代码如下：

clc;
    clear;
    
    %% 批读取NC文件的准备工作
    datadir = 'E:\yz\YZ-Date\ERA\tenperature\2013\'; %指定批量数据所在的文件夹
    filelist = dir([datadir,'*.nc']);       %列出所有满足指定类型的文件
    
    k = length(filelist);
    
    for i = 1:k  %依次读取并处理
        %% 批量读取NC文件
        ncFilePath = ['E:\yz\YZ-Date\ERA\tenperature\2013\',filelist(i).name]; %设定NC路径
        outFileName = filelist(i).name(1:7); %输出文件名，(1:33)即为第1-33个字符
          
        %% 读取变量值
         ncinf = ncinfo(ncFilePath);
        lon = ncread(ncFilePath,'longitude');          %读取经度信息（范围、精度）
        lat = ncread(ncFilePath,'latitude');          %读取维度信息
        time = ncread(ncFilePath,'time');        %读取时间序列
         pre = ncread(ncFilePath,'t2m');     %获取“温度”数据
        % pre = ncread(ncFilePath,'Rainf_f_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))]);
        R = georasterref('RasterSize', size(b),'Latlim', [double(min(lat)) double(max(lat))], 'Lonlim', [double(min(lon)) double(max(lon))]);
        geotiffwrite(['E:\yz\YZ-Date\ERA\tenperature\2013\',outFileName,'_temperature.tif'], b, R);   % 输出路径
        disp([outFileName,'done'])
        
    end
    disp('finish!')

**程序中涉及到镜像翻转的问题，我给出以下的代码，供大家可以根据不同的影像情况进行修正：**

% 逆时针旋转90°
    MonthData = rot90(MonthData);
    % 上下翻转
    MonthData = flipud（MonthData）
    % 左右翻转
    MonthData = fliplr（MonthData）

**最重要的是，你需要哪个变量（温度、降水等），要知道其属性名称，每个变量对应的属性名称如下**  
![在这里插入图片描述][watermark_type_ZHJvaWRzYW5zZmFsbGJhY2s_shadow_50_text_Q1NETiBA5L2O6LCD55qE5aSn6ICz5py15Zu-5Zu-_size_20_color_FFFFFF_t_70_g_se_x_16]  
**二、ERA5数据**  
**1.下载地址**：ERA5下载（[ERAS下载地址][ERAS]），此数据下载起来比较方便，进入网站就会下载了。但是注意一点，一定要一个月份或者一天的提交订单，否则同时选中好几个月的提交订单，所有的NC文件就会合成一个NC文件，目前这个问题我还没有解决。

**2.数据的处理**：与GLDAS数据处理的方法一样（NC-TIF），以上代码可用。那么每个变量的属性名称怎么知道呢？教大家一个方法，就是将下载的NC数据在ArcGIS中显示，看影像的名称即为属性名称！

三、**总结**  
气象产品就是拿来用着方便，但是分辨率太低了，限制了研究的深度。所以要想获取高分辨率的气象数据，就要学会插值法，获得想要的分辨率的气象数据，但是缺点数值可能不太准确。希望以上对大家有所帮助！

长路漫漫  
唯有坚持…

[GLADS]: https://disc.gsfc.nasa.gov/datasets?keywords=GLDAS&page=1
[Link 1]: https://articles.zsxq.com/id_v1n3k2wjlw7x.html
[wegt_GLDAS]: https://blog.csdn.net/weixin_44035267/article/details/111221007
[watermark_type_ZHJvaWRzYW5zZmFsbGJhY2s_shadow_50_text_Q1NETiBA5L2O6LCD55qE5aSn6ICz5py15Zu-5Zu-_size_20_color_FFFFFF_t_70_g_se_x_16]: /images/20220829/9d3df2fc0fb74cf8843958810baa5cf5.png
[ERAS]: https://cds.climate.copernicus.eu/cdsapp#!/home