forked from JBrenn/AnalyseGEOtop
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy path.Rhistory
512 lines (512 loc) · 33.4 KB
/
.Rhistory
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
krig_ord_global <- OrdKrig(wpath = "/home/jbre/Schreibtisch/krig_SMC/", datafolder = "master", rastermask = NA,
inverseDistWeigths = F, local = F, method = "ordinary",
variable = "SoilMoisture_mean", npix = 1, cutoff = 20, anis_deg = 0, anis_ax = .5,
psill = 1, nugget = 1, nmax = 20, nmin = 2, omax = 0)
plot(krig_ord_global$P2$map_pred)
krig_uni_global <- OrdKrig(wpath = "/home/jbre/Schreibtisch/krig_SMC/", datafolder = "master", rastermask = NA,
inverseDistWeigths = F, local = F, method = "universal", add_vars = "slope",
variable = "SoilMoisture_mean", npix = 1, cutoff = 20, anis_deg = 0, anis_ax = .5,
psill = 1, nugget = 1, nmax = 20, nmin = 2, omax = 0)
method
filesIN <- dir(file.path(wpath, datafolder))
if (is.null(names(cutoff)))
{
# get namezones according to input files
namezones <- c()
for (namefile in filesIN)
{
nr_char <- nchar(filesIN)
names(nr_char) <- filesIN
namezone <- strsplit(substr(namefile,1,nr_char[namefile]-4), "_")[[1]][2]
namezones <- c(namezones, namezone)
}
# list of arguments
args <- list("cutoff"=cutoff, "anis_deg"=anis_deg,"anis_ax"=anis_ax,"psill"=psill,
"nugget"=nugget,"nmax"=nmax,"nmin"=nmin,"idp"=idp)
# assing argument names to parameter value (named vector)
for (arg in 1:length(args))
{
x <- args[[arg]]
if (length(x)==length(namezones)) {
names(x) <- namezones
assign(x = names(args)[arg], value = x)
} else {
print(paste("argument", names(args)[arg], "is not given with name of investigated zone and differs in length from input .csv files in the folder", datafolder, ". Please check!", sep=" " ))
}
}
}
val_list <- list()
namefile
nr_char <- nchar(filesIN)
names(nr_char) <- filesIN
namezone <- strsplit(substr(namefile,1,nr_char[namefile]-4), "_")[[1]][2]
print(paste("processing zone", namezone, sep=" "))
# read table
worktab <- read.table(file = file.path(wpath, datafolder, namefile), header = TRUE, sep = ",", dec = ".")
if (method=="ordinary") {
var_names <- "VARIABLE"; choose <- variable
}
if (method=="universal") {
var_names <- c("VARIABLE")
for (i in 1:length(add_vars)) var_names <- c(var_names, paste("add_var",i, sep="_"))
choose <- c(variable, add_vars)
}
worktab <- cbind(worktab$x_Coord,worktab$y_Coord,worktab[,choose])
# matrix 2 data.frame
worktab <- as.data.frame(worktab)
# rename cols
names(worktab) <- c("X","Y", var_names)
# zeros
worktab[worktab$VARIABLE <= 0,"VARIABLE"] <- 0.001
str(worktab)
var_names
?krige
source('~/GitHub/SpatialInterpol/R/OrdKrig.R', echo=TRUE)
krig_uni_global <- OrdKrig(wpath = "/home/jbre/Schreibtisch/krig_SMC/", datafolder = "master", rastermask = NA,
inverseDistWeigths = F, local = F, method = "universal", add_vars = "slope",
variable = "SoilMoisture_mean", npix = 1, cutoff = 20, anis_deg = 0, anis_ax = .5,
psill = 1, nugget = 1, nmax = 20, nmin = 2, omax = 0)
source('~/GitHub/SpatialInterpol/R/OrdKrig.R', echo=TRUE)
krig_uni_global <- OrdKrig(wpath = "/home/jbre/Schreibtisch/krig_SMC/", datafolder = "master", rastermask = NA,
inverseDistWeigths = F, local = F, method = "universal", add_vars = "slope",
variable = "SoilMoisture_mean", npix = 1, cutoff = 20, anis_deg = 0, anis_ax = .5,
psill = 1, nugget = 1, nmax = 20, nmin = 2, omax = 0)
source('~/GitHub/SpatialInterpol/R/OrdKrig.R', echo=TRUE)
krig_uni_global <- OrdKrig(wpath = "/home/jbre/Schreibtisch/krig_SMC/", datafolder = "master", rastermask = NA,
inverseDistWeigths = F, local = F, method = "universal", add_vars = "slope",
variable = "SoilMoisture_mean", npix = 1, cutoff = 20, anis_deg = 0, anis_ax = .5,
psill = 1, nugget = 1, nmax = 20, nmin = 2, omax = 0)
source('~/GitHub/SpatialInterpol/R/OrdKrig.R', echo=TRUE)
krig_uni_global <- OrdKrig(wpath = "/home/jbre/Schreibtisch/krig_SMC/", datafolder = "master", rastermask = NA,
inverseDistWeigths = F, local = F, method = "universal", add_vars = "slope",
variable = "SoilMoisture_mean", npix = 1, cutoff = 20, anis_deg = 0, anis_ax = .5,
psill = 1, nugget = 1, nmax = 20, nmin = 2, omax = 0)
source('~/GitHub/SpatialInterpol/R/OrdKrig.R', echo=TRUE)
krig_uni_global <- OrdKrig(wpath = "/home/jbre/Schreibtisch/krig_SMC/", datafolder = "master", rastermask = NA,
inverseDistWeigths = F, local = F, method = "universal", add_vars = "slope",
variable = "SoilMoisture_mean", npix = 1, cutoff = 20, anis_deg = 0, anis_ax = .5,
psill = 1, nugget = 1, nmax = 20, nmin = 2, omax = 0)
source('~/GitHub/SpatialInterpol/R/OrdKrig.R', echo=TRUE)
krig_uni_global <- OrdKrig(wpath = "/home/jbre/Schreibtisch/krig_SMC/", datafolder = "master", rastermask = NA,
inverseDistWeigths = F, local = F, method = "universal", add_vars = "slope",
variable = "SoilMoisture_mean", npix = 1, cutoff = 20, anis_deg = 0, anis_ax = .5,
psill = 1, nugget = 1, nmax = 20, nmin = 2, omax = 0)
str(worktab)
coordinates(worktab) <- ~X+Y
crs(worktab) <- coordsys
coordinates(worktab) <- ~X+Y
str(worktab)
str( worktab@data)
library()
library(soilwater)
?swc
library(DataBaseAlpEnvEURAC)
?dB_readStationData
path <- "/media/alpenv/Projekte/MONALISA/05_Arbeitsbereiche/BrJ/01_data/Stations_MONALISA/NEMEF/NEMEF1500/"
header.file = "/media/alpenv/Projekte/MONALISA/05_Arbeitsbereiche/BrJ/01_data/Stations_MONALISA/NEMEF/header_NEMEF1500.txt"
station = "NEMEF1500"
files <- dir(path)
files <- dir(path)
files
header <- as.character(read.table(header.file, header=FALSE, sep=",")[,1])
header_org <- header
header
station_gen <- substr(station,1,nchar(station)-4)
station_gen
if (station_gen=="BERAT") {
skip <- 1; date_col=2; tz="Etc/GMT+1"; time_step=times("00:05:00")
} else time_step <- times("00:15:00")
if (station_gen=="SF") {
skip <- 1; date_col=2; tz="Etc/GMT+1"
}
if (station_gen=="B") {
skip <- 4; date_col=1; tz="Etc/GMT-2"
header_final <- paste(substr(header.file, 1, nchar(header.file)-16), "header_final.txt", sep="")
header_final <- as.character(read.table(header_final, header=FALSE)[,1])
}
if (station_gen=="P"| station_gen=="I") {
skip <- 4; date_col=1; tz="Etc/GMT-1"
}
if (station_gen=="M" | station_gen=="S") {
skip <- 1; date_col=2; tz="Etc/GMT+1"
header_final <- paste(substr(header.file, 1, nchar(header.file)-16), "header_final.txt", sep="")
header_final <- as.character(read.table(header_final, header=FALSE)[,1])
}
if (station_gen=="XS") {
skip <- 1; date_col=2; tz="Etc/GMT+1"
header_final <- paste(substr(header.file, 1, nchar(header.file)-17), "header_final.txt", sep="")
header_final <- as.character(read.table(header_final, header=FALSE)[,1])
}
if (station=="S2") {
skip <- 1; date_col=2; tz="Etc/GMT+2"
}
if (station_gen=="DOMEF" | station_gen=="DOMES" | station_gen=="DOPAS" | station_gen=="NEMEF") {
skip <- 5; date_col=1; tz="Etc/GMT+1"
}
skip
data <- rep(NA,length(header))
datetime <- chron(dates. = "2012-02-02", times. = "01:00:00",
format= c(dates = "y-m-d", times = "h:m:s"))
library(chron)
library(chron)
library(zoo)
data <- rep(NA,length(header))
datetime <- chron(dates. = "2012-02-02", times. = "01:00:00",
format= c(dates = "y-m-d", times = "h:m:s"))
source('~/GitHub/DataBaseAlpEnvEURAC/R/dB_readStationData.R', echo=TRUE)
data <- dB_readStationData(path = path, header.file = header.file, station = station)
str(data)
source('~/GitHub/DataBaseAlpEnvEURAC/R/dB_readStationData.R', echo=TRUE)
data <- dB_readStationData(path = path, header.file = header.file, station = station)
str(data)
smc <- dB_getSWC(path2data = path, station = "NEMEF1500", aggregation = "h", minVALUE = 0, maxVALUE = 1)
path
path2data <-"/media/alpenv/Projekte/MONALISA/05_Arbeitsbereiche/BrJ/01_data/Stations_MONALISA/NEMEF"
smc <- dB_getSWC(path2data = path2data, station = "NEMEF1500", aggregation = "h", minVALUE = 0, maxVALUE = 1)
path2data <-"/media/alpenv/Projekte/MONALISA/05_Arbeitsbereiche/BrJ/01_data/Stations_MONALISA"
smc <- dB_getSWC(path2data = path2data, station = "NEMEF1500", aggregation = "h", minVALUE = 0, maxVALUE = 1)
station
station_nr <- as.integer(substr(station, nchar(station)-3, nchar(station)))
station_nr
station_ <- substr(station, 1, nchar(station)-4)
station_
if (station_ == "XS") station_ <- "S"
path2files = file.path(path2data,station_,station)
header.file = file.path(path2data,station_,paste("header_",station,".txt",sep=""))
path2files
header.file
data_raw <- dB_readStationData(path = path2files, header.file = header.file, station = station)
smc <- dB_getSWC(path2data = path2data, station = "NEMEF1500", aggregation = "h", minVALUE = 0, maxVALUE = 1)
plot(smc)
source('~/GitHub/DataBaseAlpEnvEURAC/R/dB_getSWC.R', echo=TRUE)
smc <- dB_getSWC(path2data = path2data, station = "NEMEF1500", aggregation = "h", minVALUE = 0, maxVALUE = 1)
plot(smc)
smc_d <- dB_getSWC(path2data = path2data, station = "NEMEF1500", aggregation = "d", minVALUE = 0, maxVALUE = 1)
plot(smc)
plot(smc_d)
nemef <- dB_updatedb(stations = "NEMEF1500", variables = "TOTAL", path2data = path2data, inCloud = "/home/jbre/Schreibtisch/", write_csv = T, return_data = T)
plot(nemef)[,1:3]
plot(nemef)[,1]
plot(nemef[,5])
plot(nemef[5])
str(nemef)
str(nemef$NEMEF1500[,5])
plot(nemef$NEMEF1500[,5])
plot(nemef$NEMEF1500[,5:10])
nemef <- dB_updatedb(stations = "NEMEF1500", variables = "SWC", path2data = path2data, inCloud = "/home/jbre/Schreibtisch/", write_csv = T, return_data = T)
?dB_updatedb
nemef <- dB_updatedb(stations = "NEMEF1500", variables = "SMC", path2data = path2data, inCloud = "/home/jbre/Schreibtisch/", write_csv = T, return_data = T)
plot(nemef)
plot(nemef$NEMEF1500)
nemef <- dB_updatedb(stations = c("DOPAS2000","NEMEF1500"), variables = "TOTAL", path2data = path2data, inCloud = "/home/jbre/Schreibtisch/", write_csv = T, return_data = T)
nemef <- dB_updatedb(stations = c("DOPAS2000","NEMEF1500"), variables = "TOTAL", path2data = path2data, inCloud = "/home/jbre/Schreibtisch/", write_csv = T, return_data = T)
nemef <- dB_updatedb(stations = c("DOPAS2000","NEMEF1500"), variables = "TOTAL", path2data = path2data, inCloud = "/home/jbre/Schreibtisch/", write_csv = T, return_data = T)
nemef <- dB_updatedb(stations = c("DOMEF1500"), variables = "TOTAL", path2data = path2data, inCloud = "/home/jbre/Schreibtisch/", write_csv = T, return_data = T)
nemef <- dB_updatedb(stations = c("DOMEF1500"), variables = "SWC", path2data = path2data, inCloud = "/home/jbre/Schreibtisch/", write_csv = T, return_data = T)
?dB_getSWC
source('~/GitHub/DataBaseAlpEnvEURAC/R/dB_updatedb.R', echo=TRUE)
nemef <- dB_updatedb(stations = c("DOMEF1500"), variables = "SWC", path2data = path2data, inCloud = "/home/jbre/Schreibtisch/", write_csv = T, return_data = T)
library(Helper4me)
Helper4me::help_updateGITlib("DataBaseAlpEnvEURAC")
library(DataBaseAlpEnvEURAC)
nemef <- dB_updatedb(stations = c("DOMEF1500"), variables = "SWC", path2data = "/media/alpenv/Projekte/MONALISA/05_Arbeitsbereiche/BrJ/01_data/Stations_MONALISA/", inCloud = "/home/jbre/Schreibtisch/", write_csv = T, return_data = T)
plot(nemef)
plot(nemef$DOMEF1500)
require(RSAGA)
rsaga.get.libraries()
work_env <- rsaga.env(path = '/usr/bin',modules = '/usr/lib64/saga')
work_env <- rsaga.env(path = '/usr/bin',modules = '/usr/lib/saga')
work_env <- rsaga.env()
work_env <- rsaga.env(path = '/usr/bin')
work_env <- rsaga.env(path = '/usr/bin', modules = '/usr/lib/x86_64-linux-gnu/saga
')
work_env <- rsaga.env(path = '/usr/bin', modules = '/usr/lib/x86_64-linux-gnu/saga
')
work_env <- rsaga.env(path = '/usr/bin', modules = '/usr/lib/x86_64-linux-gnu/saga/
')
work_env <- rsaga.env(path = "/usr/bin", modules = "/usr/lib/x86_64-linux-gnu/saga/
")
work_env <- rsaga.env(modules = "/usr/lib/x86_64-linux-gnu/saga/
")
work_env <- rsaga.env(modules = "/usr/lib/x86_64-linux-gnu/saga/")
work_env <- rsaga.env(path = "/usr/bin", modules = "/usr/lib/x86_64-linux-gnu/saga/")
rsaga.get.libraries()
rsaga.get.libraries(path = work_env)
rsaga.get.libraries(env=work_env)
rsaga.get.libraries()
rsaga.get.libraries(path = "/usr/lib/x86_64-linux-gnu/saga/" )
rsaga.get.modules(lib = "ta_lighting, path = "/usr/lib/x86_64-linux-gnu/saga/" )
rsaga.get.modules(lib = "ta_lighting", path = "/usr/lib/x86_64-linux-gnu/saga/")
rsaga.get.modules(lib = "ta_lighting", env = work_env)
rsaga.get.modules(lib = "ta_hydrology", env = work_env)
?rsaga.geoprocessor
rsaga.get.modules(lib = "ta_preprocessor", env = work_env)
library(raster)
faccum <- raster("/home/jbre/Schreibtisch/Jucar_Basin_12 2.2011/GIS/flow_accum_rho8.tif")
str(values(faccum))
max(values(faccum) na.rm=T)
max(values(faccum), na.rm=T)
which.max(values(faccum), na.rm=T)
which.max(values(faccum))
values(faccum)[6271077]
which.max(values(faccum)==-1)
values(faccum)[319448]
which.max(values(faccum)==(-1))
values(faccum)[319448]
which.max(values(faccum)==(1))
which.max(values(faccum))
faccum <- NA
summary(values(faccum))
all(is.na(values(faccum)))
faccum <- raster("/home/jbre/Schreibtisch/Jucar_Basin_12 2.2011/GIS/flow_accum_rho8.tif")
values(faccum) <- NA
all(is.na(values(faccum)))
values(faccum)[6271077] <- -1
all(is.na(values(faccum)))
writeRaster(x = faccum, filename = "facc.asc")
faccum <- raster("/home/jbre/Schreibtisch/Jucar_Basin_12 2.2011/GIS/dem_filled_100.tif")
writeRaster(x = faccum, filename = "dem.asc")
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
setwd("/home/jbre/Schreibtisch/Jucar_Basin_12 2.2011/SAGA")
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
SAGAENV <- rsaga.env(path = "/usr/bin",
modules = "/usr/lib/x86_64-linux-gnu/saga/")
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "facc.asc", name_gauge = "Jucar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, SAGAENV = SAGAENV, mask_only_lc = F, GEOTOP = F, MHM = T)
SAGAENV
SAGAENV <- rsaga.env(path = "/usr/bin",
modules = "/usr/lib/x86_64-linux-gnu/saga/", version = "2.2.3")
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "facc.asc", name_gauge = "Jucar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, SAGAENV = SAGAENV, mask_only_lc = F, GEOTOP = F, MHM = T)
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "facc.asc", name_gauge = "Jucar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, SAGAENV = SAGAENV, mask_only_lc = F, GEOTOP = F, MHM = T)
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "facc.asc", name_gauge = "Jucar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, SAGAENV = SAGAENV, mask_only_lc = F, GEOTOP = F, MHM = T)
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "facc.asc", name_gauge = "Jucar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, SAGAENV = SAGAENV, mask_only_lc = F, GEOTOP = F, MHM = T)
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "facc.asc", name_gauge = "Jucar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, SAGAENV = SAGAENV, mask_only_lc = F, GEOTOP = F, MHM = T)
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "facc.asc", name_gauge = "Jucar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, SAGAENV = SAGAENV, mask_only_lc = F, GEOTOP = F, MHM = T)
SAGAENV
as.integer(SAGAENV$version)
as.integer(strsplit(SAGAENV$version,"."))
(strsplit(SAGAENV$version,"."))
(strsplit(".",SAGAENV$version))
SAGAENV$version
SAGAENV$version,1,2
substr(SAGAENV$version,1,2)
substr(SAGAENV$version,1,1)
int(substr(SAGAENV$version,1,1))
integer(substr(SAGAENV$version,1,1))
as.integer(substr(SAGAENV$version,1,1))
int.saga.version <- as.integer(substr(SAGAENV$version,1,1)) *100 + as.integer(substr(SAGAENV$version,3,3)) *10 + as.integer(substr(SAGAENV$version,5,5))
int.saga.version
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "facc.asc", name_gauge = "Jucar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, SAGAENV = SAGAENV, mask_only_lc = F, GEOTOP = F, MHM = T)
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "facc.asc", name_gauge = "Jucar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, SAGAENV = SAGAENV, mask_only_lc = F, GEOTOP = F, MHM = T)
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "facc.asc", name_gauge = "Jucar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, SAGAENV = SAGAENV, mask_only_lc = F, GEOTOP = F, MHM = T)
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "facc.asc", name_gauge = "Jucar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, SAGAENV = SAGAENV, mask_only_lc = F, GEOtop = F, MHM = T)
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "facc.asc", name_gauge = "Jucar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, SAGAENV = SAGAENV, mask_only_lc = F, GEOtop = F, MHM = T)
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "facc.asc", name_gauge = "Jucar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, SAGAENV = SAGAENV, mask_only_lc = F, GEOtop = F, MHM = T)
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "gauge_Jucar.asc", name_gauge = "Jucar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, SAGAENV = SAGAENV, mask_only_lc = F, GEOtop = F, MHM = T
efws
facc <- raster("/home/jbre/Schreibtisch/Jucar_Basin_12 2.2011/GIS/flow_accum_rho8.tif")
which(values(facc)==2043809)
which(values(facc)==2043870)
values(facc) <- NA
values(which(values(facc)==2043870)) <- -1
values(facc)[which(values(facc)==2043870)] <- -1
summary(values(facc))
unique(values(facc))
which(values(facc)==2043870)
values(facc)[6002102] <- -1
writeRaster(x = facc, filename = "gaugeJucar.asc")
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "gaugeJucar.asc", name_gauge = "Jucar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, SAGAENV = SAGAENV, mask_only_lc = F, GEOtop = F, MHM = T
)
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "gaugeJucar.asc", minArea=1000, name_gauge = "Jucar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, SAGAENV = SAGAENV, mask_only_lc = F, GEOtop = F, MHM = T)
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "gaugeJucar1.asc", minArea=100, name_gauge = "Jucar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, SAGAENV = SAGAENV, mask_only_lc = F, GEOtop = F, MHM = T)
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "gaugeJucar1.asc", minArea=100, name_gauge = "Jucar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, SAGAENV = SAGAENV, mask_only_lc = F, GEOtop = F, MHM = T)
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "gaugeJucar1.asc", minArea=100, name_gauge = "Jucar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, SAGAENV = SAGAENV, mask_only_lc = F, GEOtop = F, MHM = T)
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "gaugeJucar1.asc", minArea=100, name_gauge = "Jucar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, SAGAENV = SAGAENV, mask_only_lc = F, GEOtop = F, MHM = T)
?rsaga.sgrd.to.esri
ascii <- readLines(con =paste("./DEM/dem",res,".asc",sep=""))
res
res=100
ascii <- readLines(con =paste("./DEM/dem",res,".asc",sep=""))
header <- ascii[1:6]
str(ascii)
?vapply
sapply(acsii, function(x) gsub("-99999.00000", "-9999", x))
ascii <- readLines(con =paste("./DEM/dem",res,".asc",sep=""))
sapply(acsii, function(x) gsub("-99999.00000", "-9999", x))
sapply(ascii, function(x) gsub("-99999.00000", "-9999", x))
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
files <- dir("./mHMASC100/", full.names = T)
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "gaugeJucar1.asc", minArea=100, name_gauge = "Jucar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, SAGAENV = SAGAENV, mask_only_lc = F, GEOtop = F, MHM = T)
library(RSAGA)
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "gaugeJucar1.asc", minArea=100, name_gauge = "Jucar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, SAGAENV = SAGAENV, mask_only_lc = F, GEOtop = F, MHM = T)
SAGAENV <- rsaga.env(path = "/usr/bin",
modules = "/usr/lib/x86_64-linux-gnu/saga/", version = "2.2.3")
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "gaugeJucar1.asc", minArea=100, name_gauge = "Jucar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, SAGAENV = SAGAENV, mask_only_lc = F, GEOtop = F, MHM = T)
setwd("/home/jbre/Schreibtisch/Jucar_Basin_12 2.2011/SAGA")
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "gaugeJucar1.asc", minArea=100, name_gauge = "Jucar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, SAGAENV = SAGAENV, mask_only_lc = F, GEOtop = F, MHM = T)
files <- dir("./mHMASC100/", full.names = T)
files
for (i in files)
{
print(paste("create",i))
system(paste("sed -i -e 's/-99999.00000/-9999/g' ",i,sep=""))
}
files <- dir("./mHMASC100", full.names = T)
for (i in files)
{
print(paste("create",i))
system(paste("sed -i -e 's/-99999.00000/-9999/g' ",i,sep=""))
}
?rsaga.sgrd.to.esri
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "gaugeJucar1.asc", minArea=100, name_gauge = "Jucar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, SAGAENV = SAGAENV, mask_only_lc = F, GEOtop = F, MHM = T, rivernet = F)
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "gaugeJucar1.asc", minArea=100, name_gauge = "Jucar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, SAGAENV = SAGAENV, mask_only_lc = F, GEOtop = F, MHM = T, rivernet = F)
fdir <- read.ascii.grid(file = paste("./FlowDir/fdir",res,".asc",sep=""))
res=100
fdir <- read.ascii.grid(file = paste("./FlowDir/fdir",res,".asc",sep=""))
str(fdir)
fdir$data[fdir$data == 1] <- 128
fdir$data[fdir$data == 2] <- 1
fdir$data[fdir$data == 3] <- 2
fdir$data[fdir$data == 5] <- 8
fdir$data[fdir$data == 6] <- 16
fdir$data[fdir$data == 7] <- 32
unique(fdir$data)
summary(fdir$data)
min(fdir$data)
min(fdir$data, na.rm=T)
fdir$data[fdir$data == 0] <- 64
fdir$data[fdir$data == 0] <- 64
min(fdir$data, na.rm=T)
write.ascii.grid(data = asp, file = paste("./mHMASC",res,"/fdir.asc",sep=""))
write.ascii.grid(data = fdir, file = paste("./mHMASC",res,"/fdir.asc",sep=""))
function (..., sep = " ", collapse = NULL)
ascii <- readLines(con = paste("./mHMASC",res,"/fdir.asc",sep=""))
ascii <- readLines(con = paste("./mHMASC",res,"/fdir.asc",sep=""))
header <- ascii[1:6]
setwd("Xuquer/")
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "xuquer_gauge.asc", minArea=100, name_gauge = "Jucar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, SAGAENV = SAGAENV, mask_only_lc = F, GEOtop = F, MHM = T, rivernet = F)
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "xuquer_gauge.asc", minArea=100, name_gauge = "Jucar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, SAGAENV = SAGAENV, mask_only_lc = F, GEOtop = F, MHM = T, rivernet = F)
View(GEOtop_CreateInptsMAPS)
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
View(GEOtop_CreateInptsMAPS)
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "xuquer_gauge.asc", minArea=100, name_gauge = "Jucar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, SAGAENV = SAGAENV, mask_only_lc = F, GEOtop = F, MHM = T, rivernet = F)
facc <- raster("/home/jbre/Schreibtisch/Jucar_Basin_12 2.2011/GIS/flow_accum_rho8.tif")
library()
library(raster)
facc <- raster("/home/jbre/Schreibtisch/Jucar_Basin_12 2.2011/GIS/flow_accum_rho8.tif")
values(facc)==625095
which(values(facc)==625095)
which(values(facc)==401346)
values(facc) <- NA
values(facc)[5112563] <- -1
writeRaster(facc, "guadalaviar_gauge.asc")
values(facc) <- NA
values(facc)[3368830] <- -1
writeRaster(facc, "millars_gauge.asc")
setwd("../Guadalaviar/")
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "guadalaviar_gauge.asc", minArea=100, name_gauge = "Jucar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, SAGAENV = SAGAENV, mask_only_lc = F, GEOtop = F, MHM = T, rivernet = F)
setwd("../Millars/")
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "millars_gauge.asc", minArea=100, name_gauge = "Jucar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, SAGAENV = SAGAENV, mask_only_lc = F, GEOtop = F, MHM = T, rivernet = F)
facc <- raster("/home/jbre/Schreibtisch/Jucar_Basin_12 2.2011/GIS/flow_accum_rho8.tif")
which(values(facc)==97381)
values(facc) <- NA
values(facc)[4245470] <- -1
writeRaster(facc, "../Palancia/palancia_gauge.asc")
setwd("../Palancia")
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "palancia_gauge.asc", minArea=100, name_gauge = "Jucar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, SAGAENV = SAGAENV, mask_only_lc = F, GEOtop = F, MHM = T, rivernet = F)
facc <- raster("/home/jbre/Schreibtisch/Jucar_Basin_12 2.2011/GIS/flow_accum_rho8.tif")
which(values(facc)==158259)
values(facc) <- NA
values(facc)[9621158] <- -1
writeRaster(facc, "../Segura/segura_gauge.asc")
setwd("../Segura/")
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "segura_gauge.asc", minArea=100, name_gauge = "Jucar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, SAGAENV = SAGAENV, mask_only_lc = F, GEOtop = F, MHM = T, rivernet = F)
facc <- raster("/home/jbre/Schreibtisch/Jucar_Basin_12 2.2011/GIS/flow_accum_rho8.tif")
which(values(facc)==50861)
values(facc) <- NA
values(facc)[8699751] <- -1
writeRaster(facc, "../Campello/campello_gauge.asc")
setwd("../Campello/")
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "campello_gauge.asc", minArea=100, name_gauge = "Jucar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, SAGAENV = SAGAENV, mask_only_lc = F, GEOtop = F, MHM = T, rivernet = F)
facc <- raster("/home/jbre/Schreibtisch/Jucar_Basin_12 2.2011/GIS/flow_accum_rho8.tif")
which(values(facc)==74225)
which(values(facc)==74210)
which(values(facc)==74200)
values(facc) <- NA
values(facc)[6651939] <- -1
writeRaster(facc, "../Serpis/serpis_gauge.asc")
setwd("../Serpis/")
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "serpis_gauge.asc", minArea=100, name_gauge = "Jucar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, SAGAENV = SAGAENV, mask_only_lc = F, GEOtop = F, MHM = T, rivernet = F)
facc <- raster("/home/jbre/Schreibtisch/Jucar_Basin_12 2.2011/GIS/flow_accum_rho8.tif")
which(values(facc)==21302)
values(facc) <- NA
values(facc)[8005651] <- -1
writeRaster(facc, "../Algar/algar_gauge.asc")
setwd("../Algar/")
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "algar_gauge.asc", minArea=100, name_gauge = "Jucar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, SAGAENV = SAGAENV, mask_only_lc = F, GEOtop = F, MHM = T, rivernet = F)
facc <- raster("/home/jbre/Schreibtisch/Jucar_Basin_12 2.2011/GIS/flow_accum_rho8.tif")
which(values(facc)==21796)
values(facc) <- NA
values(facc)[8386287] <- -1
writeRaster(facc, "../Amadorio/amadorio_gauge.asc")
setwd("../Amadorio/")
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "amadorio_gauge.asc", minArea=100, name_gauge = "Jucar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, SAGAENV = SAGAENV, mask_only_lc = F, GEOtop = F, MHM = T, rivernet = F)
facc <- raster("/home/jbre/Schreibtisch/Jucar_Basin_12 2.2011/GIS/flow_accum_rho8.tif")
which(values(facc)==34699)
values(facc) <- NA
values(facc)[1311628] <- -1
writeRaster(facc, "../Vinaros/vinaros_gauge.asc")
setwd("../Vinaros/")
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "vinaros_gauge.asc", minArea=100, name_gauge = "Jucar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, SAGAENV = SAGAENV, mask_only_lc = F, GEOtop = F, MHM = T, rivernet = F)
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "vinaros_gauge.asc", minArea=100, name_gauge = "Jucar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, SAGAENV = SAGAENV, mask_only_lc = F, GEOtop = F, MHM = T, rivernet = F)
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
setwd("/home/jbre/Schreibtisch/Jucar_Basin_12 2.2011/SAGA/Algar/")
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "algar_gauge.asc", minArea=100, name_gauge = "Algar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, SAGAENV = SAGAENV, mask_only_lc = F, GEOtop = F, MHM = T, rivernet = F)
library(RSAGA)
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "algar_gauge.asc", minArea=100, name_gauge = "Algar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, SAGAENV = SAGAENV, mask_only_lc = F, GEOtop = F, MHM = T, rivernet = F)
SAGAENV <- rsaga.env(path = "/usr/bin",
modules = "/usr/lib/x86_64-linux-gnu/saga/", version = "2.2.3")
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "algar_gauge.asc", minArea=100, name_gauge = "Algar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, mask_only_lc = F, GEOtop = F, MHM = T, rivernet = F)
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "algar_gauge.asc", minArea=100, name_gauge = "Algar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, mask_only_lc = F, GEOtop = F, MHM = T, rivernet = F)
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "algar_gauge.asc", minArea=100, name_gauge = "Algar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, mask_only_lc = F, GEOtop = F, MHM = T, rivernet = F)
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "algar_gauge.asc", minArea=100, name_gauge = "Algar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, mask_only_lc = F, GEOtop = F, MHM = T, rivernet = F)
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "algar_gauge.asc", minArea=100, name_gauge = "Algar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, mask_only_lc = F, GEOtop = F, MHM = T, rivernet = F)
source('~/GitHub/AnalyseGEOtop/R/GEOtop_CreateInptsMAPS.R', echo=TRUE)
GEOtop_CreateInptsMAPS(dem = "dem.asc", res = 100, stream_gauge = "algar_gauge.asc", minArea=100, name_gauge = "Algar", preprocess = T, SkyView = F, FlowAccum = T, soil = F, landcover = F, bedrock = F, LAI = F, mask_only_lc = F, GEOtop = F, MHM = T, rivernet = F, idgauges = F)