EP2407118A2 - Methods and systems for filtering respiration noise from localization data - Google Patents
Methods and systems for filtering respiration noise from localization data Download PDFInfo
- Publication number
- EP2407118A2 EP2407118A2 EP11173450A EP11173450A EP2407118A2 EP 2407118 A2 EP2407118 A2 EP 2407118A2 EP 11173450 A EP11173450 A EP 11173450A EP 11173450 A EP11173450 A EP 11173450A EP 2407118 A2 EP2407118 A2 EP 2407118A2
- Authority
- EP
- European Patent Office
- Prior art keywords
- signal
- localization
- respiration
- acceleration
- displacement
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 230000004807 localization Effects 0.000 title claims abstract description 195
- 230000029058 respiratory gaseous exchange Effects 0.000 title claims abstract description 68
- 238000000034 method Methods 0.000 title claims abstract description 43
- 238000001914 filtration Methods 0.000 title claims abstract description 11
- 238000006073 displacement reaction Methods 0.000 claims abstract description 68
- 230000001133 acceleration Effects 0.000 claims abstract description 58
- 238000005259 measurement Methods 0.000 claims abstract description 26
- 238000012545 processing Methods 0.000 claims description 9
- 230000001131 transforming effect Effects 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 4
- 210000001015 abdomen Anatomy 0.000 description 11
- 230000000747 cardiac effect Effects 0.000 description 9
- 230000000694 effects Effects 0.000 description 7
- 238000013507 mapping Methods 0.000 description 7
- 230000000241 respiratory effect Effects 0.000 description 7
- 230000008901 benefit Effects 0.000 description 6
- 230000005684 electric field Effects 0.000 description 6
- 230000007831 electrophysiology Effects 0.000 description 6
- 238000002001 electrophysiology Methods 0.000 description 6
- 230000033001 locomotion Effects 0.000 description 6
- 210000005242 cardiac chamber Anatomy 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 210000000038 chest Anatomy 0.000 description 3
- 210000004072 lung Anatomy 0.000 description 3
- 230000004075 alteration Effects 0.000 description 2
- 239000005441 aurora Substances 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000002560 therapeutic procedure Methods 0.000 description 2
- 206010003658 Atrial Fibrillation Diseases 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 210000003484 anatomy Anatomy 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 238000009125 cardiac resynchronization therapy Methods 0.000 description 1
- 239000004020 conductor Substances 0.000 description 1
- 210000003748 coronary sinus Anatomy 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 210000005246 left atrium Anatomy 0.000 description 1
- 210000005240 left ventricle Anatomy 0.000 description 1
- 210000002414 leg Anatomy 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000036391 respiratory frequency Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 210000001562 sternum Anatomy 0.000 description 1
- 230000001225 therapeutic effect Effects 0.000 description 1
- 210000000689 upper leg Anatomy 0.000 description 1
- 210000005166 vasculature Anatomy 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/06—Devices, other than using radiation, for detecting or locating foreign bodies ; determining position of probes within or on the body of the patient
- A61B5/061—Determining position of a probe within the body employing means separate from the probe, e.g. sensing internal probe position employing impedance electrodes on the surface of the body
- A61B5/063—Determining position of a probe within the body employing means separate from the probe, e.g. sensing internal probe position employing impedance electrodes on the surface of the body using impedance measurements
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
- A61B5/7207—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise induced by motion artifacts
- A61B5/721—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise induced by motion artifacts using a separate sensor to detect motion or using motion information derived from signals other than the physiological signal to be measured
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B90/00—Instruments, implements or accessories specially adapted for surgery or diagnosis and not covered by any of the groups A61B1/00 - A61B50/00, e.g. for luxation treatment or for protecting wound edges
- A61B90/36—Image-producing devices or illumination devices not otherwise provided for
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B17/00—Surgical instruments, devices or methods, e.g. tourniquets
- A61B2017/00017—Electrical control of surgical instruments
- A61B2017/00022—Sensing or detecting at the treatment site
- A61B2017/00039—Electric or electromagnetic phenomena other than conductivity, e.g. capacity, inductivity, Hall effect
- A61B2017/00044—Sensing electrocardiography, i.e. ECG
- A61B2017/00048—Spectral analysis
- A61B2017/00053—Mapping
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
- A61B34/20—Surgical navigation systems; Devices for tracking or guiding surgical instruments, e.g. for frameless stereotaxis
- A61B2034/2046—Tracking techniques
- A61B2034/2048—Tracking techniques using an accelerometer or inertia sensor
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
- A61B34/20—Surgical navigation systems; Devices for tracking or guiding surgical instruments, e.g. for frameless stereotaxis
- A61B2034/2046—Tracking techniques
- A61B2034/2051—Electromagnetic tracking systems
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
- A61B34/20—Surgical navigation systems; Devices for tracking or guiding surgical instruments, e.g. for frameless stereotaxis
- A61B2034/2046—Tracking techniques
- A61B2034/2051—Electromagnetic tracking systems
- A61B2034/2053—Tracking an applied voltage gradient
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2562/00—Details of sensors; Constructional details of sensor housings or probes; Accessories for sensors
- A61B2562/02—Details of sensors specially adapted for in-vivo measurements
- A61B2562/0219—Inertial sensors, e.g. accelerometers, gyroscopes, tilt switches
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
- A61B5/113—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb occurring during breathing
Definitions
- the instant invention relates to the collection of data regarding a patient, such as geometric data regarding a portion of the patient's anatomy, using a localization system.
- the instant invention relates to methods and systems for processing a localization signal that describes the position, and optionally the orientation, of a device within the localization field to account for signal artifacts attributable to patient respiration.
- mapping catheter is introduced into the heart chamber of interest and moved around within the heart chamber, either randomly, pseudo-randomly, or according to one or more preset patterns.
- the three-dimensional coordinates are measured using a localization system (sometimes also referred to as a “mapping system,” “navigation system,” or “positional feedback system”).
- the localization system measures the coordinates of the mapping catheter within a localization field, typically by relating a characteristic of the localization field, such as a voltage, experienced by the mapping catheter to a location of the catheter within the field.
- Localization signals produced by localization systems often exhibit noise.
- One prevalent source of noise is patient respiration.
- Patient respiration may introduce noise by moving the generators and/or references for the localization field or by changing the patient's biological characteristics (e.g ., impedance).
- an electrode within the localization field that is, in fact, stationary may appear to be moving.
- Another object of the present invention is to provide a method of collecting localization data that ensures the collected data relates to a common fiducial point within a patient's respiration cycle, such as peak inspiration or peak expiration.
- the localization system includes a plurality of localization system patch electrodes that generate a non-ionizing localization field, and at least one of the localization system patch electrodes includes an accelerometer capable of measuring acceleration in at least one direction.
- the method includes the following steps: acquiring a localization signal of at least one position measurement sensor within the non-ionizing localization field; and acquiring an acceleration signal for the at least one localization system patch electrode; for a selected coordinate axis of the non-ionizing localization field: calculating a displacement signal of the at least one localization system patch electrode on the selected coordinate axis from the acceleration signal; transforming the displacement signal from the time domain into the frequency domain; calculating a fractional power of the transformed displacement signal indicative of patient respiration; comparing the calculated fractional power to a threshold value; and if the calculated fractional power exceeds the threshold value, applying a filter to the localization signal of the at least one position measurement sensor to filter out the respiration noise component on the selected coordinate axis.
- the steps of calculating a displacement signal of the at least one localization system patch electrode on the selected coordinate axis from the acceleration signal, transforming the displacement signal from the time domain into the frequency domain, calculating a fractional power of the transformed displacement signal indicative of patient respiration, comparing the calculated fractional power to a threshold value, and if the calculated fractional power exceeds the threshold value, applying a filter to the localization signal of the at least one position measurement sensor to filter out the respiration noise component on the selected coordinate axis may be repeated for each coordinate axis of the localization field.
- the step of calculating a fractional power of the transformed displacement signal indicative of patient respiration includes calculating a ratio of an area under the transformed displacement signal between about 0 Hz and about 0.3 Hz to an area under the transformed displacement signal as a whole.
- the at least one localization system patch electrode includes a localization field generator that defines the selected coordinate axis.
- Te step of calculating a displacement signal of the at least one localization system patch electrode on the selected coordinate axis from the acceleration signal may include: resolving a component of the acceleration signal on the selected coordinate axis; and integrating the component of the acceleration signal on the selected coordinate axis twice.
- a Fourier transform may be applied to the displacement signal in order to transform the displacement signal from the time domain into the frequency domain.
- the step of applying a filter to the localization signal of the at least one position measurement sensor to filter out the respiration noise component on the selected coordinate axis includes applying a high-pass filter to the localization signal.
- the cutoff frequency for the high-pass filter is selected to be above the patient's respiration frequency.
- the threshold value may be selected as the cutoff frequency.
- the cutoff frequency may be determined using a localization signal acquired from a pseudo-stationary reference electrode positioned within the localization field. In still other embodiments of the invention, the cutoff frequency is determined using the acquired acceleration signal.
- the patient's respiration frequency may be identified by: acquiring a displacement signal for at least one pseudo-stationary electrode on at least one coordinate axis of the localization field; transforming the acquired displacement signal from the time domain to the frequency domain; and identifying the patient respiration frequency based upon a low frequency peak of the transformed signal.
- the present invention includes a method of filtering a respiration noise component from a localization signal of a localization system that uses a plurality of localization system patch electrodes to generate a non-ionizing localization field, wherein at least one of the patch electrodes includes an accelerometer capable of measuring acceleration in at least one direction.
- the method includes the following steps: acquiring an acceleration signal for the at least one localization system patch electrode; and gating collection of data points from a localization signal of at least one position measurement sensor within a localization field based upon the acquired acceleration signal.
- the step of gating collection of data points from a localization signal includes: integrating the acquired acceleration signal to calculate a velocity signal and collecting data points from the localization signal only when the calculated velocity signal changes sign from positive to negative. In other aspects of the invention, the gating step includes: integrating the acquired acceleration signal to calculate a velocity signal and collecting data points from the localization signal only when the calculated velocity signal changes sign from negative to positive. In still other aspects of the invention, the step of gating collection of data points from a localization signal includes: integrating the acquired localization signal twice to calculate a displacement signal and collecting data points from the localization signal only when the calculated displacement signal is at a maximum value. In yet further aspects of the invention, the step of gating collection of data points from a localization signal includes: integrating the acquired localization signal twice to calculate a displacement signal and collecting data points from the localization signal only when the calculated displacement signal is at a minimum value.
- the localization system includes: a plurality of localization field generators capable of generating a non-ionizing localization field; at least one localization field reference device; at least one accelerometer capable of measuring acceleration in at least one direction coupled to at least one of the plurality of localization field generators and the at least one localization field reference device; and a signal processor that acquires at least one acceleration signal from the at least one accelerometer and that uses the acquired at least one acceleration signal to filter a respiration noise component from a localization signal of at least one position measurement sensor within the non-ionizing localization field.
- the signal processor integrates the at least one acceleration signal twice to compute at least one displacement signal and uses the computed at least one displacement signal to filter the respiration noise component from the localization signal.
- the signal processor may also use a high-pass filter having a cutoff frequency determined using the acquired at least one acceleration signal to process the localization signal.
- the signal processor may process the localization signal by gating collection of data points from the localization signal using the acquired at least one acceleration signal.
- An advantage of the present invention is that it permits the filtering of a respiration noise component from a localization signal.
- Another advantage of the present invention is that it minimizes attenuation of a localization signal when filtering that signal to remove a respiration noise component.
- Still another advantage of the present invention is that it provides a method of collecting localization data at substantially the same point in each respiration cycle.
- Figure 1 is a schematic diagram of a localization system utilized in an electrophysiology study.
- Figure 2 depicts an exemplary catheter used in an electrophysiology study.
- FIGS 3a-3e illustrate the effects of respiration on localization data, wherein:
- Figure 3 a is a representative EKG signal
- Figure 3b is a representative x-axis displacement signal for a localization system patch electrode
- Figure 3c is a representative y-axis displacement signal for a localization system patch electrode
- Figure 3d is a representative z-axis displacement signal for a localization system patch electrode
- Figure 3e is a representative three-dimensional displacement signal for a localization system patch electrode.
- Figure 4 is a flowchart of a method of filtering a respiration noise component from a localization signal according to an embodiment of the present invention.
- Figures 5a-5e illustrate a method of filtering a respiration noise component from a localization signal according to another embodiment of the present invention, wherein:
- Figure 5 a is a representative EKG signal
- Figure 5b is a representative x-axis displacement signal for a localization system patch electrode
- Figure 5c is a representative y-axis displacement signal for a localization system patch electrode
- Figure 5d is a representative z-axis displacement signal for a localization system patch electrode.
- Figure 5e is a representative three-dimensional acceleration signal for a localization system patch electrode.
- the present invention provides a method and system for filtering a respiration noise component from a localization signal.
- the term "localization signal” means a signal that describes at least the position, and optionally also the orientation, of a medical device or other object within a localization field.
- the localization field will be a non-ionizing localization field, such as an electric or magnetic field.
- Figure 1 shows a schematic diagram of a localization system 8 for conducting cardiac electrophysiology studies by navigating a cardiac catheter and measuring electrical activity occurring in a heart 10 of a patient 11 and three-dimensionally mapping the electrical activity and/or information related to or representative of the electrical activity so measured.
- System 8 can be used, for example, to create an anatomical model of the patient's heart 10 using one or more electrodes.
- System 8 can also be used to measure electrophysiology data at a plurality of points along a cardiac surface, and store the measured data in association with location information for each measurement point at which the electrophysiology data was measured, for example to create a diagnostic data map of the patient's heart 10.
- localization system 8 determines the location of objects, typically within a three-dimensional space, and expresses those locations as position information determined relative to at least one reference.
- the patient 11 is depicted schematically as an oval.
- three sets of surface electrodes e.g ., patch electrodes
- a surface of the patient 11 defining three generally orthogonal axes, referred to herein as an x-axis, a y-axis, and a z-axis.
- the electrodes could be positioned in other arrangements, for example multiple electrodes on a particular body surface.
- the electrodes do not need to be on the body surface, but could be fixed on an external apparatus, or electrodes positioned internally to the body could be used.
- the x-axis surface electrodes 12, 14 are applied to the patient along a first axis, such as on the lateral sides of the thorax region of the patient ( e.g ., applied to the patient's skin underneath each arm) and may be referred to as the Left and Right electrodes.
- the y-axis electrodes 18, 19 are applied to the patient along a second axis generally orthogonal to the x-axis, such as along the inner thigh and neck regions of the patient, and may be referred to as the Left Leg and Neck electrodes.
- the z-axis electrodes 16, 22 are applied along a third axis generally orthogonal to both the x-axis and the y-axis, such as along the sternum and spine of the patient in the thorax region, and may be referred to as the Chest and Back electrodes.
- the heart 10 lies between these pairs of surface electrodes 12/14, 18/19, and 16/22.
- An additional surface reference electrode (e.g ., a "belly patch") 21 provides a reference and/or ground electrode for the system 8.
- the belly patch electrode 21 may be an alternative to a fixed intracardiac electrode 31, described in further detail below.
- the patient 11 may have most or all of the conventional electrocardiogram ("ECG" or "EKG”) system leads in place. This ECG information is available to the system 8, although not illustrated in Figure 1 .
- a representative catheter 13 having at least one electrode 17 is also shown.
- This representative catheter electrode 17 is referred to as the "roving electrode,” “moving electrode,” “measurement electrode,” or “position measurement sensor” throughout this specification.
- multiple electrodes on catheter 13, or on multiple such catheters will be used.
- localization system 8 may comprise sixty-four electrodes on twelve catheters disposed within the heart and/or vasculature of the patient.
- this embodiment is merely exemplary, and any number of electrodes and catheters may be used within the scope of the present invention.
- catheter 13 extends into the left ventricle 50 of the patient's heart 10.
- catheter 13 includes electrode 17 on its distal tip, as well as a plurality of additional position measurement sensors (e.g ., measurement electrodes) 52, 54, 56 spaced along its length.
- additional position measurement sensors e.g ., measurement electrodes
- the spacing between adjacent electrodes will be known, though it should be understood that the electrodes may not be evenly spaced along catheter 13 or of equal size to each other. Since each of these electrodes 17, 52, 54, 56 lies within the patient, location data may be collected simultaneously for each of the electrodes by localization system 8.
- each position measurement sensor may generate a localization signal that describes the position, and optionally also the orientation, of catheter 13 within the localization field generated by patch electrodes 12/14, 18/19, and 16/22.
- an optional fixed reference electrode 31 (e.g ., attached to a wall of the heart 10) is shown on a second catheter 29.
- this electrode 31 may be stationary (e.g ., attached to or near the wall of the heart) or disposed in a fixed spatial relationship with the roving electrodes ( e.g ., electrodes 17, 52, 54, 56), and thus may be referred to as a "navigational reference" or "local reference.”
- the fixed reference electrode 31 may be used in addition or alternatively to the surface reference electrode 21 described above.
- a coronary sinus electrode or other fixed electrode in the heart 10 can be used as a reference for measuring voltages and displacements; that is, as described below, fixed reference electrode 31 may define the origin of a coordinate system for the localization field.
- Each surface electrode is coupled to the multiplex switch 24, and the pairs of surface electrodes are selected by software running on a computer 20, which couples the surface electrodes to a signal generator 25.
- the computer 20, for example, may comprise a conventional general-purpose computer, a special-purpose computer, a distributed computer, or any other type of computer.
- the computer 20 may comprise one or more processors, such as a single central processing unit (CPU), or a plurality of processing units, commonly referred to as a parallel processing environment, which may execute instructions to practice the various aspects of the present invention described herein.
- three nominally orthogonal electric fields are generated by a series of driven and sensed electric dipoles (e.g ., surface electrode pairs 12/14, 18/19, and 16/22) in order to realize catheter navigation in a biological conductor.
- these orthogonal fields can be decomposed and any pairs of surface electrodes can be driven as dipoles to provide effective electrode triangulation.
- the electrodes 12, 14, 18, 19, 16, and 22 (or any number of electrodes) could be positioned in any other effective arrangement for driving a current to or sensing a current from an electrode in the heart.
- multiple electrodes could be placed on the back, sides, and/or belly of patient 11. Additionally, such non-orthogonal methodologies add to the flexibility of the system.
- the potentials measured across the roving electrodes resulting from a predetermined set of drive (source-sink) configurations may be combined algebraically to yield the same effective potential as would be obtained by simply driving a uniform current along the orthogonal axes.
- any two of the surface electrodes 12, 14, 16, 18, 19, 22 may be selected as a dipole source and drain with respect to a ground reference, such as belly patch 21, while the unexcited electrodes measure voltage with respect to the ground reference.
- the roving electrodes 17, 52, 54, 56 placed in the heart 10 are exposed to the field from a current pulse and are measured with respect to ground, such as belly patch 21.
- the catheters within the heart may contain more or fewer electrodes than the four shown, and each electrode potential may be measured.
- at least one electrode may be fixed to the interior surface of the heart to form a fixed reference electrode 31, which is also measured with respect to ground, such as belly patch 21, and which may be defined as the origin of the coordinate system relative to which localization system 8 measures positions. Data sets from each of the surface electrodes, the internal electrodes, and the virtual electrodes may all be used to determine the location of the roving electrodes 17, 52, 54, 56 within heart 10.
- the measured voltages may be used to determine the location in three-dimensional space of the electrodes inside the heart, such as roving electrodes 17, 52, 54, 56, relative to a reference location, such as reference electrode 31. That is, the voltages measured at reference electrode 31 may be used to define the origin of a coordinate system, while the voltages measured at roving electrodes 17, 52, 54, 56 may be used to express the location of roving electrodes 17, 52, 54, 56 relative to the origin.
- the invention will be described in connection with a three-dimensional (x, y, z) Cartesian coordinate system. It should be understood, however, that other coordinate systems, such as spherical and cylindrical coordinate systems in three dimensions and polar coordinate systems in two dimensions, are within the scope of the invention.
- the data used to determine the location of the electrode(s) within the heart is measured while the surface electrode pairs impress an electric field on the heart.
- the electrode data may also be used to create a respiration compensation value used to improve the raw location data for the electrode locations as described in United States patent no. 7,263,397 , which is hereby incorporated herein by reference in its entirety.
- the electrode data may also be used to compensate for changes in the impedance of the body of the patient as described in co-pending United States application no. 11/227,580, filed 15 September 2005 , which is also incorporated herein by reference in its entirety.
- the system 8 first selects a set of surface electrodes and then drives them with current pulses. While the current pulses are being delivered, electrical activity, such as the voltages measured at least one of the remaining surface electrodes and in vivo electrodes, is measured and stored. Compensation for artifacts, such as respiration and/or impedance shifting, may be performed as indicated above.
- the localization/mapping system is the EnSite NavXTM navigation and visualization system of St. Jude Medical, Atrial Fibrillation Division, Inc., which generates the electrical fields described above.
- Other localization systems may be used in connection with the present invention, including for example, the CARTO navigation and location system of Biosense Webster, Inc., the AURORA® system of Northern Digital Inc., or Sterotaxis' NIOBE® Magnetic Navigation System, all of which utilize magnetic fields rather than electrical fields.
- the localization and mapping systems described in the following patents can also be used with the present invention: United States patent nos. 6,990,370 ; 6,978,168 ; 6,947,785 ; 6,939,309 ; 6,728,562 ; 6,640,119 ; 5,983,126 ; and 5,697,377 .
- the fields generated by localization system 8, whether an electrical field (e.g ., EnSite NavXTM), a magnetic field (e.g ., CARTO, AURORA®, NIOBE®), or another suitable field, may be referred to generically as "localization fields," while the elements generating the fields, such as surface electrodes 12, 14, 16, 18, 19, and 22 may be generically referred to as “localization field generators.” As described above, surface electrodes 12, 14, 16, 18, 19, and 22 may also function as detectors to measure the characteristics of the localization field (e.g ., the voltages measured at roving electrodes 17, 52, 54, 56, or a current from roving electrodes 17, 52, 54, 56), and thus may also be referred to as "localization elements.” Though the present invention will be described primarily in the context of a localization system that generates an electrical field, one of ordinary skill in the art will understand how to apply the principles disclosed herein in other types of localization fields, and in particular other types of non-ionizing localization
- the localization signal may include artifacts (e.g ., a noise component) attributable to patient respiration.
- artifacts e.g ., a noise component
- the measured displacement error still exceeded about one centimeter. That is, patient respiration may cause roving electrodes 17, 52, 54, and 56 to appear to be as far as 2 cm away from their actual location in the three-dimensional space of the localization field.
- this electrode exhibits considerable respiration noise on the z-axis ( Figure 3d ), and much less respiration noise on the x-axis ( Figure 3b ) and y-axis ( Figure 3c ).
- United States patent no. 7,263,397 discloses a method of compensating for these respiration artifacts.
- the present invention provides methods to filter a respiratory noise component from a localization signal acquired from a position measurement sensor (e.g ., one or more of electrodes 17, 52, 54, and 56).
- the present invention employs accelerometers attached to one or more of the localization system patch electrodes (e.g ., localization field generators 12/14, 18/19, and 16/22, as well as belly patch 21) in order to acquire a signal that describes the motion of the patch electrodes.
- each of the localization system patch electrodes is equipped with an accelerometer, but the invention could be practiced in a localization system where fewer than all of the patch electrodes includes an accelerometer. It is also desirable for the accelerometers to be capable of measuring acceleration in three dimensions, though the invention could still be practiced with accelerometers having more limited capabilities.
- FIG. 4 illustrates application of this method according to the present invention in order to filter a respiration noise component from the x-axis component of a localization signal acquired from a position measurement sensor, such as electrode 17, within the localization field.
- a position measurement sensor such as electrode 17
- an acceleration signal is acquired for at least one patch electrode.
- the acceleration signal is acquired for at least one patch electrode that defines the coordinate axis on which the respiration noise component is to be filtered (e.g ., the acceleration signal from patch electrodes 12 and/or 14 will be used to filter on the x-axis; the acceleration signal from patch electrodes 18 and/or 19 will be used to filter on the y-axis; and the acceleration signal from patch electrodes 16 and/or 22 will be used to filter on the z-axis).
- the acquired acceleration signal is used to calculate a displacement signal of the at least one patch electrode on the x-axis. That is, the acquired acceleration signal is processed in order to generate a signal that describes the movement of the at least one patch electrode.
- the acquired acceleration signal is integrated twice and resolved into vector components along the x, y, and z-axes (line 101).
- the acquired acceleration signal may be resolved into its components either prior to or after double integration.
- the x-axis displacement signal is transformed from the time domain into the frequency domain in block 104. This may be accomplished, for example, by applying a Fourier transform to the x-axis displacement signal.
- the transformed signal will typically exhibit two peaks: a low frequency peak resulting from the respiration noise component, and a high frequency peak resulting from the physiological component.
- a fractional power of the transformed signal attributable to patient respiration is calculated. This may be accomplished by calculating a ratio of two integrals of the transformed x-axis displacement signal: the integral of the transformed x-axis displacement signal between about 0 Hz and about 0.3 Hz (a range sufficient to cover most, if not all, respiratory activity), and the integral of the transformed x-axis displacement signal as a whole.
- the fractional power of the transformed signal attributable to patient respiration is compared to a threshold value 0, discussed in further detail below. If the fractional power is less than the threshold value, then there is not a significant respiration noise component on the x-axis, and no further action needs to be taken (block 110).
- the fractional power exceeds the threshold value, then there is a significant respiration noise component on the x-axis, and it would be desirable to filter this component out of the x-axis component of the localization signal.
- this is accomplished by applying a high-pass filter (block 112) to the x-axis component of the localization signal (block 114), thereby outputting a filtered x-axis component of the localization signal (block 116).
- a high-pass filter is used in order to filter out the relatively low frequency respiration noise component while preserving the relatively high frequency physiological component.
- the cutoff frequency of the high-pass filter should be set above the patient's respiratory frequency but below the frequency of the physiological signal.
- the threshold value ⁇ will typically be selected on an individualized basis for each patient.
- ⁇ is computed using patch electrode displacement data for a localization system coordinate axis that does not exhibit significant respiratory noise (referred to herein as a "quiet" axis).
- Figures 3b and 3c show that the depicted patch electrode is experiencing very little respiratory noise on the x- and y-axes. Therefore, either the x- or y-axis could be used to determine a suitable 0.
- ⁇ will be computed by transforming the displacement signal for the quiet axis into the frequency domain and selecting a value that falls above the lower frequency peak of the respiration noise signal and below the higher frequency peak of the physiological signal. As a matter of convenience, this same value may also be used as the cutoff frequency for the high pass filter.
- the threshold value, cutoff frequency, and/or patient's respiration frequency may be determined using any pseudo-stationary reference electrode.
- the term "pseudo-stationary" refers to an electrode that is not intended to be moving, although it may appear to be moving due to external influences.
- patch electrodes 12/14, 18/19, 16/22, and 21, as well as fixed intracardiac electrode 31, are all considered pseudo-stationary - they are supposed to remain stationary to provide a fixed frame of reference for the localization system. Electrodes 17, 52, 54, and 56, however, are not pseudo-stationary, as they are intended to move through the patient's heart.
- a localization signal of intracardiac electrode 31 is used to determine the threshold value, cutoff frequency, and/or patient's respiration frequency.
- this is a localization signal, it already represents displacement and need not be integrated before being transformed into the frequency domain.
- an acceleration signal from a patch electrode is used to determine the threshold value, cutoff frequency, and/or patient's respiration frequency. This signal may need to be integrated twice to compute a displacement signal prior to further processing.
- the collection of localization data is gated to the patient's respiratory cycle, using the acquired acceleration signal, such that all data is collected at the same fiducial point in each respiration cycle (e.g ., data is always collected at peak inspiration or always collected at peak expiration).
- Figures 5b-5d depict, respectively, patch electrode motion on the x-, y-, and z-axes.
- Figure 5e is the acceleration signal for the patch electrode.
- Figures 5b-5e are shown juxtaposed with an EKG signal ( Figure 5a ).
- the darker shaded bands 60 represent patient inspiration, while the lighter shaded bands 70 represent patient expiration.
- the acquired acceleration signal is integrated once to yield a velocity signal, and data points are collected only when the velocity signal on all coordinate axes changes sign. For example, when the velocity signals change sign from positive to negative, it corresponds to peak inspiration, when the patient begins to exhale. Conversely, when the velocity signals change sign from negative to positive, it corresponds to peak expiration, when the patient begins to inhale.
- the acquired acceleration signal is integrated twice to yield a displacement signal, and data points are collected only when the displacement signal on all coordinate axes is at an extreme. For example, when the displacement signals are at their maxima, it corresponds to peak inspiration, when the patient begins to exhale. Conversely, when the displacement signals are at their minima, it corresponds to peak expiration, when the patient beings to inhale.
- the methods described above may be executed by one or more computer systems (e.g ., computer system 20), and may be software implemented (e.g ., one or more software programs executed by one or more computer systems or processors), hardware implemented (e.g ., a series of instructions stored in one or more solid state devices), or a combination of both.
- the computer may be a conventional general purpose computer, a special purpose computer, a distributed computer, or any other type of computer. Further, the computer may comprise one or more processors, such as a single central processing unit or a plurality of processing units, commonly referred to as a parallel processing environment.
- processor refers to a computer microprocessor and/or a software program (e.g ., a software module or separate program) that is designed to be executed by one or more microprocessors running on one or more computer systems.
- a localization system for measuring a position of a measurement electrode within a non-ionizing localization field may include: a plurality of localization field generators capable of generating a non-ionizing localization field (e.g ., patch electrodes 12/14, 18/19, and 16/22); at least one localization field reference (e.g ., belly patch 21 and/or intracardiac electrode 31); at least one accelerometer that is capable of measuring acceleration in at least one direction coupled to at least one of the plurality of localization field generators and the at least one localization field reference; and a signal processor that acquires at least one acceleration signal from the at least one accelerometer and uses that acquired acceleration signal to filter a respiration noise component from a localization signal received from at least one position measurement sensor (e.g ., electrodes 17, 52, 54, and 56) within the localization field.
- a plurality of localization field generators capable of generating a non-ionizing localization field (e.g ., patch electrodes 12/14,
- the methods disclosed herein can be employed to good advantage in various contexts, including cardiac therapeutic contexts.
- various vector-based metrics such as the variation in distance between electrodes in various heart chambers, in connection with cardiac pacing and cardiac resynchronization therapy, such as disclosed in United States application nos. 12/621,373, filed 18 November 2009 , and 12/621,397, filed 18 November 2009 , and United States provisional application no. 61/121,737, filed 11 December 2008 , all of which are hereby incorporated by reference as though fully set forth herein.
- the present teachings can aid in minimizing the effect of respiratory oscillations in determining these desirable vector-based metrics.
- a slow-moving average filter may be applied to the displacement signal in order to smooth out the cardiac component and determine the peak-to-peak interval.
- the peak-to-peak interval will, in turn, facilitate determination of the respiration cycle, which can be used to determine the appropriate cutoff frequency for the filter (e.g ., a high-pass filter).
- All directional references e.g ., upper, lower, upward, downward, left, right, leftward, rightward, top, bottom, above, below, vertical, horizontal, clockwise, and counterclockwise
- Joinder references e.g., attached, coupled, connected, and the like
- Joinder references are to be construed broadly and may include intermediate members between a connection of elements and relative movement between elements. As such, joinder references do not necessarily infer that two elements are directly connected and in fixed relation to each other.
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Surgery (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Pathology (AREA)
- Physics & Mathematics (AREA)
- Signal Processing (AREA)
- Biophysics (AREA)
- Human Computer Interaction (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physiology (AREA)
- Psychiatry (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
Description
- The instant invention relates to the collection of data regarding a patient, such as geometric data regarding a portion of the patient's anatomy, using a localization system. In particular, the instant invention relates to methods and systems for processing a localization signal that describes the position, and optionally the orientation, of a device within the localization field to account for signal artifacts attributable to patient respiration.
- It is well known how to generate heart chamber geometry in preparation for cardiac diagnostic or therapeutic procedures. Often, a mapping catheter is introduced into the heart chamber of interest and moved around within the heart chamber, either randomly, pseudo-randomly, or according to one or more preset patterns. The three-dimensional coordinates are measured using a localization system (sometimes also referred to as a "mapping system," "navigation system," or "positional feedback system"). The localization system measures the coordinates of the mapping catheter within a localization field, typically by relating a characteristic of the localization field, such as a voltage, experienced by the mapping catheter to a location of the catheter within the field.
- Localization signals produced by localization systems often exhibit noise. One prevalent source of noise is patient respiration. Patient respiration may introduce noise by moving the generators and/or references for the localization field or by changing the patient's biological characteristics (e.g., impedance). Thus, an electrode within the localization field that is, in fact, stationary may appear to be moving.
- It is desirable to be able to remove the respiration noise component from a localization signal generated by a localization system.
- It is also desirable to be able to remove the respiration noise component from a localization signal without attenuating that signal.
- Another object of the present invention is to provide a method of collecting localization data that ensures the collected data relates to a common fiducial point within a patient's respiration cycle, such as peak inspiration or peak expiration.
- Disclosed herein is a method of filtering a respiration noise component from a localization signal of a localization system. The localization system includes a plurality of localization system patch electrodes that generate a non-ionizing localization field, and at least one of the localization system patch electrodes includes an accelerometer capable of measuring acceleration in at least one direction. The method includes the following steps: acquiring a localization signal of at least one position measurement sensor within the non-ionizing localization field; and acquiring an acceleration signal for the at least one localization system patch electrode; for a selected coordinate axis of the non-ionizing localization field: calculating a displacement signal of the at least one localization system patch electrode on the selected coordinate axis from the acceleration signal; transforming the displacement signal from the time domain into the frequency domain; calculating a fractional power of the transformed displacement signal indicative of patient respiration; comparing the calculated fractional power to a threshold value; and if the calculated fractional power exceeds the threshold value, applying a filter to the localization signal of the at least one position measurement sensor to filter out the respiration noise component on the selected coordinate axis. Optionally, the steps of calculating a displacement signal of the at least one localization system patch electrode on the selected coordinate axis from the acceleration signal, transforming the displacement signal from the time domain into the frequency domain, calculating a fractional power of the transformed displacement signal indicative of patient respiration, comparing the calculated fractional power to a threshold value, and if the calculated fractional power exceeds the threshold value, applying a filter to the localization signal of the at least one position measurement sensor to filter out the respiration noise component on the selected coordinate axis may be repeated for each coordinate axis of the localization field.
- In some embodiments of the invention, the step of calculating a fractional power of the transformed displacement signal indicative of patient respiration includes calculating a ratio of an area under the transformed displacement signal between about 0 Hz and about 0.3 Hz to an area under the transformed displacement signal as a whole.
- In some embodiments, the at least one localization system patch electrode includes a localization field generator that defines the selected coordinate axis.
- Te step of calculating a displacement signal of the at least one localization system patch electrode on the selected coordinate axis from the acceleration signal may include: resolving a component of the acceleration signal on the selected coordinate axis; and integrating the component of the acceleration signal on the selected coordinate axis twice.
- A Fourier transform may be applied to the displacement signal in order to transform the displacement signal from the time domain into the frequency domain.
- Typically, the step of applying a filter to the localization signal of the at least one position measurement sensor to filter out the respiration noise component on the selected coordinate axis includes applying a high-pass filter to the localization signal. Preferably, the cutoff frequency for the high-pass filter is selected to be above the patient's respiration frequency. For example, the threshold value may be selected as the cutoff frequency. Alternatively, the cutoff frequency may be determined using a localization signal acquired from a pseudo-stationary reference electrode positioned within the localization field. In still other embodiments of the invention, the cutoff frequency is determined using the acquired acceleration signal.
- The patient's respiration frequency may be identified by: acquiring a displacement signal for at least one pseudo-stationary electrode on at least one coordinate axis of the localization field; transforming the acquired displacement signal from the time domain to the frequency domain; and identifying the patient respiration frequency based upon a low frequency peak of the transformed signal.
- In another aspect, the present invention includes a method of filtering a respiration noise component from a localization signal of a localization system that uses a plurality of localization system patch electrodes to generate a non-ionizing localization field, wherein at least one of the patch electrodes includes an accelerometer capable of measuring acceleration in at least one direction. The method includes the following steps: acquiring an acceleration signal for the at least one localization system patch electrode; and gating collection of data points from a localization signal of at least one position measurement sensor within a localization field based upon the acquired acceleration signal.
- In some aspects of the invention, the step of gating collection of data points from a localization signal includes: integrating the acquired acceleration signal to calculate a velocity signal and collecting data points from the localization signal only when the calculated velocity signal changes sign from positive to negative. In other aspects of the invention, the gating step includes: integrating the acquired acceleration signal to calculate a velocity signal and collecting data points from the localization signal only when the calculated velocity signal changes sign from negative to positive. In still other aspects of the invention, the step of gating collection of data points from a localization signal includes: integrating the acquired localization signal twice to calculate a displacement signal and collecting data points from the localization signal only when the calculated displacement signal is at a maximum value. In yet further aspects of the invention, the step of gating collection of data points from a localization signal includes: integrating the acquired localization signal twice to calculate a displacement signal and collecting data points from the localization signal only when the calculated displacement signal is at a minimum value.
- Also disclosed herein is a localization system for measuring a position of a measurement electrode within a non-ionizing localization field. The localization system includes: a plurality of localization field generators capable of generating a non-ionizing localization field; at least one localization field reference device; at least one accelerometer capable of measuring acceleration in at least one direction coupled to at least one of the plurality of localization field generators and the at least one localization field reference device; and a signal processor that acquires at least one acceleration signal from the at least one accelerometer and that uses the acquired at least one acceleration signal to filter a respiration noise component from a localization signal of at least one position measurement sensor within the non-ionizing localization field.
- In certain embodiments of the invention, the signal processor integrates the at least one acceleration signal twice to compute at least one displacement signal and uses the computed at least one displacement signal to filter the respiration noise component from the localization signal. The signal processor may also use a high-pass filter having a cutoff frequency determined using the acquired at least one acceleration signal to process the localization signal. Alternatively, or additionally, the signal processor may process the localization signal by gating collection of data points from the localization signal using the acquired at least one acceleration signal.
- An advantage of the present invention is that it permits the filtering of a respiration noise component from a localization signal.
- Another advantage of the present invention is that it minimizes attenuation of a localization signal when filtering that signal to remove a respiration noise component.
- Still another advantage of the present invention is that it provides a method of collecting localization data at substantially the same point in each respiration cycle.
- The foregoing and other aspects, features, details, utilities, and advantages of the present invention will be apparent from reading the following description and claims, and from reviewing the accompanying drawings.
-
Figure 1 is a schematic diagram of a localization system utilized in an electrophysiology study. -
Figure 2 depicts an exemplary catheter used in an electrophysiology study. -
Figures 3a-3e illustrate the effects of respiration on localization data, wherein: -
Figure 3 a is a representative EKG signal; -
Figure 3b is a representative x-axis displacement signal for a localization system patch electrode; -
Figure 3c is a representative y-axis displacement signal for a localization system patch electrode; -
Figure 3d is a representative z-axis displacement signal for a localization system patch electrode; and -
Figure 3e is a representative three-dimensional displacement signal for a localization system patch electrode. -
Figure 4 is a flowchart of a method of filtering a respiration noise component from a localization signal according to an embodiment of the present invention. -
Figures 5a-5e illustrate a method of filtering a respiration noise component from a localization signal according to another embodiment of the present invention, wherein: -
Figure 5 a is a representative EKG signal; -
Figure 5b is a representative x-axis displacement signal for a localization system patch electrode; -
Figure 5c is a representative y-axis displacement signal for a localization system patch electrode; -
Figure 5d is a representative z-axis displacement signal for a localization system patch electrode; and -
Figure 5e is a representative three-dimensional acceleration signal for a localization system patch electrode. - The present invention provides a method and system for filtering a respiration noise component from a localization signal. As used herein, the term "localization signal" means a signal that describes at least the position, and optionally also the orientation, of a medical device or other object within a localization field. Typically, the localization field will be a non-ionizing localization field, such as an electric or magnetic field.
- Such localization fields are often employed in procedures carried out within a human body, and in particular in cardiac diagnostic and therapeutic procedures. Therefore, for purposes of illustration, the invention will be described in detail in the context of a localization system utilized in a cardiac electrophysiology procedure. It is contemplated, however, that the present invention may be practiced to good advantage in other contexts.
-
Figure 1 shows a schematic diagram of a localization system 8 for conducting cardiac electrophysiology studies by navigating a cardiac catheter and measuring electrical activity occurring in aheart 10 of a patient 11 and three-dimensionally mapping the electrical activity and/or information related to or representative of the electrical activity so measured. System 8 can be used, for example, to create an anatomical model of the patient'sheart 10 using one or more electrodes. System 8 can also be used to measure electrophysiology data at a plurality of points along a cardiac surface, and store the measured data in association with location information for each measurement point at which the electrophysiology data was measured, for example to create a diagnostic data map of the patient'sheart 10. As one of ordinary skill in the art will recognize, and as will be further described below, localization system 8 determines the location of objects, typically within a three-dimensional space, and expresses those locations as position information determined relative to at least one reference. - For simplicity of illustration, the patient 11 is depicted schematically as an oval. In the embodiment shown in
Figure 1 , three sets of surface electrodes (e.g., patch electrodes) are shown applied to a surface of the patient 11, defining three generally orthogonal axes, referred to herein as an x-axis, a y-axis, and a z-axis. In other embodiments the electrodes could be positioned in other arrangements, for example multiple electrodes on a particular body surface. Likewise, the electrodes do not need to be on the body surface, but could be fixed on an external apparatus, or electrodes positioned internally to the body could be used. - In
Figure 1 , thex-axis surface electrodes axis electrodes axis electrodes heart 10 lies between these pairs ofsurface electrodes 12/14, 18/19, and 16/22. - An additional surface reference electrode (e.g., a "belly patch") 21 provides a reference and/or ground electrode for the system 8. The
belly patch electrode 21 may be an alternative to a fixedintracardiac electrode 31, described in further detail below. It should also be appreciated that, in addition, the patient 11 may have most or all of the conventional electrocardiogram ("ECG" or "EKG") system leads in place. This ECG information is available to the system 8, although not illustrated inFigure 1 . - A
representative catheter 13 having at least one electrode 17 (e.g., a distal electrode) is also shown. Thisrepresentative catheter electrode 17 is referred to as the "roving electrode," "moving electrode," "measurement electrode," or "position measurement sensor" throughout this specification. Typically, multiple electrodes oncatheter 13, or on multiple such catheters, will be used. In one embodiment, for example, localization system 8 may comprise sixty-four electrodes on twelve catheters disposed within the heart and/or vasculature of the patient. Of course, this embodiment is merely exemplary, and any number of electrodes and catheters may be used within the scope of the present invention. - For purposes of this disclosure, an
exemplary catheter 13 is shown inFigure 2 . InFigure 2 ,catheter 13 extends into theleft ventricle 50 of the patient'sheart 10.Catheter 13 includeselectrode 17 on its distal tip, as well as a plurality of additional position measurement sensors (e.g., measurement electrodes) 52, 54, 56 spaced along its length. Typically, the spacing between adjacent electrodes will be known, though it should be understood that the electrodes may not be evenly spaced alongcatheter 13 or of equal size to each other. Since each of theseelectrodes electrodes catheter 13 within the localization field generated bypatch electrodes 12/14, 18/19, and 16/22. - Returning now to
Figure 1 , an optional fixed reference electrode 31 (e.g., attached to a wall of the heart 10) is shown on asecond catheter 29. For calibration purposes, thiselectrode 31 may be stationary (e.g., attached to or near the wall of the heart) or disposed in a fixed spatial relationship with the roving electrodes (e.g.,electrodes reference electrode 31 may be used in addition or alternatively to thesurface reference electrode 21 described above. In many instances, a coronary sinus electrode or other fixed electrode in theheart 10 can be used as a reference for measuring voltages and displacements; that is, as described below, fixedreference electrode 31 may define the origin of a coordinate system for the localization field. - Each surface electrode is coupled to the
multiplex switch 24, and the pairs of surface electrodes are selected by software running on acomputer 20, which couples the surface electrodes to asignal generator 25. Thecomputer 20, for example, may comprise a conventional general-purpose computer, a special-purpose computer, a distributed computer, or any other type of computer. Thecomputer 20 may comprise one or more processors, such as a single central processing unit (CPU), or a plurality of processing units, commonly referred to as a parallel processing environment, which may execute instructions to practice the various aspects of the present invention described herein. - Generally, three nominally orthogonal electric fields are generated by a series of driven and sensed electric dipoles (e.g., surface electrode pairs 12/14, 18/19, and 16/22) in order to realize catheter navigation in a biological conductor. Alternatively, these orthogonal fields can be decomposed and any pairs of surface electrodes can be driven as dipoles to provide effective electrode triangulation. Likewise, the
electrodes - Thus, any two of the
surface electrodes belly patch 21, while the unexcited electrodes measure voltage with respect to the ground reference. Theroving electrodes heart 10 are exposed to the field from a current pulse and are measured with respect to ground, such asbelly patch 21. In practice the catheters within the heart may contain more or fewer electrodes than the four shown, and each electrode potential may be measured. As previously noted, at least one electrode may be fixed to the interior surface of the heart to form a fixedreference electrode 31, which is also measured with respect to ground, such asbelly patch 21, and which may be defined as the origin of the coordinate system relative to which localization system 8 measures positions. Data sets from each of the surface electrodes, the internal electrodes, and the virtual electrodes may all be used to determine the location of theroving electrodes heart 10. - The measured voltages may be used to determine the location in three-dimensional space of the electrodes inside the heart, such as
roving electrodes reference electrode 31. That is, the voltages measured atreference electrode 31 may be used to define the origin of a coordinate system, while the voltages measured atroving electrodes roving electrodes - As should be clear from the foregoing discussion, the data used to determine the location of the electrode(s) within the heart is measured while the surface electrode pairs impress an electric field on the heart. The electrode data may also be used to create a respiration compensation value used to improve the raw location data for the electrode locations as described in United States patent no.
7,263,397 , which is hereby incorporated herein by reference in its entirety. The electrode data may also be used to compensate for changes in the impedance of the body of the patient as described in co-pending United States application no.11/227,580, filed 15 September 2005 - In summary, the system 8 first selects a set of surface electrodes and then drives them with current pulses. While the current pulses are being delivered, electrical activity, such as the voltages measured at least one of the remaining surface electrodes and in vivo electrodes, is measured and stored. Compensation for artifacts, such as respiration and/or impedance shifting, may be performed as indicated above.
- In a preferred embodiment, the localization/mapping system is the EnSite NavX™ navigation and visualization system of St. Jude Medical, Atrial Fibrillation Division, Inc., which generates the electrical fields described above. Other localization systems, however, may be used in connection with the present invention, including for example, the CARTO navigation and location system of Biosense Webster, Inc., the AURORA® system of Northern Digital Inc., or Sterotaxis' NIOBE® Magnetic Navigation System, all of which utilize magnetic fields rather than electrical fields. The localization and mapping systems described in the following patents (all of which are hereby incorporated by reference in their entireties) can also be used with the present invention: United States patent nos.
6,990,370 ;6,978,168 ;6,947,785 ;6,939,309 ;6,728,562 ;6,640,119 ;5,983,126 ; and5,697,377 . - The fields generated by localization system 8, whether an electrical field (e.g., EnSite NavX™), a magnetic field (e.g., CARTO, AURORA®, NIOBE®), or another suitable field, may be referred to generically as "localization fields," while the elements generating the fields, such as
surface electrodes surface electrodes roving electrodes roving electrodes electrodes - The basic localization methodology described above provides a first order indication of the location of a roving or other catheter electrode within the heart chamber. As one of ordinary skill in the art will appreciate, however, the localization signal may include artifacts (e.g., a noise component) attributable to patient respiration. For example, when measured with reference to the
belly patch 21, displacement errors exceeding two centimeters have been noted in the left atrium in data measured with aroving catheter electrode 17. When referenced to a fixed electrode, such asintracardiac electrode 31, the measured displacement error still exceeded about one centimeter. That is, patient respiration may causeroving electrodes - These errors may be caused by several different aspects of patient respiration. First, when the lungs fill with air, the current paths from the localization field generators change. This alteration of current path changes the measured potential between
roving electrode 17 and the reference electrode (e.g., belly patch 31). Respiration also moves the localization system patch electrodes (e.g.,localization field generators 12/14, 18/19, 16/22 and belly patch 31); this displacement may also change the apparent position of an object within the localization field. The resistivity of blood returning to the heart from the lungs is also altered following inspiration. Full lungs may also have a compressive effect on the heart that may cause displacement of the catheter itself. - This potential for error in the localization signal will be apparent to one of ordinary skill in the art from reviewing
Figures 3a-e , which show the movement of a localization field generator for the z-axis (e.g.,patch electrode 16 or 22) on the x, y, and z-axes (Figures 3b, 3c, and 3d , respectively), as well as in three dimensions (Figure 3e ). For the sake of understanding, the movement of the patch electrode is juxtaposed with the signal received by an EKG lead (Figure 3a ). As will be understood fromFigures 3b-3d , this electrode exhibits considerable respiration noise on the z-axis (Figure 3d ), and much less respiration noise on the x-axis (Figure 3b ) and y-axis (Figure 3c ). - As noted above, United States patent no.
7,263,397 discloses a method of compensating for these respiration artifacts. As an alternative or in addition to this respiration compensation methodology, the present invention provides methods to filter a respiratory noise component from a localization signal acquired from a position measurement sensor (e.g., one or more ofelectrodes - The present invention employs accelerometers attached to one or more of the localization system patch electrodes (e.g.,
localization field generators 12/14, 18/19, and 16/22, as well as belly patch 21) in order to acquire a signal that describes the motion of the patch electrodes. Preferably, each of the localization system patch electrodes is equipped with an accelerometer, but the invention could be practiced in a localization system where fewer than all of the patch electrodes includes an accelerometer. It is also desirable for the accelerometers to be capable of measuring acceleration in three dimensions, though the invention could still be practiced with accelerometers having more limited capabilities. - A first method of filtering a respiration noise component from a localization signal will be described with reference to the flowchart of
Figure 4. Figure 4 illustrates application of this method according to the present invention in order to filter a respiration noise component from the x-axis component of a localization signal acquired from a position measurement sensor, such aselectrode 17, within the localization field. One of ordinary skill in the art will appreciate how to extend the method ofFigure 4 , as described herein, to other coordinate axes defining the localization field. - In
block 100, an acceleration signal is acquired for at least one patch electrode. Preferably, the acceleration signal is acquired for at least one patch electrode that defines the coordinate axis on which the respiration noise component is to be filtered (e.g., the acceleration signal frompatch electrodes 12 and/or 14 will be used to filter on the x-axis; the acceleration signal frompatch electrodes 18 and/or 19 will be used to filter on the y-axis; and the acceleration signal frompatch electrodes 16 and/or 22 will be used to filter on the z-axis). - In
block 102, the acquired acceleration signal is used to calculate a displacement signal of the at least one patch electrode on the x-axis. That is, the acquired acceleration signal is processed in order to generate a signal that describes the movement of the at least one patch electrode. Thus, in moving fromblock 100 to block 102, the acquired acceleration signal is integrated twice and resolved into vector components along the x, y, and z-axes (line 101). Of course, the acquired acceleration signal may be resolved into its components either prior to or after double integration. - Once the x-axis displacement signal has been calculated, it is transformed from the time domain into the frequency domain in
block 104. This may be accomplished, for example, by applying a Fourier transform to the x-axis displacement signal. The transformed signal will typically exhibit two peaks: a low frequency peak resulting from the respiration noise component, and a high frequency peak resulting from the physiological component. - In
block 106, a fractional power of the transformed signal attributable to patient respiration is calculated. This may be accomplished by calculating a ratio of two integrals of the transformed x-axis displacement signal: the integral of the transformed x-axis displacement signal between about 0 Hz and about 0.3 Hz (a range sufficient to cover most, if not all, respiratory activity), and the integral of the transformed x-axis displacement signal as a whole. - In
decision block 108, the fractional power of the transformed signal attributable to patient respiration is compared to athreshold value 0, discussed in further detail below. If the fractional power is less than the threshold value, then there is not a significant respiration noise component on the x-axis, and no further action needs to be taken (block 110). - If, on the other hand, the fractional power exceeds the threshold value, then there is a significant respiration noise component on the x-axis, and it would be desirable to filter this component out of the x-axis component of the localization signal. Preferably, this is accomplished by applying a high-pass filter (block 112) to the x-axis component of the localization signal (block 114), thereby outputting a filtered x-axis component of the localization signal (block 116). A high-pass filter is used in order to filter out the relatively low frequency respiration noise component while preserving the relatively high frequency physiological component. Thus, the cutoff frequency of the high-pass filter should be set above the patient's respiratory frequency but below the frequency of the physiological signal.
- The threshold value θ will typically be selected on an individualized basis for each patient. Preferably, θ is computed using patch electrode displacement data for a localization system coordinate axis that does not exhibit significant respiratory noise (referred to herein as a "quiet" axis). For example,
Figures 3b and 3c show that the depicted patch electrode is experiencing very little respiratory noise on the x- and y-axes. Therefore, either the x- or y-axis could be used to determine a suitable 0. Typically, θ will be computed by transforming the displacement signal for the quiet axis into the frequency domain and selecting a value that falls above the lower frequency peak of the respiration noise signal and below the higher frequency peak of the physiological signal. As a matter of convenience, this same value may also be used as the cutoff frequency for the high pass filter. - It should be understood that the threshold value, cutoff frequency, and/or patient's respiration frequency may be determined using any pseudo-stationary reference electrode. As used herein, the term "pseudo-stationary" refers to an electrode that is not intended to be moving, although it may appear to be moving due to external influences. Thus,
patch electrodes 12/14, 18/19, 16/22, and 21, as well as fixedintracardiac electrode 31, are all considered pseudo-stationary - they are supposed to remain stationary to provide a fixed frame of reference for the localization system.Electrodes - Thus, in some embodiments of the invention, a localization signal of
intracardiac electrode 31 is used to determine the threshold value, cutoff frequency, and/or patient's respiration frequency. Of course, because this is a localization signal, it already represents displacement and need not be integrated before being transformed into the frequency domain. - In other embodiments of the invention, an acceleration signal from a patch electrode is used to determine the threshold value, cutoff frequency, and/or patient's respiration frequency. This signal may need to be integrated twice to compute a displacement signal prior to further processing.
- Another embodiment of the present invention will now be described with reference to
Figures 5a-5e . In this embodiment of the invention, the collection of localization data is gated to the patient's respiratory cycle, using the acquired acceleration signal, such that all data is collected at the same fiducial point in each respiration cycle (e.g., data is always collected at peak inspiration or always collected at peak expiration). -
Figures 5b-5d depict, respectively, patch electrode motion on the x-, y-, and z-axes.Figure 5e is the acceleration signal for the patch electrode. For the sake of comparison,Figures 5b-5e are shown juxtaposed with an EKG signal (Figure 5a ). Throughout, the darkershaded bands 60 represent patient inspiration, while the lightershaded bands 70 represent patient expiration. - There are many ways in which data collection may be gated. In some aspects of the invention, the acquired acceleration signal is integrated once to yield a velocity signal, and data points are collected only when the velocity signal on all coordinate axes changes sign. For example, when the velocity signals change sign from positive to negative, it corresponds to peak inspiration, when the patient begins to exhale. Conversely, when the velocity signals change sign from negative to positive, it corresponds to peak expiration, when the patient begins to inhale.
- In other aspects of the invention, the acquired acceleration signal is integrated twice to yield a displacement signal, and data points are collected only when the displacement signal on all coordinate axes is at an extreme. For example, when the displacement signals are at their maxima, it corresponds to peak inspiration, when the patient begins to exhale. Conversely, when the displacement signals are at their minima, it corresponds to peak expiration, when the patient beings to inhale.
- By gating the collection of data points as described above, a "fair playing field" is ensured. It should be understood that, although data points are collected on a gated basis (i.e., at the same fiducial point in each respiration cycle), the localization signal may nonetheless be constantly tracked. That is, the location of
electrodes - The methods described above may be executed by one or more computer systems (e.g., computer system 20), and may be software implemented (e.g., one or more software programs executed by one or more computer systems or processors), hardware implemented (e.g., a series of instructions stored in one or more solid state devices), or a combination of both. The computer may be a conventional general purpose computer, a special purpose computer, a distributed computer, or any other type of computer. Further, the computer may comprise one or more processors, such as a single central processing unit or a plurality of processing units, commonly referred to as a parallel processing environment. The term "processor" as used herein refers to a computer microprocessor and/or a software program (e.g., a software module or separate program) that is designed to be executed by one or more microprocessors running on one or more computer systems.
- For example, a localization system for measuring a position of a measurement electrode within a non-ionizing localization field may include: a plurality of localization field generators capable of generating a non-ionizing localization field (e.g.,
patch electrodes 12/14, 18/19, and 16/22); at least one localization field reference (e.g.,belly patch 21 and/or intracardiac electrode 31); at least one accelerometer that is capable of measuring acceleration in at least one direction coupled to at least one of the plurality of localization field generators and the at least one localization field reference; and a signal processor that acquires at least one acceleration signal from the at least one accelerometer and uses that acquired acceleration signal to filter a respiration noise component from a localization signal received from at least one position measurement sensor (e.g.,electrodes - The methods disclosed herein can be employed to good advantage in various contexts, including cardiac therapeutic contexts. For example, it has been proposed to use various vector-based metrics, such as the variation in distance between electrodes in various heart chambers, in connection with cardiac pacing and cardiac resynchronization therapy, such as disclosed in United States application nos.
12/621,373, filed 18 November 2009 12/621,397, filed 18 November 2009 61/121,737, filed 11 December 2008 - Although several embodiments of this invention have been described above with a certain degree of particularity, those skilled in the art could make numerous alterations to the disclosed embodiments without departing from the spirit or scope of this invention. For example, rather than transforming the displacement signal from the time domain into the frequency domain, a slow-moving average filter may be applied to the displacement signal in order to smooth out the cardiac component and determine the peak-to-peak interval. The peak-to-peak interval will, in turn, facilitate determination of the respiration cycle, which can be used to determine the appropriate cutoff frequency for the filter (e.g., a high-pass filter).
- All directional references (e.g., upper, lower, upward, downward, left, right, leftward, rightward, top, bottom, above, below, vertical, horizontal, clockwise, and counterclockwise) are only used for identification purposes to aid the reader's understanding of the present invention, and do not create limitations, particularly as to the position, orientation, or use of the invention. Joinder references (e.g., attached, coupled, connected, and the like) are to be construed broadly and may include intermediate members between a connection of elements and relative movement between elements. As such, joinder references do not necessarily infer that two elements are directly connected and in fixed relation to each other.
- It is intended that all matter contained in the above description or shown in the accompanying drawings shall be interpreted as illustrative only and not limiting. Changes in detail or structure may be made without departing from the spirit of the invention as defined in the appended claims.
- It is explicitly stated that all features disclosed in the description and/or the claims are intended to be disclosed separately and independently from each other for the purpose of original disclosure as well as for the purpose of restricting the claimed invention independent of the composition of the features in the embodiments and/or the claims. It is explicitly stated that all value ranges or indications of groups of entities disclose every possible intermediate value or intermediate entity for the purpose of original disclosure as well as for the purpose of restricting the claimed invention, in particular as limits of value ranges.
Claims (15)
- A method of filtering a respiration noise component from a localization signal of a localization system including a plurality of localization system patch electrodes that generate a non-ionizing localization field, wherein at least one of the localization system patch electrodes includes an accelerometer capable of measuring acceleration in at least one direction, the method comprising:acquiring a localization signal of at least one position measurement sensor within the non-ionizing localization field;acquiring an acceleration signal for the at least one localization system patch electrode; andprocessing the localization signal of the at least one position measurement sensor using the acquired acceleration signal.
- The method of claim 1, the method further comprising:for a selected coordinate axis of the non-ionizing localization field:calculating a displacement signal of the at least one localization system patch electrode on the selected coordinate axis from the acceleration signal;transforming the displacement signal from the time domain into the frequency domain;calculating a fractional power of the transformed displacement signal indicative of patient respiration;comparing the calculated fractional power to a threshold value; andif the calculated fractional power exceeds the threshold value, applying a filter to the localization signal of the at least one position measurement sensor to filter out the respiration noise component on the selected coordinate axis.
- The method according to claim 2, wherein the step of calculating a fractional power of the transformed displacement signal indicative of patient respiration comprises calculating a ratio of an area under the transformed displacement signal between about 0 Hz and about 0.3 Hz to an area under the transformed displacement signal as a whole.
- The method according to any one of claims 2 to 3, further comprising repeating the steps of:calculating a displacement signal of the at least one localization system patch electrode on the selected coordinate axis from the acceleration signal;transforming the displacement signal from the time domain into the frequency domain;calculating a fractional power of the transformed displacement signal indicative of patient respiration;comparing the calculated fractional power to a threshold value; andif the calculated fractional power exceeds the threshold value, applying a filter to the localization signal of the at least one position measurement sensor to filter out the respiration noise component on the selected coordinate axisfor each coordinate axis of the localization field.
- The method according to any one of claims 2 to 4, wherein the at least one localization system patch electrode comprises a localization field generator that defines the selected coordinate axis.
- The method according to any one of claims 2 to 5, wherein the step of calculating a displacement signal of the at least one localization system patch electrode on the selected coordinate axis from the acceleration signal comprises:resolving a component of the acceleration signal on the selected coordinate axis; andintegrating the component of the acceleration signal on the selected coordinate axis twice.
- The method according to any one of claims 2 to 6, wherein the step of transforming the displacement signal from the time domain into the frequency domain comprises applying a Fourier transform to the displacement signal.
- The method according to any one of claims 2 to 7, wherein the step of applying a filter to the localization signal of the at least one position measurement sensor to filter out the respiration noise component on the selected coordinate axis comprises applying a high-pass filter to the localization signal, wherein a cutoff frequency for the high-pass filter is selected to be above a patient respiration frequency.
- The method according to claim 8, wherein the threshold value is selected as the cutoff frequency; or
wherein the cutoff frequency is determined using a localization signal acquired from a pseudo-stationary reference electrode positioned within the localization field; or
wherein the cutoff frequency is determined using the acquired acceleration signal. - The method according to claim 8, wherein the patient respiration frequency is identified by:acquiring a displacement signal for at least one pseudo-stationary electrode on at least one coordinate axis of the localization field;transforming the acquired displacement signal from the time domain to the frequency domain; andidentifying the patient respiration frequency based upon a low frequency peak of the transformed signal.
- The method of claim 1, the method further comprising:gating collection of data points from a localization signal of at least one position measurement sensor within a localization field based upon the acquired acceleration signal.
- The method according to claim 11, wherein the step of gating collection of data points from a localization signal comprises:integrating the acquired acceleration signal to calculate a velocity signal;collecting data points from the localization signal only when the calculated velocity signal changes sign from positive to negative or changes from negative to positive; orwherein the step of gating collection of data points from a localization signal comprises:integrating the acquired localization signal twice to calculate a displacement signal; andcollecting data points from the localization signal only when the calculated displacement signal is at a maximum value or when the calculated displacement signal is at a minimum value.
- A localization system for measuring a position of a measurement electrode within a non-ionizing localization field, comprising:a plurality of localization field generators capable of generating a non-ionizing localization field;at least one localization field reference device;at least one accelerometer capable of measuring acceleration in at least one direction coupled to at least one of the plurality of localization field generators or the at least one localization field reference device; anda signal processor that acquires at least one acceleration signal from the at least one accelerometer and that uses the acquired at least one acceleration signal to filter a respiration noise component from a localization signal of at least one position measurement sensor within the non-ionizing localization field.
- The system according to claim 13, wherein the signal processor integrates the at least one acceleration signal twice to compute at least one displacement signal and uses the computed at least one displacement signal to filter the respiration noise component from the localization signal.
- The system according to claim 13 or 14, wherein the signal processor processes the localization signal using a high-pass filter having a cutoff frequency determined using the acquired at least one acceleration signal, or
wherein the signal processor processes the localization signal by gating collection of data points from the localization signal using the acquired at least one acceleration signal.
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US12/835,518 US8603004B2 (en) | 2010-07-13 | 2010-07-13 | Methods and systems for filtering respiration noise from localization data |
Publications (3)
Publication Number | Publication Date |
---|---|
EP2407118A2 true EP2407118A2 (en) | 2012-01-18 |
EP2407118A3 EP2407118A3 (en) | 2014-03-19 |
EP2407118B1 EP2407118B1 (en) | 2019-03-06 |
Family
ID=44582194
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
EP11173450.5A Active EP2407118B1 (en) | 2010-07-13 | 2011-07-11 | Systems for filtering respiration noise from localization data |
Country Status (2)
Country | Link |
---|---|
US (1) | US8603004B2 (en) |
EP (1) | EP2407118B1 (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110151156A (en) * | 2019-04-07 | 2019-08-23 | 西安电子科技大学 | A kind of Fetal Heart Rate extracting method and system based on Weighted Average Algorithm automatic in window |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9414770B2 (en) * | 2010-12-29 | 2016-08-16 | Biosense Webster (Israel) Ltd. | Respiratory effect reduction in catheter position sensing |
CN105007810B (en) * | 2013-03-15 | 2018-05-22 | 莱彻韦斯科勒公司 | For measuring the method for tissue impedance and its device |
US9237936B2 (en) | 2013-07-12 | 2016-01-19 | Pacesetter, Inc. | System and method for integrating candidate implant location test results with real-time tissue images for use with implantable device leads |
US10758133B2 (en) * | 2014-08-07 | 2020-09-01 | Apple Inc. | Motion artifact removal by time domain projection |
CN105869598B (en) * | 2016-06-07 | 2019-06-11 | 武汉华星光电技术有限公司 | The driving method and liquid crystal display of liquid crystal display |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5697377A (en) | 1995-11-22 | 1997-12-16 | Medtronic, Inc. | Catheter mapping system and method |
US6640119B1 (en) | 1992-09-23 | 2003-10-28 | Endocardial Solutions, Inc. | Method for orienting an electrode array |
US6939309B1 (en) | 1993-09-23 | 2005-09-06 | Endocardial Solutions, Inc. | Electrophysiology mapping system |
US6978168B2 (en) | 1992-09-23 | 2005-12-20 | Endocardial Solutions, Inc. | Software for mapping potential distribution of a heart chamber |
US7263397B2 (en) | 1998-06-30 | 2007-08-28 | St. Jude Medical, Atrial Fibrillation Division, Inc. | Method and apparatus for catheter navigation and location and mapping in the heart |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7697972B2 (en) * | 2002-11-19 | 2010-04-13 | Medtronic Navigation, Inc. | Navigation system for cardiac therapies |
JP4610560B2 (en) * | 2003-05-21 | 2011-01-12 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | Apparatus and method for navigating a catheter |
US7885707B2 (en) | 2005-09-15 | 2011-02-08 | St. Jude Medical, Atrial Fibrillation Division, Inc. | Method of scaling navigation signals to account for impedance drift in tissue |
US20070270686A1 (en) * | 2006-05-03 | 2007-11-22 | Ritter Rogers C | Apparatus and methods for using inertial sensing to navigate a medical device |
KR100809041B1 (en) * | 2006-06-20 | 2008-03-03 | 삼성전자주식회사 | Apparatus for sensing sleeping comdition and method for operating the apparatus |
US9220439B2 (en) * | 2006-12-29 | 2015-12-29 | St. Jude Medical, Atrial Fibrillation Division, Inc. | Navigational reference dislodgement detection method and system |
US8527049B2 (en) | 2008-12-11 | 2013-09-03 | Pacesetter, Inc. | Cardiac resynchronization therapy optimization using vector measurements obtained from realtime electrode position tracking |
JP5795576B2 (en) * | 2009-06-12 | 2015-10-14 | バード・アクセス・システムズ,インコーポレーテッド | Method of operating a computer-based medical device that uses an electrocardiogram (ECG) signal to position an intravascular device in or near the heart |
US8412327B2 (en) | 2009-11-18 | 2013-04-02 | Pacesetter, Inc. | Cardiac resynchronization therapy optimization using vector measurements obtained from realtime electrode position tracking |
-
2010
- 2010-07-13 US US12/835,518 patent/US8603004B2/en active Active
-
2011
- 2011-07-11 EP EP11173450.5A patent/EP2407118B1/en active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6640119B1 (en) | 1992-09-23 | 2003-10-28 | Endocardial Solutions, Inc. | Method for orienting an electrode array |
US6728562B1 (en) | 1992-09-23 | 2004-04-27 | Endocardial Solutions, Inc. | Method for creating a virtual electrogram |
US6978168B2 (en) | 1992-09-23 | 2005-12-20 | Endocardial Solutions, Inc. | Software for mapping potential distribution of a heart chamber |
US6990370B1 (en) | 1992-09-23 | 2006-01-24 | Endocardial Solutions, Inc. | Method for mapping heart electrophysiology |
US6939309B1 (en) | 1993-09-23 | 2005-09-06 | Endocardial Solutions, Inc. | Electrophysiology mapping system |
US6947785B1 (en) | 1993-09-23 | 2005-09-20 | Endocardial Solutions, Inc. | Interface system for endocardial mapping catheter |
US5697377A (en) | 1995-11-22 | 1997-12-16 | Medtronic, Inc. | Catheter mapping system and method |
US5983126A (en) | 1995-11-22 | 1999-11-09 | Medtronic, Inc. | Catheter location system and method |
US7263397B2 (en) | 1998-06-30 | 2007-08-28 | St. Jude Medical, Atrial Fibrillation Division, Inc. | Method and apparatus for catheter navigation and location and mapping in the heart |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110151156A (en) * | 2019-04-07 | 2019-08-23 | 西安电子科技大学 | A kind of Fetal Heart Rate extracting method and system based on Weighted Average Algorithm automatic in window |
Also Published As
Publication number | Publication date |
---|---|
US8603004B2 (en) | 2013-12-10 |
EP2407118B1 (en) | 2019-03-06 |
US20120016253A1 (en) | 2012-01-19 |
EP2407118A3 (en) | 2014-03-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP2120700B1 (en) | System and method for correction of inhomogeneous fields | |
EP2407118B1 (en) | Systems for filtering respiration noise from localization data | |
JP7093776B2 (en) | Systems and Methods for Generating Premature Ventricular Contraction Electrophysiological Maps | |
US11380029B2 (en) | System and method for mapping local activation times | |
CN108697334B (en) | Method and system for mapping cardiac activity | |
CN110381813B (en) | System and method for distinguishing adipose tissue and scar tissue during electrophysiology mapping | |
JP6914383B2 (en) | Methods and systems for mapping cardiac repolarization | |
US20220287614A1 (en) | System and method for electrophysiological mapping | |
US10349856B2 (en) | Methods and systems for mapping cardiac restitution | |
CN103220970A (en) | Navigation reference dislodgement detection method and system | |
US11103177B2 (en) | System and method for mapping cardiac activity | |
JP2022532843A (en) | Systems and methods for cardiac mapping | |
JP7410970B2 (en) | Systems and methods for cardiac mapping | |
JP2023518239A (en) | Systems, methods, and apparatus for mapping local activation times | |
EP3609400A1 (en) | System and method for electrophysiological mapping | |
CN111065318A (en) | System and method for classifying electrophysiological signals from multi-dimensional catheters | |
US20220020228A1 (en) | Systems and methods for modifying geometry surface models using electrophysiology measurements | |
CN117897097A (en) | Method and system for generating respiratory signals for use in electrophysiological procedures |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AK | Designated contracting states |
Kind code of ref document: A2 Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR |
|
AX | Request for extension of the european patent |
Extension state: BA ME |
|
PUAI | Public reference made under article 153(3) epc to a published international application that has entered the european phase |
Free format text: ORIGINAL CODE: 0009012 |
|
PUAL | Search report despatched |
Free format text: ORIGINAL CODE: 0009013 |
|
AK | Designated contracting states |
Kind code of ref document: A3 Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR |
|
AX | Request for extension of the european patent |
Extension state: BA ME |
|
RIC1 | Information provided on ipc code assigned before grant |
Ipc: A61B 19/00 20060101AFI20140213BHEP |
|
17P | Request for examination filed |
Effective date: 20140523 |
|
RBV | Designated contracting states (corrected) |
Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: EXAMINATION IS IN PROGRESS |
|
17Q | First examination report despatched |
Effective date: 20161114 |
|
REG | Reference to a national code |
Ref country code: DE Ref legal event code: R079 Ref document number: 602011056807 Country of ref document: DE Free format text: PREVIOUS MAIN CLASS: A61B0019000000 Ipc: A61B0034200000 |
|
GRAP | Despatch of communication of intention to grant a patent |
Free format text: ORIGINAL CODE: EPIDOSNIGR1 |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: GRANT OF PATENT IS INTENDED |
|
RIC1 | Information provided on ipc code assigned before grant |
Ipc: A61B 34/20 20160101AFI20180831BHEP |
|
INTG | Intention to grant announced |
Effective date: 20180927 |
|
GRAS | Grant fee paid |
Free format text: ORIGINAL CODE: EPIDOSNIGR3 |
|
GRAA | (expected) grant |
Free format text: ORIGINAL CODE: 0009210 |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: THE PATENT HAS BEEN GRANTED |
|
AK | Designated contracting states |
Kind code of ref document: B1 Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR |
|
REG | Reference to a national code |
Ref country code: GB Ref legal event code: FG4D |
|
REG | Reference to a national code |
Ref country code: CH Ref legal event code: EP Ref country code: AT Ref legal event code: REF Ref document number: 1103511 Country of ref document: AT Kind code of ref document: T Effective date: 20190315 |
|
REG | Reference to a national code |
Ref country code: DE Ref legal event code: R096 Ref document number: 602011056807 Country of ref document: DE |
|
REG | Reference to a national code |
Ref country code: IE Ref legal event code: FG4D |
|
REG | Reference to a national code |
Ref country code: NL Ref legal event code: MP Effective date: 20190306 |
|
REG | Reference to a national code |
Ref country code: LT Ref legal event code: MG4D |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: LT Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20190306 Ref country code: SE Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20190306 Ref country code: NO Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20190606 Ref country code: FI Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20190306 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: NL Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20190306 Ref country code: HR Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20190306 Ref country code: GR Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20190607 Ref country code: LV Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20190306 Ref country code: BG Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20190606 Ref country code: RS Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20190306 |
|
REG | Reference to a national code |
Ref country code: AT Ref legal event code: MK05 Ref document number: 1103511 Country of ref document: AT Kind code of ref document: T Effective date: 20190306 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: ES Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20190306 Ref country code: PT Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20190706 Ref country code: AL Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20190306 Ref country code: EE Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20190306 Ref country code: SK Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20190306 Ref country code: RO Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20190306 Ref country code: CZ Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20190306 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: PL Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20190306 Ref country code: SM Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20190306 |
|
REG | Reference to a national code |
Ref country code: DE Ref legal event code: R097 Ref document number: 602011056807 Country of ref document: DE |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: AT Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20190306 Ref country code: IS Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20190706 |
|
PLBE | No opposition filed within time limit |
Free format text: ORIGINAL CODE: 0009261 |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: NO OPPOSITION FILED WITHIN TIME LIMIT |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: DK Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20190306 |
|
26N | No opposition filed |
Effective date: 20191209 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: MC Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20190306 Ref country code: SI Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20190306 |
|
REG | Reference to a national code |
Ref country code: CH Ref legal event code: PL |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: TR Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20190306 |
|
REG | Reference to a national code |
Ref country code: BE Ref legal event code: MM Effective date: 20190731 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: BE Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20190731 Ref country code: LI Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20190731 Ref country code: LU Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20190711 Ref country code: CH Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20190731 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: IE Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20190711 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: CY Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20190306 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: HU Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT; INVALID AB INITIO Effective date: 20110711 Ref country code: MT Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20190306 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: MK Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20190306 |
|
P01 | Opt-out of the competence of the unified patent court (upc) registered |
Effective date: 20230616 |
|
PGFP | Annual fee paid to national office [announced via postgrant information from national office to epo] |
Ref country code: FR Payment date: 20230616 Year of fee payment: 13 |
|
REG | Reference to a national code |
Ref country code: DE Ref legal event code: R082 Ref document number: 602011056807 Country of ref document: DE Representative=s name: ALPSPITZ IP ALLGAYER UND PARTNER PATENTANWAELT, DE |
|
PGFP | Annual fee paid to national office [announced via postgrant information from national office to epo] |
Ref country code: IT Payment date: 20230712 Year of fee payment: 13 Ref country code: GB Payment date: 20230614 Year of fee payment: 13 |
|
PGFP | Annual fee paid to national office [announced via postgrant information from national office to epo] |
Ref country code: DE Payment date: 20230614 Year of fee payment: 13 |