# differences are minor, but lets keep using Pythons function # my_savitzky_pred = rc.savitzky_golay(y, window_size=Sav_win_size, order=sav_order) savitzky_pred = scipy.signal.savgol_filter(y, window_length= Sav_win_size, polyorder=sav_order) ############################################# ### ### find peaks ### ############################################# ################################################################################# # # savitzky # savitzky_max_min = rc.my_peakdetect(y_axis=savitzky_pred, x_axis=X, delta=delt); savitzky_max = savitzky_max_min[0]; savitzky_min = savitzky_max_min[1]; savitzky_max = rc.separate_x_and_y(m_list = savitzky_max); savitzky_min = rc.separate_x_and_y(m_list = savitzky_min); savitzky_max_DoYs_series = pd.Series(savitzky_max[0]); savitzky_max_series = pd.Series(savitzky_max[1]); savitzky_min_DoYs_series = pd.Series(savitzky_min[0]); savitzky_min_series = pd.Series(savitzky_min[1]); savitzky_max_df = pd.DataFrame({
SG_pred_93 = scipy.signal.savgol_filter(y, window_length= 9, polyorder=3) # preds_df = pd.DataFrame(data = {'SG 31':SG_pred_31, 'SG 32':SG_pred_32, # 'SG 51':SG_pred_51, 'SG 52':SG_pred_52, 'SG 53':SG_pred_53, # 'SG 71':SG_pred_71, 'SG 72':SG_pred_72, 'SG 73':SG_pred_73, # 'SG 91':SG_pred_91, 'SG 92':SG_pred_92, 'SG 93':SG_pred_93}, # index = X) # df2.set_index(keys = ['a'], drop=True, inplace=True) ############################################# ### ### find peaks ### ############################################# SG_max_min_31 = rc.my_peakdetect(y_axis=SG_pred_31, x_axis=X, delta=delt); SG_max_31 = SG_max_min_31[0]; SG_min_31 = SG_max_min_31[1]; SG_max_31 = rc.separate_x_and_y(m_list = SG_max_31); SG_min_31 = rc.separate_x_and_y(m_list = SG_min_31); SG_max_DoYs_series_31 = pd.Series(SG_max_31[0]); SG_max_series_31 = pd.Series(SG_max_31[1]); SG_min_DoYs_series_31 = pd.Series(SG_min_31[0]); SG_min_series_31 = pd.Series(SG_min_31[1]); SG_max_min_32 = rc.my_peakdetect(y_axis=SG_pred_32, x_axis=X, delta=delt); SG_max_32 = SG_max_min_32[0]; SG_min_32 = SG_max_min_32[1]; SG_max_32 = rc.separate_x_and_y(m_list = SG_max_32); SG_min_32 = rc.separate_x_and_y(m_list = SG_min_32); SG_max_DoYs_series_32 = pd.Series(SG_max_32[0]); SG_max_series_32 = pd.Series(SG_max_32[1]);