Improved myocardial tagging contrast in cine balanced SSFP images.

PURPOSE: To improve the tag persistence throughout the whole cardiac cycle by providing a constant tag-contrast throughout all the cardiac phases when using balanced steady-state free precession (bSSFP) imaging. MATERIALS AND METHODS: The flip angles of the imaging radiofrequency pulses were optimized to compensate for the tagging contrast-to-noise ratio (Tag-CNR) fading at later cardiac phases in bSSFP imaging. Complementary spatial modulation of magnetization (CSPAMM) tagging was implemented to improve the Tag-CNR. Numerical simulations were performed to examine the behavior of the Tag-CNR with the proposed method, and to compare the resulting Tag-CNR with that obtained from the more commonly used spoiled gradient echo (SPGR) imaging. A gel phantom, as well as five healthy human volunteers, were scanned on a 1.5T scanner using bSSFP imaging with and without the proposed technique. The phantom was also scanned with SPGR imaging. RESULTS: With the proposed technique, the Tag-CNR remained almost constant during the whole cardiac cycle. Using bSSFP imaging, the Tag-CNR was about double that of SPGR. CONCLUSION: The tag persistence was significantly improved when the proposed method was applied, with better Tag-CNR during the diastolic cardiac phase. The improved Tag-CNR will support automated tagging analysis and quantification methods.

Direct in vitro comparison of six three-dimensional positive contrast methods for susceptibility marker imaging.

PURPOSE: To compare different techniques for positive contrast imaging of susceptibility markers with MRI for three-dimensional visualization. As several different techniques have been reported, the choice of the suitable method depends on its properties with regard to the amount of positive contrast and the desired background suppression, as well as other imaging constraints needed for a specific application. MATERIALS AND METHODS: Six different positive contrast techniques are investigated for their ability to image at 3 Tesla a single susceptibility marker in vitro. The white marker method (WM), susceptibility gradient mapping (SGM), inversion recovery with on-resonant water suppression (IRON), frequency selective excitation (FSX), fast low flip-angle positive contrast SSFP (FLAPS), and iterative decomposition of water and fat with echo asymmetry and least-squares estimation (IDEAL) were implemented and investigated. RESULTS: The different methods were compared with respect to the volume of positive contrast, the product of volume and signal intensity, imaging time, and the level of background suppression. Quantitative results are provided, and strengths and weaknesses of the different approaches are discussed. CONCLUSION: The appropriate choice of positive contrast imaging technique depends on the desired level of background suppression, acquisition speed, and robustness against artifacts, for which in vitro comparative data are now available.

Validation of two echocardiographic indexes to improve the diagnosis of complex coarctations.

OBJECTIVES: Coarctation of the aorta is one of the most common congenital heart defects. Its diagnosis may be difficult in the presence of a patent ductus arteriosus, of other complex defects or of a poor echocardiographic window. We sought to demonstrate that the carotid-subclavian artery index (CSA index) and the isthmus-descending aorta ratio (I/D ratio), two recently described echocardiographic indexes, are effective in detection of isolated and complex aortic coarctations in children younger and older than 3 months of age. The CSA index is the ratio of the distal aortic arch diameter to the distance between the left carotid artery and the left subclavian artery. It is highly suggestive of a coarctation when it is <1.5. The I/D ratio defined as the diameter of the isthmus to the diameter of the descending aorta, suggests an aortic coarctation when it is less than 0.64. METHODS: This is a retrospective cohort study in a tertiary care children's hospital. Review of all echocardiograms in children aged 0-18 years with a diagnosis of coarctation seen at the author's institution between 1996 and 2006. An age- and sex-matched control group without coarctation was constituted. Offline echocardiographic measurements of the aortic arch were performed in order to calculate the CSA index and I/D ratio. RESULTS: Sixty-eight patients were included in the coarctation group, 24 in the control group. Patients with coarctation had a significantly lower CSA index (0.84+/-0.39 vs 2.65+/-0.82, p<0.0001) and I/D ratio (0.58+/-0.18 vs 0.98+/-0.19, p<0.0001) than patients in the control group. Associated cardiac defects and age of the child did not significantly alter the CSA index or the I/D ratio. CONCLUSIONS: A CSA index less than 1.5 is highly suggestive of coarctation independent of age and of the presence of other cardiac defects. I/D ratio alone is less specific than CSA alone at any age and for any associated cardiac lesion. The association of both indexes improves sensitivity and permits diagnosis of coarctation in all patients based solely on a bedside echocardiographic measurement.

