This notebook demonstrates a simplified and cosequently less flexible way to extract a resistivity profile from a ModEM 3D volume, along a 2D seismic line, and plot both the resistivity and seismic data as overlays. Please refer to plot_resistivity_seismic.ipynb for more flexible options.
In [1]:
%%capture
# Add mtpy folder to python path. This may not be necessary
# depending on how mtpy was installed.
import sys
#sys.path.append('/path/to/mtpy')
sys.path.append('/media/data/work/GA/ausLAMP/codes/mtGeoMtpy/')
from mtpy.modeling.modem.plot_slices import PlotSlices
%matplotlib inline
In [2]:
%%capture
# sygy_fn and stacking_velocity_fn are file names for 2D segy data
# and associated stacking velocities, respoectively.
segy_fn = '../data/seismic/seismic.sgy'
stacking_velocity_fn = '../data/seismic/stacking_velocities.txt'
# modem model and data file names
m_fn = '../data/ModEM_files/Modular_MPI_NLCG_056_im2.rho'
d_fn = '../data/ModEM_files/ModEM_Data_im2.dat'
# instantiate PlotSlices class; note that we want coordinates in metres, as in the seismic line
ps = PlotSlices(model_fn=m_fn, data_fn=d_fn, plot_yn='n', map_scale='m')
In [4]:
fig, ax = ps.plot_resistivity_on_seismic(segy_fn, velocity_model=stacking_velocity_fn, pick_every=1)