Search PubMed⌕ Search

Biomedical subjects

Daniel Rueckert

Publications and source records attributed to Daniel Rueckert.

17 recordsLinked to original sources

Automatic anatomical brain MRI segmentation combining label propagation and decision fusion.

Regions in three-dimensional magnetic resonance (MR) brain images can be classified using protocols for manually segmenting and labeling structures. For large cohorts, time and expertise requirements make this approach impractical. To achieve automation, an individual segmentation can be propagated to another individual using an anatomical correspondence estimate relating the atlas image to the target image. The accuracy of the resulting target labeling has been limited but can potentially be improved by combining multiple segmentations using decision fusion. We studied segmentation propagation and decision fusion on 30 normal brain MR images, which had been manually segmented into 67 structures. Correspondence estimates were established by nonrigid registration using free-form deformations. Both direct label propagation and an indirect approach were tested. Individual propagations showed an average similarity index (SI) of 0.754+/-0.016 against manual segmentations. Decision fusion using 29 input segmentations increased SI to 0.836+/-0.009. For indirect propagation of a single source via 27 intermediate images, SI was 0.779+/-0.013. We also studied the effect of the decision fusion procedure using a numerical simulation with synthetic input data. The results helped to formulate a model that predicts the quality improvement of fused brain segmentations based on the number of individual propagated segmentations combined. We demonstrate a practicable procedure that exceeds the accuracy of previous automatic methods and can compete with manual delineations.

Adult↗

Cerebral atrophy measurements using Jacobian integration: comparison with the boundary shift integral.

We compared two methods of measuring cerebral atrophy in a cohort of 38 clinically probable Alzheimer's disease (AD) subjects and 22 age-matched normal controls, using metrics of zero atrophy, consistency, scaled atrophy and AD/control group separation. The two methods compared were the boundary shift integral (BSI) and a technique based on the integration of Jacobian determinants from non-rigid registration. For each subject, we used two volumetric magnetic resonance (MR) scans at baseline and a third obtained 1 year later. The case of zero atrophy was established by registering the same-day baseline scan pair, which should approximate zero change. Consistency was established by registering the 1-year follow-up scan to each of the baseline scans, giving two measurements of atrophy that should be very similar, while scaled atrophy was established by reducing one of the same-day scans by a fixed amount, and rigidly registering this to the other same-day scan. Group separation was ascertained by calculating atrophy rates over the two 1-year measures for the control and AD subjects. The results showed the Jacobian integration technique was significantly more accurate in calculating scaled atrophy (P < 0.001) and was able to distinguish between control and AD subjects more clearly (P < 0.01).

Aged↗

Abnormal deep grey matter development following preterm birth detected using deformation-based morphometry.

Preterm birth is a leading risk factor for neurodevelopmental and cognitive impairment in childhood and adolescence. The most common known cerebral abnormality among preterm infants at term equivalent age is a diffuse white matter abnormality seen on magnetic resonance (MR) images. It occurs with a similar prevalence to subsequent impairment, but its effect on developing neural systems is unknown. MR images were obtained at term equivalent age from 62 infants born at 24-33 completed weeks gestation and 12 term born controls. Tissue damage was quantified using diffusion-weighted imaging, and deformation-based morphometry was used to make a non-subjective survey of the whole brain to identify significant cerebral morphological alterations associated with preterm birth and with diffuse white matter injury. Preterm infants at term equivalent age had reduced thalamic and lentiform volumes without evidence of acute injury in these regions (t = 5.81, P < 0.05), and these alterations were more marked with increasing prematurity (t = 7.13, P < 0.05 for infants born at less than 28 weeks) and in infants with diffuse white matter injury (t = 6.43, P < 0.05). The identification of deep grey matter growth failure in association with diffuse white matter injury suggests that white matter injury is not an isolated phenomenon, but rather, it is associated with the maldevelopment of remote structures. This could be mediated by a disturbance to corticothalamic connectivity during a critical period in cerebral development. Deformation-based morphometry is a powerful tool for modelling the developing brain in health and disease, and can be used to test putative aetiological factors for injury.

Adolescent↗

Tract-based spatial statistics: voxelwise analysis of multi-subject diffusion data.

