Dense deformation field estimation for atlas registration using the active contour framework

A key research area in computer vision is image segmentation. Image segmentation aims at extracting objects of interest in images or video sequences. These objects contain relevant information for a given application. For example, a video surveillance application generally requires to extract moving objects (vehicles, persons or animals) from a sequence of images in order to check that their path stays conformed to the regulation rules set for the observed scene. Image segmentation is not an easy task. In many applications, the contours of the objects of interest are difficult to delineate, even manually. The problems linked to segmentation are often due to low contrast, fuzzy contours or too similar intensities with adjacent objects. In some cases, the objects to be extracted have no real contours in the image. This kind of objects is called virtual objects. Virtual objects appear especially in medical applications. To draw them, medical experts usually estimate their position from surrounding objects. The problems related to image segmentation can be greatly simplified with information known in advance on the objects to be extracted (the prior knowledge). A widely used method consists to extract the needed prior knowledge from a reference image often called atlas. The goal of the atlas is to describe the image to be segmented like a map would describe the components of a geographical area. An atlas can contain three types of information on each object being part of the image: an estimation of its position in the image, a description of its shape and texture, and the features of its adjacent objects. The atlas-based segmentation method is rather used when the atlas can characterize a range of images. This method is thus especially adapted to medical images due to the existing consistency between anatomical structures of same type. There exist two types of atlas: the determinist atlas and the statistical atlas. The determinist atlas is an image which has been selected or computed, to be the most representative of an image category to be segmented. This image is called intensity atlas. The contours of the objects of interest (the objects to be extracted in images of the same type) have been traced manually on the intensity atlas, or by using a semi-automatic method. A label is often attributed to each one of these objects in order to differentiate them. In this way, we obtain a labeled version of the atlas called labeled atlas. The statistical atlas is an atlas created from a database of images in order to be the most representative of a certain type of images to be segmented. In this atlas, the position and the features of the objects of interest depend on statistical measures. In this thesis, we are focused on the use of determinist atlases for image segmentation. The segmentation process with a determinist atlas consists to deform the objects delineated in the atlas in order to better align them with their corresponding objects in the image to be segmented. To perform this task, we have distinguished two types of approaches in the literature. The first approach consists to reduce the segmentation problem in an image registration problem. First of all, a dense deformation field that registers (i.e. puts in point-to-point spatial correspondence) the atlas to the image to be segmented, is explicitly computed. Then, this transformation is used to project the assigned labels onto each atlas structure on the image to be segmented. The advantage of this approach is that the deformation field computed from the registration of visible contours allows to easily estimate the position of virtual objects or objects with fuzzy contours. However, the methods currently used for the atlas registration are often only based on the intensity atlas. That means that they do not exploit the object-based information that can be obtained by combining the intensity atlas with its labeled version. In the second approach, the atlas contours selected by the labeled atlas are directly deformed without using a geometrical deformation. For that, this approach is based on matching contour techniques, generally called deformable models. In this thesis, we are interested to a particular type of deformable models, which are the active contour segmentation models. The advantage of the active contour method is that this segmentation technique has been designed to exploit the image information directly linked to the object to be delineated. By using object-based information, active contour models are frequently able to extract regions where the atlas-based segmentation method by registration fails. On the other hand, the result of this local segmentation method is very sensitive to the initial atlas contour position regarding to the target contours. On the other hand, this local segmentation method is very sensitive to the initial position of the atlas contours: the closer they are to the contours to be detected, the more robust the active contour-based segmentation will be. Besides, this segmentation technique needs prior shape models to be able to estimate the position of virtual objects. The main objective of this thesis is to design an algorithm for atlas-based segmentation which combines the advantages of the dense deformation field computed by the registration algorithms, with local segmentation constraints coming from the active contour framework. This implies to design a model where the registration and segmentation by active contours are jointly performed. The atlas registration algorithm that we propose is based on a formulation allowing the integration of any segmentation or contour regularization forces derived from the theory of the active contours in a non parametric registration process. Our algorithm led us to introduce the concept of hierarchical atlas registration. Its principle is that the registration of the main image objects helps the registration of depending objects. This allows to bring progressively the atlas contours closer to their target and thus, to limit the risk to be stuck in a local minimum. Our model had been designed to be easily adaptable to various types of segmentation problems. At the end of the thesis, we present several examples of atlas registration applications in medical imaging. These applications highlight the integration of manual constraints in an atlas registration process, the modeling of a tumor growth in the atlas, the labelization of the thalamus for a statistical study on neuronal connections, the localization of the subthalamic nucleus (STN) for deep brain stimulation (DBS) and the compensation of intra-operative brain shift for neuronavigation systems

Geometry Processing of Conventionally Produced Mouse Brain Slice Images

