2026/1/14 19:51:01
网站建设
项目流程
怎样做网站推广啊视频,昆山网站建设设计,网站开发有前途么,做室内效果图的网站目录
一、区域定义与地图初始化
二、数据集加载
三、年份筛选与数据计算
四、数据导出到 Google Drive
五、代码关键特性与注意事项
六、运行结果 若觉得代码对您的研究 / 项目有帮助#xff0c;欢迎点击打赏支持#xff01;需要完整代码的朋友#xff0c;打赏后可在后…目录一、区域定义与地图初始化二、数据集加载三、年份筛选与数据计算四、数据导出到 Google Drive五、代码关键特性与注意事项六、运行结果若觉得代码对您的研究 / 项目有帮助欢迎点击打赏支持需要完整代码的朋友打赏后可在后台私信复制文章标题发给我我会尽快发您完整可运行代码感谢支持本代码基于 Google Earth EngineGEE平台利用 ECMWF欧洲中期天气预报中心发布的 ERA5-LAND 月度聚合数据集自动提取 2024 年全球指定区域北纬 90° 至南纬 60°、东经 180° 至西经 180°的三类径流数据并将其导出为 GeoTIFF 格式存储到 Google Drive 中数据单位最终转换为毫米mm具体包括总径流总和runoff_sum地下径流总和sub_surface_runoff_sum地表径流总和surface_runoff_sum一、区域定义与地图初始化var geometry ee.Geometry.Polygon( [[[-180, 90], [-180, -60], [180, -60], [180, 90]]], null, false); Map.centerObject(geometry, 3);区域几何对象创建ee.Geometry.Polygon()GEE 中用于创建多边形区域的核心函数参数为顶点坐标数组、投影信息此处为null默认采用 GEE 内部投影、是否闭合false表示按顶点顺序自动闭合。顶点坐标含义四个顶点分别定义了全球范围的矩形区域 —— 西经 180° 北纬 90°北极点附近、西经 180° 南纬 60°、东经 180° 南纬 60°、东经 180° 北纬 90°覆盖了除南极洲大部分区域外的全球陆地与海洋范围。地图居中设置Map.centerObject(geometry, 3)将地图视图中心定位到上述定义的区域第二个参数3为地图缩放级别GEE 缩放级别 1-20数值越小视野越广3 级可显示全球全貌方便在 GEE 编辑器中可视化查看研究区域。二、数据集加载var runoff ee.ImageCollection(ECMWF/ERA5_LAND/MONTHLY_AGGR);ee.ImageCollection()GEE 中用于加载影像集合多幅影像组成的数据集的函数参数为数据集的 GEE 内置 ID。数据集说明ECMWF/ERA5_LAND/MONTHLY_AGGR是 ERA5-LAND 数据集的月度聚合版本与原始逐小时数据相比已提前完成月度尺度的统计聚合包含降水、蒸发、径流等多个水文气象变量适合长时间尺度的区域水文分析。三、年份筛选与数据计算for (var year 2024; year 2024; year) { var yearly runoff.filterDate(year -01-01, year -12-31); // 总径流计算关键步骤示例 var runoffSum yearly.select(runoff_sum).sum().multiply(1000).rename(runoff_sum_ year).clip(geometry); // 地下径流、地表径流计算逻辑与总径流一致仅筛选的变量名不同 }年份循环此处为单年份循环2024 年通过for循环限定数据提取的时间范围若需扩展到多年如 2000-2024 年可修改循环条件为var year 2000; year 2024; year。时间筛选yearly runoff.filterDate(...)对加载的月度影像集合按时间筛选提取指定年份 1 月 1 日至 12 月 31 日的所有月度影像共 12 幅对应全年 12 个月。径流数据计算核心步骤select(runoff_sum)从月度影像集合中筛选出 “总径流总和” 对应的波段ERA5-LAND 数据集中不同变量以不同波段存储runoff_sum为总径流波段名sub_surface_runoff_sum为地下径流波段名surface_runoff_sum为地表径流波段名。sum()对筛选后的 12 幅月度影像进行像素级求和运算得到全年的径流总和原始单位为米m。multiply(1000)单位转换将米m转换为毫米mm1m1000mm符合水文数据常用单位习惯。rename(...)为计算后的影像波段重命名添加年份后缀如runoff_sum_2024避免多年份数据导出时波段名冲突。clip(geometry)按之前定义的geometry区域对影像进行裁剪只保留目标区域内的数据去除区域外的无效值减少数据存储量。四、数据导出到 Google Drive以总径流数据导出为例地下径流、地表径流导出逻辑完全一致仅参数细节不同Export.image.toDrive({ image: runoffSum, description: year _runoff_sum_mm, fileNamePrefix: year _runoff_sum_mm, folder: ERA5L-runoff-mm, region: geometry, scale: 5000, crs: EPSG:4326, maxPixels: 1e13 });导出函数说明Export.image.toDrive()GEE 中用于将单幅影像导出到 Google Drive 的核心函数参数以键值对形式传入定义导出的各项配置。关键参数解析image指定待导出的影像对象此处为计算后的全年总径流影像runoff_sum。description导出任务的描述名称在 GEE 任务列表中显示便于区分不同任务。fileNamePrefix导出文件的前缀名最终导出的 GeoTIFF 文件名为 “前缀名.tif”如2024_runoff_sum_mm.tif。folderGoogle Drive 中存储文件的文件夹名称需提前在 Google Drive 中创建或由导出任务自动创建。region导出的数据区域与之前裁剪的geometry保持一致确保导出范围准确。scale导出影像的空间分辨率单位为米m此处5000表示影像像素大小为 5000m×5000m约 25 平方公里 / 像素ERA5-LAND 数据集的原生空间分辨率为 0.1°约 11km此处设置为 5km 是平衡分辨率与数据量的选择。crs导出影像的投影坐标系EPSG:4326为 WGS84 地理坐标系经纬度投影是全球通用的地理坐标系便于后续在 ArcGIS、QGIS 等软件中处理。maxPixels允许导出的最大像素数量GEE 默认限制导出像素数为 1e1310 万亿此处设置为 1e13 是为了避免因全球范围数据量过大导致导出失败5km 分辨率的全球区域像素数约为 1e8 级别远小于 1e13。五、代码关键特性与注意事项关键特性自动化程度高通过循环与筛选逻辑可快速批量提取指定年份的三类径流数据无需手动逐幅影像处理。数据标准化统一完成单位转换m→mm、投影设置EPSG:4326和区域裁剪导出数据可直接用于后续水文模型、空间分析等工作。可扩展性强修改循环年份范围可扩展到多年数据提取修改geometry参数可调整研究区域如某一流域、某一国家修改scale参数可调整导出分辨率。注意事项GEE 权限与资源使用前需拥有 Google 账号并登录 GEE 平台导出全球范围数据时需确保 GEE 账户有足够的计算资源批量导出多年数据时建议分批次执行避免任务队列拥堵。文件夹创建需提前在 Google Drive 中创建名为ERA5L-runoff-mm的文件夹否则导出文件会存储在 Google Drive 根目录不利于文件管理。数据精度说明ERA5-LAND 是再分析数据集基于观测数据与数值模型融合生成径流数据为模拟结果使用时需结合实地观测数据进行验证尤其适用于大尺度区域分析小尺度区域需谨慎使用。变量名准确性ERA5-LAND 数据集的波段名如runoff_sum是固定的修改时需参考 GEE 数据集的官方文档避免因波段名错误导致数据提取失败。六、运行结果点击RUN即可下载数据对下载的数据进行后处理并可视化出图若觉得代码对您的研究 / 项目有帮助欢迎点击打赏支持需要完整代码的朋友打赏后可在后台私信复制文章标题发给我我会尽快发您完整可运行代码感谢支持