Skip to content
GitLab
Select Git revision
  • master default
Find file BlameHistoryPermalink
script.R 5.35 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("S2235610_HYDRO_QJM.txt",
27 
      # "P1712910_HYDRO_QJM.txt",
28 
      # "P0885010_HYDRO_QJM.txt",
29 
      # "A1000030_HYDRO_QJM.txt",
30 
      # "A2250310_HYDRO_QJM.txt"
31 
      # )
32
33 
## Or list selection ##
34 
# Path to the list file of BH data that will be analysed
35 
BHlistdir =
36 
    ""
37
38 
BHlistname =
39 
    # ""
40 
    "Liste-station_RRSE.docx"
41
42
43 
### NIVALE ###
44 
# Path to the directory where NV data is stored
45 
NVfiledir =
46 
    # ""
47 
    "France207"
48
49 
# Name of the file that will be analysed from the NV directory
50 
NVfilename =
51 
    # ""
52 
    "all"
53
54
55 
# Path to the list file of metadata about station that will be analysed
56 
NVlistdir =
57 
    ""
58
59 
NVlistname =
60 
    # ""
61 
    "liste_bv_principaux_global.txt"
62
63
64 
### TREND ANALYSIS ###
65 
# Time period to analyse
66 
period = c("1980-01-01","2019-12-31")
67
68
69
70 
########################
7172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140
# FILE STRUCTURE # # Set working directory setwd(computer_work_path) # Sourcing R file source('processing/extractBH.R', encoding='latin1') source('processing/extractNV.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') # 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 } ###### BHfilename = BHfilename[1:10] ###### # 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
141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210
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 # ANALYSE # # Compute gap parameters for stations # df_lac = get_lacune(df_data, df_meta) # QA TREND # res_QAtrend = get_QAtrend(df_data, period) # QMNA TREND # res_QMNAtrend = get_QMNAtrend(df_data, period) # VCN10 TREND # res_VCN10trend = get_VCN10trend(df_data, df_meta, period) # 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, # header_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), period=period, info_header=TRUE, time_header=df_data, header_ratio=2, figdir=figdir, filename_opt='') # panels_layout(list(res_QAtrend$data, res_VCN10trend$data), # layout_matrix=c(1, 2), # df_meta=df_meta, # df_trend=list(res_QAtrend$trend, # res_VCN10trend$trend), # type=list(bquote(Q[A]), bquote(V[CN10])), # missRect=list(TRUE, TRUE), # period=period, # info_header=TRUE, # time_header=df_data, # header_ratio=2, # figdir=figdir, # filename_opt='') ### /!\ Removed 185 row(s) containing missing values (geom_path) -> remove NA ###
211212213214