Studies of the effects of gravitational and inertial forces on cardiovascular and respiratory dynamics

The current status and application are described of the biplane video roentgen densitometry, videometry and video digitization systems. These techniques were developed, and continue to be developed for studies of the effects of gravitational and inertial forces on cardiovascular and respiratory dynamics in intact animals and man. Progress is reported in the field of lung dynamics and three-dimensional reconstruction of the dynamic thoracic contents from roentgen video images. It is anticipated that these data will provide added insight into the role of shape and internal spatial relationships (which is altered particularly by acceleration and position of the body) of these organs as an indication of their functional status

3D/2D Registration of Mapping Catheter Images for Arrhythmia Interventional Assistance

Radiofrequency (RF) catheter ablation has transformed treatment for tachyarrhythmias and has become first-line therapy for some tachycardias. The precise localization of the arrhythmogenic site and the positioning of the RF catheter over that site are problematic: they can impair the efficiency of the procedure and are time consuming (several hours). Electroanatomic mapping technologies are available that enable the display of the cardiac chambers and the relative position of ablation lesions. However, these are expensive and use custom-made catheters. The proposed methodology makes use of standard catheters and inexpensive technology in order to create a 3D volume of the heart chamber affected by the arrhythmia. Further, we propose a novel method that uses a priori 3D information of the mapping catheter in order to estimate the 3D locations of multiple electrodes across single view C-arm images. The monoplane algorithm is tested for feasibility on computer simulations and initial canine data.Comment: International Journal of Computer Science Issues, IJCSI, Volume 4, Issue 2, pp10-19, September 200

Semiautomated Skeletonization of the Pulmonary Arterial Tree in Micro-CT Images

We present a simple and robust approach that utilizes planar images at different angular rotations combined with unfiltered back-projection to locate the central axes of the pulmonary arterial tree. Three-dimensional points are selected interactively by the user. The computer calculates a sub- volume unfiltered back-projection orthogonal to the vector connecting the two points and centered on the first point. Because more x-rays are absorbed at the thickest portion of the vessel, in the unfiltered back-projection, the darkest pixel is assumed to be the center of the vessel. The computer replaces this point with the newly computer-calculated point. A second back-projection is calculated around the original point orthogonal to a vector connecting the newly-calculated first point and user-determined second point. The darkest pixel within the reconstruction is determined. The computer then replaces the second point with the XYZ coordinates of the darkest pixel within this second reconstruction. Following a vector based on a moving average of previously determined 3- dimensional points along the vessel\u27s axis, the computer continues this skeletonization process until stopped by the user. The computer estimates the vessel diameter along the set of previously determined points using a method similar to the full width-half max algorithm. On all subsequent vessels, the process works the same way except that at each point, distances between the current point and all previously determined points along different vessels are determined. If the difference is less than the previously estimated diameter, the vessels are assumed to branch. This user/computer interaction continues until the vascular tree has been skeletonized

Three-dimensional reconstruction of myocardial contrast perfusion from biplane cineangiograms by means of linear programming techniques

The assessment of coronary flow reserve from the instantaneous distribution of the contrast agent within the coronary vessels and myocardial muscle at the control state and at maximal flow has been limited by the superimposition of myocardial regions of interest in the two-dimensional images. To overcome these limitations, we are in the process of developing a three-dimensional (3D) reconstruction technique to compute the contrast distribution in cross sections of the myocardial muscle from two orthogonal cineangiograms. To limit the number of feasible solutions in the 3D-reconstruction space, the 3D-geometry of the endo- and epicardial boundaries of the myocardium must be determined. For the geometric reconstruction of the epicardium, the centerlines of the left coronary arterial tree are manually or automatically traced in the biplane views. Next, the bifurcations are detected automatically and matched in these two views, allowing a 3D-representation of the coronary tree. Finally, the circumference of the left ventricular myocardium in a selected cross section can be computed from the intersection points of this cross section with the 3D coronary tree using B-splines. For the geometric reconstruction of the left ventricular cavity, we envision to apply the elliptical approximation technique using the LV boundaries defined in the two orthogonal views, or by applying more complex 3D-reconstruction techniques including densitometry. The actual 3D-reconstruction of the contrast distribution in the myocardium is based on a linear programming technique (Transportation model) using cost coefficient matrices. Such a cost coefficient matrix must contain a maximum amount of a priori information, provided by a computer generated model and updated with actual data from the angiographic views. We have only begun to solve this complex problem. However, based on our first experimental results we expect that the linear programming approach with advanced cost coefficient matrices and computed model will lead to a

Analysis of myocardial motion using generalized spline models and tagged magnetic resonance images

