GEOtop_CreateInptsMAPS: Create GEOtop input maps

In JBrenn/AnalyseGEOtop: Analyse GEOtop simulations

Description

Create input maps for hydrological modeling (e.g. GEOtop, mHM) using RSAGA. Base maps are a digital elevation model and a map locating the basin outlet (defined by -1). Postprocessing includes terrain analysis (slope, aspect, curvature), flow direction and accumulation, as well as sky view factor calculation. All maps, including landcover, soil type, leaf area index (LAI) and bedrock (if provided), are stored in SAGA .srgd format. Further processing includes watershed delineation, resampling to required resolution, masking to the extent of the watershed (either for all maps or for landcover map only) and writing of GEOtop or mHM input ascii files.

Usage

GEOtop_CreateInptsMAPS(dem, res, stream_gauge, name_gauge, 
preprocess=TRUE, 
SkyView=FALSE, FlowAccum=FALSE, rivernet=FALSE, idgauges=FALSE,
soil=FALSE, landcover=FALSE, bedrock=FALSE, LAI=FALSE, geology=FALSE,
SAGAENV, mask_only_lc=FALSE, GEOtop=FALSE, MHM=FALSE, ASCII=TRUE)

Arguments

dem	ascii digital elevation model, located in working directory
res	resolution of the output maps in meters
stream_gauge	ascii, defining location of stream gauge, located in working directory (NA ascii, stream gauge location = -1)
minArea	integer, minimum size of basin (number of cells)
name_gauge	name of the stream gauge
preprocess	boolean, TRUE: preprocessing (terrain analysis); FALSE: no preprocessing, use maps created in former run; default = TRUE
SkyView	boolean, TRUE: create SkyViewFactor .asc map; default = FALSE
FlowAccum	boolean, TRUE: create flow accumulation (RHO8) .asc map; default = FALSE
rivernet	boolean, TRUE: create river network .asc map; default = FALSE
idgauges	boolean, TRUE: map of gauges IDs, describing gauge locations provided in working directory, ascii format, name: idgauges.asc; default = FALSE
soil	boolean, TRUE: soil type map provided in working directory, ascii format, name: soil.asc; FALSE: no soil type map provided; default = FALSE
landcover	boolean, TRUE: landcover map provided in working directory, ascii format, name: landcover.asc; FALSE: no landcover map provided; default = FALSE
lc_scenes	chr vector, provide names of .asc maps containing landcover scences (function will loop over them), e.g. c("lc_1990.asc", "lc_2000.asc", "lc_2006.asc")
bedrock	boolean, TRUE: bedrock map provided in working directory, ascii format, name: bedrock.asc; FALSE: no bedrock map provided; default = FALSE
LAI	boolean, TRUE: LAI map provided in working directory, ascii format, name: lai.asc; FALSE: no LAI map provided; default = FALSE
geology	boolean, TRUE: geology class map provided in working directory, ascii format, name: geology.asc; FALSE: no LAI map provided; default = FALSE
SAGAENV	SAGA environment, created with function rsaga.env of the RSAGA package
mask_only_lc	boolean, TRUE: mask only landcover map; FALSE: mask all maps and crop to data; default = FALSE
GEOtop	boolean, TRUE: create GEOtop input maps in ascii format; default = FALSE
MHM	boolean, TRUE: create mHM input maps in ascii format; default = FALSE
ASCII	boolean, TRUE: create maps in ascii format; default = TRUE

Value

GEOtop input maps, storged in working directory as ascii maps

Author(s)

Johannes Brenner, [e-mail]

JBrenn/AnalyseGEOtop documentation built on May 7, 2019, 6:48 a.m.