Online atlasing, i.e., incrementing an atlas with new images as they are acquired, is key when performing studies on very large, or still being gathered, databases. Regular approaches to atlasing however do not focus on this aspect and impose a complete reconstruction of the atlas when adding images. We propose instead a diffeomorphic online atlasing method that allows gradual updates to an atlas. In this iterative centroid approach, we integrate new subjects in the atlas in an iterative manner, gradually moving the centroid of the images towards its final position. This leads to a computationally cheap approach since it only necessitates one additional registration per new subject added. We validate our approach on several experiments with three main goals: 1- to evaluate atlas image quality of the obtained atlases with sharpness and overlap measures, 2- to assess the deviation in terms of transformations with respect to a conventional atlasing method and 3- to compare its computational time with regular approaches of the literature. We demonstrate that the transformations divergence with respect to a state-of-the-art atlas construction method is small and reaches a plateau, that the two construction methods have the same ability to map subject homologous regions onto a common space and produce images of equivalent quality. The computational time of our approach is also drastically reduced for regular updates. Finally, we also present a direct extension of our method to update spatio-temporal atlases, especially useful for developmental studies.

Multi-compartment models (MCM) are increasingly used to characterize the brain white matter microstructure from diffusion-weighted imaging (DWI). Their use in clinical studies is however limited by the inability to resample an MCM image towards a common reference frame, or to construct atlases from such brain microstructure models. We propose to solve this problem by first identifying that these two tasks amount to the same problem. We propose to tackle it by viewing it as a simplification problem, solved thanks to spectral clustering and the definition of semi-metrics between several usual compartments encountered in the MCM literature. This generic framework is evaluated for two models: the multi-tensor model where individual fibers are modeled as individual tensors and the diffusion direction imaging (DDI) model that differentiates intra- and extra-axonal components of each fiber. Results on simulated data, simulated transformations and real data show the ability of our method to well interpolate MCM images of these types. We finally present as an application an MCM template of normal controls constructed using our approach.

By shortening the acquisition time of MRI, Echo Planar Imaging (EPI) enables the acquisition of a large number of images in a short time, compatible with clinical constraints as required for diffusion or functional MRI. However such images are subject to large, local distortions disrupting their correspondence with the underlying anatomy. The correction of those distortions is an open problem, especially in regions where large deformations occur. We propose a new block-matching registration method to perform EPI distortion correction based on the acquisition of two EPI with opposite phase encoding directions (PED). It relies on new transformations between blocks adapted to the EPI distortion model, and on an adapted optimization scheme to ensure an opposite symmetric transformation. We present qualitative and quantitative results of the block-matching correction using different metrics on a phantom dataset and on in-vivo data. We show the ability of the block-matching to robustly correct EPI distortion even in strongly affected areas.

Quantifying \(T_2\) and \(T_2^*\) relaxation times from MRI becomes a standard tool to assess modifications of biological tissues over time or differences between populations. However, due to the relationship between the relaxation time and the associated MR signals such an analysis is subject to error. In this work, we provide a Bayesian analysis of this relationship. More specifically, we build posterior distributions relating the raw (spin or gradient echo) acquisitions and the relaxation time and its modifications over acquisitions. Such an analysis has three main merits. First, it allows to build hierarchical models including prior information and regularisations over voxels. Second, it provides many estimators of the parameters distribution including the mean and the \(\alpha \) -credible intervals. Finally, as credible intervals are available, testing properly whether the relaxation time (or its modification) lies within a certain range with a given credible level is simple. We show the interest of this approach on synthetic datasets and on two real applications in multiple sclerosis. These keywords were added by machine and not by the authors. This process is experimental and the keywords may be updated as the learning algorithm improves.

This paper presents a new image-based visual servoing approach to control a robotic system equipped with an ultrasound (US) imaging device. The presented method allows an automatic positioning of the probe with respect to an object of interest. Moments-based image features are computed from three orthogonal US images to servo in-plane and out-of-plane motions of the system. An efficient segmentation method, based on graph-cut strategy, is proposed to extract the object contour in each image plane. Simulation results demonstrate that this approach improves upon techniques based on a single 2-D US image in terms of probe positioning. Our method was also validated from robotic experiments performed on a US phantom with the use of a motorized 3-D probe that provides the three US images.

