Raj Shekhar

Mutual information-based rigid and nonrigid registration of ultrasound volumes

We investigated the registration of ultrasound volumes based on the mutual information measure, a technique originally applied to multimodality registration of brain images. A prerequisite for successful registration is a smooth, quasi-convex mutual information surface with an unambiguous maximum. We discuss the necessary preprocessing to create such a surface for ultrasound volumes. Abdominal and thoracic organs imaged with ultrasound typically move relative to the exterior of the body and are deformable. Consequently, four specific instances of image registration involving progressively generalized transformations were studied: rigid-body, rigid-body + uniform scaling, rigid-body + nonuniform scaling, and affine. Registration was applied to clinically acquired volumetric images. The accuracy was comparable with the voxel dimension for all transformation modes, although it degraded as the transformation grew more complex. Likewise, the capture range became narrower with the complexity of transformation. As the use of real-time three-dimensional ultrasound becomes more prevalent, the method we present should work well for a variety of applications examining serial anatomic and physiologic changes. Developers of these clinical applications would match the deformation model of their problem to one of the four transformation models presented here.

Evaluation of three-dimensional segmentation algorithms for the identification of luminal and medial-adventitial borders in intravascular ultrasound images

Intravascular ultrasound (IVUS) provides direct depiction of coronary artery anatomy, including plaque and vessel area, which is important in quantitative studies on the progression or regression of coronary artery disease. Traditionally, these studies have relied on manual evaluation, which is laborious, time consuming, and subject to large interobserver and intraobserver variability. A new technique, called active surface segmentation, alleviates these limitations and makes strides toward routine analyses. However, for three-dimensional (3-D) plaque assessment or 3-D reconstruction to become a clinical reality, methods must be developed which can analyze many images quickly. Presented is a comparison between two active surface techniques for three-dimensional segmentation of luminal and medial-adventitial borders. The force-acceleration technique and the neighborhood-search technique accurately detected both borders in vivo (r/sup 2/=0.95 and 0.99, Williams index=0.67 and 0.65, and r/sup 2/=0.95 and 0.99, WI=0.67 and 0.70, respectively). However, the neighborhood-search technique was significantly faster and required less computation. Volume calculations for both techniques (r/sup 2/=0.99 and r/sup 2/=0.99) also agreed with a known-volume phantom. Active surface segmentation allows 3-D assessment of coronary morphology and further developments with this technology will provide clinical analysis tools.

Mathematical model of corneal surface smoothing after laser refractive surgery

PURPOSEnTo construct a quantitative model of corneal surface smoothing after laser ablation for refractive correction.nnnDESIGNnExperimental study, interventional case series, and meta-analysis of literature.nnnMETHODSnA theory of epithelial smoothing in response to corneal contour change is derived from differential equations that describe epithelial migration, growth, and loss. Computer simulations calculate the effects on postoperative epithelial thickness, topography, refraction, and spherical aberration. Model parameter is matched with laser in situ keratomileusis (LASIK) outcome in literature and in a retrospective study of primary spherical myopic (77 eyes) and hyperopic (19 eyes) corrections. Surgically induced refractive change was the main outcome measure.nnnRESULTSnSimulated epithelial remodeling after myopic ablation produces central epithelial thickening, reduction in achieved correction, and induction of oblate spherical aberration. Simulation of hyperopic ablation shows peripheral epithelial thickening, a larger reduction in correction, and induction of prolate spherical aberration. Simulation using a minus cylinder laser ablation pattern shows decreased astigmatism correction and increased hyperopic shift. In the LASIK series, linear regression of achieved correction vs ablation setting in hyperopic and minus cylinder corrections shows slopes of 0.97, 0.71, and 0.74, respectively. These clinical results match model predictions when the smoothing constant is set at 0.32, 0.63, and 0.55 mm, respectively.nnnCONCLUSIONSnEpithelial thickness modulations after ablation can be modeled mathematically to explain clinically observed regression and induction of aberration. The cornea appears to smooth over ablated features smaller than approximately 0.5 mm. The model provides an approach for designing ablation patterns that precompensate for the smoothing to improve final outcome.

Texture analysis of lesions in breast ultrasound images

We investigate the use of Haralicks texture features and posterior acoustic attenuation descriptors (PAAD) for the characterization of ultrasound (US) breast lesions. 71 lesions (24 cyst, 21 benign solid mass and 26 malignant solid masses) were manually segmented on two-dimensional breast US images. 28 Haralicks descriptors and two PAAD were evaluated on these segmented lesions. Mean of Sum Average, Range of Sum Entropy and the second PAAD best discriminated cysts from noncysts. Range of Correlation and the second PAAD best discriminated solid malignant from benign lesions. Computerized analysis of breast US images can increase the specificity of breast sonography by providing a better characterization of solid lesions.

Three-dimensional segmentation of luminal and adventitial borders in serial intravascular ultrasound images

Intravascular ultrasound (IVUS) provides exact anatomy of arteries, allowing accurate quantitative analysis. Automated segmentation of IVUS images is a prerequisite for routine quantitative analyses. We present a new three-dimensional (3D) segmentation technique, called active surface segmentation, which detects luminal and adventitial borders in IVUS pullback examinations of coronary arteries. The technique was validated against expert tracings by computing correlation coefficients (range 0.83-0.97) and Williams index values (range 0.37-0.66). The technique was statistically accurate, robust to image artifacts, and capable of segmenting a large number of images rapidly. Active surface segmentation enabled geometrically accurate 3D reconstruction and visualization of coronary arteries and volumetric measurements.

Registration of real-time 3-D ultrasound images of the heart for novel 3-D stress echocardiography