Brain mapping research in most neuroanatomical laboratories relies on conventional processing techniques, which often introduce histological artifacts such as tissue tears and tissue loss. In this paper we present techniques and algorithms for automatic registration and 3D reconstruction of conventionally produced mouse brain slices in a standardized atlas space. This is achieved first by constructing a virtual 3D mouse brain model from annotated slices of Allen Reference Atlas (ARA). Virtual re-slicing of the reconstructed model generates ARA-based slice images corresponding to the microscopic images of histological brain sections. These image pairs are aligned using a geometric approach through contour images. Histological artifacts in the microscopic images are detected and removed using Constrained Delaunay Triangulation before performing global alignment. Finally, non-linear registration is performed by solving Laplace's equation with Dirichlet boundary conditions. Our methods provide significant improvements over previously reported registration techniques for the tested slices in 3D space, especially on slices with significant histological artifacts. Further, as an application we count the number of neurons in various anatomical regions using a dataset of 51 microscopic slices from a single mouse brain. This work represents a significant contribution to this subfield of neuroscience as it provides tools to neuroanatomist for analyzing and processing histological data.Comment: 14 pages, 11 figure

Imaging Informatics and the Human Brain Project: the Role of Structure, Review

(no abstract

Automated detection of brain abnormalities in neonatal hypoxia ischemic injury from MR images.

We compared the efficacy of three automated brain injury detection methods, namely symmetry-integrated region growing (SIRG), hierarchical region splitting (HRS) and modified watershed segmentation (MWS) in human and animal magnetic resonance imaging (MRI) datasets for the detection of hypoxic ischemic injuries (HIIs). Diffusion weighted imaging (DWI, 1.5T) data from neonatal arterial ischemic stroke (AIS) patients, as well as T2-weighted imaging (T2WI, 11.7T, 4.7T) at seven different time-points (1, 4, 7, 10, 17, 24 and 31 days post HII) in rat-pup model of hypoxic ischemic injury were used to assess the temporal efficacy of our computational approaches. Sensitivity, specificity, and similarity were used as performance metrics based on manual ('gold standard') injury detection to quantify comparisons. When compared to the manual gold standard, automated injury location results from SIRG performed the best in 62% of the data, while 29% for HRS and 9% for MWS. Injury severity detection revealed that SIRG performed the best in 67% cases while 33% for HRS. Prior information is required by HRS and MWS, but not by SIRG. However, SIRG is sensitive to parameter-tuning, while HRS and MWS are not. Among these methods, SIRG performs the best in detecting lesion volumes; HRS is the most robust, while MWS lags behind in both respects

Dense Deformation Field Estimation for Atlas Registration using the Active Contour Framework

In this paper, we propose a new paradigm to carry outthe registration task with a dense deformation fieldderived from the optical flow model and the activecontour method. The proposed framework merges differenttasks such as segmentation, regularization, incorporationof prior knowledge and registration into a singleframework. The active contour model is at the core of ourframework even if it is used in a different way than thestandard approaches. Indeed, active contours are awell-known technique for image segmentation. Thistechnique consists in finding the curve which minimizesan energy functional designed to be minimal when thecurve has reached the object contours. That way, we getaccurate and smooth segmentation results. So far, theactive contour model has been used to segment objectslying in images from boundary-based, region-based orshape-based information. Our registration technique willprofit of all these families of active contours todetermine a dense deformation field defined on the wholeimage. A well-suited application of our model is theatlas registration in medical imaging which consists inautomatically delineating anatomical structures. Wepresent results on 2D synthetic images to show theperformances of our non rigid deformation field based ona natural registration term. We also present registrationresults on real 3D medical data with a large spaceoccupying tumor substantially deforming surroundingstructures, which constitutes a high challenging problem

Contour-Driven Atlas-Based Segmentation

We propose new methods for automatic segmentation of images based on an atlas of manually labeled scans and contours in the image. First, we introduce a Bayesian framework for creating initial label maps from manually annotated training images. Within this framework, we model various registration- and patch-based segmentation techniques by changing the deformation field prior. Second, we perform contour-driven regression on the created label maps to refine the segmentation. Image contours and image parcellations give rise to non-stationary kernel functions that model the relationship between image locations. Setting the kernel to the covariance function in a Gaussian process establishes a distribution over label maps supported by image structures. Maximum a posteriori estimation of the distribution over label maps conditioned on the outcome of the atlas-based segmentation yields the refined segmentation. We evaluate the segmentation in two clinical applications: the segmentation of parotid glands in head and neck CT scans and the segmentation of the left atrium in cardiac MR angiography images

Deeply-Supervised CNN for Prostate Segmentation

Prostate segmentation from Magnetic Resonance (MR) images plays an important role in image guided interven- tion. However, the lack of clear boundary specifically at the apex and base, and huge variation of shape and texture between the images from different patients make the task very challenging. To overcome these problems, in this paper, we propose a deeply supervised convolutional neural network (CNN) utilizing the convolutional information to accurately segment the prostate from MR images. The proposed model can effectively detect the prostate region with additional deeply supervised layers compared with other approaches. Since some information will be abandoned after convolution, it is necessary to pass the features extracted from early stages to later stages. The experimental results show that significant segmentation accuracy improvement has been achieved by our proposed method compared to other reported approaches.Comment: Due to a crucial sign error in equation