pybert.tdip package

Time-domain induced polarization (TDIP) Data Manager.

class pybert.tdip.TDIPdata(filename=None, **kwargs)[source]

Bases: object

Class managing time-domain induced polarisation (TDIP) field data.

Methods

ensureRhoa() Make sure apparent resistivity is present in file.
filter([tmin, tmax, kmax, electrode, …]) Filter data with respect to frequencies and geometric factor.
generatePDF([rdict, mdict]) Generate a multi-page pdf file with all data as pseudosections.
individualInversion(**kwargs) Carry out individual inversion for spectral chargeability.
invertMa([nr, ma]) Invert for chargeability.
invertRhoa(**kwargs) Invert apparent resistivity values.
load(filename, **kwargs) Try loading all supported file types.
save([filename]) Save all data in some (yet-to-be-defined or -decided) format.
showDecay([nr, ab, mn, verbose]) Show decay curves.
showMa([nr]) Show apparent chargeability.
showRhoa(**kwargs) Show apparent resistivity.
ensureRhoa()[source]

Make sure apparent resistivity is present in file.

filter(tmin=0, tmax=1000000000.0, kmax=1000000.0, electrode=None, forward=False, a=None, b=None, m=None, n=None, ab=None, mn=None, corrSID=1, nr=[])[source]

Filter data with respect to frequencies and geometric factor.

Parameters:
tmin : double

minimum frequency

tmax : double

maximum frequency

kmax : double

maximum (absolute) geometric factor

electrode : int

electrode to be removed completely

a/b/m/n : int

delete data with specific current or potential dipole lengths

ab/mn : int

delete data with specific current or potential dipole lengths

corrSID: int [1]

correct sensor index (like in data files)
generatePDF(rdict=None, mdict=None, **kwargs)[source]

Generate a multi-page pdf file with all data as pseudosections.

individualInversion(**kwargs)[source]

Carry out individual inversion for spectral chargeability.

invertMa(nr=0, ma=None, **kwargs)[source]

Invert for chargeability.

invertRhoa(**kwargs)[source]

Invert apparent resistivity values.

load(filename, **kwargs)[source]

Try loading all supported file types.

save(filename=None)[source]

Save all data in some (yet-to-be-defined or -decided) format.

showDecay(nr=[], ab=None, mn=None, verbose=True)[source]

Show decay curves.

showMa(nr=0, **kwargs)[source]

Show apparent chargeability.

showRhoa(**kwargs)[source]

Show apparent resistivity.

Submodules

pybert.tdip.example module

pybert.tdip.mipmodelling module

class pybert.tdip.mipmodelling.DCIPMModelling(f, mesh, rho, verbose=False)[source]

Bases: pygimli.core.ModellingBaseMT__

DC/IP modelling class using an (FD-based) approach

Methods

