REVIEW 4 major objections 6 minor 18 references
Injection Bias Reduction Techniques in Quantitative Angiography Using Patient-Specific Phantoms of Intracranial Aneurysm
T0 review · 4 major / 6 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read By deconvolving the aneurysm dome time-density curve with the inlet arterial input function, SVD-based methods reduce the injection-induced variability in quantitative angiography parameters, cutting the mean transit time slope by 92–98%…
desk verdict Useful phantom dataset and a plausible effect, but the 'bias reduction' claim outruns the reported evidence — no ground truth, one geometry shown, no significance tests. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the impulse response function recovered by singular value decomposition of the arterial input function matrix. In standard SVD the inlet time-density curve is arranged as a Toeplitz matrix; block-circulant SVD and oscillation-index SVD use a block-circulant matrix, with the oscillation-index variant also computing an oscillation index to regularize the solution. Decomposing the AIF matrix into singular vectors and singular values, and then solving for the IRF, is what strips the injection waveform out of the aneurysm curve. The IRF then carries the hemodynamic information, and the QA parameters are read off its peak, its integral, and the ratio of integral to peak.
What would settle it
Repeat the same phantom experiment with bolus rates far outside the tested range (for example, a 2 ml bolus at 20 ml/s and a 20 ml bolus at 2 ml/s) and compare the recovered IRFs: if peak height, area, or mean transit time changes systematically with injection rate, the convolution model is not linear time-invariant and the reported invariance would not generalize beyond the tested conditions.
Extended reading notes
Core claim
The paper claims that the aneurysm dome time-density curve Q is the convolution of the inlet arterial input function Ca with a single impulse response function (IRF), and that SVD-based deconvolution inverts this relation. After deconvolution, the extracted IRF is convolved back with Ca to produce a reconvolved curve Qnew that closely tracks the measured aneurysm curve. Because the IRF is meant to encode only the aneurysm's transport behavior, its derived parameters—peak height (PHIRF), area under the curve (AUCIRF), and mean transit time (MTT)—become nearly independent of bolus volume and injection duration. The paper reports that the slope of MTT versus injection duration drops almost to zero for aneurysm geometry 1 with a 5 ml bolus, and presents this invariance as evidence that the SVD variants standardize quantitative angiography across aneurysm sizes, locations, and parent-artery conditions.
Load-bearing premise
The correction works only if one linear, time-invariant impulse response function fully describes contrast transport through the aneurysm regardless of injection rate, volume, or flow state, and only if the inlet time-density curve faithfully represents the contrast actually entering the aneurysm.
Editorial extensions
If this is right
- If SVD-deconvolved QA parameters are truly invariant to injection conditions, serial DSA runs in the same patient can be compared without matching injection volume or rate, making device-induced hemodynamic changes easier to track over time.
- The near-zero MTT slopes reported for the non-stenosed artery mean that after correction, the same aneurysm yields essentially the same transit time whether contrast is delivered over a quarter second or a full second.
- The fact that the correction also reduces MTT slope in the stenosed parent artery implies the method does not require normal flow to work, which matters for the diseased vessels where treatment decisions are made.
- Because the three SVD variants reduce variability by different amounts in the main non-stenosed case, the choice of variant changes the tightness of the resulting parameter invariance.
- The paper reports robustness across aneurysm sizes and locations, which would make the correction a general normalization step rather than a geometry-specific fix.
Reading between the lines
- Editorial inference: The linear time-invariant assumption has a testable consequence the paper does not state—if the recovered IRF is truly injection-independent, it should also remain unchanged under variations in pump stroke volume or heart rate, which alter the flow state rather than the injection.
- Editorial inference: A practical downstream gain is that SVD-normalized QA parameters could be pooled across clinical sites or compared across DSA systems that use different injector settings, since deconvolution is meant to absorb the injection waveform.
- Editorial inference: The different behavior in the stenosed artery, where standard SVD outperformed the other variants, hints that residual slope after deconvolution might itself encode information about the parent-artery flow condition; a future study could test whether the post-SVD slope separates stenosis severity from injection effects.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes SVD-based deconvolution methods (sSVD, bSVD, oSVD) as a means of reducing hand-injection variability in quantitative angiography (QA) of intracranial aneurysms. In an in-vitro patient-specific phantom setup with three aneurysm geometries at different Circle of Willis locations, the authors acquired DSA time-density curves for non-stenosed and stenosed parent arteries, with 5 ml and 10 ml boluses at four injection rates. They deconvolved the inlet TDC (AIF) from the dome TDC to obtain an impulse response function, from which they computed PH_IRF, AUC_IRF, and MTT. The central claim, stated in the Abstract, is that the SVD variants 'significantly reduced QA parameter variability due to injection techniques.' Detailed results are shown only for aneurysm geometry 1 with a 5 ml bolus, reporting slope reductions of MTT versus injection duration of 92.56% (sSVD), 98.01% (bSVD), and 98.45% (oSVD) for the non-stenosed artery, and corresponding reductions for the stenosed artery. The paper concludes that SVD-based deconvolution is a valid normalization method for QA analysis.
Significance. If the claims were fully supported, this would be a valuable contribution to neurointerventional imaging, because hand-injection variability is a known limitation of QA, and a demonstrably robust deconvolution correction would improve the clinical utility of DSA-derived hemodynamic parameters. The use of patient-specific phantoms with realistic flow waveforms, multiple injection protocols, and stenosed as well as non-stenosed conditions is a strength, and the three SVD variants are well-motivated from the perfusion-imaging literature. However, the evidence presented in the manuscript is substantially narrower than the abstract and conclusion claim: numerical results are given for only one of the three aneurysm geometries and only one bolus volume, no statistical significance tests support the word 'significantly,' and no independent ground truth is available to distinguish genuine bias reduction from regularization-induced shrinkage. The core idea remains plausible, but the current manuscript does not yet establish the claimed generality or the accuracy component of 'bias reduction.'
major comments (4)
- [Abstract and Results (Figure 3, Table 1)] The central claim of significant variability reduction 'across various aneurysm configurations' is not supported by the presented data. The Results section states that the experiment was repeated with three aneurysm geometries but 'we have provided the results for aneurysm geometry 1 ... having a 5ml bolus' only. No quantitative results, slopes, or parameter tables are shown for the other two geometries or for the 10 ml bolus condition. Since the abstract and conclusion generalize to different aneurysm sizes, locations, and injection volumes, the missing results are load-bearing; the manuscript should either include the full dataset or explicitly limit the claims to the conditions actually reported.
- [Results (Figure 3, Table 1)] The word 'significantly' in the Abstract and Conclusion is not backed by any statistical test. The table reports means and standard deviations for three repeats, but there is no t-test, ANOVA, confidence interval for the slopes, or comparison of the pre- and post-deconvolution slopes. The percentage slope reductions could be tested against the null hypothesis of no reduction, and the variability across injection rates should be compared in a formal way. As written, 'significantly' has no statistical meaning, and the claim should be reworded or supported with appropriate tests.
- [Section 2.2, Eq. (deconvolution model)] Under the assumed linear time-invariant (LTI) convolution model, the deconvolution step removes injection-duration dependence by construction: if the dome TDC is exactly the convolution of the inlet AIF with a fixed impulse response, then the recovered IRF and its MTT are independent of the input function. The observed slope reduction therefore tests consistency with LTI behavior, not the accuracy of the recovered QA parameters. The large discrepancy between raw MTT values (about 1.9-2.6 s) and oSVD-derived MTT values (about 0.08-0.11 s) in Table 1 is not discussed, and without an independent reference (e.g., a known IRF, flow measurements, or CFD-computed TDCs) it is unclear whether the near-flat slopes reflect true physiological transit times or an artifact of regularization that shrinks IRF estimates toward a near-constant value. The manuscript should address this ambiguity and temper the use of the term 'bias reduction' unless accuracy against a ground truth is demonstrated.
- [Section 2.2, Methods (regularization parameters)] The manuscript does not disclose the Tikhonov regularization parameter used in sSVD, the truncation threshold or other regularization settings for bSVD, or the oscillation-index regularization parameter and penalty strength for oSVD. These settings control the bias-variance tradeoff and directly influence the shape of the recovered IRF and hence MTT, PH_IRF, and AUC_IRF. Without these values, the results are not reproducible, and the reader cannot assess whether the reported slope reductions are sensitive to arbitrarily chosen regularization strengths. The authors should state the exact parameter values and, if possible, show sensitivity of the main results to those parameters.
minor comments (6)
- [Keywords] The keyword 'Single Value Decomposition' should be corrected to 'Singular Value Decomposition.'
- [Section 2.2] The phrase 'AUCIRF which is calculated by integrating the IRF with respect to AIF' is ambiguous; it should be 'integrating the IRF over time' or 'with respect to time,' since the area under the IRF curve is a time integral.
- [Section 2.2 and Figure 2] The reconstructed TDC Qnew is mentioned but never used in the analysis or validation; if Qnew is intended as a quality check of the deconvolution, the authors should report the residuals, correlation, or another fit metric.
- [Table 1] The table formatting is difficult to follow: the row and column alignment for DIFFSL and AVGSL is inconsistent, and some cells appear misplaced or missing. The caption also introduces SLOPESVD, SLOPE, DIFFSL, and AVGSL, but the table does not clearly show all of these quantities for every row. A clearer layout or a separate table for slope statistics would improve readability.
- [Figure 3] The axis labels and symbols in Figure 3 are not fully legible in the provided version; the authors should ensure that the units of MTT are explicitly stated on the y-axis and that stenosed versus non-stenosed curves are clearly distinguished by markers or legends.
- [References] Reference [17], 'A Survey of Singular Value Decomposition Methods for Distributed Tall/Skinny Data,' appears to be a computing-systems reference that is not directly relevant to the medical-imaging context; a standard numerical linear algebra reference (e.g., Golub and Van Loan) would be more appropriate for the SVD description.
Circularity Check
No significant circularity: the SVD deconvolution is a standard inversion applied to independent phantom measurements, and the invariance claim is an empirical result; the main limitation is absence of an external ground truth, which is a validity concern, not a circular one.
full rationale
The paper's derivation chain is not circular. Phantom DSA images are acquired under varied injection volumes and durations; inlet and dome time-density curves (TDCs) are measured directly from those images; a Toeplitz or block-circulant AIF matrix is built from the inlet TDC; SVD with regularization is applied to recover an impulse response function (IRF); QA parameters (PH_IRF, AUC_IRF, MTT) are computed from that IRF; and the reported improvement is the reduction in slope of MTT versus injection duration before and after deconvolution. This is a standard deconvolution pipeline, not a fit to the reported outcome. The load-bearing model assumption, Q = Ca * IRF, is stated explicitly in Section 2.2, and the experiment could in principle falsify it: if the system were nonlinear or flow-dependent, the recovered IRF and MTT_SVD would retain injection-duration dependence. The fact that the result aligns with the LTI model is an empirical finding from independent phantom data, not a tautology. The paper contains minor self-citations (e.g., refs. [6], [15], and [16]) that establish prior context and experimental methodology, but the central claim does not reduce to those citations. The skeptic's concern—that slope reduction could reflect regularization-induced shrinkage rather than true bias removal—is a legitimate validity limitation because no independent transit-time gold standard is used, but that is an evidentiary weakness, not circularity. No step was found where a prediction is equivalent to its input by construction or where a fitted parameter is renamed as a prediction.
Assumptions & free parameters
free parameters (1)
- Tikhonov regularization parameter (lambda) =
not reported
assumptions (3)
- domain assumption Convolution model Q(t) = Ca(t) * IRF(t) holds for the aneurysm compartment
- standard math Toeplitz and block-circulant matrix constructions implement the convolution correctly
- domain assumption The in vitro flow loop reproduces relevant hemodynamics of intracranial aneurysms
Cite this review
Pith. "Pith review of Injection Bias Reduction Techniques in Quantitative Angiography Using Patient-Specific Phantoms of Intracranial Aneurysm." pith.science (2026). https://pith.science/paper/T2KEGJSL
@misc{pith2026241114475,
author = {Pith},
title = {Pith review of: Injection Bias Reduction Techniques in Quantitative Angiography Using Patient-Specific Phantoms of Intracranial Aneurysm},
year = {2026},
howpublished = {\url{https://pith.science/paper/T2KEGJSL}},
note = {Machine review of arXiv:2411.14475}
}
read the original abstract
In intracranial aneurysm (IA) treatment, digital subtraction angiography (DSA) monitors device-induced hemodynamic changes. Quantitative angiography (QA) provides more precise assessments but is limited by hand-injection variability. This study evaluates correction methods using in vitro phantoms that mimic diverse aneurysm morphologies and locations, addressing the 2D and temporal limitations of DSA. We used a patient-specific phantom to replicate three distinct IA morphologies at various Circle of Willis points: the middle cerebral artery (MCA), anterior communicating artery (ACA), and the internal carotid artery (ICA), each varying in size and shape. The diameters of the IA at MCA, ACA and ICA are 10.1, 10 and 7 millimeters, respectively. QA parameters for both non-stenosed and stenosed conditions were measured with 5ml and 10ml boluses over various injection durations to generate time density curves (TDCs). To address the variability in injection, several singular value decomposition (SVD) variants, standard SVD (sSVD) with Tikhonov regularization, block-circulant SVD (bSVD), and oscillation index SVD (oSVD) were applied. These methods enabled the extraction of IA impulse response function (IRF), peak height (PHIRF), area under the curve (AUCIRF), and mean transit time (MTT). We evaluated the robustness of bias-reducing methods by observing the invariance of these parameters with respect to the injection conditions, and the location and size of the aneurysm. The application of SVD variants, sSVD, bSVD, and oSVD, significantly reduced QA parameter variability due to injection techniques.
Reference graph
Works this paper leans on
-
[6]
P. Mondal, S. V. Setlur Nagesh, S. Sommers-Thaler et al., “Effect of singular value decomposition on removing injection variability in 2D quantitative angiography: An in silico and in vitro phantoms study,” Med Phys, (2024)
work page 2024
-
[1]
A. M. Jersey, and D. M. Foster, [Cerebral Aneurysm], Treasure Island (FL)(2023)
work page 2023
-
[2]
Brain Vascular Imaging Techniques,
B. Lavina, “Brain Vascular Imaging Techniques,” Int J Mol Sci, 18(1), (2016)
work page 2016
-
[3]
Digital subtraction angiography (DSA): a neurosurgical perspective,
P. Butler, “Digital subtraction angiography (DSA): a neurosurgical perspective,” Br J Neurosurg , 1(3), 323-33 (1987)
work page 1987
-
[4]
C. N. Ionita, H. Suri, S. Nataranjian et al., “Angiographic imaging evaluation of patient -specific bifurcation-aneurysm phantom treatment with pre-shaped, self-expanding, flow-diverting stents: feasibility study,” Proc SPIE Int Soc Opt Eng, 7965, 79651H1-79651H9 (2011)
work page 2011
-
[5]
The principle of digital subtraction angiography and radiological protection,
K. Okamoto, J. Ito, K. Sakai et al., “The principle of digital subtraction angiography and radiological protection,” Interventional Neuroradiology, 6, 25-31 (2000)
work page 2000
- [7]
-
[8]
R. A. Rava, A. B. Allman, S. Rudin et al., “Effect of Truncated Singular Value Decomposition on Digital Subtraction Angiography Derived Angiographic Parametric Imaging Maps,” Medical Imaging 2020: Physics of Medical Imaging, 11312, (2020)
work page 2020
Show all 18 references
-
[9]
Singular value decomposition as a tool for background corrections in time -resolved XFEL scattering data,
K. Haldrup, “Singular value decomposition as a tool for background corrections in time -resolved XFEL scattering data,” Philosophical Transactions of the Royal Society B -Biological Sciences , 369(1647), (2014)
2014
-
[10]
On Tikhonov regularization, bias and variance in nonlinear system identification,
T. A. Johansen, “On Tikhonov regularization, bias and variance in nonlinear system identification,” Automatica, 33(3), 441-446 (1997)
1997
-
[11]
On singular values of block circulant matrices,
M. Niezgoda, “On singular values of block circulant matrices,” Linear Algebra and Its Applications, 563, 400-410 (2019)
2019
-
[12]
Tracer ar rival timing -insensitive technique for estimating flow in MR perfusion -weighted imaging using singular value decomposition with a block-circulant deconvolution matrix,
O. Wu, L. Ostergaard, R. M. Weisskoff et al. , “Tracer ar rival timing -insensitive technique for estimating flow in MR perfusion -weighted imaging using singular value decomposition with a block-circulant deconvolution matrix,” Magnetic Resonance in Medicine, 50(1), 164-174 (2003)
2003
-
[13]
Optimal Sparse Singular Value Decomposition for High-Dimensional High-Order Data,
A. R. Zhang, and R. G. Han, “Optimal Sparse Singular Value Decomposition for High-Dimensional High-Order Data,” Journal of the American Statistical Association, 114(528), 1708-1725 (2019)
2019
-
[14]
Initial evaluation of 2D and 3D simulated high-speed 1000 fps vascular contrast -flow image sequences using computational fluid dynamics (CFD),
A. Shields, K. Williams, S. S. Veeturi et al., “Initial evaluation of 2D and 3D simulated high-speed 1000 fps vascular contrast -flow image sequences using computational fluid dynamics (CFD),” Medical Imaging 2022: Biomedical Applications in Molecular, Structural, and Function...
2022
-
[15]
Mondal, A
P. Mondal, A. Shields, M. M. S. Bhurw ani et al., [In-silico analysis of curve fitting in angiographic parametric imaging in intracranial aneurysms to reduce patient and interventionalist -induced errors] SPIE, MI (2024)
2024
-
[16]
Effect of inject ion technique on temporal parametric imaging derived from digital subtraction angiography in patient specific phantoms,
C. N. Ionita, V. L. Garcia, D. R. Bednarek et al. , “Effect of inject ion technique on temporal parametric imaging derived from digital subtraction angiography in patient specific phantoms,” Medical Imaging 2014: Biomedical Applications in Molecular, Structural, and Functional...
2014
-
[17]
A Survey of Singular Value Decomposition Methods for Distributed Tall/Skinny Data,
D. Schmidt, “A Survey of Singular Value Decomposition Methods for Distributed Tall/Skinny Data,” Proceedings of Scala 2020: 11th Workshop on Latest Advances in Scalable Algorithms for Large - Scale Systems, 27-34 (2020)
2020
-
[18]
Application of singular value decomposition to the analysis of time-resolved macromolecular X-ray data,
M. Schmidt, S. Rajagopal, Z. Ren et al., “Application of singular value decomposition to the analysis of time-resolved macromolecular X-ray data,” Biophysical Journal, 84(3), 2112-2129 (2003)
2003
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.