REVIEW 2 major objections 1 minor 2 references
An open-source pipeline using PyPulseq and BART delivers vendor-agnostic 4D phase contrast MRI for muscle velocity, displacement, and strain.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
An open-source vendor-agnostic 4D PC-MRI pipeline using PyPulseq, BART reconstruction, and a gradient probing sequence was developed and validated for muscle velocity, displacement, and strain analysis on two Siemens systems during NMES contractions.
T0 review reviewed 2026-06-27 challenge →
load-bearing objection The paper delivers a practical open-source 4D PC-MRI pipeline with time savings and strain metrics, but the vendor-agnostic claim rests on untested extrapolation from two Siemens systems. the 2 major comments →
Vendor-agnostic 4D Phase Contrast MRI: a complete open-source pipeline for velocities, displacement, and strain analysis
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
Core claim
The pipeline demonstrates that a compressed sensing-accelerated 4D PC-MRI sequence implemented in PyPulseq, reconstructed via BART, and paired with a gradient probing sequence produces consistent velocity and strain measurements across two different 3T Siemens scanners, enabling extraction of muscle strain parameters from NMES contractions in under 11 minutes per scan.
What carries the argument
The gradient probing sequence that ensures correct velocity sign assignment across scanner orientations and vendors, together with the PyPulseq-implemented compressed sensing sequence and BART reconstruction for 4D flow data.
Load-bearing premise
The gradient probing sequence ensures correct velocity sign assignment and the results from two Siemens 3T systems generalize to other vendors and field strengths.
What would settle it
Applying the pipeline to a non-Siemens scanner or different field strength and observing mismatched velocity directions or strain values that cannot be corrected by the probing sequence.
If this is right
- Strain maps and sigmoid-fitted curves yield extractable values for peak strain, mean strain, and buildup rate during muscle contraction.
- Vastus lateralis strains reach median peak values around 0.49 compared to 0.063 in flexor digitorum superficialis.
- Compressed sensing reduces forearm and thigh acquisition times to 5 and 11 minutes respectively.
- The framework supports quantitative muscle imaging with demonstrated consistency on two distinct 3T Siemens platforms.
Where Pith is reading between the lines
- The open pipeline could support direct comparison of muscle function metrics between independent research groups without vendor-specific processing steps.
- Extending the gradient probing approach to additional field strengths might allow strain measurements in pediatric or high-field settings with minimal recalibration.
- Integration with other open-source tools for NMES modeling could link electrical stimulation parameters directly to measured strain buildup rates.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript presents a fully open-source 4D phase contrast MRI pipeline that integrates a compressed-sensing accelerated sequence implemented in PyPulseq, BART-based reconstruction, and a gradient probing sequence for velocity sign assignment. It reports validation on two Siemens 3T systems (Prisma and Vida Fit) for forearm (n=9) and thigh (n=10) muscle strain analysis during NMES-induced contractions, with acquisition time reductions from 35/80 min to 5/11 min and extraction of peak/mean strain and buildup rate via sigmoid fitting.
Significance. If the vendor-agnostic aspects hold, the work would offer a significant contribution by supplying a reproducible open-source framework for quantitative muscle imaging. The provision of a fully open-source pipeline with PyPulseq and BART integration, together with the demonstrated acceleration via compressed sensing, are explicit strengths that promote accessibility and standardization in the field.
major comments (2)
- [Abstract] Abstract and validation description: the central claim of a 'vendor-agnostic' pipeline (title and abstract) rests on data acquired exclusively on two Siemens 3T systems; no acquisitions, sign-consistency checks, or strain metrics from GE, Philips, or other field strengths are presented, directly undermining the generalization required by the title and the gradient-probing motivation.
- [Methods] Methods (gradient probing sequence) and Results: the gradient probing sequence is stated to ensure correct velocity sign assignment across vendors, yet all reported data (n=9 arm, n=10 leg) and sign validation remain confined to the same Siemens pair, leaving transferability of sequence timing and reconstruction assumptions untested.
minor comments (1)
- [Results] Results (strain metrics): the reported median peak strain values (0.49 vs. 0.063) and order-of-magnitude difference lack accompanying error bars, fit uncertainties, or statistical comparisons, which would strengthen the quantitative claims even if the primary validation scope is addressed.
Simulated Author's Rebuttal
We thank the referee for the constructive comments on the scope of the vendor-agnostic claims. We address each point below, acknowledging the limitations of our current validation dataset while defending the design rationale based on the open-source tools employed.
read point-by-point responses
-
Referee: [Abstract] Abstract and validation description: the central claim of a 'vendor-agnostic' pipeline (title and abstract) rests on data acquired exclusively on two Siemens 3T systems; no acquisitions, sign-consistency checks, or strain metrics from GE, Philips, or other field strengths are presented, directly undermining the generalization required by the title and the gradient-probing motivation.
Authors: We agree that the empirical validation is confined to two Siemens 3T systems and that this limits the strength of the generalization in the title and abstract. The vendor-agnostic framing derives from the use of PyPulseq (which generates sequences portable across vendors without proprietary pulse sequence environments), BART (vendor-independent reconstruction), and the gradient probing sequence (intended to resolve velocity sign conventions that can differ by vendor or orientation). We will revise the abstract, title if appropriate, and discussion to state explicitly that the pipeline is designed to be vendor-agnostic and demonstrated on Siemens platforms, with the gradient probing method positioned as enabling future cross-vendor use. No additional scanner data will be added. revision: partial
-
Referee: [Methods] Methods (gradient probing sequence) and Results: the gradient probing sequence is stated to ensure correct velocity sign assignment across vendors, yet all reported data (n=9 arm, n=10 leg) and sign validation remain confined to the same Siemens pair, leaving transferability of sequence timing and reconstruction assumptions untested.
Authors: The gradient probing sequence acquires short calibration data to determine velocity sign polarity independently of vendor-specific assumptions about gradient directions or coordinate systems; it is implemented entirely in PyPulseq and uses the same BART reconstruction pipeline. Sign consistency was verified on the Siemens datasets to confirm the method functions as intended. We acknowledge that transferability of timing parameters and reconstruction assumptions to other vendors has not been empirically tested. In revision we will expand the methods to detail the general design and add an explicit limitations paragraph noting the need for future multi-vendor validation. revision: partial
- Empirical validation of sequence portability, sign assignment, and strain metrics on non-Siemens vendors or other field strengths is unavailable.
Circularity Check
No circularity: empirical pipeline validation with no derivation chain
full rationale
The paper describes construction and empirical validation of an open-source 4D PC-MRI pipeline (PyPulseq sequence, BART reconstruction, gradient probing, strain analysis) on two Siemens 3T systems. No mathematical derivations, first-principles results, parameter fittings, or uniqueness theorems are claimed. All quantitative outputs (strain maps, peak/mean strain, buildup rate) are direct measurements from acquired data rather than quantities predicted from fitted inputs or self-referential definitions. The vendor-agnostic claim rests on experimental compatibility checks within the tested systems, not on any loop that reduces to the paper's own assumptions by construction. This is a standard methods-and-validation manuscript with no load-bearing self-citation chains or ansatz smuggling.
Axiom & Free-Parameter Ledger
free parameters (1)
- sigmoid curve parameters for strain
axioms (2)
- standard math Phase contrast MRI signal phase encodes velocity linearly under standard gradient conditions
- domain assumption Compressed sensing reconstruction faithfully recovers the accelerated 4D flow data
Cite this review
Pith. "Pith review of Vendor-agnostic 4D Phase Contrast MRI: a complete open-source pipeline for velocities, displacement, and strain analysis." pith.science (2026). https://pith.science/paper/UYUOEBUF
@misc{pith2026260609444,
author = {Pith},
title = {Pith review of: Vendor-agnostic 4D Phase Contrast MRI: a complete open-source pipeline for velocities, displacement, and strain analysis},
year = {2026},
howpublished = {\url{https://pith.science/paper/UYUOEBUF}},
note = {Machine review of arXiv:2606.09444}
}
read the original abstract
Phase contrast MRI (PC MRI) enables quantitative assessment of tissue motion and strain. Although it is increasingly used, standardized, vendor-agnostic pipelines for accelerated acquisitions remain scarce. We present a fully open-source 4D flow PC-MRI pipeline integrating a compressed sensing-accelerated sequence implemented in PyPulseq, BART-based reconstruction, and strain analysis. Additionally, a gradient probing sequence was developed to ensure correct velocity sign assignment across scanner orientations and vendors. The pipeline was validated across two Siemens MRI systems (3T MAGNETOM Prisma and 3T Vida Fit) in two anatomical applications: forearm (Flexor Digitorum Superficialis, n=9) and thigh (Vastus Lateralis, n=10) during Neuromuscular Electrical Stimulation (NMES)-induced contractions. Compressed sensing reduced acquisition times from 35 and 80 minutes to 5 and 11 minutes for the arm and leg acquisitions, respectively. Muscle strain maps and sigmoid-fitted strain curves enabled extraction of peak strain, mean strain, and buildup rate. Strains in the Vastus Lateralis were approximately one order of magnitude higher than in the Flexor Digitorum Superficialis (median peak strain 0.49 vs. 0.063, mean strain 0.31 vs. 0.031). The pipeline demonstrates multi-platform compatibility and provides a reproducible, open framework for quantitative muscle imaging.
Reference graph
Works this paper leans on
-
[1]
Phase-Contrast MRI: Physics, Techniques, and Clinical Applications,
D. T. Wymer, K. P. Patel, W. F. Burke, and V. K. Bhatia, “Phase-Contrast MRI: Physics, Techniques, and Clinical Applications,” RadioGraphics , vol. 40, no. 1, pp. 122–140, Jan. 2020, doi: 10.1148/rg.2020190039. [2] M. Markl, A. Frydrychowicz, S. Kozerke, M. Hope, and O. Wieben, “4D flow MRI,” J. Magn. Reson. Imaging JMRI , vol. 36, no. 5, pp. 1015–1036, N...
-
[2]
ORMIR-MIDS: an open standard for curating and sharing musculoskeletal imaging data,
F. Santini et al. , “ORMIR-MIDS: an open standard for curating and sharing musculoskeletal imaging data,” JBMR Plus , vol. 10, no. 3, p. ziag013, Mar. 2026, doi: 10.1093/jbmrpl/ziag013. [25] N. J. Pelc, M. A. Bernstein, A. Shimakawa, and G. H. Glover, “Encoding strategies for three-direction phase-contrast MR imaging of flow,” J. Magn. Reson. Imaging JMRI...
This paper was first reviewed by grok-4.3 on June 27, 2026.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.