从 EFAS 局部排水方向划定盆地

Question

我需要使用欧洲洪水意识系统 (https://confluence.ecmwf.int/display/CEMS/Auxiliary+Data) 提供的 LDD 来划定出水口部分的分水岭。我想知道如何在 QGIS 或 R 中执行此操作。

Answer 1

一般来说，你会

efas_100 <- rast('~/Downloads/flood_threshold_efas_v5.1_rl_100.0.nc')

efas_100

class       : SpatRaster 
dimensions  : 2970, 4530, 1  (nrow, ncol, nlyr)
resolution  : 0.01666667, 0.01666667  (x, y)
extent      : -25.25, 50.25, 22.75, 72.25  (xmin, xmax, ymin, ymax)
coord. ref. : lon/lat WGS 84 (CRS84) (OGC:CRS84) 
source      : flood_threshold_efas_v5.1_rl_100.0.nc 
varname     : rl_100.0 
name        : rl_100.0

efas_flow <- terrain(efas_100, 'flowdir')
# create pout point of interest as matrix or vect
shed <- watershed(efas_flow, pp) 

然后提供一个倾点，您有兴趣描绘的地方。您需要裁剪

efas_100

或

efas_flow

，因为分水岭比横贯大陆的 efas_100 小得多。