Pith. sign in

REVIEW 4 major objections 5 minor 11 references

Estimation of Blood Flow Parameters in the Left Atrial Appendage from 4DCT Dynamic Contrast Enhancement

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

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

Pith's one-line read Fit each 4DCT voxel's contrast curve with a gamma-variate model and the result is a spatial map of blood arrival and residence time in the left atrial appendage.

desk verdict A feasible CT-only pipeline for mapping LAA contrast dynamics, but the quantitative flow claims rest on unvalidated gamma-variate fits and two illustrative cases. read the letter →

arxiv 2502.02728 v1 pith:WUMWLFPX submitted 2025-02-04 eess.IV

Lauren M. Severance , Andrew M. Kahn , Juan C. del Alamo , Elliot R. McVeigh This is my paper · ORCID
classification eess.IV
keywords dynamiccontrastenhancement4DCTleftatrialappendagegamma-variatefitresidencetimefibrillationbloodstasisCTperfusion
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

The paper proposes that the time-varying CT signal from a single contrast-enhanced 4DCT scan can be read as a local indicator-dilution curve throughout the left atrial appendage (LAA). At every voxel, the 27 timepoints of contrast enhancement are fit with a gamma-variate curve after nonrigid motion correction, and the fitted curve yields two flow parameters: arrival time (when signal reaches 1% of its peak) and residence time (the intensity-weighted mean transit time). The fitted parameters are presented as smooth two-dimensional maps that show where contrast arrives late and lingers, which the authors interpret as a direct visualization of blood stasis in atrial fibrillation patients. If correct, this gives clinicians a way to assess stroke risk from an ECG-gated CT already performed for cardiac morphology, without solving fluid-dynamics equations. The claim is that these maps capture spatial-temporal flow characteristics across the LAA, not that CT measures flow directly.

What carries the argument

The load-bearing object is the gamma-variate model, a three-parameter curve widely used to describe the passage of an injected indicator bolus through a series of compartments. Written as $y(t) = y_{\max}\,(\alpha e)^{-\alpha}(t/t_{\mathrm{peak}})^\alpha - y_b$, it turns each voxel's noisy 27-point contrast time series into a smooth analytical curve. The fitted curve supplies direct parameters ($y_{\max}$, $t_{\mathrm{peak}}$, $\alpha$) and derived transit parameters ($t_a$ and $RT$) that become the values plotted on dynamic contrast enhancement maps. The model's role is to regularize the sparse, noisy temporal signal so that per-voxel flow parameters can be mapped without averaging away spatial detail.

What would settle it

Compare the residence-time maps against LAA flow velocity measured by transesophageal echocardiography or 4D-flow MRI in the same atrial-fibrillation patients: if regions flagged as high residence time show fast contrast washout on those modalities, the maps do not represent stasis. A simpler mechanical check is a flow phantom with known chamber transit times imaged with the same 27-timeframe protocol, requiring the fitted residence times to match the known values within tolerance.

Watch

Extended reading notes

Core claim

On the paper's own terms, the discovery is that voxel-wise gamma-variate fitting of 4DCT dynamic contrast enhancement produces clinically interpretable maps of LAA blood flow parameters. For each voxel, the contrast signal is baseline-corrected using the mean of the first three timepoints in the descending aorta, then fit to the gamma-variate function $y(t) = y_{\max}\,(\alpha e)^{-\alpha}(t/t_{\mathrm{peak}})^\alpha - y_b$, with fit parameters $\alpha$, $y_{\max}$, and $t_{\mathrm{peak}}$. From the fitted curve the authors derive arrival time $t_a$, defined as the time the signal reaches 1% of its peak, and residence time $RT = \int t\,y(t)\,dt / \int y(t)\,dt$. Applied across a multiplanar slice through the LAA, these parameters form smooth maps, and in two illustrative patients the maps separate a uniform-filling phenotype from a delayed-filling phenotype in the middle and distal LAA. The paper frames this as enabling quantification and visualization of spatial-temporal flow characteristics and, potentially, thrombus risk.

