System and method for noninvasive electrocardiographic imaging (ECGI)

• US 7,983,743 B2
• Filed: 07/21/2006
• Issued: 07/19/2011
• Est. Priority Date: 07/22/2005
• Status: Active Grant

• Current Assignee: Case Western Reserve University, POD Suppliers, LLC
• Original Assignee: Case Western Reserve University, Washington University In St Louis
• Inventors: Wang, Yong, Rudy, Yoram, Jia, Ping
• Primary Examiner(s): Layno; Carl H
• Assistant Examiner(s): Stice; Paula J
• Application Number:
  US11/996,441
• Publication Number:
  US 20080213128A1
• Time in Patent Office: 1,824 Days
• Field of Search: 600/509, 600/513, 600/523
• US Class Current: 600/509
• CPC Class Codes: A61B 5/0044 for the heart
• CPC Class Codes: A61B 5/055 involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
• CPC Class Codes: A61B 5/282 Holders for multiple electrodes for introduction into the body A61B5/287
• CPC Class Codes: A61B 5/318 Heart-related electrical modalities, e.g. electrocardiography [ECG]
• CPC Class Codes: A61B 6/03 Computed tomography [CT] echo-tomography A61B8/14

Abstract

Noninvasive systems and methods are provided for determining electrical activity for a heart of a living being. A processor is configured to meshlessly compute data that represents heart electrical activity from a set of noninvasively measured body surface electrical potentials. This is accomplished using data that describes a geometric relationship between a plurality of locations corresponding to where the body surface electrical potentials were measured and the heart.

Claims

1. A noninvasive system for determining electrical activity for a heart of a living being, the system comprising:

memory to store computer-readable instructions; and

a processor configured to access the memory and execute the computer readable instructions to compute heart electrical activity data that represents heart electrical activity of the living being from a set of noninvasively measured body surface electrical potentials using data that describes a geometric relationship between a plurality of locations corresponding to where the body surface electrical potentials were measured and a surface representation of the heart;

wherein the processor executes the computer readable instructions to compute the heart electrical activity data by translating the measured body surface electrical potentials to surface potentials on the surface representation of the heart without using a mesh of any heart surface.

2. The system of claim 1 wherein the heart electrical activity data comprises a set of epicardial envelope electrical potentials.

3. The system of claim 2 wherein the processor is further configured to meshlessly compute the set of epicardial envelope electrical potentials via a method of fundamental solution (MFS).

4. The system of claim 2 wherein the living being comprises a torso having a body surface, the system further comprising:

an electrode array system in communication with the processor for noninvasively measuring electrical potentials at a plurality of locations on the torso body surface via a plurality of electrodes applied to the torso to provide a set of noninvasively measured body surface potentials;

a geometry determining device in communication with the processor, the geometry determining device being configured to (1) determine a geometry of the locations on the body surface where the electrical potentials were noninvasively measured, (2) determine a geometry of a heart surface of the living being, and (3) communicate the determined geometries of the locations on the body surface and the heart surface to the processor; and

wherein the processor is further configured to compute the set of epicardial envelope electrical potentials from the set of noninvasively measured body surface potentials and the determined geometries.

5. The system of claim 4 wherein the processor is further configured to compute the set of the epicardial cardiac surface electrical potentials in near-real time.

6. The system of claim 4 wherein the processor is further configured to (1) determine a plurality of epicardial nodes that define locations on an epicardial cardiac surface for which the computed epicardial cardiac surface electrical potentials apply, (2) determine a plurality of source nodes, wherein a plurality of the source nodes define a plurality of locations along a surface outside the torso and wherein another plurality of the source nodes define a plurality of locations along a surface inside the epicardial cardiac surface, (3) determine a matrix of coefficients A that relates each electrode location to each source node location, (4) perform an inverse computation on the matrix of coefficients A and the noninvasively measured body surface potentials to compute a plurality of source node coefficients, (5) determine a matrix of coefficients B that relates each epicardial node location to each source node location, and (6) perform a forward computation using the matrix of coefficients B and the source node coefficients to compute the set of epicardial cardiac surface electrical potentials.

