diff --git a/experiment/meteo_france_data/scm_models_data/visualization/hypercube_visualization/abstract_hypercube_visualizer.py b/experiment/meteo_france_data/scm_models_data/visualization/hypercube_visualization/abstract_hypercube_visualizer.py
index 8bd9527e02cdea0706c906d1933583402dfc1743..728090c18bb9547606d45808622fabefb31c2bb0 100644
--- a/experiment/meteo_france_data/scm_models_data/visualization/hypercube_visualization/abstract_hypercube_visualizer.py
+++ b/experiment/meteo_france_data/scm_models_data/visualization/hypercube_visualization/abstract_hypercube_visualizer.py
@@ -25,7 +25,10 @@ class AbstractHypercubeVisualizer(object):
first_starting_year=None,
last_starting_year=None,
exact_starting_year=None,
- verbose=True):
+ verbose=True,
+ sigma_for_best_year=0.0):
+ assert sigma_for_best_year >= 0.0
+ self.sigma_for_best_year = sigma_for_best_year
self.reduce_strength_array = reduce_strength_array
self.verbose = verbose
self.save_to_file = save_to_file
diff --git a/experiment/meteo_france_data/scm_models_data/visualization/hypercube_visualization/altitude_hypercube_visualizer.py b/experiment/meteo_france_data/scm_models_data/visualization/hypercube_visualization/altitude_hypercube_visualizer.py
index 5db8641c49de4bef291c110af1eb2929763df068..eba9f36a612310707fd4176ffc670d1ec5e359e5 100644
--- a/experiment/meteo_france_data/scm_models_data/visualization/hypercube_visualization/altitude_hypercube_visualizer.py
+++ b/experiment/meteo_france_data/scm_models_data/visualization/hypercube_visualization/altitude_hypercube_visualizer.py
@@ -2,6 +2,7 @@ import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
from matplotlib.ticker import FormatStrFormatter, ScalarFormatter
+from scipy.ndimage import gaussian_filter
from experiment.meteo_france_data.scm_models_data.visualization.hypercube_visualization.abstract_hypercube_visualizer import \
AbstractHypercubeVisualizer
@@ -157,7 +158,8 @@ class AltitudeHypercubeVisualizer(AbstractHypercubeVisualizer):
assert isinstance(serie, pd.Series)
xlabel_values = list(serie.index)
values = list(serie.values)
- argmax_idx = np.argmax(values)
+ smooth_values = gaussian_filter(values, self.sigma_for_best_year)
+ argmax_idx = np.argmax(smooth_values)
best_year = xlabel_values[argmax_idx]
if plot_title is not None:
plot_title += ' (max reached in {})'.format(best_year)
@@ -168,6 +170,9 @@ class AltitudeHypercubeVisualizer(AbstractHypercubeVisualizer):
ax.plot([], [], label=ylabel, color=color)
linewidth = 10 if poster_plot else None
ax_reversed.plot(xlabel_values, values, label=ylabel, color=color, linewidth=linewidth)
+ if self.sigma_for_best_year > 0:
+ ax_reversed.plot(xlabel_values, smooth_values, label=ylabel + ' smooth', color=color, linestyle=':',
+ linewidth=linewidth)
fontsize = 30 if poster_plot else None
ax_reversed.set_ylabel(ylabel, color=color, fontsize=fontsize, labelpad=-20)
ax_reversed.axvline(x=best_year, color=color, linestyle='--', linewidth=linewidth)
@@ -344,6 +349,8 @@ class AltitudeHypercubeVisualizer(AbstractHypercubeVisualizer):
if self.reduce_strength_array:
title += '\nEvolution of the quantile {} every {} years'.format(AbstractGevTrendTest.quantile_for_strength,
AbstractGevTrendTest.nb_years_for_quantile_evolution)
+ else:
+ title += '\nStarting years'
if set:
plt.suptitle(title)
return title
diff --git a/experiment/paper1_steps/1 - non stationary model choice/main4_common_spatial_altitude_starting_years_impact.py b/experiment/paper1_steps/1 - non stationary model choice/main4_common_spatial_altitude_starting_years_impact.py
new file mode 100644
index 0000000000000000000000000000000000000000..ac8a63e3512088748050aa26901277f730633ae5
--- /dev/null
+++ b/experiment/paper1_steps/1 - non stationary model choice/main4_common_spatial_altitude_starting_years_impact.py
@@ -0,0 +1,58 @@
+import time
+
+from experiment.meteo_france_data.scm_models_data.visualization.hypercube_visualization.altitude_year_hypercube_visualizer import \
+ Altitude_Hypercube_Year_Visualizer, AltitudeHypercubeVisualizerWithoutTrendType
+from experiment.trend_analysis.univariate_test.gev_trend_test_one_parameter import GevScaleTrendTest, \
+ GevLocationTrendTest
+from experiment.trend_analysis.univariate_test.gev_trend_test_two_parameters import GevLocationAndScaleTrendTest
+
+"""
+Visualize the 0.99 quantile initial value and its evolution
+"""
+from experiment.paper1_steps.utils import get_full_altitude_visualizer, FULL_ALTITUDES
+
+
+def main_fast_spatial_risk_evolution():
+ vizualiser = get_full_altitude_visualizer(AltitudeHypercubeVisualizerWithoutTrendType, altitude=None,
+ reduce_strength_array=True,
+ trend_test_class=GevLocationAndScaleTrendTest,
+ offset_starting_year=28)
+ vizualiser.save_to_file = False
+ vizualiser.sigma_for_best_year = 1.0
+ res = vizualiser.visualize_year_trend_test(subtitle_specified='CrocusSwe3Days')
+ print(res)
+
+
+def main_full_spatial_risk_evolution():
+ for trend_test_class in [GevLocationAndScaleTrendTest]:
+ # Compare the risk with and without taking into account the starting year
+ vizualiser = get_full_altitude_visualizer(AltitudeHypercubeVisualizerWithoutTrendType, altitude=None,
+ reduce_strength_array=True,
+ trend_test_class=trend_test_class,
+ offset_starting_year=20)
+ vizualiser.sigma_for_best_year = 1.0
+ res = vizualiser.visualize_year_trend_test(subtitle_specified='CrocusSwe3Days')
+ best_year = res[0][1]
+ for altitude in FULL_ALTITUDES[:]:
+ # Starting Year=1958
+ vizualiser = get_full_altitude_visualizer(Altitude_Hypercube_Year_Visualizer, altitude=altitude,
+ exact_starting_year=1958, reduce_strength_array=True,
+ trend_test_class=trend_test_class)
+ vizualiser.visualize_massif_trend_test_one_altitude()
+ # Optimal common starting year
+ vizualiser = get_full_altitude_visualizer(Altitude_Hypercube_Year_Visualizer, altitude=altitude,
+ exact_starting_year=best_year, reduce_strength_array=True,
+ trend_test_class=trend_test_class)
+ vizualiser.visualize_massif_trend_test_one_altitude()
+
+
+def main_run():
+ main_full_spatial_risk_evolution()
+ # main_fast_spatial_risk_evolution()
+
+
+if __name__ == '__main__':
+ start = time.time()
+ main_run()
+ duration = time.time() - start
+ print('Full run took {}s'.format(round(duration, 1)))
diff --git a/experiment/paper1_steps/1 - non stationary model choice/main4_common_starting_years_impact.py b/experiment/paper1_steps/1 - non stationary model choice/main4_common_spatial_starting_years_impact.py
similarity index 97%
rename from experiment/paper1_steps/1 - non stationary model choice/main4_common_starting_years_impact.py
rename to experiment/paper1_steps/1 - non stationary model choice/main4_common_spatial_starting_years_impact.py
index 3f08aa4a7bcd2bd58dea631d41c16ea7851b0e79..e07bdd644fbc9fe8109eda50a4c0ce6d4611be08 100644
--- a/experiment/paper1_steps/1 - non stationary model choice/main4_common_starting_years_impact.py
+++ b/experiment/paper1_steps/1 - non stationary model choice/main4_common_spatial_starting_years_impact.py
@@ -29,7 +29,7 @@ def main_fast_spatial_risk_evolution():
def main_full_spatial_risk_evolution():
# Compare the risk with and without taking into account the starting year
for altitude in FULL_ALTITUDES[:]:
- for trend_test_class in [GevLocationAndScaleTrendTest]:
+ for trend_test_class in [GevScaleTrendTest]:
# Starting Year=1958
vizualiser = get_full_altitude_visualizer(Altitude_Hypercube_Year_Visualizer, altitude=altitude,
exact_starting_year=1958, reduce_strength_array=True,
diff --git a/experiment/paper1_steps/1 - non stationary model choice/main4_individual_starting_years_impact.py b/experiment/paper1_steps/1 - non stationary model choice/main4_individual_starting_years_impact.py
index 3c0511308167d8f9f995f6e3d492e9e9da4fa824..22368c51c58319cf5ba31404d5769c203663cc96 100644
--- a/experiment/paper1_steps/1 - non stationary model choice/main4_individual_starting_years_impact.py
+++ b/experiment/paper1_steps/1 - non stationary model choice/main4_individual_starting_years_impact.py
@@ -16,7 +16,8 @@ def main_fast_spatial_risk_evolution():
for altitude in [1800]:
vizualiser = get_full_altitude_visualizer(Altitude_Hypercube_Year_Visualizer, altitude=altitude,
reduce_strength_array=False,
- trend_test_class=GevLocationAndScaleTrendTest)
+ trend_test_class=GevLocationAndScaleTrendTest,
+ offset_starting_year=20)
vizualiser.save_to_file = False
vizualiser.visualize_massif_trend_test_one_altitude()
vizualiser.reduce_strength_array = True
@@ -25,7 +26,7 @@ def main_fast_spatial_risk_evolution():
def main_full_spatial_risk_evolution():
# Compare the risk with and without taking into account the starting year
- for altitude in FULL_ALTITUDES[:]:
+ for altitude in FULL_ALTITUDES[-2:]:
for trend_test_class in [GevLocationAndScaleTrendTest]:
vizualiser = get_full_altitude_visualizer(Altitude_Hypercube_Year_Visualizer, altitude=altitude,
exact_starting_year=1958, reduce_strength_array=True,
@@ -33,7 +34,8 @@ def main_full_spatial_risk_evolution():
vizualiser.visualize_massif_trend_test_one_altitude()
vizualiser = get_full_altitude_visualizer(Altitude_Hypercube_Year_Visualizer, altitude=altitude,
reduce_strength_array=True,
- trend_test_class=trend_test_class)
+ trend_test_class=trend_test_class,
+ offset_starting_year=20)
vizualiser.visualize_massif_trend_test_one_altitude()
vizualiser.reduce_strength_array = False
vizualiser.visualize_massif_trend_test_one_altitude()