An error occurred while loading the file. Please try again.
-
Gaetano Raffaele authored957389e4
# \\\
# Copyright 2021-2022 Louis Héraut*1,
# Éric Sauquet*2,
# Valentin Mansanarez
#
# *1 INRAE, France
# louis.heraut@inrae.fr
# *2 INRAE, France
# eric.sauquet@inrae.fr
#
# This file is part of ash R toolbox.
#
# Ash R toolbox is free software: you can redistribute it and/or
# modify it under the terms of the GNU General Public License as
# published by the Free Software Foundation, either version 3 of the
# License, or (at your option) any later version.
#
# Ash R toolbox is distributed in the hope that it will be useful, but
# WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
# General Public License for more details.
#
# You should have received a copy of the GNU General Public License
# along with ash R toolbox.
# If not, see <https://www.gnu.org/licenses/>.
# ///
#
#
# Rcode/plotting/script_layout.R
#
# Script that manages the call to the right process in order to
# realise plottings of data analyses.
# Sourcing R files
source(file.path('Rcode', 'processing', 'extract.R'), encoding='UTF-8')
source(file.path('Rcode', 'plotting', 'layout.R'), encoding='UTF-8')
## 0. SHAPEFILE LOADING ______________________________________________
# Shapefile importation in order to it only once time
df_shapefile = ini_shapefile(resources_path, df_meta,
fr_shpdir, fr_shpname,
bs_shpdir, bs_shpname,
sbs_shpdir, sbs_shpname,
cbs_shpdir, cbs_shpname, cbs_coord,
rv_shpdir, rv_shpname,
show_river=show_river)
## 1. HYDROMETRIC STATIONS LAYOUT ____________________________________
### 1.1. Flow time series for stations _______________________________
if ('station_serie_plot' %in% to_do) {
# Square root computation
df_sqrt = compute_sqrt(df_data)
# Layout
layout_panel(to_plot=c('datasheet'),
df_meta=df_meta,
df_data=list(df_data,
df_sqrt),
var=list('Q', 'sqrt(Q)'),
type=list('data', 'data'),
axis_xlim=axis_xlim,
layout_matrix=matrix(c(1, 2), ncol=1),
summary=TRUE,
info_header=df_data,
df_shapefile=df_shapefile,
figdir=figdir,
resources_path=resources_path,
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)
}