Apparatus and method for removing noise from biosignals

Disclosed are an apparatus and a method for removing noise from biosignals. An embodiment of the invention provides an apparatus for removing noise from biosignals measured for a particular period of time that includes: a similarity calculating unit configured to calculate a similarity between any one biosignal from among the biosignals and other biosignals; and a noise removal unit configured to remove noise from the any one biosignal by using the similarity. Aspects of the invention provide the advantage of effectively removing noise from biosignals without distorting the waveforms of the biosignals.

CROSS-REFERENCE TO RELATED APPLICATIONS

This application claims the benefit of Korean Patent Application No. 10-2012-0025515 filed on Mar. 13, 2012, the disclosure of which is incorporated herein by reference in its entirety.

TECHNICAL FIELD

The embodiments of the present invention relate to an apparatus and a method for removing noise from biosignals, more particularly to an apparatus and a method for removing noise from biosignals based on similarity between the biosignals.

BACKGROUND ART

Electrocardiogram (ECG) signals are electrical signals that originate from the heart and are important in diagnosing and monitoring heart conditions. With the advent of e-health and u-health programs, etc., in recent times, there is an interest in technology for measuring ECG signals in a mobile environment.

A major problem for measuring ECG signals in a mobile environment is low signal-to-noise ratio (hereinafter, “SNR”). In a typical hospital, ECG signals are measured for a stable posture in a stable environment using high-performance equipment, but in a mobile environment, the ECG measurement may be performed for a subject that is not always in a stable posture, and the equipment used may be a compact device that provides a relatively low performance.

For ECG signals obtained in such an environment, it may not be possible to sufficiently remove noise using conventional, typical noise removal methods.

That is, conventional ECG noise removal methods can be divided mainly into methods that use general high/band/low pass filters and methods that use wavelets. These conventional methods are effective when applied to signals having a somewhat high SNR, but if they are applied to ECG signals having a low SNR and noise is removed excessively, the shapes of waveforms such as the QRS-complex, P-wave, T-wave, etc., essential in diagnosing a human body, may be greatly altered.

SUMMARY

To resolve the problems of the related art, an aspect of the invention provides an apparatus and a method for effectively removing noise from biosignals measured at a mobile device.

Other objectives of the invention can be derived by those skilled in the art from the embodiments of the invention described below.

To achieve the objectives above, an embodiment of the invention provides an apparatus for removing noise from biosignals measured for a particular period of time. The apparatus includes: a similarity calculating unit configured to calculate a similarity between any one biosignal from among the biosignals and other biosignals; and a noise removal unit configured to remove noise from the any one biosignal by using the similarity.

The noise removal unit can calculate weights for the other biosignals respectively by using the similarity and can remove noise from the any one biosignal based on a weighted average of the biosignals using the weights.

The biosignal can include an ECG (electrocardiogram) signal, and the ECG signal can include a P-wave, a QRS-complex, and a T-wave waveform.

The apparatus can further include a selection unit configured to select a particular target signal for the noise removal from among signals included in the any one biosignal.

The similarity calculating unit can include: an alignment unit configured to align the biosignals with respect to an R peak of the QRS complex waveform included in the biosignals; a search window setting unit configured to set a range of a search window to include the target signal and signals within a particular range centering around the target signal from among signals included in the aligned biosignals; a unit window setting unit configured to set a size of a unit window for the calculating of similarity; and a calculating unit configured to calculate a similarity between signals included in a unit window centering around the target signal and signals included in unit windows centering around signals included in the search window.

The alignment unit can separate the biosignals based on a preset area centering around the R peak and can align the separated biosignals with respect to the R peak.

The calculating unit can calculate the similarity by using a mean absolute deviation between the signals included in a unit window centering around the target signal and the signals included in unit windows centering around signals included in the search window.

The search window setting unit can set the range of the search window differently according to a type of waveform to which the selected target signal belongs, and the unit window setting unit can set the size of the unit window differently according to a type of waveform to which the selected target signal belongs.

The noise removal unit can include: a weight calculating unit configured to calculate weights for signals included in the search window respectively by using the similarity; and a weighted averaging unit configured to perform a weighted averaging for signals included in the search window by using the weights.

Another embodiment of the invention provides a method for removing noise from biosignals measured for a particular period of time. The method includes: calculating a similarity between any one biosignal from among the biosignals and other biosignals; and removing noise from the any one biosignal by using the similarity.

The biosignal can be an ECG (electrocardiogram) signal including a P-wave, a QRS-complex, and a T-wave waveform.

