---
title: Quantitative Susceptibility Mapping (QSM)
url: https://www.emergentmind.com/topics/quantitative-susceptibility-mapping-qsm
type: topic
---

# Quantitative Susceptibility Mapping (QSM)

Quantitative Susceptibility Mapping (QSM) is a magnetic resonance imaging (MRI) phase post-processing technique that quantifies the spatial distribution of tissue magnetic susceptibility, providing critical biomarkers for iron, calcium, and myelin content. QSM has demonstrated significant potential in neurological disease research, sub-millimeter neuroimaging, and the quantification of deep grey matter changes. The primary challenge in QSM is the ill-posedness of the dipole inversion, which can lead to noise amplification and artifacts if not rigorously addressed [2406.12300][2106.00525][2501.17158].

## 1. Mathematical Model and Ill-Posed Dipole Inversion

The foundational forward model in QSM relates the measured local field map $m(r)$ to the underlying tissue susceptibility distribution $\chi(r)$ via convolution with the dipole kernel $d(r)$, with added noise $\epsilon(r)$:
\[
m(r) = (d * \chi)(r) + \epsilon(r)
\]
In the Fourier domain, this becomes:
\[
M(k) = D(k) \cdot \Chi(k) + E(k)
\]
where $D(k) = \frac{1}{3} - \frac{k_z^2}{k_x^2 + k_y^2 + k_z^2}$. The transfer function $D(k)$ is zero on a conical surface ("magic angle"), making the inversion ill-posed; direct inversion amplifies noise in the null-space, producing characteristic streaking artifacts [2406.12300][2106.00525][2501.17158].

To remedy this, various regularized inversion frameworks are employed, including Tikhonov regularization, sparsity-promoting variational models, and modern deep learning architectures. The standard variational problem takes the form:
\[
\hat\chi = \arg\min_\chi \| F^{-1}[ D \cdot F(\chi)] - m \|_2^2 + \lambda R(\chi)
\]
where $R(\chi)$ can encode sparsity, morphology, or other priors [2106.00525][2501.17158].

## 2. Deep Neural Network Approaches and Architectures

Recent advances leverage 3D convolutional neural networks (CNNs) and physics-informed models for QSM. IR2QSM [2406.12300] exemplifies modern architectures:
- **IR2QSM Cascade**: Four U-net variants (IR2U-nets) are cascaded to iteratively refine susceptibility estimates. Each stage incorporates *Reverse Concatenation* (RC), fusing semantic features from the prior decoder into the current encoder, and a *Middle Recurrent Module* (RM), inspired by gated RNNs, for capturing long-range dependencies critical to resolving the non-local nature of the dipole kernel.
- **Encoder-Decoder Blocks**: Each U-net block consists of two $3^3$ convolutions, BatchNorm, ReLU, and (up/down)sampling as appropriate.
- **Output Aggregation**: Outputs from all cascades are concatenated and fused via a $1 \times 1 \times 1$ convolution, yielding the final susceptibility volume.

The training loss is a weighted sum of MSE at each stage with empirically-decayed weights plus a terminal MSE. Regularization via noise augmentation and dropout enhances in vivo generalization.

xQSM [2106.00525] employs a U-net with octave convolutions, splitting features into high- and low-frequency streams for efficient and context-rich processing. Denoising and patch-wise training further improve robustness to field-of-view truncation and noise.

## 3. Pipeline Steps: Phase to Susceptibility

QSM reconstruction comprises several processing stages, each with algorithmic alternatives:

- **Phase Unwrapping**: Unwraps the 2π ambiguity in multi-echo phase, with 3D Laplacian or path-following algorithms.
- **Background Field Removal (BFR)**: Separates local tissue fields from background phase contributions, typically via RESHARP (spherical mean value with Tikhonov regularization), SHARP, V-SHARP, PDF, or LBV. The spatial support and mask erosion directly impact reconstruction biases, especially in reduced field-of-view or slab-coverages [2106.00525][2501.17158].
- **Dipole Inversion**: Core step, solving for $\chi$ from the preprocessed field, using inverse filtering (e.g., TKD), iterative regularized solvers (LSQR, MEDI, HEIDI), or deep learning models (U-net, xQSM, IR2QSM) [2406.12300][2106.00525][2501.17158].
- **Referencing**: Final QSM is referenced to white matter, whole brain, or CSF mean to resolve the non-uniqueness up to a constant [$\chi_{\text{ref}}(r) = \chi(r) - \langle \chi \rangle_{R_{ref}}$].

