Analysis of Left Atrial Respiratory and Cardiac Motion for Cardiac Ablation Therapy.

Cardiac ablation therapy is often guided by models built from preoperative computed tomography (CT) or magnetic resonance imaging (MRI) scans. One of the challenges in guiding a procedure from a preoperative model is properly synching the preoperative models with cardiac and respiratory motion through computational motion models. In this paper, we describe a methodology for evaluating cardiac and respiratory motion in the left atrium and pulmonary veins of a beating canine heart. Cardiac catheters were used to place metal clips within and near the pulmonary veins and left atrial appendage under fluoroscopic and ultrasound guidance and a contrast-enhanced, 64-slice multidetector CT scan was collected with the clips in place. Each clip was segmented from the CT scan at each of the five phases of the cardiac cycle at both end-inspiration and end-expiration. The centroid of each segmented clip was computed and used to evaluate both cardiac and respiratory motion of the left atrium. A total of three canine studies were completed, with 4 clips analyzed in the first study, 5 clips in the second study, and 2 clips in the third study. Mean respiratory displacement was 0.2±1.8 mm in the medial/lateral direction, 4.7±4.4 mm in the anterior/posterior direction (moving anterior on inspiration), and 9.0±5.0 mm superior/inferior (moving inferior with inspiration). At end inspiration, the mean left atrial cardiac motion at the clip locations was 1.5±1.3 mm in the medial/lateral direction, and 2.1±2.0 mm in the anterior/posterior and 1.3±1.2 mm superior/inferior directions. At end expiration, the mean left atrial cardiac motion at the clip locations was 2.0±1.5 mm in the medial/lateral direction, 3.0±1.8 mm in the anterior/posterior direction, and 1.5±1.5 mm in the superior/inferior directions.

Toward Standardized Mapping for Left Atrial Analysis and Cardiac Ablation Guidance.

In catheter-based cardiac ablation, the pulmonary vein ostia are important landmarks for guiding the ablation procedure, and for this reason, have been the focus of many studies quantifying their size, structure, and variability. Analysis of pulmonary vein structure, however, has been limited by the lack of a standardized reference space for population based studies. Standardized maps are important tools for characterizing anatomic variability across subjects with the goal of separating normal inter-subject variability from abnormal variability associated with disease. In this work, we describe a novel technique for computing flat maps of left atrial anatomy in a standardized space. A flat map of left atrial anatomy is created by casting a single ray through the volume and systematically rotating the camera viewpoint to obtain the entire field of view. The technique is validated by assessing preservation of relative surface areas and distances between the original 3D geometry and the flat map geometry. The proposed methodology is demonstrated on 10 subjects which are subsequently combined to form a probabilistic map of anatomic location for each of the pulmonary vein ostia and the boundary of the left atrial appendage. The probabilistic map demonstrates that the location of the inferior ostia have higher variability than the superior ostia and the variability of the left atrial appendage is similar to the superior pulmonary veins. This technique could also have potential application in mapping electrophysiology data, radio-frequency ablation burns, or treatment planning in cardiac ablation therapy.

Augmented environments for minimally invasive therapy: implementation barriers from technology to practice.

Augmented environments for medical applications have been explored and developed in an effort to enhance the clinician's view of anatomy and facilitate the performance of minimally invasive procedures. These environments must faithfully represent the real surgical field and require seamless integration of pre- and intra-operative imaging, surgical instrument tracking and display technology into a common framework centered around the patient. However, few image guidance environments have been successfully translated into clinical use. Several challenges that contribute to the slow progress of integrating such environments into clinical practice are discussed here in terms of both technical and clinical limitations.

Relationship of femoral neck areal bone mineral density to volumetric bone mineral density, bone size, and femoral strength in men and women.