Calculating the similarity can include: aligning the biosignals with respect to an R peak of the QRS complex waveform included in the biosignals; selecting a particular target signal for the noise removal from among signals included in the any one biosignal; setting a range of a search window to include the target signal and signals within a particular range centering around the target signal from among signals included in the aligned biosignals; setting a size of a unit window for the calculating of similarity; and calculating a similarity between signals included in a unit window centering around the target signal and signals included in unit windows centering around signals included in the search window.

Aspects of the invention provide the advantage of effectively removing noise from biosignals without distorting the waveforms of the biosignals.

DETAILED DESCRIPTION

Certain embodiments of the present invention will be described below in more detail with reference to the accompanying drawings.

Biosignals are signals created in the human body and can be used to check a person's health. Examples of biosignals include pulse signals, ECG (electrocardiogram) signals, etc.

Biosignals, which are created in the body repeatedly in accordance with the beating of the heart, tend to have similar waveforms. A principle of the invention is remove noise from biosignals by using this similarity in biosignals.

FIG. 1is a block diagram illustrating the composition of an apparatus for removing noise from biosignals according to an embodiment of the invention.

Referring toFIG. 1, the noise removal apparatus100can include a biosignal measuring unit101, a selection unit103, a similarity calculating unit105, and a noise removal unit107.

The biosignal measuring unit101may measure the biosignals, which are created repeatedly in the body and which have similar waveforms, for a particular period of time. For example, the biosignal measuring unit101can receive signals detected by sensors attached to a part of the body, in order to measure the biosignals. Of course, various other embodiments can be employed by which the user directly uses a mobile device for measuring biosignals. As this is not the key feature of the invention, further descriptions on this matter will be omitted.

The selection unit103may select any one biosignal from which to remove noise, from among the biosignals measured by the biosignal measuring unit101for a particular period of time.

The similarity calculating unit105may calculate a similarity between the any one biosignal, selected by the selection unit103from among the biosignals measured for the particular period of time, and the other biosignals.

Then, the noise removal unit107may remove noise from the any one biosignal, using the similarity between the any one biosignal and other biosignals.

For example, the noise removal unit107can calculate weights for the other biosignals, respectively, using the similarity, and can remove noise from the any one biosignal based on the weighted average of the biosignals using the weights.

Thus, by repeating the noise removal process for all of the biosignals measured for a particular period of time, the noise removal apparatus100can remove noise from the biosignals.

For convenience, it will be assumed in the descriptions below that the biosignals are ECG signals, but it would be obvious to those skilled in the art that the biosignals are not limited to ECG signals.

FIG. 2illustrates an example of an ECG signal.

Referring toFIG. 2, an ECG signal is a measured value of the action current created by the heart muscles in accordance with the beating of the heart. The ECG signal includes a P-wave, a QRS-complex, and a T-wave waveform; in the body, the ECG signals including the three types of waveforms above are created in repetition.

The QRS-complex waveform, which includes an R peak, may be measured in a short amount of time, while the P-wave and the T-wave maybe measured with relatively smaller current values in longer periods of time compared to the QRS-complex waveform.

According to an embodiment of the invention, the noise removal apparatus100can remove noise considering the waveform of the ECG signal, and to this end, the selection unit103can select any one biosignal from among the biosignals measured for a particular period of time and can select a particular target signal401for noise removal from among the signals included in the any one biosignal.

Also, the noise removal apparatus100can perform the noise removal considering which waveform the particular target signal401selected by the selection unit103belongs to. This aspect will be described later in further detail.

FIG. 3is a block diagram illustrating the detailed composition of a similarity calculating unit according to an embodiment of the invention.

Referring toFIG. 3, the similarity calculating unit105can include an alignment unit301, a search window setting unit303, a unit window setting unit305, and a calculating unit307.

The alignment unit301may align the biosignals with respect to the R peak of the QRS-complex waveform included in each biosignal. To this end, the alignment unit301can separate the biosignals measured for a particular period of time based on a preset area centering around the R peak, and can align the separated biosignals with respect to the R peak. Here, the preset area may refer to an area including the P-wave and the T-wave waveforms centering around the QRS-complex waveform.

FIG. 4illustrates an example of biosignals aligned according to an embodiment of the invention.

InFIG. 4, it can be seen that the biosignals measured by the alignment unit301for a particular period of time have been aligned with respect to the R peak, and the selection unit103can select a particular target signal401, from which to remove noise, from among the aligned biosignals.

According to an embodiment of the invention, the selection unit103can select a particular target signal401from signals that belong to any one of the P-wave, QRS-complex, and T-wave waveforms, from among the signals included in the biosignals.

Referring again toFIG. 3, the search window setting unit303may set the range of a search window, by which to calculate the similarities for the aligned biosignals.

To be more specific, the search window setting unit303may set the range of the search window to include the target signal401and the signals that lie within a particular range centering around the target signal401, from among the signals included in the aligned biosignals.

