探索Python融合地学:一文教会你下载ERA5-Land数据

探索Python融合地学:一文教会你下载ERA5-Land数据

再分析数据在气象学领域用的比较多,下载数据有的时候还挺头疼的。今天小编教你下载ERA5-Land数据的三种方式。话不多说,咱们学起来!

一、官网下载

贴出官网:https://cds.climate.copernicus.eu/datasets/reanalysis-era5-land?tab=download

相较于ERA5数据,Land数据空间分辨率提高了,为0.1°格网,时间分辨率同样是逐小时,可以说已经很方便大家开展各项研究了。

第一步:注册个人账号,在右上角的人像这里点击可以注册个人的账号,这里不多说。

第二步:点击“Download”,进入下载界面,在下载界面你可以看到很多气象变量,我们随便选择一种气象要素,比如风的U/V分量,支持多选。

第三步:选择研究时段,这里可以选择某年某个月的所有日期的数据。bug在于,它只能一个月一个月申请下载,如果你研究时间尺度很长,每一年得点12下,下载12份文件,再下载下一年的。所以,如果你是一个体验者,推荐用这种方法。如果你是一个研究者,想要快读完成数据下载,建议用后面小编罗列的代码方法。

第四步:确定研究区域,如果你的研究区域是全区,选择第一个即可,不用更改,但大多数人都是以小区域研究为主,所以可以选择“Sub-Region”,填写上经纬度的范围坐标。

第五步:确定数据格式,一般选择NetCDF4即NC文件。

第六步:确定提交申请,点击“Submit form”,再点击右上角的“Your Requests”,进去后可以看到你刚刚提交的申请链接,等读取后,点击下载即可。此过程需要一些时间,还望大家见谅。

以上是ERA5-Land数据的官网下载方式,ERA5系列数据采用这种方式也可以下载,而且也挺方便的,因为ERA5数据可以按照年份直接打包12个月的数据,不需要一个月一个月单独下载。

二、Python代码下载

如果你需要多年多月多小时的数据,那么我相信第二种方法会很适合你。

工具:配置Python环境,这里我用的是PyCharm。

第一步:一样的,先打开下载界面,找到“Terms of use”,点击“Accepted”,不然下载不了。

第二步:点击右上角头像,弹出窗口中有“Your profile”,点击进入

第三步:在界面中下滑,找到“API key”,找到ur1和key,点击右侧的“复制”符号。

第四步:选择C盘,点击用户→用户名,在里面新建一个txt文件,命名为“.cdsapirc”。将第三步复制的链接和密钥粘贴到文本里面,保存下来,这个文件非常有用。

第五步:进入编译器,搜索cdsapi包,安装最新版。

第六步:贴上代码,想要的拿去。

import cdsapi import os from pathlib import Path # 创建保存目录 save_dir = Path(r"H:\Wind_U_V") save_dir.mkdir(parents=True, exist_ok=True) # ERA5-Land数据集名称 dataset = "reanalysis-era5-land" # 初始化CDS API客户端 client = cdsapi.Client() # 循环处理2023年每个月的数据 for month in range(1, 13): # 将月份格式化为两位数字符串 month_str = f"{month:02d}" # 根据月份确定天数(考虑不同月份的天数差异) if month in [1, 3, 5, 7, 8, 10, 12]: days = [f"{day:02d}" for day in range(1, 32)] elif month in [4, 6, 9, 11]: days = [f"{day:02d}" for day in range(1, 31)] else: # 2月(2023年不是闰年) days = [f"{day:02d}" for day in range(1, 29)] # 构建请求参数 request = { "variable": [ "10m_u_component_of_wind", "10m_v_component_of_wind" ], "year": "2023", "month": month_str, "day": days, "time": [ "00:00", "01:00", "02:00", "03:00", "04:00", "05:00", "06:00", "07:00", "08:00", "09:00", "10:00", "11:00", "12:00", "13:00", "14:00", "15:00", "16:00", "17:00", "18:00", "19:00", "20:00", "21:00", "22:00", "23:00" ], "data_format": "netcdf", "download_format": "unarchived", "area": [54, 72, 4, 137] } # 设置输出文件名 output_filename = save_dir / f"ERA5_Land_wind_2023_{month_str}.nc" print(f"正在下载2023年{month_str}月的数据...") try: # 下载数据 client.retrieve(dataset, request).download(output_filename) print(f"成功下载: {output_filename}") except Exception as e: print(f"下载2023年{month_str}月数据时出错: {e}") print("所有月份数据下载完成!")