There has been much recent interest in using magnetic resonance diffusion imaging to provide information about anatomical connectivity in the brain, by measuring the anisotropic diffusion of water in white matter tracts. One of the measures most commonly derived from diffusion data is fractional anisotropy (FA), which quantifies how strongly directional the local tract structure is. Many imaging studies are starting to use FA images in voxelwise statistical analyses, in order to localise brain changes related to development, degeneration and disease. However, optimal analysis is compromised by the use of standard registration algorithms; there has not to date been a satisfactory solution to the question of how to align FA images from multiple subjects in a way that allows for valid conclusions to be drawn from the subsequent voxelwise analysis. Furthermore, the arbitrariness of the choice of spatial smoothing extent has not yet been resolved. In this paper, we present a new method that aims to solve these issues via (a) carefully tuned non-linear registration, followed by (b) projection onto an alignment-invariant tract representation (the "mean FA skeleton"). We refer to this new approach as Tract-Based Spatial Statistics (TBSS). TBSS aims to improve the sensitivity, objectivity and interpretability of analysis of multi-subject diffusion imaging studies. We describe TBSS in detail and present example TBSS results from several diffusion imaging studies.

Amyotrophic Lateral Sclerosis↗

Beyond the g-factor limit in sensitivity encoding using joint histogram entropy.

The maximum practical speed-up that can be achieved using parallel imaging methods is widely accepted to be limited by g-factor noise. An approximate expression for the g-factor noise as a function of the principal eigenvector of the inverse sensitivity matrix is derived. This formulation allows g-factor enhanced noise to be reduced by a constrained optimization procedure with joint image histogram entropy between a reference image and a SENSE image as an image quality metric. The reference image does not need to have identical resolution or contrast. The reference image may also be used for coil calibration. The limits of the method are explored using simulated and real array coil data with high g-factor using a variety of contrast and resolution combinations. The method preserves image structure, contrast, and lesions even when these were not observable in the reference data. In all cases g-factor was dramatically reduced.

Brain Mapping↗

Automatic quantification of changes in bone in serial MR images of joints.

Recent innovations in drug therapies have made it highly desirable to obtain sensitive biomarkers of disease progression that can be used to quantify the performance of candidate disease modifying drugs. In order to measure potential image-based biomarkers of disease progression in an experimental model of rheumatoid arthritis (RA), we present two different methods to automatically quantify changes in a bone in in-vivo serial magnetic resonance (MR) images from the model. Both methods are based on rigid and nonrigid image registration to perform the analysis. The first method uses segmentation propagation to delineate a bone from the serial MR images giving a global measure of temporal changes in bone volume. The second method uses rigid body registration to determine intensity change within a bone, and then maps these into a reference coordinate system using nonrigid registration. This gives a local measure of temporal changes in bone lesion volume. We detected significant temporal changes in local bone lesion volume in five out of eight identified candidate bone lesion regions, and significant difference in local bone lesion volume between male and female subjects in three out of eight candidate bone lesion regions. But the global bone volume was found to be fluctuating over time. Finally, we compare our findings with histology of the subjects and the manual segmentation of bone lesions.

Algorithms↗

Generalised overlap measures for assessment of pairwise and groupwise image registration and segmentation.

Effective validation techniques are an essential pre-requisite for segmentation and non-rigid registration techniques to enter clinical use. These algorithms can be evaluated by calculating the overlap of corresponding test and gold-standard regions. Common overlap measures compare pairs of binary labels but it is now common for multiple labels to exist and for fractional (partial volume) labels to be used to describe multiple tissue types contributing to a single voxel. Evaluation studies may involve multiple image pairs. In this paper we use results from fuzzy set theory and fuzzy morphology to extend the definitions of existing overlap measures to accommodate multiple fractional labels. Simple formulas are provided which define single figures of merit to quantify the total overlap for ensembles of pairwise or groupwise label comparisons. A quantitative link between overlap and registration error is established by defining the overlap tolerance. Experiments are performed on publicly available labeled brain data to demonstrate the new measures in a comparison of pairwise and groupwise registration.

Algorithms↗

Construction of a 4D statistical atlas of the cardiac anatomy and its use in classification.

In this paper we present a novel method for building a 4D statistical atlas describing the cardiac anatomy and how the cardiac anatomy changes during the cardiac cycle. The method divides the distribution space of cardiac shapes into two subspaces. One distribution subspace accounts for changes in cardiac shape caused by inter-subject variability. The second distribution subspace accounts for changes in cardiac shape caused by deformation during the cardiac cycle (i.e. intra-subject variability). Principal component analysis (PCA) have been performed in order to calculate the most significant modes of variation of each distribution subspace. During the construction of the statistical atlas we eliminate the need for manual landmarking of the cardiac images by using a non-rigid surface registration algorithm to propagate a set of pseudo-landmarks from an automatically landmarked atlas to each frame of all the image sequences. In order to build the atlas we have used 26 cardiac image sequences from healthy volunteers. We show how the resulting statistical atlas can be used to differentiate between cardiac image sequences from patients with hypertrophic cardiomyopathy and normal subjects.

Algorithms↗

Spatio-temporal free-form registration of cardiac MR image sequences.

