Skip to content
GitLab
Select Git revision
  • dev default
  • master
  • 136-decreased-performance-of-calibration-execution-time
  • 115-add-pcict-date-format-managment
  • 118-plot-outputsmodel-remove-case-sensitivity-to-the-which-argument
  • 96-parallelize-the-grid-screening-step-during-calibration
  • 60-migrate-transfoparam-functions-and-createcaliboptions-to-s3-methods
  • airGRplus
  • cleanAirGRplus
  • CRAN_v1.7.0
  • CRAN_v1.6.12
  • CRAN_v1.6.11
  • CRAN_v1.6.10.4
  • CRAN_v1.6.9.27
  • CRAN_v1.6.9.24
  • CRAN_v1.6.9.23
  • CRAN_v1.6.9.21
  • CRAN_v1.4.3.65
  • CRAN_v1.4.3.60
  • CRAN_v1.4.3.52
  • CRAN_v1.3.2.42
  • CRAN_v1.3.2.41
  • CRAN_v1.3.2.40
  • CRAN_v1.3.2.39
  • CRAN_v1.3.2.23
  • CRAN_v1.3.2.17
  • CRAN_v1.2.13.16
  • CRAN_v1.2.13.13
  • CRAN_v1.0.14.1
Find file BlameHistoryPermalink
Forked from HYCAR-Hydro / airGR
Source project has a limited visibility.
script.R 7.15 KiB
1 
############ A MODIFIER ############
2
3
4 
# Path to the data
5 
computer_data_path =
6 
    "/home/louis/Documents/bouleau/INRAE/CDD_stationnarite/data"
7 
    # "C:\\Users\\louis.heraut\\Documents\\CDD_stationnarite\\data"
8
9 
# Work path (it needs to end with '/ASH' directory)
10 
computer_work_path =
11 
    "/home/louis/Documents/bouleau/INRAE/CDD_stationnarite/ASH"
12 
    # "C:\\Users\\louis.heraut\\Documents\\CDD_stationnarite\\ASH"
13
14
15 
### BANQUE HYDRO ###
16 
# Path to the directory where BH data is stored from the work path
17 
filedir =
18 
    # ""
19 
    "BanqueHydro_Export2021"
20
21
22 
### MANUAL SELECTION ###
23 
# Name of the file that will be analysed from the AG directory
24 
filename =
25 
    # ""
26
27 
    # c(
28 
      # "S2235610_HYDRO_QJM.txt",
29 
      # "P1712910_HYDRO_QJM.txt",
30 
      # "P0885010_HYDRO_QJM.txt",
31 
      # "O5055010_HYDRO_QJM.txt",
32 
      # "Q7002910_HYDRO_QJM.txt"
33 
      # )
34
35 
    c("S4214010_HYDRO_QJM.txt",
36 
      "O0384010_HYDRO_QJM.txt",
37 
      "Q7002910_HYDRO_QJM.txt")
38
39
40
41 
### AGENCE EAU ADOUR GARONNE SELECTION ###
42 
# Path to the list file of AG data that will be analysed
43 
AGlistdir =
44 
    ""
45
46 
AGlistname =
47 
    ""
48 
    # "Liste-station_RRSE.docx"
49
50
51 
### NIVALE SELECTION ###
52 
# Path to the list file of metadata about station that will be analysed
53 
NVlistdir =
54 
    ""
55
56 
NVlistname =
57 
    ""
58 
    # "nival_selection.txt"
