Pith. sign in

REVIEW 4 major objections 4 minor 3 references

Theoretical derivation of blood velocity from TOF-MRA based artery centerline

T0 review · 4 major / 4 minor · reviewed 2026-08-01 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Quantitative blood velocity can be extracted from the spatial decay of standard TOF-MRA centerline signals.

desk verdict The derivation is sound and the inversion is well-posed, but the validation is an oracle experiment: regularization parameters are tuned on ground-truth velocities, so the central claim about extracting functional velocity from clinical TOF-MRA is not yet supported. read the letter →

arxiv 2607.16498 v3 pith:3BN2JWHK submitted 2026-07-17 physics.med-ph math-phmath.MP

Abrar Faiyaz This is my paper · ORCID
classification physics.med-phmath-phmath.MP
keywords TOF-MRAbloodvelocityBloch-McConnellequationsflow-relatedenhancementeffectiverelaxationrateDual-Tikhonovregularizationcerebrovascularimaging
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

Standard Time-of-Flight MRA is normally read as a purely structural image, but this paper argues that the spatial decay of its signal along a vessel contains quantitative blood velocity information. By deriving the Bloch-McConnell flow equations into a one-dimensional decay law, the author shows that a vessel's local speed can be inverted from the intensity slope at each centerline point, stabilized by a dual-penalty regularized optimization. Simulations in a tapering vessel and a focal stenosis under synthetic scanner noise recover ground-truth point-wise velocities, with errors of roughly 2 cm/s for gradual narrowing and 8 cm/s for the sharp stenotic jet. If the framework holds in vivo, velocity measurement would be available from a standard clinical MRA without phase-contrast, multi-echo, or additional scan time.

What carries the argument

The key machinery is the effective relaxation rate R_eff = 1/T1 - ln(cosα)/TR, which folds T1 relaxation and RF saturation into a single spatial decay constant, together with the inversion identity v(z) = -R_eff (S(z)-S_ss)/S'(z). The forward model maps velocity history into signal through the transit-time integral, and the dual-Tikhonov global matrix inversion (Eq. 18) replaces unstable numerical differentiation with a stabilized linear solve, allowing pointwise instead of constant velocity.

What would settle it

Measure velocity in a straight uniform vessel with phase-contrast MRI, extract the centerline TOF-MRA signal from the same subject and sequence, and check whether v(z) computed from Eq. 14 matches the independent measurement; if the recovered velocities deviate by much more than the simulated RMSE (≈2 cm/s tapered, ≈8 cm/s stenotic) at the same noise level, the idealized decay model is falsified for real data.

Watch

Extended reading notes

Core claim