In this paper we present two registration algorithms for the spatio-temporal alignment of cardiac MR image sequences. Both algorithms have the ability to correct spatial misalignment between the images sequences caused by global and local shape differences. In addition, they have the ability to correct temporal misalignment caused by differences in the length of the cardiac cycles and by differences in the dynamic properties of the hearts. The algorithms use a 4D deformable transformation model which is separated into spatial and temporal components. The first registration algorithm optimizes the spatial and temporal transformation models simultaneously, while the second registration algorithm optimizes the temporal transformation component before optimizing the spatial component. For the evaluation of the spatio-temporal registration methods we have acquired 15 MR image sequences from healthy volunteers. The registration methods were quantitatively evaluated by measuring the overlap and surface distance of anatomical regions and qualitatively by visual inspection. The results demonstrate that a significant improvement in the alignment of the image sequences is achieved by the use of the deformable spatio-temporal transformation model. We demonstrate the use of the method for the construction of a probabilistic MR cardiac atlas representing the anatomy and function of a healthy heart.

Algorithms↗

Fast generation of digitally reconstructed radiographs using attenuation fields with application to 2D-3D image registration.

Generation of digitally reconstructed radiographs (DRRs) is computationally expensive and is typically the rate-limiting step in the execution time of intensity-based two-dimensional to three-dimensional (2D-3D) registration algorithms. We address this computational issue by extending the technique of light field rendering from the computer graphics community. The extension of light fields, which we call attenuation fields (AFs), allows most of the DRR computation to be performed in a preprocessing step; after this precomputation step, DRRs can be generated substantially faster than with conventional ray casting. We derive expressions for the physical sizes of the two planes of an AF necessary to generate DRRs for a given X-ray camera geometry and all possible object motion within a specified range. Because an AF is a ray-based data structure, it is substantially more memory efficient than a huge table of precomputed DRRs because it eliminates the redundancy of replicated rays. Nonetheless, an AF can require substantial memory, which we address by compressing it using vector quantization. We compare DRRs generated using AFs (AF-DRRs) to those generated using ray casting (RC-DRRs) for a typical C-arm geometry and computed tomography images of several anatomic regions. They are quantitatively very similar: the median peak signal-to-noise ratio of AF-DRRs versus RC-DRRs is greater than 43 dB in all cases. We perform intensity-based 2D-3D registration using AF-DRRs and RC-DRRs and evaluate registration accuracy using gold-standard clinical spine image data from four patients. The registration accuracy and robustness of the two methods is virtually identical whereas the execution speed using AF-DRRs is an order of magnitude faster.

Algorithms↗

Segmentation of 4D cardiac MR images using a probabilistic atlas and the EM algorithm.

In this paper an automatic atlas-based segmentation algorithm for 4D cardiac MR images is proposed. The algorithm is based on the 4D extension of the expectation maximisation (EM) algorithm. The EM algorithm uses a 4D probabilistic cardiac atlas to estimate the initial model parameters and to integrate a priori information into the classification process. The probabilistic cardiac atlas has been constructed from the manual segmentations of 3D cardiac image sequences of 14 healthy volunteers. It provides space and time-varying probability maps for the left and right ventricles, the myocardium, and background structures such as the liver, stomach, lungs and skin. In addition to using the probabilistic cardiac atlas as a priori information, the segmentation algorithm incorporates spatial and temporal contextual information by using 4D Markov Random Fields. After the classification, the largest connected component of each structure is extracted using a global connectivity filter which improves the results significantly, especially for the myocardium. Validation against manual segmentations and computation of the correlation between manual and automatic segmentation on 249 3D volumes were calculated. We used the 'leave one out' test where the image set to be segmented was not used in the construction of its corresponding atlas. Results show that the procedure can successfully segment the left ventricle (LV) (r = 0.96), myocardium (r = 0.92) and right ventricle (r = 0.92). In addition, 4D images from 10 patients with hypertrophic cardiomyopathy were also manually and automatically segmented yielding a good correlation in the volumes of the LV (r = 0.93) and myocardium (0.94) when the atlas constructed with volunteers is blurred.

Algorithms↗

Analysis of 3-D myocardial motion in tagged MR images using nonrigid image registration.

Tagged magnetic resonance imaging (MRI) is unique in its ability to noninvasively image the motion and deformation of the heart in vivo, but one of the fundamental reasons limiting its use in the clinical environment is the absence of automated tools to derive clinically useful information from tagged MR images. In this paper, we present a novel and fully automated technique based on nonrigid image registration using multilevel free-form deformations (MFFDs) for the analysis of myocardial motion using tagged MRI. The novel aspect of our technique is its integrated nature for tag localization and deformation field reconstruction using image registration and voxel based similarity measures. To extract the motion field within the myocardium during systole we register a sequence of images taken during systole to a set of reference images taken at end-diastole, maximizing the normalized mutual information between the images. We use both short-axis and long-axis images of the heart to estimate the full four-dimensional motion field within the myocardium. We also present validation results from data acquired from twelve volunteers.

