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)
dfI_f, dfII_f, dfIII_f = dEXP.filter_ridges(dfI,
dfII,
dfIII,
minDepth=1000,
maxDepth=3000,
minlength=3,
rmvNaN=True)
df_f = dfI_f, dfII_f, dfIII_f
# df_f = dfI, dfII, dfIII
#%%
MESH.append(mesh)
LABEL.append(label_prop)
DF_F.append(df_f)
DF_FIT.append(df_fit)
XXZZ.append(xxzz)
CTm.append(CT)
#%%
plt.figure()
ax = plt.gca()
i = 0
pEXP.plot_xy(MESH[i], label=LABEL[i], ax=ax) #, ldg=)
dfI_f, dfII_f, dfIII_f = DF_F[i]
pEXP.plot_ridges_harmonic(dfI_f, dfII_f, dfIII_f, ax=ax, label=False)
pEXP.plot_ridges_sources(DF_FIT[i],
ax=ax,
z_max_source=-max_elevation * 1.2,
ridge_type=[0, 1, 2],
ridge_nb=None)
x1, x2, z1, z2 = XXZZ[i]
square([x1, x2, z1, z2])
plt.annotate(CTm[i], [(x1 + x2) / 2, (z1 + z2) / 2])
plt.figure()
ax = plt.gca()
i = 1
pEXP.plot_xy(MESH[i], label=LABEL[i], ax=ax) #, ldg=)
dfI_f, dfII_f, dfIII_f = DF_F[i]
