MORSE-PI -- Flexible and artefact-free image reconstruction for structural and functional QSM and other phase-critical imaging applications
Pith reviewed 2026-06-26 12:37 UTC · model grok-4.3
The pith
MORSE-PI generates high-SNR phase images free of singularities and artefacts for QSM by constructing a virtual reference coil.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
MORSE-PI produces artefact-free and singularity-free phase images with high SNR for both structural and functional brain imaging by using the VRC to correct phase offsets in coil sensitivity estimates, outperforming GRAPPA-ASPIRE in robustness and reproducibility at 3T and 7T, and unlike ESPIRiT or adaptive coil combination methods, avoiding singularities.
What carries the argument
The Virtual Reference Coil (VRC), a linear combination of coil sensitivity maps with correlations enhanced by the noise covariance matrix, which ensures robust signal support across the brain and corrects phase offsets without introducing singularities.
Load-bearing premise
Constructing the virtual reference coil as a linear combination of coil sensitivity maps with noise covariance enhancement will reliably ensure complete brain coverage and correct phase offsets without adding new artefacts or biases.
What would settle it
Finding phase singularities, noise amplification, or reduced reproducibility when applying MORSE-PI to the same multi-coil datasets used to benchmark GRAPPA with ASPIRE would disprove the central performance claims.
read the original abstract
Phase imaging applications such as QSM are highly sensitive to noise amplifications, phase singularities, and other artefacts, particularly in challenging scenarios such as ultra-high field (7T), under-sampled or single-echo acquisitions. We present a novel image reconstruction method, MORSE-PI, designed to produce high-SNR, artefact-free, and singularity-free phase images for both structural and functional phase-based brain imaging. MORSE-PI extends our previous approach, MORSE, by introducing a Virtual Reference Coil (VRC). The VRC is constructed as a linear combination of coil sensitivity maps, with correlations enhanced between coil elements using the noise covariance matrix. Such a VRC ensures robust signal support across the entire brain and is used to correct phase offsets in the MORSE-derived coil sensitivity estimates, resulting in artefact-free, high SNR phase. Compared to GRAPPA with ASPIRE phase correction, MORSE-PI demonstrates greater robustness to artefacts such as noise amplification and aliasing, and shows improved reproducibility in structural imaging at both 3T and 7T. Unlike ESPIRiT and GRAPPA combined with adaptive coil combination methods, MORSE-PI yields singularity-free phase maps. MORSE-PI enables high-SNR reconstructions even for the most challenging scenarios, such as single-echo EPI at 7T. Its efficient, containerised implementation using the Gadgetron framework supports deployment on the MRI scanner console during measurements. MORSE-PI offers a flexible and computationally efficient solution for generating high-SNR, artefact- and singularity-free phase images in both single- and multi-echo GRE and EPI acquisitions. This makes it particularly well-suited for structural and functional QSM, as well as other phase-based MRI applications. Its robustness and rapid computational time facilitate efficient deployment on scanners across field strengths.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript presents MORSE-PI, an extension of the prior MORSE method for MRI phase reconstruction in QSM and related applications. It introduces a Virtual Reference Coil (VRC) formed as a linear combination of coil sensitivity maps whose inter-element correlations are enhanced by the noise covariance matrix; this VRC is used to correct phase offsets in the MORSE-derived sensitivities. The central claims are that MORSE-PI yields high-SNR, artefact-free and singularity-free phase maps, exhibits greater robustness to noise amplification and aliasing than GRAPPA+ASPIRE, produces singularity-free maps unlike ESPIRiT or adaptive coil combination, and demonstrates improved reproducibility in structural imaging at both 3 T and 7 T, including single-echo EPI at 7 T. The implementation is containerised in Gadgetron for on-scanner use.
Significance. If the performance claims are supported by quantitative metrics and controlled validation in the full manuscript, MORSE-PI could offer a practically useful advance for phase-critical imaging at high fields by reducing common artefacts in under-sampled or single-echo acquisitions while remaining computationally efficient and scanner-deployable.
major comments (3)
- [Abstract] Abstract: the performance advantages ('greater robustness to artefacts such as noise amplification and aliasing', 'improved reproducibility', 'singularity-free phase maps') are asserted without any quantitative metrics, error bars, statistical tests, subject numbers, or data-exclusion criteria. This absence is load-bearing because the central claim is one of superiority over GRAPPA+ASPIRE, ESPIRiT and adaptive coil combination.
- [Abstract] Abstract (VRC paragraph): the VRC is described only declaratively as 'a linear combination of coil sensitivity maps, with correlations enhanced between coil elements using the noise covariance matrix' that 'ensures robust signal support across the entire brain' and corrects phase offsets. No explicit equations, matrix definitions, or derivation showing singularity elimination or bias-free behaviour are supplied; this construction is the load-bearing technical step for the artefact-free claim.
- [Abstract] Abstract: no validation is described for the key assumption that the covariance-weighted linear combination will maintain full-brain support in low-signal regions or avoid introducing spatially varying phase biases; without such controlled tests the singularity-free and artefact-free assertions cannot be assessed.
minor comments (1)
- [Abstract] The abstract is lengthy and repeats the list of advantages; condensing the claims while retaining the technical description of the VRC would improve readability.
Simulated Author's Rebuttal
We thank the referee for the constructive feedback on the abstract. We address each major comment below and will revise the abstract to better support the central claims with quantitative information from the full manuscript where feasible.
read point-by-point responses
-
Referee: [Abstract] Abstract: the performance advantages ('greater robustness to artefacts such as noise amplification and aliasing', 'improved reproducibility', 'singularity-free phase maps') are asserted without any quantitative metrics, error bars, statistical tests, subject numbers, or data-exclusion criteria. This absence is load-bearing because the central claim is one of superiority over GRAPPA+ASPIRE, ESPIRiT and adaptive coil combination.
Authors: We agree that the abstract would benefit from highlighting key quantitative results to support the superiority claims. The full manuscript includes such metrics from in vivo experiments at 3T and 7T (including single-echo EPI), with comparisons to the referenced methods. We will revise the abstract to incorporate representative quantitative metrics, subject numbers, and reference to the validation approach. revision: yes
-
Referee: [Abstract] Abstract (VRC paragraph): the VRC is described only declaratively as 'a linear combination of coil sensitivity maps, with correlations enhanced between coil elements using the noise covariance matrix' that 'ensures robust signal support across the entire brain' and corrects phase offsets. No explicit equations, matrix definitions, or derivation showing singularity elimination or bias-free behaviour are supplied; this construction is the load-bearing technical step for the artefact-free claim.
Authors: Abstract length constraints limit the inclusion of full equations or derivations, which are provided in the Methods section. We will revise the VRC description to include a more explicit formulation of the linear combination weighted by the noise covariance matrix and note its role in phase offset correction and singularity avoidance. revision: partial
-
Referee: [Abstract] Abstract: no validation is described for the key assumption that the covariance-weighted linear combination will maintain full-brain support in low-signal regions or avoid introducing spatially varying phase biases; without such controlled tests the singularity-free and artefact-free assertions cannot be assessed.
Authors: The full manuscript validates the VRC assumption via controlled in vivo tests at 3T and 7T demonstrating full-brain support without introduced biases or singularities. We will revise the abstract to briefly reference this validation approach. revision: yes
Circularity Check
No significant circularity detected in derivation chain
full rationale
The paper presents MORSE-PI as an extension of prior MORSE work via an explicit linear-algebra construction of the Virtual Reference Coil as a linear combination of measured coil sensitivity maps, with correlations enhanced by the noise covariance matrix. This construction is described declaratively in the abstract without any equations that reduce the resulting phase maps, singularity elimination, or performance metrics back to fitted parameters or self-referential definitions drawn from the same input data. The self-reference to the authors' prior MORSE method is not load-bearing for any circular reduction, as the new VRC step is presented as an independent linear combination rather than a tautological renaming or prediction of its own inputs. No self-definitional, fitted-input, or ansatz-smuggling patterns appear in the provided text.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption A linear combination of coil sensitivity maps weighted by the noise covariance matrix produces a virtual reference coil with robust signal support across the brain.
invented entities (1)
-
Virtual Reference Coil (VRC)
no independent evidence
Reference graph
Works this paper leans on
-
[1]
High-Resolution MR Venography at 3.0 Tesla
Reichenbach JR, Barth M, Haacke EM, Klarhöfer M, Kaiser WA, Moser E. High-Resolution MR Venography at 3.0 Tesla. J Comput Assist Tomogr. 2000;24(6):949
2000
-
[2]
Imaging iron stores in the brain using magnetic resonance imaging
Haacke EM, Cheng NYC, House MJ, et al. Imaging iron stores in the brain using magnetic resonance imaging. Magn Reson Imaging. 2005;23(1):1-25. doi:10.1016/j.mri.2004.10.001
-
[3]
Lee J, Shmueli K, Kang BT, et al. The contribution of myelin to magnetic susceptibility-weighted contrasts in high-field MRI of the brain. NeuroImage. 2012;59(4):3967-3975. doi:10.1016/j.neuroimage.2011.10.076
-
[4]
Brain magnetic resonance imaging with contrast dependent on blood oxygenation
Ogawa S, Lee TM, Kay AR, Tank DW. Brain magnetic resonance imaging with contrast dependent on blood oxygenation. Proc Natl Acad Sci. 1990;87(24):9868-9872
1990
-
[5]
Susceptibility weighted imaging (SWI)
Haacke EM, Xu Y, Cheng YCN, Reichenbach JR. Susceptibility weighted imaging (SWI). Magn Reson Med. 2004;52(3):612-618. doi:10.1002/mrm.20198
-
[6]
Magnetic susceptibility mapping of brain tissue in vivo using MRI phase data
Shmueli K, de Zwart JA, van Gelderen P, Li TQ, Dodd SJ, Duyn JH. Magnetic susceptibility mapping of brain tissue in vivo using MRI phase data. Magn Reson Med. 2009;62(6):1510-1522. doi:10.1002/mrm.22135
-
[7]
Functional quantitative susceptibility mapping (fQSM)
Balla DZ, Sanchez-Panchuelo RM, Wharton SJ, et al. Functional quantitative susceptibility mapping (fQSM). NeuroImage. 2014;100:112-124. doi:10.1016/j.neuroimage.2014.06.011
-
[8]
Liu C. Susceptibility tensor imaging. Magn Reson Med. 2010;63(6):1471-1477. doi:10.1002/mrm.22482
-
[9]
χ-separation: Magnetic susceptibility source separation toward iron and myelin mapping in the brain
Shin HG, Lee J, Yun YH, et al. χ-separation: Magnetic susceptibility source separation toward iron and myelin mapping in the brain. NeuroImage. 2021;240:118371. doi:10.1016/j.neuroimage.2021.118371
-
[10]
Chen J, Gong NJ, Chaim KT, Otaduy MCG, Liu C. Decompose quantitative susceptibility mapping (QSM) to sub-voxel diamagnetic and paramagnetic components based on gradient-echo MRI data. NeuroImage. 2021;242:118477. doi:10.1016/j.neuroimage.2021.118477
-
[11]
Robinson SD, Bredies K, Khabipova D, Dymerska B, Marques JP, Schweser F. An illustrated comparison of processing methods for MR phase imaging and QSM: combining array coil signals and phase unwrapping. NMR Biomed. 2017;30(4):e3601. doi:10.1002/nbm.3601
-
[12]
Susceptibility-weighted imaging: current status and future directions
Liu S, Buch S, Chen Y, et al. Susceptibility-weighted imaging: current status and future directions. NMR Biomed. 2017;30(4):e3552. doi:10.1002/nbm.3552
-
[13]
Bilgic B, Costagli M, Chan KS, et al. Recommended implementation of quantitative susceptibility mapping for clinical research in the brain: A consensus of the ISMRM electro-magnetic tissue properties study group. Magn Reson Med. 2024;91(5):1834-1862. doi:10.1002/mrm.30006
-
[14]
Eckstein K, Dymerska B, Bachrata B, et al. Computationally Efficient Combination of Multi-channel Phase Data From Multi-echo Acquisitions (ASPIRE). Magn Reson Med. 2018;79(6):2996-3006. doi:10.1002/mrm.26963
-
[15]
Online image reconstruction via Multiple Orthogonal Reference Sensitivity Encoding (MORSE)
Josephs O, Dymerska B, Graedel N, Balbastre Y, Corbin N, Callaghan MF. Online image reconstruction via Multiple Orthogonal Reference Sensitivity Encoding (MORSE). Magn Reson Mater Phys Biol Med. In Press 2026
2026
-
[16]
SENSE: sensitivity encoding for fast MRI
Pruessmann KP, Weiger M, Scheidegger MB, Boesiger P. SENSE: sensitivity encoding for fast MRI. Magn Reson Med. 1999;42(5):952-962
1999
-
[17]
MORSE: Multiple Orthogonal Reference Sensitivity Encoding
Josephs O, Dymerska B, Graedel NN, Balbastre Y, Corbin N, Callaghan MF. MORSE: Multiple Orthogonal Reference Sensitivity Encoding. arXiv. Preprint posted online October 10, 2025:arXiv:2510.09098. doi:10.48550/arXiv.2510.09098
-
[18]
ESPIRiT—an eigenvalue approach to autocalibrating parallel MRI: Where SENSE meets GRAPPA
Uecker M, Lai P, Murphy MJ, et al. ESPIRiT—an eigenvalue approach to autocalibrating parallel MRI: Where SENSE meets GRAPPA. Magn Reson Med. 2014;71(3):990-1001. doi:10.1002/mrm.24751
-
[19]
Gadgetron: An open source framework for medical image reconstruction
Hansen MS, Sørensen TS. Gadgetron: An open source framework for medical image reconstruction. Magn Reson Med. 2013;69(6):1768-1776. doi:10.1002/mrm.24389
-
[20]
Phase reconstruction from multiple coil data using a virtual reference coil
Parker DL, Payne A, Todd N, Hadley JR. Phase reconstruction from multiple coil data using a virtual reference coil. Magn Reson Med. 2014;72(2):563-569. doi:10.1002/mrm.24932
-
[21]
Design of a SENSE-optimized high- sensitivity MRI receive coil for brain imaging
de Zwart JA, Ledden PJ, Kellman P, van Gelderen P, Duyn JH. Design of a SENSE-optimized high- sensitivity MRI receive coil for brain imaging. Magn Reson Med. 2002;47(6):1218-1227. doi:10.1002/mrm.10169
-
[22]
Hammond KE, Lupo JM, Xu D, et al. Development of a robust method for generating 7.0 T multichannel phase images of the brain with application to normal volunteers and patients with neurological diseases. NeuroImage. 2008;39(4):1682-1692. doi:10.1016/j.neuroimage.2007.10.037
-
[23]
Published online March 3, 2026
fil-physics/gadgetron-matlab. Published online March 3, 2026. Accessed March 23, 2026. https://github.com/fil-physics/gadgetron-matlab
2026
-
[24]
Published online March 10, 2026
gadgetron/gadgetron. Published online March 10, 2026. Accessed March 23, 2026. https://github.com/gadgetron/gadgetron
2026
-
[25]
Mitigating susceptibility-induced distortions in high-resolution 3DEPI fMRI at 7T
Malekian V, Graedel NN, Hickling A, et al. Mitigating susceptibility-induced distortions in high-resolution 3DEPI fMRI at 7T. NeuroImage. 2023;279:120294. doi:10.1016/j.neuroimage.2023.120294
-
[26]
Generalized autocalibrating partially parallel acquisitions ( GRAPPA )
Griswold MA, Jakob PM, Heidemann RM, et al. Generalized autocalibrating partially parallel acquisitions (GRAPPA). Magn Reson Med. 2002;47(6):1202-1210. doi:10.1002/mrm.10171
-
[27]
Adaptive reconstruction of phased array MR imagery
Walsh DO, Gmitro AF, Marcellin MW. Adaptive reconstruction of phased array MR imagery. Magn Reson Med. 2000;43(5):682-690. doi:10.1002/(SICI)1522-2594(200005)43:5<682::AID-MRM10>3.0.CO;2-G
-
[28]
Release V3.0 - MPM_QSM for MORSE-Pi recon without phase singularities (containers >=V13.4) · fil- physics/MPM_QSM. GitHub. Accessed March 23, 2026. https://github.com/fil- physics/MPM_QSM/releases/tag/v3.0
2026
-
[29]
Eckstein K, Bachrata B, Hangel G, et al. Improved susceptibility weighted imaging at ultra-high field using bipolar multi-echo acquisition and optimized image processing: CLEAR-SWI. NeuroImage. 2021;237:118175. doi:10.1016/j.neuroimage.2021.118175
-
[30]
Phase unwrapping with a rapid opensource minimum spanning tree algorithm (ROMEO)
Dymerska B, Eckstein K, Bachrata B, et al. Phase unwrapping with a rapid opensource minimum spanning tree algorithm (ROMEO). Magn Reson Med. 2021;85(4):2294-2308. doi:10.1002/mrm.28563
-
[31]
SPM 25: open source neuroimaging analysis software
Tierney TM, Alexander NA, Avila NL, et al. SPM 25: open source neuroimaging analysis software. J Open Source Softw. 2025;10(110):8103. doi:10.21105/joss.08103
-
[32]
Ashburner J, Friston KJ. Unified segmentation. NeuroImage. 2005;26(3):839-851. doi:10.1016/j.neuroimage.2005.02.018
-
[33]
A novel background field removal method for MRI using projection onto dipole fields (PDF)
Liu T, Khalidov I, Rochefort L de, et al. A novel background field removal method for MRI using projection onto dipole fields (PDF). NMR Biomed. 2011;24(9):1129-1136. doi:10.1002/nbm.1670
-
[34]
Wei H, Dibb R, Zhou Y, et al. Streaking artifact reduction for quantitative susceptibility mapping of sources with large dynamic range. NMR Biomed. 2015;28(10):1294-1303. doi:10.1002/nbm.3383
-
[35]
SEPIA—Susceptibility mapping pipeline tool for phase images
Chan KS, Marques JP. SEPIA—Susceptibility mapping pipeline tool for phase images. NeuroImage. 2021;227:117611. doi:10.1016/j.neuroimage.2020.117611
-
[36]
Wu B, Li W, Avram AV, Gho SM, Liu C. Fast and tissue-optimized mapping of magnetic susceptibility and T2* with multi-echo and multi-shot spirals. NeuroImage. 2012;59(1):297-305. doi:10.1016/j.neuroimage.2011.07.019
-
[37]
Ultra-High Field Neuro MRI
Bloch KM, Guye M, Poser BA. Ultra-High Field Neuro MRI. Elsevier; 2023
2023
-
[38]
Nonlinear formulation of the magnetic field to source relationship for robust quantitative susceptibility mapping
Liu T, Wisnieff C, Lou M, Chen W, Spincemaille P, Wang Y. Nonlinear formulation of the magnetic field to source relationship for robust quantitative susceptibility mapping. Magn Reson Med. 2013;69(2):467-
2013
-
[39]
doi:10.1002/mrm.24272
-
[40]
Encoding and reconstruction in parallel MRI
Pruessmann KP. Encoding and reconstruction in parallel MRI. NMR Biomed. 2006;19(3):288-299. doi:10.1002/nbm.1042
-
[41]
Warrington O, Graedel NN, Callaghan MF, Kok P. Communication of perceptual predictions from the hippocampus to the deep layers of the parahippocampal cortex. Sci Adv. 2025;11(21):eads4970. doi:10.1126/sciadv.ads4970
-
[42]
Background field removal using spherical mean value filtering and Tikhonov regularization
Sun H, Wilman AH. Background field removal using spherical mean value filtering and Tikhonov regularization. Magn Reson Med. 2014;71(3):1151-1157. doi:10.1002/mrm.24765
-
[43]
Exploring linearity of deep neural network trained QSM: QSMnet+
Jung W, Yoon J, Ji S, et al. Exploring linearity of deep neural network trained QSM: QSMnet+. NeuroImage. 2020;211:116619. doi:10.1016/j.neuroimage.2020.116619
-
[44]
Moeller S, Pisharady PK, Ramanna S, et al. NOise reduction with DIstribution Corrected (NORDIC) PCA in dMRI with complex-valued parameter-free locally low-rank processing. NeuroImage. 2021;226:117539. doi:10.1016/j.neuroimage.2020.117539
-
[45]
ENLIVE: An Efficient Nonlinear Method for Calibrationless and Robust Parallel Imaging
Holme HCM, Rosenzweig S, Ong F, Wilke RN, Lustig M, Uecker M. ENLIVE: An Efficient Nonlinear Method for Calibrationless and Robust Parallel Imaging. Sci Rep. 2019;9(1):3034. doi:10.1038/s41598-019-39888- 7
-
[46]
Low-Rank Modeling of Local k-Space Neighborhoods (LORAKS) for Constrained MRI
Haldar JP. Low-Rank Modeling of Local k-Space Neighborhoods (LORAKS) for Constrained MRI. IEEE Trans Med Imaging. 2014;33(3):668-681. doi:10.1109/TMI.2013.2293974
-
[47]
Pohmann R, Speck O, Scheffler K. Signal-to-noise ratio and MR tissue parameters in human brain imaging at 3, 7, and 9.4 tesla using current receive coil arrays. Magn Reson Med. 2016;75(2):801-809. doi:10.1002/mrm.25677
-
[48]
Robust 3D Bloch‐Siegert based B1+ mapping using multi‐echo general linear modeling
Corbin N, Acosta‐Cabronero J, Malik SJ, Callaghan MF. Robust 3D Bloch‐Siegert based B1+ mapping using multi‐echo general linear modeling. Magn Reson Med. 2019;82(6):2003-2015. doi:10.1002/mrm.27851
-
[49]
A method for the dynamic correction of B0- related distortions in single-echo EPI at 7T
Dymerska B, Poser BA, Barth M, Trattnig S, Robinson SD. A method for the dynamic correction of B0- related distortions in single-echo EPI at 7T. NeuroImage. 2018;168:321-331. doi:10.1016/j.neuroimage.2016.07.009
-
[50]
Robinson SD, Bachrata B, Eckstein K, et al. Improved dynamic distortion correction for fMRI using single- echo EPI and a readout-reversed first image (REFILL). Hum Brain Mapp. 2023;44(15):5095-5112. doi:10.1002/hbm.26440
-
[51]
(ISMRM 2022) Why is 7T fMRI in the entorhinal cortex so difficult and what can we do about it? Accessed March 23, 2026
Graedel N, Warrington O, Dymerska B, et al. (ISMRM 2022) Why is 7T fMRI in the entorhinal cortex so difficult and what can we do about it? Accessed March 23, 2026. https://archive.ismrm.org/2022/1096.html
2022
-
[52]
Li J, Lin H, Liu T, et al. Quantitative susceptibility mapping (QSM) minimizes interference from cellular pathology in R2* estimation of liver iron concentration. J Magn Reson Imaging. 2018;48(4):1069-1079. doi:10.1002/jmri.26019
-
[53]
Wei H, Dibb R, Decker K, et al. Investigating magnetic susceptibility of human knee joint at 7 tesla. Magn Reson Med. Published online January 1, 2017:n/a-n/a. doi:10.1002/mrm.26596 Supplementary material Figure S 1 Noise covariance matrices illustrating the effects of whitening or increasing cross -channel correlation in a 64 -channel coil at 3T . Ψwhite...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.