Nonlinear stochastic regularization to characterize tissue residue function in bolus-tracking MRI: assessment and comparison with SVD, block-circulant SVD, and Tikhonov

Assessment of Perfusion Volumes by a New Automated Software for Computed Tomography Perfusion

INTRODUCTION

To compare the perfusion volumes assessed by a new automated CT perfusion (CTP) software iStroke with the circular singular value decomposition software RAPID and determine its predictive value for functional outcome in patients with acute ischaemic stroke (AIS) who underwent endovascular treatment (EVT).

METHODS

Data on patients with AIS were collected from four hospitals in China. All patients received CTP followed by EVT with complete recanalisation within 24 hours of symptom onset. We evaluated the agreement of CTP measures between the two softwares by Spearman's rank correlation tests and kappa tests. Bland-Altman plots were used to evaluate the agreement of infarct core volume (ICV) on CTP and ground truth on diffusion-weighted imaging (DWI). Logistic regression models were used to test the association between ICV on these two softwares and functional outcomes.

RESULTS

Among 326 patients, 228 had DWI examinations and 40 of them had infarct volume >70 mL. In all patients, the infarct core and hypoperfusion volumes on iStroke had a strong correlation with those on RAPID (ρ=0.68 and 0.66, respectively). The agreement of large infarct core (volume >70 mL) was substantial (kappa=0.73, p<0.001) between these two softwares. The ICV measured by iStroke and RAPID was significantly correlated with independent functional outcome at 90 days (p=0.009 and p<0.001, respectively). In patients with DWI examinations and those with an ICV >70 mL, the ICV of iStroke and RAPID was comparable on individual agreement with ground truth.

CONCLUSION

The automatic CTP software iStroke is a reliable tool for assessing infarct core and mismatch volumes, making it clinically useful for selecting patients with AIS for acute reperfusion therapy in the extended time window.

suMRak: a multi-tool solution for preclinical brain MRI data analysis

Introduction

Magnetic resonance imaging (MRI) is invaluable for understanding brain disorders, but data complexity poses a challenge in experimental research. In this study, we introduce suMRak, a MATLAB application designed for efficient preclinical brain MRI analysis. SuMRak integrates brain segmentation, volumetry, image registration, and parameter map generation into a unified interface, thereby reducing the number of separate tools that researchers may require for straightforward data handling.

Methods and implementation

All functionalities of suMRak are implemented using the MATLAB App Designer and the MATLAB-integrated Python engine. A total of six helper applications were developed alongside the main suMRak interface to allow for a cohesive and streamlined workflow. The brain segmentation strategy was validated by comparing suMRak against manual segmentation and ITK-SNAP, a popular open-source application for biomedical image segmentation.

Results

When compared with the manual segmentation of coronal mouse brain slices, suMRak achieved a high Sørensen-Dice similarity coefficient (0.98 ± 0.01), approaching manual accuracy. Additionally, suMRak exhibited significant improvement (p = 0.03) when compared to ITK-SNAP, particularly for caudally located brain slices. Furthermore, suMRak was capable of effectively analyzing preclinical MRI data obtained in our own studies. Most notably, the results of brain perfusion map registration to T2-weighted images were shown, improving the topographic connection to anatomical areas and enabling further data analysis to better account for the inherent spatial distortions of echoplanar imaging.

Discussion

SuMRak offers efficient MRI data processing of preclinical brain images, enabling researchers' consistency and precision. Notably, the accelerated brain segmentation, achieved through K-means clustering and morphological operations, significantly reduces processing time and allows for easier handling of larger datasets.

Multiparametric MR Evaluation of the Photoperiodic Regulation of Hypothalamic Structures in Sheep