Load-bearing premise

The premise the method depends on is that after motion correction and subtraction of a descending-aorta baseline, each voxel's CT intensity over time tracks the local contrast-agent concentration faithfully enough that the gamma-variate fit returns true arrival and residence times for LAA blood.

Editorial extensions

If this is right

  • A single ECG-gated 4DCT acquisition during a standard contrast injection can yield voxel-resolved maps of arrival time and residence time across the LAA, with no CFD simulation required.
  • Maps that are uniform across the LAA (as in Patient 1) and maps with delayed distal filling (as in Patient 2) give a direct visual readout of filling patterns that may correspond to differing thrombus risk.
  • Because the mean radiation dose is comparable to other cardiac CT exams, the approach could be added to existing clinical CT protocols for AF patients without a separate imaging study.
  • The success of fits from 27 sparse timepoints suggests the temporal sampling rate could be reduced, with the dose lowered further if fewer timeframes prove adequate.

Reading between the lines

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

  • A natural next step, not stated in the paper, is to use the ratio $RT / t_a$ or the peak-normalized curve width as a dimensionless stasis index that could be compared across patients independent of heart rate and injection timing.
  • The same per-voxel gamma-variate pipeline applied to the full left atrial volume, rather than a single 2D slice, would produce 3D stasis maps; the authors note they are working toward this with a machine-learning approach.
  • The maps' clinical value hinges on an external anchor: comparing residence-time maps with LAA emptying velocity from transesophageal echocardiography or 4D-flow MRI in the same patients would test whether 'slow filling' on CT matches 'stasis' by established measures.
  • If the gamma-variate residuals are structured (e.g., high near the LAA wall), the maps may be capturing a mixture of motion-correction error and true slow flow; inspecting residual maps would separate these.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

4 major / 5 minor

Summary. The paper presents a method to estimate blood flow parameters in the left atrial appendage (LAA) from dynamic contrast-enhanced 4DCT. A 4DCT acquisition protocol with 27 timeframes is used, followed by nonrigid registration and multiplanar reconstruction. Each voxel's time-intensity curve is fit with a gamma-variate model, and parameters including arrival time and residence time (RT) are derived from the fitted curve to create 2D parameter maps. The authors show qualitative results for two patients and state that the maps enable quantification and visualization of flow parameters across the LAA, with potential implications for stroke risk assessment in atrial fibrillation patients.

Significance. If validated, the method could provide a clinically accessible, high-spatial-resolution surrogate for blood stasis in the LAA from routine CT data, which would be valuable for stroke risk stratification. The paper demonstrates a clear technical pipeline and the use of a widely accepted indicator-dilution model. The maps shown appear visually plausible, and the acquisition protocol uses a radiation dose comparable to standard cardiac CT. However, the current evidence is preliminary: no independent validation against TEE, 4D flow MRI, CFD, or a phantom is presented, and the quantitative claim rests on an unvalidated baseline subtraction and extrapolated gamma-variate tails. The cohort of 17 patients is mentioned but results are only shown for 2, with no quantitative summary or error analysis.

