US10349852B2 - Analysing physiological electrograms - Google Patents
Analysing physiological electrograms Download PDFInfo
- Publication number
- US10349852B2 US10349852B2 US15/503,674 US201515503674A US10349852B2 US 10349852 B2 US10349852 B2 US 10349852B2 US 201515503674 A US201515503674 A US 201515503674A US 10349852 B2 US10349852 B2 US 10349852B2
- Authority
- US
- United States
- Prior art keywords
- signal
- noise
- electrogram
- model
- physiological
- 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.)
- Active
Links
Images
Classifications
-
- A61B5/04012—
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
- A61B5/346—Analysis of electrocardiograms
-
- A61B5/0402—
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
-
- 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
-
- 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/7217—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise originating from a therapeutic or surgical apparatus, e.g. from a pacemaker
-
- 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/7235—Details of waveform analysis
- A61B5/7246—Details of waveform analysis using correlation, e.g. template matching or determination of similarity
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N1/00—Electrotherapy; Circuits therefor
- A61N1/18—Applying electric currents by contact electrodes
- A61N1/32—Applying electric currents by contact electrodes alternating or intermittent currents
- A61N1/36—Applying electric currents by contact electrodes alternating or intermittent currents for stimulation
- A61N1/362—Heart stimulators
- A61N1/37—Monitoring; Protecting
- A61N1/3702—Physiological parameters
-
- A61B5/04525—
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
- A61B5/346—Analysis of electrocardiograms
- A61B5/349—Detecting specific parameters of the electrocardiograph cycle
- A61B5/35—Detecting specific parameters of the electrocardiograph cycle by template matching
-
- 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/7271—Specific aspects of physiological measurement analysis
- A61B5/7275—Determining trends in physiological measurement data; Predicting development of a medical condition based on physiological measurements, e.g. determining a risk factor
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N1/00—Electrotherapy; Circuits therefor
- A61N1/18—Applying electric currents by contact electrodes
- A61N1/32—Applying electric currents by contact electrodes alternating or intermittent currents
- A61N1/36—Applying electric currents by contact electrodes alternating or intermittent currents for stimulation
- A61N1/362—Heart stimulators
- A61N1/3627—Heart stimulators for treating a mechanical deficiency of the heart, e.g. congestive heart failure or cardiomyopathy
Definitions
- the present invention relates to the analysis of physiological electrograms, in particularly but not exclusively for identifying pathological cardiac conditions.
- the diagnostic change in electrograms consists of small deflections in the recorded electrogram following early stimulation of the heart.
- the heart is stimulated with apparatus that generates a stimulation sequence at one site in the heart and records electrograms from other sites within the heart.
- the pacing sequence comprises of a number of stimuli at a constant rate, known as S 1 stimuli. After a pre-set number of S 1 stimuli an early stimulus is introduced known as the S 2 stimulus or ‘extra-stimulus’. The sequence is repeated. Typically the interval between S 1 and S 2 stimuli is reduced on each occasion until the interval is so short that the heart is no longer able to respond to the S 2 stimulus.
- the interval between the S 2 stimulus and the following S 1 stimulus is the same as the S 1 -S 1 interval.
- the predictive method depends on demonstrating that the electrogram following an extra-stimulus becomes prolonged and contains more peaks. This effect in patients that are at high risk of sudden death becomes more pronounced as the interval between the S 1 and S 2 stimuli is reduced.
- Each individual potential of the electrogram following an extrastimulus is identified together with its delay after the extra stimulus. These data can subsequently be analyzed to predict the risk of sudden cardiac death.
- a long standing problem with this method has been the reliable detection of small individual potentials within the response to an extra stimulus. This stems from the presences of noise in the recorded electrical signals that may be created by other electrical equipment within a typical catheter laboratory. The electrical noise may vary widely between different laboratories. The problem is that of reliably distinguishing between potentials in the electrogram that are of physiological origin as opposed to spurious potentials caused by electrical interference.
- GB2439562 describes a method of processing date from electrograms to reduce noise.
- the method comprises correlating an electrogram signal with several templates to produce a correlator output associated with each template.
- the electrogram signal may be passed through a high pass filter beforehand.
- the correlator output from trace 1 is compared with the traces produced from the other templates. The selected trace that is considered most similar is used.
- the noise can be reduced further by identifying the peak-to-peak amplitude of the correlated output within a period of the electrogram when no physiological signal is presumed to occur. This is used to create a threshold in which any peak having an amplitude below this threshold is considered to be noise.
- signals of physiological origin may have amplitudes which are close to the threshold, as a consequence is if the threshold is set too low, the physiological derived peaks will be detected but many other peaks will also be detected due to noise. Conversely, if the amplitude threshold is too high, physiologically important features of the signal may not be detected.
- An object of the invention is to overcome or at least ameliorate the above problems.
- a method of analyzing an electrogram for example a cardiac electrogram, to distinguish a physiological signal from noise; the method comprising:
- the model of noise is derived from multiple portions of the electrogram in which a physiological signal is presumed to be absent. This provides a more accurate means of representing the noise thereby enabling improved detection of physiological signals from the noise.
- model of noise is derived by cross-correlation of the first portion of the electrogram in which a physiological signal is presumed to be absent with multiple templates that represent features of the presumed physiological signal to produce a number of template correlated signals.
- the number and form of the templates will depend on the signal in question and can be determined by experiment.
- a set of time dilated templates are used that correspond to different local conduction velocities in the region of the recording electrodes.
- a co-variance matrix is derived from the template correlated signals for each portion of the electrogram in which a physiological signal is presumed to be absent.
- the eigenvectors thence derived are wholly real and orthogonal (and thus independent), and the eigenvalues are real. In this way it is possible to differentiate a signal from noise because the signal will have significant components in eigenvectors where the noise is very small or non-existent.
- the method preferably comprises deriving a mean co-variance matrix from the co-variances matrices derived for each portion of the cardiac electrogram in which a physiological signal is presumed to be absent.
- the mean covariance matrix provides a better estimate of the behaviour of the noise throughout the entire recording.
- model of noise is expressed by deriving eigenvectors and eigenvalues from the mean co-variance matrix.
- the second portion of the electrogram is correlated with multiple templates that represent features of the presumed physiological signal to produce a set of template correlated signals, and additionally favorable that a vector is derived from a first time sample of each template correlated signal of the set, and further vectors from further time samples of each template correlated signal of the set.
- the vectors are represented as points in the model of noise by projecting each vector onto each eigenvector thereby representing the original signal as a trajectory in the model of noise, and so allows comparison of the signal and the noise.
- the physiological signal is preferably identified by determining points that lie outside the limits of the model of noise thereby discriminating the signal from the noise.
- the invention can also be expressed in terms of apparatus and thus according a further aspect of the invention there is provided apparatus for analyzing an electrogram to distinguish a physiological signal from noise; the apparatus comprising:
- FIG. 1 is a schematic representation of apparatus for pacing a heart, recording an electrogram and subsequent analysis
- FIG. 2 is schematic representation of a paced cardiac electrogram sequence
- FIG. 3 is schematic representation similar to FIG. 1 illustrating the regions of the signal used to evaluate noise
- FIG. 4 is a schematic illustrating the principle of the analysis technique
- FIG. 5 is a schematic of time domain representation of the templates
- FIG. 6 illustrates a correlation matrix of the templates showing correlation therebetween
- FIG. 7 is a flow diagram illustrating the steps for creation of a noise model
- FIG. 8 is a graphical illustration of the noise model for three templates
- FIG. 9 is a flow diagram illustrated the steps for identification of potentials the noise model
- FIG. 10 is a graphical illustration of the noise model for three templates showing trajectory of one noise record within the noise model illustrating that it remains within the limits of the model;
- FIG. 11 a graphical illustration of the noise model for three templates showing a trajectory of a signal that is derived from a physiological signal
- FIG. 12 is a graphically illustrates a normalized noise model for three templates
- FIG. 13 illustrates a simulated electrogram to which noise and small potentials have been added
- FIG. 14 illustrates the simulated electrogram of FIG. 13 with the time domain output signal of the model superimposed.
- the terms “a” or “an”, as used herein, are defined as one or more than one.
- the term “plurality”, as used herein, is defined as two or more than two.
- the term “another”, as used herein, is defined as at least a second or more.
- the terms “including” and “having,” as used herein, are defined as comprising i.e., open language.
- the term “coupled,” as used herein, is defined as “connected,” although not necessarily directly, and not necessarily mechanically. “Communicatively coupled” refers to coupling of components such that these components are able to communicate with one another through, for example, wired, wireless or other communications media.
- the term “communicatively coupled” or “communicatively coupling” includes, but is not limited to, communicating electronic control signals by which one element may direct or control another.
- the term “configured to” describes hardware, software or a combination of hardware and software that is adapted to, set up, arranged, commanded, altered, modified, built, composed, constructed, designed, or that has any combination of these characteristics to carry out a given function.
- the term “adapted to” describes hardware, software or a combination of hardware and software that is capable of, able to accommodate, to make, or that is suitable to carry out a given function.
- the apparatus comprises amplifiers 1 that are connected to recording electrodes 2 within the heart 3 and electronics and associated software 4 .
- a pacing and recording program 5 issues signals to a pacing signal generator 6 that is switched onto one selected electrode 2 by multiplexer 7 to stimulate the heart 3 .
- the signals sensed by the other electrodes 2 are amplified and digitized by an ADC 8 and stored in memory 9 . Subsequently the data is retrieved and analyzed by analysis program 10 .
- the use of the multiplexer 7 allows the heart to be stimulated at different sites 3 .
- FIG. 2 illustrates a pacing sequence applied at one electrode 2 , showing the constant rate stimuli S 1 and extra-stimuli S 2 .
- the intervals between S 1 -S 1 stimuli and the S 2 -S 1 stimuli remain constant, in this example with an interval of 500 ms.
- the interval between the S 1 and S 2 stimuli varies, and typically reduces by one 1 ms on each occasion.
- FIG. 3 illustrates the regions 11 preceding S 1 stimuli that are used for evaluating noise on the premise that no physiological signals occur in these regions 11 . Also shown are regions 12 that are analyzed to identify physiological signals resulting from S 2 stimuli.
- the noise regions 11 are extracted from the recorded signal and correlated with templates 13 , generated by the analysis program 10 , to create a model 14 with independent eigenvectors.
- the templates 13 are representations of potentials of varying widths that are likely to be the result of a physiological event.
- the portions 12 of the recorded signal are correlated to templates 13 and projected into model 14 that provides an output 15 that is indicative of whether the sample being analyzed is signal or noise.
- FIG. 5 illustrates the fifteen templates 13 that are correlated with the signals 11 12 .
- Template number one 13 A has the longest time duration while template number fifteen 13 B has the shortest.
- the intermediate template numbers 13 C shorten progressively.
- templates are always used in their frequency domain representations, i.e. their discrete Fourier transform, for computational efficiency.
- the correlation between templates 13 is shown as a correlation matrix, see FIG. 6 , in which each element is the correlation between template n and template m where n and m are the template numbers. This shows that the templates 13 are not independent since if they were all non-diagonal elements would be equal to zero.
- the templates 13 are computed in the frequency domain 20 .
- a record of each noise region 11 is transformed into the frequency domain by Fast Fourier transform 21 and correlated 22 with each template by multiplication. The result is summed to form the mean cross-correlation spectrum 23 .
- a mean noise co-variance matrix is derived 24 from the cross-correlation spectrum. The co-variance matrix is decomposed 25 into its eigenvectors and eigenvalues, thus forming a noise model 14 .
- FIG. 8 illustrates a noise model for three templates only (because of the difficulties of showing high dimensional models).
- the three eigenvectors ⁇ 1 , ⁇ 2 , ⁇ 3 are at rights angles, i.e. orthogonal.
- the limits of the noise model is shown by the shaded area 40 as defined by the eigenvalues.
- every signal identification region 12 is processed individually.
- the individual signal portion 26 is transformed into the frequency domain 27 and correlated 28 with templates 13 and the result expressed 29 in the time domain.
- the resultant signal template correlation records are projected 30 into the eigenvector space of the noise model 14 as a trajectory.
- FIG. 10 is an illustration of the trajectory 41 of a signal within region 11 illustrating that it remains within the limits 40 of the model 14 as defined by the eigenvectors/eigenvalues.
- FIG. 11 illustrates the trajectory 42 of a signal within the region 12 showing that it exceeds the limits 40 of the noise model 14 and thus is likely to be attributed to a physiological origin.
- a convenient method of determining whether the signal exceeds the noise is to reduce the trajectory to a single time domain signal.
- the resultant trajectory is normalized by division of each eigenvector by its eigenvalue so that the noise model becomes a spheroid 43 as illustrated in FIG. 12 .
- the norm 45 of the trajectory vector in this space 44 is computed to give the one-dimensional time domain signal.
- FIG. 13 shows a portion of a simulated electrogram 50 , corresponding to a single S 1 -S 2 interval. Noise has been added to the electrogram together with small potentials expected from a physiological response. The small potentials have the same peak-to-peak amplitude of the noise; the positions of the small potentials are indicated by arrows 51 .
- FIG. 14 shows the electrogram of FIG. 13 with a time domain signal output 52 derived using the method described above superimposed.
- the time domain signal output shows clear peaks that exceed the noise threshold 53 corresponding to the limit of the noise model, at the positions that the small potentials were inserted into the signal with a significant increase in signal to noise ratio.
- Computer program in the present context means any expression, in any language, code or notation, of a set of instructions intended to cause a system having an information processing capability to perform a particular function either directly or after either or both of the following a conversion to another language, code or, notation; and b reproduction in a different material form.
- Each computer system may include, inter alia, one or more computers and at least a computer readable medium allowing a computer to read data, instructions, messages or message packets, and other computer readable information from the computer readable medium.
- the computer readable medium may include computer readable storage medium embodying non-volatile memory, such as read-only memory ROM, flash memory, disk drive memory, CD-ROM, and other permanent storage. Additionally, a computer medium may include volatile storage such as RAM, buffers, cache memory, and network circuits.
- the computer readable medium may comprise computer readable information in a transitory state medium such as a network link and/or a network interface, including a wired network or a wireless network, that allow a computer to read such computer readable information.
- a transitory state medium such as a network link and/or a network interface, including a wired network or a wireless network, that allow a computer to read such computer readable information.
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Heart & Thoracic Surgery (AREA)
- Biomedical Technology (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- Surgery (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Pathology (AREA)
- Physiology (AREA)
- Psychiatry (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Cardiology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Hospice & Palliative Care (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
Previous research has shown that the risk of sudden death due to cardiac arrhythmias can be predicted by observing the shape of recorded endocardial electrograms in response to pacing, and in particularly detecting certain small deflections in the recorded electrogram following early stimulation of the heart. A long standing problem has been the reliable detection of these small individual potentials because of the presence of noise in the recorded electrical signals created by other electrical equipment within a typical catheter laboratory. The solution described involves deriving a model of noise from a first portion of the electrogram in which a physiological signal is presumed to be absent, and transforming a second portion of the electrogram, presumed to contain a physiological signal, into the model of noise. The physiological signal can then be identified by identifying portions of signal within the second portion of the electrogram that do not conform to the model of noise.
Description
This application is based upon and claims priority to PCT Patent Application No. PCT/GB2015/052190 filed on 29 Jul. 2015, which is based upon and claims priority to GB Patent Application No. 1414330.9 filed on 13 Aug. 2014, the collective disclosure of which being hereby incorporated by reference in its entirety.
The present invention relates to the analysis of physiological electrograms, in particularly but not exclusively for identifying pathological cardiac conditions.
Previous research has shown that the risk of sudden death due to cardiac arrhythmias can be predicted by observing the shape of recorded endocardial electrograms in response to pacing.
The diagnostic change in electrograms consists of small deflections in the recorded electrogram following early stimulation of the heart. The heart is stimulated with apparatus that generates a stimulation sequence at one site in the heart and records electrograms from other sites within the heart.
The pacing sequence comprises of a number of stimuli at a constant rate, known as S1 stimuli. After a pre-set number of S1 stimuli an early stimulus is introduced known as the S2 stimulus or ‘extra-stimulus’. The sequence is repeated. Typically the interval between S1 and S2 stimuli is reduced on each occasion until the interval is so short that the heart is no longer able to respond to the S2 stimulus.
The interval between the S2 stimulus and the following S1 stimulus is the same as the S1-S1 interval.
The predictive method depends on demonstrating that the electrogram following an extra-stimulus becomes prolonged and contains more peaks. This effect in patients that are at high risk of sudden death becomes more pronounced as the interval between the S1 and S2 stimuli is reduced.
Each individual potential of the electrogram following an extrastimulus is identified together with its delay after the extra stimulus. These data can subsequently be analyzed to predict the risk of sudden cardiac death.
A long standing problem with this method has been the reliable detection of small individual potentials within the response to an extra stimulus. This stems from the presences of noise in the recorded electrical signals that may be created by other electrical equipment within a typical catheter laboratory. The electrical noise may vary widely between different laboratories. The problem is that of reliably distinguishing between potentials in the electrogram that are of physiological origin as opposed to spurious potentials caused by electrical interference.
GB2439562 describes a method of processing date from electrograms to reduce noise. The method comprises correlating an electrogram signal with several templates to produce a correlator output associated with each template. The electrogram signal may be passed through a high pass filter beforehand.
The correlator output from trace 1 is compared with the traces produced from the other templates. The selected trace that is considered most similar is used.
A fundamental problem is that any series of templates that purport to represent a physiological signal will be correlated, and therefore the results of each correlated trace will not be independent of each other. This creates considerable difficulties in to how to combine the various correlator outputs to give optimal signal detection and avoid spurious over detection and under detection of physiological potentials within the signal.
The noise can be reduced further by identifying the peak-to-peak amplitude of the correlated output within a period of the electrogram when no physiological signal is presumed to occur. This is used to create a threshold in which any peak having an amplitude below this threshold is considered to be noise. However, signals of physiological origin may have amplitudes which are close to the threshold, as a consequence is if the threshold is set too low, the physiological derived peaks will be detected but many other peaks will also be detected due to noise. Conversely, if the amplitude threshold is too high, physiologically important features of the signal may not be detected.
An object of the invention is to overcome or at least ameliorate the above problems.
According to first aspect of the invention there is provided a method of analyzing an electrogram, for example a cardiac electrogram, to distinguish a physiological signal from noise; the method comprising:
deriving a model of noise from a first portion of the electrogram in which a physiological signal is presumed to be absent; and
transforming a second portion of the electrogram, presumed to contain a physiological signal, into the model of noise, and wherein the physiological signal is identified by identifying portions of the signal that do not conform to the model of noise.
This provides an improvement over using the amplitude method because detected potentials are more likely to be genuinely physiological in origin and so subsequent analysis is greatly simplified.
In a preferred embodiment, the model of noise is derived from multiple portions of the electrogram in which a physiological signal is presumed to be absent. This provides a more accurate means of representing the noise thereby enabling improved detection of physiological signals from the noise.
It is favorable that the model of noise is derived by cross-correlation of the first portion of the electrogram in which a physiological signal is presumed to be absent with multiple templates that represent features of the presumed physiological signal to produce a number of template correlated signals.
The number and form of the templates will depend on the signal in question and can be determined by experiment. In the context of cardiac electrograms, a set of time dilated templates are used that correspond to different local conduction velocities in the region of the recording electrodes.
It is preferred that a co-variance matrix is derived from the template correlated signals for each portion of the electrogram in which a physiological signal is presumed to be absent. By deriving covariance matrices that are inherently symmetric, the eigenvectors thence derived are wholly real and orthogonal (and thus independent), and the eigenvalues are real. In this way it is possible to differentiate a signal from noise because the signal will have significant components in eigenvectors where the noise is very small or non-existent.
The method preferably comprises deriving a mean co-variance matrix from the co-variances matrices derived for each portion of the cardiac electrogram in which a physiological signal is presumed to be absent. The mean covariance matrix provides a better estimate of the behaviour of the noise throughout the entire recording.
It is preferred that the model of noise is expressed by deriving eigenvectors and eigenvalues from the mean co-variance matrix.
It is preferred that the second portion of the electrogram is correlated with multiple templates that represent features of the presumed physiological signal to produce a set of template correlated signals, and additionally favorable that a vector is derived from a first time sample of each template correlated signal of the set, and further vectors from further time samples of each template correlated signal of the set.
It is preferred that the vectors are represented as points in the model of noise by projecting each vector onto each eigenvector thereby representing the original signal as a trajectory in the model of noise, and so allows comparison of the signal and the noise.
The physiological signal is preferably identified by determining points that lie outside the limits of the model of noise thereby discriminating the signal from the noise.
The invention can also be expressed in terms of apparatus and thus according a further aspect of the invention there is provided apparatus for analyzing an electrogram to distinguish a physiological signal from noise; the apparatus comprising:
means for deriving a model of noise from a first portion of the cardiac electrogram in which a physiological signal is presumed to be absent; and
means for transforming a second portion of the electrogram presumed to contain a physiological signal into the model of noise, and where the physiological signal is identified by identifying portions of the signal that do not conform to the model of noise.
The accompanying figures in which like reference numerals refer to identical or functionally similar elements throughout the separate views, and which together with the detailed description below are incorporated in and form part of the specification, serve to further illustrate various embodiments and to explain various principles and advantages all in accordance with the present disclosure, in which:
As required, detailed embodiments are disclosed herein; however, it is to be understood that the disclosed embodiments are merely examples and that the devices, systems and methods described herein can be embodied in various forms. Therefore, specific structural and functional details disclosed herein are not to be interpreted as limiting, but merely as a basis for the claims and as a representative basis for teaching one of ordinary skill in the art to variously employ the disclosed subject matter in virtually any appropriately detailed structure and function. Further, the terms and phrases used herein are not intended to be limiting, but rather, to provide an understandable description. Additionally, unless otherwise specifically expressed or clearly understood from the context of use, a term as used herein describes the singular and/or the plural of that term.
The terms “a” or “an”, as used herein, are defined as one or more than one. The term “plurality”, as used herein, is defined as two or more than two. The term “another”, as used herein, is defined as at least a second or more. The terms “including” and “having,” as used herein, are defined as comprising i.e., open language. The term “coupled,” as used herein, is defined as “connected,” although not necessarily directly, and not necessarily mechanically. “Communicatively coupled” refers to coupling of components such that these components are able to communicate with one another through, for example, wired, wireless or other communications media. The term “communicatively coupled” or “communicatively coupling” includes, but is not limited to, communicating electronic control signals by which one element may direct or control another. The term “configured to” describes hardware, software or a combination of hardware and software that is adapted to, set up, arranged, commanded, altered, modified, built, composed, constructed, designed, or that has any combination of these characteristics to carry out a given function. The term “adapted to” describes hardware, software or a combination of hardware and software that is capable of, able to accommodate, to make, or that is suitable to carry out a given function.
The apparatus comprises amplifiers 1 that are connected to recording electrodes 2 within the heart 3 and electronics and associated software 4. A pacing and recording program 5 issues signals to a pacing signal generator 6 that is switched onto one selected electrode 2 by multiplexer 7 to stimulate the heart 3. The signals sensed by the other electrodes 2 are amplified and digitized by an ADC 8 and stored in memory 9. Subsequently the data is retrieved and analyzed by analysis program 10.
All subsequent analysis is with sampled signals, care having been taken to conform to the Nyquist sampling theorem.
The use of the multiplexer 7 allows the heart to be stimulated at different sites 3.
The functions and arrangement described above can be derived, in conjunction with the teaching within this document by the person skilled in the art.
As illustrated in FIG. 4 , the noise regions 11 are extracted from the recorded signal and correlated with templates 13, generated by the analysis program 10, to create a model 14 with independent eigenvectors. The templates 13 are representations of potentials of varying widths that are likely to be the result of a physiological event.
Subsequently, the portions 12 of the recorded signal are correlated to templates 13 and projected into model 14 that provides an output 15 that is indicative of whether the sample being analyzed is signal or noise.
In practice the templates are always used in their frequency domain representations, i.e. their discrete Fourier transform, for computational efficiency.
The correlation between templates 13 is shown as a correlation matrix, see FIG. 6 , in which each element is the correlation between template n and template m where n and m are the template numbers. This shows that the templates 13 are not independent since if they were all non-diagonal elements would be equal to zero.
Referring to FIG. 7 , the templates 13 are computed in the frequency domain 20. A record of each noise region 11 is transformed into the frequency domain by Fast Fourier transform 21 and correlated 22 with each template by multiplication. The result is summed to form the mean cross-correlation spectrum 23. A mean noise co-variance matrix is derived 24 from the cross-correlation spectrum. The co-variance matrix is decomposed 25 into its eigenvectors and eigenvalues, thus forming a noise model 14.
Referring to FIG. 9 , every signal identification region 12 is processed individually. The individual signal portion 26 is transformed into the frequency domain 27 and correlated 28 with templates 13 and the result expressed 29 in the time domain. The resultant signal template correlation records are projected 30 into the eigenvector space of the noise model 14 as a trajectory.
A convenient method of determining whether the signal exceeds the noise is to reduce the trajectory to a single time domain signal. To achieve this, the resultant trajectory is normalized by division of each eigenvector by its eigenvalue so that the noise model becomes a spheroid 43 as illustrated in FIG. 12 . The norm 45 of the trajectory vector in this space 44 is computed to give the one-dimensional time domain signal.
Any peak in the time domain signal that is above of the noise is considered to be physiologically significant.
The present subject matter can be embedded in a computer program product, which comprises all the features enabling the implementation of the methods described herein, and which—when loaded in a computer system—is able to carry out these methods. Computer program in the present context means any expression, in any language, code or notation, of a set of instructions intended to cause a system having an information processing capability to perform a particular function either directly or after either or both of the following a conversion to another language, code or, notation; and b reproduction in a different material form.
Each computer system may include, inter alia, one or more computers and at least a computer readable medium allowing a computer to read data, instructions, messages or message packets, and other computer readable information from the computer readable medium. The computer readable medium may include computer readable storage medium embodying non-volatile memory, such as read-only memory ROM, flash memory, disk drive memory, CD-ROM, and other permanent storage. Additionally, a computer medium may include volatile storage such as RAM, buffers, cache memory, and network circuits. Furthermore, in certain embodiments of the computer readable medium, other than a computer readable storage medium as discussed above, the computer readable medium may comprise computer readable information in a transitory state medium such as a network link and/or a network interface, including a wired network or a wireless network, that allow a computer to read such computer readable information.
The Abstract is provided with the understanding that it is not intended be used to interpret or limit the scope or meaning of the claims. In addition, in the foregoing Detailed Description, various features are grouped together in a single embodiment for the purpose of streamlining the disclosure. This method of disclosure is not to be interpreted as reflecting an intention that the claimed embodiments require more features than are expressly recited in each claim. Rather, as the following claims reflect, inventive subject matter lies in less than all features of a single disclosed embodiment. Thus the following claims are hereby incorporated into the Detailed Description, with each claim standing on its own as a separately claimed subject matter.
The corresponding structures, materials, acts, and equivalents of all means or step plus function elements in the claims below are intended to include any structure, material, or act for performing the function in combination with other claimed elements as specifically claimed. The description herein has been presented for purposes of illustration and description, but is not intended to be exhaustive or limited to the examples in the form disclosed. Many modifications and variations will be apparent to those of ordinary skill in the art without departing from the scope of the examples presented or claimed. The disclosed embodiments were chosen and described in order to explain the principles of the embodiments and the practical application, and to enable others of ordinary skill in the art to understand the various embodiments with various modifications as are suited to the particular use contemplated. It is intended that the appended claims below cover any and all such applications, modifications, and variations within the scope of the embodiments.
Although specific embodiments of the subject matter have been disclosed, those having ordinary skill in the art will understand that changes can be made to the specific embodiments without departing from the scope of the disclosed subject matter. The scope of the disclosure is not to be restricted, therefore, to the specific embodiments, and it is intended that the appended claims cover any and all such applications, modifications, and embodiments within the scope of the present disclosure.
Claims (20)
1. A method of analysing an electrogram to distinguish a physiological signal from noise; the method comprising:
receiving a electrogram signal captured from a patient's heart, the electrogram signal including noise signals and physiological signals from the patient's heart, wherein:
a first portion of the electrogram signal includes a noise signal absent a physiological signal, and
a second portion of the electrogram signal includes a noise signal combined with a physiological signal;
deriving a model of noise from the first portion of the electrogram signal by cross-correlation of the first portion of the electrogram signal with multiple templates that represent features of the physiological signal to produce a number of template correlated signals, the model of noise having limits within which a captured electrogram signal is defined as a noise signal and outside of which a captured electrogram signal is defined as a physiological signal;
transforming, based on the model of noise, the second portion of the electrogram signal to a trajectory vector signal; and
identifying the physiological signal in the second portion of the electrogram signal by identifying portions of the second portion of the electrogram signal that were transformed to portions of the trajectory vector signal that do not conform to the model of noise.
2. A method according to claim 1 wherein the model of noise is derived from multiple first portions of the electrogram signal that include a noise signal absent a physiological signal.
3. A method according to claim 1 comprising deriving a co-variance matrix from the template correlated signals for a plurality of first portions of the electrogram signal that each includes a noise signal absent a physiological signal.
4. A method according to claim 3 comprising deriving a mean co-variance matrix from the co-variance matrices derived for each of a plurality of first portions of the electrogram signal that each includes a noise signal absent a physiological signal.
5. A method according to claim 4 wherein the model of noise is expressed by deriving Eigenvectors and Eigen values from the mean co-variance matrix.
6. A method according to claim 1 comprising correlating the second portion of the electrogram signal with multiple templates that represent features of the physiological signal to produce a set of template correlated signals.
7. A method according to claim 1 comprising deriving a vector from a first time sample of each template correlated signal of the set; and further vectors from further time samples of each template correlated signal of the set.
8. A method according to claim 7 comprising representing the vectors as points in the model of noise by projecting each vector onto each eigenvector thereby representing the original signal as a trajectory in the model of noise.
9. A method according to claim 8 wherein a physiological signal is identified by determining points that lie outside the limits of the model of noise.
10. A method of analyzing a cardiac electrogram according to claim 1 .
11. A method according to claim 1 comprising reducing the trajectory vector signal to a time domain signal, and outputting the time domain signal.
12. Apparatus for analysing a cardiac electrogram to distinguish a physiological signal from noise; the apparatus comprising:
a signal generator to generate a pacing signal;
an input electrode for applying the pacing signal to a patient's heart;
a receiving electrode to receive a cardiac electrogram signal from the patient's heart, the cardiac electrogram signal received based on application of the pacing signal to the patient's heart, the cardiac electrogram signal including noise signals and physiological signals from the patient's heart, wherein:
a first portion of the cardiac electrogram signal includes a noise signal absent a physiological signal, and
a second portion of the cardiac electrogram signal includes a noise signal combined with a physiological signal;
a memory store in which the received cardiac electrogram signal is stored; and
a processor communicatively coupled to the memory store, the processor being responsive to executing computer instructions, to perform operations comprising:
deriving a model of noise from the first portion of the cardiac electrogram signal by cross-correlation of the first portion of the electrogram signal with multiple templates that represent features of the physiological signal to produce a number of template correlated signals, the model of noise having limits within which a captured electrogram signal is defined as a noise signal and outside of which a captured electrogram signal is defined as a physiological signal;
transforming, based on the model of noise, the second portion of the electrogram signal to a trajectory vector signal; and
identifying the physiological signal in the second portion of the cardiac electrogram signal by identifying portions of the second portion of the cardiac electrogram signal that were transformed to portions of the trajectory vector signal that do not conform to the model of noise.
13. Apparatus according to claim 12 wherein the processor is responsive to executing computer instructions to derive the model of noise from multiple first portions of the cardiac electrogram signal that include a noise signal absent a physiological signal.
14. Apparatus according to claim 12 wherein the processor is responsive to executing computer instructions to derive a co-variance matrix from the template correlated signals for a plurality of first portions of the cardiac electrogram signal that each includes a noise signal absent a physiological signal.
15. Apparatus according to claim 14 wherein the processor is responsive to executing computer instructions to derive a mean co-variance matrix from the co-variance matrices derived for each of a plurality of first portions of the cardiac electrogram signal that each includes a noise signal absent a physiological signal.
16. Apparatus according to claim 15 wherein the processor is responsive to executing computer instructions to express the model by deriving Eigenvectors and Eigen values from the mean co-variance matrix.
17. Apparatus according to claim 12 wherein the processor is responsive to executing computer instructions to correlate the second portion of the cardiac electrogram signal with multiple templates that represent features of the physiological signal to produce a set of template correlated signals.
18. Apparatus according to claim 12 wherein the processor is responsive to executing computer instructions to derive a vector from a first time sample of each template correlated signal of the set; and further vectors from further time samples of each template correlated signal of the set.
19. Apparatus according to claim 18 wherein the processor is responsive to executing computer instructions to express the vectors as points in the model of noise by projecting each vector onto each Eigenvector thereby representing the original signal as a trajectory in the model of noise.
20. A method according to claim 11 comprising displaying the time domain signal overlaid with a noise threshold corresponding to the limits of the model of noise, portions of the time domain signal displayed above the noise threshold indicating they are physiologically significant.
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
GB1414330.9A GB2529200B (en) | 2014-08-13 | 2014-08-13 | Improvement to Analysing Physiological Electrograms |
PCT/GB2015/052190 WO2016024081A1 (en) | 2014-08-13 | 2015-07-29 | Improvement to analysing physiological electrograms |
GB1414330.9 | 2015-08-13 |
Publications (2)
Publication Number | Publication Date |
---|---|
US20170265766A1 US20170265766A1 (en) | 2017-09-21 |
US10349852B2 true US10349852B2 (en) | 2019-07-16 |
Family
ID=51629700
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US15/503,674 Active US10349852B2 (en) | 2014-08-13 | 2015-07-09 | Analysing physiological electrograms |
Country Status (9)
Country | Link |
---|---|
US (1) | US10349852B2 (en) |
EP (1) | EP3179905A1 (en) |
JP (1) | JP2017529136A (en) |
CN (1) | CN107249447A (en) |
AU (1) | AU2015303017A1 (en) |
CA (1) | CA2958026A1 (en) |
GB (1) | GB2529200B (en) |
RU (1) | RU2017107224A (en) |
WO (1) | WO2016024081A1 (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11865219B2 (en) | 2013-04-15 | 2024-01-09 | Ecolab Usa Inc. | Peroxycarboxylic acid based sanitizing rinse additives for use in ware washing |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105769173B (en) * | 2016-02-29 | 2018-08-28 | 浙江铭众科技有限公司 | A kind of cardioelectric monitor system with electrocardiosignal noise removal function |
CN105796091B (en) * | 2016-02-29 | 2018-09-28 | 浙江铭众科技有限公司 | A kind of intelligent terminal for removing electrocardiosignal vehicle movement noise |
CN105534517A (en) * | 2016-02-29 | 2016-05-04 | 浙江铭众科技有限公司 | Method for removing vehicle motion noise in three-lead electrocardiosignal |
CN110179456B (en) * | 2019-05-23 | 2021-11-02 | 中国航天员科研训练中心 | Electrocardio noise recognition model training and electrocardio noise detection method and device |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5172690A (en) | 1990-10-26 | 1992-12-22 | Telectronics Pacing Systems, Inc. | Automatic stimulus artifact reduction for accurate analysis of the heart's stimulated response |
WO1996025096A1 (en) | 1995-02-17 | 1996-08-22 | Ep Technologies, Inc. | Acquisition of endocardially or epicardially paced electrocardiograms |
US20040127949A1 (en) | 2002-12-31 | 2004-07-01 | Jaeho Kim | Method and system for detecting capture with cancellation of pacing artifact |
US20050149134A1 (en) | 2003-12-24 | 2005-07-07 | Mccabe Aaron | Method and system for removing pacing artifacts from subcutaneous electrocardiograms |
US20060253044A1 (en) | 2005-05-09 | 2006-11-09 | Yi Zhang | Arrhythmia discrimination using electrocardiograms sensed from multiple implanted electrodes |
GB2439562A (en) | 2006-06-29 | 2008-01-02 | Medilec Ltd | Analysing electrograms |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1556687A (en) * | 2001-09-21 | 2004-12-22 | Nonlinear noise reduction for magnetocardiograms using wavelet transforms | |
JP4761724B2 (en) * | 2004-04-21 | 2011-08-31 | アジレント・テクノロジーズ・インク | Method for measuring phase noise and phase noise measuring apparatus |
WO2007026266A2 (en) * | 2005-06-15 | 2007-03-08 | Koninklijke Philips Electronics N.V. | Noise model selection for emission tomography |
US9052230B2 (en) * | 2011-05-13 | 2015-06-09 | Chevron U.S.A. Inc | Industrial process monitoring and imaging |
CN102274029B (en) * | 2011-05-25 | 2013-10-02 | 中国科学院深圳先进技术研究院 | Identity recognition method and system |
CN102670188A (en) * | 2012-04-25 | 2012-09-19 | 重庆大学 | Robust adaptive estimation method for single-guide heart rate detection of fetus |
CN103006206B (en) * | 2013-01-11 | 2014-08-06 | 山东师范大学 | Method for quantitatively detecting microvolt T-wave alternans |
-
2014
- 2014-08-13 GB GB1414330.9A patent/GB2529200B/en active Active
-
2015
- 2015-07-09 US US15/503,674 patent/US10349852B2/en active Active
- 2015-07-29 CN CN201580053239.8A patent/CN107249447A/en active Pending
- 2015-07-29 RU RU2017107224A patent/RU2017107224A/en unknown
- 2015-07-29 WO PCT/GB2015/052190 patent/WO2016024081A1/en active Application Filing
- 2015-07-29 CA CA2958026A patent/CA2958026A1/en not_active Abandoned
- 2015-07-29 EP EP15747189.7A patent/EP3179905A1/en active Pending
- 2015-07-29 JP JP2017508563A patent/JP2017529136A/en active Pending
- 2015-07-29 AU AU2015303017A patent/AU2015303017A1/en not_active Abandoned
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5172690A (en) | 1990-10-26 | 1992-12-22 | Telectronics Pacing Systems, Inc. | Automatic stimulus artifact reduction for accurate analysis of the heart's stimulated response |
WO1996025096A1 (en) | 1995-02-17 | 1996-08-22 | Ep Technologies, Inc. | Acquisition of endocardially or epicardially paced electrocardiograms |
US20040127949A1 (en) | 2002-12-31 | 2004-07-01 | Jaeho Kim | Method and system for detecting capture with cancellation of pacing artifact |
US20050149134A1 (en) | 2003-12-24 | 2005-07-07 | Mccabe Aaron | Method and system for removing pacing artifacts from subcutaneous electrocardiograms |
US20060253044A1 (en) | 2005-05-09 | 2006-11-09 | Yi Zhang | Arrhythmia discrimination using electrocardiograms sensed from multiple implanted electrodes |
GB2439562A (en) | 2006-06-29 | 2008-01-02 | Medilec Ltd | Analysing electrograms |
Non-Patent Citations (2)
Title |
---|
International Search Report dated Oct. 5, 2015, and received for International Application No. PCT/GB2015/052190. |
Written Opinion of the International Search Authority dated Oct. 5, 2015, and received for International Application No. PCT/GB2015/052190. |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11865219B2 (en) | 2013-04-15 | 2024-01-09 | Ecolab Usa Inc. | Peroxycarboxylic acid based sanitizing rinse additives for use in ware washing |
Also Published As
Publication number | Publication date |
---|---|
WO2016024081A1 (en) | 2016-02-18 |
AU2015303017A1 (en) | 2017-03-23 |
GB2529200B (en) | 2019-05-01 |
GB201414330D0 (en) | 2014-09-24 |
JP2017529136A (en) | 2017-10-05 |
CA2958026A1 (en) | 2016-02-18 |
GB2529200A (en) | 2016-02-17 |
CN107249447A (en) | 2017-10-13 |
RU2017107224A (en) | 2018-09-13 |
US20170265766A1 (en) | 2017-09-21 |
EP3179905A1 (en) | 2017-06-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10349852B2 (en) | Analysing physiological electrograms | |
KR102451795B1 (en) | ECG signal detection method | |
CN108186011B (en) | Atrial fibrillation detection method, atrial fibrillation detection device and readable storage medium | |
Dhar et al. | Cross-wavelet assisted convolution neural network (AlexNet) approach for phonocardiogram signals classification | |
Safara et al. | Multi-level basis selection of wavelet packet decomposition tree for heart sound classification | |
Kao et al. | Automatic phonocardiograph signal analysis for detecting heart valve disorders | |
KR101910576B1 (en) | The apparutus and method for classify input pattern promptly using artificial neural network | |
JP5672144B2 (en) | Heart rate / respiration rate detection apparatus, method and program | |
CN108024730A (en) | Using mathematical analysis and machine learning come the method and system that diagnoses the illness | |
Panigrahy et al. | Extraction of fetal ECG signal by an improved method using extended Kalman smoother framework from single channel abdominal ECG signal | |
Phanphaisarn et al. | Heart detection and diagnosis based on ECG and EPCG relationships | |
CN104503592A (en) | Method and device for determining head gestures | |
CN110786849B (en) | Electrocardiosignal identity recognition method and system based on multi-view discriminant analysis | |
Sharma et al. | Stationary wavelet transform based technique for automated external defibrillator using optimally selected classifiers | |
Kamson et al. | Multi-centroid diastolic duration distribution based HSMM for heart sound segmentation | |
CN111358460A (en) | Arrhythmia identification method and device and electronic equipment | |
Zhu et al. | Feature extraction from a novel ECG model for arrhythmia diagnosis | |
JP2842409B2 (en) | Evoked waveform calculator | |
Liang et al. | Identification of heart sounds with arrhythmia based on recurrence quantification analysis and Kolmogorov entropy | |
CN104545890A (en) | Electrocardiosignal based identity recognition method and device | |
JP6825770B2 (en) | Machine learning equipment, classification equipment and methods and programs | |
CN115702778A (en) | Sleep stage staging method and device, electronic equipment and computer readable storage medium | |
Moukadem et al. | High Order Statistics and Time‐Frequency Domain to Classify Heart Sounds for Subjects under Cardiac Stress Test | |
Ravelo-García et al. | Cepstrum coefficients of the RR series for the detection of obstructive sleep apnea based on different classifiers | |
EP3019077A1 (en) | Classifying heart sounds |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
STPP | Information on status: patent application and granting procedure in general |
Free format text: NOTICE OF ALLOWANCE MAILED -- APPLICATION RECEIVED IN OFFICE OF PUBLICATIONS |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: PUBLICATIONS -- ISSUE FEE PAYMENT VERIFIED |
|
STCF | Information on status: patent grant |
Free format text: PATENTED CASE |
|
MAFP | Maintenance fee payment |
Free format text: PAYMENT OF MAINTENANCE FEE, 4TH YR, SMALL ENTITY (ORIGINAL EVENT CODE: M2551); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY Year of fee payment: 4 |