With our environment is set up, and the model downloaded, we can begin building the Salvus mesh. Most material models in seismic exploration are discretized onto a regular grid, and the open-source
SalvusPythonToolkit provides a simple interface to create two or three dimensional Salvus meshes based on such models. As models can be stored in a wide variety of format such as
SEG-Y, or
RSF, we use
xarray
as an intermediary to encapsulate both the model parameters and geometry. More info on
xarray
, including extensive documentation, can be found
here. Reading our Marmousi model into an
xarray.Dataset
is trivial, and the process can be inspected by opening the acompanying
getting_started_tools.py
file.
In general, if you are working working with a regularly gridded model in either two or three dimensions, a recipie similar to the function provided should be all you need to do to get your model into a form Salvus Mesh can understand. All that is required is an xarray
dataset with the following information