Automatic Localization of the Left Ventricle from Cardiac Cine Magnetic Resonance Imaging: A New Spectrum-Based Computer-Aided Tool

Traditionally, cardiac image analysis is done manually. Automatic image processing can help with the repetitive tasks, and also deal with huge amounts of data, a task which would be humanly tedious. This study aims to develop a spectrum-based computer-aided tool to locate the left ventricle using images obtained via cardiac magnetic resonance imaging. Discrete Fourier Transform was conducted pixelwise on the image sequence. Harmonic images of all frequencies were analyzed visually and quantitatively to determine different patterns of the left and right ventricles on spectrum. The first and fifth harmonic images were selected to perform an anisotropic weighted circle Hough detection. This tool was then tested in ten volunteers. Our tool was able to locate the left ventricle in all cases and had a significantly higher cropping ratio of 0.165 than did earlier studies. In conclusion, a new spectrum-based computer aided tool has been proposed and developed for automatic left ventricle localization. The development of this technique, which will enable the automatic location and further segmentation of the left ventricle, will have a significant impact in research and in diagnostic settings. We envisage that this automated method could be used by radiographers and cardiologists to diagnose and assess ventricular function in patients with diverse heart diseases.     

Electrical Wave Propagation in an Anisotropic Model of the Left Ventricle Based on Analytical Description of Cardiac Architecture

We develop a numerical approach based on our recent analytical model of fiber structure in the left ventricle of the human heart. A special curvilinear coordinate system is proposed to analytically include realistic ventricular shape and myofiber directions. With this anatomical model, electrophysiological simulations can be performed on a rectangular coordinate grid. We apply our method to study the effect of fiber rotation and electrical anisotropy of cardiac tissue (i.e., the ratio of the conductivity coefficients along and across the myocardial fibers) on wave propagation using the ten Tusscher–Panfilov (2006) ionic model for human ventricular cells. We show that fiber rotation increases the speed of cardiac activation and attenuates the effects of anisotropy. Our results show that the fiber rotation in the heart is an important factor underlying cardiac excitation. We also study scroll wave dynamics in our model and show the drift of a scroll wave filament whose velocity depends non-monotonically on the fiber rotation angle; the period of scroll wave rotation decreases with an increase of the fiber rotation angle; an increase in anisotropy may cause the breakup of a scroll wave, similar to the mother rotor mechanism of ventricular fibrillation.     

A Geometrical Approach for Automatic Shape Restoration of the Left Ventricle

This paper describes an automatic algorithm that uses a geometry-driven optimization approach to restore the shape of three-dimensional (3D) left ventricular (LV) models created from magnetic resonance imaging (MRI) data. The basic premise is to restore the LV shape such that the LV epicardial surface is smooth after the restoration and that the general shape characteristic of the LV is not altered. The Maximum Principle Curvature () and the Minimum Principle Curvature () of the LV epicardial surface are used to construct a shape-based optimization objective function to restore the shape of a motion-affected LV via a dual-resolution semi-rigid deformation process and a free-form geometric deformation process. A limited memory quasi-Newton algorithm, L-BFGS-B, is then used to solve the optimization problem. The goal of the optimization is to achieve a smooth epicardial shape by iterative in-plane and through-plane translation of vertices in the LV model. We tested our algorithm on 30 sets of LV models with simulated motion artifact generated from a very smooth patient sample, and 20 in vivo patient-specific models which contain significant motion artifacts. In the 30 simulated samples, the Hausdorff distances with respect to the Ground Truth are significantly reduced after restoration, signifying that the algorithm can restore geometrical accuracy of motion-affected LV models. In the 20 in vivo patient-specific models, the results show that our method is able to restore the shape of LV models without altering the general shape of the model. The magnitudes of in-plane translations are also consistent with existing registration techniques and experimental findings.     

Automatic Segmentation of Eight Tissue Classes in Neonatal Brain MRI

Purpose

Volumetric measurements of neonatal brain tissues may be used as a biomarker for later neurodevelopmental outcome. We propose an automatic method for probabilistic brain segmentation in neonatal MRIs.

Materials and Methods

In an IRB-approved study axial T1- and T2-weighted MR images were acquired at term-equivalent age for a preterm cohort of 108 neonates. A method for automatic probabilistic segmentation of the images into eight cerebral tissue classes was developed: cortical and central grey matter, unmyelinated and myelinated white matter, cerebrospinal fluid in the ventricles and in the extra cerebral space, brainstem and cerebellum. Segmentation is based on supervised pixel classification using intensity values and spatial positions of the image voxels. The method was trained and evaluated using leave-one-out experiments on seven images, for which an expert had set a reference standard manually. Subsequently, the method was applied to the remaining 101 scans, and the resulting segmentations were evaluated visually by three experts. Finally, volumes of the eight segmented tissue classes were determined for each patient.

