Skip to content
GitLab
Select Git revision
  • master default
Find file BlameHistoryPermalink
script.R 4.04 KiB
1
2 
###### A MODIFIER ######
3
4
5 
# Path to the data
6 
computer_data_path =
7 
    "/home/louis/Documents/bouleau/INRAE/CDD_stationnarite/data"
8 
    # "C:\\Users\\louis.heraut\\Documents\\CDD_stationnarite\\data"
9
10 
# Work path
11 
computer_work_path =
12 
    "/home/louis/Documents/bouleau/INRAE/CDD_stationnarite/ASH"
13 
    # "C:\\Users\\louis.heraut\\Documents\\CDD_stationnarite\\ASH"
14
15
16 
### BANQUE HYDRO ###
17 
# Path to the directory where BH data is stored
18 
BHfiledir =
19 
    # ""
20 
    "BanqueHydro_Export2021"
21
22 
## Manual selection ##
23 
# Name of the file that will be analysed from the BH directory
24 
BHfilename =
25 
    # ""
26 
    c("H5920011_HYDRO_QJM.txt", "K4470010_HYDRO_QJM.txt")
27 
    # "all"
28
29 
## Or list selection ##
30 
# Path to the list file of BH data that will be analysed
31 
BHlistdir =
32 
    ""
33
34 
BHlistname =
35 
    ""
36 
    # "Liste-station_RRSE.docx"
37
38
39 
### NIVALE ###
40 
# Path to the directory where NV data is stored
41 
NVfiledir =
42 
    ""
43 
    # "France207"
44
45 
# Name of the file that will be analysed from the NV directory
46 
NVfilename =
47 
    ""
48 
    # "all"
49
50
51 
# Path to the list file of metadata about station that will be analysed
52 
NVlistdir =
53 
    ""
54
55 
NVlistname =
56 
    ""
57 
    # "liste_bv_principaux_global.txt"
58
59
60 
### TREND ANALYSIS ###
61 
# Time period to analyse
62 
period = c("1960-01-01","2019-12-31")
63 
# period = c("1960-01-01","2020-01-01")
64
65
66
67 
########################
68
69
70 
# FILE STRUCTURE #
7172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140
# Set working directory setwd(computer_work_path) # Sourcing R file source('processing/extractBH.R') source('processing/extractNV.R') source('processing/format.R') source('processing/analyse.R') source('plotting/panel.R') # Usefull library # 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)) # BANQUE HYDRO # # Get only the selected station from a list station file if (BHlistname != ""){ df_selec = get_selection(computer_data_path, BHlistdir, BHlistname, cnames=c('code', 'station', 'BV_km2', 'axe_principal_concerne', 'longueur_serie', 'commentaires', 'choix'), c_num=c('BV_km2', 'longueur_serie')) BHfilename = df_selec[df_selec$ok,]$filename } # Extract metadata about selected stations df_meta_BH = extractBH_meta(computer_data_path, BHfiledir, BHfilename) # Extract data about selected stations df_data_BH = extractBH_data(computer_data_path, BHfiledir, BHfilename) # NIVALE # # Extract metadata about selected stations df_meta_NV = extractNVlist_meta(computer_data_path, NVfiledir, NVlistdir, NVlistname) ### /!\ missing station info ### # Extract data about selected stations df_data_NV = extractNV_data(computer_data_path, NVfiledir, NVfilename) # JOIN # df_join = join(df_data_BH, df_data_NV, df_meta_BH, df_meta_NV) df_data = df_join$data df_meta = df_join$meta # TIME PANEL #
141142143144145146147148149150151152153154155156157158159160161162163164165166167168169
# Plot time panel of debit by stations # panel(df_data, df_meta, figdir, "") # panel(df_data, df_meta, figdir, "", is_sqrt=TRUE) ### /!\ Removed 185 row(s) containing missing values (geom_path) -> remove NA ### # ANALYSE # # Compute gap parameters for stations # df_lac = get_lacune(df_data, df_meta) # QA TREND # res_QAtrend = get_QAtrend(df_data, period) panel(df_data=res_QAtrend$data, df_meta=df_meta, df_trend=res_QAtrend$trend, p_threshold=0.1, figdir=figdir, filedir_opt='QA', variable='QA') # QMNA TREND # # res_QMNAtrend = get_QMNAtrend(df_data, period) # VCN10 TREND # # res_VCN10trend = get_VCN10trend(df_data, df_meta, period) # panel(df_data=res_VCN10trend$data, # df_meta=df_meta, # df_trend=res_VCN10trend$trend, # figdir=figdir, filedir="VCN10")