7. The system of claim 6 wherein the processor is further configured to dynamically determine the source nodes.

8. The system of claim 6 wherein the processor is further configured to statically determine the source nodes by (1) defining each source node that is outside the body surface such that it is located a predetermined distance outward from a corresponding electrode location on a ray extending from a calculated center of the epicardial cardiac surface through a corresponding electrode location, and (2) defining each source node that is inside the epicardial cardiac surface such that it is located a predetermined distance inward from a corresponding epicardial node location on a ray extending from the calculated center of the epicardial cardiac surface through the corresponding epicardial node location.

9. The system of claim 6 wherein the processor is further configured to generate an epicardial cardiac surface potential map from the set of computed epicardial cardiac surface electrical potentials.

10. The system of claim 6 wherein the processor is further configured to compute a plurality of sets of the epicardial cardiac surface electrical potentials over a time duration from a plurality of successively noninvasively measured body surface electrical potentials.

11. The system of claim 10 wherein the processor is further configured to generate at least one selected from the group consisting of a plurality of electrograms and an isochrone from the sets of computed epicardial cardiac surface electrical potentials.

12. The system of claim 10 wherein the processor is further configured to generate, from the sets of computed epicardial cardiac surface electrical potentials, at least one selected from the group consisting of a recovery map, an integral map, and an activation-recovery interval map.

13. The system of claim 1 wherein the processor is integrated into a medical imaging platform.

14. A method for noninvasively reconstructing electrical activity of a heart of a living being, the living being having a body surface, the method comprising:

computing heart electrical activity data corresponding to heart electrical activity from a set of noninvasively measured body surface electrical potentials via an algorithm that translates the noninvasively measured body surface electrical potentials to the heart electrical activity data, the algorithm not using a mesh of any heart surface and not using a mesh of the body surface; and

generating a display based on the heart electrical activity data.

15. The method of claim 14 wherein the algorithm comprises a method of fundamental solution (MFS).

16. The method of claim 14 further comprising:

noninvasively measuring electrical potentials at a plurality of locations on the body surface;

determining a geometry of the plurality of locations on the body surface where the electrical potentials were noninvasively measured; and

determining a geometry of an epicardial envelope of the living being; and

wherein the computing comprises computing epicardial envelope electrical potential estimates from the set of noninvasively measured body surface potentials and the determined geometries of the plurality of locations on the body surface and the epicardial envelope.

17. The method of claim 16 wherein the computing further comprises:

configuring a plurality of epicardial nodes that define the locations on the epicardial envelope for which the computed epicardial envelope electrical potential estimates apply;

configuring a plurality of source nodes, wherein a plurality of the source nodes define a plurality of locations along a surface outside a torso of the living being and wherein another plurality of the source nodes define a plurality of locations along a surface inside the epicardial envelope;

determining a matrix of coefficients A that relates each electrode location to each source node location;

performing an inverse computation on the matrix of coefficients A and the noninvasively measured body surface electrical potentials to compute a plurality of source node coefficients;

determining a matrix of coefficients B that relates each epicardial node location to each source node location; and

performing a forward computation using the matrix of coefficients B and the source node coefficients to compute the epicardial envelope electrical potential estimates.

18. The method of claim 17 wherein the source node configuring comprises (1) defining each source node that is outside the torso such that it is located a predetermined distance outward from a corresponding electrode location on a ray extending from a calculated center of the epicardial envelope through the corresponding electrode location, and (2) defining each source node that is inside the heart surface such that it is located a predetermined distance inward from a corresponding epicardial node location on a ray extending from the calculated center of the epicardial envelope through the corresponding epicardial node location.

19. The method of claim 17 further comprising:

generating an epicardial cardiac surface potential map from the computed epicardial cardiac surface electrical potential estimates.

20. A non-transitory computer readable medium for use in connection with noninvasively computing electrical activity of a heart of a living being, the living being having a body surface, the computer readable medium comprising:

a code segment executable by a processor for computing data corresponding to heart electrical activity from a set of noninvasively measured body surface electrical potentials via a meshless algorithm that translates the noninvasively measured body surface electrical potentials to the heart electrical activity data;

a code segment executable by a processor for determining a geometry of locations on the body surface where the body surface electrical potentials were noninvasively measured; and

a code segment executable by a processor for determining a geometry of an epicardial envelope of the living being; and

wherein the computing code segment is configured to compute the set of epicardial envelope electrical potentials from the set of noninvasively measured body surface potentials and the determined geometries.

21. The computer readable medium of claim 20 wherein the meshless algorithm comprises a method of fundamental solution (MFS).

22. The computer readable medium of claim 20 wherein the body surface comprises the torso of the living being, and wherein the locations on the body surface where the set of noninvasively measured body surface electrical potentials were noninvasively measured correspond to locations of a plurality of electrodes that are applied to the living being'"'"'s torso.

23. The computer readable medium of claim 20 wherein the computing code segment further comprises:

a code segment executable by a processor for defining a plurality of epicardial nodes that correspond to the locations on the epicardial envelope for which the computed epicardial envelope electrical potentials apply;

a code segment executable by a processor for defining a plurality of source nodes, wherein a plurality of the source nodes correspond to a plurality of locations along a surface outside the torso and wherein another plurality of the source nodes correspond to a plurality of locations along a surface inside the epicardial envelope;

a code segment executable by a processor for determining a matrix of coefficients A that relates each electrode location to each source node location;

a code segment executable by a processor for performing an inverse computation on the matrix of coefficients A and the noninvasively measured body surface electrical potentials to compute a plurality of source node coefficients;

a code segment executable by a processor for determining a matrix of coefficients B that relates each epicardial node location to each source node location; and

a code segment executable by a processor for performing a forward computation using the matrix of coefficients B and the source node coefficients to compute the set of epicardial envelope electrical potentials.

24. The computer readable medium of claim 23 wherein the heart has an epicardial cardiac surface, wherein the epicardial envelope electrical potentials comprise a set of epicardial cardiac surface electrical potentials, and wherein the code segment for defining epicardial nodes is further configured to define the epicardial nodes such that they correspond to the locations on the epicardial cardiac surface for which the computed epicardial cardiac surface electrical potentials apply.

25. The computer readable medium of claim 24 wherein the source node configuring code segment is configured to statically define the source nodes.

26. The computer readable medium of claim 25 wherein the statically defining code segment is configured to (1) define each source node that is outside the torso surface such that it is located a predetermined distance outward from a corresponding electrode location on a ray extending from a calculated center of the epicardial cardiac surface through the corresponding electrode location, and (2) define each source node that is inside the epicardial cardiac surface such that it is located a predetermined distance inward from a corresponding epicardial node location on a ray extending from the calculated center of the epicardial cardiac surface through the corresponding epicardial node location.

27. A method for reconstructing electrical potentials on a surface inside a volume, wherein the volume comprises a field that satisfies Laplace equation, the volume having an outer surface, the method comprising:

computing data corresponding to a plurality of electrical potentials on a surface inside the volume from a set electrical potentials measured at the volume'"'"'s outer surface via a meshless algorithm that translates the measured outer surface electrical potentials to the inside surface electrical potential data without using a mesh of the surface inside the volume and without using a mesh of the outer surface of the volume; and

generating a display based on the heart electrical activity data.

28. The method of claim 27 wherein the inside surface electrical potential data comprises a plurality of electrical potentials at different locations along the inside surface.

29. The method of claim 28 further comprising:

applying a plurality of electrodes to the outer surface to measure the outer surface electrical potentials;

