else: p1 = [(northing[0]), 0] p2 = [(northing[0]), 20e3] x_axis = slicedir[1] # direction of p1p2 profil SI = 2 pEXP.plot_line(xp, yp, gravity, p1, p2, interp=True, Xaxis=x_axis) pEXP.plot_field(xp, yp, gravity, shape_grid) #%% # Upward continuation of the field data with discretisation in altitude controlled by the number of layers (nlay) and the maximum elevation desired (max_elevation) mesh, label_prop = dEXP.upwc(xp, yp, zp, gravity, shape_grid, zmin=0, zmax=max_elevation, nlayers=nlay, qorder=qorder) plt, cmap = pEXP.plot_xy(mesh, label=label_prop, Xaxis=x_axis, p1p2=np.array([p1, p2])) plt.colorbar(cmap) #%% xderiv = transform.derivx(xp, yp, gravity, shape_grid, order=qorder) yderiv = transform.derivy(xp, yp, gravity, shape_grid, order=qorder)
xp, yp, U, shape = dEXP.pad_edges(xp, yp, U, shape, pad_type=0) # reflexion=5 p1 = [min(yp), 0] p2 = [max(yp), 0] x_axis = 'y' pEXP.plot_line(xp, yp, U, p1, p2, interp=interp, Xaxis=x_axis) #%% # Upward continuation of the field data mesh, label_prop = dEXP.upwc(xp, yp, zp, U, shape, zmin=0, zmax=max_elevation, nlayers=nlay, qorder=qorder) 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')
x_resolution=interp_size, savefig=False, interp=interp, smooth=smooth, Xaxis=x_axis) plt.savefig('zderiv' + str(file) + '.png', dpi=450) #%% ------- upward continuation of the field data p = [p1_s, p2_s] mesh, label_prop = dEXP.upwc(XFs, YFs, zp, UF, shape, zmin=0, zmax=max_elevation, nlayers=nlay, qorder=qorder) plt, cmap = pEXP.plot_xy(mesh, label=label_prop, Xaxis=x_axis, Vminmax=[0, 0.0125], p1p2=p) cbar = plt.colorbar(cmap, shrink=0.25, pad=0.04) cbar.set_label('upwc voltage (V)') plt.tight_layout() plt.savefig('upwc voltage' + str(file) + '.png', dpi=450)