We present a novel approach for joint clustering and point-by-point mapping of white matter fiber pathways. Knowledge of the point correspondence along the fiber pathways is not only necessary for accurate clustering of the trajectories into fiber bundles, but also crucial for any tract-oriented quantitative analysis. We employ an expectation-maximization (EM) algorithm to cluster the trajectories in a gamma mixture model context. The result of clustering is the probabilistic assignment of the fiber trajectories to each cluster, an estimate of the cluster parameters, i.e. spatial mean and variance, and point correspondences. The fiber bundles are modeled by the mean trajectory and its spatial variation. Point-by-point correspondence of the trajectories within a bundle is obtained by constructing a distance map and a label map from each cluster center at every iteration of the EM algorithm. This offers a time-efficient alternative to pairwise curve matching of all trajectories with respect to each cluster center. The proposed method has the potential to benefit from an anatomical atlas of fiber tracts by incorporating it as prior information in the EM algorithm. The algorithm is also capable of handling outliers in a principled way. The presented results confirm the efficiency and effectiveness of the proposed framework for quantitative analysis of diffusion tensor MRI.
BACKGROUND: Use of implantable cardiac defibrillators (ICDs) in children and patients with congenital heart disease is complicated by body size and anatomy. A variety of creative implantation techniques has been used empirically in these groups on an ad hoc basis.
OBJECTIVE: To rationalize ICD placement in special populations, we used subject-specific, image-based finite element models (FEMs) to compare electric fields and expected defibrillation thresholds (DFTs) using standard and novel electrode configurations.
METHODS: FEMs were created by segmenting normal torso computed tomography scans of subjects ages 2, 10, and 29 years and 1 adult with congenital heart disease into tissue compartments, meshing, and assigning tissue conductivities. The FEMs were modified by interactive placement of ICD electrode models in clinically relevant electrode configurations, and metrics of relative defibrillation safety and efficacy were calculated.
RESULTS: Predicted DFTs for standard transvenous configurations were comparable with published results. Although transvenous systems generally predicted lower DFTs, a variety of extracardiac orientations were also predicted to be comparably effective in children and adults. Significant trend effects on DFTs were associated with body size and electrode length. In many situations, small alterations in electrode placement and patient anatomy resulted in significant variation of predicted DFT. We also show patient-specific use of this technique for optimization of electrode placement.
CONCLUSION: Image-based FEMs allow predictive modeling of defibrillation scenarios and predict large changes in DFTs with clinically relevant variations of electrode placement. Extracardiac ICDs are predicted to be effective in both children and adults. This approach may aid both ICD development and patient-specific optimization of electrode placement. Further development and validation are needed for clinical or industrial utilization.
The theories of signal sampling, filter banks, wavelets, and "overcomplete wavelets" are well established for the Euclidean spaces and are widely used in the processing and analysis of images. While recent advances have extended some filtering methods to spherical images, many key challenges remain. In this paper, we develop theoretical conditions for the invertibility of filter banks under continuous spherical convolution. Furthermore, we present an analogue of the Papoulis generalized sampling theorem on the 2-Sphere. We use the theoretical results to establish a general framework for the design of invertible filter banks on the sphere and demonstrate the approach with examples of self-invertible spherical wavelets and steerable pyramids. We conclude by examining the use of a self-invertible spherical steerable pyramid in a denoising experiment and discussing the computational complexity of the filtering framework.
We created an image-guided robot system to assist with skull base drilling by integrating a robot, a commercial navigation system, and an open source visualization platform. The objective of this procedure is to create a cavity in the skull base to allow access for neurosurgical interventions. The motivation for introducing an image-guided robot is to improve safety by preventing the surgeon from accidentally damaging critical structures during the drilling procedure. Our approach is to attach the cutting tool to the robot end-effector and operate the robot in a cooperative control mode, where robot motion is determined from the forces and torques applied by the surgeon. We employ "virtual fixtures" to constrain the motion of the cutting tool so that it remains in the safe zone that was defined on a preoperative CT scan. This paper presents the system design and the results of phantom and cadaveric experiments. Both experiments have demonstrated the feasibility of the system, with average overcut error at about 1 mm and maximum errors at 2.5 mm.
The level set method is a popular technique used in medical image segmentation; however, the numerics involved make its use cumbersome. This paper proposes an approximate level set scheme that removes much of the computational burden while maintaining accuracy. Abandoning a floating point representation for the signed distance function, we use integral values to represent the signed distance function. For the cases of 2D and 3D, we detail rules governing the evolution and maintenance of these three regions. Arbitrary energies can be implemented in the framework. This scheme has several desirable properties: computations are only performed along the zero level set; the approximate distance function requires only a few simple integer comparisons for maintenance; smoothness regularization involves only a few integer calculations and may be handled apart from the energy itself; the zero level set is represented exactly removing the need for interpolation off the interface; and evolutions proceed on the order of milliseconds per iteration on conventional uniprocessor workstations. To highlight its accuracy, flexibility and speed, we demonstrate the technique on intensity-based segmentations under various statistical metrics. Results for 3D imagery show the technique is fast even for image volumes.
The hippocampus is known to be vulnerable to hypoxia, stress, and undernutrition, all likely to be present in fetal intrauterine growth restriction (IUGR). The effect of IUGR in preterm infants on the hippocampus was studied using 3D magnetic resonance imaging at term-equivalent age Thirteen preterm infants born with IUGR after placental insufficiency were compared with 13 infants with normal intrauterine growth age matched for gestational age. The hippocampal structural differences were defined using voxel-based morphometry and manual segmentation. The specific neurobehavioral function was evaluated by the Assessment of Preterm Infants' Behavior at term and at 24 mo of corrected age by a Bayley Scales of Infant and Toddler Development. Voxel-based morphometry detected significant gray matter volume differences in the hippocampus between the two groups. This finding was confirmed by manual segmentation of the hippocampus with a reduction of hippocampal volume after IUGR. The hippocampal volume reduction was further associated with functional behavioral differences at term-equivalent age in all six subdomains of the Assessment of Preterm Infants' Behavior but not at 24 mo of corrected age. We conclude that hippocampal development in IUGR is altered and might result from a combination of maternal corticosteroid hormone exposure, hypoxemia, and micronutrient deficiency.
In this paper, we have developed a geometric-based scaling model that describes the adsorption of diblock copolymer chains from good solvents and theta-solvents onto reactive surfaces of varying curvatures. To evaluate the impact of particle size on the adsorption process, we probed the adsorption of poly(styrene-b-methymethacrylate) (PS-PMMA) diblock copolymers from solvents with different degrees of selectivity on aluminum oxide (Al(2)O(3)) surfaces belonging to particles of different sizes. When the adsorbed PMMA layer is dense enough (in the case of a theta-solvent for the PMMA block), our results show good correlation between the theory and experimental results, pointing to the formation of a PMMA adsorption layer and a brushlike PS layer. Conversely, when adsorption occurs from a nonpreferential solvent, particularly on particles with high curvature, the PMMA adsorption layer at the surface becomes less dense and the grafted PS moiety exhibits a transitional morphology consisting of several layers of increasingly sparsely spaced blobs.
BACKGROUND AND PURPOSE: Lesion volume change (LVC) assessment is essential in monitoring MS progression. LVC is usually measured by independently segmenting serial MR imaging examinations. Subtraction imaging has been proposed for improved visualization and characterization of lesion change. We compare segmentation of subtraction images (SSEG) with serial single time-point conventional segmentation (CSEG) by assessing the LVC relationship to brain atrophy and disease duration, as well as scan-rescan reproducibility and annual rates of lesion accrual.
MATERIALS AND METHODS: Pairs of scans were acquired 1.5 to 4.7 years apart in 21 patients with multiple sclerosis (MS). Scan-rescan MR images were acquired within 30 minutes in 10 patients with MS. LVC was measured with CSEG and SSEG after coregistration and normalization. Coefficient of variation (COV) and Bland-Altman analyses estimated method reproducibility. Spearman rank correlations probed associations between LVC and other measures.
RESULTS: Atrophy rate and net LVC were associated for SSEG (R = -0.446; P < .05) but not when using CSEG (R = -0.180; P = .421). Disease duration did not show an association with net lesion volume change per year measured by CSEG (R = -0.360; P = .11) but showed an inverse correlation with SSEG-derived measurements (R = -0.508; P < .05). Scan-rescan COV was lower for SSEG (0.98% +/- 1.55%) than for CSEG (8.64% +/- 9.91%).
CONCLUSION: SSEG unveiled a relationship between T2 LVC and concomitant brain atrophy and demonstrated significantly higher measurement reproducibility. SSEG, a promising tool providing detailed analysis of subtle alterations in lesion size and intensity, may provide critical outcome measures for clinical trials of novel treatments, and may provide further insight into progression patterns in MS.
BACKGROUND AND PURPOSE: White-matter hyperintensities (WMHs) detected by magnetic resonance imaging are thought to represent the effects of cerebral small-vessel disease and neurodegenerative changes. We sought to determine whether the spatial distribution of WMHs discriminates between different disease groups and healthy aging individuals and whether these distributions are related to local cerebral perfusion patterns.
METHODS: We examined the pattern of WMHs by T2/fluid-attenuated inversion recovery-weighted magnetic resonance imaging in 3 groups of subjects: cerebral amyloid angiopathy (n=32), Alzheimer disease or mild cognitive impairment (n=41), and healthy aging (n=29). WMH frequency maps were calculated for each group, and spatial distributions were compared by voxel-wise logistic regression. WMHs were also analyzed as a function of normal cerebral perfusion patterns by overlaying a single photon emission computed tomography atlas.
RESULTS: Although WMH volume was greater in cerebral amyloid angiopathy and Alzheimer disease/mild cognitive impairment than in healthy aging, there was no consistent difference in the spatial distributions when controlling for total WMH volume. Hyperintensities were most frequent in the deep periventricular WM in all 3 groups. A strong inverse correlation between hyperintensity frequency and normal perfusion was demonstrated in all groups, demonstrating that WMHs were most common in regions of relatively lower normal cerebral perfusion.
CONCLUSIONS: WMHs show a common distribution pattern and predilection for cerebral WM regions with lower atlas-derived perfusion, regardless of the underlying diagnosis. These data suggest that across diverse disease processes, WM injury may occur in a pattern that reflects underlying tissue properties, such as relative perfusion.
The accuracy and precision of segmentations of medical images has been difficult to quantify in the absence of a 'ground truth' or reference standard segmentation for clinical data. Although physical or digital phantoms can help by providing a reference standard, they do not allow the reproduction of the full range of imaging and anatomical characteristics observed in clinical data. An alternative assessment approach is to compare with segmentations generated by domain experts. Segmentations may be generated by raters who are trained experts or by automated image analysis algorithms. Typically, these segmentations differ due to intra-rater and inter-rater variability. The most appropriate way to compare such segmentations has been unclear. We present here a new algorithm to enable the estimation of performance characteristics, and a true labelling, from observations of segmentations of imaging data where segmentation labels may be ordered or continuous measures. This approach may be used with, among others, surface, distance transform or level-set representations of segmentations, and can be used to assess whether or not a rater consistently overestimates or underestimates the position of a boundary.
We have investigated the adsorption of asymmetric poly(styrene-b-methyl methacrylate) block copolymers (PS-PMMA) from a selective solvent onto alumina (Al(2)O(3)) particles having variable and controllable radii. The solvent used was a bad solvent for the PS block (block A) and a good solvent for the PMMA block (block B), which has a higher affinity of the surface. Such a case represents a new class of adsorption, where both blocks compete for the adsorption sites of the metallic surface. Two theoretical models, the modified drops model and the perforated film model, have been evaluated as appropriate representation of such an adsorption scenario. The experimental results indicated that the adsorption of the PS-PMMA block copolymer generated a patterned surface comprised of a homogeneous melt layer of the PS block perforated with holes having a variable PMMA structure, depending on the distance from the bottom of the hole (alumina surface) and the distance from walls of the hole. The density gradient of the PMMA moiety in the hole reverted to the classical brush morphology at a critical distance from the surface of the hole.
OBJECTIVE: Impairment of white matter connecting frontal and temporal cortices has been reported in schizophrenia. Yet, not much is known about the effects of age on fibers connecting these brain regions. Using diffusion tensor imaging tractography, we investigated the relationship between age and fiber integrity in patients with schizophrenia vs. healthy adults.
METHODS: DTI tractography was used to create 3D reconstructions of the cingulum, uncinate and inferior occipito-frontal fasciculi in 27 patients with schizophrenia and 34 healthy volunteers (23-56 years of age, group-matched on age). Fractional anisotropy (FA), describing fiber integrity, was then calculated along the entire length of these tracts, and correlated with subjects' age.
RESULTS: Patients revealed a significant decline in FA with age in both the cingulum and uncinate, but not in the inferior occipito-frontal fasciculi. No statistically significant correlations were found in these fiber bundles in controls.
CONCLUSIONS: These results suggest an age-associated reduction of frontal-temporal connectivity in schizophrenia, but not in healthy controls.
In this paper, we describe an approach to the problem of simultaneously enhancing image sequences and tracking the objects of interest represented by the latter. The enhancement part of the algorithm is based on Bayesian wavelet denoising, which has been chosen due to its exceptional ability to incorporate diverse a priori information into the process of image recovery. In particular, we demonstrate that, in dynamic settings, useful statistical priors can come both from some reasonable assumptions on the properties of the image to be enhanced as well as from the images that have already been observed before the current scene. Using such priors forms the main contribution of the present paper which is the proposal of the dynamic denoising as a tool for simultaneously enhancing and tracking image sequences. Within the proposed framework, the previous observations of a dynamic scene are employed to enhance its present observation. The mechanism that allows the fusion of the information within successive image frames is Bayesian estimation, while transferring the useful information between the images is governed by a Kalman filter that is used for both prediction and estimation of the dynamics of tracked objects. Therefore, in this methodology, the processes of target tracking and image enhancement "collaborate" in an interlacing manner, rather than being applied separately. The dynamic denoising is demonstrated on several examples of SAR imagery. The results demonstrated in this paper indicate a number of advantages of the proposed dynamic denoising over "static" approaches, in which the tracking images are enhanced independently of each other.
Segmentation involves separating an object from the background in a given image. The use of image information alone often leads to poor segmentation results due to the presence of noise, clutter or occlusion. The introduction of shape priors in the geometric active contour (GAC) framework has proved to be an effective way to ameliorate some of these problems. In this work, we propose a novel segmentation method combining image information with prior shape knowledge, using level-sets. Following the work of Leventon et al., we propose to revisit the use of PCA to introduce prior knowledge about shapes in a more robust manner. We utilize kernel PCA (KPCA) and show that this method outperforms linear PCA by allowing only those shapes that are close enough to the training data. In our segmentation framework, shape knowledge and image information are encoded into two energy functionals entirely described in terms of shapes. This consistent description permits to fully take advantage of the Kernel PCA methodology and leads to promising segmentation results. In particular, our shape-driven segmentation technique allows for the simultaneous encoding of multiple types of shapes, and offers a convincing level of robustness with respect to noise, occlusions, or smearing.
This paper proposes a deterministic observer framework for visual tracking based on non-parametric implicit (level-set) curve descriptions. The observer is continuous-discrete, with continuous-time system dynamics and discrete-time measurements. Its state-space consists of an estimated curve position augmented by additional states (e.g., velocities) associated with every point on the estimated curve. Multiple simulation models are proposed for state prediction. Measurements are performed through standard static segmentation algorithms and optical-flow computations. Special emphasis is given to the geometric formulation of the overall dynamical system. The discrete-time measurements lead to the problem of geometric curve interpolation and the discrete-time filtering of quantities propagated along with the estimated curve. Interpolation and filtering are intimately linked to the correspondence problem between curves. Correspondences are established by a Laplace-equation approach. The proposed scheme is implemented completely implicitly (by Eulerian numerical solutions of transport equations) and thus naturally allows for topological changes and subpixel accuracy on the computational grid.
We have previously developed a fast Monte Carlo (MC)-based joint ordered-subset expectation maximization (JOSEM) iterative reconstruction algorithm, MC-JOSEM. A phantom study was performed to compare quantitative imaging performance of MC-JOSEM with that of a triple-energy-window approach (TEW) in which estimated scatter was also included additively within JOSEM, TEW-JOSEM. We acquired high-count projections of a 5.5 cm3 sphere of 111In at different locations in the water-filled torso phantom; high-count projections were then obtained with 111In only in the liver or only in the soft-tissue background compartment, so that we could generate synthetic projections for spheres surrounded by various activity distributions. MC scatter estimates used by MC-JOSEM were computed once after five iterations of TEW-JOSEM. Images of different combinations of liver/background and sphere/background activity concentration ratios were reconstructed by both TEW-JOSEM and MC-JOSEM for 40 iterations. For activity estimation in the sphere, MC-JOSEM always produced better relative bias and relative standard deviation than TEW-JOSEM for each sphere location, iteration number, and activity combination. The average relative bias of activity estimates in the sphere for MC-JOSEM after 40 iterations was -6.9%, versus -15.8% for TEW-JOSEM, while the average relative standard deviation of the sphere activity estimates was 16.1% for MC-JOSEM, versus 27.4% for TEW-JOSEM. Additionally, the average relative bias of activity concentration estimates in the liver and the background for MC-JOSEM after 40 iterations was -3.9%, versus -12.2% for TEW-JOSEM, while the average relative standard deviation of these estimates was 2.5% for MC-JOSEM, versus 3.4% for TEW-JOSEM. MC-JOSEM is a promising approach for quantitative activity estimation in 111In SPECT.
The advancement of the computational biology field hinges on progress in three fundamental directions--the development of new computational algorithms, the availability of informatics resource management infrastructures and the capability of tools to interoperate and synergize. There is an explosion in algorithms and tools for computational biology, which makes it difficult for biologists to find, compare and integrate such resources. We describe a new infrastructure, iTools, for managing the query, traversal and comparison of diverse computational biology resources. Specifically, iTools stores information about three types of resources--data, software tools and web-services. The iTools design, implementation and resource meta-data content reflect the broad research, computational, applied and scientific expertise available at the seven National Centers for Biomedical Computing. iTools provides a system for classification, categorization and integration of different computational biology resources across space-and-time scales, biomedical problems, computational infrastructures and mathematical foundations. A large number of resources are already iTools-accessible to the community and this infrastructure is rapidly growing. iTools includes human and machine interfaces to its resource meta-data repository. Investigators or computer programs may utilize these interfaces to search, compare, expand, revise and mine meta-data descriptions of existent computational biology resources. We propose two ways to browse and display the iTools dynamic collection of resources. The first one is based on an ontology of computational biology resources, and the second one is derived from hyperbolic projections of manifolds or complex structures onto planar discs. iTools is an open source project both in terms of the source code development as well as its meta-data content. iTools employs a decentralized, portable, scalable and lightweight framework for long-term resource management. We demonstrate several applications of iTools as a framework for integrated bioinformatics. iTools and the complete details about its specifications, usage and interfaces are available at the iTools web page http://iTools.ccb.ucla.edu.
Many promising MRI approaches for research or clinical management of multiple sclerosis (MS) have recently emerged, or are under development or refinement. Advanced MRI methods need to be assessed to determine whether they allow earlier diagnosis or better identification of phenotypes. Improved post-processing should allow more efficient and complete extraction of information from images. Magnetic resonance spectroscopy should improve in sensitivity and specificity with higher field strengths and should enable the detection of a wider array of metabolites. Diffusion imaging is moving closer to the goal of defining structural connectivity and, thereby, determining the functional significance of lesions at specific locations. Cell-specific imaging now seems feasible with new magnetic resonance contrast agents. The imaging of myelin water fraction brings the hope of providing a specific measure of myelin content. Ultra-high-field MRI increases sensitivity, but also presents new technical challenges. Here, we review these recent developments in MRI for MS, and also look forward to refinements in spinal-cord imaging, optic-nerve imaging, perfusion MRI, and functional MRI. Advances in MRI should improve our ability to diagnose, monitor, and understand the pathophysiology of MS.
BACKGROUND: We have acquired dual-echo spin-echo (DE SE) MRI data of the rhesus monkey brain since 1994 as part of an ongoing study of normal aging. To analyze these legacy data for regional volume changes, we have created a reference label atlas for the Template Driven Segmentation (TDS) algorithm.
METHODS: The atlas was manually created from DE SE legacy MRI data of one behaviorally normal, young, male rhesus monkey and consisted of 14 regions of interest (ROI's). We analyzed the reproducibility and validity of the TDS algorithm using the atlas relative to manual segmentation.
RESULTS: ROI volumes were comparable between the two segmentation methodologies, except TDS overestimated the volume of basal ganglia regions. Both methodologies were highly reproducible, but TDS had lower sensitivity and comparable specificity.
CONCLUSIONS: TDS segmentation calculates accurate volumes for most ROI's. Sensitivity will be improved in future studies through the acquisition of higher quality data.