Merge branch 'main' of https://gitlab.cesbio.omp.eu/auclairj/modspa_pixel into main

SAMIR_xls/SAMIRxls_Reference_ReadAndPlot.py

@@ -12,7 +12,7 @@ import numpy as np
 
 # Modspa excel file - Reference Simmonneaux !
 # xls_file_path = './SAMIR_xls/SAMIRpixel_Reference_Simonneaux2012.xls'
-xls_file_path = '~/WORKSPACE/MODSPA/Lamasquere/SAMIRxls_Reference.FR-Lam2019.xls'
+xls_file_path = 'SAMIRxls_Reference.FR-Lam2019.xls'
 titre = 'Irrigated Maize - France (Lam-2019)'
 figname = '-Maize-Lam-2019'
 
@@ -145,8 +145,9 @@ plt.suptitle(titre)
 plt.savefig('SWC-Ks'+figname+'.1.png')
 
 
-# %%
+# %% Etude entre rapport ETRobs/ET0 et Kc
 
+plt.figure(7)
 df_samir_xls['Kcobs'] = df_samir_xls['ETobs_ET0'].rolling(window=5).mean()
 
 plt.figure()
@@ -158,3 +159,14 @@ plt.plot(dates, df_samir_xls.Kcb*df_samir_xls.Ks+(df_samir_xls.Kei+df_samir_xls.
 plt.legend(['Kcobs', 'NDVI', 'Kcobsmean', 'KcbxKs', 'KcbxKs+Ke', 'Kc_obs=ETobs/ET0'], loc='best')
 plt.xticks(rotation=20)
 plt.title(titre)
+
+# %% ETRobs vs ETRsim
+
+plt.figure()
+plt.plot(df_samir_xls.ET, df_samir_xls.ETobs, 'ok', markersize=3)
+plt.plot([0, 10], [0, 10], 'r')
+plt.xlabel("ETR sim mm/d")
+plt.ylabel("ETR obs mm/d")
+plt.text(0, 6, "rmse="+str(round(rmse_calc(df_samir_xls.ET, df_samir_xls.ETobs), 2)))
+plt.title(titre)
+plt.axis('square')