Benchmarking framework for myocardial tracking and deformation algorithms: an open access database.

In this paper we present a benchmarking framework for the validation of cardiac motion analysis algorithms. The reported methods are the response to an open challenge that was issued to the medical imaging community through a MICCAI workshop. The database included magnetic resonance (MR) and 3D ultrasound (3DUS) datasets from a dynamic phantom and 15 healthy volunteers. Participants processed 3D tagged MR datasets (3DTAG), cine steady state free precession MR datasets (SSFP) and 3DUS datasets, amounting to 1158 image volumes. Ground-truth for motion tracking was based on 12 landmarks (4 walls at 3 ventricular levels). They were manually tracked by two observers in the 3DTAG data over the whole cardiac cycle, using an in-house application with 4D visualization capabilities. The median of the inter-observer variability was computed for the phantom dataset (0.77 mm) and for the volunteer datasets (0.84 mm). The ground-truth was registered to 3DUS coordinates using a point based similarity transform. Four institutions responded to the challenge by providing motion estimates for the data: Fraunhofer MEVIS (MEVIS), Bremen, Germany; Imperial College London - University College London (IUCL), UK; Universitat Pompeu Fabra (UPF), Barcelona, Spain; Inria-Asclepios project (INRIA), France. Details on the implementation and evaluation of the four methodologies are presented in this manuscript. The manually tracked landmarks were used to evaluate tracking accuracy of all methodologies. For 3DTAG, median values were computed over all time frames for the phantom dataset (MEVIS=1.20mm, IUCL=0.73 mm, UPF=1.10mm, INRIA=1.09 mm) and for the volunteer datasets (MEVIS=1.33 mm, IUCL=1.52 mm, UPF=1.09 mm, INRIA=1.32 mm). For 3DUS, median values were computed at end diastole and end systole for the phantom dataset (MEVIS=4.40 mm, UPF=3.48 mm, INRIA=4.78 mm) and for the volunteer datasets (MEVIS=3.51 mm, UPF=3.71 mm, INRIA=4.07 mm). For SSFP, median values were computed at end diastole and end systole for the phantom dataset(UPF=6.18 mm, INRIA=3.93 mm) and for the volunteer datasets (UPF=3.09 mm, INRIA=4.78 mm). Finally, strain curves were generated and qualitatively compared. Good agreement was found between the different modalities and methodologies, except for radial strain that showed a high variability in cases of lower image quality.

A novel anisotropic fast marching method and its application to blood flow computation in phase-contrast MRI.

BACKGROUND: Phase-contrast MRI (PC MRI) can be used to assess blood flow dynamics noninvasively inside the human body. The acquired images can be reconstructed into flow vector fields. Traditionally, streamlines can be computed based on the vector fields to visualize flow patterns and particle trajectories. OBJECTIVES: The traditional methods may give a false impression of precision, as they do not consider the measurement uncertainty in the PC MRI images. In our prior work, we incorporated the uncertainty of the measurement into the computation of particle trajectories. METHODS: As a major part of the contribution, a novel numerical scheme for solving the anisotropic Fast Marching problem is presented. A computing time comparison to state-of-the-art methods is conducted on artificial tensor fields. A visual comparison of healthy to pathological blood flow patterns is given. RESULTS: The comparison shows that the novel anisotropic Fast Marching solver outperforms previous schemes in terms of computing time. The visual comparison of flow patterns directly visualizes large deviations of pathological flow from healthy flow. CONCLUSIONS: The novel anisotropic Fast Marching solver efficiently resolves even strongly anisotropic path costs. The visualization method enables the user to assess the uncertainty of particle trajectories derived from PC MRI images.

Detection of mild emphysema by computed tomography density measurements.

PURPOSE: To assess the ability of a conventional density mask method to detect mild emphysema by high-resolution computed tomography (HRCT); to analyze factors influencing quantification of mild emphysema; and to validate a new algorithm for detection of mild emphysema. MATERIAL AND METHODS: Fifty-five healthy male smokers and 34 never-smokers, 61-62 years of age, were examined. Emphysema was evaluated visually, by the conventional density mask method, and by a new algorithm compensating for the effects of gravity and artifacts due to motion and the reconstruction algorithm. Effects of the reconstruction algorithm, slice thickness, and various threshold levels on the outcome of the density mask area were evaluated. RESULTS: Forty-nine percent of the smokers had mild emphysema. The density mask area was higher the thinner the slice irrespective of the reconstruction algorithm and threshold level. The sharp algorithm resulted in increased density mask area. The new reconstruction algorithm could discriminate between smokers with and those without mild emphysema, whereas the density mask method could not. The diagnostic ability of the new algorithm was dependent on lung level. At about 90% specificity, sensitivity was 65-100% in the apical levels, but low in the rest of the lung. CONCLUSION: The conventional density mask method is inadequate for detecting mild emphysema, while the new algorithm improves the diagnostic ability but is nevertheless still imperfect.

Detection of neural activity in functional MRI using canonical correlation analysis.

A novel method for detecting neural activity in functional magnetic resonance imaging (fMRI) data is introduced. It is based on canonical correlation analysis (CCA), which is a multivariate extension of the univariate correlation analysis widely used in fMRI. To detect homogeneous regions of activity, the method combines a subspace modeling of the hemodynamic response and the use of spatial relationships. The spatial correlation that undoubtedly exists in fMR images is completely ignored when univariate methods such as as t-tests, F-tests, and ordinary correlation analysis are used. Such methods are for this reason very sensitive to noise, leading to difficulties in detecting activation and significant contributions of false activations. In addition, the proposed CCA method also makes it possible to detect activated brain regions based not only on thresholding a correlation coefficient, but also on physiological parameters such as temporal shape and delay of the hemodynamic response. Excellent performance on real fMRI data is demonstrated. Magn Reson Med 45:323-330, 2001.