FIG. 5illustrates an example of setting a search window according to an embodiment of the invention.

FIG. 5illustrates portions of the aligned biosignals, and the search window setting unit303may set the range of the search window500for calculating the similarities between the biosignals aligned as shown inFIG. 5.

According to an embodiment of the invention, the search window setting unit303can set the range of the search window, i.e. the position and size of the search window500, according to the type of waveform to which the target signal401belongs. For example, the search window setting unit303can select the range of a first search window501to include the P-wave waveform region from among the aligned biosignals if the target signal401belongs to the P-wave waveform, and can select the range of a second search window503to include the QRS-complex waveform region if the target signal401belongs to the QRS-complex waveform. Also, if the target signal401belongs to the T-wave waveform, the range of a third search window505can be set to include the T-wave waveform region.

Since the QRS-complex waveform is measured for a relatively shorter amount of time than for the P-wave and T-wave waveforms, the size of the second search window503along direction of the lateral axis can be set smaller than the lateral direction sizes of the search windows501,505for the other waveforms. That is, the search window setting unit303can set the range of the search window500according to the type of waveform, considering the length of time during which the waveform is measured.

Referring again toFIG. 3, the unit window setting unit303may set the size of a unit window for calculating similarity.

Then, the calculating unit307may calculate the similarities between the signals included in the unit window centering around the target signal401and the signals included in unit windows centering around signals included in the search window.

FIG. 6illustrates an example of calculating similarity by using a search window and unit windows according to an embodiment of the invention.

Where drawing (a) ofFIG. 6is a magnified view of the region of the third search window505, the unit window setting unit303may set the size of the unit window610for calculating similarity. According to an embodiment of the invention, the size of the unit window610can also be set differently according to the type of waveform, as is the case for the size of the search window500.

The calculating unit307may calculate the similarities between the signals included in the unit window (hereinafter referred to as the “first unit window611”) centering around the target signal401and the signals included in the unit windows (hereinafter referred to as “second windows613”) respectively centering around signals included in the search window610.

According to an embodiment of the invention, the target signal401from which to remove noise can be denoted IR, and the signals included in a search window500of a particular section centering around the target signal can be denoted I1. . . In. . . IN. Also, the k number of signals included in the first unit window611centering around the target signal IRcan be denoted BR=[BR1, BR2, . . . , BRk], while the k number of signals included in the second unit windows613centering around the remaining signals I1˜IN, respectively, can be denoted Bn=[Bn1, Bn2, . . . , Bnk], n=1 . . . N.

In this case, the calculating unit307can perform the similarity calculations between the target signal IRtargeted for noise removal and the remaining signals I1. . . In. . . INby using the sum of absolute deviation (SAD) between the first unit window611and the second unit window613for each of the remaining signals. The similarity calculation can also be performed by using MSE, correlation, etc., besides SAD.

For example, the similarity between IRand Incan be expressed as Equation 1 below.
SADn=Σ|BR−Bn|=Σi=1k|BiR−Bin|  [Equation 1]

Here, SADnrepresents the similarity between IRand In, where a smaller value means greater similarity.

Drawing (b) ofFIG. 6shows weights for the signals calculated by using the similarities; a description of the weights will be provided below with reference toFIG. 7. A darker color (i.e. nearing 0) represents low similarity, while a lighter color (i.e. nearing 1) represents high similarity.

FIG. 7is a block diagram illustrating the detailed composition of a noise removal unit according to an embodiment of the invention.

Referring toFIG. 7, the noise removal unit107can include a weight calculating unit701and a weighted averaging unit703.

Based on the similarities between the target signal401and the remaining signals as calculated by Equation 1, the weight calculating unit701may calculate the weights used for the weighted averaging, such as by Equation 2 shown below.
wn=e−SADn/β[Equation 2]

Here, wnis the weight for signal In, and β is a constant for adjusting the degree of noise removal.

The weight calculating unit701may use Equation 2 to calculate weights for all of the signals I1. . . In. . . INincluded in the search window500.

By using the weights, the weighted averaging unit703may perform a weighted averaging for all signals included in the search window500, and thereby remove noise from the target signal401. The weighted averaging can be performed using Equation 3 shown below.

Here, I′Ris the target signal with noise removed, and w is the sum of weights for all signals included in the search window500.

Thus, the noise removal apparatus100can select all of the signals in the search window as target signals and remove noise for all signals.

FIG. 8is a flowchart illustrating the overall flow of a method for removing noise according to an embodiment of the invention.

The operation at each step is described below in more detail with reference toFIG. 8.

First, in step S800, the alignment unit301may separate the biosignals, measured at the biosignal measuring unit101for a particular period of time, into certain areas centering around the R peaks, and may align the separated biosignals with respect to the R peaks.

