Retrospective optimization of ECG signal mapping by detecting discrepancies

Abstract

To provide intra-body cardiac electrocardiogram (ECG) sensing.SOLUTION: A system includes signal acquisition circuitry and a processor. The signal acquisition circuitry is configured to receive multiple intra-cardiac signals acquired by multiple electrodes of an intra-cardiac probe in a heart of a patient. The processor is configured to perform a sequence of annotation-visualization operations at subsequent times, by performing the following steps in each operation: extracting multiple annotation values from the intra-cardiac signals; selecting a group of the intra-cardiac signals; identifying in the group one or more annotation values that are statistically deviant by more than a predefined measure of deviation; and visualizing the annotation values to a user, while omitting and refraining from visualizing the statistically deviant annotation values. The processor is further configured to assess, over one or more of the annotation-visualization operations, a rate of omissions of annotation values, and to take a corrective action in response to detecting that the rate of omissions exceeds a predefined threshold.SELECTED DRAWING: Figure 1

Description

[Technical field]

[0001] (CROSS REFERENCE TO RELATED APPLICATIONS) This application is related to U.S. patent application entitled "Using Statistical Characteristics of Multiple Grouped ECG Signals to Detect Inconsistent Signals" filed on even date (Attorney Docket No. BIO6196USNP1), the disclosure of which is incorporated herein by reference.

[0002] FIELD OF THEINVENTION The present invention relates generally to internal medical procedures and devices, and more particularly to internal sensing of the electrocardiogram (ECG) of the heart. [Background technology]

[0003] When measuring and annotating intracardiac electrocardiogram (iECG) signals generated by multiple electrodes, it may be desirable to process (e.g., computer process) the signals to reduce noise embedded in the signals.

[0004] Various methods exist for such iECG signal processing, for example, U.S. Patent Application Publication No. 2009 / 0089048 describes an automated method for determining local activation time (LAT) from a multi-channel ECG signal with four or more channels, including a ventricular channel, a mapping channel, and multiple reference channels. Summary of the Invention [Means for solving the problem]

[0005] An embodiment of the present invention described herein provides a system including a signal acquisition circuit and a processor. The signal acquisition circuit is configured to receive a plurality of intracardiac signals acquired by a plurality of electrodes of an intracardiac probe in a patient's heart. The processor is configured to perform a series of annotation visualization operations at subsequent times, each operation including extracting a plurality of annotation values ​​from the intracardiac signals, selecting a group of intracardiac signals, identifying one or more annotation values ​​from the group that statistically deviate beyond a predefined measure of deviation, and omitting the statistically deviating annotation values ​​from visualization while visualizing the annotation values ​​to a user. The processor is further configured to evaluate a proportion of annotation values ​​omitted across one or more of the annotation visualization operations, and take corrective action in response to detecting a proportion of omitted values ​​exceeding a predefined threshold.

[0006] In some embodiments, the processor is configured to take corrective action by reapplying one or more of the omitted annotation values ​​and re-identifying the statistically deviant annotation values. In some embodiments, the processor is configured to take corrective action by re-extracting one or more of the annotation values ​​from the intracardiac signal.

[0007] In one embodiment, the processor is configured to define the measure of deviation using a standard score of the annotation values, hi another embodiment, the processor is configured to define the measure of deviation using one or more percentiles of the annotation values.

[0008] In an exemplary embodiment, for a given annotation visualization operation, the processor is configured to calculate the deviation of annotation values ​​for intracardiac signals acquired by a selected subset of spatially related electrodes located within a predefined distance from one another in the heart. In another embodiment, in calculating the deviation of annotation values ​​for a given annotation visualization operation, the processor is configured to average the intracardiac signals over multiple temporally related cardiac cycles occurring within a predefined time period.

[0009] In yet another embodiment, in a given annotation visualization operation, the processor is configured to modify one or more of the annotation values ​​in a given intracardiac signal acquired by a given electrode in the group to compensate for a displacement of the given electrode relative to other electrodes in the group.

[0010] In some embodiments, the annotation values ​​include local activation times (LAT). In some embodiments, the processor is configured to visualize the annotation values ​​by overlaying the annotation values, other than the statistically deviant annotation values, on a model of the heart.

[0011] According to an embodiment of the present invention, there is further provided a method including receiving a plurality of intracardiac signals acquired by a plurality of electrodes of an intracardiac probe in a patient's heart. At a subsequent time, a series of annotation visualization operations are performed by each operation (i) extracting a plurality of annotation values ​​from the intracardiac signals, (ii) selecting a group of the intracardiac signals, (iii) identifying one or more annotation values ​​within the group that are statistically deviant by more than a predefined measure of deviation, and (iv) visualizing the annotation values ​​to a user while omitting the statistically deviant annotation values ​​to prevent their visualization. An omission rate of annotation values ​​is evaluated during one or more of the annotation visualization operations. In response to detecting that the omission rate exceeds a predefined threshold, corrective action is taken.

[0012] The invention will be more fully understood from consideration of the following detailed description taken in conjunction with the drawings, in which: [Brief description of the drawings]