Algorithms↗

Construction of a statistical model for cardiac motion analysis using nonrigid image registration.

In this paper we present a new technique for tracking the movement of the myocardium using a statistical model derived from the motion fields in the hearts of several healthy volunteers. To build the statistical model we tracked the motion of the myocardium in 17 volunteers using a nonrigid registration technique based on free-form deformations and mapped the motion fields obtained into a common reference coordinate system. A principal component analysis (PCA) was then performed on the motion fields to extract the major modes of variation in the fields between the successive time frames. The modes of variation obtained were then used to parametrize the free-form deformations and build our statistical model. The results of using our model to track the motion of the heart in normal volunteers are also presented.

Algorithms↗

Automatic construction of 3-D statistical deformation models of the brain using nonrigid registration.

In this paper, we show how the concept of statistical deformation models (SDMs) can be used for the construction of average models of the anatomy and their variability. SDMs are built by performing a statistical analysis of the deformations required to map anatomical features in one subject into the corresponding features in another subject. The concept of SDMs is similar to statistical shape models (SSMs) which capture statistical information about shapes across a population, but offers several advantages over SSMs. First, SDMs can be constructed directly from images such as three-dimensional (3-D) magnetic resonance (MR) or computer tomography volumes without the need for segmentation which is usually a prerequisite for the construction of SSMs. Instead, a nonrigid registration algorithm based on free-form deformations and normalized mutual information is used to compute the deformations required to establish dense correspondences between the reference subject and the subjects in the population class under investigation. Second, SDMs allow the construction of an atlas of the average anatomy as well as its variability across a population of subjects. Finally, SDMs take the 3-D nature of the underlying anatomy into account by analysing dense 3-D deformation fields rather than only information about the surface shape of anatomical structures. We show results for the construction of anatomical models of the brain from the MR images of 25 different subjects. The correspondences obtained by the nonrigid registration are evaluated using anatomical landmark locations and show an average error of 1.40 mm at these anatomical landmark positions. We also demonstrate that SDMs can be constructed so as to minimize the bias toward the chosen reference subject.

Algorithms↗

Registration and tracking to integrate X-ray and MR images in an XMR facility.

We describe a registration and tracking technique to integrate cardiac X-ray images and cardiac magnetic resonance (MR) images acquired from a combined X-ray and MR interventional suite (XMR). Optical tracking is used to determine the transformation matrices relating MR image coordinates and X-ray image coordinates. Calibration of X-ray projection geometry and tracking of the X-ray C-arm and table enable three-dimensional (3-D) reconstruction of vessel centerlines and catheters from bi-plane X-ray views. We can, therefore, combine single X-ray projection images with registered projection MR images from a volume acquisition, and we can also display 3-D reconstructions of catheters within a 3-D or multi-slice MR volume. Registration errors were assessed using phantom experiments. Errors in the combined projection images (two-dimensional target registration error--TRE) were found to be 2.4 to 4.2 mm, and the errors in the integrated volume representation (3-D TRE) were found to be 4.6 to 5.1 mm. These errors are clinically acceptable for alignment of images of the great vessels and the chambers of the heart. Results are shown for two patients. The first involves overlay of a catheter used for invasive pressure measurements on an MR volume that provides anatomical context. The second involves overlay of invasive electrode catheters (including a basket catheter) on a tagged MR volume in order to relate electrophysiology to myocardial motion in a patient with an arrhythmia. Visual assessment of these results suggests the errors were of a similar magnitude to those obtained in the phantom measurements.

Adolescent↗

Automatic construction of multiple-object three-dimensional statistical shape models: application to cardiac modeling.

A novel method is introduced for the generation of landmarks for three-dimensional (3-D) shapes and the construction of the corresponding 3-D statistical shape models. Automatic landmarking of a set of manual segmentations from a class of shapes is achieved by 1) construction of an atlas of the class, 2) automatic extraction of the landmarks from the atlas, and 3) subsequent propagation of these landmarks to each example shape via a volumetric nonrigid registration technique using multiresolution B-spline deformations. This approach presents some advantages over previously published methods: it can treat multiple-part structures and requires less restrictive assumptions on the structure's topology. In this paper, we address the problem of building a 3-D statistical shape model of the left and right ventricle of the heart from 3-D magnetic resonance images. The average accuracy in landmark propagation is shown to be below 2.2 mm. This application demonstrates the robustness and accuracy of the method in the presence of large shape variability and multiple objects.

Algorithms↗