ax, cb = ra.showResult(cMin=min(vp), cMax=max(vp), logScale=False) pg.show(geom, ax=ax, fillRegion=False, regionMarker=False) # lines on top ############################################################################### # Note that internally the following is called # # .. code-block:: python # # ax, _ = pg.show(ra.mesh, vest, label="Velocity [m/s]", **kwargs) # ############################################################################### # Another useful tool is to show the model along with its respone on the data_ ra.showResultAndFit() ############################################################################### # Takeaway message: # A default data inversion with checking of the data consists of only few lines # (Everthing else can be looked at by introspecting the Refraction manager) # # .. code-block:: python # # from pygimli.physics import Refraction # ra = Refraction(filename) # ra.invert() # ra.showResultAndFit() pg.wait()
#!/usr/bin/env python
# -*- coding: utf-8 -*-
"""
Example refraction script applied to an example file
(shaly bedrock detection from Günther&Rücker (2006), EAGE Near Surface)
"""
import os
import pygimli as pg
from pygimli.physics import Refraction
# load example file explicitly from same directory (if called from elsewhere)
workpath=os.path.dirname(__file__)
if len(workpath) == 0:
workpath = '.'
ra = Refraction(workpath + '/example_topo.sgt')
# ra = Refraction()
# ra.loadData(os.path.dirname(__file__) + '/example_topo.sgt')
print(ra)
if False: # possible (typical) actions
ra.showData() # show only data (right after init)
ra.showVA() # show apparent velocity image
ra.createMesh(depth=15.) # pass non-default meshing options
ra.invert(lam=300, zWeight=0.1) # use vtop/vbottom for startmodel
ra.showResultAndFit() # typical output: model and data with response
pg.wait()
ax, cb = ra.showResult(cMin=min(vp), cMax=max(vp), logScale=False) pg.show(geom, ax=ax, fillRegion=False, regionMarker=False) # lines on top ############################################################################### # Note that internally the following is called # # .. code-block:: python # # ax, _ = pg.show(ra.mesh, vest, label="Velocity [m/s]", **kwargs) # ############################################################################### # Another useful tool is to show the model along with its respone on the data_ ra.showResultAndFit() ############################################################################### # Takeaway message: # A default data inversion with checking of the data consists of only few lines # (Everthing else can be looked at by introspecting the Refraction manager) # # .. code-block:: python # # from pygimli.physics import Refraction # ra = Refraction(filename) # ra.invert() # ra.showResultAndFit() pg.wait()
#!/usr/bin/env python
# -*- coding: utf-8 -*-
"""
Example refraction script applied to an example file
(shaly bedrock detection from Günther&Rücker (2006), EAGE Near Surface)
"""
import os
import pygimli as pg
from pygimli.physics import Refraction
# load example file explicitly from same directory (if called from elsewhere)
ra = Refraction(os.path.dirname(__file__) + '/example_topo.sgt')
#ra = Refraction()
#ra.loadData(os.path.dirname(__file__) + '/example_topo.sgt')
print(ra)
if False: # possible (typical) actions
ra.showData() # show only data (right after init)
ra.showVA() # show apparent velocity image
ra.createMesh(depth=15.) # pass non-default meshing options
ra.invert(lam=300, zweight=0.1) # use vtop/vbottom for startmodel
ra.showResultAndFit() # typical output: model and data with response
pg.wait()