Most organisms on earth, humans included, have developed strategies to cope with environmental day-night and seasonal cycles to survive. For most of them, their physiological and behavioral functions, including the reproductive function, are synchronized with the annual changes of day length, to ensure winter survival and subsequent reproductive success in the following spring. Sheep are sensitive to photoperiod, which also regulates natural adult neurogenesis in their hypothalamus. We postulate that the ovine model represents a good alternative to study the functional and metabolic changes occurring in response to photoperiodic changes in hypothalamic structures of the brain. Here, the impact of the photoperiod on the neurovascular coupling and the metabolism of the hypothalamic structures was investigated at 3T using BOLD fMRI, perfusion-MRI and proton magnetic resonance spectroscopy (1H-MRS). A longitudinal study involving 8 ewes was conducted during long days (LD) and short days (SD) revealing significant BOLD, rCBV and metabolic changes in hypothalamic structures of the ewe brain between LD and SD. More specifically, the transition between LD and SD revealed negative BOLD responses to hypercapnia at the beginning of SD period followed by significant increases in BOLD, rCBV, Glx and tNAA concentrations towards the end of the SD period. These observations suggest longitudinal mechanisms promoting the proliferation and differentiation of neural stem cells within the hypothalamic niche of breeding ewes. We conclude that multiparametric MRI studies including 1H-MRS could be promising non-invasive translational techniques to investigate the existence of natural adult neurogenesis in-vivo in gyrencephalic brains.

Effects of bolus injection duration on perfusion estimates in dynamic CT and dynamic susceptibility contrast MRI

Estimates of cerebral blood flow (CBF) and tissue mean transit time (MTT) have been shown to differ between dynamic CT perfusion (CTP) and dynamic susceptibility contrast MRI (DSC-MRI). This study investigates whether these discrepancies regarding CBF and MTT between CTP and DSC-MRI can be attributed to the different injection durations of these techniques. Five subjects were scanned using CTP and DSC-MRI. Region-wise estimates of CBF, MTT, and cerebral blood volume (CBV) were derived based on oscillatory index regularized singular value decomposition. A parametric model that reproduced the shape of measured time curves and characteristics of resulting perfusion parameter estimates was developed and used to simulate data with injection durations typical for CTP and DSC-MRI for a clinically relevant set of perfusion scenarios and noise levels. In simulations, estimates of CBF/MTT showed larger negative/positive bias and increasing variability for CTP when compared to DSC-MRI, especially for high CBF levels. While noise also affected estimates, at clinically relevant levels, the injection duration effect was larger. There are several methodological differences between CTP and DSC-MRI. The results of this study suggest that the injection duration is among those that can explain differences in estimates of CBF and MTT between these bolus tracking techniques.

Cumulative effect of simvastatin, L-arginine, and tetrahydrobiopterin on cerebral blood flow and cognitive function in Alzheimer's disease

BACKGROUND AND OBJECTIVES

Vascular disease is a known risk factor for Alzheimer's disease (AD). Endothelial dysfunction has been linked to reduced cerebral blood flow. Endothelial nitric oxide synthase pathway (eNOS) upregulation is known to support endothelial health. This single-center, proof-of-concept study tested whether the use of three medications known to augment the eNOS pathway activity improves cognition and cerebral blood flow (CBF).

METHODS

Subjects with mild AD or mild cognitive impairment (MCI) were sequentially treated with the HMG-CoA reductase synthesis inhibitor simvastatin (weeks 0-16), L-arginine (weeks 4-16), and tetrahydrobiopterin (weeks 8-16). The primary outcome of interest was the change in CBF as measured by MRI from baseline to week 16. Secondary outcomes included standard assessments of cognition.

RESULTS

A total of 11 subjects were deemed eligible and enrolled. One subject withdrew from the study after enrollment, leaving 10 subjects for data analysis. There was a significant increase in CBF from baseline to week 8 by ~13% in the limbic and ~15% in the cerebral cortex. Secondary outcomes indicated a modest but significant increase in the MMSE from baseline (24.2±3.2) to week 16 (26.0±2.7). Exploratory analysis indicated that subjects with cognitive improvement (reduction of the ADAS-cog 13) had a significant increase in their respective limbic and cortical CBF.

CONCLUSIONS