__call__((object)arg1, (object)model) C++ signature :
clearConstraints((object)arg1) C++ signature :
clearJacobian((object)arg1) C++ signature :
constraints((object)arg1) C++ signature :
constraintsRef((object)arg1) C++ signature :
createConstraints((object)arg1) C++ signature :
createDefaultStartModel((object)arg1) C++ signature :
createJacobian(model) create jacobian matrix using unchanged DC jacobian and m model
createMappedModel((object)arg1, …) Read only extrapolation of model values given per cell marker to values given per cell.
createRefinedForwardMesh((object)arg1 [, …) C++ signature :
createStartModel((object)arg1) C++ signature :
createStartVector((object)arg1) DEPRECATED use createStartModel
data((object)arg1) Return the associated data container.
deleteMesh((object)arg1) Delete the actual mesh.
initConstraints((object)arg1) C++ signature :
initJacobian((object)arg1) C++ signature :
initRegionManager((object)arg1) C++ signature :
jacobian((object)arg1) Return the pointer to the Jacobian matrix associated with this forward operator.
jacobianRef((object)arg1) C++ signature :
mapModel((object)arg1, (object)model [, …) C++ signature :
mesh((object)arg1) C++ signature :
multiThreadJacobian((object)arg1) Return number of threads used for Jacobian generation.
region((object)arg1, (object)marker) Syntactic sugar for this->regionManager().region(marker).
regionManager((object)arg1) C++ signature :
regionManagerRef((object)arg1) C++ signature :
response(m) return forward response as function of chargeability model
response_mt((object)arg1, (object)model [, …) C++ signature :
setConstraints((object)arg1, (object)C) C++ signature :
setData((object)arg1, (object)data) Change the associated data container
setJacobian((object)arg1, (object)J) C++ signature :
setMesh((object)arg1, (object)mesh [, …) Set new mesh to the forward operator, optionally hold region parameter for the new mesh (i.e.
setMultiThreadJacobian((object)arg1, …) Set number of threads used for brute force Jacobian generation.
setRegionManager((object)arg1, (object)reg) C++ signature :
setStartModel((object)arg1, (object)startModel) C++ signature :
setThreadCount((object)arg1, (object)nThreads) Set the maximum number of allowed threads for MT calculation.
setVerbose((object)arg1, (object)verbose) Set verbose state.
solution((object)arg1) C++ signature :
startModel((object)arg1) C++ signature :
threadCount((object)arg1) Return the maximum number of allowed threads for MT calculation
verbose((object)arg1) Get verbose state.
createJacobian_mt  
responses  
createJacobian(model)[source]

create jacobian matrix using unchanged DC jacobian and m model

response(m)[source]

return forward response as function of chargeability model

pybert.tdip.mkallpdf module

pybert.tdip.tdipdata module

class pybert.tdip.tdipdata.TDIPdata(filename=None, **kwargs)[source]

Bases: object

Class managing time-domain induced polarisation (TDIP) field data.

Methods

ensureRhoa() Make sure apparent resistivity is present in file.
filter([tmin, tmax, kmax, electrode, …]) Filter data with respect to frequencies and geometric factor.
generatePDF([rdict, mdict]) Generate a multi-page pdf file with all data as pseudosections.
individualInversion(**kwargs) Carry out individual inversion for spectral chargeability.
invertMa([nr, ma]) Invert for chargeability.
invertRhoa(**kwargs) Invert apparent resistivity values.
load(filename, **kwargs) Try loading all supported file types.
save([filename]) Save all data in some (yet-to-be-defined or -decided) format.
showDecay([nr, ab, mn, verbose]) Show decay curves.
showMa([nr]) Show apparent chargeability.
showRhoa(**kwargs) Show apparent resistivity.
ensureRhoa()[source]

Make sure apparent resistivity is present in file.

filter(tmin=0, tmax=1000000000.0, kmax=1000000.0, electrode=None, forward=False, a=None, b=None, m=None, n=None, ab=None, mn=None, corrSID=1, nr=[])[source]

Filter data with respect to frequencies and geometric factor.

Parameters:
tmin : double

minimum frequency

tmax : double

maximum frequency

kmax : double

maximum (absolute) geometric factor

electrode : int

electrode to be removed completely

a/b/m/n : int

delete data with specific current or potential dipole lengths

ab/mn : int

delete data with specific current or potential dipole lengths

corrSID: int [1]

correct sensor index (like in data files)
generatePDF(rdict=None, mdict=None, **kwargs)[source]

Generate a multi-page pdf file with all data as pseudosections.

individualInversion(**kwargs)[source]

Carry out individual inversion for spectral chargeability.

invertMa(nr=0, ma=None, **kwargs)[source]

Invert for chargeability.

invertRhoa(**kwargs)[source]

Invert apparent resistivity values.

load(filename, **kwargs)[source]

Try loading all supported file types.

save(filename=None)[source]

Save all data in some (yet-to-be-defined or -decided) format.

showDecay(nr=[], ab=None, mn=None, verbose=True)[source]

Show decay curves.

showMa(nr=0, **kwargs)[source]

Show apparent chargeability.

showRhoa(**kwargs)[source]

Show apparent resistivity.

pybert.tdip.tdipdata.importTDIPdata(filename, verbose=False)[source]

Read in TDIP data.

pybert.tdip.tdipdata1 module