The paper's central discovery is a mathematical identity linking the spatial decay of TOF-MRA signal intensity to local blood velocity. The author models a blood packet entering the imaging volume fully relaxed and then being repeatedly excited by RF pulses, which yields an effective relaxation rate R_eff = 1/T1 - ln(cosα)/TR that combines true T1 relaxation with artificial RF-driven saturation. The resulting centerline signal obeys S(z) = S_ss + (S0 - S_ss) exp(-R_eff ∫ dz'/v(z')), and differentiating the log-transformed signal gives the pointwise velocity v(z) = -R_eff (S(z)-S_ss)/S'(z). Because this derivative amplifies noise, the inversion is recast as a globally convex least-squares pro

Load-bearing premise

The load-bearing premise is that the measured TOF-MRA centerline signal is produced exactly by the 1D plug-flow Bloch-McConnell model—fully relaxed inflow at z=0, identical RF pulses, no partial-volume averaging, slice-profile, inflow-angle, or T2* effects, and known S0 and S_ss.

Editorial extensions

If this is right

  • Pointwise blood velocity becomes measurable from a standard ~4-minute TOF-MRA acquisition, eliminating the need for phase-contrast or 4D-flow scans that add 5-20 minutes.
  • Previously acquired TOF-MRA images could be re-analyzed retrospectively, since the inversion needs only the centerline signal profile and the known sequence parameters (TR, flip angle, T1).
  • The method distinguishes smooth tapering from focal stenosis by reconstructing local velocity, where global constant-velocity fits would only return a mean value.
  • The global inverse problem is convex with a closed-form solution, so no iterative non-linear fitting is required and the computation can be made fast enough for clinical workflows.
  • The in vivo centerline decay shown in the paper matches the tapering-tube model, suggesting the framework can be applied to real cerebrovascular networks once the inverse solver is run on those topologies.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • A natural next experiment is to compare the solver's output against phase-contrast MRI or transcranial Doppler in the same vessel segment; any systematic mismatch would identify which neglected physics must be added to the model.
  • The plug-flow assumption ignores the parabolic or pulsatile profile of real arterial flow; extending the transit-time integral to a distributed velocity field would test the method's sensitivity to flow shape.
  • The inversion is not limited to TOF-MRA: any inflow-weighted contrast that produces a spatial decay profile, such as arterial spin labeling, could be treated with the same log-derivative identity.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 4 minor

Summary. The paper proposes a physics-informed framework to estimate point-wise blood velocity from the 1D centerline signal of standard TOF-MRA. Starting from the conventional Bloch-McConnell expression for flow-related enhancement, the author derives a variable-velocity formulation in which the spatial decay of the signal depends on the transit-time integral (Eq. 10). A dual-Tikhonov regularized inversion is set up to recover the inverse-velocity array from the log-transformed signal (Eqs. 15–18). The method is tested on two synthetic tube geometries—a gradually narrowing tube and a focal stenosis—with ground-truth velocities derived from mass conservation. The paper reports RMSEs of 1.94 cm/s (Scenario A) and 7.71 cm/s (Scenario B), claims these recover the variable velocity profiles, and concludes that quantitative hemodynamics can be extracted from structural TOF-MRA without additional scans. An in vivo example (Figure 3) shows only signal intensity decay along traced centerlines, not recovered velocities.

Significance. If the framework were validated against independent measurements, it could offer a low-cost way to augment routine TOF-MRA with hemodynamic information, avoiding phase-contrast or multi-TE acquisitions. The analytical derivation in Sections 2.1–2.2 is transparent and appears algebraically correct, and the regularized inversion is a reasonable, convex approach to an ill-posed integral equation. The author also openly acknowledges limitations, including the need for data-driven regularization selection and the amplitude attenuation of stenotic peaks. However, the current validation is entirely self-referential: synthetic signals are generated from the same model being inverted, and the regularization weights are chosen using the ground-truth velocity. The reported error figures are thus oracle bounds rather than achievable performance. The central claim—that quantitative velocity can be recovered from real TOF-MRA—is not yet supported.

major comments (4)
  1. [Section 3.2 / Table 1 / Discussion (p. 8/13)] The validation is an oracle experiment. The regularization weights λ1 and λ2 are selected by sweeping the parameter space and choosing the combination with the lowest RMSE against the ground-truth v(z) used to synthesize the signal. The reported RMSEs (1.94 cm/s and 7.71 cm/s) therefore represent an upper bound on achievable accuracy only when the answer is already known. The paper states that 'automated, data-driven parameter selection (e.g., L-curve methods) will be necessary for large-scale clinical application,' which concedes that the current workflow is not blind. This undermines the quantitative claims in the abstract and conclusions.
  2. [Section 2.4 and Eq. 10] The signal model in Eq. 10 assumes plug flow: a coherent packet of blood moves with a single centerline velocity v(z) and experiences the same RF pulse history. In real TOF-MRA, the voxel signal is an average over a parabolic velocity distribution, partial-volume effects, slice profile, inflow angle, and T2* decay. The 'parabolic laminar flow constraint' mentioned in Section 2.4 is applied only to the velocity field used to generate ground truth; it does not enter the signal generation model. Thus the synthetic validation provides no evidence that actual TOF-MRA centerline intensities obey Eq. 10. A multi-compartment simulation or an in-vivo comparison with phase-contrast velocities would be needed.
  3. [Eqs. 14 and 16] The inversion requires S0 (fresh-blood signal) and Sss (steady-state signal) to be known. These quantities are not directly available in standard clinical TOF-MRA, and they vary spatially due to B1 inhomogeneity, coil sensitivity, and tissue properties. The paper does not propose an estimation strategy for S0 and Sss. Without such a strategy, the practical applicability of Eq. 14 and Eq. 16 is not established; the synthetic experiments assume these values are known exactly.
  4. [Figure 3 and Abstract] The in vivo illustration demonstrates only that signal intensity decreases along extracted centerlines; no velocity profile is recovered or compared with an independent reference. The abstract's claim that 'quantitative, localized functional hemodynamic metrics can be extracted from standard structural MRA imaging' is therefore not supported for real data. The paper is appropriately framed elsewhere as a 'proof-of-concept,' but the in vivo figure does not provide any quantitative validation.
minor comments (4)
  1. [Eq. 17] The superscript T on D1 and D2 is not defined; it should be stated that it denotes matrix transpose.
  2. [Table 1 / Section 3.2] For Scenario B, the lowest RMSE (7.71 cm/s) occurs at λ1=1000, λ2=10, but the text selects λ1=5000, λ2=100 as 'optimal' based on qualitative absence of ringing. The distinction between RMSE-optimal and qualitatively-optimal should be stated explicitly.
  3. [General formatting] The manuscript retains MDPI template placeholders such as 'Academic Editor: Firstname Last-name' and the DOI '10.3390/xxxxx'. These must be removed or completed before submission.
  4. [Section 3.1] The ground-truth velocities are derived from a 1D mass-conservation formula, not from a full computational fluid dynamics simulation. The term 'simulated hemodynamics' should be qualified to avoid overstating the physical realism of the test scenarios.

Circularity Check

2 steps flagged · score 7.0 of 10

Validation is an oracle experiment: synthetic signals are generated from the same Bloch–McConnell equation being inverted, and regularization parameters are selected using the known ground-truth velocities, so the reported RMSEs are in-sample rather than evidence for the central claim.

  1. fitted input called prediction [Section 3.2, Table 1, Figure 2]
    "The parameter sweep evaluated the mathematical trade-off between first-derivative noise suppression (λ1) and second-derivative structural stiffness (λ2): ... The lowest absolute RMSE (7.71 cm/s) was recorded at the weakest regularization bounds (λ1=1000, λ2=10); however, qualitative observation of the reconstructed profile revealed persistent high-frequency baseline oscillations. Moderate regularization weights (e.g., λ1=5000, λ2=100) established an optimal balance, effectively mitigating oscillatory ringing artifacts in the stable proximal and distal segments while preserving the magnitude of"

    The regularization parameters λ1 and λ2 are selected by comparing reconstruction RMSE against the known ground-truth velocity profiles. The same ground-truth velocities were used to generate the synthetic signals, so the reported errors (1.94 cm/s and 7.71 cm/s) are oracle-optimized, not blind predictions. The paper itself concedes that automated, data-driven parameter selection is future work, confirming that the current accuracy figures do not represent achievable performance without knowing the answer.

  2. other [Section 2.2 (Eq. 10) and Section 2.5 (Eqs. 15–18)]
    "Because the blood is experiencing RF pulses over this entire transmit time, the standard exponential decay equation changes from a scalar multiplication to an integral of the local decay-rate. S(z)=Sss+(S0−Sss)exp(−Reff∫0^z 1/v(z′)dz′) ... By spatially upsampling the centerline to a high-density computational grid, the Bloch equation transit time integral can be mapped as a linear system Lu≈y."

    The synthetic TOF-MRA signals are created by evaluating the same transit-time integral (Eq. 10) that the inverse solver discretizes as Lu≈y. Inverting a signal generated from this model therefore recovers the very velocities that were put into the model; it demonstrates numerical self-consistency, not that real TOF-MRA signals obey the assumed Bloch–McConnell plug-flow model. The in vivo Figure 3 shows only intensity decay, not recovered velocity, so no independent validation is provided.

full rationale

The algebraic derivation in Sections 2.1–2.3 is self-contained and not circular: Eq. 14 follows from differentiating Eq. 10, and the definition of Reff is a parameterization of the standard Bloch–McConnell exponential. No load-bearing self-citation or imported uniqueness theorem is used; the cited Nishimura reference is standard MR physics. However, the central claim — that quantitative localized velocity can be extracted from standard TOF-MRA — rests on a validation loop that closes on itself. Simulated signals are generated from the exact forward model being inverted, and the regularization hyperparameters are tuned using the same ground-truth velocities used to generate the data. The reported RMSEs are therefore in-sample, oracle bounds. The in vivo component only illustrates signal decay along centerlines and does not recover or benchmark velocities. This is a significant circular validation, though not a circular derivation, which justifies a score of 7 rather than higher.

Assumptions & free parameters 3 free parameters · 6 assumptions · 1 invented entities

The derivation is algebraically straightforward and introduces no new physics; R_eff is a composite of existing scanner parameters. The main burden is the 1D plug-flow measurement model and the ground-truth-dependent choice of regularization weights, not new entities.

free parameters (3)
  • Regularization weights (λ1, λ2) = Optimal set λ1=5000, λ2=100; swept over λ1∈[1000,100000], λ2∈[10,1000000]
    Selected by minimizing RMSE against simulated ground-truth velocity; this is fitting to the target and requires an independent selection method (e.g., L-curve) for real data.
  • Maximum physiological velocity bound = 150 cm/s
    Hard cap applied after inversion to enforce stability; not derived from data.
  • S0 and Sss (fresh-blood and steady-state signal amplitudes) = Known exactly in simulation; no estimation method provided for real images
    Eqs. 14 and 16 require both values; in real TOF-MRA they must be estimated from the image, and the paper does not specify a robust procedure.
assumptions (6)
  • domain assumption Ernst/steady-state magnetization evolution for repeated identical RF pulses (Eqs. 1–2)
    Borrowed from textbook MR physics; assumes a single blood packet experiences identical α pulses at constant TR.
  • domain assumption Continuous pulse approximation n = t(z)/TR, with t(z)=∫0^z dz'/v(z')
    Used to convert discrete pulse count into a spatial integral (Eqs. 3, 9–10); ignores discrete RF timing and pulsatile flow.
  • domain assumption 1D plug-flow centerline representation: S(z) tracks one blood packet with no partial volume, slice profile, inflow-angle, T2*, or bifurcation mixing
    The core measurement model (Eqs. 5–7); unvalidated in vivo. Figure 3 shows only intensity decay, not velocity recovery.
  • domain assumption Fully relaxed blood entry at z=0 (M_z(0)=M0)
    Required for the boundary condition in Eq. 1; may fail at vessel origins or for slow flow.
  • domain assumption Steady, incompressible flow with mass conservation v(z) ∝ 1/A(z) for simulated ground truth
    Used to generate ground-truth velocities in Section 2.4; not a property of the inversion itself.
  • standard math Fundamental Theorem of Calculus applies to the sampled, noise-corrupted signal with S(z) ≠ Sss
    Eqs. 12–14 differentiate a noisy observation; requires sufficient smoothness and non-zero denominator, which motivates the regularized inverse in Section 2.5.