Treatment of mild AD/MCI subjects with medications shown to augment the eNOS pathway was well tolerated and associated with modestly increased cerebral blood flow and cognitive improvement.

TRIAL REGISTRATION

This study is registered in https://www.

CLINICALTRIALS

gov ; registration identifier: NCT01439555; date of registration submitted to registry: 09/23/2011; date of first subject enrollment: 11/2011.

Feasibility of intravoxel incoherent motion in the assessment of tumor microvasculature and blood-brain barrier integrity: a case-based evaluation of gliomas

OBJECTIVE

To evaluate the feasibility of intravoxel incoherent motion (IVIM) in assessing blood-brain barrier (BBB) integrity and microvasculature in tumoral tissue of glioma patients.

METHODS

Images from 8 high-grade and 4 low-grade glioma patients were acquired on a 3 T MRI scanner. Acquisition protocol included pre- and post-contrast T1- and T2-weighted imaging, FLAIR, dynamic susceptibility contrast (DSC), and susceptibility-weighted imaging (SWI). In addition, IVIM was acquired with 15 b-values and fitted under the non-negative least square (NNLS) model to output the diffusion (D) and pseudo-diffusion (D*) coefficients, perfusion fraction (f), and f times D* (fD*) maps.

RESULTS

IVIM perfusion-related maps were sensitive to (1) blood flow and perfusion alterations within the microvasculature of brain tumors, in agreement with intra-tumoral susceptibility signal (ITSS); (2) enhancing areas of BBB breakdown in agreement with DSC maps as well as areas of BBB abnormality that was not detected on DSC maps; (3) enhancing perfusion changes within edemas; (4) detecting early foci of increased perfusion within low-grade gliomas.

CONCLUSION

The results suggest IVIM may be a promising approach to delineate tumor extension and progression in size, and to predict histological grade, which are clinically relevant information that characterize tumors and guide therapeutic decisions in patients with glioma.

Blood oxygen level dependent fMRI and perfusion MRI in the sheep brain

The ovine model could be an effective translational model but remains underexplored. Here, Blood Oxygen Level dependent functional MRI during visual stimulation and resting-state perfusion MRI were explored. We aimed at investigating the impact of isoflurane anesthesia during visual stimulation and evaluate resting cerebral blood flow and cerebral blood volume parameters in the lamb and adult sheep brain. BOLD fMRI and perfusion MRI after a bolus of DOTAREM were conducted in 4 lambs and 6 adult ewes at 3 T. A visual stimulation paradigm was delivered during fMRI at increasing isoflurane doses (1-3%). Robust but weak BOLD responses (0.21 ± 0.08%) were found in the lateral geniculate nucleus (LGN) up to 3% isoflurane anaesthesia. No significant differences were found beween BOLD responses in the range 1 to 3% ISO (p > 0.05). However, LGN cluster size decreased and functional localization became less reliable at high ISO doses (2.5-3% ISO). BOLD responses were weaker in adult sheep than in lambs (4.6 ± 1.5 versus 13.6 ± 8.5; p = 0.08). Relative cerebral blood volumes (rCBV) and relative cerebral blood flows (rCBF) were significantly higher (p < 0.0001) in lambs than in adult sheep for both gray and white matter. The impact of volatile anesthesia was explored for the first time on BOLD responses demonstrating increased reliability of functional localization of brain activity at low doses. Perfusion MRI was conducted for the first time in both lambs and adult ewes. Assessment of baseline cerebrovascular values are of interest for future studies of brain diseases allowing an improved interpretation of BOLD responses.

A hybrid deconvolution approach for estimation of in vivo non-displaceable binding for brain PET targets without a reference region

BACKGROUND AND AIM

Estimation of a PET tracer's non-displaceable distribution volume (VND) is required for quantification of specific binding to its target of interest. VND is generally assumed to be comparable brain-wide and is determined either from a reference region devoid of the target, often not available for many tracers and targets, or by imaging each subject before and after blocking the target with another molecule that has high affinity for the target, which is cumbersome and involves additional radiation exposure. Here we propose, and validate for the tracers [11C]DASB and [11C]CUMI-101, a new data-driven hybrid deconvolution approach (HYDECA) that determines VND at the individual level without requiring either a reference region or a blocking study.

