xuguangning1218/Typhoon_Path

关于ERA-Interim的问题

Opened this issue · 4 comments

Nianqitongs commented

作者您好，想请问一下对于ERA-Interim数据的处理方法，就是怎么得到和CMA轨迹对应经纬位置的ERA数据，我想用ERA1979年到2000年的数据，但是不知道该怎么处理ERA数据使他能够匹配CMA 1979到2000年的对应经纬。作者提供的数据ERA是与CMA里的台风轨迹位置对应的，我想问一下这是怎么处理的。谢谢作者

xuguangning1218 commented

具体步骤如下：
1、收集好CMA数据中台风中心的经纬度，由于CMA的台风中心经纬度是0.1度的，因此需要通过除以10转化为正常的经纬度
2、由于ERA-Interim目前已经暂停更新，因此建议使用CDS中的ERA5 Pressure Level数据集，里面包含了ERA-Interim的所有数据及相应更新数据
3、一般CDS数据下载有两种方式：网站下载和代码下载。其中网站下载无法选择相应的grid精度，建议使用代码下载
4、注册CDS，https://cds.climate.copernicus.eu/user/register?destination=%2F%23!%2Fhome
5、您可以登录帐号后，查看这个链接https://cds.climate.copernicus.eu/api-how-to配置您的机器，以获得代码下载权限
6、如果您不熟悉怎么配置下载参数，可以访问链接，通过网页选择您要下载的数据然后点击最下面的“Show API Request”按钮可以查看对应的生成的下载代码（通过代码可以增加grid参数选择对应的经纬度分辨率，如grid:[0.5,0.5]即表示经纬度为0.5度0.5度的空间分辨率，默认为0.25）或通过这个链接学习使用。下面为一个示例代码：
import cdsapi
c = cdsapi.Client()
c.retrieve(
'reanalysis-era5-pressure-levels',
{
'product_type': 'reanalysis',
'format': 'grib',
'variable': 'geopotential',
'pressure_level': '1000',
'year': '1979',
'month': '02',
'day': '02',
'time': [
'00:00', '01:00', '02:00',
],
'area': [
20, 80, -20,
120,
],
},
'download.grib')
7、了解如何使用代码下载CDS数据后，您可以通过第1步寻找到的台风中心经纬度，生成您想要下载的再分析数据的范围、对应时间等，并填入area及相应参数中，以实现个性化下载。示例代码中的area表示为一个20S20N,80E100E的4040的范围，其中经纬度分辨率为默认的0.25*0.25
8、关于下载的小技巧，建议您不要一次性下载所有需要的数据，分批逐个下载最佳，因为下载需要构建，生成下载。而且生成下载有效期，构建太多数据可能有效期内无法下载完成

Nianqitongs commented
xuguangning1218 commented

作者您好，关于这个代码示例我还有个疑问，在选取下载数据的经纬度时无法跨越太平洋上180度经度线，是否如果要下载180度经度线两边的数据需要分开下载。比如说下面这张图上代表的代码范围就是东经95度到东经180度，不知道我理解的对不对。 念 @.***  
…

[是的，如果您因为选择的格点很大，出现这个情况，则需要下载两份再拼接再一起，如下图所示：
WechatIMG27
WechatIMG28

Nianqitongs commented