Intelligent Imaging of Perfusion Using Arterial Spin Labelling

Arterial spin labelling (ASL) is a powerful magnetic resonance imaging technique, which can be used to noninvasively measure perfusion in the brain and other organs of the body. Promising research results show how ASL might be used in stroke, tumours, dementia and paediatric medicine, in addition to many other areas. However, significant obstacles remain to prevent widespread use: ASL images have an inherently low signal to noise ratio, and are susceptible to corrupting artifacts from motion and other sources. The objective of the work in this thesis is to move towards an "intelligent imaging" paradigm: one in which the image acquisition, reconstruction and processing are mutually coupled, and tailored to the individual patient. This thesis explores how ASL images may be improved at several stages of the imaging pipeline. We review the relevant ASL literature, exploring details of ASL acquisitions, parameter inference and artifact post-processing. We subsequently present original work: we use the framework of Bayesian experimental design to generate optimised ASL acquisitions, we present original methods to improve parameter inference through anatomically-driven modelling of spatial correlation, and we describe a novel deep learning approach for simultaneous denoising and artifact filtering. Using a mixture of theoretical derivation, simulation results and imaging experiments, the work in this thesis presents several new approaches for ASL, and hopefully will shape future research and future ASL usage

ExploreASL: An image processing pipeline for multi-center ASL perfusion MRI studies

Arterial spin labeling (ASL) has undergone significant development since its inception, with a focus on improving standardization and reproducibility of its acquisition and quantification. In a community-wide effort towards robust and reproducible clinical ASL image processing, we developed the software package ExploreASL, allowing standardized analyses across centers and scanners. The procedures used in ExploreASL capitalize on published image processing advancements and address the challenges of multi-center datasets with scanner-specific processing and artifact reduction to limit patient exclusion. ExploreASL is self-contained, written in MATLAB and based on Statistical Parameter Mapping (SPM) and runs on multiple operating systems. To facilitate collaboration and data-exchange, the toolbox follows several standards and recommendations for data structure, provenance, and best analysis practice. ExploreASL was iteratively refined and tested in the analysis of >10,000 ASL scans using different pulse-sequences in a variety of clinical populations, resulting in four processing modules: Import, Structural, ASL, and Population that perform tasks, respectively, for data curation, structural and ASL image processing and quality control, and finally preparing the results for statistical analyses on both single-subject and group level. We illustrate ExploreASL processing results from three cohorts: perinatally HIV-infected children, healthy adults, and elderly at risk for neurodegenerative disease. We show the reproducibility for each cohort when processed at different centers with different operating systems and MATLAB versions, and its effects on the quantification of gray matter cerebral blood flow. ExploreASL facilitates the standardization of image processing and quality control, allowing the pooling of cohorts which may increase statistical power and discover between-group perfusion differences. Ultimately, this workflow may advance ASL for wider adoption in clinical studies, trials, and practice

Measurement of vascular water transport in human subjects using time-resolved pulsed arterial spin labelling.

Most approaches to arterial spin labelling (ASL) data analysis aim to provide a quantitative measure of the cerebral blood flow (CBF). This study, however, focuses on the measurement of the transfer time of blood water through the capillaries to the parenchyma (referred to as the capillary transfer time, CTT) as an alternative parameter to characterise the haemodynamics of the system. The method employed is based on a non-compartmental model, and no measurements need to be added to a common time-resolved ASL experiment. Brownian motion of labelled spins in a potential was described by a one-dimensional general Langevin equation as the starting point, and as a Fokker-Planck differential equation for the averaged distribution of labelled spins at the end point, which takes into account the effects of flow and dispersion of labelled water by the pseudorandom nature of the microvasculature and the transcapillary permeability. Multi-inversion time (multi-TI) ASL data were acquired in 14 healthy subjects on two occasions in a test-retest design, using a pulsed ASL sequence and three-dimensional gradient and spin echo (3D-GRASE) readout. Based on an error analysis to predict the size of a region of interest (ROI) required to obtain reasonably precise parameter estimates, data were analysed in two relatively large ROIs, i.e. the occipital lobe (OC) and the insular cortex (IC). The average values of CTT in OC were 260 ± 60 ms in the first experiment and 270 ± 60 ms in the second experiment. The corresponding IC values were 460 ± 130 ms and 420 ± 139 ms, respectively. Information related to the water transfer time may be important for diagnostics and follow-up of cerebral conditions or diseases characterised by a disrupted blood-brain barrier or disturbed capillary blood flow. Copyright © 2015 John Wiley & Sons, Ltd

ICA-based denoising for ASL perfusion imaging