METHODS

HYDECA requires the tracer metabolite-corrected concentration curve in blood plasma and uses a singular value decomposition to estimate the impulse response function across several brain regions from measured time activity curves. HYDECA decomposes each region's impulse response function into the sum of a parametric non-displaceable component, which is a function of VND, assumed common across regions, and a nonparametric specific component. These two components differentially contribute to each impulse response function. Different regions show different contributions of the two components, and HYDECA examines data across regions to find a suitable common VND. HYDECA implementation requires determination of two tuning parameters, and we propose two strategies for objectively selecting these parameters for a given tracer: using data from blocking studies, and realistic simulations of the tracer. Using available test-retest data, we compare HYDECA estimates of VND and binding potentials to those obtained based on VND estimated using a purported reference region.

RESULTS

For [11C]DASB and [11C]CUMI-101, we find that regardless of the strategy used to optimize the tuning parameters, HYDECA provides considerably less biased estimates of VND than those obtained, as is commonly done, using a non-ideal reference region. HYDECA test-retest reproducibility is comparable to that obtained using a VND determined from a non-ideal reference region, when considering the binding potentials BPP and BPND.

CONCLUSIONS

HYDECA can provide subject-specific estimates of VND without requiring a blocking study for tracers and targets for which a valid reference region does not exist.

Stable spline deconvolution for dynamic susceptibility contrast MRI

PURPOSE

To present the stable spline (SS) deconvolution method for the quantification of the cerebral blood flow (CBF) from dynamic susceptibility contrast MRI.

METHODS

The SS method was compared with both the block-circulant singular value decomposition (oSVD) and nonlinear stochastic regularization (NSR) methods. oSVD is one of the most popular deconvolution methods in dynamic susceptibility contrast MRI (DSC-MRI). NSR is an alternative approach that we proposed previously. The three methods were compared using simulated data and two clinical data sets.

RESULTS

The SS method correctly reconstructed the dispersed residue function and its peak in presence of dispersion, regardless of the delay. In absence of dispersion, SS performs similarly to oSVD and does not correctly reconstruct the residue function and its peak. SS and NSR better differentiate healthy and pathologic CBF values compared with oSVD in all simulated conditions. Using acquired data, SS and NSR provide more clinically plausible and physiological estimates of the residue function and CBF maps compared with oSVD.

CONCLUSION

The SS method overcomes some of the limitations of oSVD, such as unphysiological estimates of the residue function and NSR, the latter of which is too computationally expensive to be applied to large data sets. Thus, the SS method is a valuable alternative for CBF quantification using DSC-MRI data. Magn Reson Med 78:1801-1811, 2017. © 2017 International Society for Magnetic Resonance in Medicine.

Model-free quantification of dynamic PET data using nonparametric deconvolution

Dynamic positron emission tomography (PET) data are usually quantified using compartment models (CMs) or derived graphical approaches. Often, however, CMs either do not properly describe the tracer kinetics, or are not identifiable, leading to nonphysiologic estimates of the tracer binding. The PET data are modeled as the convolution of the metabolite-corrected input function and the tracer impulse response function (IRF) in the tissue. Using nonparametric deconvolution methods, it is possible to obtain model-free estimates of the IRF, from which functionals related to tracer volume of distribution and binding may be computed, but this approach has rarely been applied in PET. Here, we apply nonparametric deconvolution using singular value decomposition to simulated and test-retest clinical PET data with four reversible tracers well characterized by CMs ([(11)C]CUMI-101, [(11)C]DASB, [(11)C]PE2I, and [(11)C]WAY-100635), and systematically compare reproducibility, reliability, and identifiability of various IRF-derived functionals with that of traditional CMs outcomes. Results show that nonparametric deconvolution, completely free of any model assumptions, allows for estimates of tracer volume of distribution and binding that are very close to the estimates obtained with CMs and, in some cases, show better test-retest performance than CMs outcomes.

