Chapter 5 快速开始

本章通过一个完整的操作示例,带领你从零开始完成一次 GHDC 建模任务:准备边界文件、填写参数、提交任务、确认邮件、等待处理,直到下载结果。

本章使用四川彭州市龙漕沟作为示例流域——面积约 \(16\ \mathrm{km}^2\) 的山地小流域,建模年份 2017–2018,气象驱动选用 GLDAS。

整个流程共五步,如图 5.1 所示:

GHDC 任务完整流程(用户视角)。图为英文原版,各步骤含义:Step 1 提交任务,Step 2 邮件确认,Step 3 下载结果。

Figure 5.1: GHDC 任务完整流程(用户视角)。图为英文原版,各步骤含义:Step 1 提交任务,Step 2 邮件确认,Step 3 下载结果。


5.1 第一步:准备流域边界文件

GHDC 需要你提供研究流域的矢量边界,格式为 Shapefile(ZIP 压缩包)

5.1.1 边界文件格式要求

  • 压缩包扩展名为 .zip
  • ZIP 内必须包含以下文件(文件名可不同,但扩展名必须完整):
文件 说明
.shp 几何图形
.shx 空间索引
.dbf 属性表
.prj 坐标参考系
  • 几何类型必须为 Polygon(面),不支持线或点
  • 坐标系建议使用地理坐标系(WGS 84 / EPSG:4326)
  • 文件大小上限为 50 MB

提示:如果你没有现成的边界文件,也可以直接在网站地图上手动绘制多边形,系统会自动将绘制结果用作流域边界(见第 5.2 节)。


5.2 第二步:访问网站并提交任务

5.2.1 打开网站

在浏览器中访问 https://ghdc.ac.cn

页面左侧为交互地图,右侧为任务表单(图 5.2)。

GHDC v2 网站主界面。左侧为高德地图(可绘制或上传流域边界),右侧为任务表单。

Figure 5.2: GHDC v2 网站主界面。左侧为高德地图(可绘制或上传流域边界),右侧为任务表单。

5.2.2 划定流域边界

你可以通过以下两种方式之一提供边界:

方式一:上传边界文件

点击地图左上角的「上传边界」按钮,选择已准备好的 .zip 文件。上传成功后,地图会自动缩放显示流域范围,并提示边界面积和中心坐标(图 5.3)。

上传边界文件后,地图自动高亮显示流域轮廓,右上角提示面积和中心坐标。

Figure 5.3: 上传边界文件后,地图自动高亮显示流域轮廓,右上角提示面积和中心坐标。

方式二:在地图上绘制

点击地图左上角的「绘制边界」按钮,在地图上点击若干顶点以勾勒流域轮廓,双击结束绘制。绘制完成后可点击「清除」重新绘制(图 5.4)。

在地图上逐点绘制流域多边形边界,双击关闭多边形。

Figure 5.4: 在地图上逐点绘制流域多边形边界,双击关闭多边形。

5.2.3 填写任务表单

确认边界后,填写右侧表单中的各项参数:

基本信息

字段 说明 示例
邮箱 用于接收确认邮件和结果通知 your@email.com
项目名称 仅限字母、数字、下划线,3–20 字符 lcg_2017

模型参数

字段 说明 本示例值
气象驱动 驱动模型的大气再分析数据集 GLDAS
起始年份 模拟开始年份 2017
结束年份 模拟结束年份 2018
高程数据 DEM 来源 ASTER 30m
最小单元数 流域划分的最少网格单元数 400
最大单元面积 单个网格单元的面积上限(km²) 0.1
含水层厚度 地下含水层厚度(m) 10
土壤数据 土壤质地来源 HWSD
土地利用 土地利用/覆被数据来源 MODIS

各参数的详细含义和选择建议见第 3 章;可用数据集说明见第 2.3 章。

5.2.4 提交任务

确认所有参数无误后,点击「提交任务」按钮(图 5.5)。

填写完成的任务表单。确认参数无误后点击底部「提交任务」按钮。

Figure 5.5: 填写完成的任务表单。确认参数无误后点击底部「提交任务」按钮。

提交成功后,页面会弹出确认对话框,显示任务编号,并提示系统已向你的邮箱发送一封确认邮件(图 5.6)。

提交成功对话框:显示任务编号,并提示确认邮件已发送。

Figure 5.6: 提交成功对话框:显示任务编号,并提示确认邮件已发送。

注意:此时任务尚未正式进入处理队列,需要你完成邮件确认(第三步)后才会开始处理。


5.3 第三步:邮件确认

5.3.1 查收确认邮件

提交后几分钟内,你会收到一封来自 GHDC 的邮件,主题类似:

【GHDC 任务确认】lcg_2017 / [GHDC Task Confirmation] lcg_2017

邮件正文中包含:

  • 任务编号
  • 你所设定的建模参数摘要
  • 一个确认链接

如图 5.7 所示:

收到的确认邮件(示例)。邮件包含建模参数摘要和蓝色确认按钮。

Figure 5.7: 收到的确认邮件(示例)。邮件包含建模参数摘要和蓝色确认按钮。

