Differentiable Deconvolution for Improved Stroke Perfusion Analysis

Perfusion imaging is the current gold standard for acute ischemic stroke analysis. It allows quantification of the salvageable and non-salvageable tissue regions (penumbra and core areas respectively). In clinical settings, the singular value decomposition (SVD) deconvolution is one of the most accepted and used approaches for generating interpretable and physically meaningful maps. Though this method has been widely validated in experimental and clinical settings, it might produce suboptimal results because the chosen inputs to the model cannot guarantee optimal performance. For the most critical input, the arterial input function (AIF), it is still controversial how and where it should be chosen even though the method is very sensitive to this input. In this work we propose an AIF selection approach that is optimized for maximal core lesion segmentation performance. The AIF is regressed by a neural network optimized through a differentiable SVD deconvolution, aiming to maximize core lesion segmentation agreement with ground truth data. To our knowledge, this is the first work exploiting a differentiable deconvolution model with neural networks. We show that our approach is able to generate AIFs without any manual annotation, and hence avoiding manual rater’s influences. The method achieves manual expert performance in the ISLES18 dataset. We conclude that the methodology opens new possibilities for improving perfusion imaging quantification with deep neural networks.

[1]  Paul Suetens,et al.  Perfusion parameter estimation using neural networks and data augmentation , 2018, BrainLes@MICCAI.

[2]  Cristian Sminchisescu,et al.  Training Deep Networks with Structured Layers by Matrix Backpropagation , 2015, ArXiv.

[3]  Nitish Srivastava,et al.  Dropout: a simple way to prevent neural networks from overfitting , 2014, J. Mach. Learn. Res..

[4]  K Scheffler,et al.  Analysis of input functions from different arterial branches with gamma variate functions and cluster analysis for quantitative blood volume measurements. , 2000, Magnetic resonance imaging.

[5]  Cristian Sminchisescu,et al.  Matrix Backpropagation for Deep Networks with Structured Layers , 2015, 2015 IEEE International Conference on Computer Vision (ICCV).

[6]  F. Calamante Arterial input function in perfusion MRI: a comprehensive review. , 2013, Progress in nuclear magnetic resonance spectroscopy.

[7]  K. Murase,et al.  Determination of arterial input function using fuzzy clustering for quantification of cerebral blood flow with dynamic susceptibility contrast‐enhanced MR imaging , 2001, Journal of magnetic resonance imaging : JMRI.

[8]  Manolis I. A. Lourakis,et al.  Estimating the Jacobian of the Singular Value Decomposition: Theory and Applications , 2000, ECCV.

[9]  Geoffrey E. Hinton,et al.  ImageNet classification with deep convolutional neural networks , 2012, Commun. ACM.

[10]  T-Y Lee,et al.  Serial changes in CT cerebral blood volume and flow after 4 hours of middle cerebral occlusion in an animal model of embolic cerebral ischemia. , 2007, AJNR. American journal of neuroradiology.

[11]  Søren Christensen,et al.  Automatic selection of arterial input function using cluster analysis , 2006, Magnetic resonance in medicine.

[12]  Alessandra Bertoldo,et al.  Automatic selection of arterial input function on dynamic contrast-enhanced MR images , 2011, Comput. Methods Programs Biomed..

[13]  M. Reiser,et al.  Deconvolution of bolus-tracking data: a comparison of discretization methods , 2007, Physics in medicine and biology.

[14]  et al.,et al.  ISLES 2015 ‐ A public evaluation benchmark for ischemic stroke lesion segmentation from multispectral MRI , 2017, Medical Image Anal..

[15]  Lin Shi,et al.  Automatic detection of arterial input function in dynamic contrast enhanced MRI based on affinity propagation clustering , 2014, Journal of magnetic resonance imaging : JMRI.

[16]  D. Altman,et al.  STATISTICAL METHODS FOR ASSESSING AGREEMENT BETWEEN TWO METHODS OF CLINICAL MEASUREMENT , 1986, The Lancet.

[17]  Rebecca Fahrig,et al.  Deconvolution-Based CT and MR Brain Perfusion Measurement: Theoretical Model Revisited and Practical Implementation Details , 2011, Int. J. Biomed. Imaging.

[18]  Paul Suetens,et al.  Optimization with soft Dice can lead to a volumetric bias , 2019, BrainLes@MICCAI.

[19]  Qi Yang,et al.  An Automatic Estimation of Arterial Input Function Based on Multi-Stream 3D CNN , 2019, Front. Neuroinform..

[20]  Matus Straka,et al.  A benchmarking tool to evaluate computer tomography perfusion infarct core predictions against a DWI standard , 2016, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[21]  Michael Kistler,et al.  The Virtual Skeleton Database: An Open Access Repository for Biomedical Research and Collaboration , 2013, Journal of medical Internet research.

[22]  M. Moseley,et al.  Automated method for generating the arterial input function on perfusion-weighted MR imaging: validation in patients with stroke. , 2005, AJNR. American journal of neuroradiology.

[23]  M. Wintermark,et al.  Automated CT perfusion imaging for acute ischemic stroke , 2019, Neurology.

[24]  Mark W Parsons,et al.  Whole-Brain CT Perfusion to Quantify Acute Ischemic Penumbra and Core. , 2016, Radiology.

[25]  Roland Bammer,et al.  Ischemic core and hypoperfusion volumes predict infarct size in SWIFT PRIME , 2016, Annals of neurology.