Results

The Dice similarity coefficients of the segmented tissue classes, except myelinated white matter, ranged from 0.75 to 0.92. Myelinated white matter was difficult to segment and the achieved Dice coefficient was 0.47. Visual analysis of the results demonstrated accurate segmentations of the eight tissue classes. The probabilistic segmentation method produced volumes that compared favorably with the reference standard.

Conclusion

The proposed method provides accurate segmentation of neonatal brain MR images into all given tissue classes, except myelinated white matter. This is the one of the first methods that distinguishes cerebrospinal fluid in the ventricles from cerebrospinal fluid in the extracerebral space. This method might be helpful in predicting neurodevelopmental outcome and useful for evaluating neuroprotective clinical trials in neonates.     

Dynamic Dyssynchrony and Impaired Contractile Reserve of the Left Ventricle in Beta-Thalassaemia Major: An Exercise Echocardiographic Study

Background

Performance of the left ventricle during exercise stress in thalassaemia patients is uncertain. We aimed to explore the phenomenon of dynamic dyssynchrony and assess contractile reserve in patients with beta-thalassaemia major and determine their relationships with myocardial iron load.

Methods and Results

Thirty-two thalassaemia patients (16 males), aged 26.8±6.9 years, without heart failure and 17 healthy controls were studied. Their left ventricular (LV) volumes, ejection fraction, systolic dyssynchrony index (SDI), and myocardial acceleration during isovolumic LV contraction (IVA) were determined at rest and during submaximal bicycle exercise testing using 3-dimensional and tissue Doppler echocardiography. Myocardial iron load as assessed by T2* cardiac magnetic resonance in patients were further related to indices of LV dyssynchrony and contractile reserve. At rest, patients had significantly greater LV SDI (p<0.001) but similar IVA (p = 0.22) compared with controls. With exercise stress, the prevalence of mechanical dyssynchrony (SDI>4.6%, control+2SD) increased from baseline 25% to 84% in patients. Δ SDI_{exercise-baseline} correlated with exercise-baseline differences in LV ejection fraction (p<0.001) and stroke volume (p = 0.006). Compared with controls, patients had significantly less exercise-induced increase in LV ejection fraction, cardiac index, and IVA (interaction, all p<0.05) and had impaired contractile reserve as reflected by the gentler IVA-heart rate slope (p = 0.018). Cardiac T2* in patients correlated with baseline LV SDI (r = −0.44, p = 0.011) and IVA-heart rate slope (r = 0.36, p = 0.044).

Conclusions

Resting LV dyssynchrony is associated with myocardial iron load. Exercise stress further unveils LV dynamic dyssynchrony and impaired contractile reserve in patients with beta-thalassaemia major.     

A Tissue-Level Electromechanical Model of the Left Ventricle: Application to the Analysis of Intraventricular Pressure

The ventricular pressure profile is characteristic of the cardiac contraction progress and is useful to evaluate the cardiac performance. In this contribution, a tissue-level electromechanical model of the left ventricle is proposed, to assist the interpretation of left ventricular pressure waveforms. The left ventricle has been modeled as an ellipsoid composed of twelve mechano-hydraulic sub-systems. The asynchronous contraction of these twelve myocardial segments has been represented in order to reproduce a realistic pressure profiles. To take into account the different energy domains involved, the tissue-level scale and to facilitate the building of a modular model, multiple formalisms have been used: Bond Graph formalism for the mechano-hydraulic aspects and cellular automata for the electrical activation. An experimental protocol has been defined to acquire ventricular pressure signals from three pigs, with different afterload conditions. Evolutionary Algorithms have been used to identify the model parameters in order to minimize the error between experimental and simulated ventricular pressure signals. Simulation results show that the model is able to reproduce experimental ventricular pressure. In addition, electro-mechanical activation times have been determined in the identification process. For example, the maximum electrical activation time is reached, respectively, 96.5, 139.3 and 131.5 ms for the first, second, and third pigs. These preliminary results are encouraging for the application of the model on non-invasive data like ECG, arterial pressure or myocardial strain.     

Stress Distribution in the Canine Left Ventricle during Diastole and Systole

