CN114159045A - High-robustness automatic respiratory signal measuring system - Google Patents
High-robustness automatic respiratory signal measuring system Download PDFInfo
- Publication number
- CN114159045A CN114159045A CN202010947437.1A CN202010947437A CN114159045A CN 114159045 A CN114159045 A CN 114159045A CN 202010947437 A CN202010947437 A CN 202010947437A CN 114159045 A CN114159045 A CN 114159045A
- Authority
- CN
- China
- Prior art keywords
- signal
- respiratory
- ppg
- analysis
- signals
- 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
- 230000000241 respiratory effect Effects 0.000 title claims abstract description 53
- 230000029058 respiratory gaseous exchange Effects 0.000 claims abstract description 38
- 238000004458 analytical method Methods 0.000 claims abstract description 34
- 238000000605 extraction Methods 0.000 claims abstract description 29
- 238000005259 measurement Methods 0.000 claims abstract description 11
- 230000036391 respiratory frequency Effects 0.000 claims abstract description 11
- 230000036387 respiratory rate Effects 0.000 claims description 52
- 230000004927 fusion Effects 0.000 claims description 21
- 238000001514 detection method Methods 0.000 claims description 19
- 238000010183 spectrum analysis Methods 0.000 claims description 17
- 238000001228 spectrum Methods 0.000 claims description 8
- 238000011156 evaluation Methods 0.000 claims description 7
- 239000000284 extract Substances 0.000 claims description 6
- 238000012795 verification Methods 0.000 claims description 5
- 239000008280 blood Substances 0.000 claims description 4
- 210000004369 blood Anatomy 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 3
- 238000013441 quality evaluation Methods 0.000 claims description 3
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 claims description 2
- 230000001121 heart beat frequency Effects 0.000 claims description 2
- 229910052760 oxygen Inorganic materials 0.000 claims description 2
- 239000001301 oxygen Substances 0.000 claims description 2
- 238000013186 photoplethysmography Methods 0.000 abstract description 63
- 238000000034 method Methods 0.000 abstract description 22
- 238000012544 monitoring process Methods 0.000 abstract description 10
- 206010051379 Systemic Inflammatory Response Syndrome Diseases 0.000 abstract description 4
- 238000012545 processing Methods 0.000 abstract description 4
- 208000001647 Renal Insufficiency Diseases 0.000 abstract description 2
- 201000006370 kidney failure Diseases 0.000 abstract description 2
- 208000008784 apnea Diseases 0.000 abstract 1
- 238000011161 development Methods 0.000 abstract 1
- 230000006806 disease prevention Effects 0.000 abstract 1
- 238000000691 measurement method Methods 0.000 abstract 1
- 230000000474 nursing effect Effects 0.000 abstract 1
- 230000003449 preventive effect Effects 0.000 abstract 1
- 238000013459 approach Methods 0.000 description 4
- 238000000354 decomposition reaction Methods 0.000 description 4
- 208000010496 Heart Arrest Diseases 0.000 description 3
- 238000013528 artificial neural network Methods 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 239000002184 metal Substances 0.000 description 3
- 230000000737 periodic effect Effects 0.000 description 3
- 230000001174 ascending effect Effects 0.000 description 2
- 230000000747 cardiac effect Effects 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 230000000670 limiting effect Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000000513 principal component analysis Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000001931 thermography Methods 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- 239000013598 vector Substances 0.000 description 2
- 206010038669 Respiratory arrest Diseases 0.000 description 1
- 208000004301 Sinus Arrhythmia Diseases 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000000537 electroencephalography Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000007500 overflow downdraw method Methods 0.000 description 1
- 230000010412 perfusion Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000001303 quality assessment method Methods 0.000 description 1
- 230000003252 repetitive effect Effects 0.000 description 1
- 208000024584 respiratory abnormality Diseases 0.000 description 1
- 230000000284 resting effect Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
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/08—Detecting, measuring or recording devices for evaluating the respiratory organs
- A61B5/0816—Measuring devices for examining respiratory frequency
-
- 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, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
- A61B5/1455—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters
- A61B5/14551—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, 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
- A61B5/7207—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise induced by motion artifacts
-
- 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/725—Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
-
- 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/7257—Details of waveform analysis characterised by using transforms using Fourier transforms
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Public Health (AREA)
- Surgery (AREA)
- Veterinary Medicine (AREA)
- General Health & Medical Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Physiology (AREA)
- Signal Processing (AREA)
- Psychiatry (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Pulmonology (AREA)
- Mathematical Physics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Optics & Photonics (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
The application provides a high robustness respiratory signal automatic measurement system, including Photoelectric Plethysmograph (PPG) signal acquisition module, and Photoelectric Plethysmograph (PPG) signal extraction and analysis module. The invention obtains photoplethysmography pulse wave signals through photoplethysmography (PPG) equipment, and carries out algorithm processing analysis with high robustness on the photoplethysmography (PPG) signals to obtain accurate respiration signals and give out respiration frequency (RR). The method improves the accuracy and robustness of respiratory signal measurement, and the measurement method is full-automatic and has accurate and reliable measurement data. The invention is used for disease prevention of heart and apnea, Systemic Inflammatory Response Syndrome (SIRS), renal failure and the like which take respiratory signals and respiratory frequency as important monitoring indexes, and has important clinical significance. And because the automatic respiratory measuring equipment provided by the invention has low cost, the automatic respiratory measuring equipment can be widely applied to places such as general wards, community hospitals, nursing homes and the like, and the development of preventive medical health careers is promoted.
Description
Technical Field
The application relates to the technical field of physiological information processing, in particular to a system and a device for monitoring and implementing respiratory signals on line.
Background
Respiratory signal and Respiratory Rate (RR) are among the indices of hospital prevention of abnormal conditions, such as cardiac and respiratory arrest, Systemic Inflammatory Response Syndrome (SIRS), renal failure, etc. The normal respiration rate in adults is 8-20 beats per minute (bpm). We found that in a study on respiratory abnormalities, 54% of patients with cardiac arrest had a respiratory frequency (RR) >27bpm at least three days before cardiac arrest. Therefore, it is crucial for monitoring the breathing of the patient. Unlike other vital signs, in most cases the measurement of the Respiratory Rate (RR) is performed manually, since currently widely available automated devices meeting clinical requirements are used for the measurement of respiratory rate. The importance of Respiratory Rate (RR) and the lack of technical means at present indicate that finding an automatic and continuous respiratory monitoring method which meets clinical requirements is of great clinical significance. Currently, methods used clinically include capnography (capnography), impedance plethysmography (impedance plethysmography), and flow thermography (flow thermography), which are expensive and invasive and not suitable for patients in general wards.
With the advancement of photoplethysmography, photoplethysmography (PPG) is widely used for respiratory monitoring both inside and outside of hospitals. Scientists have given a number of methods to estimate the respiratory rate from photoplethysmography (PPG) signals, such as the earliest Wavelet decomposition (Wavelet decomposition), digital filtering (digital filters), fourier transform (FFT), and complex demodulation (complex decomposition), classical modal decomposition and autoregressive models (AR), to the recent Principal Component Analysis (PCA) and artificial Neural Networks (NN), which have been successfully applied to various photoplethysmography (PPG) databases, with good estimates. There are studies using neural networks to analyze the three modulation signals to select the best waveform for the final algorithm design. There are also studies using data fusion in combination with estimates of multiple modulation signals. Despite such many advances, the use of pulse oximeters to measure respiratory signals has only recently been adopted for commercial use due to the lack of accuracy with which the above methods estimate respiratory signals. It is therefore important to provide a reliable method. A more common approach today is to introduce a Signal Quality Index (SQI) to evaluate the signal quality, and not to perform an algorithmic estimation of the photoplethysmography (PPG) signal if it does not carry meaningful physiological information. Lack of quality indicators may lead to serious clinical errors, and introduction of evaluation indicators may improve accuracy and reliability. That is, there is no respiratory signal monitoring method that is safe, accurate and convenient for wide application in clinic.
Disclosure of Invention
The application aims to provide a high-robustness automatic respiratory signal measuring system, and provides a method and equipment which are accurate, convenient to use and capable of monitoring on line for measuring by using a photoplethysmograph (PPG) device and obtaining a respiratory signal and Respiratory Rate (RR) with high accuracy from a PPG signal.
To achieve the above object, the present invention provides a highly robust respiratory signal automatic measurement system, comprising:
a PPG signal acquisition module 100, which acquires a photoplethysmograph (PPG) signal by a photoplethysmograph device;
the PPG signal extraction and analysis module 200 extracts three respiratory modulation signals, namely, peak sequence amplitude (rriiv), continuous peak time difference (RIFV) and pulse height (RIAV), performs spectral analysis on the respiratory modulation signals, performs data fusion on the signals and data obtained by the spectral analysis, and finally obtains accurate and reliable respiratory signals and Respiratory Rate (RR).
Preferably, the PPG signal acquisition module comprises a photoplethysmograph device employing a wearable pulse oximeter, estimating the blood oxygen saturation (SpO2) according to Beer-Lambert's law, using the photoplethysmogram to display the change of blood volume in the finger over time, including a pulse component and a constant component, the respiratory and heartbeat frequencies being included in the pulse component of the photoplethysmograph (PPG) signal, and transmitting the signal to a PPG signal extraction and analysis module;
preferably, the PPG signal extraction and analysis module comprises: a respiration modulation signal extraction unit which extracts and evaluates three respiration modulation signals of peak sequence amplitude (RRIV), continuous peak time difference (RIFV) and pulse height (RIAV); the signal analysis unit is used for carrying out spectrum analysis on the three modulation signals; and the data fusion unit is used for performing data fusion on the result of the frequency spectrum analysis.
Preferably, the respiration signal modulation extraction unit includes: three respiration modulation signals of peak value sequence amplitude (RRIV), continuous peak value time difference (RIFV) and pulse height (RIAV) are obtained by adopting peak value detection, peak value verification and trough detection, and signal quality evaluation is carried out on the signals, wherein the evaluation method comprises an increment merging subdivision algorithm (IMS), double pulse detection and Signal Quality Index (SQI) calculation.
Preferably, the signal analyzing unit includes: performing spectrum analysis on the three modulation signals, including high-efficiency Fourier transform (FFT), autocorrelation Analysis (AC) and establishing an autoregressive model to obtain three respiratory signals and Respiratory Rate (RR) pre-estimated values;
preferably, the data fusion unit includes: and performing data fusion on the estimated signal and the data obtained by the frequency spectrum analysis to obtain a high-accuracy respiratory signal and Respiratory Rate (RR).
The high-robustness automatic respiratory signal measuring system provided by the invention is simple in equipment operation, low in cost, accurate and reliable in result, and suitable for online monitoring of respiratory frequency in various places.
The invention relates to a high robustness respiration signal automatic measuring system which is characterized in that:
1) the device can be used for monitoring the respiratory frequency on line in various places such as hospitals, community hospitals and families;
2) the equipment is simple and safe to operate, high in automation degree and accurate and reliable in result;
3) the equipment cost is low, the limitation by fund is small, and the device can be used in places such as common medical institutions and the like.
Drawings
In order to more clearly illustrate the technical solution of the embodiment of the present application, a structural block diagram of a high-robustness automatic respiratory signal measurement system in fig. 1 is provided. The drawings illustrate only certain embodiments of the application and are therefore not to be considered limiting of scope, for those skilled in the art will appreciate that other related drawings may be derived therefrom without the exercise of inventive faculty.
Detailed Description
The technical solutions in the embodiments of the present application will be described below with reference to the drawings in the embodiments of the present application, and it should be noted that the described embodiments are only intended to facilitate understanding of the present invention, and do not have any limiting effect.
The high-robustness respiratory signal automatic measurement system in the embodiment of the application is composed of a PPG signal acquisition module 100 and a PPG signal extraction and analysis module 200.
The PPG signal acquisition module is a wearable pulse oximeter 100. Wearable pulse oximeter 100 collects the photoplethysmograph (PPG) signals, transmits them to PPG signal extraction and analysis module 200 via wifi.
The PPG signal extraction and analysis module 200 is composed of a respiratory modulation signal extraction unit 210, a signal analysis unit 220, and a data fusion unit 230. The respiratory modulation signal extraction unit 210 extracts three respiratory modulation signals, namely peak sequence amplitude (RRIV), continuous peak time difference (RIFV) and pulse height (RIAV), from a photoplethysmograph (PPG) signal; the signal analysis unit 220 analyzes the extracted respiration modulation signal to obtain three estimated respiration signals and estimated respiration frequency (RR) values; the data fusion unit 230 performs data fusion on the obtained estimated respiratory signal and the estimated Respiratory Rate (RR) value to finally obtain an accurate and reliable respiratory signal and Respiratory Rate (RR).
The following detailed description of embodiments of the invention:
The signal acquisition module 100 of the photoplethysmography (PPG) is worn at the finger end by the wearable pulse oximeter 111, driven by a firmware program, and transmits the acquired signal of the photoplethysmography (PPG) to the PPG signal extraction and analysis module 200 through Wifi via the data acquisition and transmission unit 110.
(II) PPG Signal extraction and analysis Module 200
The PPG signal extraction and analysis module 200 runs on a computer and is composed of a respiration modulation signal extraction unit 210, a signal analysis unit 220, and a data fusion unit 230.
The respiratory modulation signal extraction unit 210 extracts three respiratory modulation signals, namely peak sequence amplitude (RRIV), continuous peak time difference (RIFV) and pulse height (RIAV), from a photoplethysmograph (PPG) signal; the signal analysis unit 220 analyzes the extracted respiration modulation signal to obtain an estimated respiration signal and an estimated respiration frequency (RR); the data fusion unit 230 performs data fusion on the obtained result to finally obtain an accurate and reliable respiratory signal and Respiratory Rate (RR). The specific analysis method is as follows:
1. extraction of respiratory modulation signals
A. Method idea
a. The respiratory frequency (RR) is extracted from a photoplethysmograph (PPG), which is generally segmented using sliding windows, each window yielding a respective respiratory frequency (RR). The first step of the algorithm is typically to extract the respiration modulation signal from a photoplethysmograph (PPG) signal, which has many methods, the most common being the detection of peaks and troughs of the photoplethysmograph (PPG) signal. We define the photoplethysmograph (PPG) signal peak time series as a set of points tpk,i,yPk,i}i=1...NpkDefining a trough time series as a set of points { t }tr,i,ytr,i}i=1...NtrNpk, Ntr are the number of peaks and troughs respectively. Due to noise in the signal or false detection in the detection algorithm, Npk≠Ntr。
b. The peak-to-valley time series will be used to generate three different breath modulation signals representing three different pieces of information for a breath.
1) Peak sequence amplitude (RRIV): the breath-guided intensity variation is the amplitude of the sequence of peaks. This is due to changes in intrathoracic pressure resulting in changes in perfusion baseline, which are manifested as changes in amplitude of the photoplethysmograph (PPG) peak. Therefore, the temperature of the molten metal is controlled,
2) pulse height (RIAV): the respiratory-guided amplitude variation, i.e. the height of the photoplethysmograph (PPG) pulse, is defined as the height difference between two adjacent peak troughs. This is due to respiration-induced cardiac output changes, resulting in changes in the absolute amplitude of the photoplethysmography (PPG) signal. Therefore, the temperature of the molten metal is controlled,
3) continuous peak time difference (RIFV): respiratory induction frequency changes are changes in the instantaneous HR value of the respiratory cycle, a phenomenon known as Respiratory Sinus Arrhythmia (RSA). Defined as the time difference between successive peaks. Therefore, the temperature of the molten metal is controlled,still other respiratory modulation signals are also given, such as pulse width variations. These modulation signals can be used to estimate the Respiratory Rate (RR).
c. Preprocessing operations are typically required before data processing can take place. High frequency noise is first filtered out using low pass filtering. The photoplethysmograph (PPG) signal is then captured using an incremental merging and subdividing algorithm (IMS). The respiratory modulation signal is then extracted using the three approaches mentioned above. Due to the non-uniform time sequence, the signal needs to be resampled at fs-4 Hz. To unify the metrics of the three signals, normalization with zero mean-one variance was used. In order to increase the reliability of the signal, a signal evaluation index (signal quality index SQI) is added to evaluate the modulation signal.
B. For Respiratory Rate (RR) estimation, the most critical is the extraction of the respiratory modulation signal. The main methods for extracting the respiration modulation signal are peak detection and signal quality evaluation. The present invention presents an improved peak detection method and signal evaluation method, followed by estimation of respiratory signal and frequency using spectral analysis and data fusion.
a. Peak detection
The general principle of peak detection is that any singular point of a differentiable signal corresponds to a zero crossing or two inflection points in its derivative signal. The left side of the peak value is raised and the right side is lowered intuitively. The general approach to peak detection is to have the first derivative be zero and the second derivative be less than zero. There is also a method to provide that the maximum value of the original signal is the peak value between two inflection points of the first derivative, i.e. the subscripts corresponding to the maximum value and the minimum value. The invention provides a new method for peak detection without a second derivative and without judging the inflection point of the first derivative. In order to avoid the influence of the amplitude change of the front and back of the photoplethysmograph (PPG) signal on peak detection and verification, windowing is carried out on the photoplethysmograph (PPG) signal, a 10s time window is used, and the overlapping time is 5s (avoiding missing part of peaks).
Peak verification: peak verification is primarily initiated from two aspects, the amplitude threshold and the time interval threshold. In terms of amplitude threshold setting, let thresh1 be the ninety-th quantile, thresh2 be the first ten-decile, and then,then orderWherein peaks is a peak point. If the highdiff is less than middlediff&If highdiff < lowdiff, the peak point is recorded. In terms of time interval, since the pulse wave is mainly controlled by the heartbeat, the normal resting heart rate is 60-100, and the peak value in this range is detected.
Wave trough detection: when a Photoelectric Plethysmograph (PPG) signal is observed, the minimum value between two peak values is the trough.
b. Signal quality assessment (incremental Merge Subdivision Algorithm IMS, Dual pulse recognition, Signal Quality Index (SQI))
Since motion artifacts (artifacts) and noise that cannot be filtered out exist in the photoplethysmograph (PPG) signal, an evaluation of the quality of the photoplethysmograph (PPG) signal is required. The method for evaluating the signal quality comprises the steps of dividing a pulse of a photoplethysmography (PPG) into line segments (convenient to process and judge) by using an IMS (IP multimedia subsystem) algorithm, and judging the proportion of an artifact and a clip in the double pulse as a Signal Quality Index (SQI) to judge the signal quality after the double pulse is processed.
IMS (Incremental-Merge Segmentation increment merging subdivision) algorithm: the IMS algorithm has a sliding window structure and can be used for real-time processing. The algorithm only needs to set a parameter m (the number of points moved each time, mainly related to the sampling rate). The principle is that the Photoelectric Plethysmograph (PPG) signal is divided into n segments with the length of m, the slopes of the segments are calculated, if the slopes have the same slope, the segments are combined (the starting point of the first line segment and the end point of the second line segment are connected into a straight line), and if the slopes have different slopes, the segments are divided into new line segments.
Double pulse detection: by judging the slope as a negative line segment and a subsequent line segment, if a certain condition is met, the pulse is double pulse. After the double pulse is judged, the three line segments of the double pulse are merged.
Signal Quality Index (SQI): the rise segments are in one-to-one correspondence with the fall segments, and therefore, the rise segments are analyzed separately. If the threshold and the slope of the ascending line segment both exceed the threshold, judging the ascending line segment to be artifact; if the slope is zero, determining that the slope is clip; the segment immediately following the clip is also artifact; and calculating the ratio of artifact to clip in a photoplethysmography (PPG) signal, calculating a Signal Quality Index (SQI), and judging the signal quality.
2. Signal analysis (FFT, autocorrelation, autoregressive, autocorrelation + autoregressive)
After obtaining the peak-to-valley sequence from the peak detection, three respiratory modulation signals, peak sequence amplitude (RRIV), continuous peak time difference (RIFV), and pulse height (RIAV), were calculated using the methods described above. Since a non-uniform time series is obtained, it is resampled f by linear interpolation for the purpose of subsequent spectral analysiss4 Hz. And meanwhile, zero-mean square difference standardization is used, so that the amplitudes of the three modulation signals are unified to the same range. Then, a high-pass filter is used to filter out the low-frequency signal interference, and a moving average filter is used to smooth the signal.
FFT: performing efficient Fourier transform on the three modulation signals, and selecting the frequency with the maximum amplitude from a frequency domain as a respiratory frequency (RR); because the number of modulation signal points is less, zero padding processing is carried out on the modulation signals in order to improve the frequency spectrum resolution.
AC: autocorrelation analysis is a mathematical tool that looks for repetitive patterns, such as looking for periodic signals that are masked by noise. The breathing signal can be regarded as a noisy periodic signal. The autocorrelation formula isWherein xiFor time series of signals, xi+τFor time series shifted by τ units, μ being the mean, σ2Is the variance. An autocorrelation sequence can be obtained from τ -0 to τ -n-1. The abscissa is the number of translation points tau, the relationship between the frequency and the number of translation points is f ═ fs/tau, and fs is the sampling frequency. In the autocorrelation signal, each peak (except the first one) represents a period of strong autocorrelation, and the period with the largest correlation is taken as the Respiratory Rate (RR).
AR: the AR model is an alternative method to Discrete Fourier Transform (DFT) and is also one of the methods for short-time series high-resolution spectral estimation. In biomedical engineering, AR models are widely used for spectral analysis of heart rate variability and electroencephalography. In the AR model, each value of the time series is a regression to its past values. The number of past values used is referred to as the order of the model. The AR model can be viewed as a filter that divides the time series into a predictable time series and a prediction error series. It provides a smoother, more easily understood power spectrum than the discrete fourier transform, but also has the disadvantage of complex model identification. The AR model is defined as:
wherein M is the model order, aiAs a weight, ε [ n ]]Is a prediction error term and obeys ε -N (0, σ). Using least squares to make the prediction error epsilon n]Minimizing to obtain optimal parameter ai。
The above equation is developed
x[n]=a1x[n-1]+a2x[n-2]+...+aMx[n-M]+ε[n] (2)
The first N predictions are:
m unknowns require M equations that can be solved using the least squares method. Matrixing the above formula:
x=Xa+ε (4)
wherein X is a matrix with M rows and M columns, and a and epsilon are M-dimensional column vectors.
The optimum parameter a is requiredoptI.e. the minimum prediction error n is required]Namely:
ε=x-Xaopt0; (5) i.e. by
(XTX)-1,
(XTX)-1(XTX)aopt=aopt=(XTX)-1XTx (7)
At the same time, due to XTX is very close to the autocorrelation matrix of X,
XTx is very close to the autocorrelation vector of x,
obtained by synthesizing (7), (8) and (9), aopt=R-1r (10), this equation is called Yule-Walker equation, also called autocorrelation. The use of a recursive approach greatly facilitates solving for aopt。
Another important point of the autoregressive model is the choice of the order M of the model, with different order fits of the model being good or bad. In practical application, a plurality of models are fitted to the sequence, and different models are selected according to effects. The most common selection criteria is Akaike's Information Criterion (AIC), expressed as:
whereinIs the prediction error variance associated with M. The M that minimizes AIC is the best model order.
Meanwhile, according to a Yule-Walker equation:
σ2=[r(0) r(1) r(2)…r(M)][1 a1 a2…aM]T。
from the z field, formula (1) can be written as:
And because z is ejwTTo obtainThen, the frequency spectrum R (e) of the time seriesjw) Can be obtained by multiplying the square of the transfer function by the variance of the prediction error. Namely, it is
AC + AR: as can be seen from the characteristics of the autocorrelation method, autocorrelation can remove part of the noise in the periodic signal. Each peak of the autocorrelation coefficients represents a period of strong autocorrelation, so that the period of autocorrelation corresponds to the period of the respiration signal. Therefore, the respiratory frequency (RR) can be estimated using the autocorrelation signal of the respiratory modulation signal as an input signal for AR spectral analysis.
3. Data fusion
For the estimated respiratory signal and Respiratory Rate (RR) obtained by the three respiratory modulation signals, a data fusion can be performed in order to improve the accuracy and reliability. The common fusion method is to directly average the frequency spectrums of the three modulation signals, and take the maximum value as the Respiratory Rate (RR). The present invention provides two improved methods:
first, from the above spectral analysis, for each kind of respiration modulation signal, a Respiration Rate (RR) is obtained, which is respectively a Respiration Rate (RR) _ peak sequence amplitude (RRIV), a Respiration Rate (RR) _ pulse height (RIAV), a Respiration Rate (RR) _ continuous peak time difference (RIFV), and if the absolute value of the difference between the three Respiration Rates (RR) is within 2bpm, the average value of the three signal spectra is determined; if the absolute value of the difference value of any two respiratory frequencies (RR) is within 2bpm, taking the average value of the frequency spectrum; otherwise, the Respiration Rate (RR) _ continuous peak time difference (RIFV) is taken (because it is observed from the experiment that the waveform of the continuous peak time difference (RIFV) is closer to the respiration waveform).
Secondly, as the autocorrelation signal of the modulation signal can be used as a respiration signal to be analyzed, and the waveform is more regular, the autocorrelation signal is segmented by using an IMS algorithm (as the autocorrelation signal of the normal respiration signal is approximate to a sine wave, the variance value is smaller, and the mean value is close to 1, namely the waveform is more stable), variance analysis is carried out on the decomposed segment, the Signal Quality Index (SQI) ═ var/mean is used as an index for measurement, and the modulation signal of which the Signal Quality Index (SQI) is smaller than a certain range is taken as frequency spectrum average.
To sum up, the high-robustness respiration signal automatic measurement system provided by the embodiment of the application acquires the PPG signal through the PPG signal acquisition module, and analyzes the signal through the PPG signal extraction and analysis module to finally obtain the accurate and reliable respiration signal and the respiration frequency. Accurate and effective respiratory signal online monitoring data can be provided for clinical use.
In the embodiments provided in the present application, it should be understood that the disclosed apparatus and method may be implemented in other ways. The above-described embodiments of the apparatus are merely illustrative, and for example, the division of the units is only one logical division, and there may be other divisions in actual implementation, and for example, a plurality of units may be combined or integrated into another module, or a unit of a module may further refine more units.
The above description is only an example of the present application and is not intended to limit the scope of the present application, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, improvement and the like made within the spirit and principle of the present application shall be included in the protection scope of the present application.
Claims (6)
1. A highly robust respiratory signal automatic measurement system, comprising:
a PPG signal acquisition module 100, which acquires a photoplethysmograph (PPG) signal by a photoplethysmograph device;
a PPG signal extraction and analysis module 200; the PPG signal extraction and analysis module 200 comprises the following units: a respiratory signal extraction unit 210 that extracts three respiratory modulation signals, namely, a peak sequence amplitude (RRIV), a continuous peak time difference (RIFV), and a pulse height (RIAV); a signal analysis unit 220 that performs spectrum analysis on the respiration modulation signal; the data fusion unit 230 performs data fusion on the signals and data obtained by the spectrum analysis, and finally obtains accurate and reliable respiratory signals and Respiratory Rate (RR).
2. The system as claimed in claim 1, wherein the PPG signal acquisition module comprises:
the photoplethysmograph device employs a wearable pulse oximeter, estimates the blood oxygen saturation (SpO2) according to Beer-Lambert's law, uses a photoplethysmogram to display the change in blood volume in the finger over time, including a pulse component and a constant component, with the respiratory and heartbeat frequencies included in the pulse component of the photoplethysmograph (PPG) signal, and transmits the signal to a PPG signal extraction and analysis module.
3. The system as claimed in claim 1, wherein the PPG signal extraction and analysis module comprises:
a respiration modulation signal extraction unit which extracts and evaluates three respiration modulation signals of peak sequence amplitude (RRIV), continuous peak time difference (RIFV) and pulse height (RIAV);
the signal analysis unit is used for carrying out spectrum analysis on the three modulation signals;
and the data fusion unit is used for performing data fusion on the result of the frequency spectrum analysis.
4. The PPG signal extraction and analysis module of claim 3, wherein the respiratory signal modulation extraction unit comprises:
three respiration modulation signals of peak value sequence amplitude (RRIV), continuous peak value time difference (RIFV) and pulse height (RIAV) are obtained by adopting peak value detection, peak value verification and trough detection, and signal quality evaluation is carried out on the signals, wherein the evaluation method comprises an increment merging subdivision algorithm (IMS), double pulse detection and Signal Quality Index (SQI) calculation.
5. The PPG signal extraction and analysis module of claim 3, wherein the signal analysis unit comprises:
and performing spectrum analysis on the three modulation signals, including high-efficiency Fourier transform (FFT), autocorrelation Analysis (AC) and establishing an autoregressive model to obtain three respiratory signals and Respiratory Rate (RR) pre-estimated values.
6. The PPG signal extraction and analysis module of claim 3, wherein the data fusion unit further comprises:
and performing data fusion on the estimated signal and data obtained by the spectrum analysis, performing spectrum evaluation or autocorrelation signal quality analysis, and obtaining a high-accuracy respiratory signal and Respiratory Rate (RR).
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010947437.1A CN114159045A (en) | 2020-09-10 | 2020-09-10 | High-robustness automatic respiratory signal measuring system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010947437.1A CN114159045A (en) | 2020-09-10 | 2020-09-10 | High-robustness automatic respiratory signal measuring system |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114159045A true CN114159045A (en) | 2022-03-11 |
Family
ID=80475641
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010947437.1A Pending CN114159045A (en) | 2020-09-10 | 2020-09-10 | High-robustness automatic respiratory signal measuring system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114159045A (en) |
-
2020
- 2020-09-10 CN CN202010947437.1A patent/CN114159045A/en active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Pimentel et al. | Toward a robust estimation of respiratory rate from pulse oximeters | |
US11103155B2 (en) | Method of determining the frequency of a periodic physiological process of a subject, and a device and system for determining the frequency of a periodic physiological process of a subject | |
Todd et al. | The identification of peaks in physiological signals | |
Fleming et al. | A comparison of signal processing techniques for the extraction of breathing rate from the photoplethysmogram | |
Shah et al. | Respiratory rate estimation during triage of children in hospitals | |
Bansal et al. | A review of measurement and analysis of heart rate variability | |
Birrenkott et al. | Robust estimation of respiratory rate via ECG-and PPG-derived respiratory quality indices | |
Mirmohamadsadeghi et al. | Respiratory rate estimation from the ECG using an instantaneous frequency tracking algorithm | |
Jeyhani et al. | Comparison of simple algorithms for estimating respiration rate from electrical impedance pneumography signals in wearable devices | |
Sharma et al. | ECG-derived respiration using Hermite expansion | |
Nayan et al. | Development of respiratory rate estimation technique using electrocardiogram and photoplethysmogram for continuous health monitoring | |
Reddy et al. | Automated hilbert envelope based respiration rate measurement from PPG signal for wearable vital signs monitoring devices | |
C. Dela Cruz et al. | Deriving heart rate and respiratory rate from ecg using wavelet transform | |
Kozia et al. | ECG-derived respiration using a real-time QRS detector based on empirical mode decomposition | |
Surrel et al. | Low-power wearable system for real-time screening of obstructive sleep apnea | |
Chakraborty et al. | A robust PPG onset and systolic peak detection algorithm based on Hilbert transform | |
Dong et al. | An integrated framework for evaluation on typical ECG-derived respiration waveform extraction and respiration | |
CN114580477B (en) | Wearable dynamic respiratory rate estimation system based on multi-time sequence fusion | |
Modak et al. | A novel adaptive multilevel thresholding based algorithm for QRS detection | |
Alikhani et al. | Spectral fusion-based breathing frequency estimation; experiment on activities of daily living | |
Nayan et al. | Breathing rate estimation from a single-lead electrocardiogram acquisition system | |
Protopsaltis et al. | Continuous measurement of respiratory rate via single-wavelength reflective photoplethysmography | |
CN114159045A (en) | High-robustness automatic respiratory signal measuring system | |
CN109528214A (en) | A kind of Multifunctional wrist BOLD contrast | |
Nyamukuru et al. | Extracting Fractional Inspiratory Time from Electrocardiograms |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20220311 |
|
WD01 | Invention patent application deemed withdrawn after publication |