This invention relates to the determination and representation of physiological information relating to a heart surface using, e.g., a non-contact catheter.
Use of minimally invasive procedures, such as catheter ablation, to treat a variety of heart conditions, such as supraventricular and ventricular arrhythmias, is becoming increasingly more prevalent. Such procedures involve the mapping of electrical activity in the heart (e.g., based on cardiac signals), such as at various locations on the endocardium surface (“cardiac mapping”), to identify the site of origin of the arrhythmia followed by a targeted ablation of the site. To perform such cardiac mapping a catheter with one or more electrodes can be inserted into the patient's heart chamber.
Conventional 3D mapping techniques include contact mapping and non-contact mapping. In contact mapping techniques one or more catheters are advanced into the heart. Physiological signals resulting from the electrical activity of the heart are acquired with one or more electrodes located at the catheter distal tip after determining that the tip is in stable and steady contact with the endocardium surface of a particular heart chamber. Location and electrical activity is usually measured sequentially on a point-by-point basis at about 50 to 200 points on the internal surface of the heart to construct an electro-anatomical depiction of the heart. The generated map may then serve as the basis for deciding on a therapeutic course of action, for example, tissue ablation, to alter the propagation of the heart's electrical activity and to restore normal heart rhythm. On the other hand, in non-contact-based mapping systems a multiple electrodes catheter is percutaneously placed in the heart chamber of interest. Once in the chamber, the catheter is deployed to assume a 3D shape. Using the signals detected by the non-contact electrodes and information on chamber anatomy and relative electrode location, the system provides physiological information regarding the endocardium of the heart chamber.
In some aspects a method includes measuring physiological signals at one or more electrodes in response to electrical activity in a heart cavity. The hear cavity having a surface. The method also includes determining, based at least in part on Laplace's equation, bipolar physiological information at multiple locations of a surface based on the measured physiological signals and positions of the one or more electrodes with respect to the surface.
Embodiments can include one or more of the following.
The physiological signals can be unipolar signals such as unipolar potential signals.
Measuring the physiological signals can include measuring a potential between a first electrode and a second electrode. The first electrode can be located within the heart cavity and the second electrode can be located at a distance from the heart cavity such that the electrode is not affected by local tissue activation in a heart cavity.
Measuring the physiological signals can include measuring a potential between the first electrode and the second electrode that are separated by a distance greater than 5 cm. Measuring the physiological signals can include measuring a potential between the first electrode and the second electrode that are separated by a distance greater than 10 cm. Measuring the physiological signals can include measuring a potential between the first electrode and the second electrode that are separated by a distance greater than 15 cm.
The bipolar physiological information can be current density information. The current density information can include a normal component of the current density. The current density information can include a magnitude of a current density vector. The current density information can include a magnitude of a tangential component of the current density.
Measuring the physiological signals can include measuring the physiological signals for multiple different catheter positions in the heart cavity. The number of catheter positions at which the signals are measured and used to determine the bipolar physiological information at the multiple locations of the endocardium surface can be more than three.
The signals can be measured for at least one electrical heart cycle for each catheter position.
The determination of the physiological information at the multiple locations of the endocardium surface can include synchronizing the signals measured at the different catheter positions with one another according to an electrical heart beat cycle. The measured signals can be synchronized based on physiological data including at least one of: ECG and intercardiac electrograms.
The determination of the physiological information at the multiple locations of the endocardium surface can include processing the synchronized signals as though they were obtained at one time from all of the positions sampled by the catheter electrodes for the different positions of the catheter in the heart cavity.
Determining the bipolar physiological information can include applying a transformation function to the measured physiological signals. Applying the transformation function can include solving directly for current. Solving directly for current can include using finite element analysis. Applying the transformation function can include solving simultaneously for current and potential. Solving simultaneously for current and potential can include using finite element analysis. Solving simultaneously for current and potential can include solving a transformation function which can be expressed as one or more matrices. One or more matrices can be represented by the following equation:
The one or more matrices can be further represented by one or more regularization terms.
Determining the bipolar physiological information at multiple locations of the endocardium surface can include applying a transformation function to the physiological signals, wherein the transformation functions relates the physiological signals to the physiological information at the multiple locations of the endocardium surface. Determining the bipolar physiological information can include determining the transformation function by calculating a forward transformation for relating the physiological information at the multiple locations to the signals measured and inverting the forward transformation. Inverting the forward transformation can include reformulating an underdetermined matrix inversion by regularization. The regularization can be based on a physical relationship. The physical relationship can be a relationship between current and potential. The physical relationship can be represented by the following equation: [Ka,
The one or more electrodes can include multiple spatially distributed electrodes mounted on one or more catheters that are placed inside the heart cavity.
During the measurement of the physiological signals, at least some of the electrodes can be spaced apart from the endocardium surface and/or during the measurement of the physiological signals, at least some of the electrodes can be in contact with the endocardium surface.
The one or more electrodes can include one or more body-surface electrodes. The one or more electrodes can include both electrodes mounted on one or more catheters that are placed inside the body and body-surface electrodes. The one or more electrodes can include electrodes mounted on the one or more catheters that can be moved and positioned at multiple locations in an organ.
The method can also include displaying at least a portion of the determined physiological information. The physiological information can include electrical information.
The method can also include using the determined physiological information to guide treatment of the heart cavity. The treatment can include ablation of one or more selected regions of the heart. The method can also include repeating the measurement of catheter electrode signals and the determination of the physiological information after the ablation treatment. The treatment can include cell therapy, gene therapy, or the application of other biological agents.
The method can also include determining the positions of the electrodes with respect to the endocardium surface. Determining the positions of the electrodes with respect to the endocardium surface can include using at least one of electric fields, magnetic fields, fluoroscopy, and ultrasound to determine a position of the electrodes in a first coordinate system. Determining the positions of the electrodes with respect to the endocardium surface can include registering a representation of the endocardium surface with the first coordinate system. Determining the position of the electrodes can include measuring information about at least one of a position and orientation of the electrodes within the heart cavity. Determining the position of the electrodes can include using a tracking system to track the position of the electrodes. The tracking system can be a system that uses a magnetic field for tracking. The tracking system can be a system that uses injected currents for tracking.
The one or more electrodes can include multiple spatially distributed electrodes.
In some additional aspects, a system includes one or more electrodes configured to measure unipolar signals in response to electrical activity in a heart cavity having a surface and a processing unit configured to determine bipolar physiological information at multiple locations of the surface based on unipolar signals measured by the one or more electrodes and positions of the electrodes with respect to the surface.
Embodiments can include one or more of the following.
The unipolar signals can be unipolar potential signals.
The one or more electrodes can include at least a first electrode and a second electrode. The first electrode can be located within the heart cavity and the second electrode can be located at a distance from the heart cavity such that the electrode is not affected by local tissue activation in a heart cavity. The first and second electrodes can be configured to measure a potential between the first electrode and the second electrode.
The bipolar physiological information can include current density information. The current density information can include a normal component of the current density. The current density information can include a magnitude of a current density vector. The current density information can include a magnitude of a tangential component of the current density.
The processing unit can be further configured to synchronize signals measured at the different catheter positions with one another according to an electrical heart beat cycle. The processing unit can be further configured to synchronize the measured signals based on physiological data including at least one of: ECG and intercardiac electrograms.
The processing unit can be further configured to process the synchronized signals as though they were obtained at one time from all of the positions sampled by the catheter electrodes for the different positions of the catheter in the heart cavity.
The processing unit is further configured to apply a transformation function to the measured unipolar signals. The processing unit can be further configured to solve directly for current. The processing unit can be further configured to solve simultaneously for current and potential. The processing unit can be further configured to solve a transformation function which can be expressed as one or more matrices. The one or more matrices can include one or more regularization terms. The processing unit can be further configured to apply a transformation function to the unipolar signals. The transformation functions can relate the unipolar signals to the bipolar physiological information at the multiple locations of the endocardium surface.
The processing unit can be further configured to determine the transformation function by calculating a forward transformation for relating the physiological information at the multiple locations to the signals measured and inverting the forward transformation. The processing unit can be further configured to reformulate an underdetermined matrix inversion by regularization based on a physical relationship. The processing unit can be further configured to reformulate an underdetermined matrix inversion by a volumetric regularization.
The system can also include a display system configured to display at least a portion of the determined physiological information.
The electrodes can be mounted on one or more catheters that are placed inside the body. The electrodes can include one or more body-surface electrodes. The electrodes can include both electrodes mounted on one or more catheters that are placed inside the body and body-surface electrodes.
The system can also include a tracking system configured to obtain the positions of one or more electrodes. The tracking system can include a system using a magnetic field for tracking. The tracking system can include a system using injected currents for tracking.
The one or more electrodes can be multiple spatially distributed electrodes.
In some additional aspects, a method includes measuring electrical potentials at one or more electrodes in response to electrical activity in a heart cavity having a surface and determining current information at multiple locations of the surface based on the measured electrical potentials and positions of the one or more electrodes with respect to the surface.
Embodiments can include one or more of the following.
The current information can include a normal component of the current density. The current information can include a magnitude of a current density vector. The current information can include a magnitude of a tangential component of the current density.
Measuring the electrical potentials can include measuring the electrical potentials for multiple different catheter positions in a heart cavity having an endocardium surface. The electrical potentials can be measured for at least one electrical heart cycle for each of multiple catheter positions.
The determination of the current information at the multiple locations of the endocardium surface can include synchronizing the signals measured at the different catheter positions with one another according to an electrical heart beat cycle. The determination of the current information at the multiple locations of the endocardium surface further can include processing the synchronized signals as though they were obtained at one time from all of the positions sampled by the catheter electrodes for the different positions of the catheter in the heart cavity. Determining the current information can include applying a transformation function to the measured potential signals. Applying the transformation function can include solving directly for current. Solving directly for current can include using finite element analysis. Applying the transformation function can include solving simultaneously for current and potential. Solving simultaneously for current and potential can include using finite element analysis. Solving simultaneously for current and potential can include solving a transformation function which can be expressed as one or more matrices. The one or more matrices can be represented by the following equation:
The one or more matrices can be further represented by one or more regularization terms. Determining the current information at multiple locations of the endocardium surface can include applying a transformation function to the potential signals. The transformation function can relate the potential signals to the current information at the multiple locations of the endocardium surface. Determining the current information further can include determining the transformation function by calculating a forward transformation for relating the current information at the multiple locations to the potential signals measured and inverting the forward transformation. Inverting the forward transformation can include reformulating an underdetermined matrix inversion by regularization. The regularization can be based on a physical relationship. The physical relationship can be relationship between current and potential. The physical relationship can be represented by the following equation: [Ka,
The regularization can also include limiting the current magnitude. The regularization can be a volumetric regularization.
In some additional aspects, a system includes one or more electrodes configured to measure electrical potentials in response to electrical activity in a heart cavity and a processing unit configured to determine current information at multiple locations of a surface based on unipolar signals measured by the one or more electrodes and positions of the electrodes with respect to the surface.
In some additional aspects, a method can include generating bipolar information based on three or more unipolar measurements.
Embodiments can include one or more of the following.
The unipolar measurements can be unipolar potential signals. The bipolar information can include current information. The bipolar information can current density information. The current density information can include a normal component of the current density. The current density information can include a magnitude of a current density vector. The current density information can include a magnitude of a tangential component of the current density.
The method can also include measuring the unipolar signals by measuring a potential between a first electrode and a second electrode, the first electrode being located within the heart cavity and the second electrode being located at a distance from the heart cavity such that the electrode is not affected by local tissue activation in a heart cavity.
The method can also include measuring the unipolar signals at one or more electrodes for multiple different catheter positions in the heart cavity in response to electrical activity in a heart cavity.
Generating bipolar information based on three or more unipolar measurements can include determining bipolar information at multiple locations of an endocardium surface based on the unipolar signals and positions of one or more electrodes used to measure the unipolar signals with respect to the endocardium surface.
Generating bipolar information can include applying a transformation function to the unipolar signals. Applying the transformation function can include solving directly for current. Applying the transformation function can include solving simultaneously for current and potential.
In some additional aspects, a system includes a processing unit configured to receive unipolar measurements and generate bipolar information based on three or more of the unipolar measurements.
Embodiments of the system may also include devices, software, components, and/or systems to perform any features described above in connection with the first method and/or described below in connection with the second method.
Embodiments of the methods and systems generally disclosed herein can be applied to determining the position of any object within an organ in a patient's body such as the patient's heart, lungs, brain, or liver.
As used herein, the “position” of an object means information about one or more of the 6 degrees of freedom that completely define the location and orientation of a three-dimensional object in a three-dimensional coordinate system. For example, the position of the object can include: three independent values indicative of the coordinates of a point of the object in a Cartesian coordinate system and three independent values indicative of the angles for the orientation of the object about each of the Cartesian axes; or any subset of such values.
As used herein, “heart cavity” means the heart and surrounding tissue.
Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. In case of conflict with documents incorporated herein by reference, the present document controls.
The details of one or more embodiments of the invention are set forth in the accompanying drawings and the description below. Other features, objects, and advantages of the invention will be apparent from the description and drawings, and from the claims.
Like reference symbols in the various drawings indicate like elements.
Embodiments disclosed herein include a method for generating non-contact electro-anatomical maps (“EAM”) with bipolar or current density information. Prior methods for generating non-contact EAMs rely on generating EAMs with unipolar information.
Embodiments disclosed herein further describe a method for generating contact EAM with bipolar or current density information in locations that were not probed with a contact electrode.
In both methods (e.g., the non-contact method and the contact method), further post-processing operations can be performed on the reconstructed physiological information (e.g., the bipolar or current density information) to extract and display useful features of the bipolar or current density information.
In some embodiments, as described in more detail below, Laplace's Equation can be used for deriving the physiological information, sometimes using approximations or discretization methods. In addition, a volumetric regularization scheme that regularizes electrode measurements in addition to surface potentials can be used (e.g., as opposed to the use of surface regularization alone).
While the examples described herein focus on intracardiac endocardial mapping, the advantages of generating a bipolar non-contact map also hold for non-invasive body surface mapping which generate a bipolar epicardial EAM.
In general, when generating an EAM the physician is interested in relating electrical properties of tissue to anatomy. Two types of electrical measurements are typically performed to investigate electrical properties of tissue. The first is unipolar or potential measurement, Vu. In this case a voltage is measured between an investigated electrode and an “indifferent” electrode assumed to be far enough away so that it is not affected by local tissue activation (e.g., at least about 5 cm away from the investigated electrode, at least about 10 cm away from the investigated electrode, at least about 15 cm away from the investigated electrode, at least about 20 cm away from the investigated electrode, at least about 25 cm away from the investigated electrode, etc.). The “indifferent” electrode can be positioned inside the body, or derived from a surface potential such as Wilson's Central Terminal. The second type of measurement is bipolar measurement, Vb. In this case a voltage is measured between two close electrodes, usually a few millimeters (e.g., less than 3-5 mm) apart. In the bipolar case, both electrodes are in the heart and affected by local tissue activation. The bipolar measurement can be viewed as an estimation of the current density in the direction of the bipolar electrode pair.
Referring to
Both unipolar and bipolar measurements can be used in catheter ablation procedures. Bipolar measurements are believed to have a number of advantages over unipolar measurements. See, for example, Fred M. Kusumoto, “Unipolar Recording in Cardiac Electrophysiologic studies”, Journal of Interventional Cardiac Electrophysiology 1999; 3.
One advantage of bipolar measurements is far field rejection. A goal of electro-anatomical mapping is to characterize the properties of underlying tissue in the vicinity of the measuring electrode. The far field is a field generated by tissue that is further from the point of interest, but is still large enough to affect measurement. Electrodes measure potential which results from a summation of all fields generated by tissue. In the unipolar case, it is estimated that the amplitude of measured signal is proportional to the reciprocal value of the square of the distance between tissue and the electrode. In the bipoloar case it is estimated that the effect of distance is cubed. Therefore, bipolar measurements are better at measuring local tissue excitation and rejecting far field. Far field rejection can be of particular importance when mapping in the atria and when mapping in infracted tissue.
Another advantage of bipolar measurements is rejection of interference and noise. In the bipolar case, because the electrodes are spaced closely and their signal carrying wires are run down the catheter in close proximity, the effect of noise on the two electrodes has a large common component. Since electrode values are subtracted, much of the noise common to both electrodes can be rejected. Due to the larger distance between electrodes in the uniploar case, a smaller portion of the noise collected by electrodes is common, and less of it may be rejected. This leads to better signal to noise performance in the bipolar case.
Another advantage of bipolar measurements is the experience level and accumulated knowledge about such measurements. Due to some of the advantages of bipolar measurements much experience and knowledge has been accumulated with bipolar measurements. For example, algorithms and methodologies that depend on specific bipolar threshold and signal behavior have been developed. To benefit from this body of knowledge and user experience it is important to provide bipolar EAM in a non-contact mapping system.
Methods for generating estimated unipolar measurements using non-contact mapping are described, for example, in U.S. Pat. No. 7,515,954 entitled “Non-contact cardiac mapping, including moving catheter and multi-beat integration” and filed Jun. 13, 2006 which is hereby incorporated in its entirety by reference. In contrast to generating estimated unipolar measurements, embodiments disclosed herein include systems and methods of directly and indirectly computing an estimation of bipolar measurement based on measured unipolar potential values (e.g., an estimation of a current density value). It is believed that one advantage of computing an estimation of bipolar measurement based on measured unipolar potential values is making available the advantages of bipolar mapping in a non-contact methodology. In addition, due to the nature of non-contact mapping and algorithms presented herein, it is also possible to provide current density maps, which can be viewed as mathematically ideal bipolar maps.
Embodiments disclosed herein include methods of performing regularization. In prior systems for generating EAM maps, due to the underdetermined ill-posed nature of the inverse problem surface regularization is applied such that a-priori assumptions regarding the distribution of the potential on the heart surface are used to constrain the estimated surface potential, for example, as described in U.S. Pat. No. 7,515,954 entitled “Non-contact cardiac mapping, including moving catheter and multi-beat integration” and filed Jun. 13, 2006 which is hereby incorporated in its entirety by reference. These a-priori assumptions are not governed by physical or biological laws, and therefore, add error to the solution in cases where actual surface distribution varies from the a-priori assumptions. Examples of such surface regularization schemes that employ these a-priory assumptions include Tikhonov 0 and 1 which limits signal or gradient magnitude. Embodiments disclosed herein include a regularization method which takes advantage of the physical relationship between current and potential and regularizes the relationship between those two estimated values. Since this method relies on the physical relationship between current and potential it is believed to greatly improve EAM accuracy.
It should be appreciated that the term non-contact as used herein is considered for any measurement that is not required to be on the surface of the cardiac chamber (e.g., measurements where the catheter is spaced apart from the endocardium surface and/or measurements where the catheter touches the endocardium surface but such contact is not required for subsequent calculations based on the measurement). The same governing equations and calculations can be applied for measurements that are taken at the surface of the chamber in order to determine information at locations on the surface that were not directly probed. In other words, it is possible to collect multiple contact unipolar or potential measurements and use the same formulation to derive bipolar or current density EAM of the entire chamber. Similarly, a combination of contact and non-contact measurements can be used and the method disclosed would still apply.
The relationship between bipolar measurement and current density is described below. In addition, two methods for providing an estimation of bipolar measurement are described herein. The first method uses a unipolar potential based inverse engine to estimate bipolar signal. The second method provides a means to compute bipolar measurement or current density directly from measured electrode potentials. Each of these methods is described in more detail below following a description of the relationship between bipolar measurements and current density.
Current Density and Bipolar Signal
As the distance between electrodes and their physical size diminish in a bipolar catheter, the bipolar measurement becomes proportional the electric field vector (which in turn is proportional to the current density vector) in the given catheter orientation.
As Se1, Se2 and de approach zero, the bipolar measurement is directly proportional to the current density. In a typical catheter, these dimensions are a few millimeters, such that a bipolar catheter can be viewed as estimating an electric field, which in a homogenous conductor such as blood is directly proportional to the current vector. It therefore follows that a current density map can be viewed as an ideal bipolar map. The map is ideal in that its value does not depend on electrode size and spacing which create inconsistencies across real physical bipolar catheters.
Indirect Solution for Bipolar Measurement
One approach to obtaining bipolar measurements using a non-contact system is to subtract unipolar data from neighboring locations. Referring to
We note that the orientation of the bipolar electrodes relative to the heart will affect the bipolar measurement amplitude. For example, a reversed orientation (180°) will provide the negative value of the original (e.g., in
In order to obtain an estimated bipolar signal from a unipolar potential inverse engine, we generate an inverse solution for surface potentials {circumflex over (V)}e. As shown in
As shown in
As shown in
While the above provides a method of expanding potential inversion to provide a bipolar measurement, there is a reduction in accuracy due to regularization. When surface Tikhonov regularization is used in potential inversion, the estimation of potential distribution on the surface is smoothed by the regularization operator. The smoothing operator reduces the difference in measurement between neighboring locations, which in turn reduces and adds error to bipolar amplitude obtained in the above scheme. Therefore, in some embodiments, it can be advantageous to use a direct solver for the bipolar signal as described below.
Forward Solution Formulations
In the example provided below, the direct solution for current density is demonstrated using the finite element method. It should be appreciated that a number of other numerical methods such as boundary element, finite volume, finite difference, etc. could be used to accomplish the same goal. In the example below, after introducing some governing equations, a forward finite element formulation for the unipolar case is described, followed by a more generalized forward formulation for the potential and current density case. Additionally, and inverse operation which adds a volumetric regularization is also described below.
Governing Equations for Endocardial and Epicardial Problems
A stationary conduction forward problem is formulated by three Maxwell's equations:
∇×E=0 (0.1)
∇·J=0 (0.2)
J=σ(x)E (0.3)
Equation (0.1) expresses that the electric field E is conservative. Thus, E can be described by a scalar electric potential Φ as follows:
E=−∇Φ (0.4)
Equation (0.2) states that the current density vector J is source free in the intracavitary domain and (0.3) is the differential Ohm's law. In the following description, it is assumed that the electric conductivity σ(x) is spatially varying, thereby treating both the endocardial and epicardial problems simultaneously. The endocardial solution to this problem is useful for intracardiac mapping which can be used in minimally invasive procedures such as catheter ablation. The epicardial problem is useful for non-invasive body surface based EAM which can be used to diagnose myocardial infarction and other heart disease. For the sake of simplicity, as used herein, σ denotes general spatial variation.
−∇·(σ∇Φ)=0 in Ω (0.5)
By prescribing the electric scalar unipolar potential on the endocardium,
Φ=Ve on Γe (endocardium) (0.6)
a boundary value problem (BVP) is obtained, the so-called Dirichlet problem, which has a well-posed unique solution.
The other problem is the epicardial problem of electrocardiography, which is depicted in
where Γe and Ve now denote the epicardial surface and the epicardial unipolar potential, respectively. Equation (0.8) expresses the fact that no current flows through the torso surface Γt and it is called the homogeneous Neumann boundary condition. It is well established that the BVP (0.5), (0.7) and (0.8) also has a unique and stable solution.
As shown below in the following two subsections, the boundary condition (0.8) does not explicitly appear in the finite element formulation (so comes the alternative name natural boundary condition). Therefore, the solution of endocardial problems with spatially varying electric conductivity, can be applied without modification to epicardial body surface problems as well.
On the endocardium or epicardium—depending on the type of the problem—different types of EAM's can be generated. One way is to display the electric potential, Ve. As explained above, it can be advantageous to also compute current density based EAM's which are proportional to the bipolar measurement. Three options of expressing the current density are described below. One is to display the magnitude of the current density vector
J=|J|=√{square root over (Jx2+Jy2+Jz2)} (0.9)
A second option is to consider only the normal component of the current density
Jn=J·n (0.10)
where n is the unit normal vector of the endocardial or epicardial surfaces. A third option is to represent the magnitude of the tangential component of the current density
Jt=|J−J·n| (0.11)
The total and tangential magnitudes, J and Jt, loose “sign” information while Jn retains it.
Electric Scalar Potential (Unipolar) Formulation
Weak Form of the Forward Problems
In order to derive the solution algorithm of the Dirichlet problem by the finite element method, the solution below begins by introducing the weak forms of the endocardial and epicardial BVP's (0.5), (0.6) and (0.5), (0.7), (0.8), respectively. Beginning with the endocardial problem, multiply equation (0.5) by an arbitrary admissible weighting function w with vanishing value on the endocardium and integrate the product over the whole domain Ω
By applying the vector identity
−w∇·a=∇w·a−∇·(wa) (0.13)
with a=−σ∇Φ and the Gauss's law, the following equation is obtained
where the second term is zero due to the vanishing weighting functions on the endocardial surface.
The weak form for the Laplace epicardial problem (0.5), (0.7), (0.8) can be derived similarly
where the first surface term is zero due to the vanishing weighting functions on the epicardial surface similarly to the endocardial case in (0.14). The natural boundary condition (0.8) can be “weakly” enforced by setting the second surface term in (0.15) to zero. We can conclude that the weak forms of the endocardial and epicardial problems (0.14), (0.15) are in the same form if spatially varying electric conductivity were assumed. Thus, the examples below discuss the finite element discretization and the inverse solutions of the endocardial problem only and the results can be applied to the epicardial problem straightforwardly by simply substituting endocardial surface with epicardial surface.
In order to handle non-zero Dirichlet boundary conditions, the electric vector potential can be split into two parts
Φ=Φe+ΦaΦe=Ve, Φa=0 on Γe (0.16)
where the only requirement for the volumetric function Φe is to have the same value as the prescribed surface potential Ve on the endocardium but arbitrary otherwise. By substituting (0.16) into (0.14) the following equation is obtained:
Equation (0.17) is called the weak form of the Dirichlet problem (0.5), (0.6). This means that if (0.17) is satisfied for any admissible weighting function w then Φ=Φe+Φa is the solution of the BVP (0.5), (0.6).
Finite Element Discretization of Laplace Equation
Equation (0.17) provides a starting point for the finite element formulation. By interpolating Φ:
where αn is a linear nodal interpolation function associated with vertex n as illustrated in
αn=λn (0.19)
The split (0.16) can easily be realized by the finite element interpolation scheme
where e and
a are the sets of endocardial and intracavitary vertex numbers, respectively. By substituting (0.20) into (0.17) and using the intracavitary finite element interpolation functions as weighting functions, the following equation is obtained:
By changing the order of the integration and summation the following equation is generated:
By choosing the intracavitary interpolation functions as weighting functions in (0.21) the constraint for the weighting functions in (0.14) is satisfied automatically. Equation (0.22) is a symmetric algebraic equation system and it can be written in a compact form
[Kaa][Φa]=−[Kae][Φe] (0.23)
where
and [Φa] and [Φe] are column vectors with the intracavitary and endocardial nodal potentials. With such nodal finite element representation, [Φe]=[Ve] thus (0.23) becomes
[Kaa][Φa]=−[Kae][Ve] (0.26)
Equations system (0.26) is a positive definite linear system for the intracavitary nodal potential values and can be solved uniquely.
Direct Solution for Current Density—Scalar Potential Mixed Formulation
Governing Equations for the Mixed Formulation
There are several alternative ways to formulate the electric conduction problem (0.1)-(0.3). Our goal is to use a formulation where the current density J is explicitly computed together with the scalar potential. In order to achieve that, we write the conduction problem (0.1)-(0.3) in a “mixed” form that includes two equations
Equations (0.27), (0.28) with the boundary condition (0.6) is also a well-posed BVP with a unique solution. The advantage of this formulation can be that it provides both the current density and the electric potential without the need of any further differentiation.
Weak Form of the Mixed Formulation
To obtain the finite element discretization, we have to first obtain the weak form of the system (0.27), (0.28). Let us multiply (0.27) with an arbitrary vector weighting function W and integrate over Ω
where L2(Ω) is the function space of square integrable functions in Ω. If we utilize the vector identity (0.13) again, we arrive at
Similarly, multiply (0.28) by a scalar weighting function w and integrate over the domain
If equations (0.30), (0.31) are satisfied for all admissible vector and scalar weighting functions W and w, J and Φ are the solutions of the BVP (0.27), (0.28). This means that the weak forms (0.30), (0.31) are equivalent to the strong problem (0.27), (0.28), respectively. Note that the same formulation can be used for the epicardial problem.
Finite Element Discretization
In order to obtain the finite element equations, we proceed similarly to the scalar potential formulation. We interpolate the scalar potential by linear nodal finite elements as in (0.20). In order to obtain a unique solution, we have to use linear facet finite element functions
where
βf=βmno=2(λm∇λn×∇λo+λn∇λo×∇λm+λo∇λm×∇λn)
and Jf are unknown facet current density coefficients. Note that the interpolation form (0.33) ensures the desired property of weighting functions in (0.29). The facet interpolation functions are illustrated in
Let us plug (0.32), (0.20) into (0.30), (0.31), and use the vector and nodal interpolation functions as weighting functions. After changing the order of the summations and integrations, we obtain
Equations (0.34), (0.35) can be written in compact form
Equation system (0.36) has a unique solution and it provides the current density and scalar potential simultaneously.
Inverse Algorithm for Current Density and Electric Scalar Potential Simultaneously
Inversion in Terms of Electric Scalar Potential
In order to illustrate how the inverse conduction problems of electrocardiography can be formulated when finite element discretization is utilized, the electric scalar potential formulation (0.22) provides a starting point. As described above, only the endocardial inverse problem shown in c and the rest of the volume vertex numbers by
a=
c∪
The inverse finite element equation can be formed by rearranging the equation system (0.22) to have the catheter vertices on the right hand side
As above, we can write this equation in a compact form:
[Ka,
Equation (2.2) is severely ill-posed and introduction of regularization is needed to solve it.
We have stated previously that all the results for the endocardial forward problem are valid for the epicardial problem without modification. This is also true for the inverse discussion of the epicardial inverse problems. Then c denotes the set of vertex numbers on the torso surface where measurements are available and
Inversion by Mixed Formulation
In this subsection, a formulation of the solution of the endocardial inverse problem in terms of the mixed formulation (0.27), (0.28) is shown. One advantage of the mixed formulation is that it provides both the endocardial current density and scalar unipolar potential simultaneously. Furthermore, this formulation couples the two physical quantities as part of its regularization scheme.
Similarly to the rearrangement in the scalar potential case, (0.34), (0.35) can be rearranged as
These equations can also be written in a compact form as:
The difficulty is that equation (2.5) is ill-posed and needs regularization to obtain a solution.
It can be beneficial to utilize the unipolar potential formulation 2.2) as regularization thereby coupling the estimated current to scalar unipolar potential. This is substituted into the zero block of the system matrix (2.5) as a regularization matrix multiplied by a regularization parameter
An additional regularization for the current density J has to be added to limit the current density magnitude
where the parameter μ is another regularization parameter. The regularization matrix [MTT] is also a volumetric regularization matrix. One exemplary form can be represented as:
where SK=[volume(ΩK)]2/3 is a weighting factor and ΩK denotes the tetrahedrons in which the integration takes place. The solution of (2.7) is well-posed and unique for J and Φ if we use reasonable values for the regularization parameters μ and λ. It is believed that, values between 0.4 and 1.0 for both regularization parameters provide increased inversion quality.
Post Processing
Various types of post-processing can be used. In some embodiments, the post processing may involve selecting a format for outputting (e.g., displaying) the reconstructed bipolar or current density information to a user. In other embodiments, the post-processing may involve significant further mathematical manipulation of the reconstructed potentials to provide additional types of physiological information.
Some of the post-processing operations performed on the reconstructed set(s) of physiological information include the generation of a resolution map. Such a resolution map indicates the spatial resolution of reconstructed physiological information at points on the endocardium surface, thereby providing a measure of the reliability and accuracy of the information at various points on the endocardium surface. Resolution maps may be used with any form of post-processing operation including all modes listed below. The resolution map may be superimposed with physiological information to highlight those areas that are accurate and those that are not. Strictly speaking, information about the resolution maps can be determined prior to obtaining the reconstructed information; however, herein we generally refer to the generation and display of the resolution map as “post-processing” because such information is typically presented to the user with reconstructed physiological information.
Another type of post-processing operation that may be performed includes the generation of isocurrent maps. Particularly, where the reconstructed physiological information pertains to bipolar or current density, the magnitude of reconstructed information may be color coded and superimposed on the 3D endocardial representation. Isocurrent maps are the reconstructed current densities computed for every sampled set of data over a single or multiple heart beats. For the magnitude either the normal, tangential or total bipolar or current density may be computed.
Another type of post-processing operation that may be performed includes the generation of vector isocurrent maps. Particularly, where the reconstructed physiological information pertains to bipolar or current density, the reconstructed information may be represented by arrows on vertices in the 3D endocardial representation. Vector isocurrent maps are the reconstructed current densities computed for every sampled set of data over a single or multiple heart beats. In this representation the length of the arrow represents magnitude while its direction is the direction of the current or bipolar vector. Yet another type of post-processing operation includes the generation of timing maps (such as activation time maps). The timing maps provide information on the time-dependent behavior of the heart's electrical activity. Particularly, the activation map indicates at what point in time particular points on the endocardium surface experience a change in their electrical activity. For example, the activation map could identify the point in time at which particular cells on the endocardium surface experienced depolarization. Another type of timing map may be an iso-duration map where the amount of time certain tissue has been active for is detected. Timing maps may be computed from the reconstructed bipolar or current density information over a single or multiple heart beats. Timing maps may be determined and displayed for one or more points on the endocardium surface representation.
Another type of post processing operation is the generation of amplitude maps. Amplitude maps can be used to display characteristics of bipolar or current amplitude in a given area. The amplitude maps may be computed from the reconstructed bipolar or current density information over a single or multiple heart beats. Useful amplitude map information that may be determined and displayed for one or more points on the endocardium surface representation includes the maximum amplitude, or root mean square values.
Another type of post-processing operation is the generation of a difference map. The difference map provides information regarding the effectiveness of the clinical procedure (e.g., ablation) performed on the patient to ameliorate the symptoms of arrhythmias. The difference map compares the electrical behavior of the heart, as reflected from two or more voltage maps generated before and after the performance of the particular clinical procedure.
A further type of post processing operation is the generation of frequency maps. Frequency mapping, and more generally spectral analysis, are used to identify on the endocardium surface localized sites of high-frequency activity during fibrillation. Frequency maps are computed by acquiring multiple sets of reconstructed information over a particular time interval which includes a single or multiple heart beats. The acquired raw data is then used to obtain the frequency representation of that data. Specific information (e.g., dominant frequency components) from the frequency representation is subsequently identified, and that identified information may be displayed.
Other types of post-processing information may likewise be performed.
In order to generate the exemplary physiological data shown in
The dataset that was collected was used to generate several maps for comparing different computational methodologies including generating of a unipolar measurement, an indirect solution for bipolar measurement, and a direct solution for bipolar measurement.
Representative System
In some embodiments the distal end of the catheter 110 is fitted with multiple electrodes spread somewhat uniformly over the catheter. For example, the electrodes may be mounted on the catheter 110 following a 3D olive shape. The electrodes are mounted on a device capable of deploying the electrodes into the desired shape while inside the heart, and retracting the electrodes when the catheter is removed from the heart. To allow deployment into a 3D shape in the heart, electrodes may be mounted on a balloon, or shape memory material such as Nitinol.
At each of the locations to which the catheter 110 is moved, the catheter's multiple electrodes acquire signals resulting from the electrical activity in the heart cavity. Consequently, reconstructing and presenting to a user (such as a doctor and/or technician) physiological data pertaining to the heart's electrical activity may be based on information acquired at multiple locations, thereby providing a more accurate and faithful reconstruction of physiological behavior of the endocardium surface. The acquisition of signals at multiple catheter locations in the heart chamber enables the catheter to effectively act as a “mega-catheter” whose effective number of electrodes and electrode span is proportional to the product of the number of locations in which signal acquisition is performed and the number of electrodes the catheter has.
To enhance the quality of the reconstructed physiological information at the endocardium surface, in some embodiments the catheter 110 is moved to more than three locations (for example, more than 5, 10, or even 50 locations) within the heart chamber. Further, the spatial range over which the catheter is moved may be larger than one third (⅓) of the diameter of the heart cavity (for example, larger than 35%, 40%, 50% or even 60% of the diameter of the heart cavity). Additionally, in some embodiments the reconstructed physiological information is computed based on signals measured over several heart beats, either at a single catheter location within the heart chamber or over several locations. In circumstances where the reconstructed physiological information is based on multiple measurements over several heart beats, the measurements are synchronized with one another so that the measurement are performed at approximately the same phase of the heart cycle. The signal measurements over multiple beats can be synchronized based on features detected from physiological data such as surface ECG or intracardiac electrograms.
Mapping system 100 further includes the processing unit 120 which performs several of the operations pertaining to the mapping procedure, including the reconstruction procedure to determine the physiological information at the surface (e.g., as described above). To expedite the computational operations performed by the mapping system 100, the processing unit 120 can compute, generally prior to the insertion of the catheter into the heart chamber and/or before signal acquisition by the catheter's electrodes has commenced, transformation functions that can be used in real-time to facilitate the reconstruction process. Once the catheter 110 is inserted and is displaced to a particular location in the heart chamber, the mapping procedure can be performed expeditiously by computing in real-time those transformation components that were not computed ahead of the signal acquisition stage, and combining those components with the appropriate pre-processed transformation components to obtain the overall transformation function(s). That overall transformation function is applied to the acquired raw data to perform the inverse reconstruction operation.
The processing unit 120 also performs a catheter registration procedure. The location of the catheter 110 inserted into the heart chamber can be determined using a conventional sensing and tracking system (not shown) that provide the 3D spatial coordinates of the catheter and/or its multiple electrodes with respect to the catheter's coordinate system as established by the sensing and tracking system. However, to perform the mapping procedure and reconstruct physiological information on the endocardium surface, it is necessary to align the coordinate system of the catheter 110 with the endocardium surface's coordinate system. The processing unit 120 (or some other processing module of system 100) determines a coordinate system transformation function that transforms the 3D spatial coordinates of the catheter's locations into coordinates expressed in terms of the endocardium surface's coordinate system, or vice-versa.
The processing unit 120 also performs post-processing operations on the reconstructed physiological information to extract and display useful features of the information to the operator of the system 100 and/or other persons (e.g., a physician).
As further shown in
As further shown in
The methods and systems described herein are not limited to a particular hardware or software configuration, and may find applicability in many computing or processing environments. The methods and systems can be implemented in hardware, or a combination of hardware and software, and/or can be implemented from commercially available modules applications and devices. Where the implementation of the systems and methods described herein is at least partly based on use of microprocessors, the methods and systems can be implemented in one or more computer programs, where a computer program can be understood to include one or more processor executable instructions. The computer program(s) can execute on one or more programmable processors, and can be stored on one or more storage medium readable by the processor (including volatile and non-volatile memory and/or storage elements), one or more input devices, and/or one or more output devices. The processor thus can access one or more input devices to obtain input data, and can access one or more output devices to communicate output data. The input and/or output devices can include one or more of the following: Random Access Memory (RAM), Redundant Array of Independent Disks (RAID), floppy drive, CD, DVD, magnetic disk, internal hard drive, external hard drive, memory stick, or other storage device capable of being accessed by a processor as provided herein, where such aforementioned examples are not exhaustive, and are for illustration and not limitation.
The computer program(s) can be implemented using one or more high level procedural or object-oriented programming languages to communicate with a computer system; however, the program(s) can be implemented in assembly or machine language, if desired. The language can be compiled or interpreted. The device(s) or computer systems that integrate with the processor(s) can include, for example, a personal computer(s), workstation (e.g., Sun, HP), personal digital assistant (PDA), handheld device such as cellular telephone, laptop, handheld, or another device capable of being integrated with a processor(s) that can operate as provided herein. Accordingly, the devices provided herein are not exhaustive and are provided for illustration and not limitation.
References to “a microprocessor” and “a processor”, or “the microprocessor” and “the processor,” can be understood to include one or more microprocessors that can communicate in a stand-alone and/or a distributed environment(s), and can thus be configured to communicate via wired or wireless communications with other processors, where such one or more processor can be configured to operate on one or more processor-controlled devices that can be similar or different devices. Furthermore, references to memory, unless otherwise specified, can include one or more processor-readable and accessible memory elements and/or components that can be internal to the processor-controlled device, external to the processor-controlled device, and can be accessed via a wired or wireless network using a variety of communications protocols, and unless otherwise specified, can be arranged to include a combination of external and internal memory devices, where such memory can be contiguous and/or partitioned based on the application. Accordingly, references to a database can be understood to include one or more memory associations, where such references can include commercially available database products (e.g., SQL, Informix, Oracle) and also proprietary databases, and may also include other structures for associating memory such as links, queues, graphs, trees, with such structures provided for illustration and not limitation.
Accordingly, other embodiments are within the scope of the following claims.
Number | Name | Date | Kind |
---|---|---|---|
4649924 | Taccardi | Mar 1987 | A |
4674518 | Salo | Jun 1987 | A |
4840182 | Carlson | Jun 1989 | A |
4920490 | Isaacson | Apr 1990 | A |
5156151 | Imran | Oct 1992 | A |
5284142 | Goble et al. | Feb 1994 | A |
5297549 | Beatty et al. | Mar 1994 | A |
5300068 | Rosar et al. | Apr 1994 | A |
5341807 | Nardella | Aug 1994 | A |
5381333 | Isaacson et al. | Jan 1995 | A |
5469858 | Osborne | Nov 1995 | A |
5480422 | Ben-haim | Jan 1996 | A |
5500011 | Desai | Mar 1996 | A |
5553611 | Budd et al. | Sep 1996 | A |
5568809 | Ben-Haim | Oct 1996 | A |
5577502 | Darrow et al. | Nov 1996 | A |
5588429 | Isaacson et al. | Dec 1996 | A |
5634469 | Bruder et al. | Jun 1997 | A |
5662108 | Budd et al. | Sep 1997 | A |
5687737 | Branham et al. | Nov 1997 | A |
5697377 | Wittkampf | Dec 1997 | A |
5704365 | Albrecht et al. | Jan 1998 | A |
5722402 | Swanson et al. | Mar 1998 | A |
5840025 | Ben-Haim | Nov 1998 | A |
5840031 | Crowley | Nov 1998 | A |
5846198 | Killmann | Dec 1998 | A |
5848972 | Triedman et al. | Dec 1998 | A |
5921982 | Lesh et al. | Jul 1999 | A |
5954665 | Ben-Haim | Sep 1999 | A |
5971933 | Gopakumaran et al. | Oct 1999 | A |
5983126 | Wittkampf | Nov 1999 | A |
6014581 | Whayne et al. | Jan 2000 | A |
6050267 | Nardella et al. | Apr 2000 | A |
6095150 | Panescu et al. | Aug 2000 | A |
6167296 | Shahidi | Dec 2000 | A |
6226542 | Reisfeld | May 2001 | B1 |
6236886 | Cherepenin et al. | May 2001 | B1 |
6240307 | Beatty et al. | May 2001 | B1 |
6246898 | Vesely et al. | Jun 2001 | B1 |
6254536 | DeVito | Jul 2001 | B1 |
6278894 | Salo et al. | Aug 2001 | B1 |
6298257 | Hall et al. | Oct 2001 | B1 |
6308093 | Armoundas et al. | Oct 2001 | B1 |
6314310 | Ben-Haim et al. | Nov 2001 | B1 |
6317619 | Boernert et al. | Nov 2001 | B1 |
6318375 | Plicchi et al. | Nov 2001 | B1 |
6360123 | Kimchi et al. | Mar 2002 | B1 |
6368285 | Osadchy et al. | Apr 2002 | B1 |
6400981 | Govari | Jun 2002 | B1 |
6516807 | Panescu et al. | Feb 2003 | B1 |
6547082 | Babini | Apr 2003 | B1 |
6556695 | Packer et al. | Apr 2003 | B1 |
6574492 | Ben-Haim et al. | Jun 2003 | B1 |
6574498 | Gilboa | Jun 2003 | B1 |
6593884 | Gilboa et al. | Jul 2003 | B1 |
6600948 | Ben-Haim et al. | Jul 2003 | B2 |
6603996 | Beatty et al. | Aug 2003 | B1 |
6631290 | Guck et al. | Oct 2003 | B1 |
6640119 | Budd et al. | Oct 2003 | B1 |
6650927 | Keidar | Nov 2003 | B1 |
6690963 | Ben-Haim et al. | Feb 2004 | B2 |
6701176 | Halperin et al. | Mar 2004 | B1 |
6728562 | Budd et al. | Apr 2004 | B1 |
6773402 | Govari et al. | Aug 2004 | B2 |
6807439 | Edwards et al. | Oct 2004 | B2 |
6839588 | Rudy | Jan 2005 | B1 |
6847839 | Ciaccio et al. | Jan 2005 | B2 |
6872428 | Yang et al. | Mar 2005 | B2 |
6892090 | Verard et al. | May 2005 | B2 |
6892091 | Ben-Haim et al. | May 2005 | B1 |
6939309 | Beatty et al. | Sep 2005 | B1 |
6957101 | Porath et al. | Oct 2005 | B2 |
6978168 | Beatty et al. | Dec 2005 | B2 |
6990370 | Beatty et al. | Jan 2006 | B1 |
7016719 | Rudy et al. | Mar 2006 | B2 |
7043292 | Tarjan et al. | May 2006 | B2 |
7263397 | Hauck et al. | Aug 2007 | B2 |
7433730 | Berrier | Oct 2008 | B1 |
7505810 | Harlev et al. | Mar 2009 | B2 |
7515954 | Harlev et al. | Apr 2009 | B2 |
7729752 | Harlev et al. | Jun 2010 | B2 |
8137343 | Harlev et al. | Mar 2012 | B2 |
20020151807 | Goldin | Oct 2002 | A1 |
20030018251 | Solomon | Jan 2003 | A1 |
20030065271 | Khoury | Apr 2003 | A1 |
20030076277 | Muramatsu et al. | Apr 2003 | A1 |
20030078509 | Panescu | Apr 2003 | A1 |
20030216630 | Jersey-Willuhn et al. | Nov 2003 | A1 |
20040077942 | Hall et al. | Apr 2004 | A1 |
20040097806 | Hunter et al. | May 2004 | A1 |
20040243015 | Smith et al. | Dec 2004 | A1 |
20040254437 | Hauck et al. | Dec 2004 | A1 |
20050038337 | Edwards | Feb 2005 | A1 |
20050054918 | Sra | Mar 2005 | A1 |
20050107834 | Freeman et al. | May 2005 | A1 |
20050154282 | Li et al. | Jul 2005 | A1 |
20050288599 | MacAdam et al. | Dec 2005 | A1 |
20060085049 | Cory et al. | Apr 2006 | A1 |
20060116575 | Willis | Jun 2006 | A1 |
20060122526 | Berenfeld et al. | Jun 2006 | A1 |
20060173251 | Govari et al. | Aug 2006 | A1 |
20060178587 | Khoury | Aug 2006 | A1 |
20060241401 | Govari et al. | Oct 2006 | A1 |
20070016007 | Govari et al. | Jan 2007 | A1 |
20070038078 | Osadchy | Feb 2007 | A1 |
20070049821 | Willis | Mar 2007 | A1 |
20070197929 | Porath et al. | Aug 2007 | A1 |
20070265539 | Hastings et al. | Nov 2007 | A1 |
20070270703 | He et al. | Nov 2007 | A1 |
20070287902 | Fuimaono et al. | Dec 2007 | A1 |
20070299351 | Harlev et al. | Dec 2007 | A1 |
20070299352 | Harlev et al. | Dec 2007 | A1 |
20070299353 | Harlev et al. | Dec 2007 | A1 |
20080190438 | Harlev et al. | Aug 2008 | A1 |
20080221566 | Krishnan | Sep 2008 | A1 |
20080234588 | Feldman et al. | Sep 2008 | A1 |
20080249424 | Harlev et al. | Oct 2008 | A1 |
20090171274 | Harlev et al. | Jul 2009 | A1 |
20090177072 | Harlev et al. | Jul 2009 | A1 |
20090253976 | Harlev et al. | Oct 2009 | A1 |
20100286551 | Harlev et al. | Nov 2010 | A1 |
Number | Date | Country |
---|---|---|
2002-51998 | Feb 2002 | JP |
2007-185516 | Jul 2007 | JP |
5703288 | Apr 2015 | JP |
Entry |
---|
Adams et al., “Seeded Region Growing”, IEEE Transactions on Pattern Analysis and Machine Intelligence, 16(6):641-647, 1994. |
Baan, Jan et al., “Continuous Measurement Of Left Ventricular Volume In Animals And Humans By Conductance Catheter”, Circulation, vol. 70,pp. 812-823. |
Ben-Haim et al., “Nonfluoroscopic, in Vivo Navigation and Mapping Technology”, Nature Medicine, 2(12):1393-1395, 1996. |
Besl et al., “A Method for Registration of 3-D Shapes”, IEEE Transactions on Pattern Analysis and Machine Intelligence, 14(2):239-256, 1992. |
Blomström-Lundqvist et al., “ACC/AHA/ESC Guidelines for the Management of Patients With Supraventricular Arrhythmias-Executive Summary”, Journal of the American College of Cardiology, 42(8):1493-1531, 2003. |
Breithardt et al., “AHA Medical/Scientific Statement—Special Report: Standards for Analysis of Ventricular Late Potentials Using High-Resolution or Signal-Averaged Electrocardiography”, Circulation, 83(4):1481-1488, 1991. |
Brooks et al., “Electrical Imaging of the Heart”, IEEE Signal Processing Magazine, pp. 24-42, 1997. |
Caspi et al., “Stem Cell Research: Regenerating the Heart Using Human Embryonic Stem Cells—from Cell to Bedside”, IMAJ, 8:208-214, 2006. |
De Groot et al., “Three-Dimensional Catheter Positioning During Radiofrequency Ablation in Patients: First Application of a Real-Time Position Management System”, Journal of Cardiovascular Electrophysiology, 11:1183-1192, 2000. |
Donahue et al., “Focal Modification of Electrical Conduction in the Heart by Viral Gene Transfer”, Nature Medicine, 6(12):1395-1398, 2000. |
Dong et al., “Integrated Electroanatomic Mapping With Three-Dimensional Computed Tomographic Images for Real-Time Guided Ablations”, Circulation, 113:186-194, 2006. |
Durrer et al., “Total Excitation of the Isolated Human Heart”, Circulation, vol. XLI, pp. 899-912, 1970. |
Ector, Joris et al., “Cardiac Three-Dimensional Magnetic Resonance Imaging And Fluoroscopy Merging”, Circulation, (Dec. 13, 2005). |
Friedman, “Catheter Cryoablation of Cardiac Arrhythmias”, Current Opinion in Cardiology, 20:48-54, 2005. |
Friedman, “Novel Mapping Techniques for Cardiac Electrophysiology”, Heart, 87:575-582, 2002. |
Geddes, L.A. et al., “Criteria for the Selection of Materials for Implanted Electrodes”, Annals of Biomedical Engineering, vol. 31, pp. 879-890 (2003). |
Gepstein et al., “A Novel Method for Nonfluoroscopic Catheter-Based Electroanatomical Mapping of the Heart”, Circulation 95:1611-1622, 1997. |
Huang, Yi-Chih et al., “Development Of A Third Generation Intraventricular Impedance Imaging (Iii) System Evaluation Of Hardware Design”, Engineering in Medicine and Biology Society,. Proceedings of the 19th Annual International Conference of the IEEE, Oct. 30-Nov. 2, 1997 vol. 6,(1997). |
Jain et al., “Cell Therapy Attenuates Deleterious Ventricular Remodeling and Improves Cardiac Performance After Myocardial Infarction”, Circulation, 103:1920-1927, 2001. |
Jalife, “Rotors and Spiral Waves in Atrial Fibrillation”, Journal of Cardiovascular Electrophysiology, 14:776-780, 2003. |
Jané et al., “Alignment Methods for Averaging of High-Resolution Cardiac Signals: A Comparative Study of Performance”, IEEE Transactions on Biomedical Engineering, 38(6):571-579, 1991. |
Jia et al., “Electrophysiologic Endocardial Mapping From a Noncontact Nonexpandable Catheter: A Validation Study of a Geometry-Based Concept”, Journal of Cardiovascular Electrophysiology, 11:1238-1251, 2000. |
Kikuchi et al., “Targeted Modification of Atrial Electrophysiology by Homogeneous Transmural Atrial Gene Transfer”, Circulation, 111:264-270, 2005. |
Kistler et al., “Validation of Three-Dimensional Cardiac Image Integration: Use of Integrated CT Image into Electroanatomic Mapping System to Performa Catheter Ablation of Atrial Fibrillation”, Journal of Cardiovascular Electrophysiology, 17:341-348, 2006. |
Kun, Stevan et al., “Conductance Volumetric Model Of An Eccentrically Positioned Catheter Within A Three-Compartment Ellipsoidal Ventricle”, U, IEEE Transactions on, Jun. 1993, vol. 40, Issue: 6. |
Laciar et al., “Improved Alignment Method for Noisy High-Resolution ECG and Holter Records Using Multiscale Cross-Correlation”, IEEE Transactions on Biomedical Engineering, 50(3):344-353, 2003. |
Liu et al., “Endocardial Potential Mapping From a Noncontact Nonexpandable Catheter: A Feasibility Study”, Annals of Biomedical Engineering, 26:994-1009, 1998. |
Lorensen et al. “Marching Cubes: A High Resolution 3D Surface Construction Algorithm”. Computer Graphics 21(4):163-169, Jul. 1987. |
Mäkelä et al., “A Review of Cardiac Image Registration Methods”, IEEE Transactions on Medical Imaging, 21(9):1011-1021, 2002. |
Malladi, R. et al., “A Geometric Approach To Segmentation And Analysis Of 3D Medical Images”, Mathematical Methods in Biomedical Image Analysis, Proceedings of the Workshop on, Jun. 21-22, 1996, pp, 244-252, (1996). |
Mangan, Alan et al., “Partitioning 3D Surface Meshes Using Watershed Segmentation”, IEEE Transactions on Visualization and Computer Graphics, vol. 05, No. 4, pp. 308-321, (Oct.-Dec. 1999). |
Meininger et al., “Initial Experience with a Novel Focused Ultrasound Ablation System for Ring Ablation Outside the Pulmonary Vein”, Journal of Interventional Cardiac Electrophysiology, 8:141-148, 2003. |
Merrill, Daniel R. et al., “Electrical stimulation of excitable tissue: design of efficacious and safe protocols”, Journal of Neuroscience Methods, vol. 141, pp. 171-198 (2005). |
Miller, “Editor's Forum—Application of Registration for Ablation: A Marriage of Technologies”, Journal of Interventional Cardiac Electrophysiology, 11:87-89, 2004. |
Nademanee et al., “A New Approach for Catheter Ablation of Atrial Fibrillation: Mapping of the Electrophysiologic Substrate”, Journal of the American College of Cardiology, 43(11):2044-2053, 2004. |
Noseworthy et al., “The Impact of Respiration on Left Atrial and Pulmonary Venous Anatomy: Implications for Image-Guided Intervention”, Heart Rhythm, 2:1173-1178, 2005. |
Notification of Transmittal of the International Search Report and the Written Opinion of the International Searching Authority, or the Declaration, PCT/US2009/036099, Dated Apr. 28, 2009, 21 pages. |
Pappone et al., “Robotic Magnetic Navigation for Atrial Fibrillation Ablation”, Journal of the American College of Cardiology, 47(7):1390-1400, 2006. |
Paragios, “A Level Set Approach for Shape-Driven Segmentation and Tracking of the Left Ventricle”, IEEE Transactions on Medical Imaging, 22(6):773-776, 2003. |
Persson et al., “A Simple Mesh Generator in MATLAB”, SIAM Review, 46(2):329-345, 2004. |
Persson, “Mesh Generation for Implicit Geometries”, Massachusetts Institute of Technology—Thesis, Feb. 5, 2006. |
Pham, Dzung et al., “Current Methods In Medical Image Segmentation”, Annu. Rev. Biomed. Eng., 02: pp. 315-337, (2000). |
Rao et al., “Novel Noncontact Catheter System for Endocardial Electrical and Anatomical Imaging”, Annals of Biomedical Engineering, 32(4):573-584, 2004. |
Reddy et al., “Integration of Cardiac Magnetic Resonance Imaging With Three-Dimensional Electroanatomic Mapping to Guide Left Ventricular Catheter Manipulation—Feasibility in a Porcine Model of Healed Myocardial Infarction”, Journal of the American College of Cardiology, 44(11):2202-2213, 2004. |
Reddy et al., “Use of a Diode Laser Balloon Ablation Catheter To Generate Circumferential Pulmonary Venous Lesions in an Open-Thoracotomy Caprine Model”, PACE, 27:52-57, 2004. |
Sanders et al., “Spectral Analysis Identifies Sites of High-Frequency Activity Maintaining Atrial Fibrillation in Humans”, Circulation, 112:789-797, 2005. |
Sethian. “Level Set Methods and Fast Marching Methods: Evolving Interfaces in Computational Geometry, Fluid Mechanics, Computer Vision, and Materials Science”. Department of Mathematics—University of California, Berkeley. Cambridge University Press, 1999. |
Simon et al., “Electroanatomic Mapping of the Right Atrium With a Right Atrial Basket Catheter and Three-Dimensional Intracardiac Echocardiography”, PACE, 27:318-326, 2004. |
Smits et al., “Catheter-Based Intramyocarial Injection of Autologous Skeletal Myoblasts as a Primary Treatment of Ischemic Heart Failure”, Journal of the American College of Cardiology, 42(12):2063-2069, 2003. |
Solomon et al., “Real-Time Cardiac Catheter Navigation on Three-Dimensional CT Images”, Journal of Interventional Cardiac Electrophysiology, 8:27-36, 2003. |
Sra, Jasbir et al, “Registration Of 3D Computed Tomographic Images With Interventional Systems: Implications For Catheter Ablation Of Atrial Fibrillation”, J Intery Card Electrophysiol, 16: pp. 141-148, (2006). |
Sra et al., “Registration of Three-Dimensional Left Atrial Computed Tomographic Images With Projection Images Obtained Using Fluoroscopy”, Circulation, 112:3763-3768, 2005. |
Stevenson, “Radiofrequency Catheter Ablation of Ventricular Tachycardia After Myocardial Infarction”, Circulation, 98:308-314, 1998. |
Taccardi et al., “A New Intracavitary Probe for Detecting the Site of the Origin of Ectopic Ventricular Beats During One Cardiac Cycle”, Circulation, 75(1):272-281, 1987. |
Thal et al., “Novel Applications in Catheter Ablation”, Journal of Interventional Cardiac Electrophysiology, 13:17-21, 2005. |
Thiagalingam et al., “Noncontact Mapping of the Left Ventricle: Insights from Validation With Transmural Contact Mapping”, PACE, 27:570-578, 2004. |
Voth, “The Inverse Problem of Electrocardiography: Industrial Solutions and Simulations”, BEM and NFSI Conference Proceedings, Minneapolis, MN, May 12-15, 2005, pp. 191-194. |
Wittkampf et al., “LocaLisa: New Technique for Real-Time 3-Dimensional Localization of Regular Intracardiac Electrodes”, Circulation, 99:1312-1317, 1999. |
Yezzi, Anthony et al., “A Geometric Snake Model for Segmentation”, IEEE Transactions on Medical Imaging, vol. 16, No. 2, (Apr. 1997). |
International Search Report and the Written Opinion, PCT/US07/70854, Sep. 12, 2008, 15 pages. |
International Search Report and the Written Opinion, PCT/US08/52385, Aug. 8, 2008, 11 pages. |
Arthur, “Clinical Use of Intracardiac Impedance: Current Applications and Furture Perspectives”, PACE, vol. 24:500-506, Apr. 2001. |
Authorized officer Simin Baharlou, International Preliminary Report on Patentability in PCT/US2010/027568, mailed Nov. 3, 2011, 5 pages. |
Authorized officer Sung Chan Chung, International Search Report and Written Opinion in PCT/US2010/027568, mailed Nov. 4, 2010, 8 pages. |
Badics, “Real-Time Reconstruction of Endocardial Potential Maps in Non-Contact Cardiac Mapping”, International Journal for computation and Mathematics in Electrical Engineering (COMPEL), vol. 28, No. 4, 2009. |
Cheney et al., “Electrical Impedance Tomography”, SIAM Review 41:85-101, 1999. |
E. J. Haug, K. K. Choi, V. Komkov: Design Sensitivity Analysis of Structural Systems, Mathematics in Science and Engineering, vol. 177 (1986). |
Kuklik et al., “The reconstruction, from a set of points, and analysis of the interior surface of the heart chamber”, Physiol. Meas. 25:617-627, 2004. |
L. Piegl, W. Tiller: The NURBS Book, 2nd Edition, Springer (1997). |
Japanese Office Action, with English Translation, JP Application No. 2012-507231, Jun. 3, 2014, 9 pages. |
European Search Report issued in EP Application No. 10767478.0 mailed Jun. 10, 2015, 6 pages. |
Authorized officer Sung Chan Chung, Search Report and Written Opinion in PCT/US2010/027568, mailed Nov. 4, 2010, 8 pages. |
Number | Date | Country | |
---|---|---|---|
20100274150 A1 | Oct 2010 | US |