An active contour-based atlas registration model applied to automatic subthalamic nucleus targeting on MRI: method and validation.

This paper presents a new non parametric atlas registration framework, derived from the optical flow model and the active contour theory, applied to automatic subthalamic nucleus (STN) targeting in deep brain stimulation (DBS) surgery. In a previous work, we demonstrated that the STN position can be predicted based on the position of surrounding visible structures, namely the lateral and third ventricles. A STN targeting process can thus be obtained by registering these structures of interest between a brain atlas and the patient image. Here we aim to improve the results of the state of the art targeting methods and at the same time to reduce the computational time. Our simultaneous segmentation and registration model shows mean STN localization errors statistically similar to the most performing registration algorithms tested so far and to the targeting expert's variability. Moreover, the computational time of our registration method is much lower, which is a worthwhile improvement from a clinical point of view.

Active Contour-Based Segmentation of Head and Neck with Adaptive Atlas Selection

This paper presents automated segmentation of structuresin the Head and Neck (H\&N) region, using an activecontour-based joint registration and segmentation model.A new atlas selection strategy is also used. Segmentationis performed based on the dense deformation fieldcomputed from the registration of selected structures inthe atlas image that have distinct boundaries, onto thepatient's image. This approach results in robustsegmentation of the structures of interest, even in thepresence of tumors, or anatomical differences between theatlas and the patient image. For each patient, an atlasimage is selected from the available atlas-database,based on the similarity metric value, computed afterperforming an affine registration between each image inthe atlas-database and the patient's image. Unlike manyof the previous approaches in the literature, thesimilarity metric is not computed over the entire imageregion; rather, it is computed only in the regions ofsoft tissue structures to be segmented. Qualitative andquantitative evaluation of the results is presented.

Fast Geodesic Active Fields for Image Registration Based on Splitting and Augmented Lagrangian Approaches.

In this paper, we present an efficient numerical scheme for the recently introduced geodesic active fields (GAF) framework for geometric image registration. This framework considers the registration task as a weighted minimal surface problem. Hence, the data-term and the regularization-term are combined through multiplication in a single, parametrization invariant and geometric cost functional. The multiplicative coupling provides an intrinsic, spatially varying and data-dependent tuning of the regularization strength, and the parametrization invariance allows working with images of nonflat geometry, generally defined on any smoothly parametrizable manifold. The resulting energy-minimizing flow, however, has poor numerical properties. Here, we provide an efficient numerical scheme that uses a splitting approach; data and regularity terms are optimized over two distinct deformation fields that are constrained to be equal via an augmented Lagrangian approach. Our approach is more flexible than standard Gaussian regularization, since one can interpolate freely between isotropic Gaussian and anisotropic TV-like smoothing. In this paper, we compare the geodesic active fields method with the popular Demons method and three more recent state-of-the-art algorithms: NL-optical flow, MRF image registration, and landmark-enhanced large displacement optical flow. Thus, we can show the advantages of the proposed FastGAF method. It compares favorably against Demons, both in terms of registration speed and quality. Over the range of example applications, it also consistently produces results not far from more dedicated state-of-the-art methods, illustrating the flexibility of the proposed framework.

High-resolution three-dimensional aortic magnetic resonance angiography and quantitative vessel wall characterization of different atherosclerotic stages in a rabbit model

PURPOSE: Atherosclerosis results in a considerable medical and socioeconomic impact on society. We sought to evaluate novel magnetic resonance imaging (MRI) angiography and vessel wall sequences to visualize and quantify different morphologic stages of atherosclerosis in a Watanabe hereditary hyperlipidemic (WHHL) rabbit model. MATERIAL AND METHODS: Aortic 3D steady-state free precession angiography and subrenal aortic 3D black-blood fast spin-echo vessel wall imaging pre- and post-Gadolinium (Gd) was performed in 14 WHHL rabbits (3 normal, 6 high-cholesterol diet, and 5 high-cholesterol diet plus endothelial denudation) on a commercial 1.5 T MR system. Angiographic lumen diameter, vessel wall thickness, signal-/contrast-to-noise analysis, total vessel area, lumen area, and vessel wall area were analyzed semiautomatically. RESULTS: Pre-Gd, both lumen and wall dimensions (total vessel area, lumen area, vessel wall area) of group 2 + 3 were significantly increased when compared with those of group 1 (all P < 0.01). Group 3 animals had significantly thicker vessel walls than groups 1 and 2 (P < 0.01), whereas angiographic lumen diameter was comparable among all groups. Post-Gd, only diseased animals of groups 2 + 3 showed a significant (>100%) signal-to-noise ratio and contrast-to-noise increase. CONCLUSIONS: A combination of novel 3D magnetic resonance angiography and high-resolution 3D vessel wall MRI enabled quantitative characterization of various atherosclerotic stages including positive arterial remodeling and Gd uptake in a WHHL rabbit model using a commercially available 1.5 T MRI system.