[0013] [Figure 1] FIG. 1 is a schematic, pictorial illustration of an electro-anatomical system for multi-channel measurement of intracardiac ECG signals, in accordance with an embodiment of the present invention. [Diagram 2] FIG. 2 shows a schematic diagram of signal acquisition by multiple electrodes over multiple cardiac cycles according to an embodiment of the present invention. [Figure 3A] 1 is a flow chart that illustrates a schematic diagram of a first method for increasing the reliability of an annotation value, according to an embodiment of the present invention; [Figure 3B] 4 is a flow chart that illustrates a schematic diagram of a second method for increasing the reliability of an annotation value, according to an embodiment of the present invention; [Figure 4] 1 is a flow chart that generally illustrates an improved method for increasing the reliability of annotation values, according to an embodiment of the present invention; [Figure 5A] FIG. 13 is a schematic diagram illustrating the annotation of a series of LAT values ​​from a group of spatially related electrodes in an electroanatomical system, according to an embodiment of the present invention. [Figure 5B] FIG. 5B illustrates a schematic diagram of a statistical characteristic calculation of the LAT values ​​received by the processor (in FIG. 5A) when no corrective action is taken, according to an embodiment of the present invention. [Figure 5C] FIG. 2 illustrates, in a simplified manner, the avoidance of error propagation by reapplying default values ​​according to an embodiment of the present invention; [Figure 6A] FIG. 2A is a schematic diagram illustrating probabilities of LAT values ​​for a series of signals acquired by spatially related electrodes in an electro-anatomical system, according to an embodiment of the present invention. [Figure 6B] FIG. 2 shows a schematic diagram of annotation of LAT values ​​for group-1 signals according to an embodiment of the present invention. [Figure 6C]FIG. 2 illustrates a schematic diagram of annotation of LAT values ​​for group-2 signals according to an embodiment of the present invention. [Figure 6D] FIG. 2 illustrates a schematic diagram of re-annotation of Group-1 signals according to an embodiment of the present invention. [Figure 6E] FIG. 13 is a schematic diagram illustrating the recalculation of annotation values ​​for group-2 after re-annotation according to an embodiment of the present invention. [Figure 7] 1 is a flow chart that illustrates generally a method for avoiding error propagation through re-annotation according to an embodiment of the present invention; DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

[0014] Overview Intracardiac probe-based (e.g., catheter-based) cardiac diagnostic and therapeutic systems may measure multiple intracardiac signals, such as an electrocardiogram (ECG), during an invasive procedure. Such systems may acquire multiple intracardiac signals using electrodes attached to the distal end of a probe (hereinafter also referred to as "distal electrodes"). The measured signals may be used to provide a physician with visual information of the heart, such as 3D mapping of the source of pathological electrical patterns within a patient's heart, to aid in corrective medical procedures, such as ablation.

[0015] The measured signals are typically weak and have a low signal-to-noise ratio (SNR). Furthermore, the galvanic connection between some electrodes and the tissue may be poor or non-existent. Meanwhile, because many electrodes are used, there may be some redundancy in the data the system receives from the electrodes.

[0016] Embodiments of the invention disclosed herein provide an intracardiac probe-based electroanatomical measurement and analysis system and method that uses statistical characteristics of the signals collected by the distal electrodes to improve the quality and reliability of the collected data.

[0017] In the following description, reference is made herein to annotation values ​​of local activation time (LAT). However, the disclosed technology is not limited to LAT, and various embodiments of the invention may use annotation values ​​of various other suitable signal parameters.