invented entities (1)
  • Effective relaxation rate R_eff = 1/T1 − ln(cos α)/TR (Eq. 4) independent evidence
    purpose: Consolidates sequence/scanner parameters into a single decay constant so that the spatial TOF signal decay can be inverted for velocity.
    Not a new physical entity: it is an algebraic combination of existing MR parameters. It is independently measurable from scanner settings or phantom experiments, so it adds no hidden freedom beyond the assumption that the Ernst/flow model holds.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Theoretical derivation of blood velocity from TOF-MRA based artery centerline." pith.science (2026). https://pith.science/paper/3BN2JWHK

@misc{pith2026260716498,
  author       = {Pith},
  title        = {Pith review of: Theoretical derivation of blood velocity from TOF-MRA based artery centerline},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/3BN2JWHK}},
  note         = {Machine review of arXiv:2607.16498}
}
read the original abstract

Time-of-flight magnetic resonance angiography (TOF-MRA) is widely used for structural vascular imaging, but extracting functional hemodynamics like blood velocity typically requires supplementary phase-contrast scans. This study proposes a novel, physics-informed computational framework to extract variable fluid velocity directly from standard TOF-MRA signal profiles. We analytically expand the Bloch equations into Bloch-McConnell flow equations, establishing a mathematical relationship between the spatial decay of longitudinal magnetization and fluid velocity. To validate this derivation and overcome the limitations of constant-velocity assumptions, a MATLAB simulation framework was developed to model fluid flow in two variable-geometry flowing tube cases i.e continuous tapering and focal stenosis -under synthetic scanner noise. A global inverse optimization approach utilizing Dual-Tikhonov regularization was deployed to stably invert the ill-posed transit time integral, actively penalizing high-frequency numerical ringing while preserving structural curve stiffness. The computational sim-ulations successfully recovered ground-truth point-wise velocities, accurately tracking gradual hemodynamic accelerations and sharp stenotic jets. This theoretical framework provides a robust mathematical proof-of-concept that quantitative, localized functional hemodynamic metrics can be extracted from standard structural MRA imaging, estab-lishing a foundation for advanced flow quantification without requiring additional scan time.