Model-based Segmentation and Image Fusion of 3D Computed Tomography and 3D Ultrasound of the Eye for Radiotherapy Planning

For radiotherapy treatment planning of retinoblastoma inchildhood, Computed Tomography (CT) represents thestandard method for tumor volume delineation, despitesome inherent limitations. CT scan is very useful inproviding information on physical density for dosecalculation and morphological volumetric information butpresents a low sensitivity in assessing the tumorviability. On the other hand, 3D ultrasound (US) allows ahigh accurate definition of the tumor volume thanks toits high spatial resolution but it is not currentlyintegrated in the treatment planning but used only fordiagnosis and follow-up. Our ultimate goal is anautomatic segmentation of gross tumor volume (GTV) in the3D US, the segmentation of the organs at risk (OAR) inthe CT and the registration of both. In this paper, wepresent some preliminary results in this direction. Wepresent 3D active contour-based segmentation of the eyeball and the lens in CT images; the presented approachincorporates the prior knowledge of the anatomy by usinga 3D geometrical eye model. The automated segmentationresults are validated by comparing with manualsegmentations. Then, for the fusion of 3D CT and USimages, we present two approaches: (i) landmark-basedtransformation, and (ii) object-based transformation thatmakes use of eye ball contour information on CT and USimages.

Direct three-dimensional myocardial strain tensor quantification and tracking using zHARP.

Images of myocardial strain can be used to diagnose heart disease, plan and monitor treatment, and to learn about cardiac structure and function. Three-dimensional (3D) strain is typically quantified using many magnetic resonance (MR) images obtained in two or three orthogonal planes. Problems with this approach include long scan times, image misregistration, and through-plane motion. This article presents a novel method for calculating cardiac 3D strain using a stack of two or more images acquired in only one orientation. The zHARP pulse sequence encodes in-plane motion using MR tagging and out-of-plane motion using phase encoding, and has been previously shown to be capable of computing 3D displacement within a single image plane. Here, data from two adjacent image planes are combined to yield a 3D strain tensor at each pixel; stacks of zHARP images can be used to derive stacked arrays of 3D strain tensors without imaging multiple orientations and without numerical interpolation. The performance and accuracy of the method is demonstrated in vitro on a phantom and in vivo in four healthy adult human subjects.

Efficient total variation algorithm for fetal brain MRI reconstruction.

Fetal MRI reconstruction aims at finding a high-resolution image given a small set of low-resolution images. It is usually modeled as an inverse problem where the regularization term plays a central role in the reconstruction quality. Literature has considered several regularization terms s.a. Dirichlet/Laplacian energy, Total Variation (TV)- based energies and more recently non-local means. Although TV energies are quite attractive because of their ability in edge preservation, standard explicit steepest gradient techniques have been applied to optimize fetal-based TV energies. The main contribution of this work lies in the introduction of a well-posed TV algorithm from the point of view of convex optimization. Specifically, our proposed TV optimization algorithm or fetal reconstruction is optimal w.r.t. the asymptotic and iterative convergence speeds O(1/n2) and O(1/√ε), while existing techniques are in O(1/n2) and O(1/√ε). We apply our algorithm to (1) clinical newborn data, considered as ground truth, and (2) clinical fetal acquisitions. Our algorithm compares favorably with the literature in terms of speed and accuracy.

Functional compensation of motor function in pre-symptomatic Huntington's disease.