Arterial Spin Labelling (ASL) imaging derives a perfusion image by tracing the accumulation of magnetically labeled blood water in the brain. As the image generated has an intrinsically low signal to noise ratio (SNR), multiple measurements are routinely acquired and averaged, at a penalty of increased scan duration and opportunity for motion artefact. However, this strategy alone might be ineffective in clinical settings where the time available for acquisition is limited and patient motion are increased. This study investigates the use of an Independent Component Analysis (ICA) approach for denoising ASL data, and its potential for automation.72 ASL datasets (pseudo-continuous ASL; 5 different post-labeling delays: 400, 800, 1200, 1600, 2000 m s; total volumes = 60) were collected from thirty consecutive acute stroke patients. The effects of ICA-based denoising (manual and automated) where compared to two different denoising approaches, aCompCor, a Principal Component-based method, and Enhancement of Automated Blood Flow Estimates (ENABLE), an algorithm based on the removal of corrupted volumes. Multiple metrics were used to assess the changes in the quality of the data following denoising, including changes in cerebral blood flow (CBF) and arterial transit time (ATT), SNR, and repeatability. Additionally, the relationship between SNR and number of repetitions acquired was estimated before and after denoising the data.The use of an ICA-based denoising approach resulted in significantly higher mean CBF and ATT values (p [less than] 0.001), lower CBF and ATT variance (p [less than] 0.001), increased SNR (p [less than] 0.001), and improved repeatability (p [less than] 0.05) when compared to the raw data. The performance of manual and automated ICA-based denoising was comparable. These results went beyond the effects of aCompCor or ENABLE. Following ICA-based denoising, the SNR was higher using only 50% of the ASL-dataset collected than when using the whole raw data.The results show that ICA can be used to separate signal from noise in ASL data, improving the quality of the data collected. In fact, this study suggests that the acquisition time could be reduced by 50% without penalty to data quality, something that merits further study. Independent component classification and regression can be carried out either manually, following simple criteria, or automatically

Combined Denoising and Suppression of Transient Artifacts in Arterial Spin Labeling MRI Using Deep Learning

Background: Arterial spin labeling (ASL) is a useful tool for measuring cerebral blood flow (CBF). However, due to the low signal-to-noise ratio (SNR) of the technique, multiple repetitions are required, which results in prolonged scan times and increased susceptibility to artifacts. Purpose: To develop a deep-learning-based algorithm for simultaneous denoising and suppression of transient artifacts in ASL images. Study Type: Retrospective. Subjects: 131 pediatric neuro-oncology patients for model training and 11 healthy adult subjects for model evaluation. Field Strength/Sequence: 3T / pseudo-continuous and pulsed ASL with 3D gradient-and-spin-echo readout. Assessment: A denoising autoencoder (DAE) model was designed with stacked encoding/decoding convolutional layers. Reference standard images were generated by averaging 10 pairwise ASL subtraction images. The model was trained to produce perfusion images of a similar quality using a single subtraction image. Performance was compared against Gaussian and non-local means (NLM) filters. Evaluation metrics included SNR, peak SNR (PSNR), and structural similarity index (SSIM) of the CBF images, compared to the reference standard. Statistical Tests: One-way analysis of variance (ANOVA) tests for group comparisons. Results: The DAE model was the only model to produce a significant increase in SNR compared to the raw images (P < 0.05), providing an average SNR gain of 62%. The DAE model was also effective at suppressing transient artifacts, and was the only model to show a significant improvement in accuracy in the generated CBF images, as assessed using PSNR values (P < 0.05). In addition, using data from multiple inflow time acquisitions, the DAE images produced the best fit to the Buxton kinetic model, offering a 75% reduction in the fitting error compared to the raw images. Data Conclusion: Deep-learning-based algorithms provide superior accuracy when denoising ASL images, due to their ability to simultaneously increase SNR and suppress artifactual signals in raw ASL images. Level of Evidence: 3. Technical Efficacy Stage: 1

Improved quantification of perfusion in patients with cerebrovascular disease.

In recent years measurements of cerebral perfusion using bolus-tracking MRI have become common clinical practice in the diagnosis and management of patients with stroke and cerebrovascular disease. An active area of research is the development of methods to identify brain tissue that is at risk of irreversible damage, but amenable to salvage using reperfusion treatments, such as thrombolysis. However, the specificity and sensitivity of these methods are limited by the inaccuracies in the perfusion data. Accurate measurements of perfusion are difficult to obtain, especially in patients with cerebrovascular diseases. In particular, if the bolus of MR contrast is delayed and/or dispersed due to cerebral arterial abnormalities, perfusion is likely to be underestimated using the standard analysis techniques. The potential for such underestimation is often overlooked when using the perfusion maps to assess stroke patients. Since thrombolysis can increase the risk of haemorrhage, a misidentification of 'at-risk' tissue has potentially dangerous clinical implications. This thesis presents several methodologies which aim to improve the accuracy and interpretation of the analysed bolus-tracking data. Two novel data analysis techniques are proposed, which enable the identification of brain regions where delay and dispersion of the bolus are likely to bias the perfusion measurements. In this way true hypoperfusion can be distinguished from erroneously low perfusion estimates. The size of the perfusion measurement errors are investigated in vivo, and a parameterised characterisation of the bolus delay and dispersion is obtained. Such information is valuable for the interpretation of in vivo data, and for further investigation into the effects of abnormal vasculature on perfusion estimates. Finally, methodology is presented to minimise the perfusion measurement errors prevalent in patients with cerebrovascular diseases. The in vivo application of this method highlights the dangers of interpreting perfusion values independently of the bolus delay and dispersion