major comments (4)
  1. [Section II.D and Eq. (2)] Residence time is computed as an integral over the fitted gamma-variate curve using Eq. (2). The 27 timepoints are acquired over 40–60 s plus one at ~100 s, so for voxels with slow washout—precisely the stasis regions of interest—the tail of the curve is extrapolated and not constrained by measured data. This makes the derived RT potentially an artifact of the model's exponential decay rather than a measured physiological quantity. Please provide per-voxel fit uncertainty estimates (e.g., confidence intervals from residuals) or restrict the analysis to voxels where the washout portion is actually sampled.
  2. [Section II.D, baseline subtraction] The baseline y0 is computed as the mean of timepoints 1–3 in a descending aorta ROI. This assumes that the LAA voxel's pre-contrast CT intensity is identical to the aortic baseline, which is not validated. A constant offset error in y0 will bias the fitted y_max, t_arrival, and RT. The authors should validate this assumption by comparing the pre-contrast LAA voxel values to the aortic baseline, or by using a phantom with known contrast concentrations to establish the accuracy of the baseline correction.
  3. [Section III.B and Fig. 4] The paper reports that 17 patients underwent imaging, but quantitative results are shown for only 2 patients. There is no summary of the range of RT values, inter-patient variability, or the reproducibility of the maps, nor any comparison with independent flow measurements (e.g., TEE LAA emptying velocity, 4D flow MRI, or CFD). The claim that the method 'enables quantification and visualization of spatial-temporal characteristics of flow parameters' is therefore not supported by the evidence presented. Please include a quantitative validation study or a statistically meaningful analysis across the full cohort, and compare the derived parameters against an independent reference.
  4. [Table I and Section II.D] The gamma-variate fit uses loose parameter bounds (alpha unbounded, y_max and t_arrival within ±40 HU and ±5 s of the data maximum). No assessment of fit quality (e.g., R^2, chi-square, or residual analysis) is reported per voxel or per patient. It is unclear whether the smooth maps reflect physiological variation or simply the smoothness of the fitting function. Please report fit quality metrics and consider adding a goodness-of-fit criterion to exclude poorly fitted voxels.
minor comments (5)
  1. [Abstract] The word '4-dimentional' should be 'four-dimensional'.
  2. [Section II.B] The threshold of cumulative sum <100 HU for voxel exclusion is introduced without justification; please explain how this value was chosen and its sensitivity.
  3. [Section II.D] The parameter t_arrival is not clearly introduced before it is used; please define it explicitly and consistently (the text appears to have formatting artifacts for the symbol).
  4. [Fig. 1] The term '2D fit over time' is ambiguous; please clarify whether this refers to fitting the temporal signal at each pixel or to a spatial 2D fit.
  5. [References] In reference [10], a comma is missing after 'Z.M. Lin'; please verify the reference formatting.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: all reported parameters are explicit functions of the gamma-variate fit to measured 4DCT contrast curves, and the only author-overlapping citation is non-load-bearing background.

full rationale

The derivation chain is self-contained. Section II.D fits Eq. (1) to each voxel's measured contrast-versus-time curve (raw 4DCT signal minus a descending-aorta baseline), then computes t_arrival as the time the fitted curve reaches 1% of its maximum and RT as the first moment of the analytic gamma-variate curve via Eq. (2). These are explicit algebraic definitions in terms of the fitted function; the paper transparently labels them 'derived parameters,' so there is no self-definitional inversion in which a target quantity is secretly an input. The gamma-variate model is supported by external classical references (Thompson 1964; Madsen 1992; Bindschadler et al.; Lin et al.), not by a self-citation chain. The only reference with overlapping authorship is [2] (Garcia-Villalba et al., which includes McVeigh, Kahn, and del Alamo), and it is cited merely as an example of prior CFD-based flow estimation in the introduction; it is not load-bearing for the present fitting or parameter equations, and no uniqueness theorem or ansatz is imported from the authors' prior work. Concerns that the 27-sample acquisition, aortic baseline subtraction, or extrapolated washout tail limit accuracy are real robustness/validation risks—the paper reports no comparison with 4D-flow MRI, TEE, CFD, or phantom data—but they are not circularity: the reported maps are exactly what the model produces from the data, and the paper does not claim an independent measurement of flow against which the fit is checked. Under the circularity rubric, the result is an internally consistent model-based estimate, not a reduction of a prediction to its own inputs.

Assumptions & free parameters 5 free parameters · 4 assumptions · 0 invented entities

The paper adds a per-voxel gamma-variate fit of an assumed contrast dilution curve; it postulates that the fitted curve and its first moment represent blood residence time, and relies on registration and baseline subtraction to make the signal meaningful. No new physical entities are postulated.

