Patent Application
20160175023
1. Field of the Invention
This invention relates to cardiac physiology. More particularly, this invention relates to the evaluation of electrical propagation in the heart.
2. Description of the Related Art
The meanings of certain acronyms and abbreviations used herein are given in Table 1.
Cardiac arrhythmias such as atrial fibrillation are an important cause of morbidity and death. Commonly assigned U.S. Pat. No. 5,546,951, and U.S. Pat. No. 6,690,963, both issued to Ben Haim and PCT application WO 96/05768, all of which are incorporated herein by reference, disclose methods for sensing an electrical property of heart tissue, for example, local activation time, as a function of the precise location within the heart. Data are acquired with one or more catheters having electrical and location sensors in their distal tips, which are advanced into the heart. Methods of creating a map of the electrical activity of the heart based on these data are disclosed in commonly assigned U.S. Pat. No. 6,226,542, and U.S. Pat. No. 6,301,496, both issued to Reisfeld, which are incorporated herein by reference. As indicated in these patents, location and electrical activity is typically initially measured on about 10 to about 20 points on the interior surface of the heart. These data points are then generally sufficient to generate a preliminary reconstruction or map of the cardiac surface. The preliminary map is often combined with data taken at additional points in order to generate a more comprehensive map of the heart's electrical activity. Indeed, in clinical settings, it is not uncommon to accumulate data at 100 or more sites to generate a detailed, comprehensive map of heart chamber electrical activity. The generated detailed 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.
Catheters containing position sensors may be used to determine the trajectory of points on the cardiac surface. These trajectories may be used to infer motion characteristics such as the contractility of the tissue. As disclosed in U.S. Pat. No. 5,738,096, issued to Ben Haim, which is incorporated herein in its entirety by reference, maps depicting such motion characteristics may be constructed when the trajectory information is sampled at a sufficient number of points in the heart.
Electrical activity at a point in the heart is typically measured by advancing a multiple-electrode catheter to measure electrical activity at multiple points in the heart chamber simultaneously. A record derived from time varying electrical potentials as measured by one or more electrodes is known as an electrogram. Electrograms may be measured by unipolar or bipolar leads, and are used, e.g., to determine onset of electrical propagation at a point, known as local activation time.
Sensors in a cardiac chamber may detect far-field electrical activity, i.e., the ambient electrical activity originating away from the sensors, which can distort or obscure local electrical activity, i.e., signals originating at or near the sensor location. Commonly assigned U.S. Patent Application Publication No. 2014/0005664 of Govari et al., which is herein incorporated by reference, discloses distinguishing a local component in an intracardiac electrode signal, due to the tissue with which the electrode is in contact from a remote-field contribution to the signal, and explains that a therapeutic procedure applied to the tissue can be controlled responsively to the distinguished local component.
U.S. Patent Application Publication No. 2014/0187991 of Thakur et al. proposes a method for mapping a cardiac chamber by sensing activation signals of intrinsic physiological activity with a plurality of electrodes disposed in or near the cardiac chamber. The method includes isolating R-wave events in the activation signals, generating a far-field activation template representative of a far-field activation signal component based on the R-wave events, and filtering the far-field activation template from the activation signals to identify near-field activation signal components in the activation signals.
There is provided according to embodiments of the invention a method of far field reduction, which is carried out by detecting peaks of a QRS complex of a cardiac electrogram, extracting unipolar beats of an intracardiac electrogram that occur within a predetermined time interval that includes respective detected peaks, constructing a first mean unipolar beat by averaging the extracted unipolar beats, cross correlating the extracted unipolar beats with the first mean unipolar beat by determining respective time-lags therebetween, and accepting ones of the extracted unipolar beats wherein a cross correlation thereof exceeds a predetermined correlation threshold. The method is further carried out by constructing a second mean unipolar beat from the accepted unipolar beats, determining a ventricular far field component by subtracting the second mean unipolar beat from the extracted unipolar beats, subtracting the ventricular far field component from the intracardiac electrogram to distinguish a local component of the intracardiac electrogram, and controlling a medical procedure responsively to the local component.
According to an aspect of the method, the predetermined correlation threshold includes a time lag that does not exceed 10 ms.
According to another aspect of the method, the predetermined correlation threshold includes a difference between peak amplitudes of the extracted unipolar beats and the first mean unipolar beat that does not exceed 10 per cent.
According to a further aspect of the method, constructing a first mean unipolar beat includes computing a median derivative beat by differentiating the extracted unipolar beats, and cross correlating is performed between the extracted unipolar beats and the median derivative beat.
According to one aspect of the method, constructing a second mean unipolar beat includes aligning the extracted unipolar beats.
According to yet another aspect of the method, constructing a second mean unipolar beat includes dividing the accepted unipolar beats into a plurality of clusters, and constructing respective cluster-specific mean unipolar beats for each of the clusters.
According to an additional aspect of the method, dividing the accepted unipolar beats is performed by k-means clustering according to Euclidean distances between the accepted unipolar beats.
According to still another aspect of the method, determining a ventricular far field component includes selecting an optimum one of the cluster-specific mean unipolar beats for respective accepted unipolar beats, and subtracting the selected cluster-specific mean unipolar beats from the extracted unipolar beats.
There is further provided according to embodiments of the invention an apparatus, including an intra-body probe, the probe being adapted to contact tissue in a heart, a display, and a processor. The processor is configured to receive an electrical signal from an electrode of the probe, to distinguish a local component of the signal that is due to the tissue with which the electrode is in contact from a remote-field contribution to the signal by performing the steps of: detecting peaks of a QRS complex of a cardiac electrogram, extracting unipolar beats of an intracardiac electrogram received via the electrode, wherein the extracted unipolar beats occur within a predetermined time interval that includes respective detected peaks, constructing a first mean unipolar beat by averaging the extracted unipolar beats, cross correlating the extracted unipolar beats with the first mean unipolar beat by determining respective time-lags therebetween, accepting ones of the extracted unipolar beats wherein a cross correlation thereof exceeds a predetermined correlation threshold, constructing a second mean unipolar beat from the accepted unipolar beats, determining a ventricular far field component by subtracting the second mean unipolar beat from the extracted unipolar beats, subtracting the ventricular far field component from the intracardiac electrogram to distinguish the local component of the intracardiac electrogram, and outputting the local component to the display.
There is further provided according to embodiments of the invention a method which is carried out by inserting a probe into a heart of a subject. A unipolar intracardiac electrode is disposed on the distal segment of the probe. The method is further carried out by concurrently receiving an intracardiac electrogram from a location in the heart with the intracardiac electrode and receiving another electrogram with a second electrode, forming a first template by averaging beats of the intracardiac electrogram in alignment with R-wave complexes of the other electrogram, computing respective correlations between the beats of the intracardiac electrogram and the first template, selecting beats of the intracardiac electrogram, wherein the respective correlations thereof exceed a predetermined value, forming a second template from the selected beats of the intracardiac electrogram, generating, based on the second template, a signal representing a far field component of the intracardiac electrogram, removing the far field component from the intracardiac electrogram to yield a filtered intracardiac electrogram, and reporting the filtered intracardiac electrogram.
According to yet another aspect of the method, forming a second template includes dividing the selected beats of the intracardiac electrogram into a plurality of clusters, and constructing a plurality of cluster-specific second templates for each of the clusters.
According to still another aspect of the method, dividing the selected beats of the intracardiac electrogram is performed by k-means clustering according to Euclidean distances between the selected beats.
According to a further aspect of the method, generating includes selecting an optimum one of the cluster-specific second templates for respective selected beats, and subtracting the optimum template from the selected beats, respectively.
For a better understanding of the present invention, reference is made to the detailed description of the invention, by way of example, which is to be read in conjunction with the following drawings, wherein like elements are given like reference numerals, and wherein:
In the following description, numerous specific details are set forth in order to provide a thorough understanding of the various principles of the present invention. It will be apparent to one skilled in the art, however, that not all these details are necessarily needed for practicing the present invention. In this instance, well-known circuits, control logic, and the details of computer program instructions for conventional algorithms and processes have not been shown in detail in order not to obscure the general concepts unnecessarily.
Turning now to the drawings, reference is initially made to
The system 10 may comprise a general purpose or embedded computer processor, which is programmed with suitable software for carrying out the functions described hereinbelow. Thus, although portions of the system 10 shown in other drawing figures herein are shown as comprising a number of separate functional blocks, these blocks are not necessarily separate physical entities, but rather may represent, for example, different computing tasks or data objects stored in a memory that is accessible to the processor. These tasks may be carried out in software running on a single processor, or on multiple processors. The software may be provided to the processor or processors on tangible non-transitory media, such as CD-ROM or non-volatile memory. Alternatively or additionally, the system 10 may comprise a digital signal processor or hard-wired logic.
The catheter 14 typically comprises a handle 20, having suitable controls on the handle to enable the operator 16 to steer, position and orient the distal end of the catheter as desired to the ablation. To aid the operator 16, the distal portion of the catheter 14 contains position sensors (not shown) that provide signals to a positioning processor 22, located in a console 24. The catheter 14 may be adapted, mutatis mutandis, from the ablation catheter described in commonly assigned U.S. Pat. No. 6,669,692, whose disclosure is herein incorporated by reference. The console 24 typically contains an ECG processor 26 and a display 30.
The positioning processor 22 measures location and orientation coordinates of the catheter 14. In one embodiment, the system 10 comprises a magnetic position tracking system that determines the position and orientation of the catheter 14. The system 10 typically comprises a set of external radiators, such as field generating coils 28, which are located in fixed, known positions external to the patient. The coils 28 generate electromagnetic fields in the vicinity of the heart 12. These fields are sensed by magnetic field sensors located in the catheter 14.
Typically, the system 10 includes other elements, which are not shown in the figures for the sake of simplicity. For example, the system 10 may include an electrocardiogram (ECG) monitor, coupled to receive signals from one or more body surface electrodes, to provide an ECG synchronization signal to the console 24. The system 10 typically also includes a reference position sensor, either on an externally-applied reference patch attached to the exterior of the subject's body, or on an internally-placed catheter, which is inserted into the heart 12 maintained in a fixed position relative to the heart 12. Conventional pumps and lines for circulating liquids through the catheter 14 for cooling an ablation site may be provided.
One system that embodies the above-described features of the system 10 is the CARTO® 3 System, available from Biosense Webster, Inc., 3333 Diamond Canyon Road, Diamond Bar, Calif. 91765. This system may be modified by those skilled in the art to embody the principles of the invention described herein. Multi-electrode basket and spline catheters are known that are suitable for obtaining unipolar and bipolar electrograms. An example of such a spline catheter is the Pentaray® NAV catheter, available from Biosense Webster.
One embodiment of a time ensemble averaging (TEA) algorithm for removing ventricular far field effects may be summarized as follows:
9. Reduce estimated VFF from UEGMn channel by subtracting VFF′est from unipolar channel.
Reference is now made to
The sequence shown in
Reference is now made to
Block 48 accepts the time-stamp signal 50 and the UEGMn input 36. Beats of the UEGMn channel are extracted, using the time-stamp signal 50 as a reference. Block 48 outputs an intracardiac electrode signal 52, which forms an input to blocks 54, 56.
Several intracardiac beats are averaged in block 54, which produces a template signal 58 that consists of a moving average of the last m intracardiac beats that represents the moving average. The value of m may be configured by the operator. Averaging the beats over a 30 sec time interval is suitable.
The beats in the intracardiac electrode signal 52 are correlated with the template signal 58 in block 56, which outputs an intracardiac electrode signal 60 that includes correlation information. Signals having a high correlation with the template signal 58 are retained and aligned with the template to obtain the best correlation in a range of −3 to +3 ms. Signals failing to show a high correlation are discarded.
The intracardiac electrode signal 60 is input to block 62, which aligns and selects intracardiac beats for further processing. In block 62 beats failing to show ventricular activity, i.e., ventricular far field activity are rejected as they correlate poorly with the template signal 58. The rejected beats typically only show atrial activity. Block 62 outputs an intracardiac beat signal 64 that is highly correlated with the template signal 58.
The intracardiac beat signal 64 is an input to block 66, which produces a moving average of the last n intracardiac beats and outputs a template signal 68 that represents a moving average of the beats that were selected in block 62. The value of n is configurable. Typically the value n represents beats selected from 30 sec of cardiac activity.
A time-varying signal that estimates the ventricular far field component of the UEGMn input 36 is generated in block 70 and is output as the estimated ventricular far field output 38.
Another embodiment of a time ensemble averaging (TEA) algorithm for removing ventricular far field effects may be summarized as follows:
1. Extract unipolar beats recorded on unipolar IC channel (UEGMN) and its neighboring electrodes around R peak detection timestamps.
2. Mean average all extracted beats from electrode N to a single beat to generate base template.
3. Calculate median filter on all extracted beats derivatives from Electrode N and its neighboring electrodes.
4. Calculate the Cross correlation between all extracted beats derivatives and the median derivative template.
5. Reject all beats that are blow correlation threshold.
6. Operate k-Means algorithm to divide the remaining beats into k clusters based on the Euclidean distance between them.
7. Mean average each cluster beats to k templates.
8. Select best template for VFF reduction per each beat in Electrode N based on cross correlation with template & minimum RM after reduction with template.
9. Use the selected template to construct the estimated (VFF′est) signal for reduction.
10. Reduce estimated VFF from UEGMN channel by subtracting VFF′est from unipolar channel.
Reference is now made to
A version of the time ensemble averaging is performed in block 72, which accepts the input 36 and the time-stamp signal 50 (
Reference is now made to
The signal 82 is subjected to mean averaging in block 86 to form a template signal 88. The template signal 88 is processed in block 110, which is described below. Typically the template signal 88 represents an average of all beats of one electrode during a recording interval, and typically comprises about 500 beats.
The signal 84 is differentiated with respect to time and the derivatives of each of the nine electrodes subjected to median averaging in block 90. A median template signal 92 generated from the median averaging process is output from the block 90.
The derivatives that are output as the median template signal 92 are correlated with respective intracardiac beats of the electrode signal 84 in block 94. The correlation involves squaring the signals from each of the 9 electrodes. Correlation results are output as signal 96 to block 98, where intracardiac beats having a correlation coefficient that is below a predefined threshold are dropped. Such poorly correlated beats typically have a high component of atrial activity and are undesirable for obtaining an isolated ventricular far field component. The remaining intracardiac beats of the signal 84 are output as signal 100.
The signal 100 is an input to block 102 where the filtered intracardiac beats are classified into k groups using k-means clustering, each group being 50 ms in duration. The value of k is 3 in a current embodiment. K-means clustering is a partitioning method that partitions observed data into k mutually exclusive clusters using the Euclidean distance as a metric, and returns the index of the cluster to which it has assigned each observation. The k groups are output as a k-cluster signal 104. A mean is calculated for each of the k groups in block 106. K templates and the individual filtered beats are output as a signal 108.
The k-template signal 108 and the template signal 88 along with the individual beats and their cluster assignments are input to a selection block 110 where each filtered beat is associated with the most suitable one of the k templates of the k-template signal 108 or the template signal 88 and output as signal 112. A ventricular far field estimate is developed on a per beat basis in block 114 using the signal 112 and output as the estimated ventricular far field signal 74, where it is processed in block 76 (
Reference is now made to
It will be appreciated by persons skilled in the art 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 that are not in the prior art, which would occur to persons skilled in the art upon reading the foregoing description.