Automatic determination of the arterial input function in dynamic susceptibility contrast MRI: comparison of different reproducible clustering algorithms

Introduction

Arterial input function (AIF) plays an important role in the quantification of cerebral hemodynamics. The purpose of this study was to select the best reproducible clustering method for AIF detection by comparing three algorithms reported previously in terms of detection accuracy and computational complexity.

Methods

First, three reproducible clustering methods, normalized cut (Ncut), hierarchy (HIER), and fast affine propagation (FastAP), were applied independently to simulated data which contained the true AIF. Next, a clinical verification was performed where 42 subjects participated in dynamic susceptibility contrast MRI (DSC-MRI) scanning. The manual AIF and AIFs based on the different algorithms were obtained. The performance of each algorithm was evaluated based on shape parameters of the estimated AIFs and the true or manual AIF. Moreover, the execution time of each algorithm was recorded to determine the algorithm that operated more rapidly in clinical practice.

Results

In terms of the detection accuracy, Ncut and HIER method produced similar AIF detection results, which were closer to the expected AIF and more accurate than those obtained using FastAP method; in terms of the computational efficiency, the Ncut method required the shortest execution time.

Conclusion

Ncut clustering appears promising because it facilitates the automatic and robust determination of AIF with high accuracy and efficiency.

Estimation of arterial arrival time and cerebral blood flow from QUASAR arterial spin labeling using stable spline

PURPOSE

QUASAR arterial spin labeling (ASL) permits the application of deconvolution approaches for the absolute quantification of cerebral perfusion. Currently, oscillation index regularized singular value decomposition (oSVD) combined with edge-detection (ED) is the most commonly used method. Its major drawbacks are nonphysiological oscillations in the impulse response function and underestimation of perfusion. The aim of this work is to introduce a novel method to overcome these limitations.

METHODS

A system identification method, stable spline (SS), was extended to address ASL peculiarities such as the delay in arrival of the arterial blood in the tissue. The proposed framework was compared with oSVD + ED in both simulated and real data. SS was used to investigate the validity of using a voxel-wise tissue T1 value instead of using a single global value (of blood T1 ).

RESULTS

SS outperformed oSVD + ED in 79.9% of simulations. When applied to real data, SS exhibited a physiologically realistic range for perfusion and a higher mean value with respect to oSVD + ED (55.5 ± 9.5 SS, 34.9 ± 5.2 oSVD + ED mL/100 g/min).

CONCLUSION

SS can represent an alternative to oSVD + ED for the quantification of QUASAR ASL data. Analysis of the retrieved impulse response function revealed that using a voxel wise tissue T1 might be suboptimal.

Perfusion magnetic resonance imaging: a comprehensive update on principles and techniques

Perfusion is a fundamental biological function that refers to the delivery of oxygen and nutrients to tissue by means of blood flow. Perfusion MRI is sensitive to microvasculature and has been applied in a wide variety of clinical applications, including the classification of tumors, identification of stroke regions, and characterization of other diseases. Perfusion MRI techniques are classified with or without using an exogenous contrast agent. Bolus methods, with injections of a contrast agent, provide better sensitivity with higher spatial resolution, and are therefore more widely used in clinical applications. However, arterial spin-labeling methods provide a unique opportunity to measure cerebral blood flow without requiring an exogenous contrast agent and have better accuracy for quantification. Importantly, MRI-based perfusion measurements are minimally invasive overall, and do not use any radiation and radioisotopes. In this review, we describe the principles and techniques of perfusion MRI. This review summarizes comprehensive updated knowledge on the physical principles and techniques of perfusion MRI.

Modeling and correction of bolus dispersion effects in dynamic susceptibility contrast MRI

PURPOSE

Bolus dispersion in DSC-MRI can lead to errors in cerebral blood flow (CBF) estimation by up to 70% when using singular value decomposition analysis. However, it might be possible to correct for dispersion using two alternative methods: the vascular model (VM) and control point interpolation (CPI). Additionally, these approaches potentially provide a means to quantify the microvascular residue function.