59
60
61 
### TREND ANALYSIS ###
62 
# Time period to analyse
63 
periodAll = c("1800-01-01", "2019-12-31")
64 
periodSub = c("1968-01-01", "2019-12-31")
65 
trend_period = list(periodAll, periodSub)
66
67 
# Time period to mean
68 
period1 = c("1968-01-01", "1994-12-31")
69 
period2 = c("1995-01-01", "2019-12-31")
70 
mean_period = list(period1, period2)
7172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140
# p value p_thresold = 0.1 #c(0.01, 0.05, 0.1) ### MAP ### fr_shpdir = 'map/france' fr_shpname = 'gadm36_FRA_0.shp' bs_shpdir = 'map/bassin' bs_shpname = 'BassinHydrographique.shp' rv_shpdir = 'map/river' rv_shpname = 'CoursEau_FXX.shp' #################################### # FILE STRUCTURE # # Set working directory setwd(computer_work_path) # Sourcing R file source('processing/extract.R', encoding='latin1') source('processing/format.R', encoding='latin1') source('processing/analyse.R', encoding='latin1') source('plotting/panel.R', encoding='latin1') source('plotting/layout.R', encoding='latin1') # Result directory resdir = file.path(computer_work_path, 'results') if (!(file.exists(resdir))) { dir.create(resdir) } print(paste('resdir :', resdir)) # Figure directory figdir = file.path(computer_work_path, 'figures') if (!(file.exists(figdir))) { dir.create(figdir) } print(paste('figdir :', figdir)) # Initialization of null data frame if there is no data selected df_data_AG = NULL df_data_NV = NULL df_meta_AG = NULL df_meta_NV = NULL # AGENCE EAU ADOUR GARONNE SELECTION # if (AGlistname != ""){ # Get only the selected station from a list station file df_selec_AG = get_selection_AG(computer_data_path, AGlistdir, AGlistname, cnames=c('code', 'station', 'BV_km2', 'axe_principal_concerne', 'longueur_serie', 'commentaires', 'choix'), c_num=c('BV_km2', 'longueur_serie'))
141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210
# Get filenames of the selection filename = df_selec_AG[df_selec_AG$ok,]$filename ##### # filename = filename[(1+30):(16+30)] ##### # Extract metadata about selected stations df_meta_AG = extract_meta(computer_data_path, filedir, filename) # Extract data about selected stations df_data_AG = extract_data(computer_data_path, filedir, filename) } # NIVALE SELECTION # if (NVlistname != ""){ # Get only the selected station from a list station file df_selec_NV = get_selection_NV(computer_data_path, NVlistdir, NVlistname) # Get filenames of the selection filename = df_selec_NV[df_selec_NV$ok,]$filename ##### # filename = filename[(1+20):(16+20)] ##### # Extract metadata about selected stations df_meta_NV = extract_meta(computer_data_path, filedir, filename) # Extract data about selected stations df_data_NV = extract_data(computer_data_path, filedir, filename) } # MANUAL SELECTION # if (AGlistname == "" & NVlistname == "") { # Extract metadata about selected stations df_meta_AG = extract_meta(computer_data_path, filedir, filename) # Extract data about selected stations df_data_AG = extract_data(computer_data_path, filedir, filename) } # JOIN # # Make the join df_join = join(df_data_AG, df_data_NV, df_meta_AG, df_meta_NV) df_data = df_join$data df_meta = df_join$meta # ANALYSE # # Compute gap parameters for stations # df_lac = get_lacune(df_data, df_meta) # QA TREND # res_QAtrend = get_QAtrend(df_data, period=trend_period, p_thresold=p_thresold) # QMNA TREND # res_QMNAtrend = get_QMNAtrend(df_data, period=trend_period, p_thresold=p_thresold) # VCN10 TREND # res_VCN10trend = get_VCN10trend(df_data, df_meta, period=trend_period,
211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275
p_thresold=p_thresold) # df_break = get_break(res_QAtrend$data, df_meta) # df_break = get_break(res_QMNAtrend$data, df_meta) # df_break = get_break(res_VCN10trend$data, df_meta) # histogram(df_break$Date, df_meta, # figdir=figdir) # cumulative(df_break$Date, df_meta, dyear=8, # figdir=figdir) # TIME PANEL # # Plot time panel of debit by stations # panels_layout(list(df_data, df_data), # layout_matrix=c(1, 2), # df_meta=df_meta, # missRect=list(TRUE, TRUE), # type=list('Q', 'sqrt(Q)'), # info_header=TRUE, # time_header=NULL, # var_ratio=3, # figdir=figdir, # filename_opt='time') panels_layout(list(res_QAtrend$data, res_QMNAtrend$data, res_VCN10trend$data), layout_matrix=c(1, 2, 3), df_meta=df_meta, df_trend=list(res_QAtrend$trend, res_QMNAtrend$trend, res_VCN10trend$trend), type=list(bquote(Q[A]), bquote(Q[MNA]), bquote(V[CN10])), missRect=list(TRUE, TRUE, TRUE), trend_period=trend_period, mean_period=mean_period, info_header=TRUE, time_header=df_data, info_ratio=2, time_ratio=2, var_ratio=3, computer_data_path=computer_data_path, fr_shpdir=fr_shpdir, fr_shpname=fr_shpname, bs_shpdir=bs_shpdir, bs_shpname=bs_shpname, rv_shpdir=rv_shpdir, rv_shpname=rv_shpname, figdir=figdir, filename_opt='') # map_panel(computer_data_path, # fr_shpdir, fr_shpname, # bs_shpdir, bs_shpname, # rv_shpdir, rv_s hpname) ### /!\ Removed 185 row(s) containing missing values (geom_path) -> remove NA ###