Involuntary choreiform movements are a clinical hallmark of Huntington's disease. Studies in clinically affected patients suggest a shift of motor activations to parietal cortices in response to progressive neurodegeneration. Here, we studied pre-symptomatic gene carriers to examine the compensatory mechanisms that underlie the phenomenon of retained motor function in the presence of degenerative change. Fifteen pre-symptomatic gene carriers and 12 matched controls performed button presses paced by a metronome at either 0.5 or 2 Hz with four fingers of the right hand whilst being scanned with functional magnetic resonance imaging. Subjects pressed buttons either in the order of a previously learnt 10-item finger sequence, from left to right, or kept still. Error rates ranged from 2% to 7% in the pre-symptomatic gene carriers and from 0.5% to 4% in controls, depending on the condition. No significant difference in task performance was found between groups for any of the conditions. Activations in the supplementary motor area (SMA) and superior parietal lobe differed with gene status. Compared with healthy controls, gene carriers showed greater activations of left caudal SMA with all movement conditions. Activations correlated with increasing speed of movement were greater the closer the gene carriers were to estimated clinical diagnosis, defined by the onset of unequivocal motor signs. Activations associated with increased movement complexity (i.e. with the pre-learnt 10-item sequence) decreased in the rostral SMA with nearing diagnostic onset. The left superior parietal lobe showed reduced activation with increased movement complexity in gene carriers compared with controls, and in the right superior parietal lobe showed greater activations with all but the most demanding movements. We identified a complex pattern of motor compensation in pre-symptomatic gene carriers. The results show that preclinical compensation goes beyond a simple shift of activity from premotor to parietal regions involving multiple compensatory mechanisms in executive and cognitive motor areas. Critically, the pattern of motor compensation is flexible depending on the actual task demands on motor control.

Use of a dynamic foot pressure index to monitor the effects of treatment for equinus gait in children with cerebral palsy.

The purpose of this study is to introduce and describe a newly developed index using foot pressure analysis to quantify the degree of equinus gait in children with cerebral palsy before and after injection with botulinum toxin. Data were captured preinjection and 12 weeks postinjection. Ten children aged 2(1/2) to 6(1/2) years took part (5 boys and 5 girls). Three of them had a diagnosis of spastic diplegia and 7 of congenital hemiplegia. In total, 13 limbs were analyzed. After orientation and segmentation of raw pedobarographic data, we determined a dynamic foot pressure index graded 0 to 100 that quantified the relative degree of heel and forefoot contact during stance. These data were correlated (Pearson correlation) with clinical measurements of dorsiflexion at the ankle (on a slow and fast stretch) and video observation (using the Observational Gait Scale). Pedobarograph data were strongly correlated with both the Observational Gait Scale scores (R = 0.79, P < 0.005) and clinical measurements of dorsiflexion on a fast stretch, which is reflective of spasticity (R = 0.70, P < 0.005). We demonstrated the index's sensitivity in detecting changes in spasticity and good correlation with video observations seems to indicate this technique's potential validity. When manipulated and segmented appropriately, and with the development of a simple ordinal index, we found that foot pressure data provided a useful tool in tracking changes in patients with spastic equinus.

Reconfigurable Architecture To Estimate The Motion Of An Underwater Vehicle

This work proposes a parallel architecture for a motion estimation algorithm. It is well known that image processing requires a huge amount of computation, mainly at low level processing where the algorithms are dealing with a great numbers of data-pixel. One of the solutions to estimate motions involves detection of the correspondences between two images. Due to its regular processing scheme, parallel implementation of correspondence problem can be an adequate approach to reduce the computation time. This work introduces parallel and real-time implementation of such low-level tasks to be carried out from the moment that the current image is acquired by the camera until the pairs of point-matchings are detected

Nontraumatic subarachnoid hemorrhage management: evaluation with reduced iodine volume at CT angiography.

PURPOSE: To evaluate the technical quality and the diagnostic performance of a protocol with use of low volumes of contrast medium (25 mL) at 64-detector spiral computed tomography (CT) in the diagnosis and management of adult, nontraumatic subarachnoid hemorrhage (SAH). MATERIALS AND METHODS: This study was performed outside the United States and was approved by the institutional review board. Intracranial CT angiography was performed in 73 consecutive patients with nontraumatic SAH diagnosed at nonenhanced CT. Image quality was evaluated by two observers using two criteria: degree of arterial enhancement and venous contamination. The two independent readers evaluated diagnostic performance (lesion detection and correct therapeutic decision-making process) by using rotational angiographic findings as the standard of reference. Sensitivity, specificity, and positive and negative predictive values were calculated for patients who underwent CT angiography and three-dimensional rotational angiography. The intraclass correlation coefficient was calculated to assess interobserver concordance concerning aneurysm measurements and therapeutic management. RESULTS: All aneurysms were detected, either ruptured or unruptured. Arterial opacification was excellent in 62 cases (85%), and venous contamination was absent or minor in 61 cases (84%). In 95% of cases, CT angiographic findings allowed optimal therapeutic management. The intraclass correlation coefficient ranged between 0.93 and 0.95, indicating excellent interobserver agreement. CONCLUSION: With only 25 mL of iodinated contrast medium focused on the arterial phase, 64-detector CT angiography allowed satisfactory diagnostic and therapeutic management of nontraumatic SAH.