这段代码能够实现对单年逐月逐小时ERA5-Land双变量数据的下载。

运行后,控制台状态如下:

三、Google Earth Engine平台下载

官网贴上:https://earthengine.google.com/

第一步:在GEE平台上找到数据集,这里可以查看数据集的基本信息,调用代码,波段信息等。

我们的目标是下载1970-2025年逐小时的风速数据,所以需要先选择u、v两个波段,按照风速的公式计算,完成后再导出。

非常喜欢GEE,这样我们就不用下载两个变量了。

第二步,打开GEE的编译器,贴上代码,这里我是按照月度数据合成的。

代码如下:

// 导入ERA5 Land数据集 var dataset = ee.ImageCollection('ECMWF/ERA5_LAND/HOURLY'); // 设置时间范围(1970年到2025年) var start_date = '1970-01-01'; var end_date = '2025-12-31'; // 选择u和v风速分量 var u_component = dataset.select('u_component_of_wind_10m'); // 10m u-component of wind var v_component = dataset.select('v_component_of_wind_10m'); // 10m v-component of wind // 计算风速(m/s) // 使用 map() 遍历每个图像并对每个图像应用平方、相加、开方操作 var wind_speed = u_component.map(function(image) { var u_squared = image.pow(2); // 计算u分量的平方 var v_squared = v_component.filterDate(image.date()).first().pow(2); // 同样计算v分量的平方 return u_squared.add(v_squared).sqrt(); // 合并并计算风速 }); // 定义中国范围,坐标范围为:东经73.4°~135°,北纬18°~53° var china_region = ee.Geometry.Polygon([ [[73.4, 18], [135, 18], [135, 53], [73.4, 53]] ]); // 计算1970年到2025年的逐月逐小时风速数据 var years = ee.List.sequence(1970, 2025); // 按年和月合成风速数据 years.getInfo().forEach(function(year) { var start = ee.Date.fromYMD(year, 1, 1); // 当前年份的起始日期 var end = ee.Date.fromYMD(year, 12, 31); // 当前年份的结束日期 // 按月份合成(逐月合成) for (var month = 1; month <= 12; month++) { var monthly_wind_speed = wind_speed .filterDate(start.advance(month-1, 'month'), start.advance(month, 'month')) // 过滤出每月的数据 .mean() // 计算该月的风速均值 .clip(china_region); // 限定区域为中国 // 导出每月的风速数据 Export.image.toDrive({ image: monthly_wind_speed, description: 'wind_speed_' + year + '_' + month, scale: 1000, // 设置分辨率,ERA5 Land数据分辨率是31km,设置为1000m region: china_region, fileFormat: 'GeoTIFF', maxPixels: 1e8 // 限制最大像素数 }); } }); // 计算逐小时风速数据并导出 var hourly_wind_speed = wind_speed .filterDate(start_date, end_date) .filterBounds(china_region) .map(function(image) { var hour = ee.Date(image.get('system:time_start')).get('hour'); // 获取时间中的小时 return image.set('hour', hour); // 为每个影像添加小时标签 }); // 导出逐小时数据(每小时导出为一个GeoTIFF文件) hourly_wind_speed.evaluate(function(result) { result.features.forEach(function(feature) { var hour = feature.properties.hour; var image = ee.Image(feature.id); // 获取对应的影像 var file_name = 'wind_speed_' + ee.Date(feature.properties['system:time_start']).format('YYYYMMdd_HH').getInfo(); // 导出每小时的风速数据 Export.image.toDrive({ image: image, description: file_name, scale: 1000, region: china_region, fileFormat: 'GeoTIFF', maxPixels: 1e8 }); }); }); 

综上,如果需要效率比较高的话,我推荐第二种方式,当然,如果你是GEE老用户,用第三种方式也行。祝大家下载数据愉快!

还有一点需要提醒,如果你是要做中国区域的分析,记得把下载出来的数据转成中国区域的时间再做分析哈,ERA5数据用的是UTC时间!!要加8小时才能变成中国时间!不然就是个big mistake!

