3
\ begingroup美元

我有下载ERA5 netCDF文件为了得到降水每月/每日房价为11年。但我已经意识到,没有选项来选择我感兴趣的领域。似乎这是一个很大的文件来处理。

但是,我目前使用R语言,但仍然在基本水平。

我的问题是关于如果有提取一个特定区域的降水值吗?

\ endgroup美元
3
  • \ begingroup美元 请编辑你的标题反映实际问题。回答,几乎肯定会有一条路在R阅读只是netCDF文件的一部分,但我没有足够的R用户告诉你它是什么。希望有人会过来,反潜战。 \ endgroup美元
    - - - - - -半日西蒙
    2019年8月28日20:37
  • \ begingroup美元 @SemidiurnalSimon,我做了一个编辑它看起来如何? \ endgroup美元
    - - - - - -uhoh
    2019年8月28日20:58
  • \ begingroup美元 我建议类似“如何netCDF文件的只读部分R”或类似的。 \ endgroup美元
    - - - - - -半日西蒙
    2019年8月28日21:26

1回答1

2
\ begingroup美元

有两种方法来做到这一点:

1——你可以选择感兴趣的特定领域,当下载文件。这是高度优先当你需要下载大量的变量和年。ERA5 netCDF提供在两个不同的格式:每小时的时间步骤,或者每月平均(见)。如果你需要一个sub-monthly时间步长时期全球许多变量文件非常重的磁盘空间。你可以设置一个python脚本下载后只选择域的数据程序的建议里特•Stauffer的好文章(当然还有其他的方法可以做到这一点)。

2 -您可以下载的全球netCDF变量和年在考虑广告使用软件来提取您感兴趣的领域。的一个最有效的系统就是使用气候数据操作符CDO。你可以用sellonlatbox这样的一个命令

cdo sellonlatbox、LON1 LON2, LAT1 LAT2 Input_file_Name。数控Output_file_Name.nc

或者你可以使用R至于你的情况下,做一些像:

库(光栅)库(rgdal) #加载您感兴趣的区域setwd (“path_domain_shapefile”) domain_shp < -readOGR (shapefile.shp) #读netCDF setwd (“path_to_ERA5_files_folder”) pr_data <堆栈(ERA5_file.nc) #提取数据为你感兴趣的领域pr_data_domain <提取(pr_data domain_shp) #你可以改变这一数据帧和把它写成一个csv

空间数据的基本操作提供了参考依据R我建议介绍指南

希望这个有帮助。

\ endgroup美元

    你的答案

    通过点击“发布你的答案”,你同意我们服务条款并承认您已阅读并理解我们的隐私政策的行为准则

    不是你要找的答案?浏览其他问题标记问你自己的问题