Skip to content
GitLab
Select Git revision
  • master default
Find file BlameHistoryPermalink
script_layout.R 5.18 KiB
1 
# \\\
2 
# Copyright 2021-2022 Louis Héraut*1,
3 
#                     Éric Sauquet*2,
4 
#                     Valentin Mansanarez
5 
#
6 
# *1   INRAE, France
7 
#      louis.heraut@inrae.fr
8 
# *2   INRAE, France
9 
#      eric.sauquet@inrae.fr
10 
#
11 
# This file is part of ash R toolbox.
12 
#
13 
# Ash R toolbox is free software: you can redistribute it and/or
14 
# modify it under the terms of the GNU General Public License as
15 
# published by the Free Software Foundation, either version 3 of the
16 
# License, or (at your option) any later version.
17 
#
18 
# Ash R toolbox is distributed in the hope that it will be useful, but
19 
# WITHOUT ANY WARRANTY; without even the implied warranty of
20 
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
21 
# General Public License for more details.
22 
#
23 
# You should have received a copy of the GNU General Public License
24 
# along with ash R toolbox.
25 
# If not, see <https://www.gnu.org/licenses/>.
26 
# ///
27 
#
28 
#
29 
# Rcode/plotting/script_layout.R
30 
#
31 
# Script that manages the call to the right process in order to
32 
# realise plottings of data analyses.
33
34
35 
# Sourcing R files
36 
source(file.path('Rcode', 'processing', 'extract.R'), encoding='UTF-8')
37 
source(file.path('Rcode', 'plotting', 'layout.R'), encoding='UTF-8')
38
39
40 
## 0. SHAPEFILE LOADING ______________________________________________
41 
# Shapefile importation in order to it only once time
42 
df_shapefile = ini_shapefile(resources_path, df_meta,
43 
                             fr_shpdir, fr_shpname,
44 
                             bs_shpdir, bs_shpname,
45 
                             sbs_shpdir, sbs_shpname,
46 
                             cbs_shpdir, cbs_shpname, cbs_coord,
47 
                             rv_shpdir, rv_shpname,
48 
                             show_river=show_river)
49
50
51 
## 1. HYDROMETRIC STATIONS LAYOUT ____________________________________
52 
### 1.1. Flow time series for stations _______________________________
53 
if ('station_serie_plot' %in% to_do) {
54 
    # Square root computation
55 
    df_sqrt = compute_sqrt(df_data)
56 
    # Layout
57 
    layout_panel(to_plot=c('datasheet'),
58 
                 df_meta=df_meta,
59 
                 df_data=list(df_data,
60 
                              df_sqrt),
61 
                 var=list('Q', 'sqrt(Q)'),
62 
                 type=list('data', 'data'),
63 
                 axis_xlim=axis_xlim,
64 
                 layout_matrix=matrix(c(1, 2), ncol=1),
65 
                 summary=TRUE,
66 
                 info_header=df_data,
67 
                 df_shapefile=df_shapefile,
68 
                 figdir=figdir,
69 
                 resources_path=resources_path,
70 
                 logo_dir=logo_dir,
7172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140
PRlogo_file=PRlogo_file, AEAGlogo_file=AEAGlogo_file, INRAElogo_file=INRAElogo_file, FRlogo_file=FRlogo_file) } ### 1.2. Analyses layout _____________________________________________ if ('station_trend_plot' %in% to_do) { layout_panel(to_plot=to_plot_station, map_to_plot=map_to_plot, df_meta=df_meta, df_data=DF_DATA, df_trend=DF_TREND, var=var, type=type, glose=glose, layout_matrix=matrix(c(1, 2, 3, 4, 5), ncol=1), missRect=TRUE, trend_period=trend_period, mean_period=mean_period, colorForce=TRUE, summary=TRUE, info_header=df_data, time_header=df_data, foot_note=TRUE, info_height=2.8, time_ratio=2, var_ratio=3, foot_height=1.25, df_shapefile=df_shapefile, figdir=figdir, filename_opt='', resdir=resdir, resources_path=resources_path, logo_dir=logo_dir, PRlogo_file=PRlogo_file, AEAGlogo_file=AEAGlogo_file, INRAElogo_file=INRAElogo_file, FRlogo_file=FRlogo_file) } ## 2. BREAK LAYOUT ___________________________________________________ if ('station_break_plot' %in% to_do) { # For all the variable for (v in var) { # Gets the break results for the variable df_break = DF_BREAK[[v]] histogram(df_break, df_meta, title=v, figdir=figdir) cumulative(df_break, df_meta, title=v, dyear=8, figdir=figdir) } } ## 3. CLIMATE LAYOUT _________________________________________________ if ('climate_trend_plot' %in% to_do) { layout_panel( to_plot=c('datasheet'), df_meta=df_climate_meta, df_data=list( res_PAtrend$data, res_TAtrend$data, res_ETPAtrend$data ), df_trend=list( res_PAtrend$trend, res_TAtrend$trend, res_ETPAtrend$trend ),
141142143144145146147148149150151152153154155156157
var=var_climate, type=type_climate, glose=glose_climate, layout_matrix=matrix(c(1, 2, 3), ncol=1), missRect=TRUE, trend_period=trend_period, mean_period=mean_period, colorForce=TRUE, linetype_per=c('longdash', 'solid'), info_header='code', info_height=0.5, time_header=NULL, foot_note=FALSE, paper_size=c(21, 18), figdir=figdir) }