[0018] In some embodiments according to the invention, the processor extracts the annotation values ​​(e.g., LAT) of the signals and then calculates statistical properties of the LAT values ​​of the signals acquired by the corresponding electrode group (which may include all or a portion of the electrodes). In one embodiment, the statistical properties are the average (e.g.,

[0019]

number

[0020]

number

[0021] In another embodiment, the statistical characteristics include quartiles of the LAT values. The processor calculates a first quartile Q1 and a third quartile Q3, and then ignores all values ​​below Q1 or above Q3 (wherein the first quartile (Q1) is defined as the halfway number between the minimum and median of the data set, and the third quartile (Q3) is the halfway value between the median and maximum of the data set). Alternatively, the processor may define the measure of deviation of the LAT values ​​by any other suitable percentile (or percentiles) of the LAT values. As a further alternative, any other suitable processing method for discarding outlying LAT values ​​can be used.

[0022] The techniques disclosed herein assume the absence of noise and poor galvanic connections, so that the electrodes in each group exhibit similar annotation values. Typically, if the electrodes are far apart, the annotation values ​​obtained by those electrodes may be substantially different. Moreover, the signals from each electrode may be periodically annotated on a beat-by-beat basis ("cardiac cycle"), and if the cardiac cycles are far apart in time, the annotation values ​​derived from those cardiac cycles may be different. In one embodiment, the signals within a group are interrelated. In some embodiments, a tracking system measures the geometrical location of the electrodes, and the groups include annotation values ​​derived only from electrodes that are adjacent to each other ("spatially related", i.e., electrodes located within a predetermined distance of each other). In other embodiments, the group includes annotation values ​​derived only from cardiac cycles that are adjacent to each other ("temporally related", i.e., cardiac cycles that all occur within a predefined continuous time period), and in one embodiment, the group includes values ​​that are related to each other in both space and time (abbreviated as "related values").

[0023] In some embodiments, the processor calculates statistical properties of a group of related LAT values ​​and then discards LAT values ​​within the group that are statistically outliers, e.g., substantially different from the average value of the values ​​within the group (the remaining group of LAT values ​​is referred to as the group of valid LAT values). Thus, LAT values ​​corresponding to electrodes with poor connections or that experience excessive noise may be excluded from the group of valid LAT values.

[0024] In an embodiment, to determine whether a LAT value is statistically deviant from the average LAT of the group of signals, the processor measures the deviation of the annotated LAT value from the average of the group of LAT values. In one embodiment, the measure of deviation is the standard score of the value (defined as the difference between the value and the average divided by the standard deviation), which is compared to a pre-set limit value. For example, values ​​more than 3.5 standard deviations above the average (standard score = 3.5) or more than 1.5 standard deviations below the average (standard score = -1.5) may be considered statistically deviant and therefore may be omitted. In another embodiment, the processor omits values ​​below the first quartile or above the third quartile.

[0025] In some embodiments of the present invention, the processor may mitigate discrepancies between LAT values ​​of spatially related electrodes due to different time delays of cardiac signal propagation within the heart. According to embodiments, the processor may modify the LAT annotation acquired by a given electrode by compensating for the displacement of the given electrode relative to other electrodes to account for differences in propagation delays.

[0026] An ECG signal may sometimes be ambiguously interpreted, and two LAT values ​​may be annotated with different probabilities (a correct value and an incorrect value). For example, a signal may have two peaks close to each other. As a result, in some cases, the calculated mean value of a subset may be closer to the incorrect LAT value than to the correct LAT value.

[0027] In some embodiments, the processor repeats the selection of groups over time (e.g., once every 32 cardiac cycles), with partial overlap between subsequent groups. The processor forms new groups by removing some (but not all) of the old values ​​(e.g., values ​​extracted in the oldest 16 cardiac cycles) and adding new values ​​(e.g., values ​​extracted in the most recent 16 cardiac cycles) (in this context, "old" and "new" refer to the consecutive numbers of cardiac cycles in which the signals were acquired). If the processor calculates an incorrect average of the first group (e.g., due to signal ambiguity), the processor may omit the overlapping correct LAT values ​​from the next group and discontinue the average calculation, so that the newly annotated LAT values ​​will be correct, although the error in the first average calculation may propagate to the subsequent groups.

[0028] The inventive embodiments disclosed herein avoid such error propagation. In one embodiment, the processor monitors the number of omitted LAT values, and depending on the proportion of values ​​omitted, the processor may decide to consider the omitted values ​​in the calculation of the new average, so that the error does not propagate beyond the erroneous group. In another embodiment, if the processor detects a large number of omitted values, the processor may re-annotate some of the values ​​by alternatively interpreting the corresponding signal, and then count the omitted values ​​again. The processor selects the alternative annotation if the number of omitted values ​​decreases.

[0029] In summary, a processor according to embodiments of the present invention may improve the quality and reliability of a set of annotation values ​​of spatially and / or temporally related intercardiac signals by calculating statistical properties of the annotation values, comparing the annotation values ​​to the mean value of the set, and excluding values ​​far from the mean value from the set of valid values. In some embodiments, before calculating the statistical properties, the processor may modify the set of annotation values ​​to correct for signal propagation delays. To avoid error propagation, the processor may monitor the number of omitted values ​​and, depending on the percentage of omitted values, re-evaluate the set with omitted values ​​restored. In other embodiments, depending on the percentage of omitted values, the processor may re-annotate the values ​​and look for alternative LAT interpretations.

[0030] System Description 1 is a schematic, pictorial illustration of an electroanatomical system 21 for multi-channel measurement of intracardiac ECG signals, in accordance with one embodiment of the present invention. In some embodiments, the system 21 is used for electroanatomical mapping of the heart.

[0031] Figure 1 depicts a physician 22 using an electroanatomical catheter 23 to perform electroanatomical mapping of a heart 24 of a patient 25. The catheter 23 includes one or more arms 26, which may be mechanically flexible, at its distal end, with one or more distal electrodes 27 coupled to each of the arms. As will be appreciated, although Figure 1 shows a catheter having five arms, other types of catheters may be used in alternative embodiments according to the invention. The electrodes are coupled to a processor 34 via an interface 32.

[0032] During the electroanatomical mapping procedure, a tracking system may be used to track the intracardiac location of the distal electrodes 27, thereby associating each of the acquired electrophysiological signals with a known intracardiac location. One example of a tracking system is Active Current Location (ACL), described in U.S. Pat. No. 8,456,182. In the ACL system, a processor estimates the location of each of the distal electrodes 27 based on the impedance measured between each of the distal electrodes 27 and a number of surface electrodes 28 coupled to the skin of the patient 25 (for ease of illustration, only one surface electrode is shown in FIG. 1). The processor may then associate any electrophysiological signals received from the distal electrodes 27 with the location at which the signal was acquired.

[0033] In some embodiments, the plurality of distal electrodes 27 acquire intracardiac ECG signals from atrioventricular tissue of the heart 24. The processor comprises signal acquisition circuitry 36 coupled to receive the intracardiac signals from the interface 32, a memory 38 for storing data and / or instructions, and a processing unit 42 (e.g., a CPU or other processor).

[0034] The signal acquisition circuitry 36 digitizes the intracardiac signals to generate a plurality of digital signals. The acquisition circuitry then communicates the digitized signals to a processing unit 42 included in the processor 34.

[0035] Processing unit 42 is configured, among other tasks, to extract annotation parameters from the signals and to calculate statistical properties, such as the average of the annotated parameters of adjacent signals that may be similar to each other (note that in this context adjacent signals refer to signals from electrodes located close to each other ("spatially related") and / or annotation values ​​extracted from cardiac cycles that are close to each other in time ("temporally related").

[0036] After calculating the statistical properties, the processing unit is further configured to filter out (i.e., omit) from the set annotation values ​​that are likely to be invalid (e.g., annotation values ​​from electrodes with poor galvanic connections or subject to strong temporal noise, etc.) The remaining annotation values ​​shall hereafter be referred to as "valid annotation values."

[0037] Processor 34 visualizes to the user the valid annotation values, i.e., annotation values ​​excluding omitted statistically inconsistent annotation values. In some embodiments, processor 34 visualizes the valid annotation values, for example, by overlaying the valid annotation values ​​on an electroanatomical map 50 of the heart and displaying the map on screen 52 to physician 22. Alternatively, processor 34 may visualize the valid annotation values ​​(after omitting invalid annotation values) in any other suitable manner.

[0038] The exemplary diagram shown in Figure 1 has been chosen purely for conceptual clarity. For example, in alternative embodiments of the invention, position measurements may be made by applying a voltage gradient between pairs of surface electrodes 28 and measuring the resulting potential with distal electrodes 27 (i.e., using CARTO® 4 technology manufactured by Biosense-Webster, Inc., Irvine, Calif.). Thus, embodiments of the invention apply to any position sensing method.

[0039] Other types of catheters, such as a Lasso® Catheter (manufactured by Biosense-Webster) or a basket catheter, may equally be used. A contact sensor may be attached to the distal end of the electroanatomical catheter 23. Other types of electrodes, such as those used for ablation, may also be used similarly to electrode 27 to acquire electrophysiological signals within the heart.

[0040] 1 primarily illustrates the portions of the system that are relevant to an embodiment of the present invention, omitting other system elements such as external ECG recording electrodes and their connections, as well as various ECG recording system elements such as filtering, digitization, circuit protection, and other elements.

[0041] In any embodiment, a readout application specific integrated circuit (ASIC) is used to measure the intracardiac ECG signal. Various elements of the signal acquisition circuitry 36 routing may be implemented in hardware using, for example, one or more discrete components, a field programmable gate array (FPGA) or an ASIC. In some embodiments, some elements of the signal acquisition circuitry 36 and / or processing unit 42 may be implemented in software or using a combination of software and hardware elements.

[0042] Processing unit 42 typically comprises a general-purpose processor having software programmed to perform the functions described herein. The software may be downloaded in electronic form, for example over a network, or it may alternatively or additionally be provided on and / or stored on a non-transitory, tangible medium, such as, for example, a magnetic, optical, or electronic memory.

[0043] Associated annotation value Associated annotation values ​​are derived from spatially related electrodes (e.g., electrodes that are geometrically close to each other, i.e., located no more than a predefined distance away from each other) and / or from temporally related signals (e.g., values ​​extracted from cardiac cycles that are close to each other, i.e., values ​​occurring within a predefined duration). More precisely, associated annotation values ​​are annotation values ​​whose combined distance, including the geometric distance between electrodes and the temporal distance between cardiac cycles, is below some predefined threshold.

[0044] Fig. 2 is a diagram 200 that shows a schematic of signal acquisition by multiple electrodes over multiple cardiac cycles. The horizontal axis 202 shows the cardiac cycle (each vertical line is one cardiac cycle) and the vertical axis 204 shows the distance of the electrode from a reference point (only one spatial dimension is shown). (As will be appreciated, two or three dimensions may actually be used, but are not shown for clarity). According to the exemplary embodiment shown in Fig. 2, there are electrodes on all horizontal lines and LAT annotation values ​​are registered for all intersections of horizontal and vertical lines (each intersection shall hereafter be referred to as a LAT point).

[0045] Curve 206 is an iso-LAT line showing the location of the indicated LAT value, and the electrodes likely measure values ​​interpolated from adjacent iso-LAT curves during the corresponding cardiac cycle. For example, the expected registration value of LAT point 208 (which is vertically midway between iso-LAT line 714 and iso-LAT line 716) is 715, while the expected registration value of LAT point 210 is 708.5.

[0046] As can be seen, adjacent vertical lines have similar LAT values, and adjacent horizontal lines have similar LAT values. Circles 212 represent groups of related LAT values ​​214 that are close to each other in terms of geometric (vertical) and temporal (horizontal) distance.

[0047] The exemplary diagram shown in Figure 2 is simplified for purposes of conceptual clarity only - for example, in alternative embodiments, the distance between electrodes may not be uniform and the constellation of associated signals may not be circular.

[0048] 3A is a flow chart 300 that generally illustrates a first method for increasing the reliability of annotation values, according to an embodiment of the present invention. The flow is executed by processor 34 (FIG. 1). The flow begins with a signal recording step 302, in which the processor records ECG signals monitored by electrodes 27 and acquired by acquisition circuitry 36 (FIG. 1). Next, in an annotation value extraction step 304, the processor calculates annotation values ​​for each electrode and each cardiac cycle.

[0049] The processor then proceeds to an electrode position acquisition step 306, where the positions of the electrodes are acquired (e.g., using ACL techniques) and the spatial location of each electrode is registered. The processor then proceeds to a group selection step 308.

[0050] In step 308, the processor selects groups of related annotation values. As previously explained, groups include annotation values ​​from spatially and / or temporally related signals that are likely to be similar to each other.

[0051] Next, the processor calculates the mean and standard deviation of all annotation values ​​in the group, in a mean and standard deviation calculation step 310. Any suitable type of mean value can be used in this context, such as, for example, the arithmetic mean, geometric mean, median, root mean square (RMS) value, centroid, or any other average value.

[0052] The processor then iterates for each annotation value in the set, proceeding sequentially to step 312, step 314, and either step 316 or step 318. In standard score calculation step 312, the processor calculates a standard score for the annotation value (e.g., by dividing the difference between the annotation value and the mean by the standard deviation). In standard score comparison step 314, the processor compares the standard score calculated in step 312 with a pre-set limit value. If the standard score exceeds the pre-set limit value, the process proceeds to value exclusion step 316, where the processor excludes the statistically deviant annotation value. If the standard score is within the pre-set limit value, the process proceeds to value addition step 318, where the processor adds the annotation value to the set of valid annotation values.

[0053] The processor repeats the sequence of steps 312, 314, and either step 316 or step 318 for all annotation values ​​in the group. The flowchart may then repeat the steps (from step 308) for other groups of associated electrodes.

[0054] When the flow is finished, a set of valid annotation values ​​replaces the original set; the new set is more reliable because it omits extreme values ​​(e.g., values ​​from electrodes with poor galvanic connections).

[0055] Figure 3B is a flow chart 350 that illustrates in schematic form a second method for increasing the reliability of annotation values ​​according to an embodiment of the present invention. The method illustrated in Figure 3B differs from the method illustrated in Figure 3A only in the statistical properties and the selection of omitted values. Thus, steps 302-318 illustrated in Figure 3A are identical to steps 352-368, respectively, in Figure 3B, with the exception of steps 360 and 364 in Figure 3B, which are different from steps 310 and 314, respectively, in Figure 3A and are therefore described below.

[0056] In a quartile calculation step 362, the processing unit 42 (FIG. 1) calculates the first quartile (Q1) and third quartile (Q3) of the group of LAT values ​​(where Q1 is defined as the halfway number between the minimum and the median of the group of LAT values, and Q3 is the halfway number between the median and the maximum of the group of LAT values).

[0057] In annotation value comparison step 364, the processing unit compares the annotated LAT value with Q1 and Q3, respectively. If the value is less than Q1 or greater than Q3, the processing unit proceeds to annotation value removal step 366, and if the value is between Q1 and Q3, the processing unit proceeds to value addition step 368.

[0058] The exemplary flow charts shown in Figures 3A and 3B are chosen purely for the purpose of conceptual clarity. For example, in alternative embodiments, annotation values ​​may be extracted as the signals are acquired (rather than after they are recorded). In one embodiment, the selection of the group of signals may be done by a physician, while in other embodiments, the processor selects the group according to a spatial and / or temporal range prescribed by the physician.

[0059] In some embodiments, step 318 (step 368 in FIG. 3B) is not required and the processor removes extreme values ​​from the group in step 316 (366) so that only good values ​​remain when the flow is complete. In other embodiments, all annotation values ​​are initially marked as invalid and step 316 (366) is not required.

[0060] In some embodiments, other statistical properties used are different from those already described. For example, in one embodiment, octiles may be used instead of quartiles, and the processing unit may omit values ​​below the first octile or above the last octile. Furthermore, alternatively, any other suitable percentiles may be used.

[0061] Any other suitable statistical method for detecting and discarding extreme values ​​may be used in alternative embodiments.

[0062] Propagation Delay Compensation In some embodiments, the previously described techniques may be improved by correcting the extracted LAT values ​​prior to statistical characteristic calculations for expected variations in values ​​due to different spatial locations of the electrodes. For example, waves through the heart may be assumed to travel at a given speed (e.g., 1 m / sec). Using the known locations of the electrodes acquiring the signals, theoretical differences in LAT values ​​may be applied in calculating the average value.

[0063] 4 is a flow chart 400 that generally illustrates an improved method for increasing the reliability of annotation values ​​according to an embodiment of the present invention. The flow is executed by processor 34 (FIG. 1). The flow starts with a signal recording step 402, followed by an annotation value calculation step 404, followed by an electrode position acquisition step 406, followed by a group selection step 408, although steps 402, 404, 406, and 408 may be identical to steps 302, 304, 306, and 308 (FIG. 3), respectively.

[0064] The processor then proceeds to a LAT value correction step 410 where, for each LAT value in the group, the processor calculates and applies an estimated correction value according to the spatial location of the electrodes and the assumed wave travel speed. After step 410, flow returns to Figure 3 at mean and standard deviation calculation step 310.

[0065] In this way, estimates of deviations caused by propagation delays can be removed from the group, further improving the reliability of the annotation signal.

[0066] The exemplary flow chart shown in Figure 4 has been chosen purely for the purposes of conceptual clarity. For example, in alternative embodiments, the correction of expected signal delays can be integrated into the mean and standard deviation calculation steps. In other embodiments, the correction is performed before the groups are selected (so that groups may contain more relevant LAT values).

[0067] Error Propagation and Prevention In some embodiments according to the invention, processor 34 successively selects groups of related LAT values ​​in response to receiving signals from successive heartbeats. In some cases, the processor may select groups of values ​​corresponding to signals that overlap in time. For example, if a first group includes LAT values ​​extracted from signals corresponding to cardiac cycles n through m, and a second group includes values ​​corresponding to cycles x through y, the groups will overlap if x is greater than m (x>m) and y is greater than n (y>n).

[0068] In the preceding description, reference is made to LAT values ​​annotated from different cardiac cycles with reference to Figures 5A-5C, 6A-6E, and 7. Each of the LAT values ​​may in fact include multiple spatially related LAT values ​​(including LAT values ​​corrected for signal propagation delays). However, for simplicity, one LAT value per cardiac cycle is shown.

[0069] 5A-5C are timing diagrams illustrating events along a horizontal time axis (not shown). Fig. 5A is a diagram 500 that generally illustrates the annotation of a series of LAT values ​​in an electroanatomical system, according to an embodiment of the present invention. In the exemplary embodiment shown in Fig. 5A, each group includes 10 LAT values ​​corresponding to 10 consecutive cardiac cycles, and each pair of consecutive groups shares five common LAT values.

[0070] For clarity, in the example shown in Figure 5A, the annotated LAT values ​​can be either close to "X" or close to "Y" (for short, we will refer to them as X and Y values). As will be understood, this simplification is not meant to limit the scope of the present disclosure in any way. Embodiments of the present invention may annotate any suitable set of LAT values.

[0071] In the example shown in Figure 5A, the first five LAT values ​​are "X", but when a subsequent intracardiac signal is received, the majority of the LAT values ​​are "Y", meaning that the initial annotation and decision taken for the first group was probably incorrect. As will be explained below (with reference to Figure 5B), in some embodiments according to the invention, no corrective action is taken and the incorrect decision taken for the first group may propagate to further groups. As will be explained further (with reference to Figure 5C), in other embodiments, corrective action is taken to prevent error propagation.

[0072] 5B is a diagram 510 that generally illustrates the computation of statistical characteristics of LAT values ​​received by a processor (in FIG. 5A) when no corrective action is taken, in accordance with an embodiment of the present invention. Diagram 510 includes an averaging scheme for group-1 512, an averaging scheme for group-2 514, and an averaging scheme for group-3 516 (as shown, the averaging for group-4 is identical to groups-2 and-3 and is not shown).

[0073] Group-1 averaging scheme 512 represents the averaging of the first 10 LAT values. Since there are 7 X values ​​and 3 Y values, the average value will be close to X. Therefore, the Y values ​​of the last 5 samples will be omitted from the subsequent groups.

[0074] Group-2 averaging scheme 514, which represents the averaging of LAT values ​​6-15, includes four X values ​​and three Y values ​​(with the two original Y values ​​omitted). This average is still close to X. Therefore, the Y values ​​of the last five samples are omitted from the subsequent group. As can be seen, the error from the first group has propagated to the second group.

[0075] Group-3 averaging scheme 516, representing the averaging of LAT values ​​11-20, also contains four X and three Y values ​​because two original Y values ​​were omitted, and again the average value is close to X, leading to the omission of a Y value. Thus, the error from the first group will propagate to this group and continue to further groups. (According to the example shown in FIG. 5A, 5B, the propagation will stop only if the majority (e.g., 80%) of the subsequent LAT values ​​are Y.)

[0076] 5C is a diagram 520 that illustrates, in schematic form, avoiding error propagation by reapplying omitted values ​​according to an embodiment of the present invention. Diagram 520 includes an averaging scheme 522 for group-1, a preliminary averaging scheme 524 for group-2, a re-evaluated averaging scheme 526 for group-2, and an averaging scheme 528 for group-3.

[0077] Group-1 averaging scheme 522 represents the averaging of the first 10 LAT values ​​and is identical to group-1 averaging scheme 512 (FIG. 5B). Again, the average value will be close to X, and Y values ​​will be omitted from the next group.

[0078] The preliminary averaging scheme 524 for group-2 is the same as the averaging scheme 514 for group-2 (FIG. 5B), and again the average value will be close to X, and Y values ​​will be omitted from the next group. However, according to the exemplary embodiment shown in FIG. 5C, the frequency with which LAT values ​​are omitted is monitored. Now, since there are 6 omitted values ​​out of the 10 sample groups, the processor decides to reapply the omitted values ​​in the statistical characteristic calculation. The decision to reapply the omitted values ​​(e.g., group 2 in FIG. 5C designated by reference numeral 524) is the "corrective action" to be taken.

[0079] The recalculated averaging scheme 526 for group-2 includes the omitted Y value in the statistical calculations. Now, in group-2 there are 6 Y values ​​and 4 X values, the average value is close to Y, and the X value will be omitted from the next group. Now, the averaging scheme 528 for group-3 includes 2 X values ​​and 6 Y values, the average value is close to Y, and the X value will be omitted. This pattern will continue for group-4.

[0080] Thus, according to the exemplary embodiment shown in FIG. 5C, the processor can stop the propagation of errors that may result from the omitted values ​​by monitoring how often values ​​are omitted, and, depending on how often values ​​are omitted, can take corrective action, which may include recalculating the average without the omitted values.

[0081] Other embodiments according to the invention can prevent error propagation by re-annotating at least some of the LAT values ​​depending on how often values ​​are omitted, which can be useful when intracardiac signals are ambiguous and can be interpreted in more than one way.

[0082] 6A-6E are timing diagrams illustrating events along a horizontal time axis (not shown).

[0083] FIG. 6A is a schematic diagram illustrating the probability of LAT values ​​for a set of signals acquired by spatially related electrodes in an electroanatomical system according to an embodiment of the present invention. Each rectangle in FIG. 6A represents the probability that the LAT value of the corresponding signal is X (top number) or Y (bottom number). The first five signals have a 60% probability of X value as the LAT value and a 40% probability of Y value. Later signals (signals related to later cardiac cycles) are mostly more likely to be Y value than X value. The signals are divided into partially overlapping groups (1-3), each with 10 samples and with an overlap of 5 samples between subsequent groups.

[0084] 6B is a schematic diagram illustrating annotation of LAT values ​​for signals in group-1 according to an embodiment of the present invention. The processor annotates X values ​​for six of the samples and Y values ​​for four. The average value is close to X, and the processor omits four Y values ​​(three of the omitted signals, designated by the dotted rectangle, overlap with group-2).

[0085] 6C is a schematic diagram illustrating annotation of LAT values ​​for signals in group-2 according to an embodiment of the present invention. From the five signals common to group-1, the processor has omitted three signals with higher y probability. Of the remaining seven values, four are X values ​​and three are Y values. Again, the average value is close to X, and the processor may omit Y values ​​from the next group. However, the processor now determines that there are many omitted values ​​(e.g., above a pre-set limit) and attempts to re-annotate the LAT values ​​to reduce the number of omitted values.

[0086] 6D is a schematic diagram illustrating re-annotation of group-1 signals according to an embodiment of the present invention, where the processor annotates the second most probable value for the first five samples (Y, Y, Y, X, Y for signals 1-5, respectively). The annotated LAT values ​​for the first group will now include three X values ​​and seven Y values, with the average value close to Y. The seventh and ninth values ​​will now be omitted from the next group.

[0087] 6E is a schematic diagram illustrating the re-annotation of values ​​of group-2 following the re-annotation of group-1 according to an embodiment of the present invention. Since two X values ​​from five overlapping signals from group-1 are not omitted, group-2 contains six Y values ​​and two X values, where the average value is close to Y, preventing error propagation.

[0088] The error propagation prevention methods described with reference to Figures 5A-5C and 6A-6E are exemplary methods cited merely for conceptual clarity. Embodiments in accordance with the invention may employ alternative methods for error propagation prevention. For example, the group length and overlap length may be any suitable number, including group sizes and / or overlaps that vary over time. In some embodiments, the threshold may be dynamic. In one embodiment, the processor may not annotate the LAT value with the second highest probability if its probability is below a pre-set limit or if the probability is below a limit that is a function of how often LAT values ​​are omitted.

[0089] Finally, in some embodiments, a combination of restoring omitted values ​​and re-evaluating annotations may be used.

[0090] 7 is a flow chart 700 that generally illustrates a method for avoiding error propagation through re-annotation, according to an embodiment of the present invention. The flow chart is executed by processor 34 (FIG. 1).

[0091] The flow begins with receiving overlap annotations step 702, where a processor retrieves from memory annotations of the intracardiac ECG signal that are common to the previous and current groups (as previously mentioned, the processor may omit values ​​that deviate substantially from the average value of the previous group).

[0092] Next, in receive new ECG signals step 704, the processor receives ECG signals from the electrodes (via acquisition circuitry 36 (FIG. 1)) relating to the subsequent cardiac cycle and annotates the corresponding LAT values. The processor then proceeds to calculate statistics step 706, where it calculates statistical properties (e.g., mean and standard deviation) of the group of values, including duplicate values ​​from the previous group and the newly annotated values.

[0093] The processor then moves to a group scan step 708, in which the processor compares each of the group's LAT values ​​to the group mean (calculated in step 706) and omits any values ​​that deviate from the mean by more than a pre-set threshold (or alternatively, values ​​where the standard score is not within pre-set limits). The processor also counts the number of omitted LAT values.

[0094] The processor then proceeds to a counter comparison step 710, where the processor compares the number of omitted values ​​to a preset limit value. If the number of omitted values ​​does not exceed the limit value, the processor proceeds to the next group (e.g., returns to step 702).

[0095] If, at step 710, the number of omitted values ​​exceeds a preset limit, the processor assumes that error propagation occurred as a result of a previous incorrect annotation of the LAT value. The processor then proceeds to receive duplicate signal step 712, where it retrieves the duplicate signal from memory. The processor then proceeds to extract alternative values ​​step 714, where it re-annotates the LAT value from the signal. Since the processor has already annotated the LAT value (at step 704 in the previous group), the processor now looks for an alternative annotation (e.g., the second most probable LAT value for each signal).

[0096] After step 714, the processor transitions back to step 706 and recalculates the statistics using the alternative annotation values. If the number of omitted values ​​is now less than the previous count in step 708, the new values ​​will take effect.

[0097] In summary, according to the exemplary embodiment shown in Figure 7, the LAT values ​​of successive groups of related intracardiac signals are annotated. If an error occurs in the first group, the processor may detect that values ​​in subsequent groups are frequently omitted. The processor then tries alternative annotations of the signals and selects the LAT value corresponding to the alternative annotation if it results in fewer omitted annotations.

[0098] The exemplary flowchart shown in FIG. 7 is chosen purely for conceptual clarity. For example, in an alternative embodiment, the processor extracts a first set of most probable annotation values ​​and a second set of alternative annotation values, so that step 714 is replaced by an alternative value selection step. In some embodiments, more than two LAT value candidates are extracted, and the flowchart is modified accordingly. In one embodiment, step 720, in which the number of omitted values ​​is compared to a threshold, may be replaced by other means of comparing frequencies rather than counts. In an embodiment, instead of extracting alternative annotations, the processor recalculates statistical properties of the group using the omitted values. In yet another embodiment, the processor attempts both statistical recalculation and alternative annotation values.

[0099] It will be understood that the above embodiments are given by way of example, and the present invention is not limited to what has been particularly shown and described above. Rather, the scope of the present invention includes both combinations and subcombinations of the various features described above, as well as variations and modifications thereof not disclosed in the prior art that would occur to one skilled in the art upon reading the above description. Documents incorporated by reference into this patent application are to be considered as part of this application, except that if any term is defined in such incorporated documents in a manner that is inconsistent with the definition expressly or impliedly given herein, then only the definition in this specification shall be considered.

[0100] [Embodiment] (1) a signal acquisition circuit configured to receive a plurality of intracardiac signals acquired by a plurality of electrodes of an intracardiac probe within a patient's heart; 1. A processor comprising: At a subsequent point in time, a series of annotation visualization operations are performed, with each operation (i) extracting a plurality of annotation values ​​from the intracardiac signal; (ii) selecting said group of intracardiac signals; (iii) identifying one or more annotation values ​​within the set that are statistically deviant from a predefined measure of deviance; and (iv) visualizing the annotation values ​​to a user while omitting the statistically deviant annotation values ​​to prevent them from being visualized; By carrying out assessing an omission rate of annotation values ​​across one or more of the annotation visualization operations; taking corrective action in response to detecting that the omission rate exceeds a predefined threshold a processor configured to: A system comprising: (2) The processor: reapplying one or more of said omitted annotation values; re-identifying the statistically deviant annotation values; The system of claim 1, wherein the system is configured to take the corrective action by: (3) The system of embodiment 1, wherein the processor is configured to perform the corrective action by re-extracting one or more of the annotation values ​​from the intracardiac signal. (4) The system of embodiment 1, wherein the processor is configured to define the measure of deviation using a standard score of the annotation value. (5) The system of embodiment 1, wherein the processor is configured to define the measure of deviation using one or more percentiles of the annotation values.

[0101] (6) The system of embodiment 1, wherein, for a given annotation visualization operation, the processor is configured to calculate a deviation of the annotation value relative to intracardiac signals acquired by a selected subset of spatially related electrodes located within a predetermined distance from each other within the heart. (7) The system of embodiment 1, wherein when calculating the deviation of the annotation value for a given annotation visualization operation, the processor is configured to average the intracardiac signal over multiple temporally related cardiac cycles occurring within a predefined time period. (8) The system of embodiment 1, wherein in a given annotation visualization operation, the processor is configured to modify one or more of the annotation values ​​in a given intracardiac signal acquired by a given electrode in the group to compensate for a displacement of the given electrode relative to other electrodes in the group. (9) The system of embodiment 1, wherein the annotation value includes local excitation time (LAT). (10) The system of embodiment 1, wherein the processor is configured to visualize the annotation values ​​by overlaying the annotation values ​​other than the statistically deviant annotation values ​​on a model of the heart.

[0102] (11) receiving a plurality of intracardiac signals acquired by a plurality of electrodes of an intracardiac probe within the patient's heart; At a subsequent point in time, a series of annotation visualization operations are performed, with each operation (i) extracting a plurality of annotation values ​​from the intracardiac signal; (ii) selecting said group of intracardiac signals; (iii) identifying one or more annotation values ​​within the set that are statistically deviant from a predefined measure of deviance; and (iv) visualizing the annotation values ​​to a user while omitting the statistically deviant annotation values ​​to prevent them from being visualized; and assessing an omission rate of annotation values ​​across one or more of the annotation visualization operations; and taking corrective action in response to detecting that the skip rate exceeds a predefined threshold; The method includes: (12) Taking such corrective action is reapplying one or more of the omitted annotation values; and re-identifying the statistically deviant annotation values; and 12. The method of embodiment 11, comprising: (13) The method of embodiment 11, wherein taking the corrective action includes re-extracting one or more of the annotation values ​​from the intracardiac signal. (14) The method of embodiment 11, wherein the measure of deviation is defined using a standard score of the annotation values. (15) The method of embodiment 11, wherein the measure of deviation is defined using one or more percentiles of the annotation values.

[0103] (16) The method of embodiment 11, comprising, for a given annotation visualization operation, calculating deviation of the annotation value relative to intracardiac signals acquired by a selected subset of spatially related electrodes located within a predefined distance from one another within the heart. (17) The method of embodiment 11, further comprising averaging the intracardiac signal over multiple temporally related cardiac cycles occurring within a predefined time period when calculating the deviation of the annotation value for a given annotation visualization operation. (18) The method of embodiment 11, comprising, in a given annotation visualization operation, modifying one or more of the annotation values ​​in a given intracardiac signal acquired by a given electrode in the group to compensate for a displacement of the given electrode relative to other electrodes in the group. (19) The method of embodiment 11, wherein the annotation value comprises a local excitation time (LAT). (20) The method of embodiment 11, wherein visualizing the annotation values ​​includes superimposing the annotation values ​​other than the statistically deviant annotation values ​​on the model of the heart.

Claims

1. a signal acquisition circuit configured to receive a plurality of intracardiac signals acquired by a plurality of electrodes of an intracardiac probe within the patient's heart; 1. A processor comprising: At a subsequent point in time, a series of annotation visualization operations are performed, with each operation (i) extracting a plurality of annotation values ​​from the intracardiac signal; (ii) selecting said group of intracardiac signals; (iii) identifying one or more annotation values ​​within the set that are statistically deviant by more than a predefined measure of deviance; (iv) visualizing the annotation values ​​to a user while omitting the statistically deviant annotation values ​​to prevent them from being visualized; By carrying out assessing an omission rate of annotation values ​​across one or more of the annotation visualization operations; In response to detecting that the omission rate exceeds a predetermined threshold, re-identifying the statistically deviant annotation values ​​using one or more of the non-abbreviated annotation values ​​and the abbreviated annotation values. take corrective action by a processor configured to: A system comprising:

2. The system of claim 1 , wherein the processor is configured to define the measure of deviation using a standard score of the annotation values.

3. The system of claim 1 , wherein the processor is configured to define the measure of deviation using one or more percentiles of the annotation values.

4. 2. The system of claim 1, wherein, for a given annotation visualization operation, the processor is configured to calculate a deviation of the annotation value relative to intracardiac signals acquired by a selected subset of spatially related electrodes located within a predefined distance from one another within the heart.

5. The system of claim 1 , wherein when calculating the deviation of the annotation value for a given annotation visualization operation, the processor is configured to average the intracardiac signal over multiple cardiac cycles that occur within a predefined time period and are close in time to one another.

6. 2. The system of claim 1, wherein in a given annotation visualization operation, the processor is configured to modify one or more of the annotation values ​​in a given intracardiac signal acquired by a given electrode in the group to compensate for a displacement of the given electrode relative to other electrodes in the group.

7. The system of claim 1 , wherein the annotation value comprises a local excitation time (LAT).

8. 2. The system of claim 1, wherein the processor is configured to visualize the annotation values ​​by overlaying the annotation values ​​other than the statistically deviant annotation values ​​on the model of the heart.

9. A method of operating a processor, comprising: receiving a plurality of intracardiac signals acquired by a plurality of electrodes of an intracardiac probe within a patient's heart; At a subsequent point in time, a series of annotation visualization operations are performed, with each operation (i) the processor extracting a plurality of annotation values ​​from the intracardiac signal; (ii) the processor selecting the group of intracardiac signals; and (iii) the processor identifying one or more annotation values ​​within the group that are statistically deviant by more than a predetermined measure of deviance; (iv) visualizing the annotation values ​​to a user while the processor omits and prevents the statistically deviant annotation values ​​from being visualized; and evaluating, by the processor, a rate of omission of annotation values ​​over one or more of the annotation visualization operations; In response to detecting that the omission rate exceeds a predetermined threshold, the processor re-identifying the statistically deviant annotation values ​​using one or more of the non-abbreviated annotation values ​​and the omitted annotation values; and taking corrective action by The method of operation includes:

10. The method of claim 9 , wherein the measure of deviation is defined using a standard score of the annotation values.

11. The method of claim 9 , wherein the measure of deviation is defined using one or more percentiles of the annotation values.

12. 10. The method of claim 9, further comprising: calculating, for a given annotation visualization operation, deviations of the annotation values ​​for intracardiac signals acquired by a selected subset of spatially associated electrodes located within a predefined distance from one another within the heart.

13. The method of claim 9, further comprising the processor averaging the intracardiac signal over multiple cardiac cycles occurring within a predefined time period and in close temporal proximity to one another when calculating the deviation of the annotation value for a given annotation visualization operation.

14. 10. The method of claim 9, further comprising: in a given annotation visualization operation, the processor modifies one or more of the annotation values ​​in a given intracardiac signal acquired by a given electrode in the group to compensate for a displacement of the given electrode relative to other electrodes in the group.

15. The method of claim 9 , wherein the annotation value comprises a local activation time (LAT).

16. The operating method of claim 9, wherein the processor visualizing the annotation values ​​includes the processor overlaying the annotation values ​​other than the statistically deviant annotation values ​​onto a model of the heart.

Images