We present a new automatic method for segmentation of multiple sclerosis (MS) lesions in magnetic resonance images. The method performs tissue classification using a model of intensities of the normal appearing brain tissues. In order to estimate the model, a trimmed likelihood estimator is initialized with a hierarchical random approach in order to be robust to MS lesions and other outliers present in real images. The algorithm is first evaluated with simulated images to assess the importance of the robust estimator in presence of outliers. The method is then validated using clinical data in which MS lesions were delineated manually by several experts. Our method obtains an average Dice similarity coefficient (DSC) of 0.65, which is close to the average DSC obtained by raters (0.66).

In image processing, restoration is expected to improve the qualitative inspection of the image and the performance of quantitative image analysis techniques. In this paper, an adaptation of the nonlocal (NL)-means filter is proposed for speckle reduction in ultrasound (US) images. Originally developed for additive white Gaussian noise, we propose to use a Bayesian framework to derive a NL-means filter adapted to a relevant ultrasound noise model. Quantitative results on synthetic data show the performances of the proposed method compared to well-established and state-of-the-art methods. Results on real images demonstrate that the proposed method is able to preserve accurately edges and structural details of the image.

A critical issue in image restoration is the problem of noise removal while keeping the integrity of relevant image information. Denoising is a crucial step to increase image quality and to improve the performance of all the tasks needed for quantitative imaging analysis. The method proposed in this paper is based on a 3-D optimized blockwise version of the nonlocal (NL)-means filter (Buades, et al., 2005). The NL-means filter uses the redundancy of information in the image under study to remove the noise. The performance of the NL-means filter has been already demonstrated for 2-D images, but reducing the computational burden is a critical aspect to extend the method to 3-D images. To overcome this problem, we propose improvements to reduce the computational complexity. These different improvements allow to drastically divide the computational time while preserving the performances of the NL-means filter. A fully automated and optimized version of the NL-means filter is then presented. Our contributions to the NL-means filter are: 1) an automatic tuning of the smoothing parameter; 2) a selection of the most relevant voxels; 3) a blockwise implementation; and 4) a parallelized computation. Quantitative validation was carried out on synthetic datasets generated with BrainWeb (Collins, et al., 1998). The results show that our optimized NL-means filter outperforms the classical implementation of the NL-means filter, as well as two other classical denoising methods [anisotropic diffusion (Perona and Malik, 1990)] and total variation minimization process (Rudin, et al., 1992) in terms of accuracy (measured by the peak signal-to-noise ratio) with low computation time. Finally, qualitative results on real data are presented .

Automated computer classification (ACC) techniques are needed to facilitate physician's diagnosis of complex diseases in individual patients. We provide an example of ACC using computational techniques within the context of cross-sectional analysis of magnetic resonance images (MRI) in neurodegenerative diseases, namely Alzheimer's dementia (AD). In this paper, the accuracy of our ACC methodology is assessed when presented with real life, imperfect data, i.e., cohorts of MRI with varying acquisition parameters and imaging quality. The comparative methodology uses the Jacobian determinants derived from dense deformation fields and scaled grey-level intensity from a selected volume of interest centered on the medial temporal lobe. The ACC performance is assessed in a series of leave-one-out experiments aimed at separating 75 probable AD and 75 age-matched normal controls. The resulting accuracy is 92% using a support vector machine classifier based on least squares optimization. Finally, it is shown in the Appendix that determinants and scaled grey-level intensity are appreciably more robust to varying parameters in validation studies using simulated data, when compared to raw intensities or grey/white matter volumes. The ability of cross-sectional MRI at detecting probable AD with high accuracy could have profound implications in the management of suspected AD candidates.

This paper examines a novel approach for temporal calibration of a three-dimensional (3-D) freehand ultrasound system. A localization system fixed on the probe gives the position and orientation of the probe. For quantitative use, calibration is needed to correctly localize a B-scan in four-dimensional (4-D) (3-D+t) space. Temporal latency estimation is defined in a general robust formulation using no specific probe motion constraints. Experiments were performed on synthetic and real data using a 3-D freehand ultrasound system. The achieved precision is lower than the image acquisition rate (40 ms). A validation study using a calibration phantom has been performed to evaluate the influence of incorrect latency estimation on the 3-D reconstruction procedure. We showed that for latency estimation errors less than 40 ms, the 3-D reconstruction errors are negligible for volume estimation.