A model is proposed for stress analysis of the left ventricular wall (LV wall) based on the realistic assumption that the myocardium is essentially composed of fiber elements which carry only axial tension and vary in orientation through the wall. Stress analysis based on such a model requires an extensive study of muscle fiber orientation and curvature through the myocardium. Accordingly, the principal curvatures were studied at a local site near the equator in ten dog hearts rapidly fixed in situ at end diastole and end systole; the fiber orientation for these hearts had already been established in a previous study. The principal radii of curvature were (a) measured by fitting templates to the endocardial and epicardial wall surfaces in the circumferential and longitudinal directions and (b) computed from measured lengths of semiaxes of ellipsoids of revolution representing the LV wall (“ellipsoid” data). The wall was regarded as a tethered set of nested shells, each having a unique fiber orientation. Results indicate the following. (a) Fiber curvature, k, is maximum at midwall at end systole; this peak shifts towards endocardium at end diastole. (b) The pressure or radial stress through the wall decreases more rapidly near the endocardium than near the epicardium at end diastole and at end systole when a constant tension is assumed for each fiber through the wall. (c) At end diastole the curve for the circumferential stress vs. wall thickness is convex with a maximum at midwall. In the longitudinal direction the stress distribution curve is concave with a minimum at midwall. Similar distributions are obtained at end systole when a constant tension is assumed for each fiber through the wall. (d) The curvature and stress distributions obtained by direct measurements at a selected local site agree well with those computed from “ellipsoid” data.     

Cerebral Blood Volume Analysis in Glioblastomas Using Dynamic Susceptibility Contrast-Enhanced Perfusion MRI: A Comparison of Manual and Semiautomatic Segmentation Methods

Purpose

To compare the reproducibilities of manual and semiautomatic segmentation method for the measurement of normalized cerebral blood volume (nCBV) using dynamic susceptibility contrast-enhanced (DSC) perfusion MR imaging in glioblastomas.

Materials and Methods

Twenty-two patients (11 male, 11 female; 27 tumors) with histologically confirmed glioblastoma (WHO grade IV) were examined with conventional MR imaging and DSC imaging at 3T before surgery or biopsy. Then nCBV (means and standard deviations) in each mass was measured using two DSC MR perfusion analysis methods including manual and semiautomatic segmentation method, in which contrast-enhanced (CE)-T1WI and T2WI were used as structural imaging. Intraobserver and interobserver reproducibility were assessed according to each perfusion analysis method or each structural imaging. Interclass correlation coefficient (ICC), Bland-Altman plot, and coefficient of variation (CV) were used to evaluate reproducibility.

Results

Intraobserver reproducibilities on CE-T1WI and T2WI were ICC of 0.74–0.89 and CV of 20.39–36.83% in manual segmentation method, and ICC of 0.95–0.99 and CV of 8.53–16.19% in semiautomatic segmentation method, repectively. Interobserver reproducibilites on CE-T1WI and T2WI were ICC of 0.86–0.94 and CV of 19.67–35.15% in manual segmentation method, and ICC of 0.74–1.0 and CV of 5.48–49.38% in semiautomatic segmentation method, respectively. Bland-Altman plots showed a good correlation with ICC or CV in each method. The semiautomatic segmentation method showed higher intraobserver and interobserver reproducibilities at CE-T1WI-based study than other methods.

Conclusion

The best reproducibility was found using the semiautomatic segmentation method based on CE-T1WI for structural imaging in the measurement of the nCBV of glioblastomas.     

Quantitative Estimation of Renal Function with Dynamic Contrast-Enhanced MRI Using a Modified Two-Compartment Model

Objective

To establish a simple two-compartment model for glomerular filtration rate (GFR) and renal plasma flow (RPF) estimations by dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI).

Materials and Methods

A total of eight New Zealand white rabbits were included in DCE-MRI. The two-compartment model was modified with the impulse residue function in this study. First, the reliability of GFR measurement of the proposed model was compared with other published models in Monte Carlo simulation at different noise levels. Then, functional parameters were estimated in six healthy rabbits to test the feasibility of the new model. Moreover, in order to investigate its validity of GFR estimation, two rabbits underwent acute ischemia surgical procedure in unilateral kidney before DCE-MRI, and pixel-wise measurements were implemented to detect the cortical GFR alterations between normal and abnormal kidneys.

Results

The lowest variability of GFR and RPF measurements were found in the proposed model in the comparison. Mean GFR was 3.03±1.1 ml/min and mean RPF was 2.64±0.5 ml/g/min in normal animals, which were in good agreement with the published values. Moreover, large GFR decline was found in dysfunction kidneys comparing to the contralateral control group.