UNLABELLED: Using combined dual-energy X-ray absorptiometry (DXA) and quantitative computed tomography, we demonstrate that men matched with women for femoral neck (FN) areal bone mineral density (aBMD) have lower volumetric BMD (vBMD), higher bone cross-sectional area, and relatively similar values for finite element (FE)-derived bone strength. INTRODUCTION: aBMD by DXA is widely used to identify patients at risk for osteoporotic fractures. aBMD is influenced by bone size (i.e., matched for vBMD, larger bones have higher aBMD), and increasing evidence indicates that absolute aBMD predicts a similar risk of fracture in men and women. Thus, we sought to define the relationships between FN aBMD (assessed by DXA) and vBMD, bone size, and FE-derived femoral strength obtained from quantitative computed tomography scans in men versus women. METHODS: We studied men and women aged 40 to 90 years and not on osteoporosis medications. RESULTS: In 114 men and 114 women matched for FN aBMD, FN total cross-sectional area was 38% higher (P < 0.0001) and vBMD was 16% lower (P < 0.0001) in the men. FE models constructed in a subset of 28 women and 28 men matched for FN aBMD showed relatively similar values for bone strength and the load-to-strength ratio in the two groups. CONCLUSIONS: In this cohort of young and old men and women from Rochester, MN, USA who are matched by FN aBMD, because of the offsetting effects of bone size and vBMD, femoral strength and the load-to-strength ratio tended to be relatively similar across the sexes.

An event-driven distributed processing architecture for image-guided cardiac ablation therapy.

Medical imaging data is becoming increasing valuable in interventional medicine, not only for preoperative planning, but also for real-time guidance during clinical procedures. Three key components necessary for image-guided intervention are real-time tracking of the surgical instrument, aligning the real-world patient space with image-space, and creating a meaningful display that integrates the tracked instrument and patient data. Issues to consider when developing image-guided intervention systems include the communication scheme, the ability to distribute CPU intensive tasks, and flexibility to allow for new technologies. In this work, we have designed a communication architecture for use in image-guided catheter ablation therapy. Communication between the system components is through a database which contains an event queue and auxiliary data tables. The communication scheme is unique in that each system component is responsible for querying and responding to relevant events from the centralized database queue. An advantage of the architecture is the flexibility to add new system components without affecting existing software code. In addition, the architecture is intrinsically distributed, in that components can run on different CPU boxes, and even different operating systems. We refer to this Framework for Image-Guided Navigation using a Distributed Event-Driven Database in Real-Time as the FINDER architecture. This architecture has been implemented for the specific application of image-guided cardiac ablation therapy. We describe our prototype image-guidance system and demonstrate its functionality by emulating a cardiac ablation procedure with a patient-specific phantom. The proposed architecture, designed to be modular, flexible, and intuitive, is a key step towards our goal of developing a complete system for visualization and targeting in image-guided cardiac ablation procedures.

Application of magnetic resonance imaging to measure fasting and postprandial volumes in humans.

Our aims were to measure the gastric volume response in excess of ingested meal volume (i.e. gastric accommodation), contribution of swallowed air to this excess, day-to-day variability of gastric volumes measured by MRI and their relationship to volumes measured by single-photon-emission computed tomography (SPECT). In 20 healthy volunteers, fasting and postprandial gastric volumes were measured after technetium(99m)-pertechnetate labeling of the gastric mucosa by SPECT and separately by MRI, using 3D gradient echo and 2D half-Fourier acquisition single-shot turbo spin echo (HASTE) sequences. Ten of these subjects had a second MRI exam to assess intra-individual variation. Thereafter, another 10 subjects had two MRI studies during which they ingested the nutrient in 30 or 150 mL aliquots. During MRI, the postprandial gastric volume change exceeded the ingested meal volume by 106 +/- 12 mL (Mean +/- SEM). The HASTE and gradient echo sequences distinguished air from fluid under fasting and postprandial conditions respectively. This postprandial excess mainly comprised air (61 +/- 5 mL), which was not significantly different when ingested as 30 or 150 mL aliquots. Fasting and postprandial gastric volumes measured by MRI were generally reproducible within subjects. During SPECT, postprandial volumes increased by 158 +/- 18 mL; gastric volumes measured by SPECT were higher than MRI. MRI measures gastric volumes with acceptable performance characteristics; the postprandial excess primarily consists of air, which is not affected by the mode of ingestion. Gastric volumes are technique specific and differ between MRI and SPECT.

Patient specific physical anatomy models.

The advent of small footprint stereo-lithographic printers and the ready availability of segmentation and surface modeling software provide a unique opportunity to create patient-specific physical models of anatomy, validation of image guided intervention applications against phantoms that exhibit naturally occurring anatomic variation. Because these models can incorporate all structures relevant to a procedure, this allows validation to occur under realistic conditions using the same or similar techniques as would be used in a clinical application. This in turn reduces the number of trials and time spent performing in-vivo validation experiments. In this paper, we describe our general approach for the creation of both non-tissue and tissue-mimicking patient-specific models as part of a general-purpose patient emulation system used to validate image guided intervention applications.

