Python下载ERA5波谱数据

Never灬more

已于 2022-04-06 19:35:47 修改

阅读量1.6k

点赞数 1

CC 4.0 BY-SA版权

文章标签： python

于 2022-04-06 16:12:53 首次发布

本文链接：https://blog.csdn.net/never_moreee/article/details/123991244

本文档介绍了如何下载和处理用于验证MPI算法准确性的ERA5波谱数据，以及如何处理哨兵1SAR数据。首先，详细阐述了在ECMWF网站上创建账户并使用API下载netcdf格式的波谱数据的过程。接着，针对哨兵1数据，提到了使用eoreader库进行处理，并提供了在Windows环境下安装相关库的注意事项和步骤。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

1.下载准备

为了验证MPI算法的准确性，下载ERA5波谱数据，作为海浪的初猜谱进行输入。下载哨兵1 SAR数据，作为SAR图像谱数据。
ECMWF波谱数据链接

1.1 创建ECMWF账户

下载该数据首先要创建账户。点击下面链接，网页上方创建新账户。
ECMWF创建
ECMWF数据有两种下载方式，一是可以在网站直接请求下载，这种下载方式适用于直接存储在CDS目录中的数据。具体下载方式看以下链接
CDS在线目录数据下载
另一种下载方式是利用python，通过API对其进行请求。可以下载未在CDS在线目录中的数据。本文需要下载的波谱数据就在CDS未在线目录中。
下载过程如下：
ECMWF官方文档: 链接
注意，需要先利用python配置CDS API。配置该文件时，首先要进行注册，得到一个API和key，需要保存该信息到usr/你的账户名/.cdsapirc 文件中。
配置具体教程如：
配置官方文件
下载出现的常见错误如下：
常见错误
注意：我们需要下载的是netcdf格式，需要在API下载时，加上分辨率信息(grid)和下载格式信息(netcdf)，如下

import cdsapi
c = cdsapi.Client()
c.retrieve('reanalysis-era5-complete', {
     'class': 'ea',
     'date': '2021-01-01/to/2021-01-28',
     'direction': '1/2/3/4/5/6/7/8/9/10/11/12/13/14/15/16/17/18/19/20/21/22/23/24',
     'domain': 'g',
     'expver': '1',
     'frequency': '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',
     'param': '251.140',
     'stream': 'wave',
     'time': '00:00:00/01:00:00/02:00:00/03:00:00/04:00:00/05:00:00/06:00:00/07:00:00/08:00:00/09:00:00/10:00:00/11:00:00/12:00:00/13:00:00/14:00:00/15:00:00/16:00:00/17:00:00/18:00:00/19:00:00/20:00:00/21:00:00/22:00:00/23:00:00',
     'type': 'an',
     'area':'50/-50/-25/0',
     'grid':'0.01/0.01',
     'format': 'netcdf',
 }, 'era5_wave_spectra_2021_01.nc')

2.数据处理

对哨兵1数据进行处理。下载eoreader库，可以对哨兵1下载的数据直接进行处理。但是安装过程有很多坑。下面进行记录以避免踩雷。
首先，安装eoreader库需要geopandas rasterior ioxarray三个库。这三个库分别又有各自的依赖包，在windows直接利用pip install进行安装会出现一系列问题。因此需要到该网站链接下进行wheel形式的安装。该网站包含众多wheel包。可以直接进行下载。
安装顺序为:
GDAL → Fiona → Geopandas
Geopandas可以利用pip3 install进行安装。
随后，在安装rasterio，利用wheel安装。最后pip3 install eoreader即可。