5.3.2 点击确认链接

点击邮件中的确认链接,浏览器会打开任务状态页面,提示「任务已确认,正在排队处理」(图 5.8)。

确认链接打开的任务状态页面,提示任务已成功确认并进入队列。

Figure 5.8: 确认链接打开的任务状态页面,提示任务已成功确认并进入队列。

确认链接有效期为 24 小时。超时后链接失效,需要重新提交任务。

如果你在 24 小时内没有收到邮件,请检查垃圾邮件文件夹,或联系管理员。


5.4 第四步:等待处理

点击确认后,任务进入处理队列,系统将自动完成以下工作:

  1. 从公开数据源裁剪流域范围内的高程、土壤、土地利用数据
  2. 提取流域水系与子流域边界
  3. 生成非结构三角网格
  4. 准备气象强迫数据
  5. 生成 SHUD 模型所需的全套输入文件
  6. 将所有结果打包为 ZIP 文件

处理时间因流域面积和参数不同而异:

流域面积 预计时间
< 100 km² 15–60 分钟
100–1000 km² 1–4 小时
> 1000 km² 数小时至数天

处理期间你无需保持网页开启,系统完成后会自动发送完成通知邮件


5.5 第五步:下载与使用结果

5.5.1 查收完成邮件

任务处理完成后,你会收到一封完成通知邮件,主题类似:

【GHDC 任务完成】lcg_2017 / [GHDC Task Completed] lcg_2017

邮件中包含结果下载链接,点击即可下载结果 ZIP 文件(图 5.9)。

完成通知邮件(示例)。包含下载按钮和任务信息摘要。

Figure 5.9: 完成通知邮件(示例)。包含下载按钮和任务信息摘要。

5.10 展示了 GHDC 提取的 DEM 与河网叠加效果(以龙漕沟为例):

下载结果中的 DEM 与自动提取的河网(蓝线)叠加图(龙漕沟,16 km²)。颜色表示高程,蓝色线为提取的水系。

Figure 5.10: 下载结果中的 DEM 与自动提取的河网(蓝线)叠加图(龙漕沟,16 km²)。颜色表示高程,蓝色线为提取的水系。

也可以使用命令行工具下载:

wget -c "https://ghdc.ac.cn/download/<你的下载链接>"

5.5.2 结果 ZIP 包含的内容

解压下载的 ZIP 文件后,目录结构如下:

lcg_2017_result/
├── Modeling/
│   ├── input/         # SHUD 模型输入数据(.sp, .ic, .para 等)
│   └── output/        # 空文件夹,供本地运行模型后存放输出
├── StaticFiles/
│   ├── GIS/           # 流域 GIS 数据(边界、DEM、水系、网格等 Shapefile)
│   ├── SHUD_model/    # SHUD 模型可执行程序与分析脚本
│   └── Forcing/       # 气象驱动数据(逐小时 .csv 文件)
└── README.txt         # 任务说明与使用指引

主要文件说明:

路径 内容
Modeling/input/ SHUD 模型可直接读取的输入文件
StaticFiles/GIS/mesh.shp 流域三角网格(可在 GIS 软件中查看)
StaticFiles/GIS/river.shp 提取的河网
StaticFiles/GIS/watershed.shp 子流域边界
StaticFiles/SHUD_model/ 预编译的 SHUD 可执行程序(Linux)
StaticFiles/SHUD_model/rAnalysis/ R 分析脚本(需要 rSHUD 包)

5.5.3 运行 SHUD 模型(可选)

如果你希望在本地完成水文模拟,进入 StaticFiles/SHUD_model/,执行编译与运行:

chmod +x configure
./configure
make clean
make shud
cp shud ../../../Modeling/input/
cd ../../../Modeling/input/
./shud lcg_2017

模型输出结果保存在 Modeling/output/lcg_2017.out/ 中。

5.5.4 使用 rSHUD 分析结果(可选)

StaticFiles/SHUD_model/rAnalysis/ 目录中,用 RStudio 打开 rAnalysis.Rproj,运行示例分析脚本,可绘制出口径流过程线、地下水空间分布等图表:

library(rSHUD)

prjname <- 'lcg_2017'
shud.env(prjname = prjname,
         inpath  = file.path('../../../Modeling/input/', prjname),
         outpath = file.path('../../../Modeling/output/', paste0(prjname, '.out')))

# 出口径流过程线
oid  <- getOutlets()
qdown <- readout('rivqdown')
qs   <- ts2Daily(qdown[, oid[1]])
autoplot(qs)

# 地下水空间分布
gw    <- readout('eleygw')
mgw   <- colMeans(gw)
sp.mesh <- sp.mesh2Shape(dbf = data.frame(gw = mgw))
plot_sp(sp.mesh, zcol = 'gw')

关于 rSHUD 的详细使用方法,请参考 rSHUD 文档


小结:一次完整的 GHDC 建模任务只需五步:上传边界 → 填写参数 → 确认邮件 → 等待处理 → 下载结果。接下来的章节将详细介绍各参数的含义(第 3 章)与可用数据集(第 2.3 章)。