In this paper, three different calibration methods for three-dimensional (3-D) freehand ultrasound (US) are evaluated. Calibration is the process of estimating the rigid transformation from US image coordinates to the coordinate system of the tracking sensor mounted onto the probe. Calibration accuracy has an important impact on quantitative studies. Geometrical precision can also be crucial in many interventions and surgery. The proposed evaluation framework relies on a single point phantom and a 3-D US phantom which mimics the US characteristics of human liver. Four quality measures are used: 3-D point localization criterion, distance and volume measurements, and shape based criterion. Results show that during the acquisition procedure, volumetric measurements and shapes of the reconstructed object depend on probe motion used, particularly fan motions for which errors are larger. It is also shown that accurate calibration is essential to obtain reliable quantitative information.

A new registration method for ultrasound volumes relying on on a statistical texture-based similarity measure is investigated. Texture information is given by spatial Gabor filters and represented by statistical kernel-based distributions. The registration similarity measure is then defined as a probabilistic distance, derived from Bhattacharyya coefficient, between two statistical distributions. Given this similarity measure, parametric ultrasound image registration is stated as a robust minimization issue. We also exploit frequency properties of spatial Gabor filters to propose a multiresolution approach to perform this minimization. We provide a preliminary evaluation of the new registration technique on clinical data.

In this paper, we investigate the introduction of cortical constraints for non rigid intersubject brain registration. We extract sulcal patterns with the active ribbon method, presented by Le Goualher et al. (1997). An energy based registration method (Hellier et al., 2001), which will be called photometric registration method in this paper, makes it possible to incorporate the matching of cortical sulci. The local sparse similarity and the photometric similarity are, thus, expressed in a unified framework. We show the benefits of cortical constraints on a database of 18 subjects, with global and local assessment of the registration. This new registration scheme has also been evaluated on functional magnetoencephalography data. We show that the anatomically constrained registration leads to a substantial reduction of the intersubject functional variability.

Although numerous methods to register brains of different individuals have been proposed, no work has been done, as far as we know, to evaluate and objectively compare the performances of different nonrigid (or elastic) registration methods on the same database of subjects. In this paper, we propose an evaluation framework, based on global and local measures of the relevance of the registration. We have chosen to focus more particularly on the matching of cortical areas, since intersubject registration methods are dedicated to anatomical and functional normalization, and also because other groups have shown the relevance of such registration methods for deep brain structures. Experiments were conducted using 6 methods on a database of 18 subjects. The global measures used show that the quality of the registration is directly related to the transformation's degrees of freedom. More surprisingly, local measures based on the matching of cortical sulci did not show significant differences between rigid and non rigid methods.

In this paper we address the problem of registering 3D scattered data by the mean of a statistical shape model. This model is built from a training set on which a principal component analysis (PCA) is applied. A local system of reference is computed for each sample shape of the learning set, which enables to align the training set. PCA then reveals the main modes of deformation of the class of objects of interest. Furthermore, the deformation field obtained between a given shape and a reference shape is extended to a local neighborhood of these shapes by using an interpolation based on the thin-plate splines. It is then used to register objects associated with these shapes in a local and non-linear way. The data involved here are cerebral data, both anatomical (cortical sulci) and functional (MEG dipoles).

A new method for medical image registration is formulated as a minimization problem involving robust estimators. We propose an efficient hierarchical optimization framework which is both multiresolution and multigrid. An anatomical segmentation of the cortex is introduced in the adaptive partitioning of the volume on which the multigrid minimization is based. This allows to limit the estimation to the areas of interest, to accelerate the algorithm, and to refine the estimation in specified areas. At each stage of the hierarchical estimation, we refine current estimate by seeking a piecewise affine model for the incremental deformation field. The performance of this method is numerically evaluated on simulated data and its benefits and robustness are shown on a database of 18 magnetic resonance imaging scans of the head.

In this paper we describe a new method for medical image registration. The registration is formulated as a minimization problem involving robust estimators. We propose an efficient hierarchical optimization framework which is both multiresolution and multigrid. An anatomical segmentation of the cortex is introduced in the adaptive partitioning of the volume on which the multigrid minimization is based. This allows to limit the estimation to the areas of interest, to accelerate the algorithm, and to refine the estimation in specified areas. Furthermore we introduce a methodology to constrain the registration with landmarks such as anatomical structures. The performances of this method are objectively evaluated on simulated data and its benefits are demonstrated on a large database of real acquisitions.

