CN113616203A - A method for detecting underwater blood oxygen of divers based on wavelet filtering algorithm - Google Patents
A method for detecting underwater blood oxygen of divers based on wavelet filtering algorithm Download PDFInfo
- Publication number
- CN113616203A CN113616203A CN202111025637.2A CN202111025637A CN113616203A CN 113616203 A CN113616203 A CN 113616203A CN 202111025637 A CN202111025637 A CN 202111025637A CN 113616203 A CN113616203 A CN 113616203A
- Authority
- CN
- China
- Prior art keywords
- wavelet
- signal
- decomposition
- blood oxygen
- filtering
- 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.)
- Pending
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/145—Measuring characteristics of blood in vivo, e.g. gas concentration or pH-value ; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid or cerebral tissue
- A61B5/1455—Measuring characteristics of blood in vivo, e.g. gas concentration or pH-value ; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid or cerebral tissue using optical sensors, e.g. spectral photometrical oximeters
- A61B5/14551—Measuring characteristics of blood in vivo, e.g. gas concentration or pH-value ; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid or cerebral tissue using optical sensors, e.g. spectral photometrical oximeters for measuring blood gases
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7253—Details of waveform analysis characterised by using transforms
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Public Health (AREA)
- Medical Informatics (AREA)
- Veterinary Medicine (AREA)
- Signal Processing (AREA)
- General Health & Medical Sciences (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Physiology (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Optics & Photonics (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
The invention discloses a method for detecting underwater blood oxygen of a diver based on a wavelet filtering algorithm, which relates to the technical field of underwater, and comprises the steps of obtaining filtering parameters by utilizing sample photoelectric volume pulse waves when the diver is in a plurality of different states under the same oxygen supply, respectively carrying out wavelet decomposition on an original signal by utilizing a preset wavelet base according to the decomposition layer number in the obtained filtering parameters to obtain a plurality of signal components, setting coefficients of interference layers except for effective signal layers in all the signal components corresponding to the obtained filtering parameters to zero, carrying out wavelet reconstruction to obtain corresponding reconstruction signals, and then detecting to obtain the blood oxygen saturation; the wavelet transform filtering algorithm adopted by the method has the localization characteristics in the time domain and the frequency domain, has a good filtering effect on the photoplethysmogram which is interfered by underwater strong motion artifacts and belongs to weak signals, and therefore the blood oxygen saturation can be extracted more accurately.
Description
Technical Field
The invention relates to the field of underwater technology, in particular to a diver underwater blood oxygen detection method based on a wavelet filtering algorithm.
Background
People can perform underwater exploration, salvage, rescue or entertainment activities for a long time, which provides challenges for underwater safety of divers. In the past half century, experts or scholars research and discuss underwater safety of divers, and blood oxygen saturation refers to the percentage of oxyhemoglobin in arterial blood in the capacity of all combinable hemoglobin, and is an important index for people to directly judge life safety in life and movement, so that monitoring the life safety of divers by collecting blood oxygen signals is a mainstream method at present. Relevant data show that, at present, there are two main ways of collecting blood oxygen signals: transmissive optical path detection and reflective optical path detection. The difficulty of detecting the tail end of a human body by using a transmission mode under water is high, signals are more easily influenced, and from practical application, the reflection type optical path detection is not limited by a detection part and is more suitable for detecting the blood oxygen signals under water.
However, the signal amplitude detected by the reflection-type optical path is small, and the underwater operation can hardly be in a static state, so that the signal is easy to introduce motion artifact interference except the interference of the external environment. Because the motion artifacts generated by different motion states are different, the strong motion artifact mixed in the signal is more difficult to remove particularly during severe motion, the accuracy of blood oxygen saturation calculation is seriously influenced, the removal of the motion artifact for collecting vital sign signals underwater is always a difficult point of life safety monitoring of current underwater operation, and great challenge is brought to accurate real-time blood oxygen saturation calculation under different motion states.
Disclosure of Invention
The invention provides a diver underwater blood oxygen detection method based on a wavelet filtering algorithm aiming at the problems and the technical requirements, and the technical scheme of the invention is as follows:
a diver underwater blood oxygen detection method based on wavelet filtering algorithm comprises the following steps:
collecting sample photoplethysmography pulse waves of a diver in a plurality of different states under the same oxygen supply, respectively performing wavelet decomposition by using a preset wavelet basis, and performing data analysis on wavelet decomposition results of the plurality of different states to determine filtering parameters, wherein the filtering parameters comprise the decomposition layer number and an effective signal layer;
collecting a photoelectric volume pulse wave to be detected of a diver as an original signal;
respectively carrying out wavelet decomposition on the original signal by utilizing a preset wavelet basis according to the decomposition layer number in the filtering parameters to obtain a plurality of signal components;
setting coefficients corresponding to interference layers except the effective signal layer in the filtering parameters in all the signal components to zero, and performing wavelet reconstruction to obtain corresponding reconstructed signals;
and extracting characteristic points of the reconstructed signal and detecting to obtain the blood oxygen saturation.
The further technical scheme is that the data analysis is carried out on wavelet decomposition results in a plurality of different states to determine filtering parameters, and the method comprises the following steps:
and performing wavelet decomposition on each sample photoplethysmogram pulse wave by using a preset wavelet base according to the current filtering parameters and reconstructing to obtain corresponding reconstruction signals, adjusting the decomposition layer number and/or effective signal layer until the reconstruction errors between all sample photoplethysmogram pulse waves and the respective corresponding reconstruction signals are within a preset error range, and determining to obtain the filtering parameters.
The further technical scheme is that a db5 wavelet is used as a preset wavelet base, and the number of decomposition layers included in the filtering parameters is 5.
The further technical scheme is that the effective signal layer included by the filtering parameters is a third layer signal component and a fourth layer signal component, and the interference layer includes a first layer signal component, a second layer signal component and a fifth layer signal component.
The further technical scheme is that photoelectric volume pulse waves of a diver under two different wavelengths are collected, and one path of photoelectric volume pulse waves with higher amplitude is selected as an original signal.
The beneficial technical effects of the invention are as follows:
the application discloses a diver underwater blood oxygen detection method based on a wavelet filtering algorithm, the method adopts the wavelet transformation filtering algorithm to filter photoplethysmographic pulse waves, the wavelet transformation filtering algorithm has localization characteristics in a time domain and a frequency domain, the window size is fixed but the shape is variable, the method has self-adaptability to signals, the defects of the traditional filtering algorithm are overcome, and the method has a better filtering effect on the photoplethysmographic pulse waves which are interfered by underwater strong motion artifacts and belong to weak signals, so that the blood oxygen saturation can be extracted more accurately.
Drawings
Fig. 1 is a schematic flow chart of the underwater blood oxygen detection method for divers disclosed in the present application.
Fig. 2 is a comparison diagram of waveforms of collected photoplethysmographic waves of a diver at two different wavelengths.
Fig. 3 is a schematic diagram of a waveform reconstructed from the photoplethysmogram in fig. 2 by the method of the present application.
Detailed Description
The following further describes the embodiments of the present invention with reference to the drawings.
The application discloses a method for underwater blood oxygen detection of divers based on wavelet filtering algorithm, please refer to the flow chart shown in fig. 1, the method includes:
1. collecting sample photoplethysmography pulse waves of a diver in a plurality of different states under the same oxygen supply, respectively performing wavelet decomposition by using a preset wavelet basis, and performing data analysis on wavelet decomposition results of the plurality of different states to determine filtering parameters, wherein the filtering parameters comprise the decomposition layer number and an effective signal layer.
The signal acquisition end of the signal acquisition equipment is fixed at the arm of a diver to acquire photoplethysmography, and different states of the diver at least comprise a static state and a working state.
Performing wavelet decomposition on each sample photoplethysmogram pulse wave by using a preset wavelet base according to current filtering parameters and reconstructing to obtain corresponding reconstruction signals, adjusting the number of decomposition layers and/or effective signal layers until reconstruction errors between all sample photoplethysmogram pulse waves and the corresponding reconstruction signals are within a preset error range, taking the current filtering parameters as finally determined filtering parameters, and calculating the reconstruction errors by adopting the existing algorithm, which is not repeated in the application.
According to the method, a db5 wavelet is used as a preset wavelet base, and the number of decomposition layers included in finally selected filtering parameters is 5 according to effect comparison and debugging. The filter parameters comprise effective signal layers of a third layer signal component and a fourth layer signal component, whereby the interference layer comprises a first layer signal component, a second layer signal component and a fifth layer signal component. The first and second layer signal components are mainly power frequency interference and environment white noise, and the five layer signal components are mainly baseline drift.
2. Collecting the photoelectric volume pulse waves to be detected of the diver as original signals, actually collecting the photoelectric volume pulse waves of the diver under two different wavelengths, and selecting one path of the photoelectric volume pulse waves with higher amplitude as the original signals. For example, in the present application, it can be found through comparison that the change frequencies of the photoplethysmography of the two wavelengths are the same and the positions of the peak-valley points are substantially the same, so that the photoplethysmography of the infrared light with a higher amplitude is selected as the original signal, for example, fig. 2 takes the number of sampling points as 200.
3. And respectively carrying out wavelet decomposition on the original signal by utilizing a preset wavelet basis according to the decomposition layer number in the filtering parameters to obtain a plurality of signal components. I.e., in this application, a 5-level wavelet decomposition of the original signal is performed using the db5 wavelet as the predetermined wavelet basis.
4. And setting coefficients corresponding to interference layers except the effective signal layer in the filtering parameters in all the signal components to zero, and performing wavelet reconstruction to obtain corresponding reconstructed signals. In the application, the coefficients of the first layer, the second layer and the layer are set to zero and reconstructed, so that the reconstructed signal can be obtained as shown in fig. 3, and the peak and trough positions of the reconstructed signal are clear, so that the blood oxygen signal characteristics are better saved, and the characteristic value is more convenient to obtain.
5. Carry out the characteristic point extraction and detect and obtain oxyhemoglobin saturation to the reconstruction signal, mainly through carrying out data peak valley point detection to the reconstruction signal and draw the characteristic point, carry out the specific algorithm that the characteristic point extracted and finally obtained oxyhemoglobin saturation based on the reconstruction signal can refer to current algorithm, and this application is no longer repeated.
What has been described above is only a preferred embodiment of the present application, and the present invention is not limited to the above embodiment. It is to be understood that other modifications and variations directly derivable or suggested by those skilled in the art without departing from the spirit and concept of the present invention are to be considered as included within the scope of the present invention.
Claims (5)
1. A diver underwater blood oxygen detection method based on wavelet filtering algorithm is characterized by comprising the following steps:
collecting sample photoplethysmography pulse waves of a diver in a plurality of different states under the same oxygen supply, respectively performing wavelet decomposition by using a preset wavelet basis, and performing data analysis on wavelet decomposition results of the plurality of different states to determine filtering parameters, wherein the filtering parameters comprise the decomposition layer number and an effective signal layer;
collecting a photoelectric volume pulse wave to be detected of a diver as an original signal;
respectively carrying out wavelet decomposition on the original signal by utilizing the preset wavelet basis according to the decomposition layer number in the filtering parameter to obtain a plurality of signal components;
setting coefficients corresponding to interference layers except the effective signal layer in the filtering parameters in all the signal components to zero, and performing wavelet reconstruction to obtain corresponding reconstructed signals;
and extracting characteristic points of the reconstructed signal and detecting to obtain the blood oxygen saturation.
2. The method according to claim 1, wherein the data analysis of the wavelet decomposition results of several different states to determine the filtering parameters comprises:
and performing wavelet decomposition on each sample photoplethysmogram pulse wave by using a preset wavelet base according to the current filtering parameters and reconstructing to obtain corresponding reconstruction signals, adjusting the decomposition layer number and/or effective signal layer until the reconstruction errors between all sample photoplethysmogram pulse waves and the respective corresponding reconstruction signals are within a preset error range, and determining to obtain the filtering parameters.
3. The method according to claim 1 or 2, wherein a db5 wavelet is used as the predetermined wavelet basis, and the filter parameters include 5 decomposition levels.
4. The method of claim 3, wherein the filter parameters include effective signal layers comprising third and fourth layer signal components and the interference layers comprise first, second and fifth layer signal components.
5. The method of claim 1, wherein photoplethysmography pulse waves of the diver at two different wavelengths are acquired, and one of the photoplethysmography pulse waves with a higher amplitude is selected as the original signal.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111025637.2A CN113616203A (en) | 2021-09-02 | 2021-09-02 | A method for detecting underwater blood oxygen of divers based on wavelet filtering algorithm |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111025637.2A CN113616203A (en) | 2021-09-02 | 2021-09-02 | A method for detecting underwater blood oxygen of divers based on wavelet filtering algorithm |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113616203A true CN113616203A (en) | 2021-11-09 |
Family
ID=78388868
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111025637.2A Pending CN113616203A (en) | 2021-09-02 | 2021-09-02 | A method for detecting underwater blood oxygen of divers based on wavelet filtering algorithm |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113616203A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN120167910A (en) * | 2025-05-22 | 2025-06-20 | 广东海洋大学 | Diver status monitoring method based on biofeedback |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101897591A (en) * | 2010-08-03 | 2010-12-01 | 秦皇岛市康泰医学系统有限公司 | Digital separation integrated blood oxygen probe |
| US20160361041A1 (en) * | 2015-06-15 | 2016-12-15 | The Research Foundation For The State University Of New York | System and method for infrasonic cardiac monitoring |
| CN107223036A (en) * | 2017-04-10 | 2017-09-29 | 深圳市汇顶科技股份有限公司 | Wearable device, method and device for adaptive filtering motion interference |
| CN107784278A (en) * | 2017-10-18 | 2018-03-09 | 清华大学 | With the prior-constrained sparse graph that improves of structuring as reconstruction accuracy reduces complexity method |
| CN108444713A (en) * | 2018-05-09 | 2018-08-24 | 济南大学 | A kind of Rolling Bearing Fault Character extracting method based on DShi wavelet energy bases |
| CN108920420A (en) * | 2018-03-23 | 2018-11-30 | 同济大学 | A kind of Wavelet noise-eliminating method suitable for driving evaluation test data processing |
| CN111603174A (en) * | 2020-05-11 | 2020-09-01 | 浙江荷清柔性电子技术有限公司 | Blood oxygen detection method and device, electronic equipment and storage medium |
| CN112089405A (en) * | 2020-08-26 | 2020-12-18 | 上海大学 | Pulse wave characteristic parameter measuring and displaying device |
| CN112998704A (en) * | 2021-03-04 | 2021-06-22 | 研和智能科技(杭州)有限公司 | Wearable device blood oxygen saturation calculation method |
-
2021
- 2021-09-02 CN CN202111025637.2A patent/CN113616203A/en active Pending
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101897591A (en) * | 2010-08-03 | 2010-12-01 | 秦皇岛市康泰医学系统有限公司 | Digital separation integrated blood oxygen probe |
| US20160361041A1 (en) * | 2015-06-15 | 2016-12-15 | The Research Foundation For The State University Of New York | System and method for infrasonic cardiac monitoring |
| CN107223036A (en) * | 2017-04-10 | 2017-09-29 | 深圳市汇顶科技股份有限公司 | Wearable device, method and device for adaptive filtering motion interference |
| CN107784278A (en) * | 2017-10-18 | 2018-03-09 | 清华大学 | With the prior-constrained sparse graph that improves of structuring as reconstruction accuracy reduces complexity method |
| CN108920420A (en) * | 2018-03-23 | 2018-11-30 | 同济大学 | A kind of Wavelet noise-eliminating method suitable for driving evaluation test data processing |
| CN108444713A (en) * | 2018-05-09 | 2018-08-24 | 济南大学 | A kind of Rolling Bearing Fault Character extracting method based on DShi wavelet energy bases |
| CN111603174A (en) * | 2020-05-11 | 2020-09-01 | 浙江荷清柔性电子技术有限公司 | Blood oxygen detection method and device, electronic equipment and storage medium |
| CN112089405A (en) * | 2020-08-26 | 2020-12-18 | 上海大学 | Pulse wave characteristic parameter measuring and displaying device |
| CN112998704A (en) * | 2021-03-04 | 2021-06-22 | 研和智能科技(杭州)有限公司 | Wearable device blood oxygen saturation calculation method |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN120167910A (en) * | 2025-05-22 | 2025-06-20 | 广东海洋大学 | Diver status monitoring method based on biofeedback |
| CN120167910B (en) * | 2025-05-22 | 2025-08-01 | 广东海洋大学 | Diver state monitoring method based on biofeedback |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107157492B (en) | Embedded human physiological information noninvasive detection system and data processing method | |
| CN106037694B (en) | A kind of continuous blood pressure measurer based on pulse wave | |
| Castillo et al. | Efficient wavelet-based ECG processing for single-lead FHR extraction | |
| KR102532764B1 (en) | Apparatus and method for estimating biophysiological rates | |
| TWI505816B (en) | Detecting method and apparatus for blood oxygen saturation | |
| CN109864713B (en) | Heart rate monitoring method based on multi-channel parallel filtering and spectral peak weighting selection algorithm | |
| Tang et al. | PPG signal reconstruction using a combination of discrete wavelet transform and empirical mode decomposition | |
| CN113925482B (en) | Heart rate calculating method, wearable electronic device and storage medium | |
| Raghuram et al. | Dual-tree complex wavelet transform for motion artifact reduction of PPG signals | |
| CN101732050A (en) | Photoelectric volume wave-based breathing rate monitoring method | |
| CN106618542A (en) | Denoising heart rate detecting device and method | |
| CN108937878A (en) | A kind of method that pulse wave signal motion artifacts are eliminated | |
| CN102885616A (en) | Method for removing baseline drift from pulse wave signal | |
| CN102631198A (en) | Dynamic spectrum data processing method based on difference value extraction | |
| CN102274014A (en) | Korotkoff sound blood pressure measurement and signal analysis method based on wavelet denoising | |
| CN106264505A (en) | A kind of heart rate spectral peak system of selection based on support vector machine | |
| CN112998704A (en) | Wearable device blood oxygen saturation calculation method | |
| Seljuq et al. | Selection of an optimal mother wavelet basis function for ECG signal denoising | |
| CN101683260A (en) | Rhythmicity physiological signal optimization method and device | |
| CN108403094A (en) | A method for identifying pulse wave peaks | |
| CN111493821B (en) | PPG signal real-time denoising method based on MODWT and median filtering | |
| Park et al. | Energy-efficient photoplethysmogram compression to estimate heart and respiratory rates simultaneously | |
| CN113616203A (en) | A method for detecting underwater blood oxygen of divers based on wavelet filtering algorithm | |
| Elbuni et al. | ECG parameter extraction algorithm using (DWTAE) algorithm | |
| CN110432885A (en) | A kind of photoplethysmographic noise remove method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20211109 |