Conclusion

Results in our study demonstrate that measurement of renal kinetic parameters based on the proposed model is feasible and it has the ability to discriminate GFR changes in healthy and diseased kidneys.     

Quantitative Study of the Effect of Tissue Microstructure on Contraction in a Computational Model of Rat Left Ventricle

Tissue microstructure, in particular the alignment of myocytes (fibre direction) and their lateral organisation into sheets, is fundamental to cardiac function. We studied the effect of microstructure on contraction in a computational model of rat left ventricular electromechanics. Different fibre models, globally rule-based or locally optimised to DT-MRI data, were compared, in order to understand whether a subject-specific fibre model would enhance the predictive power of our model with respect to the global ones. We also studied the impact of sheets on ventricular deformation by comparing: (a) a transversely isotropic versus an orthotropic material law and (b) a linear model with a bimodal model of sheet transmural variation. We estimated ejection fraction, wall thickening and base-to-apex shortening and compared them with measures from cine-MRI. We also evaluated Lagrangian strains as local metrics of cardiac deformation. Our results show that the subject-specific fibre model provides little improvement in the metric predictions with respect to global fibre models while material orthotropy allows closer agreement with measures than transverse isotropy. Nonetheless, the impact of sheets in our model is smaller than that of fibres. We conclude that further investigation of the modelling of sheet dynamics is necessary to fully understand the impact of tissue structure on cardiac deformation.     

Evaluation of Global Longitudinal Strain of Left Ventricle and Regional Longitudinal Strain in the Region of Left Ventricular Leads Predicts the Response to Cardiac Resynchronization Therapy in Patients with Ischemic Heart Failure

Myocardium viability in ischemic heart failure (HF) may affect the effect of cardiac resynchronization therapy (CRT). We hypothesized that longitudinal strain of 2D-STE, which reflects myocardium viability, can predict the response to CRT in patients with ischemic HF. 2D-STE was performed in 42 patients with HF, 1 week before and 1 year after CRT. GLS, RLS, and the LV synchrony index (SI), defined as the difference in timing to peak radial strain between LV anterior septal and posterior wall in LV short axis view, were calculated. A decrease in the LV end-systolic volume (LVESV) value of ≥15 % 1 year after CRT was defined as response to CRT. Twenty-nine patients responded to CRT (CRT-R group), while 13 patients did not respond and were assigned as CRT-NR group. Pre-CRT RLS and GLS were higher, while SI is lower, in CRT-R patients compared with CRT-NR group (p < 0.001). The ROC curve revealed that RLS of ?11.5 % predicted response to CRT with sensitivity of 80.0 % and specificity of 77.9 % (AUC = 0.84, p < 0.001). Further, GLS of ?13 % predicted response to CRT with sensitivity of 73.0 % and specificity of 73.4 % (AUC = 0.79, p < 0.001). In conclusion, LV dyssynchrony, GLS, and RLS calculated by 2D-STE can predict long-term response to CRT in patients with ischemic HF.     

A Model of the Structural and Functional Development of the Normal Human Fetal Left Ventricle Based on a Global Growth Law

The purpose of this research is to study the growth of the normal human left ventricle (LV) during the fetal period from 14 to 40 weeks of gestation. A new constitutive law for the active myocardium describing the mechanical properties of the active muscle during the whole cardiac cycle has been proposed. The LV model is a thick-walled, incompressible, hyperelastic cylinder, with families of helicoidal fibers running on cylindrical surfaces [1] . Based on the works of Lin and Taber [2] done on the embryonic chick heart, we use for the human fetal heart a growth law in which the growth rate depends on the wall stresses. The parameters of the growth law are adapted to agree with sizes and volumes inferred from two dimensional ultrasound measurements performed on 18 human fetuses. Then calculations are performed to extrapolate the cardiac performance during normal growth of the fetal LV. The results presented support the idea that a growth law in which the growth rate depends linearly on the mean wall stresses averaged through the space and during whole cardiac cycle, is adapted to the normal human fetal LV development.     

Phenotyping of Left and Right Ventricular Function in Mouse Models of Compensated Hypertrophy and Heart Failure with Cardiac MRI

Background

Left ventricular (LV) and right ventricular (RV) function have an important impact on symptom occurrence, disease progression and exercise tolerance in pressure overload-induced heart failure, but particularly RV functional changes are not well described in the relevant aortic banding mouse model. Therefore, we quantified time-dependent alterations in the ventricular morphology and function in two models of hypertrophy and heart failure and we studied the relationship between RV and LV function during the transition from hypertrophy to heart failure.