Figures

Figures reproduced from arXiv: 2607.16498 by the authors.

Figure 1
Figure 1. Scenario A) Narrowing/Tapering Flow. B) Stenotic Flow; velocity field shown in perula colormap for both cases. Simulation conducted for an assumed artery segment of 12cm. • Scenario A (Narrowing Tube): The estimation achieved high fidelity across most of the parameter space, with root-mean-square error (RMSE) ranging from 1.94 cm/s to 4.21 cm/s. Because the true velocity undergoes a gradual, low-frequency transition… view at source ↗
Figure 2
Figure 2. Evaluation of Dual-Tikhonov Regularization for Variable Velocity Estimation. (Top Row) Reconstructed spatial velocity profiles for the continuous tapering tube (Scenario A, left) and the focal stenotic tube (Scenario B, right). The estimations are derived across a defined sweep of first derivative (λ1) and second-derivative (λ2) penalty weights and compared against both the theoret￾ical ground-truth variable velocit… view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

3 extracted references

  1. [1]

    After absorbing an integer history of n successive radiofrequency pulses, the moving packet's transient ap-proach toward the steady state is described by (9)- 𝑀!(𝑛)=𝑀##+(𝑀

    Introduction To maintain metabolic and neural function, the human brain requires continuous and tightly regulated blood flow. Consequently, the quantitative assessment of cerebro-vascular hemodynamics is important for the diagnosis and longitudinal management of neurovascular pathologies, including cerebral small vessel disease (CSVD) (1,2). In stand-ard ...

    arXiv 2026

  2. [3]

    Leung J, Behpour A, Sokol N, Mohanta A, Kassner A. Assessment of intracranial blood flow velocities using a computer controlled vasoactive stimulus: A comparison between phase contrast magnetic resonance angiography and transcranial doppler ultraso-nography. J Magn Reson Imaging. 2013;38(3):733–8. doi:10.1002/jmri.23911 4. Lim RP, Shapiro M, Wang EY, Law ...

    arXiv 2013

  3. [17]

    Transcranial Doppler Ultrasonography as a Diagnostic Tool for Cerebrovascular Disorders

    Pan Y, Wan W, Xiang M, Guan Y. Transcranial Doppler Ultrasonography as a Diagnostic Tool for Cerebrovascular Disorders. Front Hum Neurosci. 2022 Apr 29;16. doi:10.3389/fnhum.2022.841809 18. Haller S, Zaharchuk G, Thomas DL, Lovblad KO, Barkhof F, Golay X. Arterial Spin Labeling Perfusion of the Brain: Emerging Clinical Applications. Radiology. 2016 Nov;28...

    arXiv 2022

Pith tools

Reviewed August 1, 2026 · model on record in the stance chip above.