1. Technical Field
The invention generally relates to a device and a method for signal processing, and more particularly, to a device and a method for signal processing which can selectively initiate a reduction process according to the characteristic of noise components present within an input signal.
2. Description of Related Art
With the development of biomedical engineering, doctors are able to make convincing diagnoses according to different sorts of medical reports produced by biomedical inspection instruments. For example, the instruments may obtain a physiologic signal, such as an electrocardiogram (ECG) signal, through electrodes directly attached to a patient's chest, arm, or leg. Then, the physiological signals can be detected and transmitted to processing device for further processing.
Typically, the amplitude of the ECG signal transmitted from electrodes is usually measured around millivolts, so that noise components within the ECG signal may significantly affect the actual waveform of the ECG signal. Thus, it is almost meaningless to analyze the ECG signal with the presence of noise components in the ECG waveforms. In order to deal with this sorts of dilemma, high pass filters or band pass filters are generally used to conduct a noise reduction process to reduce the noise components caused by environment, such as power-line interference, ambient electromagnetic activity or movements of electrodes. Besides, signal processing techniques are also considered useful for reducing the noise components of the ECG signal, such as baseline wander effect and motion artifact. As reducing noise components within the ECG signal may be one of the most critical biomedical issues nowadays, it is worthy to develop a signal processing device and a signal processing method which can filter out the noise components.
According to an embodiment of the invention, a signal processing device is disclosed. The signal processing device includes a sampling module, a first segmentation module, a second segmentation module, and a detection module. The sampling module samples an input signal to generate a sample signal. The first segmentation module calculates a first segment value according to the sample signal during a first time interval. The second segmentation module calculates a second segment value according to the sample signal during a second time interval, wherein the length of the first time interval is different from that of the second time interval. The detection module is coupled to the first segmentation module and the second segmentation module to generate a detection signal according to the determination of whether the first segment value lies out of a first range, and whether the second segment value lies out of a second range.
According to another embodiment of the invention, a signal processing method is disclosed. The signal processing method includes the following steps. First, a sample signal is generated by sampling an input signal. Then, a first segment value is calculated according to the sample signal during a first time interval. And, a second segment value is calculated according to the sample signal during a second time interval, wherein the length of the first time interval is different from that of the second time interval. Further, a detection signal is generated according to the determination of whether the first segment value lies out of a first range, and whether the second segment value lies out of a second range.
According to another embodiment of the invention, a signal processing device is disclosed. The signal processing device includes a sampling module, a first segmentation module, a second segmentation module, and a reduction module. The sampling module samples an input signal to generate a sample signal. The first segmentation module calculates a first segment value according to the sample signal during a first time interval. The second segmentation module calculates a second segment value according to the sample signal during a second time interval, wherein the length of the first time interval is different from that of the second time interval. The reduction module is coupled to the sampling module, the first segmentation module and the second segmentation module to generate an output signal according to the sample signal, the first segment value and the second segment value.
According to still another embodiment of the invention, a signal processing method is disclosed. The signal processing method includes the following steps. First, a sample signal is generated by sampling an input signal. Then, a first segment value is calculated according to the sample signal during a first time interval. And, a second segment value is calculated according to the sample signal during a second time interval, wherein the length of the first time interval is different from that of the second time interval. Further, an output signal is generated according to the sample signal, the first segment value and the second segment value.
In order to further the understanding regarding the invention, the following embodiments are provided along with illustrations to facilitate the disclosure of the invention.
The aforementioned illustrations and following detailed description are exemplary for the purpose of illustrating the invention, but not limiting the scope of the invention. Some objectives and advantages related to the invention will be illustrated in the subsequent description and appended drawings.
As an example for illustrating the operation of the sampling module 10, please refer to
S2={S1[δt], S1[2δt], S1[3δt], S1[4δt], . . . } (1)
Note that the value of δt should be chosen small enough to preserve the characteristic of the input signal S1. For instance, if the input signal S1 is an ECG signal, the at might be within several milliseconds.
Please refer to
In one embodiment, the first segment value S3 is generated by averaging the plurality of sample data of the sample signal S2 within the first time interval. Of course, the first segment value S3 may be calculated in other mathematical methods as long as the amplitude characteristics of the sample signal S2 is preserved. In addition, the period of the first time interval may be set by considering the characteristic of the input signal. For example, if the input signal S1 is a ECG signal with baseline wander effect, the first time interval ΔT may be approximately 100˜200 ms. As such, during the first time interval ΔT, there might be tens of samples of the sample signal S2.
Please refer to
In one embodiment, the second segment value S4 is generated by averaging the sample data of the sample signal S2 within the second time interval. Of course, the second segment value S4 may be calculated in other mathematical methods as long as the amplitude characteristics of the sample signal S2 is preserved.
Please refer to
As a common application, the input signal S1 can be an ECG signal containing different types of noise. Noise, such as motion artifact effect, might have a frequency range above 0.5 Hz while baseline wander might have a lower frequency and may affect the ECG signal for more than half a second. On the other hand, for patients with certain heart disease, the morphology of the ECG might have irregular variation, usually having higher frequency than the baseline wander effect. By properly selecting the first time interval, the second time interval, the first range and the second range, it might be helpful in distinguishing noise from the morphology variation through the first segment value and the second segment value. Generally speaking, the low frequency noise will last a longer period of time. Therefore, if a low frequency noise, e.g. baseline wander effect, exists in the ECG signal, both of the short term first segment value S3 and the long term segment value S4 will be laid out of certain amplitude ranges. Thus, the noise existence can be detected by determining whether the short term first segment value S3 lies out of the first range and the corresponding contemporary long term second segment value S4 lies out the second range. In other words, if the short term first segment value S3 lies out of the first range and the corresponding contemporary long term second segment value S4 lies out the second range, e.g. detection signal at “1”, it is possible that the input signal S1 may bear a low frequency noise-like effect. If the short term first segment value S3 lies out of the first range but the corresponding contemporary long term second segment value S4 lies within the second range, e.g. detection signal at “0”, it is possible that the input signal S1 may bear a disease-like effect rather than a noise-like effect. Further, if the short term first segment value S3 lies within the first range and the corresponding contemporary long term second segment value S4 lies within the second range, e.g. the detection signal at “−1”, it is likely that there might be no or little noise inside the input signal S1.
Please refer to
For better understanding the operation of the reduction module 28, please refer to
In practice, the baseline tracking unit 280 may generate the estimated baseline signal BS according to at least one of the first segment values S3 and the second segment values S4, or in combination. For example, the baseline tracking unit 280 may perform interpolation over the plurality of first segment values S3 to generate the estimated baseline signal BS.
In one embodiment, please refer to
As mentioned earlier, when the characteristics of the input signal S1 changes, the value of the detection signal might change. For example, when there's less noise inside the input signal S1, either the first segment values S3 lies within the first range or the corresponding contemporary of the second segment values S4 lies within the second range. In this situation, the detection signal will bear a corresponding value. The reduction module 28 then enters into a second state. And, a second step size will be adopted by the adaptive filter 282 to update the variable weights 286_1˜286_n. In other words, the adaptive filter 282 then generates the output signal S5 by filtering the estimated baseline signal BS according to the second step size and the sample signal S2 while the reduction module 28 being in the second state.
To be noted, when there's no or little noise inside the input signal S1, the reduction module 28 might not enter into the first state nor the second state. In this case, the adaptive filter 282 might just bypass the sample signal S2 as the output signal S5.
As another embodiment, the estimated baseline signal BS may be directly subtracted from the sample signal S2 to generate the output signal S5. In this situation, the adaptive filter 282 is not needed and may be entirely replaced by an adder, resulting in a lower cost. Under certain operating conditions, experiments show that this approach also bear competitive performance in reducing the noise inside the sample signal S2.
Please refer back to
In practice, the schematic diagram of the reduction module 528 may be almost the same as that shown in
Please refer to
As an example showing operation of the calibration module 39, considering the input signal S1 is an ECG signal. Upon startup of the signal processing device 3, the plurality of first segment values S3 and the plurality of second segment values S4 are monitored by the calibration module 39. Typically, during this startup procedure, the patients may be asked to stay inactive physically, e.g. lie down on a bed, to get a resting ECG signal. In this situation, the calibration module 39 can set up initial values of the first range, the second range, the first time interval and the second time interval to the signal processing device 3.
Secondly, the calibration module 39, during a normal operation procedure, continues monitoring the plurality of first segment values S3, the plurality of second segment values S4 and the detection signal D1 so as to adjust the first range, the second range, the first time interval and/or the second time interval. During this normal operation procedure, there might be some noise or morphology variation affecting the input signal. For example, if the detection signal indicates the input signal S1 has a disease-like pattern frequently, the first time interval and the second time interval may be adjusted by the calibration module 39 to have longer period than the original settings. In this way, the plurality of first segment values and the plurality of second segment values are not likely to exceed the first range and the second range too easily. The detection module 36 will be less sensitive to the morphology variation, and the reduction module 38 may at least not tend to enter into the first state to perform noise reduction filtering. Thus, the morphology of the first signal S1 can be preserved from distortion.
On the other hand, if the input signal S1 is indicated to have a noise-like pattern, the first time interval and the second time interval may be adjusted to have shorter period than the original settings, so that the detection module 36 can be more sensitive to the noise. Therefore, the reduction module 38 might be more likely to enter into the first or second state to perform noise reduction filtering and the noise component might be more easily removed from the input signal S1.
In order to explain how the signal processing device 1 works, please refer to
To sum up, the embodiments of the invention disclose a method and a device for signal processing which can reduce the noise components of a specific signal, especially the ECG signal. One of the advantages of the approach is that users can easily and correctly analyze the waveform profile between successive heartbeats in the filtered ECG signal, and can avoid misdiagnosing the level of cardiac stability.
The descriptions illustrated supra set forth simply the preferred embodiments of the present invention; however, the characteristics of the present invention are by no means restricted thereto. All changes, alternations, or modifications conveniently considered by those skilled in the art are deemed to be encompassed within the scope of the present invention delineated by the following claims.
Number | Name | Date | Kind |
---|---|---|---|
4906941 | Kato et al. | Mar 1990 | A |
6441661 | Aoki | Aug 2002 | B1 |
6701445 | Majos | Mar 2004 | B1 |
7054402 | Muellner | May 2006 | B2 |
8542068 | Doyle et al. | Sep 2013 | B2 |
20050008112 | Takeo et al. | Jan 2005 | A1 |
20090310730 | Zhidkov et al. | Dec 2009 | A1 |
20100318151 | Saha | Dec 2010 | A1 |
20130162309 | Tamura | Jun 2013 | A1 |
20140266824 | Lowney | Sep 2014 | A1 |
Number | Date | Country |
---|---|---|
1032288 | Apr 1989 | CN |
1596825 | Mar 2005 | CN |
102215040 | Oct 2011 | CN |
103 58 397 | Jul 2005 | DE |
1 995 685 | Nov 2008 | EP |
1995685 | Nov 2008 | EP |
2 294 978 | Mar 2011 | EP |
2012129413 | Sep 2012 | WO |
Entry |
---|
Vessela Krasteva et al., Assessment of ECG frequency and morphology parameters for automatic classification of life-threatening cardiac arrhythmias, Physiological Measurement, Jun. 17, 2005, pp. 707-723, vol. 26, XP020092224, Institute of Physics Publishing, UK. |
Antti Koski, Primitive coding of structural ECG features, Pattern Recognition Letters, 1996, pp. 1215-1222, vol. 17, XP004025305, Elsevier Science B.V. |
Chris F. Zhang et al., VLSI Friendly ECG QRS Complex Detector for Body Sensor Networks, IEEE Journal on Emerging and Selected Topics in Circuits and Systems, Mar. 2012, pp. 52-59, vol. 2, No. 1, XP011441042, IEEE. |
Number | Date | Country | |
---|---|---|---|
20150061758 A1 | Mar 2015 | US |