Improvements In computed tomography perfusion output using complex singular value decomposition and the maximum slope algorithm

OBJECTIVE: Determine if complex singular value decomposition (cSVD) used as preprocessing in the maximum slope algorithm reduces image noise of resultant physiologic parametric images. Noise will be decreased in the parametric maps of cerebral blood flow (CBF), cerebral blood volume (CBV) as compared to the same algorithm and data set with no cSVD applied. MATERIALS AND METHODS: A set of 10 patients (n=15) underwent a total combined 15 CT perfusion studies upon presenting with stroke symptoms. It was determined these patients suffered from occlusions resulting in a prolonged arrival time of blood to the brain. DICOM data files of these patients scans were selected based on this increased arrival delay. We compared the output of estimation calculations for cerebral blood flow (CBF), and cerebral blood volume (CBV), using preprocessing cSVD against the same scan data with no preprocessing cSVD. Image noise was assessed through the calculation of the standard deviation within specific regions of interest copied to specific areas of grey and white matter as well as CSF space. A decrease in the standard deviation values will indicate improvement in the noise level of the resultant images.. Results for the mean value within the regions of interest are expected to be similar between the groups calculated using cSVD and those calculated under the standard method. This will indicate the presence of minimal bias. RESULTS: Between groups of the standard processing method and the cSVD method standard deviation (SD) reductions were seen in both CBF and CBV values across all three ROIs. In grey matter measures of CBV, SD was reduced an average of 0.0034 mL/100g while measures of CBF saw SD reduced by an average of 0.073 mL/100g/min. In samples of white matter, standard deviations of CBV values were reduced on average by 0.0041mL/100g while CBF SD's were reduced by 0.073 mL/100g/min. CSF ROIs in CBV calculations saw SD reductions averaging 0.0047 mL/100g and reductions of 0.074 mL/100g/min in measures of CBF. Bias within CBV calculations was at most minimal as determined by no significant changes in mean calculated values. Calculations of CBF saw large downward bias in the mean values. CONCLUSIONS: The application of the cSVD method to preprocessing of CT perfusion imaging studies produces an effective method of noise reduction. In calculations of CBV, cSVD noise reduction results in overall improvement. In calculations of CBF, cSVD, while effective in noise reduction, caused mean values to be statistically lower than the standard method. It should be noted that there is currently no evaluation of which values can be considered more accurate physiologically. Simulations of the effect of noise on CBF showed a positive correlation suggesting that the CBF algorithm itself is sensitive to the level of noise

Parallel perfusion imaging processing using GPGPU

AbstractBackground and purposeThe objective of brain perfusion quantification is to generate parametric maps of relevant hemodynamic quantities such as cerebral blood flow (CBF), cerebral blood volume (CBV) and mean transit time (MTT) that can be used in diagnosis of acute stroke. These calculations involve deconvolution operations that can be very computationally expensive when using local Arterial Input Functions (AIF). As time is vitally important in the case of acute stroke, reducing the analysis time will reduce the number of brain cells damaged and increase the potential for recovery.MethodsGPUs originated as graphics generation dedicated co-processors, but modern GPUs have evolved to become a more general processor capable of executing scientific computations. It provides a highly parallel computing environment due to its large number of computing cores and constitutes an affordable high performance computing method. In this paper, we will present the implementation of a deconvolution algorithm for brain perfusion quantification on GPGPU (General Purpose Graphics Processor Units) using the CUDA programming model. We present the serial and parallel implementations of such algorithms and the evaluation of the performance gains using GPUs.ResultsOur method has gained a 5.56 and 3.75 speedup for CT and MR images respectively.ConclusionsIt seems that using GPGPU is a desirable approach in perfusion imaging analysis, which does not harm the quality of cerebral hemodynamic maps but delivers results faster than the traditional computation

Automatic quantification of the microvascular density on whole slide images, applied to paediatric brain tumours

Angiogenesis is a key phenomenon for tumour progression, diagnosis and treatment in brain tumours and more generally in oncology. Presently, its precise, direct quantitative assessment can only be done on whole tissue sections immunostained to reveal vascular endothelial cells. But this is a tremendous task for the pathologist and a challenge for the computer since digitised whole tissue sections, whole slide images (WSI), contain typically around ten gigapixels. We define and implement an algorithm that determines automatically, on a WSI at objective magnification $40\times$, the regions of tissue, the regions without blur and the regions of large puddles of red blood cells, and constructs the mask of blur-free, significant tissue on the WSI. Then it calibrates automatically the optical density ratios of the immunostaining of the vessel walls and of the counterstaining, performs a colour deconvolution inside the regions of blur-free tissue, and finds the vessel walls inside these regions by selecting, on the image resulting from the colour deconvolution, zones which satisfy a double-threshold criterion. A mask of vessel wall regions on the WSI is produced. The density of microvessels is finally computed as the fraction of the area of significant tissue which is occupied by vessel walls. We apply this algorithm to a set of 186 WSI of paediatric brain tumours from World Health Organisation grades I to IV. The segmentations are of very good quality although the set of slides is very heterogeneous. The computation time is of the order of a fraction of an hour for each WSI on a modest computer. The computed microvascular density is found to be robust and strongly correlates with the tumour grade. This method requires no training and can easily be applied to other tumour types and other stainings

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

Augmented breast tumor classification by perfusion analysis

Magnetic resonance and computed tomography imaging aid in the diagnosis and analysis of pathologic conditions. Blood flow, or perfusion, through a region of tissue can be computed from a time series of contrast-enhanced images. Perfusion is an important set of physiological parameters that reflect angiogenesis. In cancer, heightened angiogenesis is a key process in the growth and spread of tumorous masses. An automatic classification technique using recovered perfusion may prove to be a highly accurate diagnostic tool. Such a classification system would supplement existing histopathological tests, and help physicians to choose the most optimal treatment protocol. Perfusion is obtained through deconvolution of signal intensity series and a pharmacokinetic model. However, many computational problems complicate the accurate-consistent recovery of perfusion. The high time-resolution acquisition of images decreases signal-to-noise, producing poor deconvolution solutions. The delivery of contrast agent as a function of time must also be determined or sampled before deconvolution can proceed. Some regions of the body, such as the brain, provide a nearby artery to serve as this arterial input function. Poor estimates can lead to an over or under estimation of perfusion. Breast tissue is an example of one tissue region where a clearly defined artery is not present. This proposes a new method of using recovered perfusion and spatial information in an automated classifier. This classifier grades suspected lesions as benign or malignant. This method can be integrated into a computer-aided diagnostic system to enhance the value of medical imagery