Heart wall motion abnormalities are the very sensitive indicators of common heart diseases, such as myocardial infarction and ischemia. Regional strain analysis is especially important in diagnosing local abnormalities and mechanical changes in the myocardium. In this work, we present a complete method for the analysis of cardiac motion and the evaluation of regional strain in the left ventricular wall. The method is based on the generalized spline models and tagged magnetic resonance images (MRI) of the left ventricle. The whole method combines dynamical tracking of tag deformation, simulating cardiac movement and accurately computing the regional strain distribution. More specifically, the analysis of cardiac motion is performed in three stages. Firstly, material points within the myocardium are tracked over time using a semi-automated snake-based tag tracking algorithm developed for this purpose. This procedure is repeated in three orthogonal axes so as to generate a set of one-dimensional sample measurements of the displacement field. The 3D-displacement field is then reconstructed from this sample set by using a generalized vector spline model. The spline reconstruction of the displacement field is explicitly expressed as a linear combination of a spline kernel function associated with each sample point and a polynomial term. Finally, the strain tensor (linear or nonlinear) with three direct components and three shear components is calculated by applying a differential operator directly to the displacement function. The proposed method is computationally effective and easy to perform on tagged MR images. The preliminary study has shown potential advantages of using this method for the analysis of myocardial motion and the quantification of regional strain

Mapping Molecular Agents Distributions in Whole Mice Hearts Using Born-Normalized Optical Projection Tomography

To date there is a lack of tools to map the spatio-temporal dynamics of diverse cells in experimental heart models. Conventional histology is labor intensive with limited coverage, whereas many imaging techniques do not have sufficiently high enough spatial resolution to map cell distributions. We have designed and built a high resolution, dual channel Born-normalized near-infrared fluorescence optical projection tomography system to quantitatively and spatially resolve molecular agents distribution within whole murine heart. We validated the use of the system in a mouse model of monocytes/macrophages recruitment during myocardial infarction. While acquired, data were processed and reconstructed in real time. Tomographic analysis and visualization of the key inflammatory components were obtained via a mathematical formalism based on left ventricular modeling. We observed extensive monocyte recruitment within and around the infarcted areas and discovered that monocytes were also extensively recruited into non-ischemic myocardium, beyond that of injured tissue, such as the septum

Three-dimensional reconstruction of stenosed coronary artery segments with assessment of the flow impedance

In this paper preliminary results of a study about the diagnostic benefits of 3D visualization and quantitation of stenosed coronary artery segments are presented. As is well known, even biplane angiographic images do not provide enough information for binary reconstruction. Therefore,a priori information about the slice to be reconstructed must be incorporated into the reconstruction algorithm. One approach is to assume a circular cross-section of the coronary artery. Hence, the diameter is estimated from the contours of the vessels in both projections. Another approach is to search for a solution of the reconstruction problem close to the previously reconstructed adjacent slice. In this paper we follow the first method based on contour information. The reconstructed coronary segment is visualized in three dimensions. Based on the obtained geometry of the obstruction the pertinent blood flow impedance is estimated on the basis of fluid dynamic principles. The results of applying the reconstruction algorithms to clinical coronary biplane exposures are presented with an indication of the assessed flow impedance

4D Flow Patterns and Relative Pressure Distribution in a Left Ventricle Model by Shake-the-Box and Proper Orthogonal Decomposition Analysis

Purpose: Intraventricular blood flow dynamics are associated with cardiac function. Accurate, noninvasive, and easy assessments of hemodynamic quantities (such as velocity, vortex, and pressure) could be an important addition to the clinical diagnosis and treatment of heart diseases. However, the complex time-varying flow brings many challenges to the existing noninvasive image-based hemodynamic assessments. The development of reliable techniques and analysis tools is essential for the application of hemodynamic biomarkers in clinical practice. Methods: In this study, a time-resolved particle tracking method, Shake-the-Box, was applied to reconstruct the flow in a realistic left ventricle (LV) silicone model with biological valves. Based on the obtained velocity, 4D pressure field was calculated using a Poisson equation-based pressure solver. Furthermore, flow analysis by proper orthogonal decomposition (POD) of the 4D velocity field has been performed. Results: As a result of the Shake-the-Box algorithm, we have extracted: (i) particle positions, (ii) particle tracks, and finally, (iii) 4D velocity fields. From the latter, the temporal evolution of the 3D pressure field during the full cardiac cycle was obtained. The obtained maximal pressure difference extracted along the base-to-apex was about 2.7 mmHg, which is in good agreement with those reported in vivo. The POD analysis results showed a clear picture of different scale of vortices in the pulsatile LV flow, together with their time-varying information and corresponding kinetic energy content. To reconstruct 95% of the kinetic energy of the LV flow, only the first six POD modes would be required, leading to significant data reduction. Conclusions: This work demonstrated Shake-the-Box is a promising technique to accurately reconstruct the left ventricle flow field in vitro. The good spatial and temporal resolutions of the velocity measurements enabled a 4D reconstruction of the pressure field in the left ventricle. The application of POD analysis showed its potential in reducing the complexity of the high-resolution left ventricle flow measurements. For future work, image analysis, multi-modality flow assessments, and the development of new flow-derived biomarkers can benefit from fast and data-reducing POD analysis.</p