This page was generated from jupyter-notebook/nb_3d-surveys.ipynb. Interactive online version: Binder badge. Download notebook.

3D visualization

EMagPy can also display 3D inverted survey (inverted based on 1D inversion). It uses the Python package pyvista to render the .vtk produced.

[1]:
import sys
sys.path.append('../src/')
from emagpy import Problem
import pyvista as pv
[6]:
k = Problem()
k.createSurvey('../src/examples/cover-crop/coverCrop.csv')
k.setInit(depths0=[0.5], fixedDepths=[False])
k.invert(method='ROPE', njobs=1)
k.show3D()
Removing 1 NaN from survey
120/120 inverted

EMagPy can also smooth the 1D inversion in 3D, similarly to lateral smoothing for 2D. For this set the threed=True argument and specify a smoothing value using the beta parameter. EMagPy will penalize the inversion from the closest 1D survey (in 3D distance).

[7]:
k = Problem()
k.createSurvey('../src/examples/saprolite/regolith.csv')
k.convertFromNMEA()
k.invert(threed=True, beta=0.5)
k.show3D()
31/31 inverted
[ ]: