CLUSTERING OF ELECTROPHYSIOLOGICAL (EP) SIGNALS USING SIMILARITIES AMONG ARRHYTHMOGENIC ACTIVATIONS

Abstract

A method includes defining multiple different types of arrhythmias. Similarity measures are defined for the types of arrhythmias. A set of EP signals is received, the set acquired in a heart of a patient. Using the similarity measures, the set of EP signals is partitioned into at least two clusters, each cluster containing EP signals complying with a respective similarity measure for a respective type of arrhythmia. The clusters are indicated to a user.

Claims

1. A method for clustering of electrophysiological signals using similarities among arrhythmogenic activations, the method comprising: defining multiple different types of arrhythmias; defining similarity measures for the types of arrhythmias; receiving a set of EP signals acquired in a heart of a patient; using the similarity measures, partitioning the set of EP signals into at least two clusters, each cluster containing EP signals complying with a respective similarity measure for a respective type of arrhythmia; and indicating the clusters to a user.

2. The method according to claim 1, wherein indicating the clusters comprises presenting a representative EP signal of each cluster.

3. The method according to claim 1, wherein partitioning the set of EP signals comprises assigning the EP signals to the clusters based on a number, N, of activations in each EP signal over a given window of interest.

4. The method according to claim 1, wherein partitioning the set of EP signals comprises assigning the EP signals to the clusters based on a number of cross iso-peaks within an activation in each EP signal.

5. The method according to claim 1, wherein partitioning the set of EP signals comprises assigning the EP signals to the clusters based on a width of an activation in each EP signal.

6. The method according to claim 1, wherein partitioning the set of EP signals comprises assigning the EP signals to the clusters based on a peak-to-peak bi-polar voltage of an activation in each EP signal.

7. The method according to claim 1, wherein partitioning the set of EP signals comprises assigning the EP signals to the clusters based on a sharpest slope of an activation in each EP signal.

8. The method according to claim 1, wherein the similarity measure is one of a linear sum function and a quadratic sum function.

9. The method according to claim 1, wherein the EP signals are electrograms.

10. A system for clustering of electrophysiological signals using similarities among arrhythmogenic activations, the system comprising: a memory configured to store: definitions multiple different types of arrhythmias; and definitions of similarity measures for the types of arrhythmias; and a processor, which is configured to: receive a set of EP signals acquired in a heart of a patient; using the similarity measures, partition the set of EP signals into at least two clusters, each cluster containing EP signals complying with a respective similarity measure for a respective type of arrhythmia; and indicate the clusters to a user.

11. The system according to claim 10, wherein the processor is configured to indicate the clusters by presenting a representative EP signal of each cluster.

12. The system according to claim 10, wherein the processor is configured to partition the set of EP signals by assigning the EP signals to the clusters based on a number, N, of activations in each EP signal over a given window of interest.

13. The system according to claim 10, wherein the processor is configured to partition the set of EP signals by assigning the EP signals to the clusters based on a number of cross iso-peaks within an activation in each EP signal.

14. The system according to claim 10, wherein the processor is configured to partition the set of EP signals by assigning the EP signals to the clusters based on a width of an activation in each EP signal.

15. The system according to claim 10, wherein the processor is configured to partition the set of EP signals by assigning the EP signals to the clusters based on a peak-to-peak bi-polar voltage of an activation in each EP signal.

16. The system according to claim 10, wherein the processor is configured to partition the set of EP signals by assigning the EP signals to the clusters based on a sharpest slope of an activation in each EP signal.

17. The system according to claim 10, wherein the similarity measure is one of a linear sum function and a quadratic sum function.

18. The system according to claim 10, wherein the EP signals are electrograms.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0016] The present invention will be more fully understood from the following detailed description of the embodiments thereof, taken together with the drawings in which:

[0017] FIG. 1 is a schematic, pictorial illustration of a system for electrophysiological (EP) mapping and ablation, in accordance with an exemplary embodiment of the present invention;

