cedalion.imagereco.solver

Solver for the image reconstruction problem.

Functions

pseudo_inverse_stacked(Adot[, alpha, Cmeas, ...])

Calculate the pseudo-inverse of a stacked sensitivity matrix.
cedalion.imagereco.solver.pseudo_inverse_stacked(
Adot: DataArray,
alpha: float = 0.01,
Cmeas: _SupportsArray[dtype[Any]] | _NestedSequence[_SupportsArray[dtype[Any]]] | bool | int | float | complex | str | bytes | _NestedSequence[bool | int | float | complex | str | bytes] | None = None,
alpha_spatial: float | None = None,
)[source]

Calculate the pseudo-inverse of a stacked sensitivity matrix.

Parameters:

  • Adot – Stacked matrix

  • alpha – Thikonov regularization parameter

  • Cmeas – Optional measurement regularization parameters. If specified this can be either a vector of size nchannel or a matrix of size nchannelxnchannel.

  • alpha_spatial – Optional spatial regularization parameter. Suggested default is 1e-3, or 1e-2 when spatial basis functions are used.

Returns:

Pseudo-inverse of the stacked matrix.

Return type:

xr.DataArray