METHODS

VM and CPI were extended to correct for dispersion by means of a vascular transport function. Simulations were performed at multiple dispersion levels and an in vivo analysis was performed on a healthy subject and two patients with carotid atherosclerotic disease.

RESULTS

Simulations showed that methods that could not address dispersion tended to underestimate CBF (ratio in CBF estimation, CBFratio = 0.57-0.77) in the presence of dispersion; whereas modified CPI showed the best performance at low-to-medium dispersion; CBFratio = 0.99 and 0.81, respectively. The in vivo data showed trends in CBF estimation and residue function that were consistent with the predictions from simulations.

CONCLUSION

In patients with atherosclerotic disease the estimated residue function showed considerable differences in the ipsilateral hemisphere. These differences could partly be attributed to dispersive effects arising from the stenosis when dispersion corrected CPI was used. It is thus beneficial to correct for dispersion in perfusion analysis using this method.

The 39 steps: evading error and deciphering the secrets for accurate dynamic susceptibility contrast MRI

Dynamic susceptibility contrast (DSC) MRI is the most commonly used MRI method to assess cerebral perfusion and other related haemodynamic parameters. Although the technique is well established and used routinely in clinical centres, there are still many problems that impede accurate perfusion quantification. In this review article, we present 39 steps which guide the reader through the theoretical principles, practical decisions, potential problems, current limitations and latest advances in DSC-MRI. The 39 steps span the collection, analysis and interpretation of DSC-MRI data, expounding issues and possibilities relating to the contrast agent, the acquisition of DSC-MRI data, data pre-processing, the contrast concentration-time course, the arterial input function, deconvolution, common perfusion parameters, post-processing possibilities, patient studies, absolute versus relative quantification and automated analysis methods.

Bayesian hemodynamic parameter estimation by bolus tracking perfusion weighted imaging

A delay-insensitive probabilistic method for estimating hemodynamic parameters, delays, theoretical residue functions, and concentration time curves by computed tomography (CT) and magnetic resonance (MR) perfusion weighted imaging is presented. Only a mild stationarity hypothesis is made beyond the standard perfusion model. New microvascular parameters with simple hemodynamic interpretation are naturally introduced. Simulations on standard digital phantoms show that the method outperforms the oscillating singular value decomposition (oSVD) method in terms of goodness-of-fit, linearity, statistical and systematic errors on all parameters, especially at low signal-to-noise ratios (SNRs). Delay is always estimated sharply with user-supplied resolution and is purely arterial, by contrast to oSVD time-to-maximum TMAX that is very noisy and biased by mean transit time (MTT), blood volume, and SNR. Residue functions and signals estimates do not suffer overfitting anymore. One CT acute stroke case confirms simulation results and highlights the ability of the method to reliably estimate MTT when SNR is low. Delays look promising for delineating the arterial occlusion territory and collateral circulation.

Interventional 4-D C-arm CT perfusion imaging using interleaved scanning and partial reconstruction interpolation

Tissue perfusion measurement during catheter-guided stroke treatment in the interventional suite is currently not possible. In this work, we present a novel approach that uses a C-arm angiography system capable of computed tomography (CT)-like imaging (C-arm CT) for this purpose. With C-arm CT one reconstructed volume can be obtained every 4-6 s which makes it challenging to measure the flow of an injected contrast bolus. We have developed an interleaved scanning (IS) protocol that uses several scan sequences to increase temporal sampling. Using a dedicated 4-D reconstruction approach based on partial reconstruction interpolation (PRI) we can optimally process our data. We evaluated our combined approach (IS-PRI) with simulations and a study in five healthy pigs. In our simulations, the cerebral blood flow values (unit: ml/100 g/min) were 60 (healthy tissue) and 20 (pathological tissue). For one scan sequence the values were estimated with standard deviations of 14.3 and 2.9, respectively. For two interleaved sequences the standard deviations decreased to 3.6 and 1.5, respectively. We used perfusion CT to validate the in vivo results. With two interleaved sequences we achieved promising correlations ranging from r=0.63 to r=0.94. The results suggest that C-arm CT tissue perfusion imaging is feasible with two interleaved scan sequences.