今天的内容就分享到这里,希望对大家有帮助!

Read more

嘿嘿 解决了Dev C++ 中文乱码(有效版)

嘿嘿 解决了Dev C++ 中文乱码(有效版)

这是博主第一篇博客!记录一下博主的小小小小解决史! 很早就下载用了Dev c++ ,但现在隔了很长时间没去用过了再次打开发现出现中文乱码的现象!在网站上翻阅了许久!终于解决了问题!困扰了许久! ——————————————————————— 这个中文乱码看着是真烦得慌!!! tips:不要急不要急,事情慢慢都能解决掉滴! 还有不要保存在C盘哦!最好都保存在D盘内!本博客示范的未命名1.c 保存于C盘桌面上是为了演示方便! ———————————————————————————  图1 这是我们原来出现中文乱码的界面 编译的时候会出现这个窗口   图一 (再说一遍!这个中文乱码在之前没解决掉问题的时候一看到这个就很烦! ) 图二是编译过后(中文乱码版) 图二          ————————————————————————— 第一种方法(也是博主强推亲测有效法) ·第一步         请点击左上角<控制台界面>左上角                选中<默认值D> 图三   操作第一步     ·第二步          将下方“使用旧版本

By Ne0inhk
【 C/C++ 算法】入门动态规划-----路径问题(以练代学式)

【 C/C++ 算法】入门动态规划-----路径问题(以练代学式)

>每日激励:“不设限和自我肯定的心态:I can do all things。 — Stephen Curry” 绪论 : 本章是动态规划的第二篇,本章将开始二维的动态规划,在二维中的动态规划本质和一维的分析来说差不太多,只不过状态表示从一维变成了二维,而在二维上所能管理的状态就从一维的两个变成了二维的三个,也就是x轴,y轴,数组中的值。若没看了解过动规算法,我强烈建议先看第一篇blog,因为当你看完第一篇你就对动规基本认识了,其中也就能认识到它的五步骤分析法,这里也就不扩充说明而是直接使用了 ———————— 早关注不迷路,话不多说安全带系好,发车啦(建议电脑观看)。 路径问题🛣️ 本章主要还是在二维数组中的进行的动态规划: 同样还是五步走:状态表示、状态方程、初始化、移动方向、返回结果 1. 其中在二维中状态表示就会和一位略有不同,不同本质一样: 从以 i 结尾.,… ==》从左上角到达 i j 位置,… 1. 当然在最后一题中发现上面这种常规方法实现不通,因为状态方程会受后面状态影响 2.

By Ne0inhk
c++好用的刷题网址(学习c++的必看系列)

c++好用的刷题网址(学习c++的必看系列)

作为学习有很多方向,我按照分类给出。 按照刷题专用,初级,高级,API参考等进行分类,建议先收藏防止找不到 一、刷题专用 以下内容供刷题使用。 1、 LeetCode - C++ Just a moment... leetcode.com/problemset/?topicSlugs=cpp 简介: 提供大量C++算法和数据结构题目,支持在线评测。 推荐理由: 面试准备必备,提升C++编码能力。 2、HackerRank - C++ https://www.hackerrank.com/domains/cpp www.hackerrank.com/domains/cpp 简介: 提供C++编程挑战和竞赛,涵盖基础到高级题目。 推荐理由:

By Ne0inhk
【算法竞赛】C/C++ 的输入输出你真的玩会了吗?

【算法竞赛】C/C++ 的输入输出你真的玩会了吗?

🔭 个人主页:散峰而望 《C语言:从基础到进阶》《编程工具的下载和使用》《C语言刷题》《算法竞赛从入门到获奖》《人工智能AI学习》《AI Agent》 愿为出海月,不做归山云 🎬博主简介 文章目录 * 前言 * 1. OJ(online judge)题目输入情况汇总 * 1.1 单组测试用例 * 1.2 多组测试用例 * 1.2.1 测试数据组数已知 * 1.2.2 测试数据组未知 * 1.2.3 特殊值结束测试数据 * 2. 输入时特殊技巧 * 2.1 含空格字符串的特殊处理方式 * 2.2 数字的特殊处理方式 * 3. scanf/printf 和

By Ne0inhk