Expand description
Tikhonov regularization for QSM
Tikhonov regularization adds an L2 penalty term to stabilize the inversion:
χ = argmin_x ||Dx - f||₂² + λ||Γx||₂²
This has a closed-form solution in k-space: χ̂ = D* · f̂ / (|D|² + λ|Γ|²)
Reference: Bilgic, B., et al. (2014). “Fast image reconstruction with L2-regularization.” Journal of Magnetic Resonance Imaging, 40(1):181-191. https://doi.org/10.1002/jmri.24365
Reference implementation: https://github.com/kamesy/QSM.jl
Structs§
- Tikhonov
Params - Tikhonov algorithm parameters
Enums§
- Regularization
- Regularization type for Tikhonov
Functions§
- tikhonov
- Tikhonov regularization for dipole inversion
- tikhonov_
default - Tikhonov with default parameters (identity regularization, λ=0.01; matches QSM.jl)