Methods

MRI was used to quantify RV and LV function and morphology in healthy (n = 4) and sham operated (n = 3) C57BL/6 mice, and animals with a mild (n = 5) and a severe aortic constriction (n = 10).

Results

Mice subjected to a mild constriction showed increased LV mass (P<0.01) and depressed LV ejection fraction (EF) (P<0.05) as compared to controls, but had similar RVEF (P>0.05). Animals with a severe constriction progressively developed LV hypertrophy (P<0.001), depressed LVEF (P<0.001), followed by a declining RVEF (P<0.001) and the development of pulmonary remodeling, as compared to controls during a 10-week follow-up. Myocardial strain, as a measure for local cardiac function, decreased in mice with a severe constriction compared to controls (P<0.05).

Conclusions

Relevant changes in mouse RV and LV function following an aortic constriction could be quantified using MRI. The well-controlled models described here open opportunities to assess the added value of new MRI techniques for the diagnosis of heart failure and to study the impact of new therapeutic strategies on disease progression and symptom occurrence.     

Distinctive Profile of IsomiR Expression and Novel MicroRNAs in Rat Heart Left Ventricle

MicroRNAs (miRNAs) are single-stranded non-coding RNAs that negatively regulate target gene expression through mRNA cleavage or translational repression. There is mounting evidence that they play critical roles in heart disease. The expression of known miRNAs in the heart has been studied at length by microarray and quantitative PCR but it is becoming evident that microRNA isoforms (isomiRs) are potentially physiologically important. It is well known that left ventricular (patho)physiology is influenced by transmural heterogeneity of cardiomyocyte phenotype, and this likely reflects underlying heterogeneity of gene expression. Given the significant role of miRNAs in regulating gene expression, knowledge of how the miRNA profile varies across the ventricular wall will be crucial to better understand the mechanisms governing transmural physiological heterogeneity. To determinine miRNA/isomiR expression profiles in the rat heart we investigated tissue from different locations across the left ventricular wall using deep sequencing. We detected significant quantities of 145 known rat miRNAs and 68 potential novel orthologs of known miRNAs, in mature, mature* and isomiR formation. Many isomiRs were detected at a higher frequency than their canonical sequence in miRBase and have different predicted targets. The most common miR-133a isomiR was more effective at targeting a construct containing a sequence from the gelsolin gene than was canonical miR-133a, as determined by dual-fluorescence assay. We identified a novel rat miR-1 homolog from a second miR-1 gene; and a novel rat miRNA similar to miR-676. We also cloned and sequenced the rat miR-486 gene which is not in miRBase (v18). Signalling pathways predicted to be targeted by the most highly detected miRNAs include Ubiquitin-mediated Proteolysis, Mitogen-Activated Protein Kinase, Regulation of Actin Cytoskeleton, Wnt signalling, Calcium Signalling, Gap junctions and Arrhythmogenic Right Ventricular Cardiomyopathy. Most miRNAs are not expressed in a gradient across the ventricular wall, with exceptions including miR-10b, miR-21, miR-99b and miR-486.     

The Nature of the Disorder of Function in Chronic Postinfarction Aneurysm of the Left Ventricle

Assessment of left ventricular function in five patients with chronic postinfarction left ventricular aneurysm was carried out at the time of left heart catheterization and compared with that in six normal subjects. One patient was investigated before and after surgical resection of the aneurysm. The presence of the aneurysm placed the left ventricle at a mechanical disadvantage in systole and increased the resistance to diastolic filling (impedance). This was true even in patients with normal cardiac indices who were not badly disabled. Resection of the aneurysm corrected both these abnormalities, and, as well, lowered the time-tension index at a time when calculated left ventricular work was much increased. These differences between normals and patients with aneurysms, and the changes occurring as a result of resection of an aneurysm, show that the presence of the aneurysm places the left ventricle at a mechanical disadvantage in systole as well as altering its diastolic filling characteristics.     

Repeatability of Cerebral Perfusion Using Dynamic Susceptibility Contrast MRI in Glioblastoma Patients

OBJECTIVES

This study evaluates the repeatability of brain perfusion using dynamic susceptibility contrast magnetic resonance imaging (DSC-MRI) with a variety of post-processing methods.

METHODS