[0018] FIG. 2 is a set of schematic graphs of EP signals, each graph representing a different type of cardiac EP signal, in accordance with an exemplary embodiment of the present invention; and

[0019] FIG. 3 is a flow chart that schematically illustrates a method of clustering EP signals using similarities among arrhythmogenic EP signals, in accordance with an exemplary embodiment of the present invention.

DETAILED DESCRIPTION OF EMBODIMENTS

Overview

[0020] The electrical activity of tissue of an organ of a patient, such as a cardiac chamber, can be mapped, e.g., electroanatomically and/or electrophysiologically, using a mapping catheter having suitable electrodes fitted at its distal end. An EP mapping system may acquire and analyze the signals output by the catheter.

[0021] In practice, during a cardiac arrhythmia analysis procedure, a very large number, for example, up to few tens of thousands or even more EP signals (e.g., bi-polar electrograms, may be acquired. Within these numerous EP signals, a physician may be interested in looking only at a small subset of the EP signals, for example electrograms showing some special fractionation, or electrograms characterized by a particular cycle length or local activation time (LAT). However, because of the large number of acquired EP signals, it is both difficult and time consuming to locate members of a subset among the entire set of acquired EP signals.

[0022] Embodiments of the present invention that are described hereinafter provide algorithms that automatically divide a set of EP signals into clusters (EP signals may be electrograms (EGM) or electrocardiograms (ECGs), for example). In an exemplary embodiment, a set of electrograms is acquired from a patient with arrhythmia. The electrograms are then analyzed to find parameters such as activation times, activation duration, number of peaks in an electrogram, slopes of the peaks, and the like. An algorithm uses these parameters to associate each electrogram with one of multiple types of predefined EP activations, such as normal conduction, atrial fibrillations, flutter, tachycardia, or other type of arrhythmia.

[0023] For each type of EP signal, a “similarity measure” is predefined, being the distance relating EP signals (e.g., electrograms) within the defined type of EP signal (e.g., type of arrhythmia). This mathematical similarity measure, for example one defined by a similarity metric as described below, delineates the bounds of the parameters defining “neighboring” signals. As an example, electrograms for a particular type of arrhythmia may be defined as those having two activations, each peak having a maximum voltage within some specified range, e.g., between 2 mV and 5 mV. An electrogram is deemed in the neighborhood of a given electrogram of this type if its two activations are within a given similarity difference (e.g., 10%) of those of the given electrogram.

[0024] In some embodiments, a processor runs a clustering algorithm that uses the predefined types and their respective distances to decide if a particular EP signal is clustered within one of the defined types. Typical clustering algorithms that may be used include the DBSCAN and the OPTICS algorithms, which were previously described in the literature. The two algorithms are described, for example, in a review paper entitled “Density-based clustering,” by Kriegel et al., WIREs Data Mining and Knowledge Discovery, Volume 1, Issue 3, May/June 2011, pages 231-240 of, which is incorporated herein by reference.

[0025] In some embodiments, the processor clusters the EP signals automatically without user intervention. Subdividing the EP signals into clustered groups (i.e., types) of signals, for example ten groups, reduces the amount of time a physician may need to examine the EP signals. Thus, the physician may only need to inspect ten cluster-type representative signals, rather than many or all of the acquired signals.

System Description

[0026] FIG. 1 is a schematic, pictorial illustration of a system 21 for electrophysiological (EP) mapping and ablation, in accordance with an embodiment of the present invention. FIG. 1 depicts a physician 27 using an EP mapping catheter 29 to perform an EP mapping of a heart 23 of a patient 25. Catheter 29 comprises, at its distal end, a multi-channel electrode array 50 comprising one or more arms 20, each of which is coupled to mapping-electrodes 22. During the mapping procedure, electrodes 22 acquire and/or inject signals from and/or to the tissue of heart 23. In particular, electrodes 22 acquire intracardiac EP signals, such as unipolar and/or bipolar electrograms.

[0027] Catheter 29 may be further used to perform an ablation, such as a radiofrequency (RF) or irreversible electroporation (IRE).

[0028] The respective locations of mapping-electrodes 22 inside heart 23 (i.e., where the electrograms are measured) are also tracked, such that a processor 28 may link each acquired electrogram with the location at which the signal was acquired. System 21 externally senses electrical position signals using a plurality of external electrodes 24 coupled to the body of patient 25; for example, three external electrodes 24 may be coupled to the patient's chest, and another three external electrodes may be coupled to the patient's back. For ease of illustration, only one external electrode is shown in FIG. 1.

[0029] An example of a system capable of using the sensed electrical position signals to track the locations of mapping-electrodes 22 inside heart 23 of the patient is the CARTO® 3 system (produced by Biosense Webster, Irvine, Calif.). The CARTO® 3 system uses a tracking method known as Advanced Current Location (ACL), which is described in detail in U.S. Pat. No. 8,456,182 whose disclosure is incorporated herein by reference.

[0030] A data acquisition module 33 receives the multiple electrograms conveyed to an electrical interface 35 over a wire link that runs in catheter 29. Processor 28 of system 21 receives these cardiac signals via the electrical interface 35, and uses the disclosed clustering algorithm based on the aforementioned similarity function to cluster these signals according to different predefined arrhythmia types, e.g., atrial fibrillations, tachycardias, flutter and more, which may amount, as in the example shown in FIG. 2, to about ten EP signal types (e.g., types of electrograms 40 shown on a display 26 of system 21). In this way physician 27 may need to inspect only a small number (e.g., ten) of the EP signals (e.g., signals 40), rather than many more of the signals acquired, in order to clinically assess the mapping results. Based on the information, physician 27 may, for example, decide to perform an ablation.

[0031] The definitions of the various types of arrhythmia, and the corresponding definitions of the similarity measures, are typically stored in a memory 31 for use by processor 28. Processor 28 also stores the clustered electrograms in memory 31 for further analysis, such as for constructing an EP map.

[0032] The exemplary embodiment shown in FIG. 1 is chosen purely for the sake of conceptual clarity. Other types of multi-electrode sensing geometries, such as of the Lasso® catheter (produced by Biosense Webster) or a basket catheter may also be employed. Additionally, contact sensors may be fitted at the distal end of electro-anatomical catheter 29 and transmit data indicative of the physical quality of electrode contact with tissue. In an exemplary embodiment, measurements of some electrodes 22 may be discarded because their physical contact quality is poor, and the measurements of other electrodes may be regarded as valid because their contact quality is high.

[0033] Processor 28 typically comprises a general-purpose computer with software programmed to carry out the functions described herein. The software may be downloaded to the computer in electronic form, over a network, for example, or it may, alternatively or additionally, be provided and/or stored on non-transitory tangible media, such as magnetic, optical, or electronic memory. In particular, processor 28 runs a dedicated algorithm that enables processor 28 to perform the steps described in FIG. 3.

Clustering of Ep Signals Using Similarities Among Arrhythmogenic Activations

[0034] As noted above, a clustering algorithm that is provided analyzes acquired EP signals to find parameters of the EP signals such as activation times, activation duration, number of peaks in an EP signal, and slopes of the peaks. The processor utilizes an algorithm that uses these parameters in a similarity measure to associate each EP signal, if found aberrative, with different predefined arrhythmia types.

[0035] In an exemplary embodiment, by-type clustering parameters for atrial flutter arrhythmia (for the similarity function) are:

[0036] Number of activations within the window of interest (WOI): N

[0037] Number of cross iso-peaks within an activation: P

[0038] Width of an activation: W

[0039] Peak to peak bipolar voltage of an activation: V

[0040] Sharpest slope of an activation: S

[0041] Between a pair of EP signals, a normalized value of the similarity function within the range of [0,1] is correspondingly calculated using the differences: dP, dW, dV, dS.

[0042] Assuming a stage clustering model, a first stage is characterized solely by a number N of activations within the WOI, and produce four clusters (flat: 0 activations; single potential, double potential, other: more than two activations).

[0043] The second stage of the model is calculated based on a similarity function F(EP_signal_1, EP_signal_2) on the rest of the parameters:

F(EP_sig_1,EP_sig_2)=0.5.Math.dP+0.25.Math.dS+0.125.Math.dV+0.125.Math.dW

[0044] The above linear similarity function is brought by way of example. Other functions, or norms, for example, quadratic ones in dP, dW, dV, and dS, may be used.

[0045] FIG. 2 is a set 41 of schematic graphs of EP signals, each graph representing a different cardiac activation type, in accordance with an exemplary embodiment of the present invention. FIG. 2 shows ten such schematic graphs representative of respective predefined clusters of cardiac EP signal types, such as types of arrhythmogenic activations. The EP signals shown are bipolar electrograms.

[0046] In set 41, different graphs may have a different number of activations 42, N, within a common WOI 202. For example, EP signals types 2 and 4 have N=2 activations 42 (e.g., peaks, complexes), whereas EP signals types 6 and 9 have N=3 peaks. In that regard, EP signal type 3 is unique by having N>4 peaks.

[0047] As another example, EP signals type 7 shows a large width of an activation, W, whereas EP signals types 4 and 10 show a narrower W.

[0048] A physician may need to inspect only set 41 of the ten-cluster representative EP signals types 1 to 10 in order to assess the cardiac activity and to determine if and which of the types has clinical significance. For example, most types may represent normal cardiac activities with only a small number of types representing suspicious electrical activity, abnormal substrate, possible ablation target, or other morphology of interest to a physician.

[0049] The graphs of FIG. 2 are schematic and brought purely by way of example to demonstrate the technique. In particular, set 41 of the shown graphs is not intended to represent a particular use case (e.g., to represent a diagnostic value case). Actual EP signals of actual arrhythmias may vary from the ones shown here.

[0050] FIG. 3 is a flow chart that schematically illustrates a method of clustering EP signals using similarities among arrhythmogenic EP signals, in accordance with an exemplary embodiment of the present invention. The algorithm, according to the presented exemplary embodiment, carries out a process that begins when a user, such as physician 27, defines, (e.g., selects using processor 28 with installed algorithm) a group of different arrhythmia types, at an arrhythmia type defining step 302.

[0051] Next, the algorithm selects a similarity measure defined for each arrhythmia type, with physician 27 optionally adjusting parameters of the selected one or more similarity measures, at similarity measures defining step 304.

[0052] At an EP data receiving step 306, processor 28 receives a set of thousands of EP signals acquired by a multi-electrode catheter (e.g., catheter 29) that mapped (or is mapping in real-time) a cardiac chamber.

[0053] Using the one or more similarity measures, at a clustering step 308, processor 28 partitions the set of EP signals into at least two clusters, each cluster characterized by its EP signals complying with the respective similarity measure for the type of arrhythmia.

[0054] Finally, processor 28 presents on display 26 representative graphs, such as set 41, of the different clusters, at EP signals presentation step 310. Typically, the processor selects a median signal from each cluster of signals as a representative signal for the cluster. Alternatively, however, any other suitable selection can be used.

[0055] Although the embodiments described herein mainly address cardiac applications, the methods and systems described herein can also be used in other applications, such as in mapping of electrical activity in the brain.

[0056] It will thus be appreciated that the embodiments described above are cited by way of example, and that the present invention is not limited to what has been particularly shown and described hereinabove. Rather, the scope of the present invention includes both combinations and sub-combinations of the various features described hereinabove, as well as variations and modifications thereof which would occur to persons skilled in the art upon reading the foregoing description and which are not disclosed in the prior art. Documents incorporated by reference in the present patent application are to be considered an integral part of the application except that to the extent any terms are defined in these incorporated documents in a manner that conflicts with the definitions made explicitly or implicitly in the present specification, only the definitions in the present specification should be considered.