PyHydroGeophysX.analysis package#

Submodules#

PyHydroGeophysX.analysis.sensitivity module#

Sensitivity and resolution analysis utilities.

PyHydroGeophysX.analysis.sensitivity.compute_cumulative_sensitivity(J)[source]#: Compute cumulative sensitivity as absolute column sums of the Jacobian.

PyHydroGeophysX.analysis.sensitivity.compute_depth_of_investigation(inv_class, data: Any, mesh: Any, scale_low: float = 0.8, scale_high: float = 1.2) → Tuple[ndarray, Dict[str, ndarray]][source]#

Estimate DOI using two inversions with different reference/initial models.

This follows the Oldenburg-Li style idea of quantifying model sensitivity to reference-model choice.

PyHydroGeophysX.analysis.sensitivity.compute_resolution_matrix(J, Wd, Wm, lam)[source]#

Compute model resolution matrix:

R = (J^T Wd Wd J + lam Wm^T Wm)^(-1) J^T Wd Wd J

PyHydroGeophysX.analysis.sensitivity.plot_sensitivity_map(sensitivity, mesh, ax=None)[source]#: Plot a sensitivity map on a mesh (or as a simple line plot fallback).

Module contents#

Sensitivity and resolution analysis utilities.

PyHydroGeophysX.analysis.compute_cumulative_sensitivity(J)[source]#: Compute cumulative sensitivity as absolute column sums of the Jacobian.

PyHydroGeophysX.analysis.compute_depth_of_investigation(inv_class, data: Any, mesh: Any, scale_low: float = 0.8, scale_high: float = 1.2) → Tuple[ndarray, Dict[str, ndarray]][source]#

Estimate DOI using two inversions with different reference/initial models.

This follows the Oldenburg-Li style idea of quantifying model sensitivity to reference-model choice.

PyHydroGeophysX.analysis.compute_resolution_matrix(J, Wd, Wm, lam)[source]#

Compute model resolution matrix:

R = (J^T Wd Wd J + lam Wm^T Wm)^(-1) J^T Wd Wd J

PyHydroGeophysX.analysis.plot_sensitivity_map(sensitivity, mesh, ax=None)[source]#: Plot a sensitivity map on a mesh (or as a simple line plot fallback).