Thirty-two patients with newly diagnosed glioblastoma were recruited. On a 3-T MRI using a dual-echo, gradient-echo spin-echo DSC-MRI protocol, the patients were scanned twice 1 to 5 days apart. Perfusion maps including cerebral blood volume (CBV) and cerebral blood flow (CBF) were generated using two contrast agent leakage correction methods, along with testing normalization to reference tissue, and application of arterial input function (AIF). Repeatability of CBV and CBF within tumor regions and healthy tissues, identified by structural images, was assessed with intra-class correlation coefficients (ICCs) and repeatability coefficients (RCs). Coefficients of variation (CVs) were reported for selected methods.

RESULTS

CBV and CBF were highly repeatable within tumor with ICC values up to 0.97. However, both CBV and CBF showed lower ICCs for healthy cortical tissues (up to 0.83), healthy gray matter (up to 0.95), and healthy white matter (WM; up to 0.93). The values of CV ranged from 6% to 10% in tumor and 3% to 11% in healthy tissues. The values of RC relative to the mean value of measurement within healthy WM ranged from 22% to 42% in tumor and 7% to 43% in healthy tissues. These percentages show how much variation in perfusion parameter, relative to that in healthy WM, we expect to observe to consider it statistically significant. We also found that normalization improved repeatability, but AIF deconvolution did not.

CONCLUSIONS

DSC-MRI is highly repeatable in high-grade glioma patients.     

Simulation of Left Atrial Function Using a Multi-Scale Model of the Cardiovascular System

During a full cardiac cycle, the left atrium successively behaves as a reservoir, a conduit and a pump. This complex behavior makes it unrealistic to apply the time-varying elastance theory to characterize the left atrium, first, because this theory has known limitations, and second, because it is still uncertain whether the load independence hypothesis holds. In this study, we aim to bypass this uncertainty by relying on another kind of mathematical model of the cardiac chambers. In the present work, we describe both the left atrium and the left ventricle with a multi-scale model. The multi-scale property of this model comes from the fact that pressure inside a cardiac chamber is derived from a model of the sarcomere behavior. Macroscopic model parameters are identified from reference dog hemodynamic data. The multi-scale model of the cardiovascular system including the left atrium is then simulated to show that the physiological roles of the left atrium are correctly reproduced. This include a biphasic pressure wave and an eight-shaped pressure-volume loop. We also test the validity of our model in non basal conditions by reproducing a preload reduction experiment by inferior vena cava occlusion with the model. We compute the variation of eight indices before and after this experiment and obtain the same variation as experimentally observed for seven out of the eight indices. In summary, the multi-scale mathematical model presented in this work is able to correctly account for the three roles of the left atrium and also exhibits a realistic left atrial pressure-volume loop. Furthermore, the model has been previously presented and validated for the left ventricle. This makes it a proper alternative to the time-varying elastance theory if the focus is set on precisely representing the left atrial and left ventricular behaviors.     

Characterization of Respiratory and Cardiac Motion from Electro-Anatomical Mapping Data for Improved Fusion of MRI to Left Ventricular Electrograms

Accurate fusion of late gadolinium enhancement magnetic resonance imaging (MRI) and electro-anatomical voltage mapping (EAM) is required to evaluate the potential of MRI to identify the substrate of ventricular tachycardia. However, both datasets are not acquired at the same cardiac phase and EAM data is corrupted with respiratory motion limiting the accuracy of current rigid fusion techniques. Knowledge of cardiac and respiratory motion during EAM is thus required to enhance the fusion process. In this study, we propose a novel approach to characterize both cardiac and respiratory motion from EAM data using the temporal evolution of the 3D catheter location recorded from clinical EAM systems. Cardiac and respiratory motion components are extracted from the recorded catheter location using multi-band filters. Filters are calibrated for each EAM point using estimates of heart rate and respiratory rate. The method was first evaluated in numerical simulations using 3D models of cardiac and respiratory motions of the heart generated from real time MRI data acquired in 5 healthy subjects. An accuracy of 0.6–0.7 mm was found for both cardiac and respiratory motion estimates in numerical simulations. Cardiac and respiratory motions were then characterized in 27 patients who underwent LV mapping for treatment of ventricular tachycardia. Mean maximum amplitude of cardiac and respiratory motion was 10.2±2.7 mm (min = 5.5, max = 16.9) and 8.8±2.3 mm (min = 4.3, max = 14.8), respectively. 3D Cardiac and respiratory motions could be estimated from the recorded catheter location and the method does not rely on additional imaging modality such as X-ray fluoroscopy and can be used in conventional electrophysiology laboratory setting.