Abstract-

Abstract. In order to improve the robustness of rigid registration algorithms in various medical imaging problems, we propose in this article a general framework built on block matching strategies. This framework combines two stages in a multi-scale hierarchy. The first stage consists in finding for each block (or subregion) of the first image, the most similar subregion in the other image, using a similarity criterion which depends on the nature of the images. The second stage consists in finding the global rigid transformation which best explains most of these local correspondances. This is done with a robust procedure which allows up to 50% of false matches. We show that this approach, besides its simplicity, provides a robust and efficient way to rigidly register images in various situations. This includes for instance the alignment of 2D histological sections for the 3D reconstructions of trimmed organs and tissues, the automatic computation of the mid-sagittal plane in multimodal 3D images of the brain, and the multimodal registration of 3D CT and MR images of the brain. A quantitative evaluation of the results is provided for this last example, as well as a comparison with the classical approaches involving the minimization of a global measure of similarity based on Mutual Information or the Correlation Ratio. This shows a significant improvement of the robustness, for a comparable final accuracy. Although slightly more expensive in terms of computational requirements, the proposed approach can easily be implemented on a parallel architecture, which opens potentialities for real time applications using a large number of processors.

Diffusion-Weighted MRI (DW-MRI) is subject to random noise yielding measures that are different from their real values, and thus biasing the subsequently estimated tensors. The Non-Local Means (NLMeans) filter has recently been proposed to denoise MRI with high signal-to-noise ratio (SNR). This filter has been shown to allow the best restoration of image intensities for the estimation of diffusion tensors (DT) compared to state-of-the-art methods. However, for DW-MR images with high b-values (and thus low SNR), the noise, which is strictly Rician-distributed, can no longer be approximated as additive white Gaussian, as implicitly assumed in the classical formulation of the NLMeans. High b-values are typically used in high angular resolution diffusion imaging (HARDI) or q-space imaging (QSI), for which an optimal restoration is critical. In this paper, we propose to adapt the NLMeans filter to Rician noise corrupted data. Validation is performed on synthetic data and on real data for both conventional MR images and DT images. Our adaptation outperforms the original NLMeans filter in terms of peak-signal-to-noise ratio (PSNR) for DW-MRI.

Background and Purpose Diffusion tensor imaging tractography reconstruction of white matter pathways can help guide brain tumor resection. However, DTI tracts are complex mathematical objects and the validity of tractography-derived information in clinical settings has yet to be fully established. To address this issue, we initiated the DTI Challenge, an international working group of clinicians and scientists whose goal was to provide standardized evaluation of tractography methods for neurosurgery. The purpose of this empirical study was to evaluate different tractography techniques in the first DTI Challenge workshop. Methods Eight international teams from leading institutions reconstructed the pyramidal tract in four neurosurgical cases presenting with a glioma near the motor cortex. Tractography methods included deterministic, probabilistic, filtered, and global approaches. Standardized evaluation of the tracts consisted in the qualitative review of the pyramidal pathways by a panel of neurosurgeons and DTI experts and the quantitative evaluation of the degree of agreement among methods. Results The evaluation of tractography reconstructions showed a great inter-algorithm variability. Although most methods found projections of the pyramidal tract from the medial portion of the motor strip, only a few algorithms could trace the lateral projections from the hand, face, and tongue area. In addition, the structure of disagreement among methods was similar across hemispheres despite the anatomical distortions caused by pathological tissues. Conclusions The DTI Challenge provides a benchmark for the standardized evaluation of tractography methods on neurosurgical data. This study suggests that there are still limitations to the clinical use of tractography for neurosurgical decision-making.

Abstract. In this paper we study the impact of denoising the raw high angular resolution diffusion imaging (HARDI) data with the Non-Local Means filter adapted to Rician noise (NLMr). We first show that NLMr filtering improves robustness of apparent diffusion coefficient (ADC) and orientation distribution function (ODF) reconstructions from synthetic HARDI datasets. Our results suggest that the NLMr filtering improve the quality of anisotropy maps computed from ADC and ODF and improve the coherence of q-ball ODFs with the underlying anatomy while not degrading angular resolution. These results are shown on a biological phantom with known ground truth and on a real human brain dataset. Most importantly, we show that multiple measurements of diffusion-weighted (DW) images and averaging these images along each direction can be avoided because NLMr filtering of the individual DW images produces better quality generalized fractional anisotropy maps and more accurate ODF fields than when computed from the averaged DW datasets.

Abstract. Diffusion tensor imaging (DT-MRI) is very sensitive to corrupting noise due to the non linear relationship between the diffusionweighted image intensities (DW-MRI) and the resulting diffusion tensor. Denoising is a crucial step to increase the quality of the estimated tensor field. This enhanced quality allows for a better quantification and a better image interpretation. The methods proposed in this paper are based on the Non-Local (NL) means algorithm. This approach uses the natural redundancy of information in images to remove the noise. We introduce three variations of the NL-means algorithms adapted to DW-MRI and to DT-MRI. Experiments were carried out on a set of 12 diffusion-weighted images (DW-MRI) of the same subject. The results show that the intensity based NL-means approaches give better results in the context of DT-MRI than other classical denoising methods, such as Gaussian Smoothing, Anisotropic Diffusion and Total Variation.

We present a new method to automatically compute, reorient, and recenter the mid-sagittal plane in anatomical and functional three-dimensional (3-D) brain images. This iterative approach is composed of two steps. At first, given an initial guess of the mid-sagittal plane (generally, the central plane of the image grid), the computation of local similarity measures between the two sides of the head allows to identify homologous anatomical structures or functional areas, by way of a block matching procedure. The output is a set of point-to-point correspondences: the centers of homologous blocks. Subsequently, we define the mid-sagittal plane as the one best superposing the points on one side and their counterparts on the other side by reflective symmetry. Practically, the computation of the parameters characterizing the plane is performed by a least trimmed squares estimation. Then, the estimated plane is aligned with the center of the image grid, and the whole process is iterated until convergence. The robust estimation technique we use allows normal or abnormal asymmetrical structures or areas to be treated as outliers, and the plane to be mainly computed from the underlying gross symmetry of the brain. The algorithm is fast and accurate, even for strongly tilted heads, and even in presence of high acquisition noise and bias field, as shown on a large set of synthetic data. The algorithm has also been visually evaluated on a large set of real magnetic resonance (MR) images. We present a few results on isotropic as well as anisotropic anatomical (MR and computed tomography) and functional (single photon emission computed tomography and positron emission tomography) real images, for normal and pathological subjects.

No abstract