plt, cmap = pEXP.plot_xy(mesh, label=label_prop, Xaxis=x_axis) plt.colorbar(cmap) #%% # ridges identification # dEXP.ridges_minmax_plot(xp, yp, mesh, p1, p2, # label=label_prop, # fix_peak_nb=2, # method_peak='find_peaks') # or find_peaks or peakdet or spline_roots dfI, dfII, dfIII = dEXP.ridges_minmax(xp, yp, mesh, p1, p2, label=label_prop, fix_peak_nb=2, method_peak='find_peaks', showfig=True, Xaxis=x_axis) #%% # Plot ridges over continuated section fig = plt.figure() ax = plt.gca() pEXP.plot_xy(mesh, label=label_prop, ax=ax, Xaxis=x_axis) pEXP.plot_ridges_harmonic(dfI, dfII, dfIII, ax=ax) #%% # Filter ridges regionally constrainsted)
# plt, cmap = pEXP.plot_xy(mesh, label=label_prop) # plt.colorbar(cmap) #%% # Ridges identification # dEXP.ridges_minmax_plot(xp, yp, mesh, p1, p2, # label=label_prop, # fix_peak_nb=2, # method_peak='find_peaks') # or find_peaks or peakdet or spline_roots dfI, dfII, dfIII = dEXP.ridges_minmax(xp, yp, mesh, p1, p2, label=label_prop, fix_peak_nb=2, method_peak='find_peaks') #%% # Plot ridges over continuated section # fig = plt.figure() # ax = plt.gca() # pEXP.plot_xy(mesh, label=label_prop, ax=ax) #, ldg=) # pEXP.plot_ridges_harmonic(dfI,dfII,dfIII,ax=ax) #%% # Filter ridges (regionally constrainsted)
label=label_prop, method_peak='find_peaks', showfig=True, interp=True, smooth=True, Xaxis=x_axis) #%% # Ridges identification at all levels: plot extremas obtained via find_peaks function (numpy) for all 3 types of extremas familly RI, RII and RIII D = dEXP.ridges_minmax(xp, yp, mesh, p1, p2, label=label_prop, method_peak='find_peaks', qorder=qorder, interp=True, smooth=True, fix_peak_nb=1, returnAmp=True, showfig=True, Xaxis=x_axis) dfI, dfII, dfIII = D[0:3] hI, hII, hIII = D[3:6] H = D[3:6] #%% # filter ridges using a minimum length criterium and and filter for a specific range of altitude # a =2.25 # if x_axis=='y': # xf_min = a*x1
# plt, cmap = pEXP.plot_xy(mesh, label=label_prop) # plt.colorbar(cmap) #%% # Ridges identification # dEXP.ridges_minmax_plot(xp, yp, mesh, p1, p2, # label=label_prop, # fix_peak_nb=2, # method_peak='find_peaks') # or find_peaks or peakdet or spline_roots dfI, dfII, dfIII = dEXP.ridges_minmax(xp, yp, mesh, p1, p2, label=label_prop, fix_peak_nb=fix_peak_nb, method_peak='find_peaks', smooth=smooth, showfig=True) #%% # Plot ridges over continuated section # fig = plt.figure() # ax = plt.gca() # pEXP.plot_xy(mesh, label=label_prop, ax=ax) #, ldg=) # pEXP.plot_ridges_harmonic(dfI,dfII,dfIII,ax=ax) #%% # Filter ridges (regionally constrainsted)
# or find_peaks or peakdet or spline_roots # dfI,dfII, dfIII = dEXP.ridges_minmax(Xs, Ys, mesh, p1_s, p2_s,interp=interp,x_resolution= interp_size, # label=label_prop,fix_peak_nb=2, # smooth=smooth, # true = low pass, otherwise specify the filter to apply # method_peak='find_peaks', # showfig=True, # Xaxis=x_axis) D = dEXP.ridges_minmax(Xs, Ys, mesh, p1_s, p2_s, label=label_prop, method_peak='find_peaks', fix_peak_nb=2, returnAmp=True, showfig=True, Xaxis=x_axis, interp=interp, x_resolution=interp_size, smooth=smooth, qorder=qorder) dfI, dfII, dfIII = D[0:3] hI, hII, hIII = D[3:6] heights = D[3:6] #%% ------------------------------- plot ridges over continuated section fig = plt.figure()
# %% ridges identification # dEXP.ridges_minmax_plot(xp, yp, mesh, p1, p2, # label=label_prop, # fix_peak_nb=4, # interp=interp,smooth=smooth, # method_peak='find_peaks') # or find_peaks or peakdet or spline_roots dfI, dfII, dfIII = dEXP.ridges_minmax(Xs, Ys, mesh, p1_s, p2_s, interp=interp, label=label_prop, fix_peak_nb=2, smooth=smooth, method_peak='find_peaks', showfig=True, Xaxis=x_axis) # dfI, dfII, dfIII = dEXP.ridges_minmax(xp, yp, mesh, p1, p2, # label=label_prop, # # minAlt_ridge=minAlt_ridge, # maxAlt_ridge=maxAlt_ridge, # fix_peak_nb=None) #%% ------------------------------- plot ridges over continuated section fig = plt.figure()