WO2007133920A2 - Method and system for real-time digital filtering for electrophysiological and hemodynamic amplifiers - Google Patents
Method and system for real-time digital filtering for electrophysiological and hemodynamic amplifiers Download PDFInfo
- Publication number
- WO2007133920A2 WO2007133920A2 PCT/US2007/067754 US2007067754W WO2007133920A2 WO 2007133920 A2 WO2007133920 A2 WO 2007133920A2 US 2007067754 W US2007067754 W US 2007067754W WO 2007133920 A2 WO2007133920 A2 WO 2007133920A2
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- digital data
- digital
- notch
- filtering
- data
- Prior art date
Links
Classifications
-
- 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/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/333—Recording apparatus specially adapted therefor
-
- 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/25—Bioelectric electrodes therefor
- A61B5/279—Bioelectric electrodes therefor specially adapted for particular uses
- A61B5/28—Bioelectric electrodes therefor specially adapted for particular uses for electrocardiography [ECG]
- A61B5/283—Invasive
-
- 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
- A61N1/3704—Circuits specially adapted therefor, e.g. for sensitivity control
Definitions
- the subject matter described herein relates to the field of invasive cardiology, and more particularly, to the field of signal filtering and conditioning.
- analog cardiac data is collected with heart catheters by electrophysiological and hemodynamic amplifiers and requires signal filtering and conditioning before it may be used.
- Current electrophysiological and hemodynamic amplifiers use cascaded analog circuits to perform signal filtering and conditioning. This signal filtering and conditioning is performed after analog anti-aliasing, but before analog-to- digital conversion.
- This (analog circuit based) signal filtering and conditioning process can be problematic, specifically, the amplified analog signals are consistently noisy and have poor signal resolution due to imprecise analog components. In fact, the precision of the analog components is approximately 1%.
- the systems utilized to process these analog signals are often times inflexible due to the unadj us table nature of the analog circuits.
- a method and system for real-time digital filtering for electrophysiological and hemodynamic amplifiers replaces the analog circuits currently used in such systems for signal filtering and conditioning with digital filters that may be implemented in a software application.
- the method and system includes digitizing the analog signal collected from the patient prior to performing the signal filtering and conditioning.
- the method and system also includes removing stimulus artifacts, as well as performing sample rate conversion and scaling on the digital signal.
- the processed digital signals may be used, displayed, saved and converted to analog signal through digital-to- analog conversion.
- a method of real-time digital filtering for electrophysiological and hemodynamic amplifiers comprises inputting a set of digital data from a patient, removing a stimulus artifact from the set of digital data when a pacing function is on, filtering the set of digital data with a first order high-pass filter and a second order low-pass filter or a special resonator, and performing a sample rate conversion and a data scaling on the set of digital data when a notch function is off.
- This method further comprises filtering the set of digital data with three first-order Butterworth notch filters at base, 2nd and 3rd harmonics of the power line frequency when the notch function is on, and the notch function is fixed, and filtering the set of digital data with three amplitude adaptive notch filters at base, 2nd and 3rd harmonics of the power line frequency when the notch function is on, and the notch function is not fixed, and collecting a set of analog cardiac data from the patient, and converting the analog cardiac data to the set of digital data wherein the first order high-pass filter and the second order low-pass filter are Butterworth type digital filers.
- the method further comprises outputting the set of digital data to an output device, converting the outputted set of digital data to an analog signal.
- the system comprises an acquisition device configured to collect a set of analog cardiac data from a patient, wherein the acquisition device is further configured to convert the set of analog cardiac data to a set of digital data, a storage media for storing a computer application and a processing unit coupled to the acquisition device and the storage media, the processing unit configured to execute the computer application, and further configured to receive the set of digital data from the acquisition device, wherein when the computer application is executed, a stimulus artifact is removed from the set of digital data when a pacing function is on, and the set of digital data is filtered with a first order high-pass filter and a second order low-pass filter or a special resonator, and a sample rate conversion and a data scaling is performed on the set of digital data when a notch function is off.
- the set of digital data is filtered with three first-order Butterworth notch filters at base, 2nd and 3rd harmonics of the power line frequency when the notch function is on, and the notch function is fixed, and the set of digital data is filtered with three amplitude adaptive notch filters at base, 2nd and 3rd harmonics of the power line frequency when the notch function is on, and the notch function is not fixed.
- the first order high-pass filter and the second order low-pass filter are Butterworth type digital filers.
- the system further comprises an output device configured to receive an outputted set of digital data wherein the outputted set of digital data is converted to an analog signal.
- Figure 1 shows a flow chart of an embodiment of a method of real-time digital filtering.
- Figure 2 illustrates a block diagram of an embodiment of a system for real-time digital filtering.
- the method and system improves current electrophysiological and hemodynamic amplifiers by replacing cascaded analog filtering circuits with cascaded digital filters, which may be implemented in computer software for signal filtering and conditioning based on mathematical algorithms.
- the analog high pass filter will be replaced with a digital high pass filter
- the analog low pass filter with a digital low pass filter
- the analog fixed notch filters with digital fixed notch filters.
- users will have the option to select amplitude adaptive notch filtering; and in order to compensate for signal attenuation from analog anti-aliasing and analog-to-digital conversion, a digital special resonator will be used for signal boosting.
- filter coefficients are calculated and updated.
- the sample data is removed of pacing stimulus artifacts if pacing is on, and processed by the high pass filter, then the low pass filter, or the special digital resonator depending on the user-selected low pass cutoff frequency, which determines signal attenuation from analog anti-aliasing and analog- to -digital conversion.
- the signal is then filtered by the notch filters, the fixed notch filters or amplitude adaptive notch filters depending on user's choice if notch is on, and finally the sample rate conversion and data scaling are performed before the sample is output.
- New electrophysiological and hemodynamic amplifiers clearly solve the above three problems because they provide cleaner signals and better signal resolution since digital algorithms have better precision (0.0001% round-off error), flexible systems since digital filter coefficients can be easily calculated, and they have potential in addressing application problems such as removal of pacing stimulus artifacts by adding additional digital filtering algorithms.
- a filtering method 10 of the present invention is depicted in the flow chart. After collecting the analog cardiac data from a patient with a catheter, the signal is digitized and the digital data is inputted in step 12. If a pacing function is on in step 14, the stimulus artifact is removed from the digital signal in step 16. Step 16 is performed by using exponential curve fitting for the pace noise, and then subtracting that portion from the digital signal. Before describing the subsequent steps of the filtering method 10, it is important to describe the various Butterworth filters utilized in this invention.
- the transfer functions H (z) for the nth order Butterworth low-pass and high-pass digital filters may be obtained by applying the following bilinear transformation:
- f c is the 3 dB cutoff frequency in Hz of the low-pass or high-pass digital filter and f s is the sampling frequency in Hz.
- N are a coefficients, b coefficients and the length or number of a or b coefficients of the filters ( a o ⁇ 0 ).
- step 18 the digital data is filtered with a first order Butterworth high-pass filter.
- the specific operation of the first order Butterworth high-pass filter is described below: 1 st Order Rntterworth High-pass ⁇ TR Filter
- step 20 the digital data is then filtered with a second-order Butterworth low-pass filter or a special digital resonator.
- a second-order Butterworth low-pass filter or a special digital resonator.
- a detailed description of the 2 n order Butterworth low-pass filter is included below as well.
- a sample rate conversion and data scaling of the digital data is performed in step 30.
- the system and method will utilize the Shannon's sampling theorem to reconstruct the signal from its digital original samples, at the original sample rate, and then resample it at the output sampling rate.
- This type of sample rate conversion will incoiporate anti-aliasing before down sampling and remove duplicate spectra after up-sampling after step 30, the digital data will be outputted in step 32.
- the digital data may be outputted to a digital to analog converter so that it may be relayed onto a monitor, monitoring system, or any type of physician workstation or physician device. Likewise, the outputted digital signal may be utilized without conversion to analog.
- step 24 it will be determined whether the notch is fixed. If the notch is fixed, then the digital data will be filtered with a first order Butterworth notch filter at the base in step 26. The specifics of this filtering method are shown below.
- step 28 the digital data is filtered with adaptive notch filters at the base.
- step 28 A detailed description of step 28 is outlined below.
- These filters have two portions: 1) detection of base power line and its harmonic frequencies; 2) adaptive amplitude filter using modeling technique.
- the frequency of power line is not very stable due to poor load control, where the frequency can be higher when load is low, and frequency can be lower when load is high. In this way, a fixed power line frequency notch filter may not be able to reduce the noise to acceptable level.
- the new filter scan a segment of the data and perform power spectrum analysis to detect the spectrum peak around 50 Hz and 60 Hz depend on the presetting of the regional code.
- a threshold is set to a 3 times of standard deviation of the power spectrum of the surrounding region. If the power spectrum is exceeding this threshold, the location of the peak value is detected as notch filter frequency.
- the 2 nd and 3 rd harmonic spectrum peaks are compared to the surrounding spectrum with the same 3 times standard deviation rule. If the spectrum peak exceeds the threshold, the harmonic notch filter(s) are set for the detected frequency. Otherwise, the harmonic filters are turned off to save the computation time.
- the adaptive amplitude filter is to model the sinusoid signal first and then subtract the filter from the signal.
- signal s(n+l) need to be adapted to correct modeling error after the first estimation.
- a small constant is added to or subtracted from the estimated s(n+l) depending on the polarity of the signal difference.
- K is the amplitude of power line interference noise at notch frequency, i.e., the power line base, 2 n harmonic or 3 r harmonic frequency. This K is calculated from its past two points using the following formula:
- ⁇ a is the pre-configured amplitude adaptive step size.
- an embodiment of the system 50 of the present invention is depicted in the block diagram. It should be noted that the method as discussed above in Fig. 1 may be implemented as a software application and executed by the system 50, herein described with respect to Fig. 2.
- an embodiment of the system 50 of the present invention includes an acquisition device 54 configured to obtain a set of cardiac data from a patient 52.
- the acquisition device 54 is further configured to convert the analog data into digital data and relay the digital data to the processor 56.
- the storage media 58 in Fig. 2 includes the computer code embodying the software configured to Ca 1 Ty out the method as described above.
- the processor executes the method on the collected digital data from the acquisition device 54 and provides the digital data output as described above in the method of Fig. 1 a technical effect of the system 50 is output of the digital data to the output device 62 for display on a graphical user interface 64 if desired.
- the output device 62 may be any digital or analog device used to view, measure, manipulate or record such cardiac data.
- a digital to analog converter will be utilized to convert the output data from the processor 56.
- the electrophysiological and hemodynamic amplifiers of the present invention have three major advantages in that the digital filtering algorithms are much less bulky than analog filter circuits, there is a much lower cost in manufacturing with free duplication of digital algorithms versus the costly analog component, and the present invention is a flexible and extendable system capable of adding features and addressing application problems by using additional digital algorithms.
Abstract
The present invention includes a method and system for real-time digital filtering for electrophysiological and hemodynamic amplifiers. The invention replaces the analog circuits currently used for signal filtering and conditioning in such systems with digital filters that may be implemented in a software application. The method and system includes digitizing the analog signal collected from the patient prior to performing the signal filtering and conditioning. The method and system also includes removing stimulus artifacts (16), as well as performing sample rate conversion and scaling (30) on the digital signal. The processed digital signals may be used, displayed, saved and converted to analog signal thru digital-to-analog conversion.
Description
METHOD AND SYSTEM FOR REAL-TIME DIGITAL FILTERING FOR ELECTROPHYSIOLOGICAL AND HEMODYNAMIC AMPLIFIERS
RELATED APPLICATIONS
The present application is based on and claims priority to U.S. Provisional Patent Application Serial No. 60/800,372, filed May 15, 2006 and to U.S. Non-Provisional Patent Application Serial No. 11/467,361, filed August 25, 2006.
FIELD OF THE INVENTION
The subject matter described herein relates to the field of invasive cardiology, and more particularly, to the field of signal filtering and conditioning.
BACKGROUND OF THE INVENTION
In the field of invasive cardiology, analog cardiac data is collected with heart catheters by electrophysiological and hemodynamic amplifiers and requires signal filtering and conditioning before it may be used. Current electrophysiological and hemodynamic amplifiers use cascaded analog circuits to perform signal filtering and conditioning. This signal filtering and conditioning is performed after analog anti-aliasing, but before analog-to- digital conversion. This (analog circuit based) signal filtering and conditioning process can be problematic, specifically, the amplified analog signals are consistently noisy and have poor signal resolution due to imprecise analog components. In fact, the precision of the analog components is approximately 1%. Furthermore, the systems utilized to process these analog signals are often times inflexible due to the unadj us table nature of the analog circuits. Lastly, the analog amplifiers currently used make it difficult or even impossible to address application problems such as removal of pacing stimulus artifacts in the signal. What is needed is a method and system of real-time digital filtering for these amplifiers that produce better quality signals, flexibility in system structure, and additional functionality such as removing pacing stimulus artifacts from the signal.
. l .
BRIEF DESCRIPTION OF THE INVENTION
A method and system for real-time digital filtering for electrophysiological and hemodynamic amplifiers is provided. The system replaces the analog circuits currently used in such systems for signal filtering and conditioning with digital filters that may be implemented in a software application. The method and system includes digitizing the analog signal collected from the patient prior to performing the signal filtering and conditioning. The method and system also includes removing stimulus artifacts, as well as performing sample rate conversion and scaling on the digital signal. The processed digital signals may be used, displayed, saved and converted to analog signal through digital-to- analog conversion.
In one embodiment a method of real-time digital filtering for electrophysiological and hemodynamic amplifiers comprises inputting a set of digital data from a patient, removing a stimulus artifact from the set of digital data when a pacing function is on, filtering the set of digital data with a first order high-pass filter and a second order low-pass filter or a special resonator, and performing a sample rate conversion and a data scaling on the set of digital data when a notch function is off. This method further comprises filtering the set of digital data with three first-order Butterworth notch filters at base, 2nd and 3rd harmonics of the power line frequency when the notch function is on, and the notch function is fixed, and filtering the set of digital data with three amplitude adaptive notch filters at base, 2nd and 3rd harmonics of the power line frequency when the notch function is on, and the notch function is not fixed, and collecting a set of analog cardiac data from the patient, and converting the analog cardiac data to the set of digital data wherein the first order high-pass filter and the second order low-pass filter are Butterworth type digital filers. The method further comprises outputting the set of digital data to an output device, converting the outputted set of digital data to an analog signal.
Another embodiment of a system for real-time digital filtering for electrophysiological and hemodynamic amplifiers is provided. The system comprises an acquisition device configured to collect a set of analog cardiac data from a patient, wherein the acquisition
device is further configured to convert the set of analog cardiac data to a set of digital data, a storage media for storing a computer application and a processing unit coupled to the acquisition device and the storage media, the processing unit configured to execute the computer application, and further configured to receive the set of digital data from the acquisition device, wherein when the computer application is executed, a stimulus artifact is removed from the set of digital data when a pacing function is on, and the set of digital data is filtered with a first order high-pass filter and a second order low-pass filter or a special resonator, and a sample rate conversion and a data scaling is performed on the set of digital data when a notch function is off. The set of digital data is filtered with three first-order Butterworth notch filters at base, 2nd and 3rd harmonics of the power line frequency when the notch function is on, and the notch function is fixed, and the set of digital data is filtered with three amplitude adaptive notch filters at base, 2nd and 3rd harmonics of the power line frequency when the notch function is on, and the notch function is not fixed. The first order high-pass filter and the second order low-pass filter are Butterworth type digital filers. The system further comprises an output device configured to receive an outputted set of digital data wherein the outputted set of digital data is converted to an analog signal.
BRIEF DESCRIPTION OF THE DRAWINGS
Figure 1 shows a flow chart of an embodiment of a method of real-time digital filtering.
Figure 2 illustrates a block diagram of an embodiment of a system for real-time digital filtering.
DETAILED DESCRIPTION
The method and system improves current electrophysiological and hemodynamic amplifiers by replacing cascaded analog filtering circuits with cascaded digital filters, which may be implemented in computer software for signal filtering and conditioning based on mathematical algorithms. Specifically, the analog high pass filter will be replaced with a digital high pass filter, the analog low pass filter with a digital low pass filter, and the analog fixed notch filters with digital fixed notch filters.
In addition, in order to improve removal of power line noise while preserving signal contents around the power line frequency and the corresponding 2n and 3r harmonics, users will have the option to select amplitude adaptive notch filtering; and in order to compensate for signal attenuation from analog anti-aliasing and analog-to-digital conversion, a digital special resonator will be used for signal boosting.
At the start of the electrophysiological and hemodynamic application, or when the user changes system configuration parameters that affect the digital filtering, filter coefficients are calculated and updated. Once a digital sample is output from the analog-to-digital converter, the sample data is removed of pacing stimulus artifacts if pacing is on, and processed by the high pass filter, then the low pass filter, or the special digital resonator depending on the user-selected low pass cutoff frequency, which determines signal attenuation from analog anti-aliasing and analog- to -digital conversion. The signal is then filtered by the notch filters, the fixed notch filters or amplitude adaptive notch filters depending on user's choice if notch is on, and finally the sample rate conversion and data scaling are performed before the sample is output.
New electrophysiological and hemodynamic amplifiers clearly solve the above three problems because they provide cleaner signals and better signal resolution since digital algorithms have better precision (0.0001% round-off error), flexible systems since digital filter coefficients can be easily calculated, and they have potential in addressing application problems such as removal of pacing stimulus artifacts by adding additional digital filtering algorithms.
Referring to Fig. 1, an embodiment of a filtering method 10 of the present invention is depicted in the flow chart. After collecting the analog cardiac data from a patient with a catheter, the signal is digitized and the digital data is inputted in step 12. If a pacing function is on in step 14, the stimulus artifact is removed from the digital signal in step 16. Step 16 is performed by using exponential curve fitting for the pace noise, and then
subtracting that portion from the digital signal. Before describing the subsequent steps of the filtering method 10, it is important to describe the various Butterworth filters utilized in this invention.
Bntferworfh High-pfiss. T .ow-pass and Notch Filters
The transfer function H(s) for the nth order Butterworth low-pass analog filter is generally given in Eq. Al below:
Where [n/2] is the integer part of n/2, P - 0 or 1 depending on whether n is even or odd, Ωc is the analog 3 dB cutoff frequency. This general formula gives an expression whose denominator is the product of 2nd order polynomials in s. When n is odd, there is a 1st order polynomials also.
Replacing s/Ωc withΩf /s m above Eq. Al, we obtain the transfer function H(s) for the nth order Butterworth high-pass analog filter:
From Equations A 1-2, the transfer functions H (z) for the nth order Butterworth low-pass and high-pass digital filters may be obtained by applying the following bilinear transformation:
Where fc is the 3 dB cutoff frequency in Hz of the low-pass or high-pass digital filter and fs is the sampling frequency in Hz.
Where «/ , fy (1 = 0 to N-I) and N are a coefficients, b coefficients and the length or number of a or b coefficients of the filters (ao ≠ 0 ).
The difference Eq. describing the output response yn to input sample xn of the filter is readily obtained (same as Eq. 1) below:
Finally, the transfer function H(z) of the 1st order Butterworth band-stop or notch digital filter is obtained by applying the following transformation to the transfer function of the Is order Butterworth low-pass digital filter:
Where,
(Zr = the 3 dB cutoff frequency in Hz of the original low-pass filter, /„ , /, = upper, lower 3 dB cutoff frequencies in Hz of the new notch filter and Is = sampling rate in Hz).
Referring back to Fig. 1, in step 18, the digital data is filtered with a first order Butterworth high-pass filter. The specific operation of the first order Butterworth high-pass filter is described below:
1 st Order Rntterworth High-pass ΪTR Filter
The mathematical formula for lsl order Butterworth high-pass IIR digital filter is illustrated below:
Where,
( fc = 3 dB cutoff frequency in Hz and = sampling rate in Hz).
In step 20, the digital data is then filtered with a second-order Butterworth low-pass filter or a special digital resonator. A detailed description of the 2n order Butterworth low-pass filter is included below as well.
?.nd Order Bπtfe.rworth ϊ ,nw-pass TTR Filter
The mathematical formula for 2nd order Butterworth low-pass IIR digital filter is illustrated below:
Where,
( /,. = 3 dB cutoff frequency in Hz and Js = sampling rate in Hz).
Still refeπing to Fig. 1 , if the notch function is off in step 22, then a sample rate conversion and data scaling of the digital data is performed in step 30. The system and method will utilize the Shannon's sampling theorem to reconstruct the signal from its digital original samples, at the original sample rate, and then resample it at the output sampling rate. This type of sample rate conversion will incoiporate anti-aliasing before down sampling and remove duplicate spectra after up-sampling after step 30, the digital data will be outputted in step 32. The digital data may be outputted to a digital to analog converter so that it may be relayed onto a monitor, monitoring system, or any type of physician workstation or physician device. Likewise, the outputted digital signal may be utilized without conversion to analog. Referring back to step 22 of Fig. 1, if the notch function is on, then in step 24, it will be determined whether the notch is fixed. If the notch is fixed, then the digital data will be filtered with a first order Butterworth notch filter at the base in step 26. The specifics of this filtering method are shown below.
1st Orrle.r Rntterworth Rand-stop (Nntr.η) TTR Filter
The mathematical formula for lsl order Butterworth band-stop (notch) IIR digital filter is illustrated below:
Where,
(/„ , /, = upper, lower 3 dB cutoff frequencies in Hz and Js = sampling rate in Hz). After the filtering step of step 26, the method will continue to step 30.
Referring back to step 24 of Fig. 1, if the notch is not fixed, then the method continues on to step 28, where the digital data is filtered with adaptive notch filters at the base. A detailed description of step 28 is outlined below.
Amplitude Adaptive Nntπh Filters with Power T .ine Frequency Dp.tp.rti on
These filters have two portions: 1) detection of base power line and its harmonic frequencies; 2) adaptive amplitude filter using modeling technique.
In quite a few countries, the frequency of power line is not very stable due to poor load control, where the frequency can be higher when load is low, and frequency can be lower when load is high. In this way, a fixed power line frequency notch filter may not be able to reduce the noise to acceptable level. The new filter scan a segment of the data and perform power spectrum analysis to detect the spectrum peak around 50 Hz and 60 Hz depend on the presetting of the regional code. A threshold is set to a 3 times of standard deviation of the power spectrum of the surrounding region. If the power spectrum is exceeding this threshold, the location of the peak value is detected as notch filter frequency.
After the base frequency /p is detected, the 2nd and 3rd harmonic spectrum peaks are compared to the surrounding spectrum with the same 3 times standard deviation rule. If the spectrum peak exceeds the threshold, the harmonic notch filter(s) are set for the detected frequency. Otherwise, the harmonic filters are turned off to save the computation time.
The adaptive amplitude filter is to model the sinusoid signal first and then subtract the filter from the signal. A trigonometric identity sin(a+b) = 2sin(a)cos(b) - sin(a-b). That lead to the time series sinusoid signal model as: s(n+l) = 2*C*s(n) - s(n-l), where C =
cos(2πf/fs), f = sinusoid signal frequency, fs is the sample frequency. In the real-time filtering, signal s(n+l) need to be adapted to correct modeling error after the first estimation. A small constant is added to or subtracted from the estimated s(n+l) depending on the polarity of the signal difference.
The mathematical formula for Amplitude Adaptive Notch Filter is illustrated below:
Where K is the amplitude of power line interference noise at notch frequency, i.e., the power line base, 2n harmonic or 3r harmonic frequency. This K is calculated from its past two points using the following formula:
Where,
is notch frequency dependent coefficient ( fp = notch frequency in Hz andΛ = sampling rate in Hz).
Finally, the signal amplitude difference between the current and last sample is checked in order to provide a better estimate of the amplitude of the power line interference noise for next point using the following formula:
Where δa is the pre-configured amplitude adaptive step size.
After the filtering step is completed in step 28, the method continues on to step 30.
Referring now to Fig. 2, an embodiment of the system 50 of the present invention is depicted in the block diagram. It should be noted that the method as discussed above in Fig. 1 may be implemented as a software application and executed by the system 50, herein described with respect to Fig. 2. Referring now to Fig. 2, an embodiment of the system 50 of the present invention includes an acquisition device 54 configured to obtain a set of cardiac data from a patient 52. The acquisition device 54 is further configured to convert the analog data into digital data and relay the digital data to the processor 56. The storage media 58 in Fig. 2 includes the computer code embodying the software configured to Ca1Ty out the method as described above. Executing the computer code, the processor performs the method on the collected digital data from the acquisition device 54 and provides the digital data output as described above in the method of Fig. 1 a technical effect of the system 50 is output of the digital data to the output device 62 for display on a graphical user interface 64 if desired. As noted above, the output device 62 may be any digital or analog device used to view, measure, manipulate or record such cardiac data. In the case of the output device 62 being an analog type device, a digital to analog converter will be utilized to convert the output data from the processor 56.
The electrophysiological and hemodynamic amplifiers of the present invention have three major advantages in that the digital filtering algorithms are much less bulky than analog filter circuits, there is a much lower cost in manufacturing with free duplication of digital algorithms versus the costly analog component, and the present invention is a flexible and extendable system capable of adding features and addressing application problems by using additional digital algorithms.
By replacing analog filtering circuits with real-time digital filtering, new electrophysiological and hemodynamic amplifiers provide cleaner signals and better signal resolution, flexibility and extensibility in signal filtering and conditioning, and makes solving some tough technical problems such as removal of pacing stimulus artifacts possible.
The present invention has been described in terms of specific embodiments incorporating details to facilitate the understanding of the principals of construction and operation of the invention. Such reference herein to specific embodiments and details thereof is not intended to limit the scope of the claims appended hereto. It will be apparent to those skilled in the art that modifications may be made in the embodiment chosen for illustration without departing from the spirit and scope of the invention.
Claims
1. A method of real-time digital filtering for electrophysiological and hemodynamic amplifiers, the method comprising: inputting a set of digital data from a patient; removing a stimulus artifact from the set of digital data when a pacing function is on; filtering the set of digital data with a first order high-pass filter and a second order low-pass filter or a special resonator; and performing a sample rate conversion and a data scaling on the set of digital data when a notch function is off.
2. The method as claimed in claim 1, further comprising filtering the set of digital data with three first-order Butterworth notch filters at base, 2nd and 3rd harmonics of the power line frequency when the notch function is on, and the notch function is fixed.
3. The method as claimed in claim 1, further comprising filtering the set of digital data with three amplitude adaptive notch filters at base, 2nd and 3rd harmonics of the power line frequency when the notch function is on, and the notch function is not fixed.
4. The method as claimed in claim 1, further comprising collecting a set of analog cardiac data from the patient, and converting the analog cardiac data to the set of digital data.
5. The method as claimed in claim 1, wherein the first order high-pass filter and the second order low-pass filter are Butterwoith type digital filters.
6. The method as claimed in claim 1, further comprising outputting the set of digital data to an output device.
7. The method as claimed in claim 6, further comprising converting the outputted set of digital data to an analog signal.
8. The method as claimed in claim 1, further comprising storing the set of digital data in a storage medium.
9. A system for real-time digital filtering for electrophysiological and hemodynamic amplifiers, the system comprising: an acquisition device configured to collect a set of analog cardiac data from a patient, wherein the acquisition device is further configured to convert the set of analog cardiac data to a set of digital data; a storage media for storing a computer application; and a processing unit coupied to the acquisition device and the storage media, the processing unit configured to execute the computer application, and further configured to receive the set of digital data from the acquisition device, wherein when the computer application is executed, a stimulus artifact is removed from the set of digital data when a pacing function is on, and the set of digital data is filtered with a first order high-pass filter and a second order low-pass filter or a special resonator, and a sample rate conversion and a data scaling is performed on the set of digital data when a notch function is off.
10. The system as claimed in claim 9, wherein the set of digital data is filtered with three first-order Butterworth notch filters at base, 2nd and 3rd harmonics of the power line frequency when the notch function is on, and the notch function is fixed.
11. The system as claimed in claim 9, wherein the set of digital data is filtered with three amplitude adaptive notch filters at base, 2nd and 3rd harmonics of the power line frequency when the notch function is on, and the notch function is not fixed.
12. The system as claimed in claim 9, wherein the first order high-pass filter and the second order low-pass filter are Butterworth type digital filers.
13. The system as claimed in claim 9, further comprising an output device configured to receive an outputted set of digital data.
14. The system as claimed in claim 13, wherein the outputted set of digital data is converted to an analog signal.
15. A method of real-time digital filtering, the method comprising: collecting a set of analog cardiac data from a patient; converting the analog cardiac data to a set of digital data; removing a stimulus artifact from the set of digital data when a pacing function is on; filtering the set of digital data with a first order Butterworth high-pass filter and a second order Butterworth low-pass filter or a special resonator; performing a sample rate conversion and a data scaling on the set of digital data when a notch function is off; filtering the set of digital data with three first-order Butterworth notch filters at base,
2nd and 3rd harmonics of the power line frequency when the notch function is on, and the notch function is fixed; filtering the set of digital data with three amplitude adaptive notch filters at base, 2nd and 3rd harmonics of the power line frequency when the notch function is on, and the notch function is not fixed; and outputting an analog signal derived from the set of digital data.
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US80037206P | 2006-05-15 | 2006-05-15 | |
US60/800,372 | 2006-05-15 | ||
US11/467,361 | 2006-08-25 | ||
US11/467,361 US7783339B2 (en) | 2006-05-15 | 2006-08-25 | Method and system for real-time digital filtering for electrophysiological and hemodynamic amplifers |
Publications (2)
Publication Number | Publication Date |
---|---|
WO2007133920A2 true WO2007133920A2 (en) | 2007-11-22 |
WO2007133920A3 WO2007133920A3 (en) | 2008-10-02 |
Family
ID=38686029
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/US2007/067754 WO2007133920A2 (en) | 2006-05-15 | 2007-04-30 | Method and system for real-time digital filtering for electrophysiological and hemodynamic amplifiers |
Country Status (2)
Country | Link |
---|---|
US (1) | US7783339B2 (en) |
WO (1) | WO2007133920A2 (en) |
Families Citing this family (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8201000B2 (en) * | 2009-02-24 | 2012-06-12 | International Business Machines Corporation | Computing load management method and system |
US8205106B2 (en) * | 2009-02-24 | 2012-06-19 | International Business Machines Corporation | Energy load management method and system |
CN103532412A (en) * | 2012-07-05 | 2014-01-22 | 台达电子工业股份有限公司 | Feedback control circuit of power converter, and power converter system |
US10085662B2 (en) | 2012-10-30 | 2018-10-02 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Device and method for digital-to-analog transformations and reconstructions of multi-channel electrocardiograms |
CN103040459B (en) * | 2013-01-05 | 2014-06-04 | 西安交通大学 | Method of high-fidelity filtering for power frequency interferences in multichannel feeble physiological information recording system |
EP4218559A1 (en) | 2014-02-25 | 2023-08-02 | ICU Medical, Inc. | Patient monitoring system with gatekeeper signal |
JP6674553B2 (en) | 2015-10-19 | 2020-04-01 | アイシーユー・メディカル・インコーポレーテッド | Hemodynamic monitoring system with detachable display unit |
EP4278982A3 (en) | 2016-07-19 | 2024-02-21 | NuVera Medical, Inc. | Medical devices and methods of use |
US10525264B2 (en) | 2016-11-16 | 2020-01-07 | Samsung Electronics Co., Ltd. | Stimulator and method of controlling stimulator |
CA3057946A1 (en) | 2017-03-30 | 2018-10-04 | Shifamed Holdings, Llc | Medical tool positioning devices, systems, and methods of use and manufacture |
US10783223B2 (en) * | 2017-10-12 | 2020-09-22 | Pacesetter, Inc. | Hybrid signal processing circuit for implantable medical devices and methods |
US10686715B2 (en) | 2018-05-09 | 2020-06-16 | Biosig Technologies, Inc. | Apparatus and methods for removing a large-signal voltage offset from a biomedical signal |
IL281028B2 (en) | 2018-08-23 | 2024-02-01 | Nuvera Medical Inc | Medical tool positioning devices, systems, and methods of use and manufacture |
EP3880285B1 (en) | 2018-11-13 | 2024-02-21 | NuVera Medical, Inc. | Medical tool positioning devices |
US11576601B2 (en) * | 2019-04-18 | 2023-02-14 | X Development Llc | Artifact identification in EEG measurements |
US11419556B2 (en) * | 2020-07-02 | 2022-08-23 | Vital Connect, Inc. | Method and system for pacing pulse detection and pacing artifact rejection |
CN112754495A (en) * | 2021-01-19 | 2021-05-07 | 湖南万脉医疗科技有限公司 | Portable wearable non-inductive electrocardiogram monitoring system based on cloud computing technology |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4944299A (en) * | 1989-08-08 | 1990-07-31 | Siemens-Pacesetter, Inc. | High speed digital telemetry system for implantable device |
US5197467A (en) * | 1992-06-22 | 1993-03-30 | Telectronics Pacing Systems, Inc. | Multiple parameter rate-responsive cardiac stimulation apparatus |
US5400795A (en) * | 1993-10-22 | 1995-03-28 | Telectronics Pacing Systems, Inc. | Method of classifying heart rhythms by analyzing several morphology defining metrics derived for a patient's QRS complex |
US5494042A (en) * | 1994-01-28 | 1996-02-27 | Ep Technologies, Inc. | Systems and methods for deriving electrical characteristics of cardiac tissue for output in iso-characteristic displays |
US5647379A (en) * | 1994-11-22 | 1997-07-15 | Ventritex, Inc. | Correlator based electromagnetic interference responsive control system useful in medical devices |
US5795303A (en) * | 1995-02-17 | 1998-08-18 | Ep Technologies, Inc. | Systems and methods for making time-sequential measurements of biopotentials sensed in myocardial tissue |
WO2003043494A1 (en) * | 2001-11-23 | 2003-05-30 | Medit As | A cluster system for remote monitoring and diagnostic support |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5660184A (en) * | 1995-05-15 | 1997-08-26 | Johnson & Johnson Medical, Inc. | Pacemaker pulse detection and artifact rejection |
US6351664B1 (en) * | 1999-11-12 | 2002-02-26 | Ge Medical Systems Information Technologies, Inc. | Method of removing signal interference from sampled data and apparatus for effecting the same |
US20030194030A1 (en) * | 2002-04-15 | 2003-10-16 | Cory Samuel Modlin | Power line interference compensation in PCM upstream transmissions |
US20040199082A1 (en) * | 2003-04-03 | 2004-10-07 | Ostroff Alan H. | Selctable notch filter circuits |
-
2006
- 2006-08-25 US US11/467,361 patent/US7783339B2/en active Active
-
2007
- 2007-04-30 WO PCT/US2007/067754 patent/WO2007133920A2/en active Application Filing
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4944299A (en) * | 1989-08-08 | 1990-07-31 | Siemens-Pacesetter, Inc. | High speed digital telemetry system for implantable device |
US5197467A (en) * | 1992-06-22 | 1993-03-30 | Telectronics Pacing Systems, Inc. | Multiple parameter rate-responsive cardiac stimulation apparatus |
US5400795A (en) * | 1993-10-22 | 1995-03-28 | Telectronics Pacing Systems, Inc. | Method of classifying heart rhythms by analyzing several morphology defining metrics derived for a patient's QRS complex |
US5494042A (en) * | 1994-01-28 | 1996-02-27 | Ep Technologies, Inc. | Systems and methods for deriving electrical characteristics of cardiac tissue for output in iso-characteristic displays |
US5647379A (en) * | 1994-11-22 | 1997-07-15 | Ventritex, Inc. | Correlator based electromagnetic interference responsive control system useful in medical devices |
US5795303A (en) * | 1995-02-17 | 1998-08-18 | Ep Technologies, Inc. | Systems and methods for making time-sequential measurements of biopotentials sensed in myocardial tissue |
WO2003043494A1 (en) * | 2001-11-23 | 2003-05-30 | Medit As | A cluster system for remote monitoring and diagnostic support |
Non-Patent Citations (1)
Title |
---|
DATABASE INSPEC [Online] THE INSTITUTION OF ELECTRICAL ENGINEERS, STEVENAGE, GB; January 2006 (2006-01), RAMASHRI T ET AL: "Data compression using FFT and digital filters in PCM telemetry applications" XP002480512 Database accession no. 8918983 & Journal of the Institution of Engineers (India) Electronics and Telecommunication Engineering Division Instn. Eng. (India) India, vol. 86, January 2006 (2006-01), pages 63-67, XP002480511 ISSN: 0251-1096 * |
Also Published As
Publication number | Publication date |
---|---|
US20070265537A1 (en) | 2007-11-15 |
US7783339B2 (en) | 2010-08-24 |
WO2007133920A3 (en) | 2008-10-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7783339B2 (en) | Method and system for real-time digital filtering for electrophysiological and hemodynamic amplifers | |
JP5191791B2 (en) | Elimination of coherent signals in ECG | |
TW201034624A (en) | ECG device with impulse and channel switching ADC noise filter and error corrector for derived leads | |
Pandey et al. | High frequency noise removal from ECG using moving average filters | |
CN106687047B (en) | Sensor interface device providing digital processing of intravascular flow and pressure data | |
RU2677007C2 (en) | Apparatus and method for ecg motion artifact removal | |
KR20140139564A (en) | Systems and methods for ecg monitoring | |
CN110680308A (en) | Electrocardiosignal denoising method based on fusion of improved EMD and threshold method | |
CN107303177A (en) | The detection method and system of a kind of ECG T wave, P ripples | |
WO2014080571A1 (en) | Recursive least squares adaptive acoustic signal filtering for physiological monitoring system | |
Kumar et al. | Performance comparison of windowing techniques for ECG signal enhancement | |
CN110236528B (en) | Method and device for acquiring respiratory information | |
El Bouny et al. | Removal of 50Hz PLI from ECG signal using undecimated wavelet transform | |
Vityazeva et al. | Accuracy Loss in Multi-Rate Processing of Biomedical Signals | |
Gupta et al. | Hardware implementation and reduction of artifacts from ECG signal | |
Uma et al. | ECG Noise Removal Using Modified Distributed Arithmetic Based Finite Impulse Response Filter | |
Toan et al. | Design an open DSP-based system to acquire and process the bioelectric signal in realtime | |
Zhang et al. | A Method for Reduction of Noise in the ECG | |
H Abbas | Removing 0.5 Hz baseline wander from ECG signal using multistage adaptive filter | |
Maghfiroh et al. | Design of Single Lead Electrocardiography using Filter Order 3 to Reduce Noise Using Spektrum Analysis Based on Fast Fourier Transform | |
Mahdiani et al. | A review of transient suppression methods of IIR notch filters used for power-line interference rejection in ECG measurement | |
Sonal et al. | Comparative analysis of ECG baseline wander correction techniques | |
CN114617544A (en) | Breathing abnormity detection method based on frequency domain analysis | |
Huang et al. | Muscle Artifacts Cancellation Framework for ECG Signals Combining Convolution Auto-encoder and Average Beat Subtraction | |
Patel et al. | Nyquist Filter Design to Remove Power Line Interference From Bio-Electric Signals |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 07761561 Country of ref document: EP Kind code of ref document: A2 |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 07761561 Country of ref document: EP Kind code of ref document: A2 |