An integrated system for real-time image guided cardiac catheter ablation.

Minimally invasive cardiac catheter ablation procedures require effective visualization of the relevant heart anatomy and electrophysiology (EP). In a typical ablation procedure, the visualization tools available to the cardiologist include bi-plane fluoroscopy, real-time ultrasound, and a coarse 3D model which gives a rough representation of cardiac anatomy and electrical activity. Recently, there has been increased interest in incorporating detailed, patient specific anatomical data into the cardiac ablation procedure. We are currently developing a prototype system which both integrates a patient specific, preoperative data model into the procedure as well as fuses the various visualization modalities (i.e. fluoroscopy, ultrasound, EP) into a single display. In this paper, we focus on two aspects of the prototype system. First, we describe the framework for integrating the various system components, including an efficient communication protocol. Second, using a simple two-chamber phantom of the heart, we demonstrate the ability to integrate preoperative data into the ablation procedure. This involves the registration and visualization of tracked catheter points within the cardiac chambers of the preoperative model.

Quantitative characterization of lung disease.

The increase in prevalence, incidence and variety of pulmonary diseases has precipitated the need for more non-invasive quantitative assessment of structure/function relationships in the lung. This need requires concise description not only of lung anatomy but also of associated underlying mechanics of pulmonary function, as well as deviation from normal in specific diseases. This can be facilitated through the use of adaptive deformable surface models of the lung at end inspiratory and expiratory volumes. Lung surface deformation may be used to represent tissue excursion, which can characterize both global and regional lung mechanics. In this paper, we report a method for robust determination and visualization of pulmonary structure and function using clinical CT scans. The method provides both intuitive 3D parametric visualization and objective quantitative assessment of lung structure and associated function, in both normal and pathological cases.

Patient specific dynamic geometric models from sequential volumetric time series image data.

Generating patient specific dynamic models is complicated by the complexity of the motion intrinsic and extrinsic to the anatomic structures being modeled. Using a physics-based sequentially deforming algorithm, an anatomically accurate dynamic four-dimensional model can be created from a sequence of 3-D volumetric time series data sets. While such algorithms may accurately track the cyclic non-linear motion of the heart, they generally fail to accurately track extrinsic structural and non-cyclic motion. To accurately model these motions, we have modified a physics-based deformation algorithm to use a meta-surface defining the temporal and spatial maxima of the anatomic structure as the base reference surface. A mass-spring physics-based deformable model, which can expand or shrink with the local intrinsic motion, is applied to the metasurface, deforming this base reference surface to the volumetric data at each time point. As the meta-surface encompasses the temporal maxima of the structure, any extrinsic motion is inherently encoded into the base reference surface and allows the computation of the time point surfaces to be performed in parallel. The resultant 4-D model can be interactively transformed and viewed from different angles, showing the spatial and temporal motion of the anatomic structure. Using texture maps and per-vertex coloring, additional data such as physiological and/or biomechanical variables (e.g., mapping electrical activation sequences onto contracting myocardial surfaces) can be associated with the dynamic model, producing a 5-D model. For acquisition systems that may capture only limited time series data (e.g., only images at end-diastole/end-systole or inhalation/exhalation), this algorithm can provide useful interpolated surfaces between the time points. Such models help minimize the number of time points required to usefully depict the motion of anatomic structures for quantitative assessment of regional dynamics.

Parametric display of myocardial function.