free parameters (5)
  • Per-voxel gamma-variate fit parameters (y_max, t_arrival, alpha) = Fit to each voxel; values not reported
    The central flow parameter maps are direct outputs of these per-voxel fits. They are estimated from data, with bounds chosen by the authors (Table I).
  • Cumulative signal threshold (100 HU) = 100 HU
    Voxels with cumulative sum across time below 100 HU are excluded from fitting; no sensitivity analysis or justification is given.
  • Gaussian smoothing sigma = 1 mm
    Applied to all timeframes before fitting; chosen by authors and not varied.
  • Baseline timepoints and descending aorta ROI = timepoints 1-3
    Baseline subtraction uses the mean of the first three timepoints in the descending aorta; this choice is applied to all LAA voxels and is not validated.
  • Gamma-variate fit parameter bounds = y_max: peak +/- 40 HU; t_arrival: peak time +/- 5 s; alpha: 0 to infinity
    Constraints chosen by the authors, not derived; they can bias parameter estimates in noisy voxels.
assumptions (4)
  • domain assumption The gamma-variate model describes the contrast concentration time course through the LAA.
    The model is standard for indicator dilution, but it is applied here without verification for LAA flow; Section II.C.
  • domain assumption CT intensity is linearly proportional to local contrast concentration after baseline subtraction.
    Needed for the fitted curve to represent contrast transit; stated implicitly in Section II.D but not calibrated.
  • domain assumption Nonrigid registration accurately aligns all 27 timeframes so each voxel tracks the same tissue location.
    Relies on GPNLPerf (Ref. 7) without assessing residual misregistration in the LAA; Section II.B.
  • domain assumption Blood residence time and arrival time derived from the gamma-variate fit correspond to actual blood flow properties.
    This is the interpretive bridge from curve fits to flow parameters; no comparison with flow measurements is provided.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Estimation of Blood Flow Parameters in the Left Atrial Appendage from 4DCT Dynamic Contrast Enhancement." pith.science (2026). https://pith.science/paper/WUMWLFPX

@misc{pith2026250202728,
  author       = {Pith},
  title        = {Pith review of: Estimation of Blood Flow Parameters in the Left Atrial Appendage from 4DCT Dynamic Contrast Enhancement},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/WUMWLFPX}},
  note         = {Machine review of arXiv:2502.02728}
}
read the original abstract

Cardiac CT is often used clinically in electrophysiology to evaluate cardiac morphology. One such case is to evaluate patients with Atrial Fibrillation (AF). AF can cause regions of slow blood flow and blood stasis through the left atrial appendage (LAA), and therefore, it may be preferable to evaluate blood flow through the LAA in addition to morphology. Although CT cannot measure flow directly, CT data has been used to estimate flow using modeling approaches such as Computational Fluid Dynamics, which take into account the cardiac geometry to simulate flow. Advances in CT technology now enable high-resolution imaging of the whole heart with low radiation doses. With multi-heartbeat imaging during a contrast injection, we can obtain 4-dimentional CT (4DCT) images to measure dynamic contrast enhancement directly. In this study, we use high-resolution 4DCT to acquire images of contrast enhancement across the LAA over multiple heartbeats. The CT contrast signal at each voxel over time is used to create dynamic contrast enhancement maps of parameters derived from a gamma-variate fit. These contrast enhancement maps enable quantification and visualization of spatial-temporal characteristics of flow parameters across the LAA.

Figures

Figures reproduced from arXiv: 2502.02728 by the authors.

