US8498696B2 - Health monitoring device and human electric signal processing method - Google Patents
Health monitoring device and human electric signal processing method Download PDFInfo
- Publication number
- US8498696B2 US8498696B2 US12/796,634 US79663410A US8498696B2 US 8498696 B2 US8498696 B2 US 8498696B2 US 79663410 A US79663410 A US 79663410A US 8498696 B2 US8498696 B2 US 8498696B2
- Authority
- US
- United States
- Prior art keywords
- electric signal
- monitoring device
- human
- health monitoring
- processing method
- 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.)
- Expired - Fee Related, expires
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/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/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/346—Analysis of electrocardiograms
-
- 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/7253—Details of waveform analysis characterised by using transforms
- A61B5/726—Details of waveform analysis characterised by using transforms using Wavelet transforms
Definitions
- the present invention relates to health monitoring devices and human electric signal processing methods, and in particular relates to health monitoring devices and human electric signal processing methods using wavelet transform and spectral estimation.
- every heart beat period includes a P wave, a Q wave, a R wave, a S wave, a T wave and a U wave. Additionally, every wave has individual physiological representations. Therefore, whether a human being is healthy may be determined by the characteristics of the waves.
- Republic of Chinese Patent discloses a method for locating desired points within an electrocardiogram signal.
- the method applies wavelet transform and multi-scale differential operations to locate at least a desired point in at least a frequency band within the electrocardiogram signal for health determination. Additionally, after the wavelet transform is performed, coefficients thereof are further used in the multi-scale differential operation.
- the present invention provides a health monitoring device, which comprises a detecting unit and a processing unit.
- the detecting unit detects a first electric signal from a human body; and the processing unit performs the following steps of: receiving the first electric signal; performing wavelet transform on the first electric signal, to filter frequency bands thereof which include undesired noises, and then perform reverse wavelet transform to obtain a second electric signal; combining frequency bands of the first electric signal which include primary characteristics and then perform reverse wavelet transform on the combined frequency bands to obtain a third electric signal; obtaining period information of the second electric signal from the third electric signal, and then re-sampling the second electric signal according to the period information to obtain a fourth electric signal; and comparing the fourth electric signal with a plurality of electric signal patterns in a database to determine whether the fourth electric signal is an electric signal for healthiness.
- the present invention also provides a human electric signal processing method, which comprises: receiving a first electric signal of a human body; performing wavelet transform on the first electric signal to filter frequency bands thereof which include undesired noises, and then perform reverse wavelet transform to obtain a second electric signal; combining frequency bands of the first electric signal which include primary characteristics and then performing reverse wavelet transform on the combined frequency bands to obtain a third electric signal; obtaining period information of the second electric signal from the third electric signal, and then re-sampling the second electric signal according to the period information to obtain a fourth electric signal.
- FIG. 1 is a block diagram of a health monitor device according to the present invention
- FIG. 2 is a flow chart of a human electric signal processing method performed by the processing unit of FIG. 1 ;
- FIG. 3 is a schematic diagram illustrating a plurality of decomposed waveforms after performing a wavelet transform
- FIG. 4 a is a schematic diagram of a QRS waveband after combining the signals d 3 and d 4 and then performing a reverse wavelet transform on the combined signals;
- FIG. 4 b is a schematic diagram of a P, T, and U waveband after combining the signals d 6 , d 7 and d 8 and then performing a reverse wavelet transform on the combined signals.
- FIG. 1 is a block diagram of a health monitor device 10 according to the present invention.
- the health monitor device 10 is disposed on a human body part and detect the health signal of a human body.
- the health monitor device 10 may be designed to have a watch-like shape to be worn on a patient's wrist, but the invention is not limited thereto.
- the health monitor device 10 comprises a detecting unit 102 , a processing unit 104 and a wireless unit 106 .
- the detecting unit 102 detects a first electric signal of a human body, and sends the first electric signal to the processing unit 104 .
- the first electric signal may be an electrocardiogram signal.
- the processing unit 104 is coupled to the detecting unit 102 , and receives the first electric signal, determines whether the first electric signal is an electric signal for healthiness, by performing a human electric signal processing method 20 (as shown in FIG. 2 ), and then sends the result from the determination to the wireless unit 106 . After the wireless unit 106 receives the result from the processing unit 104 , the result is sent by wireless communications to an external source.
- the determination result may be sent to a remote health control center, and when the result indicates that the patient under observation is not healthy, the remote health control center would assign medical professionals to help the patient.
- various communication methods may be used in the wireless unit 106 , such as a Wi-Fi, WLAN and WiMAX communication method and others.
- FIG. 2 is a flow chart of a human electric signal processing method 20 performed by the processing unit 104 of FIG. 1 .
- step S 200 The process starts from step S 200 .
- step S 201 the processing unit 104 receives the first electric signal (electrocardiogram signal) from the detecting unit 102 .
- step S 202 the processing unit 104 performs wavelet transform on the first electric signal, to filter the frequency bands thereof which comprise undesired noises, and then performs reverse wavelet transform on to generate a second electric signal.
- the frequency bands which include undesired noises in the first electric signal are low frequency drift bands and high frequency noise bands.
- FIG. 3 is a schematic diagram illustrating a plurality of decomposed waveforms after performing a db3-based wavelet transform.
- the wavelet transform technique is well-known in the prior art and will not be discussed hereinafter. In FIG.
- the signal S is the first electric signal, and after 11 levels of wavelet transform, the signals d 1 , d 2 , d 3 , d 4 , d 5 , d 6 , d 7 , d 8 , d 9 , d 10 , d 11 and a 11 are generated.
- S d 1 +d 2 +d 3 +d 4 +d 5 +d 6 +d 7 +d 8 +d 9 +d 10 +d 11 +a 11 , wherein the signals d 1 and d 2 belong to high frequency noise bands, and the signals d 9 , d 10 , d 11 belong to low frequency drift bands. Therefore, in step S 202 , the signals may be eliminated and the de-noised second electric signal will be constructed after the wavelet transform is completed.
- step S 203 the processing unit 104 performs normalization to the second electric signal so that magnitude range of the electric signal is consistent with the magnitude range of the electric signal patterns in a database.
- the waveform, the wavelength, and the amplitude of the electric signal patterns are all important parameters to be estimated, thus, all of the signals in the database would be adjusted to be within the same scale for comparison. Therefore, in step S 203 , the magnitude of the second electric signal may be normalized so that the magnitude range of the electric signal is consistent with the magnitude range of the electric signal patterns in the database. For example, when the voltages of the electric signal patterns in the database are within ⁇ 1 ⁇ +1 mV, the voltage of the second electric signal would be normalized to be within the rage ( ⁇ 1 ⁇ +1 mV).
- step S 204 the processing unit 104 combines frequency bands of the first electric signal which include primary characteristics and then performs reverse wavelet transform on the combined frequency bands to obtain a third electric signal.
- the frequency bands which include primary characteristics further comprise a QRS waveband, a P waveband, a T waveband, and a U waveband.
- the primary characteristics of a QRS waveband are located in the signals d 3 and d 4
- the primary characteristics of the P, T and U wavebands are located in the signals d 6 , d 7 and d 8 .
- FIG. 4 a shows a QRS waveband after combining the signals d 3 and d 4 and then performing a reverse wavelet transform on the combined signals
- FIG. 4 b shows a P, T, and U waveband after combining the signals d 6 , d 7 and d 8 and then performing a reverse wavelet transform on the combined signals.
- step S 205 the processing unit 104 obtains period information of the second electric signal from the third electric signal, and then re-samples the second electric signal according to the period information to obtain a fourth electric signal.
- the processing unit 104 utilizes the QRS waveband or the P, T, and U wavebands obtained from step S 204 , to obtain a heart beat period, and then obtains the second electric signal which has a single period according to the heart beat period Then, the processing unit 104 re-samples the second electric signal with the single period to obtain a fourth electric signal.
- the number of data points in a single heart beat period are usually different from one to one, which causes the characteristics in frequency domain unable to be compared with the other ones under the same standards, thus, a step for re-sampling the data points is necessary,
- the number of re-sampled data points is 300, but the number is not limited thereto in other embodiments.
- step S 206 the processing unit 104 performs spectral estimation on the fourth electric signal to obtain a plurality of characteristic coefficients.
- spectral estimation is performed by using Autoregressive Model (AR Model) to obtain the characteristic coefficients, and the characteristic coefficients may be three coefficients of a 3-order ARM.
- AR Model Autoregressive Model
- y[n] is the fourth electric signal
- y[n ⁇ 1], y[n ⁇ 2], y[n ⁇ 3] are electric signals previously measured
- u[n] is noise
- a1, a2 and a3 are correlation coefficients.
- the characteristic coefficients obtained in step S 206 are the correlation coefficients described above.
- the 3-order ARM is taken as an example, but those skilled in the art may use other ARMs with higher orders (5-order, 6-order, etc.) to obtain more characteristic coefficients in other embodiments.
- step S 207 the processing unit 104 compares the fourth electric signal with a plurality of electric signal patterns in a database to determine whether the fourth electric signal is an electric signal for healthiness.
- Step S 207 may be performed by a Support Vector Machine (SVM).
- SVM Support Vector Machine
- the processing unit 104 may couple the fourth electric signal to the characteristic coefficients (a1, a2, and a3) in series to generate a new characteristic ⁇ right arrow over (x) ⁇ ; and then, output the new characteristic ⁇ right arrow over (x) ⁇ to the SVM.
- the SVM is performed by the x following equation (B):
- ⁇ i * is Lagrange multipliers
- ⁇ right arrow over (x) ⁇ j is generated by coupling the electric signal patterns in the database to their corresponding characteristic coefficients in series
- y i indicates that ⁇ right arrow over (x) ⁇ j is health patterns or ill health patterns. For example, when ⁇ right arrow over (x) ⁇ j is a pattern for health, y i is smaller than zero; otherwise, y i is greater than zero.
- exp ⁇ ( - 1 ⁇ 2 ⁇ ⁇ x ⁇ i - x ⁇ j ⁇ 2 ) is a Radial-basis function (RBF) Kernel function, and the SVM is performed by using RBF Kernel function to map data to a high dimension space to improve accuracy for data comparison.
- b is a correlation value.
- the SVM may compare the new characteristic ⁇ right arrow over (x) ⁇ with the numerous patterns ⁇ right arrow over (x) ⁇ j in the database. When the result of the equation (B) is greater than zero, the fourth electric signal may be determined as a signal for ill health; otherwise, when the result of the equation (B) is smaller than zero, the fourth electric signal may be determined as a signal for healthiness.
- step S 208 the steps end.
- the present invention firstly performs the wavelet transform technique to the electrocardiogram to eliminate undesired noises, normalizes and re-samples. Secondly, spectral estimation is performed to obtain characteristic coefficients and the characteristic coefficients are coupled to the processed electrocardiogram signal in series to generate new characteristics. Finally, the present invention compares the new characteristics by utilizing SVM to determine whether the fourth electric signal is an electric signal for healthiness.
- the present invention is not limited to application only to an electrocardiogram signal.
- Other electric signals such as electroencephalogram signals, electromyogram signals or blood sugar level concentration signals may also be applied.
- the normalization may also be performed to the first electric signal in step S 201 or to the fourth electric signal in step S 205 .
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Cardiology (AREA)
- Surgery (AREA)
- Medical Informatics (AREA)
- Veterinary Medicine (AREA)
- Physics & Mathematics (AREA)
- Public Health (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Animal Behavior & Ethology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Physiology (AREA)
- Signal Processing (AREA)
- Artificial Intelligence (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
Abstract
Description
y[n]=a1*y[n−1]+a2*y[n−2]+a3*y[n−3]+u[n] (A)
is a Radial-basis function (RBF) Kernel function, and the SVM is performed by using RBF Kernel function to map data to a high dimension space to improve accuracy for data comparison. Meanwhile, b is a correlation value. By using the function (B), the SVM may compare the new characteristic {right arrow over (x)} with the numerous patterns {right arrow over (x)}j in the database. When the result of the equation (B) is greater than zero, the fourth electric signal may be determined as a signal for ill health; otherwise, when the result of the equation (B) is smaller than zero, the fourth electric signal may be determined as a signal for healthiness.
Claims (20)
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| TW98119888A | 2009-06-15 | ||
| TWTW098119888 | 2009-06-15 | ||
| TW098119888A TWI365062B (en) | 2009-06-15 | 2009-06-15 | Health monitoring device and human electric signal processing method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| US20100317934A1 US20100317934A1 (en) | 2010-12-16 |
| US8498696B2 true US8498696B2 (en) | 2013-07-30 |
Family
ID=43307006
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| US12/796,634 Expired - Fee Related US8498696B2 (en) | 2009-06-15 | 2010-06-08 | Health monitoring device and human electric signal processing method |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US8498696B2 (en) |
| TW (1) | TWI365062B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105232032A (en) * | 2015-11-05 | 2016-01-13 | 福州大学 | Remote electrocardiograph monitoring and early warning system and method based on wavelet analysis |
| KR20200120494A (en) | 2019-04-11 | 2020-10-21 | 삼성전자주식회사 | A drift, noise, and motion artifact correction method for photoplethysmogram(ppg) signal and apparatus thereof |
| US11607178B2 (en) | 2019-04-11 | 2023-03-21 | Samsung Electronics Co., Ltd. | Drift, noise, and motion artifact correction method for photoplethysmogram (PPG) signals and system thereof |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103750835A (en) * | 2013-11-12 | 2014-04-30 | 天津工业大学 | Electrocardiosignal characteristic detection algorithm |
| CN104458173B (en) * | 2014-11-27 | 2017-04-12 | 广东电网有限责任公司中山供电局 | Steel framework structure mutational damage recognition method and system |
| CN107341769A (en) * | 2016-05-03 | 2017-11-10 | 中国科学院微电子研究所 | Electrocardiosignal denoising method and system |
| US9642577B1 (en) * | 2016-08-16 | 2017-05-09 | American Reliance, Inc. | Methods and systems for disease analysis based on transformations of diagnostic signals |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| TW200822904A (en) | 2006-11-30 | 2008-06-01 | Taipei Veterans General Hospital Vac | Method for orientating characteristic points of electrocardio-signal and storage media and electronic apparatus and electrocardio-signal analysis system used by the method |
| US7751873B2 (en) * | 2006-11-08 | 2010-07-06 | Biotronik Crm Patent Ag | Wavelet based feature extraction and dimension reduction for the classification of human cardiac electrogram depolarization waveforms |
| US7809433B2 (en) * | 2005-08-09 | 2010-10-05 | Adidas Ag | Method and system for limiting interference in electroencephalographic signals |
-
2009
- 2009-06-15 TW TW098119888A patent/TWI365062B/en active
-
2010
- 2010-06-08 US US12/796,634 patent/US8498696B2/en not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7809433B2 (en) * | 2005-08-09 | 2010-10-05 | Adidas Ag | Method and system for limiting interference in electroencephalographic signals |
| US7751873B2 (en) * | 2006-11-08 | 2010-07-06 | Biotronik Crm Patent Ag | Wavelet based feature extraction and dimension reduction for the classification of human cardiac electrogram depolarization waveforms |
| TW200822904A (en) | 2006-11-30 | 2008-06-01 | Taipei Veterans General Hospital Vac | Method for orientating characteristic points of electrocardio-signal and storage media and electronic apparatus and electrocardio-signal analysis system used by the method |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105232032A (en) * | 2015-11-05 | 2016-01-13 | 福州大学 | Remote electrocardiograph monitoring and early warning system and method based on wavelet analysis |
| KR20200120494A (en) | 2019-04-11 | 2020-10-21 | 삼성전자주식회사 | A drift, noise, and motion artifact correction method for photoplethysmogram(ppg) signal and apparatus thereof |
| US11607178B2 (en) | 2019-04-11 | 2023-03-21 | Samsung Electronics Co., Ltd. | Drift, noise, and motion artifact correction method for photoplethysmogram (PPG) signals and system thereof |
Also Published As
| Publication number | Publication date |
|---|---|
| TW201043195A (en) | 2010-12-16 |
| TWI365062B (en) | 2012-06-01 |
| US20100317934A1 (en) | 2010-12-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US8498696B2 (en) | Health monitoring device and human electric signal processing method | |
| US6389308B1 (en) | System and device for multi-scale analysis and representation of electrocardiographic data | |
| US7343197B2 (en) | Multi-scale analysis and representation of physiological and health data | |
| US9760536B2 (en) | Method and apparatus for low complexity spectral analysis of bio-signals | |
| EP3366203B1 (en) | Method and system for cuffless blood pressure estimation using photoplethysmogram features and pulse transit time | |
| İşler et al. | Combining classical HRV indices with wavelet entropy measures improves to performance in diagnosing congestive heart failure | |
| Chawla | A comparative analysis of principal component and independent component techniques for electrocardiograms | |
| Li et al. | Design of a continuous blood pressure measurement system based on pulse wave and ECG signals | |
| US10912464B2 (en) | System and method for monitoring vascular system health | |
| EP3003140B1 (en) | Electrocardiogram analysis | |
| US20190175042A1 (en) | Method and apparatus for determining blood pressure | |
| US11576601B2 (en) | Artifact identification in EEG measurements | |
| JP6434432B2 (en) | Improved ECG calculation method for use in generating 12-lead ECG measurements from a device having less than 10 electrodes | |
| CN112107303A (en) | Apparatus and method for calibrating biological information estimation model, and biological information estimation apparatus | |
| JP2021115429A (en) | Signal measuring system, measurement signal processing device and program | |
| Keshtkar et al. | Discriminant analysis between myocardial infarction patients and healthy subjects using Wavelet Transformed signal averaged electrocardiogram and probabilistic neural network | |
| Degirmenci et al. | Synchrosqueezing transform in biomedical applications: A mini review | |
| CN115721285B (en) | Device and method for extracting heartbeat data based on wireless radar signals | |
| Rajeshkumar et al. | UTO-LAB model: USRP based touchless lung anomaly detection model with optimized machine learning classifier | |
| WO2018018570A1 (en) | Device and method for measuring electrocardiogram | |
| JP6706996B2 (en) | Biological signal processing device, abnormality determination method and program | |
| KR101464066B1 (en) | apparatus and method for measuring living body information | |
| CN118749933A (en) | A non-contact heart rate detection method based on ultra-wideband bio-radar | |
| Orhan et al. | Heart rate detection on single-arm ECG by using dual-median approach | |
| Paul et al. | Mental stress detection using multimodal characterization of PPG signal for personal healthcare applications |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| AS | Assignment |
Owner name: NATIONAL TAIWAN UNIVERSITY, TAIWAN Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:YEN, JIA-YUSH;YEN, FAN-CHE;LEE, CHIH-KUNG;AND OTHERS;SIGNING DATES FROM 20100218 TO 20100303;REEL/FRAME:024519/0652 |
|
| STCF | Information on status: patent grant |
Free format text: PATENTED CASE |
|
| FPAY | Fee payment |
Year of fee payment: 4 |
|
| MAFP | Maintenance fee payment |
Free format text: PAYMENT OF MAINTENANCE FEE, 8TH YR, SMALL ENTITY (ORIGINAL EVENT CODE: M2552); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY Year of fee payment: 8 |
|
| FEPP | Fee payment procedure |
Free format text: MAINTENANCE FEE REMINDER MAILED (ORIGINAL EVENT CODE: REM.); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY |
|
| LAPS | Lapse for failure to pay maintenance fees |
Free format text: PATENT EXPIRED FOR FAILURE TO PAY MAINTENANCE FEES (ORIGINAL EVENT CODE: EXP.); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY |
|
| STCH | Information on status: patent discontinuation |
Free format text: PATENT EXPIRED DUE TO NONPAYMENT OF MAINTENANCE FEES UNDER 37 CFR 1.362 |
|
| FP | Lapsed due to failure to pay maintenance fee |
Effective date: 20250730 |