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()
