Abstract:
Devices, systems, and methods for localizing and/or treating arrhythmias of a patient&#39;s heart, which are particularly useful for localizing focal atrial fibrillation, allow locating arrhythmogenic regions of a chamber of the heart using heart cycle signals measured from a body surface of the patient. Non-invasive localization of the ectopic origin or exit site allows focal, circular, and/or perimeter treatment to be directed so as to inhibit complex arrhythmias without having to rely on wide-spread and time consuming sequential searches and/or on massively invasive simultaneous electrocardial sensors. The invention recognizes that effective localization of these complex arrhythmias can be significantly enhanced by techniques and structures which separate heart cycle signals originating from differing chambers and/or regions of the heart tissue.

Description:
CROSS-REFERENCES TO RELATED APPLICATIONS  
       [0001]     This patent application is a continuation of U.S. patent application Ser. No. 09/808,728 filed on Mar. 14, 2001, which is a continuation-in-part of U.S. patent application Ser. No. 09/611,179 filed on Jul. 6, 2000, which claims benefit under 35 USC 119(e) from U.S. Provisional Patent Application No. 60/200,965 filed on May 1, 2000; and U.S. Provisional Patent Application Nos. 60/189,610 and 60/189,611, both filed on Mar. 15, 2000. The subject matter of this application is related to that of concurrently filed applications entitled:  Continuous Localization and Guided Treatment of Cardiac Arrhythmias,  co-assigned with the present application, and  QRST Subtraction Using an Adaptive Template for Analysis of T - Wave Obscured Atrial Activity;  and is also related to that of U.S. patent application Ser. No. 09/724,947 filed Nov. 28, 2000, and U.S. Provisional Patent Application Nos. 60/196,204 filed Apr. 11, 2000, and 60/189,513 filed Mar. 15, 2000. The full disclosures of all of these applications are incorporated herein by reference. 
     
    
     STATEMENT AS TO RIGHTS TO INVENTIONS MADE UNDER FEDERALLY SPONSORED RESEARCH OR DEVELOPMENT  
       [0002]     NOT APPLICABLE  
       REFERENCE TO A “SEQUENCE LISTING,” A TABLE, OR A COMPUTER PROGRAM LISTING APPENDIX SUBMITTED ON A COMPACT DISK  
       [0003]     NOT APPLICABLE  
       BACKGROUND OF THE INVENTION  
       [0004]     I. Field of Invention  
         [0005]     The present invention generally relates to devices, systems, and methods for diagnosing and/or treating of the heart. In a particular embodiment, the invention provides techniques for localizing and/or treating atrial fibrillation and other arrhythmias.  
         [0006]     Significant progress has recently been made toward effective treatments of many cardiac arrhythmias. Contraction of a healthy human heart generally propagates through the heart tissue from the sinus node in the right atrium, and eventually the associated ventricles. This normal propagation of contraction forces blood to flow from the atria to the ventricles in a synchronized pumping action. Focal or re-entrant arrhythmias of the heart often originate at, and propagate from alternative heart tissue locations, resulting in irregular contractions of some or all of the heart tissues. Radiofrequency intracardiac catheter ablation of the alternative ectopic origin is now used to effectively treat a variety of arrhythmias, including ventricular tachycardia (VT).  
         [0007]     Although quite effective, current catheter ablation of arrhythmogenic sites has significant disadvantages. A particular challenge in an effective catheter ablation treatment is the time required for proper identification of the treatment site. As it is generally desirable to limit the size of the ablation, significant time is often spent testing candidate ablation sites. These candidate sites are often tested sequentially by positioning the intracardiac catheter against a site within (for example) the right ventricle, identifying the engaged tissue location within the ventricle, sensing and/or pacing the heart at the candidate site, repositioning the intracardiac catheter to a new candidate site, and repeating this process until the ectopic origin has been identified.  
         [0008]     As fluoroscopy is often used to identify the location of the engaged tissue, this sequential iterative process can result in significant exposure of the patient and treating personnel to potentially harmful radiation. While alternative (and more complex) intracardiac catheter probe structures have been proposed to allow more rapid identification of the ectopic origin(s) of VTs and other focal arrhythmias, the size and cost of these complex structures may limit their acceptability, particularly for treatment of (for example) the left atrium of the heart, which is often accessed from the right atrium by a puncture through the atrial septum.  
         [0009]     To overcome the disadvantages associated with these known time consuming and/or invasive intracardiac arrhythmia sensing and localization techniques, researchers have been working on alternative arrhythmia localization techniques which rely on body surfacing mapping often during pacing. Pacing often comprises initiating the arrhythmia by applying a small electrical pulse from a catheter. Electrocardiograms (ECG) may be recorded during abnormal atrial or ventricular activity and compared with ECGs taken during pacing at different sites within the heart, optionally using a standard 12-lead ECG system. More detailed information regarding ectopic sites can be obtained by recording heart cycle signals at the body surface using a more comprehensive sensor array. These heart cycle signals, which generally comprise small amplitude variations in electrical potential along the anterior and/or posterior torso, can be manipulated and/or mapped so as to provide an indication of the origin of the arrhythmia within the heart. Much of this work has concentrated on VT. More recent work has begun to investigate the possibility of localizing certain atrial arrhythmias, such as right atrial tachycardia. While the initial results of this research appear quite promising for treatment of selected individuals, significant advancements would be beneficial to allow widespread treatment of patients suffering from cardiac arrhythmias.  
         [0010]     The most common form of cardiac arrhythmia may be atrial fibrillation (AFib). Atrial fibrillation is often paroxysmal in nature, which may contribute to the significant risks of the disorder. Atrial fibrillation may result in twice as many hospitalizations annually as VT, and may cause significant morbidity and/or mortality, leading not only to heart failure, but associated risks of thrombo-embolism and stroke.  
         [0011]     Many current AFib patients are managed using antiarrhythmic drugs. Unfortunately, existing drug treatments are merely palliative, since they are aimed at suppression of the arrhythmia and not at curing the underlying disease. Many researchers are directing resources to development of therapeutic catheters to treat atrial fibrillation, attempting to build on several years of successful ablation for treatment of other arrhythmias. Early indications are that when accurately identified, ablation of ectopic origins of focal AFib may provide an effective treatment for the disorder. Hence, there would be significant benefits to extending the new body surface localization techniques to atrial fibrillation. The nature of AFib, however, represents a significant barrier to the direct application of known mapping techniques used with other arrhythmias.  
         [0012]     Atrial fibrillation is generally more complex and difficult to localize than other arrhythmias. Focal AFib often exhibits an infrequent, irregular occurrence, and may be difficult to induce with known catheter mapping techniques. Even when atrial fibrillation is ongoing and/or successfully induced in the lab, AFib may exhibit prolonged occurrences in many patients, possibly requiring repetitive direct current shock cardioversion to convert the patient back into a normal sinus rhythm. Atrial fibrillation may also have multiple focal arrhythmia sources, possibly leading to detailed catheter mapping and unacceptably long procedures. Procedure times in general may be excessively long, particularly when conducted under prolonged fluoroscopic imaging, leading to excessive x-ray exposure to the patient, physician, and nursing staff. The current invasive options for AFib mapping also have significant disadvantages, particularly when they involve extended and/or traumatic catheter manipulation in the left atrium.  
         [0013]     In light of the above, it would be desirable to provide improved devices, systems, and methods for localizing and/or treating AFib and other arrhythmias within a heart of a patient. The present invention provides such improvements, mitigating and/or overcoming at least some of the disadvantages of known approaches for diagnosing and treating arrhythmias.  
         [0014]     II. Related Art  
         [0015]     The following patents may be relevant to the subject matter of the present invention, and their full disclosures incorporated herein by reference: U.S. Pat. No. 5,311,873; and U.S. Pat. No. 5,634,469. Peeters, H. A. P., SippensGroenewegen, A. and others described “ Clinical Application of an Integrated  3- Phase Mapping Technique for Localization of the Site of Origin of Idiopathic Ventricular Tachycardia ”in  Circulation,  99:1300-1311 (1999). SippensGroenewegen, A. et al. also described “ Body Surface Mapping of A trial Arrhythmias: Atlas of Paced P wave Integral Maps to Localize the Focal Origin of Right Atrial Tachycardia ”, in  J. Electrocardiol.,  31(Supp.): 85-91 (1998). Related work was described by SippensGroenewegen, A. et al. in, “ Value of Body Surface Mapping in Localizing the Site of Origin of Ventricular Tachycardia in Patients with Previous Myocardial Infarction”, J. Am. Coll. Cardiol.  24:1708-1724 (1994). Each of these references is incorporated herein by reference.  
       BRIEF SUMMARY OF THE INVENTION  
       [0016]     The present invention provides improved devices, systems, and methods for localizing and/or treating arrhythmias of a heart. The techniques of the present invention are particularly useful for localizing atrial fibrillation, and allow locating arrhythmogenic regions of a chamber of the heart using heart cycle signals measured from a body surface of the patient. Non-invasive localization of the ectopic origin allows focal treatment to be quickly targeted to effectively inhibit these complex arrhythmias without having to rely on widespread and time consuming sequential searches, and/or on massively invasive simultaneous intracardiac sensor techniques. The invention recognizes that effective localization of these complex arrhythmias can be significantly enhanced by techniques and structures which separate superimposed heart cycle signals originating from differing chambers and/or regions of the heart tissue. In the exemplary embodiment, P wave signal portions are separated from superimposed QRST wave complex signals so as to isolate signals originating in an atrium from concurrent activity in the ventricle. The P wave signals may be measured by a thoracic array of electrical sensors distributed along the patient&#39;s skin across the torso. The invention allows invasive pace-mapping to be limited to a predetermined arrhythmogenic region within a particular chamber of the heart, often followed by ablation of the ectopic origin to inhibit the arrhythmia.  
         [0017]     In a first aspect, the invention provides a method for treating fibrillation in a heart of a patient. The patient has an exposed body surface, and the method comprises measuring the fibrillation from the body surface. An arrhythmogenic region of the heart is located in response to the measured fibrillation. Treatment is directed at or near the arrhythmogenic region so that fibrillation is inhibited.  
         [0018]     Preferably, heart cycle signals are sensed while no intracardiac probe is present in the heart, with the arrhythmogenic region being determined (at least in part) using the non-invasively sensed heart cycle signals. In many cases, an array having more than about 20 sensing locations will be coupled to thoracic skin of the patient. The sensed heart signals will often include signals originating in an arrhythmogenic chamber of the heart, superimposed with signals from other chambers. Typically, an atrial signal will be superimposed with a ventricular signal. These superimposed signals are separable by signal separators. The methods of the present invention often include separating the atrial and ventricular signals with such a signal separator.  
         [0019]     In the exemplary embodiment of the present method, at least one reference cycle is selected from among a plurality of heart cycles measured by the sensor array. The arrhythmogenic region can be determined from the measured signals during this reference cycle, typically by selecting a time portion of the reference cycle, by integrating the separated signals from each sensor location within the selected time portion, and by arranging the resulting integral values within a data matrix according to the locations at which the signals are sensed along the body surface. Such a data matrix may be graphically plotted with the plots often including lines of constant integral values. These plots can be used to identify the arrhythmogenic region, most commonly by comparing the data matrix and/or plots to a database having a plurality of known arrhythmia cycles. Such databases will often have associated known arrhythmogenic regions for each known arrhythmia cycle.  
         [0020]     Advantageously, the locating of the arrhythmogenic region can be performed using measurements of spontaneous fibrillation. The located arrhythmogenic region may have a surface area of less than about five square centimeters, often having an outer radius of less than about 2.5 cm, and ideally having an outer radius of about 1.0 cm or less. An ectopic site or exit site within such a limited arrhythmogenic region may be more precisely identified by subsequently introducing a mapping/pacing probe into the arrhythmogenic atrium. Alternatively, where locating can be performed to sufficient accuracy solely using the sensor array, ablation may directly proceed based on the non-invasively identified arrhythmogenic region. Such ablation will often be performed using radiofrequency energy, cryogenic cooling or electrosurgical energy, buy may alternatively be effected by focal or circular delivery of cryogenic cooling, ablative compounds, ultrasound energy, microwave energy, laser energy, or the like. Circular or perimeter lesions may be particularly beneficial for isolating arrhythmogenic regions in or near pulmonary veins and the like.  
         [0021]     In another aspect, the invention provides a method for treating arrhythmia in a heart of a patient. The patient has an accessible body surface, and the heart has a left atrium, a left ventricle, a right atrium, and a right ventricle. Heart signals at the body surface include atrial signals superimposed with ventricular signals. The atrial and ventricular signals are separable by a separator. A database has information regarding a plurality of known arrhythmia cycles, each known arrhythmia cycle having an associated known arrhythmogenic region. The method comprises sensing signals during an arrhythmia initiation cycle from the body surface. Alternatively, one can also sense premature atrial beats with focal AFib, and/or persistent AFib. The atrial signals are separated from the sensed signals with the signal separator. An arrhythmogenic region (or in some embodiments, an insertion point of a concealed accessory pathway) of an atrium of the heart is located by comparing the separated signals to the database. A treatment is directed at or near an ectopic site or exit site within the arrhythmogenic region so that the arrhythmia is inhibited.  
         [0022]     In another aspect, the present invention provides a system for treating arrhythmia in a heart of a patient. The patient has an accessible body surface and the heart has an atrium and a ventricle. A sensor array can be coupled to the body surface for sensing heart signals, the heart signals including atrial signals superimposed with ventricular signals. The atrial and ventricular signals are separable by a signal separator. A database includes information regarding a plurality of known cycles, each known arrhythmia cycle having an associated arrhythmogenic region. The system comprises a processor coupled to the sensor array and the database. The processor derives an arrhythmogenic region of the atrium from the heart cycle signals by separating the atrial and ventricular signals with the signal separator, and by comparing the separated signals to the database. A probe directs focal treatment at or near the arrhythmogenic region so that the arrhythmia is inhibited.  
         [0023]     In yet another aspect, the invention provides a kit for use with a probe and a sensor array to treat arrhythmia in an atrium of a heart of a patient. The sensor array is coupleable to the body surface for sensing heart cycle signals. The probe has a treatment delivery surface. The kit comprises a processor coupleable with the sensor array. The processor generates an output in response to sensed heart cycle signals from the array. The kit also includes instructions for locating the arrhythmogenic region within the atrium by sensing the heart cycle signals from the body surface when no intracardiac probe is present in the heart, and by comparing the output of the processor with the database.  
         [0024]     In yet another aspect, the invention provides a system for localizing an arrhythmia. The system comprises an input for body surface signals and a processor coupled to the input. The processor derives an arrhythmogenic region of a heart in response to the body surface signals. An output graphically indicates the arrhythmogenic region. 
     
    
     BRIEF DESCRIPTION OF THE DRAWINGS  
       [0025]      FIG. 1A  schematically illustrates a sensor system having an array of sensing locations distributed across a patient&#39;s torso.  
         [0026]      FIG. 1B  graphically illustrates the method for calculating an integral value across a selected time portion of a heart signal cycle from a single sensor location.  
         [0027]      FIG. 1C  illustrates a plot of a data matrix generated by mapping the integral values with positions corresponding to the locations of the sensors across the patient&#39;s torso.  
         [0028]      FIG. 2  schematically illustrates a method and computer program for separating a portion of an electrocardiogram signal relevant to a first portion of a heart from signals relevant to alternative portions of the heart.  
         [0029]     FIGS.  3 A-D graphically illustrate signal portion separation method steps employed by the program and method of  FIG. 2 .  
         [0030]      FIG. 4A  graphically illustrates a database of known atrial paced heart cycles as 17 mean P wave integral maps.  
         [0031]      FIGS. 4B and 4C  illustrate 17 known right atrial ectopic origins associated with the 17 mean P wave integral maps of  FIG. 4A .  
         [0032]     FIGS.  5 A-C illustrate a database of mean P wave integral maps and associated locations of ectopic origins of the left atrium.  
         [0033]     FIGS.  6 A-F illustrate correlations between integral maps of paced heart signal cycles obtained in different patients at a common region in the left atrium.  
         [0034]      FIG. 7  illustrates a database of QRS integral maps and associated ectopic origins within the right ventricle.  
         [0035]      FIG. 8  illustrates a database of QRS integral maps and associated ectopic origins within the left ventricle.  
         [0036]      FIG. 9  illustrates a diagnosis and treatment methodology using an arrhythmia location database and signal separator to help locate and treat the origins of focal AFib and other arrhythmias.  
         [0037]      FIGS. 9A and 9B  illustrate arrhythmogenic regions of the left and right ventricles, respectively.  
         [0038]      FIG. 10  schematically illustrates a system and/or kit for diagnosing and/or treating focal AFib and other arrhythmias, according to the principles of the present invention.  
         [0039]      FIGS. 11A-11H  illustrate exposed and skin-engaging surfaces of four panels supporting heart cycle sensors in an exemplary sensor array structure.  
         [0040]      FIGS. 12A and 12B  are computer screen prints of an AFib localization program showing the localization output. 
     
    
     DETAILED DESCRIPTION OF THE INVENTION  
       [0041]     While the following description is largely directed to localization and/or treatment of focal atrial fibrillation, the methods, devices, and systems of the present invention may be used for a variety of arrhythmias, including both focal and re-entrant arrhythmias (such as those resulting from scars). When used for treatment of re-entrant arrhythmia, treatment may be directed at or near an exit site or insertion point of a pathway. The invention is well suited for use with pulmonary vein isolation therapies now being developed (in which linear, circumferential, and/or perimeter lesions may isolate one or more pulmonary veins to inhibit propagation from triggers or exit sites in or near the veins) by allowing selection of target veins and/or indicating whether vein isolation should be utilized. The invention is particularly useful for localizing paroxysmal AFib, but may also find application for persistent and chronic AFib, localizing an insertion point of a concealed accessory pathway, and the like.  
         [0042]     Referring now to  FIG. 1A , the techniques of the present invention will generally make use of an array  10  of sensors  12  distributed across anterior and posterior skin surfaces of torso T on patient P. Array  10  provides multilead electrocardiogram (ECG) data at a plurality of sensing locations distributed across torso T, typically at over 20 sensing locations, more preferably at over 40 sensing locations, and ideally at 62 or more sensing locations.  
         [0043]     Sensors  12  generally comprise unipolar or bipolar electrodes coupled to the patient&#39;s skin, or to an alternative accessible body surface (for example via a transesophageal approach) suitable for measuring electrical body surface potential. Suitable electrode structures may include those described in U.S. Pat. Nos. 5,311,873 and 5,634,496, previously incorporated herein by reference. Exemplary arrays for use in locations having large amounts of electromagnetic noise (such as an electrophysiology lab or other location in which electrosurgery or electrical stimulation of tissues for intracardiac pacing is performed) was described by Metting van Rijn, A. C. et al. in  IEEE Trans. Biomed. Eng ., BME-40:302-308; (1993). Alternative sensor array structures and associated data acquisition and manipulations components were described by SippensGroenewegen, A. et al. in an article entitled, “ Body Surface Mapping During Pacing at Multiple Sites in the Human Atrium: P wave Morphology of Ectopic Right Atrial Activation”, Circulation,  97:369-380 (1998); and by Linnenbank, A. C. in a 1996 thesis for the University of Amsterdam entitled, “ On - Site Recording, Analysis, and Presentation of Multi - channel ECG Data”.    
         [0044]     Referring now to  FIG. 1B , ECG data is preferably acquired simultaneously from each sensor  12  of array  10  at a sampling rate of over about 500 Hz, ideally at a sampling rate of about 1,000 Hz or more. In some embodiments, sequential sampling of sensor  12  from array  10  may alternatively be used, and higher or lower sampling rates are also feasible. When a lower sampling rate is used, the data may be upsampled using multi-rate filter banks.  
         [0045]     Preferably, signals which are absent, for example, due to electrode obscurement by defibrillator patches or lead dislodgment, may be deleted. Poor quality signals may also be visually and/or automatically identified and rejected. Such rejected signals may be replaced using interpolation of adjacent lead recording data. Interpolation techniques may also be utilized to correct for offset variation among electrodes, and for linear baseline drifting.  
         [0046]     Graph  14  includes an ECG signal tracing  16  representing the variation in voltage over time, as sensed by sensors  12 , optionally at about 1 to 2 ms intervals. Signal tracing  16  may be used to evaluate heart cycle signals from the heart of patient P. In general, one or more reference heart cycles will be selected for manipulation and comparison. The reference heart cycle will typically be the heart cycle coinciding with initiation of the arrhythmia for localizing focal AFib, or a premature heart cycle that has the same morphology as the cycle that led to initiation of focal AFib. ECG Tracing  16  can be used to determine a beginning  18  and end  20  of a time portion  22  of the reference heart signal cycle which is of particular interest for evaluating one or more regions of the heart. In the example illustrated in  FIG. 1B , a P wave onset may be determined by the time at which the voltage progresses beyond 30 μV while termination of the P wave may be defined at the atrial J-point, as is generally understood in the field of electrocardiography. Alternative criteria for P wave onset and offset might also be utilized, and automated detection of time portion  22  is also feasible. Alternative time portions may also be selected.  
         [0047]     Referring to  FIGS. 1B and 1C , measurements made at each sensor  12  are preferably mapped onto a data matrix  24  according to the locations of the associated sensor. In the exemplary embodiment, a P wave integral numerical value  24  may be calculated based on heart cycle signals  16  within selected time portion  22  for a particular sensor location N 10 . This calculated P wave integral value reflects the time/amplitude area of ECG signal at that sensor location within the selected time portion. Similar integral values are calculated for each sensor location, and the sensor values are mapped within data matrix  24  continuously from a portion of the data matrix associated with a front F of torso T, across a side of the patient P, and to a back B portion of torso T. As shown in  FIG. 1C , the data matrix will often be presented graphically by calculating lines of constant integral values  28  based on the individual discrete integral values and their associated positions within the data matrix. In some embodiments, this information can be summarized by presenting a single line  30  of zero integral value between a region of positive integral values  32  and a region of negative integral values  34 . In much of the description which follows, the region of positive integral values  32  is presented as a shaded region within a graphically depiction of data matrix  24 . Exemplary alternative data matrices may be presented with shades of a first color (red, for example) for positive values, a second color (blue, for example) for negative values, and optionally a third color (such as green) for zero.  
         [0048]     For localizing of certain arrhythmias, possibly including certain ventricular tachycardias and some types of atrial tachycardia, directly using measurements from sensors  12  to calculate integral values  26  for the selected time portion  22  maybe sufficient to identify an arrhythmogenic region (which may be relatively large) of a particular ventricle, and in some cases, a particular atrium. Localizing directly from the sensed heart cycle signals is significantly facilitated when the signals within the time portion of interest are predominantly indicative of activity within a candidate ectopic region of the heart. For example, when localizing ventricular tachycardia (VT), selecting a time portion dominated by the QRS complex in the signal can effectively localize arrhythmogenic foci or exit site, as more fully described in the  J. Am. Coll. Cardiol.,  24:1708-1724 (1994), the full disclosure of which is incorporated herein by reference. This localizing of tachycardia foci within the ventricle may be facilitated by the domination of the QRS complex in the signal of the overall body surface potential.  
         [0049]     Unfortunately, when localizing fibrillation foci within an atrium, the P wave (which can be indicative of activity within the atrium) will often be superimposed, either partially or completely, by the T-U wave. Physiologically speaking, the atrial activity of interest may coincide with ventricular recovery of the preceding cardiac cycle.  
         [0050]     To accurately localize focal triggers during the initiation of paroxysmal or persistent AFib, the present invention makes use of systems and methods for effectively separating a signal portion of interest from a superimposed signal portion, with the two signal portions often being separated from a single signal sensed from at least one single sensor location. These signal separation techniques are particularly advantageous when used to isolate the P wave from a simultaneously occurring T-U wave. It may be possible in some circumstances to artificially separate these waves by active overdrive pacing using an intracardiac catheter with a pacing period selected to avoid superimposition of these two signal portions during artificially initiated arrhythmia. For the reason described above, non-invasive electrocardiographic localization of atrial arrhythmias, particularly atrial premature beats before or during invasive mapping procedures is highly advantageous. As will be understood with reference to  FIGS. 2-3D , a QRST subtraction program helps to isolate and preserve the P wave morphology so as to enable trigger localization of focal AFib and other arrhythmias. The application of similar subtraction methodologies may also enhance the ability of body surface mapping systems to localize triggers, exit sites, pathway insertion points, flutter waves, and/or fibrillation waves of other atrial arrhythmias such as atrial flutter, chronic AFib, and the like. QRST subtraction program  40  may also enable application of inverse problem techniques to analysis of atrial arrhythmias, for example, when atrial depolarization is obscured by the preceding ventricular repolarization. It should be understood that alternative signal separation methods and systems might also be used, including those described in the following references, which are incorporated herein by reference: Slocum, J. et al., “ Computer Detection of Atrioventricular Dissociation from Surface Electrocardiograms During Wide QRS Complex Tachycardias,” Circulation,  72:1028-1036 (1985); Slocum, J. et al., “ Diagnosis of A trial Fibrillation from Surface Electrocardiograms Based on Computer - Detected Atrial Activity,” J. Electrocardial.,  25:1-8 (1992); Holm, M. et al., “ Non - Invasive Assessment of the Atrial Cycle Length During Atrial Fibrillation in Man: Introducing, Validating and Illustrating a New ECG Method,” Cardiovasc. Res.,  38:69-81 (1998); Bollman, A. et al., “ Frequency Analysis of Human Atrial Fibrillation Using the Surface Electrocardiogram and its Response to Ibutilide,” Am. J. Cardiol.,  81:1439-1445 (1998); and Ingemansson, M. P. et al., “ Autonomic Modulation of the Atrial Cycle Length by the Head Up Tilt Test: Non - Invasive Evaluation in Patients with Chronic Atrial Fibrillation,” Heart,  80:71-76 (1998).  
         [0051]     Referring now to  FIG. 2 , automated QRST subtraction program  40  uses an adaptive QRST template constructed from averaged QRST complexes combined with ECG body surface measurements to enable isolation of the otherwise obscured ectopic atrial activity. Generally, this approach allows the surface ECG measurements to retain their intricate spatial and temporal detail within the P wave morphology. Subtraction program  40  is capable of unmasking and preserving subtle heart signal details of relatively low voltage P wave signal portions despite the obscuring superimposed relatively high voltage QRST complex. The QRST subtraction method of  FIG. 2  is described in more detail in U.S. Patent Application No. 60/189,513 filed Mar. 15, 2000, previously incorporated herein by reference.  
         [0052]     As described above, the method of program  40  generally includes recording of unipolar ECG data from the array of torso sites in step  42 . The measured signal will include both the P wave (which is of interest for AFib) and a superimposed QRST signal portion  44 .  
         [0053]     In the exemplary embodiment, about 100 cardiac cycles of 62-channel ECG data are measured during sinus rhythm or atrial overdrive pacing. Optionally, ECG signals can be acquired during both sinus rhythm and atrial pacing. Fewer cycles may be used if the spatial and temporal variations of the QRST complex are relatively low. Typically, more than ten (10) cycles will be used, often more than 50 cycles for construction of the QRST template. As illustrated in  FIG. 3A , off-line digital filtering of the data designated for template construction may be performed using a 0.5 Hz high-pass filter, such as an IIR Chebychev Type-1 filter. This can help to correct for respiration-related baseline drifting. Additionally, a 100 Hz low-pass filter, such as a IIR Chebychev Type-1 filter, may be used to remove high-frequency signal artifacts. Additionally, a 50-60 Hz notch filter may be used to remove line-frequency interference. Similar filtering may be employed on the superimposed signals to be separated.  
         [0054]     Each filtered QRS complex  48  may be identified using a complex-resonator/comb filter, together with a dual-edge threshold detection technique similar to that described by Ruha, A. et al., in an article entitled “ A Real - Time Microprocessor QRS Detector System with a  1- ms Timing Accuracyfor the Measurement of Ambulatory HR V”, IEEE Trans. Biomed. Eng.,  44:159-167 (1997), the disclosure of which is incorporated herein by reference. Alternative QRS detection methods might also be used.  
         [0055]     R wave fiducial points are marked and the average R—R interval is computed. The dominant QRS morphology is identified, optionally using visual identification from pooled data, automated statistical methods, or the like. This dominant QRS morphology is used to select complexes for template creation.  
         [0056]     QRST template construction  50  may be understood with reference to  FIG. 3B . The selection of complexes for template creation may be based on two criteria: QRS pattern, and R—R interval length.  
         [0057]     Regarding QRS pattern criteria, complexes of each cycle are compared with the template using parametric cross-correlations. In the exemplary embodiment, Pearson&#39;s coefficients are computed for a fixed time window sliding over a 20 ms time period. Once again, a variety of alternative cross-correlations may be used. Each newly selected complex (for example, those having r≧0.98) is aligned and averaged with previously selected cycles, and the QRS template is updated. Each complex which does not have adequate correlation with the template is excluded, with these complexes often being ectopic or aberrantly conducted ventricular beats.  
         [0058]     Minimum R-R duration threshold is computed and only complexes having an R-R duration above the computed threshold are used to create the QRST template. For each complex, the Q wave fiducial point is identified and the average Q-R interval is computed. The average Q-T interval (the length of the QRST complex) is computed, optionally using a modified Bazett&#39;s formula. Additional correction for baseline drifting may be performed using linear interpolation after averaging a window prior to the Q wave. An adaptive template may be constructed from the selected complexes by averaging their QRST intervals, optionally with the complexes aligned by a window surrounding the R wave fiducial point.  
         [0059]     In step  52 , the above-described fiducial window in the QRS complexes, together with an additional fiducial time window around the peak of the T wave, are marked automatically. These fiducial windows are marked in both the QRST template, and in the superimposed signal  48  containing a QRST complex together with a superimposed P wave.  
         [0060]     In step  54 , the template  56  and the superimposed signals are aligned. Typically, T wave window onset for the template is estimated as 0.64*QT. The T wave window is adjusted for each complex based on the QT interval length. Optionally, the operator can manually adjust any fiducial marker, such as by manipulating a mouse, joystick, by putting a numerical value, or the like, the template&#39;s QRS fiducial points and T wave windows are aligned with the respective windows of the superimposed QRST complex. Alignment of superimposed signal  48  and template  56  may be performed manually, or automatically by sliding windows over each other in 10-ms increments and calculating cross-correlation coefficients.  
         [0061]     After alignment of the fiducial windows, the template is resampled as well as amplitude modulated  58 , and then QRST subtraction is carried out as can be understood with reference to  FIG. 3D  and step  60 . Template resampling and modulation are performed to compensate for discrepancies in duration (as a result of rate-related differences in the QRST interval) and amplitude (as a consequence of variations in peak R-T wave voltage, which may be predominantly caused by respiratory variation). This process is particularly geared toward obtaining optimal subtraction performance of the T-U wave complex. Optionally, cubic spline interpolation methods may be used for template resampling, although other multi-rate processing methods can also be used.  
         [0062]     The subtraction procedure from channel-to-channel may be performed according to electrode position. This facilitates maintaining sliding intervals as small as 10 ms while having the Q-T interval dispersion accessed in multiple leads. Additional low-pass filtering can be carried out to smooth possible QRS leakage after QRST subtraction. The remaining ECG signal after QRST subtraction features a P wave which is effectively isolated from the previously superimposed T-U wave.  
         [0063]     After the QRST segment is effectively subtracted from each measured signal containing a P wave of interest, the morphology of the isolated P wave can be analyzed as described above. Specifically, an integral map may be computed of the separated, previously superimposed P wave. This integral map can be compared with a database of P wave maps created by pacing, so that atrial tachycardia, atrial accessory pathway insertion sites, and focal triggers of paroxysmal or persistent AFib may be localized using data from the surface ECG array  10 , as shown in  FIG. 1A . Advantageously, the method described herein above may obtain high performance in the T-U wave range by correcting specifically for differences in both the QRST duration, and in the voltage of the T wave. Additionally, the above-described QRST subtraction methodology makes use of a separate data set than the superimposed wave to be separated, with the separate data set optionally being obtained during sinus rhythm or atrial overdrive pacing to help insure that atrial and ventricular activity are clearly separated. As was also mentioned above, the superimposed wave or reference cycle of interest will often comprise a single ectopic atrial heartbeat which can be readily separated using the above-described method.  
         [0064]     Once the signal separation has been successfully completed for the reference heartbeat at each sensor location  12  of array  10 , a P wave integral map may be plotted from the data matrix  26  as described above with reference to FIGS.  1 A-C. Similar QRS integral data matrices and plots may be generated for localizing ventricular arrhythmias.  
         [0065]     Referring now to FIGS.  4 A-C, a graphical plot of a particular patient&#39;s P wave integral may be used to localize an arrhythmogenic region in an atrium by comparing the P wave integral plot for the patient to a database of P wave integral plots having associated known ectopic foci within the right atrium. Each of the 17 plots of database  70  has an associated ectopic region (identified by the encircled numbers illustrated in  FIGS. 4B and 4C ), the database having been gathered using pacing.  
         [0066]     The anterior-posterior view AP shown in  FIG. 4B  and the posterior-anterior PA view of  FIG. 4C  illustrate the right atrial cavity. Anatomical landmarks included in these diagrams include the superior vena cava SVC, and inferior vena cava IVC; the right atrial appendage RAA; the smooth right atrium SRA; the trabeculated right atrium TRA; the crista terminalis CT; the fossa ovalis FO; the left atrium LA; the Eustachian valve EV; the coronary sinus os CSO; the tricuspid valve TV; the right pulmonary artery RPA; and the left pulmonary artery LPA.  
         [0067]     Methods for assembling a right atrial database are described in detail in the  J. Electrocardiol.,  31 (Supp.): 85-91 (1998), previously incorporated herein by reference. A similar left atrial database  80   a, b , and  c  (generally referred to as left atrial database  80 ) is illustrated in FIGS.  5 A-C. Once again, the encircled numbers relate mean P wave integral plots to specific endocardial regions of pacing shown in the anatomical diagrams.  
         [0068]     The mean P wave integral maps of left and right atrial databases  70 ,  80  feature extreme positions and zero line contours without positive and negative integral contour lines. Alternative plot formats, such as three-dimensional or chest anatomy-based formats, map displays using various color schemes, and the like, may also be used.  
         [0069]     The anatomical diagrams illustrated in FIGS.  5 A-C present a posterior-anterior PA view and anterior-posterior AP view, and a left posterior-oblique LPO view of the left atrium LA. Once again, major anatomical landmarks are highlighted including the left and right pulmonary arteries LPA, RPA; the superior and inferior vena cava SVC, IVC; the left atrial appendage LAA; the right atrium RA; the coronary sinus CS; the left ventricle LV; the left upper pulmonary vein LUPV; the left lower pulmonary vein LLPV; the right upper pulmonary vein RUPV; the right lower pulmonary vein RLPV; Bachmann&#39;s bundle BB; the mitral annulus MA; the anterior mitral valve leaflet AMVL; and the fossa ovalis FO.  
         [0070]     Referring now to FIGS.  6 A-F, six individual P wave integral maps included within group  6  were each obtained during pacing at the left upper and left lower pulmonary veins. These six similarly located pacing sites were grouped together within group  6  of left atrial database  80 , and these plots were averaged to produce the sixth numbered mean plot of the left atrial database. The spatial compatibility of these patterns can be clearly seen, particularly with reference to the location and orientations of both the highest positive and negative integral values, as well as with reference to the zero line contour separating the shaded from unshaded regions. While each of these six patient-specific maps were generated using intracardiac pacing, naturally occurring ectopic origins may be identified by comparing reference heart cycle signals measured during premature atrial beats, the onset of AFib, and/or atrial tachycardia (and optionally separated from superimposed signals described above) to these mean paced plots.  
         [0071]     Referring now to  FIGS. 7 and 8 , a right ventricular database  82  and a left ventricular database  84  each include mean QRS integral maps for paced ectopic origins in the right and left ventricles, respectively. These ventricular databases are more fully described in an article by Peeters, H. A. P. et al. entitled “ Clinical Application of an Integrated  3- Phase Mapping Technique for Localization of the Site of Origin of Idiopathic Ventricular Tachycardia”, Circulation  99:1300-1311 (1999) the disclosure of which is incorporated herein by reference.  
         [0072]     Referring now to  FIG. 9 , a localization/treatment method  100  may be performed by first establishing a database of arrhythmogenic regions and associated heart cycle signal characteristics, as described with reference to  FIGS. 4A-8  above. Body surface measurements are taken in step  102 , typically using sensor array  10  described with reference to  FIG. 1A . As the P wave will often overlap other heart signal portions, localization of AFib and other arrhythmias will benefit from a signal separation program as described above regarding  FIGS. 2-3D .  
         [0073]     Once signal separation has been effected, P wave integrals (or other selected indicia) are determined for each sensor  12 , and data matrices assembled as described with reference to  FIGS. 1B and 1C . Once again, it should be understood that these P wave integrals are calculated for a reference heart cycle signal of interest, usually for a premature atrial beat or for the heart cycle signal at the onset of an arrhythmia, for focal arrhythmias. Optionally, the data matrix may be compared with the plots of the established database in step  104 . This database comparison method may simply involve visually selecting the mean-paced integral map which appears most similar to the plot for a specific patient. Alternatively, statistical correlation coefficients may be generated between the data matrix for the patient and each of the mean-paced groups of the database. In some embodiments, different patient plots may be prepared for comparison with different databases, for example, a P wave integral may be calculated for comparison with atrial database groups, while a QRS integral plot may be prepared for comparison with ventricular database groups.  
         [0074]     Once a mean-paced integral plot from the database has been selected as the closest correlation to the reference cycled plot for a particular patient, an arrhythmogenic region  106  associated with the corresponding mean-paced plot has effectively been identified. Arrhythmogenic regions  106  associated with mean-paced plots  1 - 25  of left ventricular database  84  and plots  1 - 13  of right ventricular database  82  are illustrated in  FIGS. 9A and 9B , respectively. In many embodiments, these arrhythmogenic regions will be discrete locations based on the information within the associated database. Preferably, arrhythmogenic regions  106  will have surface areas of less than about 5 cm 2 . Optionally, the arrhythmogenic regions may have an outer radius which is less than about 2.5 cm, ideally about 1.0 cm or less. In some embodiments, the arrhythmogenic regions identified by sensor array  10  on the patient&#39;s torso may be small enough that no further localization is needed, and ablation of the ectopic site or exit site within the arrhythmogenic region may proceed without excessive collateral damage.  
         [0075]     Referring now to  FIGS. 9 and 10 , once an arrhythmogenic region  106  has been identified, it will often be advantageous to further localize an ectopic origin or exit site  108  within arrhythmogenic region  106  using a pace mapping catheter  110 . Advantageously, this mapping of ectopic origin or exit site  108  may proceed rapidly within the limited confines of the arrhythmogenic region  106 , thereby reducing fluoroscopy time and radiation exposure to patient and attending personnel, decreasing the trauma associated with accessing alternative portions of heart H, and the like. In general, pace mapping is effected by electrical stimulation of candidate ectopic origins within arrhythmogenic region  106  using a distal electrode pair of catheter  110 . This can induce ectopic heart beats which can also be measured by array  10  (or a similar array adapted for use in a high electromagnetic noise environment). Positioning of the catheter tip may be monitored using biplane x-ray imaging.  
         [0076]     The surface ECG corresponding to the paced beats is recorded, and the desired integrals and associated data matrix is generated, as described above. By comparing the data matrix plot of the induced arrhythmia to the database and/or original arrhythmia recordings, the focal origin or exit site of the arrhythmia relative to the catheter position can be estimated, optionally using the method described in U.S. Pat. No. 5,311,873, previously incorporated herein by reference. The catheter may be moved to the indicated alternative site, and the pacing and measurement steps repeated iteratively until the ectopic site is found where the paced data matrix plot best correlates with the data matrix plot of the arrhythmia. Once again, this iterative process is greatly expedited by concentrating the ectopic origin or exit site search within arrhythmogenic region  106  identified using sensor array  10 .  
         [0077]     Once the ectopic origin or exit site has been sufficiently localized, ablation of the ectopic origin  109  is effected, often using an ablation electrode of pacing catheter  110 . A variety of alternative tissue treatment modalities might be applied to the ectopic origin or exit site, including radiofrequency ablation, cryogenic cooling, ultrasound, laser, microwave, bioactive agents, and the like. Similarly, a variety of intracardiac localization techniques might be used in place of intracardiac pace mapping  108  under fluoroscopy. Suitable three-dimensional electro-anatomical point-by-point mapping systems may be commercially available for localization of an ectopic origin or exit site within an arrhythmogenic region from BIOSENSE-WEBSTER, INC. under the trademark CARTO®, and a related Real-Time Position Management™ system may be available from CARDIAC PATHWAYS CORPORATION. Alternative multi-electrode catheters may be commercially available from CARDIMA, INC., BIOSENSE-WEBSTER, INC., CARDIAC PATHWAYS CORPORATION, BARD, INC. and/or EP TECHNOLOGIES, INC. A still further alternative for localizing of the ectopic origin within an arrhythmogenic region maybe provided using a three-dimensional non-contact multi-electrode mapping system under development by ENDOCARDIAL SOLUTIONS, INC. Mapping/ablation catheters may also be available from both ENDOCARDIAL SOLUTIONS, INC. and CARDIMA, INC. Exemplary cryogenic systems may be available from CRYOCATH, INC. and from CRYOGEN, INC. A suitable cooled radiofrequency ablation catheter is sold commercially as the CHILLI®-Cooled Ablation System from CARDIAC PATHWAYS CORPORATION. Pulmonary vein isolation systems for use with the invention are now being developed by ATRIONIX (ultrasound) and CARDIOFOCUS (laser ablation).  
         [0078]     Referring to  FIG. 10 , a kit  120  for localization of an arrhythmogenic region of heart H may include a tangible media having a machine-readable code embodying any of the methods described herein above. Media  122  will often be used in a general-purpose computer PC coupleable to array  10 , with the PC typically having a processor for effecting the method steps embodied in media  122 , as well as input devices such as a mouse, keyboard, an Internet, Ethernet, and/or Intranet, as well as output devices such as monitor, a printer, an I/O port, and/or the like. The PC will often manipulate the data in response to heart cycle signals sensed by array  10 , and also in response to a heart cycle signal database  124 , as described herein above. Instructions for use  126  will often be included within kit  120 , with the instructions and at least one other component of the kit often being packaged together, instructions for use  126  optionally being embodied as printed information (optionally appearing at least in part on the packing material), a VCR tape, media embodying a machine readable code, or the like.  
         [0079]     In many embodiments a system  130  will include one or more components of kit  120 , and may also include array  10  and/or a probe such as catheter  110  for localization and/or treatment of the arrhythmias.  
         [0080]     In some embodiments, continuous localization of local triggers or exit sites may be provided using correlations such as those described in an article by Potse, M. et al., entitled “ Continuous Localization of Cardiac Activation Sites Using a Database of Multichannel ECG Recordings,” IEEE Trans. Biomed. Eng.,  47:682-689 (2000), the full disclosure of which is incorporated herein by reference.  
         [0081]     FIGS.  11 A-H illustrate skin-engaging and outer surfaces of four flexible sensor panels. These panel structures and the use of their array or sensors are more fully described in U.S. patent application Ser. No. 09/611,179, filed on Jul. 6, 2000, the disclosure of which is incorporated herein by reference.  
         [0082]     Referring now to  FIGS. 12A and 12B , screen prints  150   a ,  150   b  are at least a portion of the output of a localization computer program as described herein. Integral maps  152  are multi-color graphical representations of the morphology derived from body surface potential maps  154  (also displayed in a multi-color format). Time portion  156  of reference cycle  158  is shown graphically, and may be separated from a superimposed signal by inputting a “subtract” command  160 , ideally using a graphical user interface. Identified arrhythmogenic sites  12  are output with their associated probabilities, and are graphically illustrated  164  relative to the adjacent anatomy. In the exemplary output, candidate sites of a selected database are also shown, with one or more of the identified sites being highlighted using color, blinking, an enhanced font or icon, or the like. Optionally, cycle signals from different recorded segments  166  may be selected, and the cycles of interest may be time scaled or zoomed  168  to show the desired cycle intervals. Signals from selected leaps or sensor panels may be displayed, and a variety of additional outputs may be provided. In the output illustrated in screen prints  150   a ,  150   b , two different reference cycles from different arrhythmia events result in similar identified arrhythmogenic sites. Hence, the present invention can also provide confirmation of an arrhythmogenic site identification from one or more confirmation reference cycle.  
         [0083]     While the exemplary embodiments have been described in some detail, by way of example and for clarity of understanding, a variety of changes, adaptations, and modifications will be obvious to those of skill in the art. Hence, the scope of the present invention is limited solely by the appended claims.