Automatic selection of arterial input function on dynamic contrast-enhanced MR images

Dynamic susceptibility contrast-magnetic resonance imaging (DSC-MRI) data analysis requires the knowledge of the arterial input function (AIF) to quantify the cerebral blood flow (CBF), volume (CBV) and the mean transit time (MTT). AIF can be obtained either manually or using automatic algorithms. We present a method to derive the AIF on the middle cerebral artery (MCA). The algorithm draws a region of interest (ROI) where the MCA is located. Then, it uses a recursive cluster analysis on the ROI to select the arterial voxels. The algorithm had been compared on simulated data to literature state of art automatic algorithms and on clinical data to the manual procedure. On in silico data, our method allows to reconstruct the true AIF and it is less affected by partial volume effect bias than the other methods. In clinical data, automatic AIF provides CBF and MTT maps with a greater contrast level compared to manual AIF ones. Therefore, AIF obtained with the proposed method improves the estimate reliability and provides a quantitatively reliable physiological picture.

Deconvolution-Based CT and MR Brain Perfusion Measurement: Theoretical Model Revisited and Practical Implementation Details

Deconvolution-based analysis of CT and MR brain perfusion data is widely used in clinical practice and it is still a topic of ongoing research activities. In this paper, we present a comprehensive derivation and explanation of the underlying physiological model for intravascular tracer systems. We also discuss practical details that are needed to properly implement algorithms for perfusion analysis. Our description of the practical computer implementation is focused on the most frequently employed algebraic deconvolution methods based on the singular value decomposition. In particular, we further discuss the need for regularization in order to obtain physiologically reasonable results. We include an overview of relevant preprocessing steps and provide numerous references to the literature. We cover both CT and MR brain perfusion imaging in this paper because they share many common aspects. The combination of both the theoretical as well as the practical aspects of perfusion analysis explicitly emphasizes the simplifications to the underlying physiological model that are necessary in order to apply it to measured data acquired with current CT and MR scanners.

Assessment of clinical data of nonlinear stochastic deconvolution versus block-circulant singular value decomposition for quantitative dynamic susceptibility contrast magnetic resonance imaging

Dynamic susceptibility contrast magnetic resonance imaging (DSC-MRI) allows the noninvasive assessment of brain hemodynamics alterations by quantifying, via deconvolution, the cerebral blood flow (CBF) and mean transit time (MTT). Singular value decomposition (SVD) and block-circulant SVD (cSVD) are the most widely adopted deconvolution method, although they bear some limitations, including unphysiological oscillations in the residue function and bias in the presence of delay and dispersion between the tissue and the arterial input function. A nonlinear stochastic regularization (NSR) has been proposed, which performs better than SVD and cSVD on simulated data both in the presence and absence of dispersion. Moreover, NSR allows to quantify the dispersion level. Here, cSVD and NSR are compared for the first time on a group of nine patients with severe atherosclerotic unilateral stenosis of internal carotid artery before and after carotid stenting to investigate the effect of arterial dispersion. According to region of interest-based analysis, NSR characterizes the pathologic tissue more accurately than cSVD, thus improving the quality of the information provided to physicians for diagnosis. In fact, in 7 (78%) of the 9 subjects, CBF and MTT maps provided by NSR allow to correctly identify the pathologic hemisphere to the physician. Moreover, by emphasizing the difference between pathologic and healthy tissues, NSR may be successfully used to monitor the subject's recovery after the treatment and/or surgery. NSR also generates dispersion level and non-dispersed CBF and MTT maps. The dispersion level provides information on CBF and MTT estimates reliability and may also be used as a clinical indicator of pathological tissue state complementary to CBF and MTT, thus increasing the clinical information provided by DSC-MRI analysis.