defining a plurality of inside surface nodes that correspond to the locations on the inside surface for which the computed inside surface electrical potentials apply;

defining a plurality of source nodes, wherein a plurality of the source nodes correspond to a plurality of locations along a surface outside the outer surface and wherein another plurality of the source nodes correspond to a plurality of locations along a surface inside the inside surface;

determining a matrix of coefficients A that relates each outer surface location where an outer surface electrical potential was measured to each source node location;

performing an inverse computation on the matrix of coefficients A and the measured outer surface electrical potentials to compute a plurality of source node coefficients;

determining a matrix of coefficients B that relates each inside surface node location to each source node location; and

performing a forward computation using the matrix of coefficients B and the source node coefficients to compute the inside surface electrical potentials.

30. The method of claim 27 wherein the volume comprises a human torso.

31. The method of claim 27 wherein the inside surface comprises a human epicardial cardiac surface, and wherein the meshless algorithm comprises a method of fundamental solution (MFS).

32. A system for reconstructing heart electrical activity for a heart of a living being, the living being having a torso, the torso having an outer surface, the system comprising:

memory to store computer-readable instructions; and

a processor configured to access the memory and execute the computer readable instructions to (1) receive data representing a plurality of electrical potentials measured on the torso outer surface at a plurality locations along the torso outer surface, (2) receive spatial relationship data that describes a spatial relationship between the torso outer surface and an epicardial envelope, (3) determine the torso outer surface locations where the electrical potentials were measured, (4) define a plurality of locations along the epicardial envelope, (5) define a plurality of source node locations outside the torso outer surface, (6) define a plurality of source node locations inside the epicardial envelope, (7) based at least in part upon the spatial relationship data, the determined torso outer surface locations, and the defined source node locations, compute a matrix of coefficients A that spatially relates each determined torso outer surface location to each defined source node location, (8) perform an inverse computation on A and the received electrical potential data to compute a plurality of source node location coefficients, (9) based at least in part upon the spatial relationship data, the defined epicardial envelope locations, and the defined source node locations, compute a matrix of coefficients B that relates each defined epicardial envelope location to each defined source node location, (10) perform a forward computation using B and the source node coefficients to compute electrical activity data that represents electrical activity on the epicardial envelope, and (11) store the electrical activity data in the memory.

33. The system of claim 32 further comprising:

a plurality of electrodes for application to the torso outer surface to measure electrical potentials at the plurality of torso outer surface locations; and

a signal acquisition and processing device, wherein the electrodes are in communication with the processor via the signal acquisition and processing device, the signal acquisition and processing device being configured to acquire and process the measured electrical potentials to a data format suitable for processing by the processor.

34. The system of claim 33 further comprising a geometry determining device in communication with the processor, the geometry determining device being configured to acquire and provide the spatial relationship data to the processor.

35. The system of claim 33 wherein the processor is further configured to automatically define the plurality of epicardial envelope locations.

36. The system of claim 35 wherein the processor is further configured to automatically define the plurality of epicardial envelope locations by evenly distributing the defined epicardial envelope locations along the epicardial envelope.

37. A non-transitory computer readable medium having computer executable instructions for performing a method that comprises:

computing heart electrical activity data that represents heart electrical activity from a set of noninvasively measured body surface electrical potentials using data that describes a geometric relationship between a plurality of locations corresponding to where the body surface electrical potentials were measured and a representation of a surface of a heart of a living being, wherein the computing comprising;

performing a first computation on a first transfer matrix and the noninvasively measured body surface electrical potentials to compute a first plurality of coefficients, the first computation being performed without using a mesh of any heart surface; and

performing a second computation using a second transfer matrix and the first plurality of coefficients to compute a set of cardiac surface electrical potentials for the representation of the surface of the heart, the second computation being performed without using a mesh of any heart surface,

the set of cardiac surface electrical potentials for the representation of the surface of the heart being stored in memory.