Stress echocardiography is a routinely used clinical procedure to diagnose cardiac dysfunction by comparing wall motion information in prestress and poststress ultrasound images. Incomplete data, complicated imaging protocols and misaligned prestress and poststress views, however, are known limitations of conventional stress echocardiography. We discuss how the first two limitations are overcome via the use of real-time three-dimensional (3-D) ultrasound imaging, an emerging modality, and have called the new procedure 3-D stress echocardiography: We also show that the problem of misaligned views can be solved by registration of prestress and poststress 3-D image sequences. Such images are misaligned because of variations in placing the ultrasound transducer and stress-induced anatomical changes. We have developed a technique to temporally align 3-D images of the two sequences first and then to spatially register them to rectify probe placement error while preserving the stress-induced changes. The 3-D spatial registration is mutual information-based. Image registration used in conjunction with 3-D stress echocardiography can potentially improve the diagnostic accuracy of stress testing.

Three-dimensional reconstruction of the coronary artery wall by image fusion of intravascular ultrasound and bi-plane angiography

Background: Intravascular ultrasound (IVUS) is becoming increasingly accepted for assessing coronary anatomy. However, its utility in visualizing and quantifying coronary morphology has been limited by its 2D tomographic nature. This study presents a 3D reconstruction technique that accurately preserves 3D geometric information. Methods and Results: Images obtained from manual IVUS pullbacks and continuous bi-plane angiography were fused, using angiography to reconstruct the transducer trajectory and aid in solving for the correct rotational orientation. A novel 3D active surface method automatically identified the luminal and medial–adventitial borders which, when superimposed on the transducer trajectory, could be surface-rendered for visualization and morphometry. Segmentation agreed well with manual assessment, and 3D luminal shape matched that of angiography when projected to 2D. Conclusions: We conclude that this method provides an accurate reconstruction of the vessels anatomy, which accounts for the true curvature of the vessel.

Assessment of coronary compensatory enlargement by three-dimensional intravascular ultrasound.

Several techniques have been used to demonstrate that human arteries respond to atherosclerosis by increasing their total arterial area to prevent a decrease in blood flow. Three-dimensional reconstructions of coronary arteries can document this compensatory response accurately and specifically. Seven human coronary arteries were reconstructed using intravascular ultrasound and biplane angiography, and vessel geometries were quantified. In all seven vessels, as plaque area increased, overall vessel area increased (R = 0.986, 0.933, 0.984, 0.678, 0.763, 0.963, and 0.830), but luminal cross-sectional area did not significantly decrease. Focal compensatory enlargement was identified in each vessel, and in some cases this response appeared to occur until the vessel was 65% occluded. Luminal enlargement near the proximal ends was attributed to the natural taper of the vessel. The semi-automated, three-dimensional segmentation technique used in this study allows reproducible quantification, as there is no subjective manual tracing involved. Following the intravascular ultrasound transducer in time and space with biplane angiography allows for accurate reconstruction with or without automated pullback devices. Information on the rate of change of vessel measurements is also presented, which, when combined with visualization of accurate 3D geometry, provides a unique assessment of coronary compensatory enlargement. This reconstruction technique can be applied in a clinical environment with no major modification.

Mutual information-based multimodality registration of cardiac ultrasound and SPECT images: a preliminary investigation.

Background: Ultrasound (US) and single photon emission computed tomography (SPECT) are the two most commonly prescribed procedures for diagnosing coronary artery disease (CAD). We have demonstrated the feasibility of multimodality registration of two-dimensional (2D) and three-dimensional (3D) cardiac US images with cardiac SPECT images with an aim to simultaneously present the complementary anatomical and perfusion information from the two modalities. We have also tested the clinicians assessment of the clinical adequacy of the registered images. Methods and Results: We have demonstrated temporal and spatial registration for nine sets of cardiac US and SPECT cine loops covering the entire cardiac cycle. Temporal alignment was performed by interpolation of existing SPECT images at cardiac phases corresponding to available US images. Spatial registration was performed in 3D image space using a mutual information-based approach. Experts from echocardiography and nuclear medicine determined the clinical utility of the registration by rating each registration on a scale of 1 to 5, a rating of 3 or above indicating clinical utility. 2DUS-SPECT registration (five cases) received an average rating of 4.2, whereas 3DUS-SPECT registration (four cases) received an average rating of 2.85. By one-sample t test, the overall evaluations (mean 3.58) were greater than the pre-specified clinical cut-off of 3 with p < 0.05, indicating likelihood of clinical utility. Conclusion: Our method demonstrates the feasibility of registering cardiac US and SPECT images in their present as well as possible future forms. Such registration has the potential to provide a more accurate and powerful tool for diagnosing CAD.

Cine MPR: interactive multiplanar reformatting of four-dimensional cardiac data using hardware-accelerated texture mapping

Four-dimensional (4-D) imaging to capture the three-dimensional (3-D) structure and motion of the heart in real time is an emerging trend. We present here our method of interactive multiplanar reformatting (MPR), i.e., the ability to visualize any chosen anatomical cross section of 4-D cardiac images and to change its orientation smoothly while maintaining the original heart motion. Continuous animation to show the time-varying 3-D geometry of the heart and smooth dynamic manipulation of the reformatted planes, as well as large image size (100-300 MB), make MPR challenging. Our solution exploits the hardware acceleration of 3-D texture mapping capability of high-end commercial PC graphics boards. Customization of volume subdivision and caching concepts to periodic cardiac data allows us to use this hardware effectively and efficiently. We are able to visualize and smoothly interact with real-time 3-D ultrasound cardiac images at the desired frame rate (25 Hz). The developed methods are applicable to MPR of one or more 3-D and 4-D medical images, including 4-D cardiac images collected in a gated fashion.