pybert.joint package

Submodules

pybert.joint.dcrefraction module

Created on Fri Jan 23 12:51:21 2015

@author: Marcus

class pybert.joint.dcrefraction.DCRefraction(filename_dc=None, filename_ra=None, verbose=True, **kwargs)

Bases: object

Class for managing a structurally coupled DC resistivity and refraction inversion. Takes care of the standard setup logistics.

Methods

createFOP([refine])	create forward operator working on refined mesh
createInv([verbose, dosave])	create inversion instance
estimateError([absoluteError, relativeError])	estimate error composed of an absolute and a relative part
getDepth()	Return the maximum (roughly estimated) depth of investigation of the two methods.
load(filename_dc, filename_ra)	Load data from files.
makeMesh([depth, quality, max_cell_area])	create (inversion)
run([zweight_dc, zweight_ra, lambda_dc, …])	Run a structurally coupled inversion with an optional separate preparatory inversion.
saveResults([folder, fig])	Saves the results in the specified folder.
setMesh(mesh)	set mesh from outside (e.g.
showData([ax])	Show data in form of apparent resistivities and apparent velocoties.
showMesh([ax])	show mesh in given axes or in a new figure
showResult([c_min_dc, c_max_dc, c_min_ra, …])	Show resulting resistivity and velocity vectors.
showResults([folder, size])	show all results in one figure
standardizedCoverage([threshhold])	Return standardized coverage vector (0|1) using threshholding.

createFOP(refine=True)

create forward operator working on refined mesh

createInv(verbose=True, dosave=False)

create inversion instance

estimateError(absoluteError=0.001, relativeError=0.03)

estimate error composed of an absolute and a relative part

getDepth()

Return the maximum (roughly estimated) depth of investigation of the two methods.

load(filename_dc, filename_ra)

Load data from files.

makeMesh(depth=None, quality=34.3, max_cell_area=0.0)

create (inversion)

run(zweight_dc=0.7, zweight_ra=0.7, lambda_dc=30.0, lambda_ra=30.0, a=0.1, b=0.1, c=1.0, max_iter=10, sep_iter=2, vel_top=500.0, vel_bottom=5000.0, verbose=False)

Run a structurally coupled inversion with an optional separate preparatory inversion.

saveResults(folder=None, fig=None, **kwargs)

Saves the results in the specified folder.

Saved items are:

Inverted profiles Resistivity vector Velocity vector Coverage vectors Standardized coverage vectors

Joint standardized coverage vector Mesh (bms and vtk with results)

setMesh(mesh)

set mesh from outside (e.g. load from bert mesh folder)

showData(ax=None)

Show data in form of apparent resistivities and apparent velocoties.

showMesh(ax=None)

show mesh in given axes or in a new figure

showResult(c_min_dc=None, c_max_dc=None, c_min_ra=None, c_max_ra=None, log_scale_dc=True, log_scale_ra=False, **kwargs)

Show resulting resistivity and velocity vectors.

showResults(folder=None, size=(16, 20), **kwargs)

show all results in one figure

standardizedCoverage(threshhold=0.01)

Return standardized coverage vector (0|1) using threshholding.

Coverage is where we have cells adequately covered by both methods.

pybert.joint.dcrefraction.getcweight(vec, a=0.1, b=0.1, c=1.0)

Structural constraint weight as function of roughness.