Quantitative assessment of regional heart motion has significant potential to provide more specific diagnosis of cardiac disease and cardiac malfunction than currently possible. Local heart motion may be captured from various medical imaging scanners. In this study, 3-D reconstructions of pre-infarct and post-infarct hearts were obtained from the Dynamic Spatial Reconstructor (DSR)[Ritman EL, Robb RA, Harris LD. Imaging physiological functions: experience with DSR. Philadelphia: Praeger, 1985; Robb RA, Lent AH, Gilbert BK, Chu A. The dynamic spatial reconstructor: a computed tomography system for high-speed simultaneous scanning of multiple cross sections of the heart. J Med Syst 1980;4(2):253-88; Jorgensen SM, Whitlock SV, Thomas PJ, Roessler RW, Ritman EL. The dynamic spatial reconstructor: a high speed, stop action, 3-D, digital radiographic imager of moving internal organs and blood. Proceedings of SPIE, Ultrahigh- and High-speed Photography, Videography, Photonics, and Velocimetry 1990;1346:180-91.] (DSR). Using functional parametric mapping of disturbances in regional contractility and relaxation, regional myocardial motion during a cardiac cycle is color mapped onto a deformable heart model to facilitate appreciation of the structure-to-function relationships in the myocardium, such as occurs in regional patterns of akinesis or dyskinesis associated with myocardial ischemia or infarction resulting from coronary artery occlusion.

Measurement and display of regional myocardial motion during post infarct treatment.

Quantitative assessment of 3-D regional heart motion has significant potential to provide more specific diagnosis of cardiac malfunction than currently possible. Using functional parametric mapping, regional myocardial motion during a cardiac cycle can be color-mapped onto a deformable heart model to provide better understanding of the structure-to-function relationships in the myocardium, including regional patterns of akinesis or dyskinesis associated with ischemia or infarction. In this study, 3-D reconstructions of human hearts were obtained from Electron-Beam Computed Tomography [1] (EB-CT), comparing stages of treatment after myocardial infarction.

3D localization of implanted radioactive sources in the prostate using trans-urethral ultrasound.

Prostate cancer is the most common cancer diagnosed in men in the United States. Many techniques have been developed to diagnose and treat prostate cancer. 3D Trans-Urethral Ultrasound (TUUS) is a new technique for the diagnosis and treatment of prostate disease. This research focuses on the potential of TUUS for therapy-guidance during and after transperineal interstitial permanent prostate brachytherapy (TIPPB). Computed tomography (CT) is currently used to determine the source locations and the effective radiation dose distribution throughout the tissue after the completion of the procedure. TUUS may be a viable alternative to CT for determining source locations within the prostate. Placement of the TUUS catheter into the urethra provides excellent 2D images of the prostate. In addition, the catheter can be used to acquire 3D volumetric data for 3D analysis of the prostate, associated tissue, and radioactive sources. Initial work on source localization was conducted on an ultrasound-equivalent prostate phantom. Cylindrical dummy radiation sources with diameter of .8 mm and length of 4.5 mm were placed into the prostate phantom for assessment with TUUS. The TUUS imaging device is a 10 fr catheter with a linear array of ultrasonic crystals at one end. The ultrasound catheter operates at 10 MHz and is controlled by the Acuson Sequoia ultrasound workstation. The catheter was placed in the phantom urethra and 3D scans were acquired. A corresponding CT was acquired for comparative purposes. Segmentation of the prostate capsule was done semi-automatically. Dummy radiation seed segmentation was conducted manually. Additional processing was necessary to account for image artifact and correctly reconstruct the seeds. The prostate shell and radioactive source reconstructions provide an excellent 3D representation. Comparison to the CT data suggests that the TUUS data provided: 1) greater resolution and 2) better soft tissue differentiation. The reconstructed sources were measured and corresponded to the physical dimensions of those placed in the phantoms. The method is now being evaluated on cadavers and patients.

Patient specific physics-based model for interactive visualization of cardiac dynamics.

Cardiac disorders result mainly from defects in cardiac structure or failure to generate and regulate electrical impulses. Knowledge of the structure, motion patterns, local deformation, and associated electrical activation patterns of the heart is necessary for precise diagnosis and treatment. Electrical and mechanical performance of the heart is strongly influenced by the anisotropic nature of myocardial tissue. Diffusion-encoded MR imaging provides in vivo myocardial fiber track information that can be used for precise simulation of cardiac conduction and contraction. We propose a method that incorporates such fiber track information with a physics-based deformable model to realistically simulate cardiac contraction and subsequent relaxation. The simulation aims to reproduce the myocardial deformation during the heartbeat. The system allows interactive visualization of dynamic 3-D heart structures during the cardiac cycle. In procedures such as catheter ablation, the interactive 4-D model provides updated anatomy and physiology of the patient's heart simultaneously, and can be used to guide the procedure for efficient targeting of the treatment regions.

An axial skeleton based surface deformation algorithm for patient specific anatomic modeling.

