pygimli

name	pygimli
description	Multi-method geophysical modelling and inversion framework. Use when Claude needs to: (1) Perform electrical resistivity tomography (ERT) inversion, (2) Run seismic refraction tomography (SRT), (3) Model induced polarization (IP) data, (4) Simulate ground penetrating radar (GPR), (5) Create finite element meshes for geophysical problems, (6) Perform joint inversions of multiple datasets, (7) Forward model geophysical responses, (8) Analyze time-lapse monitoring data.
version	1.0.0
author	Geoscience Skills
license	MIT
tags	["Geophysical Inversion","ERT","Seismic Refraction","IP","Near-Surface","FEM"]
dependencies	["pygimli>=1.4.0","numpy","matplotlib"]
complements	["simpeg","verde","pyvista"]
workflow_role	modelling

pyGIMLi - Geophysical Inversion

Quick Reference

import pygimli as pg
from pygimli.physics import ert, srt

# Load ERT data
data = ert.load("survey.ohm")

# Invert
mgr = ert.ERTManager(data)
model = mgr.invert(lam=20, verbose=True)

# View result
mgr.showResult()

Key Classes

Class	Purpose
pg.Mesh	Finite element meshes
pg.DataContainer	Survey data and geometry
pg.Inversion	Base inversion framework
ert.ERTManager	ERT processing and inversion
srt.SRTManager	Seismic refraction inversion

Essential Operations

Load and View ERT Data

import pygimli as pg
from pygimli.physics import ert

data = ert.load("survey.ohm")
print(f"Measurements: {data.size()}")
ert.showData(data)  # Pseudosection

ERT Inversion

from pygimli.physics import ert

mgr = ert.ERTManager(data)
model = mgr.invert(
    lam=20,          # Regularization
    verbose=True
)
mgr.showResult()
resistivity = mgr.model

Seismic Refraction

from pygimli.physics import srt

data = srt.load("traveltimes.sgt")
mgr = srt.SRTManager(data)
model = mgr.invert(lam=, zWeight=)
mgr.showResult()