update wind_wave_dataset

efekhari27 · Apr 1, 2023 · 07271ad · 07271ad
1 parent 30c1fdf
commit 07271ad
Show file tree

Hide file tree

Showing 7 changed files with 9,998 additions and 1,005 deletions.
diff --git a/examples/data/wind_waves_ANEMOC_1H.csv b/examples/data/wind_waves_ANEMOC_1H.csv
diff --git a/examples/figures/wind_waves.jpg b/examples/figures/wind_waves.jpg
diff --git a/examples/figures/wind_waves_contours.jpg b/examples/figures/wind_waves_contours.jpg
diff --git a/examples/figures/wind_waves_simulated_10000.jpg b/examples/figures/wind_waves_simulated_10000.jpg
diff --git a/examples/figures/wind_waves_woutput.jpg b/examples/figures/wind_waves_woutput.jpg
diff --git a/examples/wind_waves_dataset.py b/examples/wind_waves_dataset.py
@@ -13,6 +13,7 @@
 
 # 
 data = pd.read_csv("data/wind_waves_ANEMOC_1H.csv", index_col=0)
+data = data.iloc[:1000]
 
 # Draw static copulogram on data
 data.columns = ["$\\theta_{wind}$", "$U$", "$\\theta_{wave}$", "$H_s$", "$T_p$"]

diff --git a/examples/wind_waves_nonparametric_fit.py b/examples/wind_waves_nonparametric_fit.py
@@ -12,17 +12,9 @@
 import copulogram as cp
 
 
-df_processed = pd.read_csv('data/df_ANEMOC_1H.csv')
-all_indexes = df_processed.index.values
-saved_columns = ['θ_wind (deg)', 'U_hub (m/s)', 'θ_wave_new (deg)', 'Hs (m)', 'Tp (s)']
-df = df_processed.loc[:, saved_columns]
-
-# Monte Carlo sub-sample to make the plots and the fits faster
-N = 10000
-subsampled_indexes = np.random.choice(all_indexes, N, replace=False)
-data = df.loc[subsampled_indexes, saved_columns]
-new_columns = ["$\\theta_{wind}$", "$U$", "$\\theta_{wave}$", "$H_s$", "$T_p$"]
-data.columns = new_columns
+data = pd.read_csv('data/wind_waves_ANEMOC_1H.csv', index_col=0)
+data.columns = ["$\\theta_{wind}$", "$U$", "$\\theta_{wave}$", "$H_s$", "$T_p$"]
+N = data.shape[0]
 # Draw static copulogram on data
 copulogram = cp.Copulogram(data, latex=True)
 alpha = 0.1
@@ -43,7 +35,7 @@
 
 WaveDir_data = ot.Sample(data["$\\theta_{wave}$"].values.reshape(-1, 1))
 kernel = ot.KernelSmoothing()
-mixed_kw = kernel.computeMixedBandwidth(WaveDir_data[:1000])
+mixed_kw = kernel.computeMixedBandwidth(WaveDir_data[:2000])
 WaveDir_kde = kernel.build(WaveDir_data[:2000], mixed_kw)
 WaveDir_kde = ot.TruncatedDistribution(WaveDir_kde, 0., 365.)