Figure 1
Figure 1. Pipeline to generate dynamic contrast enhancement maps. The 4- dimensional CT (4DCT) image consists of 27 3D volumes, each collected during a 170 ms window around end-systole. In Step 1, the 27 timeframes are aligned using a nonrigid registration motion correction to produce registered 4DCT images. In Step 2, multiplanar reconstruction is used to create a 2- dimensional (2D) slice through the LAA, creating a set of … view at source ↗
Figure 2
Figure 2. Example gamma-variate fit and model parameters. CT signal contrast as a function of time (blue dots) is fit using the gamma-variate model (blue line). CT contrast is calculated as the CT image intensity at a given pixel minus a baseline signal. The baseline signal is calculated as the average CT signal at the first 3 timepoints in a region of interest in the descending aorta. The gamma￾variate fit parameters t$%"& a… view at source ↗
Figure 3
Figure 3. CT images of contrast enhancement over time for 2 patients. The post [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

11 extracted references · 11 canonical work pages

  1. [1]

    Computational Fluid Dynamic Analysis of the Left Atrial Appendage to Predict Thrombosis Risk,

    G. M. Bosi, A. Cook, R. Rai, L. J. Menezes, S. Schievano, R. Tori, G. Burriesci, “Computational Fluid Dynamic Analysis of the Left Atrial Appendage to Predict Thrombosis Risk,” Front. Cardiovasc. Med., vol. 5, Apr. 2018, DOI: https://doi.org/10.3389/fcvm.2018.00034

    work page arXiv 2018

  2. [2]

    Demonstration of Patient-Specific Simulations to Assess Left Atrial Appendage Thrombogenesis Risk,

    M. Garcia-Villalba, L. Rossini, A. Gonzalo, D. Vigneault, P. Martinez-Legazpi, E. Duran, O. Flores, J. Bermejo, E. McVeigh, A.M. Kahn, J.C. del Alamo, “Demonstration of Patient-Specific Simulations to Assess Left Atrial Appendage Thrombogenesis Risk,” Front. Physiol., vol.12, Feb. 2021, DOI: https://doi.org/10.3389/fphys.2021.596596

    work page arXiv 2021

  3. [7]

    GPNLPerf: Robust 4d Non-rigid Motion Correction for Myocardial Perfusion Analysis,

    S. Thiruvenkadam, K.S. Shriram, B. Patil, G. Nicholas, M. Teisseire, C. Cardon, J. Knoplioch, N. Subramanian, S. Kaushik, R. Mullick, “GPNLPerf: Robust 4d Non-rigid Motion Correction for Myocardial Perfusion Analysis,” Proc. of MICCAI 2016, pp. 255-263, Oct. 2016, DOI: https://doi.org/10.1007/978-3-319-46726-9_30

    work page doi:10.1007/978-3-319-46726-9_30 2016

  4. [10]

    CT perfusion: comparison of gamma-variate fit and deconvolution,

    Z.M. Lin S. Pohlman, A. Cook, S. Chandra, “CT perfusion: comparison of gamma-variate fit and deconvolution,” Proc. SPIE 4683, pp. 102-109, April 2002, DOI: 10.1117/12.463572

    work page doi:10.1117/12.463572 2002

  5. [11]

    Comparison of Blood Flow Models and Acquisitions for Quantitative Myocardial Perfusion Estimation from Dynamic CT,

    M. Bindschadler, D. Modgil, K. R. Branch, P.J. La Riviere, A.M. Alessio, “Comparison of Blood Flow Models and Acquisitions for Quantitative Myocardial Perfusion Estimation from Dynamic CT,” Phys. Med. Biol., vol. 59, no. 7, pp. 1533-1556. DOI: 10.1088/0031-9155/59/7/1533

    work page doi:10.1088/0031-9155/59/7/1533

  6. [1964]

    DOI: 10.1161/01.res.14.6.502

    work page doi:10.1161/01.res.14.6.502

  7. [1992]

    DOI: 10.1088/0031-9155/37/7/010

    work page doi:10.1088/0031-9155/37/7/010

  8. [1993]

    DOI: 10.1016/0002-9149(93)90737-w

    work page doi:10.1016/0002-9149(93)90737-w

Show all 11 references
  1. [1996]

    DOI: 10.1016/0003-4975(95)00887-X

    doi

  2. [2016]

    DOI: 10.1161/CIRCIMAGING.116.004984

    doi

  3. [2020]

    DOI: 10.1016/j.jcct.2019.09.003

    2019 doi

Pith tools

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