Iterative Projection Onto Convex Sets for Quantitative Susceptibility Mapping

Overview

Journal Magn Reson Med

Publisher Wiley

Specialty Radiology

Date 2014 Mar 8

PMID 24604410

Citations 1

Authors

Weiran Deng

Fernando Boada

Benedikt A Poser

Claudiu Schirda

Victor Andrew Stenger

Affiliations

Soon will be listed here.

Abstract

Purpose: Quantitative susceptibility map (QSM) reconstruction is ill posed due to the zero values on the "magic angle cone" that make the maps prone to streaking artifacts. We propose projection onto convex sets (POCS) in the method of steepest descent (SD) for QSM reconstruction.

Methods: Two convex projections, an object-support projection in the image domain and a projection in k-space were used. QSM reconstruction using the proposed SD-POCS method was compared with SD and POCS alone as well as with truncated k-space division (TKD) for numerically simulated and 7 Tesla (T) human brain phase data.

Results: The QSM reconstruction error from noise-free simulated phase data using SD-POCS is at least two orders of magnitude lower than using SD, POCS, or TKD and has reduced streaking artifacts. Using the l1 -TV reconstructed susceptibility as a gold standard for 7T in vivo imaging, SD-POCS showed better image quality comparing to SD, POCS, or TKD from visual inspection.

Conclusion: POCS is an alternative method for regularization that can be used in an iterative minimization method such as SD for QSM reconstruction.

Citing Articles

The 2016 QSM Challenge: Lessons learned and considerations for a future challenge design.

Milovic C, Tejos C, Acosta-Cabronero J, Ozbay P, Schwesser F, Marques J Magn Reson Med. 2020; 84(3):1624-1637.

PMID: 32086836 PMC: 7526054. DOI: 10.1002/mrm.28185.

References

Bilgic B, Pfefferbaum A, Rohlfing T, Sullivan E, Adalsteinsson E . MRI estimates of brain iron concentration in normal aging using quantitative susceptibility mapping. Neuroimage. 2011; 59(3):2625-35. PMC: 3254708. DOI: 10.1016/j.neuroimage.2011.08.077. View

Wharton S, Bowtell R . Whole-brain susceptibility mapping at high field: a comparison of multiple- and single-orientation methods. Neuroimage. 2010; 53(2):515-25. DOI: 10.1016/j.neuroimage.2010.06.070. View

de Rochefort L, Liu T, Kressler B, Liu J, Spincemaille P, Lebon V . Quantitative susceptibility map reconstruction from MR phase data using bayesian regularization: validation and application to brain imaging. Magn Reson Med. 2009; 63(1):194-206. DOI: 10.1002/mrm.22187. View

Liu C . Susceptibility tensor imaging. Magn Reson Med. 2010; 63(6):1471-7. PMC: 2990786. DOI: 10.1002/mrm.22482. View

Abdul-Rahman H, Gdeisat M, Burton D, Lalor M, Lilley F, Moore C . Fast and robust three-dimensional best path phase unwrapping algorithm. Appl Opt. 2007; 46(26):6623-35. DOI: 10.1364/ao.46.006623. View

Liu C, Li W, Wu B, Jiang Y, Johnson G . 3D fiber tractography with susceptibility tensor imaging. Neuroimage. 2011; 59(2):1290-8. PMC: 3235503. DOI: 10.1016/j.neuroimage.2011.07.096. View

Liu C, Li W, Johnson G, Wu B . High-field (9.4 T) MRI of brain dysmyelination by quantitative mapping of magnetic susceptibility. Neuroimage. 2011; 56(3):930-8. PMC: 3085608. DOI: 10.1016/j.neuroimage.2011.02.024. View

Li L . Magnetic susceptibility quantification for arbitrarily shaped objects in inhomogeneous fields. Magn Reson Med. 2001; 46(5):907-16. DOI: 10.1002/mrm.1276. View

Liu T, Khalidov I, de Rochefort L, Spincemaille P, Liu J, Tsiouris A . A novel background field removal method for MRI using projection onto dipole fields (PDF). NMR Biomed. 2011; 24(9):1129-36. PMC: 3628923. DOI: 10.1002/nbm.1670. View

10.

Raj A, Zhang H, Prince M, Wang Y, Zabih R . Automatic algorithm for correcting motion artifacts in time-resolved two-dimensional magnetic resonance angiography using convex projections. Magn Reson Med. 2006; 55(3):649-58. DOI: 10.1002/mrm.20806. View

11.

Schweser F, Deistung A, Lehr B, Reichenbach J . Quantitative imaging of intrinsic magnetic tissue properties using MRI signal phase: an approach to in vivo brain iron metabolism?. Neuroimage. 2010; 54(4):2789-807. DOI: 10.1016/j.neuroimage.2010.10.070. View

12.

Sezan M, Stark H . Image restoration by the method of convex projections: part 2 applications and numerical results. IEEE Trans Med Imaging. 1982; 1(2):95-101. DOI: 10.1109/TMI.1982.4307556. View

13.

HALLGREN B, SOURANDER P . The non-haemin iron in the cerebral cortex in Alzheimer's disease. J Neurochem. 1960; 5:307-10. DOI: 10.1111/j.1471-4159.1960.tb13369.x. View

14.

Liu T, Spincemaille P, de Rochefort L, Kressler B, Wang Y . Calculation of susceptibility through multiple orientation sampling (COSMOS): a method for conditioning the inverse problem from measured magnetic field map to susceptibility source image in MRI. Magn Reson Med. 2008; 61(1):196-204. DOI: 10.1002/mrm.21828. View

15.

Smith S, Jenkinson M, Woolrich M, Beckmann C, Behrens T, Johansen-Berg H . Advances in functional and structural MR image analysis and implementation as FSL. Neuroimage. 2004; 23 Suppl 1:S208-19. DOI: 10.1016/j.neuroimage.2004.07.051. View

16.

Youla D, Webb H . Image restoration by the method of convex projections: part 1 theory. IEEE Trans Med Imaging. 1982; 1(2):81-94. DOI: 10.1109/TMI.1982.4307555. View

17.

Samsonov A, Kholmovski E, Parker D, Johnson C . POCSENSE: POCS-based reconstruction for sensitivity encoded magnetic resonance imaging. Magn Reson Med. 2004; 52(6):1397-406. DOI: 10.1002/mrm.20285. View

18.

Robinson S, Grabner G, Witoszynskyj S, Trattnig S . Combining phase images from multi-channel RF coils using 3D phase offset maps derived from a dual-echo scan. Magn Reson Med. 2011; 65(6):1638-48. DOI: 10.1002/mrm.22753. View

19.

Samsonov A, Velikina J, Jung Y, Kholmovski E, Johnson C, Block W . POCS-enhanced correction of motion artifacts in parallel MRI. Magn Reson Med. 2010; 63(4):1104-10. PMC: 2852198. DOI: 10.1002/mrm.22254. View

20.

Shmueli K, de Zwart J, van Gelderen P, Li T, Dodd S, Duyn J . Magnetic susceptibility mapping of brain tissue in vivo using MRI phase data. Magn Reson Med. 2009; 62(6):1510-22. PMC: 4275127. DOI: 10.1002/mrm.22135. View