Published August 1, 2026 | Version 2

QSM-CI method: WaveSep (v2)

Description

WaveSep — a wavelet-based iterative approach for susceptibility source separation. It splits net susceptibility (χ_total/QSM) into paramagnetic (χ+, iron) and diamagnetic (χ−, myelin·calcium) sources using an R2' map to break the para/dia degeneracy, under two voxel-wise data-fidelity terms (χ+ + χ− ≈ χ_total, χ+ − χ− ≈ R2'/Dr) and a wavelet-domain L1 sparsity prior solved by proximal gradient (soft-thresholding, db4). Pure Python on CPU — no network, no pretrained weights, no GPU.

QSM-CI reconstruction method wavesep. Browse and run it at https://qsmxt.github.io/QSM-CI/submission.html?method=wavesep.

Files

wavesep.zip

Files (6.2 kB)

Name Size Download all
md5:ec290e6aff736d397ba77d56a363ecf4
6.2 kB Preview Download

Additional details