Traditionally, finite element analysis or mass-spring systems are used to calculate deformations of geometric surfaces. Patient-specific geometric models can be comprised of tens of thousands, even hundreds of thousands of polygons, making finite element analysis and mass-spring systems computationally demanding. Simulations using deformable patient specific models at real time rates are prohibitive under such a computational burden. This paper presents a method for simulating deformable surfaces by deforming a skeletal representation of the surface, rather than the surface itself, yielding an efficient method for interactive simulation with models.

Voxel significance mapping using local image variances in subtraction ictal SPET.

Subtraction ictal SPET co-registered to MRI (SISCOM) has been shown to aid epileptogenic localization and improve surgical outcome in partial epilepsy patients. This paper reports a method of identifying significant areas of epileptogenic activation in the SISCOM subtraction image, taking into account normal variation between sequential 99Tcm-ethyl cysteinate diethylester SPET scans of single individuals, and attempts to assess the clinical value of statistical mapping in subtraction SPET. Non-linear inter-subject registration is used to combine a group of subtraction images into a common anatomical framework. A map of the pixel intensity standard deviation values in the subtraction images is created, and this map is non-linearly registered to a patient's SISCOM subtraction image. Pixels in the patient subtraction image were then evaluated based upon the statistical characteristics of corresponding pixels in the atlas. SISCOM images created with the voxel variance method were rated higher in quality than the conventional image variance method in 15 patients. No difference in localization rate was observed between the voxel variance mapping and image variance methods. The voxel significance mapping method was shown to improve the quality of clinical SISCOM images.

Virtual endoscopy: development and evaluation using the Visible Human datasets.

Virtual endoscopy (VE) is a new method of diagnosis using computer processing of 3D image datasets (such as CT or MRI scans) to provide simulated visualizations of patient specific organs similar or equivalent to those produced by standard endoscopic procedures. Conventional endoscopy is invasive and often uncomfortable for patients. It sometimes has serious side effects such as perforation, infection and hemorrhage. VE visualization avoids these risks and can minimize difficulties and decrease morbidity when used before actual endoscopic procedures. In addition, there are many body regions not compatible with real endoscopy that can be explored with VE. Eventually, VE may replace many forms of real endoscopy. There remains a critical need to refine and validate VE visualizations for routine clinical use. We have used the Visible Human Dataset from the National Library of Medicine to develop and test these procedures and to evaluate their use in a variety of clinical applications. We have developed specific clinical protocols to compare virtual endoscopy with real endoscopy. We have developed informative and dynamic on-screen navigation guides to help the surgeon or physician interactively determine body orientation and precise anatomical localization while performing the VE procedures. Additionally, the adjunctive value of full 3D imaging (e.g. looking "outside" of the normal field of view) during the VE exam is being evaluated. Quantitative analyses of local geometric and densitometric properties obtained from the virtual procedures ("virtual biopsy") are being developed and compared with other direct measures. Preliminary results suggest that these virtual procedures can provide accurate, reproducible and clinically useful visualizations and measurements. These studies will help drive improvements in and lend credibility to VE procedures and simulations as routine clinical tools. VE holds significant promise for optimizing endoscopic diagnostic procedures, minimizing patient risk and morbidity, and reducing health care costs.

Subtraction ictal SPECT coregistered to MRI for seizure focus localization in partial epilepsy.

Peri-ictal single-photon emission computed tomography (SPECT) of the brain is increasingly used in localizing the seizure focus in presurgical evaluation of patients with partial epilepsy. However, traditional side-by-side visual interpretation of ictal and interictal SPECT films is hampered by differences in slice location and tracer activity. Precise correlation of the seizure focus with a high-quality image of the underlying brain anatomy can improve the physician's understanding of seizure neurophysiology and assist in surgical planning. Computer-based methods have been developed for aligning, normalizing, and subtracting digital ictal and interictal SPECT images of the patient's brain to produce a map of the blood flow changes occurring between the seizure and resting states. These maps are then aligned with a high-resolution magnetic resonance image (MRI) of the patient's brain anatomy and fused to identify anatomical regions involved in the seizure. The purpose of this article is to review the technical components and clinical implementation of subtraction ictal SPECT, as well as to discuss recent technological advances that could extend and improve the diagnostic and localizing capacity of this method.