In step S805, the selection unit103may select a particular target signal, for which to perform noise removal, from among the signals included in any one biosignal.

Then, in step S810, the search window setting unit303may set the range of a search window. To be more specific, the search window setting unit303may set the range of the search window to include the target signal and signals within a particular range centering around the target signal, from among the signals included in the aligned biosignals.

In step S815, the unit window setting unit305may set the size of a unit window by which to perform similarity calculations.

In step S820, the calculating unit307may calculate the similarities between the signals included in the unit window centering around the target signal and the signals included in the unit windows centering around the signals included in the search window.

Then, in step S825, the weight calculating unit701may calculate the weights for signals included in the search window by using the similarities calculated by the calculating unit307.

Lastly, in step S830, the weighted averaging unit703may remove noise from the target signal by a weighted averaging of the signals included in the search window using the weights calculated by the weight calculating unit701.

FIG. 9throughFIG. 11show graphs representing the results of removing noise from ECG signals generated by simulations, using an existing technique and a technique based on an embodiment of the invention.

InFIG. 9, graph (a) represents multiple ECG signals generated by a simulation, graph (b) shows the signals obtained by passing the ECG signals through a band pass filter to remove noise, graph (c) shows the signals obtained by passing the ECG signals through a bilateral filter to remove noise, and graph (d) shows the signals obtained by subjecting the ECG signals to a noise removal method according to an embodiment of the invention.

InFIG. 9, comparing the signals of graph (a) with the signals of graphs (b) through (d), it can be seen that noise has been removed from the signals most effectively in graph (d), where a noise removal method according to an embodiment of the invention has been applied.

Referring toFIG. 10, it can be seen that the signal in graph (c), where the noise removal method according to an embodiment of the invention has been applied, has the smoothest curve compared with the other signals and has the noise removed with the least distortion to the P-wave, QRS-complex, and T-wave waveforms of the ECG signal.

Graph (a) ofFIG. 11shows an average error signal of the signal shown in graph (a) ofFIG. 10, graph (b) ofFIG. 11shows the error rate of signal shown in graph (b) ofFIG. 10, and graph (c) ofFIG. 11shows the error rate of the signal shown in graph (c) ofFIG. 10.

In the graphs ofFIG. 11, values closer to 0 represent smaller error. In graph (a) ofFIG. 11, it can be seen that the error is uniformly distributed over all areas, and in graph (b) ofFIG. 11, it can be seen that the error is especially large in the area of the QRS-complex waveform.

However, in graph (c) ofFIG. 11, where the noise removal method according to an embodiment of the invention has been applied, it can be seen that there is smaller error compared with the other graphs.

FIG. 12shows graphs representing the results of removing noise from actual measured ECG signals, using an existing technique and a technique based on an embodiment of the invention.

Graph (a) ofFIG. 12shows multiple actually measured ECG signals, graph (b) ofFIG. 12shows the signals with noise removed by passing the ECG signals through a band pass filter, graph (c) ofFIG. 12shows the signals with noise removed by passing the ECG signals through a bilateral pass filter, and graph (d) ofFIG. 12shows the signals with noise removed by subjecting the ECG signals to a noise removal method according to an embodiment of the invention.

Referring toFIG. 12, it can be seen that in graph (d), where a noise removal method according to an embodiment of the invention has been applied, the noise has been removed from the signals most effectively with the least amount of distortion to the waveforms of the ECG signals, compared with the other signals.

The embodiments of the present invention can be implemented in the form of program instructions that may be performed using various computer means and can be recorded in a computer-readable medium. Such a computer-readable medium can include program instructions, data files, data structures, etc., alone or in combination. The program instructions recorded on the medium can be designed and configured specifically for the present invention or can be a type of medium known to and used by the skilled person in the field of computer software. Examples of a computer-readable medium may include magnetic media such as hard disks, floppy disks, magnetic tapes, etc., optical media such as CD-ROM's, DVD's, etc., magneto-optical media such as floptical disks, etc., and hardware devices such as ROM, RAM, flash memory, etc. Examples of the program of instructions may include not only machine language codes produced by a compiler but also high-level language codes that can be executed by a computer through the use of an interpreter, etc. The hardware mentioned above can be made to operate as one or more software modules that perform the actions of the embodiments of the invention, and vice versa.

While the present invention has been described above using particular examples, including specific elements, by way of limited embodiments and drawings, it is to be appreciated that these are provided merely to aid the overall understanding of the present invention, the present invention is not to be limited to the embodiments above, and various modifications and alterations can be made from the disclosures above by a person having ordinary skill in the technical field to which the present invention pertains. Therefore, the spirit of the present invention must not be limited to the embodiments described herein, and the scope of the present invention must be regarded as encompassing not only the claims set forth below, but also their equivalents and variations.