Presents a cooperative strategy between volumetric registration and segmentation. The segmentation method is based on the level set formalism. Starting from an initial position, a closed 3D surface propagates towards the desired boundaries through the evolution of a 4D implicit function. We show that the number of iterations required for convergence is significantly reduced by using a registration process to initialize the surface. Furthermore it makes the segmentation fully automatic. The registration is achieved through a robust multiresolution and multigrid minimization scheme appropriate to our problem. In addition, a bidirectional propagation force depending on local intensity values has been designed for the evolution of the surface. Finally, an adaptive iteration step is automatically computed at each iteration in order to improve the robustness and the efficiency of the algorithm. Results on volumetric brain MR images are presented and discussed.

Systematic mapping of the variability in cortical sulcal anatomy is an area of increasing interest which presents numerous methodological challenges. To address these issues, we have implemented sulcal extraction and assisted labeling (SEAL) to automatically extract the two-dimensional (2-D) surface ribbons that represent the median axis of cerebral sulci and to neuroanatomically label these entities. To encode the extracted three-dimensional (3-D) cortical sulcal schematic topography (CSST) we define a relational graph structure composed of two main features: vertices (representing sulci) and arcs (representing the relationships between sulci). Vertices contain a parametric representation of the surface ribbon buried within the sulcus. Points on this surface are expressed in stereotaxic coordinates (i.e., with respect to a standardized brain coordinate system). For each of these vertices, we store length, depth, and orientation as well as anatomical attributes (e.g., hemisphere, lobe, sulcus type, etc.). Each arc stores the 3-D location of the junction between sulci as well as a list of its connecting sulci. Sulcal labeling is performed semiautomatically by selecting a sulcal entity in the CSST and selecting from a menu of candidate sulcus names. In order to help the user in the labeling task, the menu is restricted to the most likely candidates by using priors for the expected sulcal spatial distribution. These priors, i.e., sulcal probabilistic maps, were created from the spatial distribution of 34 sulci traced manually on 36 different subjects. Given these spatial probability maps, the user is provided with the likelihood that the selected entity belongs to a particular sulcus. The cortical structure representation obtained by SEAL is suitable to extract statistical information about both the spatial and the structural composition of the cerebral cortical topography. This methodology allows for the iterative construction of a successively more complete statistical models of the cerebral topography containing spatial distributions of the most important structures, their morphometrics, and their structural components.

We propose a method for the 3D segmentation and representation of cortical folds with a special emphasis on the cortical sulci. These cortical structures are represented using "active ribbons". Active ribbons are built from active surfaces, which represent the median surface of a particular sulcus filled by CSF. Sulci modeling is obtained from MRI acquisitions (usually T1 images). The segmentation is performed using an automatic labeling procedure to separate gyri from sulci based on curvature analysis of the different iso-intensity surfaces of the original MRI volume. The outer parts of the sulci are used to initialize the convergence of the active ribbon from the outer parts of the brain to the interior. This procedure has two advantages: first, it permits the labeling of voxels belonging to sulci on the external part of the brain as well as on the inside (which is often the hardest point) and secondly, this segmentation allows 3D visualization of the sulci in the MRI volumetric environment as well as showing the sophisticated shapes of the cortical structures by means of isolated surfaces. Active ribbons can be used to study the complicated shape of the cortical anatomy, to model the variability of these structures in shape and position, to assist nonlinear registrations of human brains by locally controlling the warping procedure, to map brain neurophysiological functions into morphology or even to select the trajectory of an intra-sulci (virtual) endoscope.

Previous extensions of two-dimensional ultrasonic imaging to three dimensions used lattice diagrams which give measurement information, but no anatomic detail. The authors conducted three sets of experiments to test the hypothesis that complete acoustic backscatter data should be retained to produce useful information about heart structure and function. First, in vitro compound B-scans were taken under ideal conditions; second, in vitro rotating conventional sector scans were taken to test clinically applicable methods; and third, clinical in vivo rotating conventional sector scans were taken of a human volunteer. It is concluded that the resulting images show details of cardiac anatomy and have great clinical promise. Interactive analysis and surface and volume displays give context and perspective information which should improve diagnostic accuracy, communication with noncardiologists and yield more precise measurements of anatomical structure and function.>