Table: Impact of Processing Choices on DGM Susceptibility Error [2501.17158][2106.00525]

| Step                 | Algorithm           | Impact on DGM Accuracy      |
|----------------------|--------------------|----------------------------|
| Background Field     | RESHARP (BFR)      | Minimal bias, stable error  |
|                      | PDF, LBV           | More artifacts near edges   |
| Dipole Inversion     | IR2QSM, xQSM       | Lowest error, robust to FOV |
|                      | iLSQR, MEDI        | Higher error, more artifacts|
| Referencing          | WM, WB             | Lower CV, higher sensitivity|
|                      | CSF                | Increased variability       |

## 4. Quantitative Performance and Clinical Sensitivity

In both simulation and in vivo studies, IR2QSM and xQSM deliver leading performance. On COSMOS-based simulated brains [2406.12300]:
- IR2QSM achieves NRMSE 27.59% (versus 35–56% for competing methods), with regression slopes in deep grey matter ROIs approximating 0.94–1.01, outperforming iLSQR, MEDI, and other network baselines in both accuracy and artifact suppression.
- In real-world protocols focused on deep grey matter, xQSM achieves <5% error in globus pallidus with only 48 mm axial coverage, compared to >40% underestimation with classical approaches unless much larger SOVs (>112 mm) are acquired [2106.00525].

Systematic clinical pipeline evaluations using 378 variants [2501.17158] show that the choice of BFR, inversion, and referencing can affect reproducibility error and sensitivity by up to two orders of magnitude. RESHARP BFR with LSQR, HEIDI, or AMP-PE inversion, and white matter/whole brain referencing, provide optimal trade-offs between sensitivity to physiological changes and reproducibility.

## 5. Limitations and Algorithmic Trade-Offs

- **Ill-posedness**: All dipole inversions are fundamentally ill-posed; regularization or data-driven priors are mandatory.
- **Finite Field-of-View**: Reduced spatial coverage degrades background field removal efficacy and can induce DGM underestimation unless specifically addressed by network design and training (xQSM) [2106.00525].
- **Computational Complexity**: Networks like IR2QSM incur elevated computational load (∼16.87 GFLOPS/patch) and increased inference time compared to single-shot U-nets, but with substantial gains in error reduction (∼25% NRMSE decrease from T=1→4 cascades) [2406.12300].
- **Generalization**: Both IR2QSM and xQSM were trained primarily on 1 mm isotropic brain MRI. Their generalizability to other resolutions, anatomies, or pathologies remains underexplored.
- **Sensitivity to Algorithmic Details**: Sensitivity and reproducibility are strongly pipeline-dependent. In multi-center or longitudinal studies, standardization is critical [2501.17158].

## 6. Future Directions

Emergent directions in QSM research include:
- **Lightweight Networks**: Architectural simplifications of RC/RM modules for real-time and deployment constraints [2406.12300].
- **Multi-Resolution Training**: Networks that generalize across spatial resolutions and anatomical targets.
- **Joint BFR + Inversion**: Physics-informed or end-to-end models that subsume BFR, improving through-plane fidelity in high-resolution or truncated SOV data [2406.12300][2106.00525].
- **Uncertainty Quantification**: Bayesian networks, dropout-based confidence estimation, and robust harmonization strategies for clinical deployment [2501.17158].
- **Open Benchmarking**: Large-scale comparative studies across varied pathologies and acquisition protocols, with open-source reference pipelines and trained models.

## 7. Significance and Clinical Impact

QSM provides a direct in vivo measure of tissue microstructure relevant to neurodegeneration, demyelination, iron overload, and microbleeds. Advancements in robust inversion (IR2QSM, xQSM) facilitate accurate, time-efficient, and artifact-suppressed mapping even under clinical scan time or coverage constraints. Optimal pipeline customization—BFR method, inversion approach, referencing—is essential for high-sensitivity and reproducibility in translational research [2501.17158]. The field continues to evolve rapidly toward interpretable, scalable, and generalizable QSM solutions [2406.12300][2106.00525].

Source: https://www.emergentmind.com/topics/quantitative-susceptibility-mapping-qsm