Learning-based Algorithms for Inverse Problems in MR Image Reconstruction and Quantitative Perfusion Imaging

Medical imaging has become an integral part of the clinical pipeline through its widespread use in the diagnosis, prognosis and treatment planning of several diseases. Magnetic Resonance Imaging (MRI) is particularly useful because it is free from ionizing radiation and is able to provide excellent soft tissue contrast. However, MRI suffers from drawbacks like long scanning durations that increase the cost of imaging and render the acquired images vulnerable to artifacts like motion. In modalities like Arterial Spin Labeling (ASL), which is used for non-invasive and quantitative perfusion imaging, low signal-to-noise ratio and lack of precision in parameter estimates also present significant problems. In this thesis, we develop and present algorithms whose focus can be divided into two broad categories. First, we investigate the reconstruction of MR images from fewer measurements, using data-driven machine learning to fill in the gaps in acquisition, thereby reducing the scan duration. Specifically, we first combine a supervised and an unsupervised (blind) learned dictionary in a residual fashion as a spatial prior in MR image reconstruction, and then extend this framework to include deep supervised learning. The latter, called blind primed supervised (BLIPS) learning, proposes that there exists synergy between features learned using shallower dictionary-based methods or traditional prior-based image reconstruction and those learned using newer deep supervised learning-based approaches. We show that this synergy can be exploited to yield reconstructions that are approx. 0.5-1 dB better in PSNR (in avg. across undersampling patterns). We also observe that the BLIPS algorithm is more robust to a scarcity of available training data, yielding reconstructions that are approx. 0.8 dB better (in terms of avg. PSNR) compared to strict supervised learning reconstruction when training data is very limited. Secondly, we aim to provide more precise estimates for multiple physiological parameters and tissue properties from ASL scans by estimation-theory-based optimization of ASL scan design, and combination with MR Fingerprinting. For this purpose, we use the Cramer-Rao Lower Bound (CRLB) for optimizing the scan design, and deep learning for regression-based estimation. We also show that regardless of the estimator used, optimization improves the precision in parameter estimates, and enables us to increase the available ‘useful’ information obtained in a fixed scanning duration. Specifically, we successfully improve the theoretical precision of perfusion estimates by 4.6% compared to a scan design where the repetition times are randomly chosen (a popular choice in literature) thereby yielding a 35.2% improvement in the corresponding RMSE in our in-silico experiments. This improvement is also visually evident in our in-vivo studies on healthy human subjects.PHDElectrical and Computer EngineeringUniversity of Michigan, Horace H. Rackham School of Graduate Studieshttp://deepblue.lib.umich.edu/bitstream/2027.42/169819/1/anishl_1.pd

ExploreASL: an image processing pipeline for multi-center ASL perfusion MRI studies

Arterial spin labeling (ASL) has undergone significant development since its inception, with a focus on improving standardization and reproducibility of its acquisition and quantification. In a community-wide effort towards robust and reproducible clinical ASL image processing, we developed the software package ExploreASL, allowing standardized analyses across centers and scanners. The procedures used in ExploreASL capitalize on published image processing advancements and address the challenges of multi-center datasets with scanner-specific processing and artifact reduction to limit patient exclusion. ExploreASL is self-contained, written in MATLAB and based on Statistical Parameter Mapping (SPM) and runs on multiple operating systems. To facilitate collaboration and data-exchange, the toolbox follows several standards and recommendations for data structure, provenance, and best analysis practice. ExploreASL was iteratively refined and tested in the analysis of >10,000 ASL scans using different pulse-sequences in a variety of clinical populations, resulting in four processing modules: Import, Structural, ASL, and Population that perform tasks, respectively, for data curation, structural and ASL image processing and quality control, and finally preparing the results for statistical analyses on both single-subject and group level. We illustrate ExploreASL processing results from three cohorts: perinatally HIV-infected children, healthy adults, and elderly at risk for neurodegenerative disease. We show the reproducibility for each cohort when processed at different centers with different operating systems and MATLAB versions, and its effects on the quantification of gray matter cerebral blood flow. ExploreASL facilitates the standardization of image processing and quality control, allowing the pooling of cohorts which may increase statistical power and discover between-group perfusion differences. Ultimately, this workflow may advance ASL for wider adoption in clinical studies, trials, and practice