CN115486833A - Respiration state detection method and device, computer equipment and storage medium - Google Patents
Respiration state detection method and device, computer equipment and storage medium Download PDFInfo
- Publication number
- CN115486833A CN115486833A CN202211004088.5A CN202211004088A CN115486833A CN 115486833 A CN115486833 A CN 115486833A CN 202211004088 A CN202211004088 A CN 202211004088A CN 115486833 A CN115486833 A CN 115486833A
- Authority
- CN
- China
- Prior art keywords
- respiratory
- sequence
- feature sequence
- heart rate
- signal
- 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.)
- Granted
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/08—Detecting, measuring or recording devices for evaluating the respiratory organs
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/024—Detecting, measuring or recording pulse rate or heart rate
- A61B5/02405—Determining heart rate variability
-
- 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/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
- A61B5/1102—Ballistocardiography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/48—Other medical applications
- A61B5/4806—Sleep evaluation
-
- 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/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
- A61B5/7267—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems involving training the classification device
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Surgery (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Animal Behavior & Ethology (AREA)
- Physiology (AREA)
- Cardiology (AREA)
- Artificial Intelligence (AREA)
- Pulmonology (AREA)
- Fuzzy Systems (AREA)
- Evolutionary Computation (AREA)
- Dentistry (AREA)
- Mathematical Physics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Signal Processing (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
The application relates to a breathing state detection method, a breathing state detection device, a computer device and a storage medium, wherein the method comprises the following steps: acquiring a sign signal of a user at a preset first time scale; extracting a ballistocardiogram signal and a respiration signal from the sign signal; determining a first feature sequence corresponding to the ballistocardiogram signal and the respiration signal; determining a second sequence of features of the heart rate feature sequence, the heart rate variability feature sequence, the respiratory rate feature sequence, and the sleep feature sequence at a number of different second time scales; the second feature sequence comprises a reference feature sequence, a fluctuation feature sequence, a maximum feature sequence and a mode feature sequence; and inputting the second characteristic sequences under a plurality of different second time scales into the trained respiratory state detection model to obtain a respiratory state detection result, so that the cost is reduced, and the accuracy of respiratory state detection is improved.
Description
Technical Field
The present application relates to the field of breath detection, and in particular to a breath state detection, apparatus, computer device and storage medium.
Background
Chronic Obstructive Pulmonary Disease (COPD) is a respiratory Disease that can be prevented and treated, and is mainly characterized by the existence of continuous airflow limitation, and clinically, the COPD is manifested by symptoms of dyspnea, cough or expectoration and the like in different degrees, and the Disease progress can be accelerated without timely intervention, so that the life quality of patients is seriously affected.
At present, early screening for COPD can be assisted by detecting the respiratory state of a patient. The existing respiratory state detection method has short monitoring time and is mostly limited to daytime monitoring, so that the respiratory state detection accuracy is low and the price is high.
Disclosure of Invention
Based on this, it is an object of the present application to provide a breathing state detection, an apparatus, a computer device and a storage medium, which may reduce costs and improve the accuracy of the breathing state detection.
According to a first aspect of embodiments of the present application, there is provided a respiratory state detection method, including the steps of:
acquiring a physical sign signal of a user at a preset first time scale;
extracting ballistocardiogram signals and respiration signals from the sign signals;
determining a first feature sequence corresponding to the ballistocardiogram signal and the respiration signal; the first characteristic sequence comprises a heart rate characteristic sequence, a heart rate variability characteristic sequence, a respiratory rate characteristic sequence and a sleep characteristic sequence;
determining a second signature sequence of the heart rate signature sequence, the heart rate variability signature sequence, the respiratory rate signature sequence, and the sleep signature sequence at a number of different second time scales; the second feature sequence comprises a reference feature sequence, a fluctuation feature sequence, a maximum feature sequence and a mode feature sequence;
and inputting the second characteristic sequences under a plurality of different second time scales into a trained respiratory state detection model to obtain a respiratory state detection result.
According to a second aspect of embodiments of the present application, there is provided a respiratory state detection apparatus comprising:
the sign signal acquisition module is used for acquiring sign signals of a user under a plurality of preset different time scales;
the signal extraction module is used for extracting ballistocardiogram signals and respiratory signals under a plurality of different time scales from the physical sign signals;
the first characteristic sequence determining module is used for determining corresponding first characteristic sequences of the ballistocardiogram signal and the respiration signal under a plurality of different time scales; the first feature sequence comprises a heart rate feature sequence, a heart rate variability feature sequence, a respiratory rate feature sequence and a sleep feature sequence;
a second feature determination module, configured to determine second features corresponding to the heart rate feature sequence, the heart rate variability feature sequence, the respiratory rate feature sequence, and the sleep feature sequence at a plurality of different time scales; the second features comprise reference features, fluctuation features, most-valued features and mode features;
and the detection result obtaining module is used for inputting the second characteristics under a plurality of different time scales into the trained respiratory state detection model to obtain a plurality of respiratory state detection results under different time scales.
According to a third aspect of embodiments of the present application, there is provided a computer apparatus comprising: a processor and a memory; wherein the memory stores a computer program adapted to be loaded by the processor and to perform the breathing state detection method according to any of the above.
According to a fourth aspect of embodiments of the present application, there is provided a computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, implements the respiratory state detection method according to any one of the above.
According to the embodiment of the application, the sign signals of the user under the preset first time scale are obtained; extracting a ballistocardiogram signal and a respiration signal from the sign signal; determining a first feature sequence corresponding to the ballistocardiogram signal and the respiration signal; the first feature sequence comprises a heart rate feature sequence, a heart rate variability feature sequence, a respiratory rate feature sequence and a sleep feature sequence; determining a second sequence of features of the heart rate feature sequence, the heart rate variability feature sequence, the respiratory rate feature sequence, and the sleep feature sequence at a number of different second time scales; the second feature sequence comprises a reference feature sequence, a fluctuation feature sequence, a most-valued feature sequence and a mode feature sequence; and inputting the second characteristic sequences under a plurality of different second time scales into the trained respiratory state detection model to obtain a respiratory state detection result, so that the cost is reduced, and the accuracy of respiratory state detection is improved.
It is to be understood that both the foregoing general description and the following detailed description are exemplary and explanatory only and are not restrictive of the application.
For a better understanding and practice, the invention is described in detail below with reference to the accompanying drawings.
Drawings
Fig. 1 is a schematic flow chart of a respiratory state detection method according to an embodiment of the present application;
fig. 2 is a schematic flowchart of step S20 according to an embodiment of the present application;
fig. 3 is a schematic flowchart of step S30 according to an embodiment of the present application;
fig. 4 is a schematic flowchart of step S40 according to an embodiment of the present application;
fig. 5 is a block diagram of a respiratory state detection apparatus according to an embodiment of the present application;
fig. 6 is a block diagram illustrating a structure of an electronic device according to an embodiment of the present application.
Detailed Description
In order to make the objects, technical solutions and advantages of the present application more clear, embodiments of the present application will be described in further detail below with reference to the accompanying drawings.
It should be understood that the embodiments described are only a few embodiments of the present application, and not all embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present application.
The terminology used in the embodiments of the present application is for the purpose of describing particular embodiments only and is not intended to be limiting of the embodiments of the present application. As used in the examples of this application and the appended claims, the singular forms "a", "an", and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise. It should also be understood that the term "and/or" as used herein refers to and encompasses any and all possible combinations of one or more of the associated listed items.
When the following description refers to the accompanying drawings, like numbers in different drawings represent the same or similar elements unless otherwise indicated. The implementations described in the following exemplary examples do not represent all implementations consistent with the present application. Rather, they are merely examples of apparatus and methods consistent with certain aspects of the application, as detailed in the claims that follow. In the description of the present application, it is to be understood that the terms "first," "second," "third," and the like are used solely to distinguish one from another and are not necessarily used to describe a particular order or sequence, nor are they to be construed as indicating or implying relative importance. The specific meaning of the above terms in the present application can be understood by those of ordinary skill in the art as appropriate.
Further, in the description of the present application, "a plurality" means two or more unless otherwise specified. "and/or" describes the association relationship of the associated object, indicating that there may be three relationships, for example, a and/or B, which may indicate: a exists alone, A and B exist simultaneously, and B exists alone. The character "/" generally indicates that the former and latter associated objects are in an "or" relationship.
Example 1
Please refer to fig. 1, which is a flowchart illustrating a respiratory state detection method according to an embodiment of the present application. The respiratory state detection method provided by the embodiment of the application comprises the following steps:
s10: and acquiring a physical sign signal of the user at a preset first time scale.
In the embodiment of the present application, an execution subject of the respiratory state detection method is a respiratory state detection device (hereinafter referred to as a detection device), and the detection device may be one computer device, a server, or a server cluster formed by combining a plurality of computer devices.
The detection device can acquire the sign signals of the user under a plurality of preset different time scales by inquiring in a preset database. The detection device can also adopt an undisturbed sensor to acquire vital sign signals of the user under a plurality of preset different time scales in real time. The undisturbed sensor comprises a signal acquisition module and a data storage module.
Specifically, the signal acquisition module is placed below a user pillow, the user sleeps on the pillow, the user generates micro vibration due to heart activity, respiratory activity and the like, the center of gravity shifts, and force signals are generated. The signal acquisition module can convert the force signal into an analog electric signal, and then the analog electric signal is filtered, amplified and A/D converted into a digital signal with the sampling rate of 1000Hz, namely a physical sign signal, through a built-in filter circuit, an amplifying circuit and an A/D conversion circuit.
The preset first time scale is an average time length of the user using the detection device, for example, the preset first time scale is 100 days.
S20: extracting a ballistocardiogram signal and a respiration signal from the sign signal.
In the embodiment of the application, the sign signals are aliasing signals of the ballistocardiogram signals and the respiration signals, and the independent ballistocardiogram signals and the independent respiration signals are respectively extracted from the sign signals.
S30: determining a first feature sequence corresponding to the ballistocardiogram signal and the respiration signal; the first signature sequence includes a heart rate signature sequence, a heart rate variability signature sequence, a respiratory rate signature sequence, and a sleep signature sequence.
The Heart Rate (HR) is the number of beats per minute of a normal person in a resting state. Heart Rate Variability (HRV) refers to the variation of the difference between successive Heart cycles. The Respiratory Rate (RR) is the number of breaths per minute.
In an embodiment of the application, several cardiac intervals may be obtained from the ballistocardiogram signal and several respiratory intervals may be obtained from the respiratory signal. A heart rate characteristic sequence and a heart rate variability characteristic sequence can be determined according to the plurality of cardiac intervals, a respiratory frequency characteristic sequence can be determined according to the plurality of respiratory intervals, and a sleep characteristic sequence can be determined through the plurality of cardiac intervals and the plurality of respiratory intervals.
S40: determining a second sequence of features of the heart rate feature sequence, the heart rate variability feature sequence, the respiratory rate feature sequence, and the sleep feature sequence at a number of different second time scales; the second feature sequence comprises a reference feature sequence, a fluctuation feature sequence, a maximum feature sequence and a mode feature sequence.
In the embodiment of the present application, the second time scale is a continuous time scale, and the plurality of different second time scales may be 3 consecutive days, 5 consecutive days, 7 consecutive days, 10 consecutive days, 14 consecutive days, and 30 consecutive days. After obtaining the heart rate feature sequence, the heart rate variability feature sequence, the respiratory rate feature sequence and the sleep feature sequence at the preset first time scale, a reference feature sequence, a fluctuation feature sequence, a most significant feature sequence and a mode feature sequence at a plurality of different second time scales can be further obtained.
S50: and inputting the second characteristic sequences under a plurality of different second time scales into a trained respiratory state detection model to obtain a respiratory state detection result.
In the embodiment of the present application, the breathing state of the user includes two states, normal and abnormal, which can be represented by labels 0 and 1. And inputting the second characteristic sequences under a plurality of different second time scales into a trained respiratory state detection model, wherein the trained respiratory state detection model can output one of a normal state and an abnormal state.
By applying the embodiment of the application, the sign signals of the user under the preset first time scale are obtained; extracting a ballistocardiogram signal and a respiration signal from the sign signal; determining a first feature sequence corresponding to the ballistocardiogram signal and the respiration signal; the first feature sequence comprises a heart rate feature sequence, a heart rate variability feature sequence, a respiratory rate feature sequence and a sleep feature sequence; determining a second signature sequence of the heart rate signature sequence, the heart rate variability signature sequence, the respiratory rate signature sequence, and the sleep signature sequence at a number of different second time scales; the second feature sequence comprises a reference feature sequence, a fluctuation feature sequence, a most-valued feature sequence and a mode feature sequence; and inputting the second characteristic sequences under a plurality of different second time scales into the trained respiratory state detection model to obtain a respiratory state detection result. This application acquires user's sign signal through undisturbed sensor acquisition equipment, can monitor for a long time to the cost is reduced has improved the degree of accuracy that respiratory state detected.
In an alternative embodiment, referring to fig. 2, the step S20 includes steps S21 to S25, which are as follows:
s21: and filtering the sign signals by adopting a filter to obtain filtered sign signals.
In the embodiment of the application, a second-order Butterworth band-pass filter with the cut-off frequency of 2-20 Hz is adopted to remove noise in sign signals of a plurality of different time scales. Further, the power frequency interference and Gaussian additive noise are removed through a fourth-order Butterworth band-pass filter with the cut-off frequency of 20 Hz. Further, taking 15s as a fixed time window, taking the sliding scale as 1s, detecting the peak-to-average ratio of the signal of the time window to detect the body motion signal, when the peak-to-average ratio of the signal is greater than a threshold value, considering that the body motion signal exists in the time window, emptying the signal in the time window and re-entering the signal into a new queue signal.
S22: carrying out corrosion-first and expansion-second opening operation on the filtered sign signals to obtain first sign signals;
s23: performing closed operation of expansion and corrosion on the filtered sign signal to obtain a second sign signal;
s24: averaging the first volume sign signal and the second volume sign signal, and taking an average result as a respiratory signal;
s25: and subtracting the respiratory signal from the filtered physical sign signal to obtain a ballistocardiogram signal.
In the embodiment of the application, on the basis of eliminating the body motion signal, due to the difference of the main frequency ranges of the respiration signal and the cardiac shock signal, the contribution of the respiration signal to the cardiac shock signal is filtered out through morphological filtering, and then the respiration signal is extracted. Since the respiration signal and the ballistocardiogram signal are in a superimposed relationship, the ballistocardiogram signal can be obtained by subtracting the respiration signal from the filtered sign signal.
In an alternative embodiment, referring to fig. 3, the step S30 includes steps S31 to S36, which are as follows:
s31: and performing time sequence bidirectional detection on the ballistocardiogram signal according to a preset ballistocardiogram signal template, determining a peak point coordinate of the ballistocardiogram signal, and obtaining a plurality of cardiac intervals according to the peak point coordinate.
In the embodiment of the application, a preset ballistocardiogram signal template is established according to morphological characteristics of ballistocardiogram signals, and then the ballistocardiogram signals are subjected to time sequence bidirectional detection according to the preset ballistocardiogram signal template to locate the J peak position. The position of the positioned J peak is further refined and calibrated through setting the refractory period and interval stability, and the final position of the J peak is obtained. The difference calculation of a plurality of front and back J peaks can obtain the corresponding value of the JJ interval, namely a plurality of cardiac intervals.
S32: obtaining a plurality of heart rates according to the plurality of cardiac intervals; acquiring heart rate data parameters corresponding to a plurality of heart rates, and taking the heart rate data parameters as a heart rate characteristic sequence; wherein the heart rate data parameters include a mode of the heart rate, a mean of the heart rate, a maximum of the heart rate, and a minimum of the heart rate.
In an embodiment of the present application, several heart rates may be calculated over several cardiac intervals. Wherein, the calculation formula of the heart rate per minute is as follows:
where HR represents heart rate and T _ BCG represents cardiac intervals.
After obtaining the plurality of heart rates, a mode, a mean, a maximum, and a minimum of the plurality of heart rates may be calculated, thereby obtaining a heart rate feature sequence. Specifically, taking the sleep time period of 5 hours corresponding to one night as an example, 5 minutes is one unit time period, and 5 hours is divided into 60 unit time periods. An average cardiac interval may be obtained every 5 minutes so that 60 average cardiac intervals may be obtained all night, and the respective heart rates are calculated for the 60 average cardiac intervals, thereby obtaining the mode, mean, maximum and minimum of the heart rates.
S33: calculating standard deviations of the cardiac intervals, and obtaining a plurality of power spectrums according to the standard deviations and a preset autoregressive spectrum estimation method; comparing the plurality of power spectrums with a preset frequency spectrum range, and determining the frequency spectrum range of each power spectrum; obtaining heart rate deceleration force according to the plurality of cardiac intervals; obtaining a heart rate variability characteristic sequence according to the average value of all power spectrums in the same frequency spectrum range and the heart rate deceleration force;
in the embodiment of the present application, according to a plurality of cardiac intervals, a corresponding standard deviation may be calculated, and the calculation formula is as follows:
wherein N represents the total number of cardiac intervals,
representing the mean of N cardiac intervals, T _ BCG i Representing the ith cardiac interval. Specifically, the mean is calculated for 60 average cardiac intervals, which in turn is calculated to obtain the standard deviation.
The autoregressive spectrum estimation method is an important spectrum estimation method, has the advantage of high frequency resolution, and improves the phenomena of spectral line splitting and frequency shift. And performing frequency domain analysis on each standard deviation by a preset autoregressive spectrum estimation method to obtain a plurality of power spectrums. According to a preset Frequency spectrum range, dividing a plurality of power spectrums into four Frequency bands including a High Frequency band (HF): 0.15 Hz-0.4 Hz, low Frequency band (LF): 0.04 Hz-0.15 Hz, very Low Frequency (Very Low Frequency, VLF): 0.003 Hz-0.04 Hz, ultra Low Frequency (ULF): <0.003Hz. And taking the power of the four frequency ranges as the HRV frequency domain characteristics.
According to a plurality of cardiac intervals, heart Rate Deceleration Capacity (DC) is obtained. Specifically, the heart rate deceleration force calculation formula is as follows:
DC＝[X(0)+X(1)-X(-1)-X(-2)]*1/4
wherein X (0) is the mean of the cardiac intervals of all the deceleration center points; x (1) is the average value of the first cardiac intervals closely adjacent to the right side of the deceleration central point; x (-1) is the mean of the first cardiac interval adjacent to the left of the deceleration center point; x (-2) is the mean of all cardiac intervals of the second adjacent left hand side of the deceleration center point.
The frequency band powers LF, HF, ULF, VLF, LF to HF ratios and DC are taken as the heart rate variability signature sequence.
S34: a number of breathing intervals of the breathing signal are determined according to a preset zero crossing detection method.
In the embodiment of the present application, a zero-crossing point detection method is adopted, in which a zero-crossing point of a respiratory signal in two adjacent rising edges is detected as a respiratory interval, that is, an end point of a respiratory signal is regarded as a start point of a next respiratory signal.
S35: obtaining a plurality of respiratory frequencies according to the plurality of respiratory intervals; acquiring breathing data parameters corresponding to a plurality of breathing frequencies, and taking the breathing data parameters as a breathing frequency characteristic sequence; the respiratory data parameters comprise a mode of respiratory frequency, a mean value of respiratory frequency, a maximum value of respiratory frequency, a minimum value of respiratory frequency, a mean value of respiratory frequency which is greater than or equal to a first preset threshold value and a mean value of respiratory frequency which is less than a second preset threshold value.
In the embodiment of the present application, the calculation formula of the breathing rate per minute is as follows:
where RR denotes respiratory frequency and T _ RR denotes respiratory interval.
After obtaining several respiratory frequencies, the mode RR1, the mean RR2, the maximum RR3, and the minimum RR4 of the respiratory frequencies may be calculated. And simultaneously comparing the plurality of respiratory frequencies with a first preset threshold and a second preset threshold so as to determine all respiratory frequencies which are greater than or equal to the first preset threshold and all respiratory frequencies which are less than the second preset threshold. Averaging all the respiratory frequencies greater than or equal to the first preset threshold value to obtain an average result RR5. Wherein the first preset threshold is 21 times/min, and RR5 is used to indicate that the breathing rate is too fast. Averaging all the respiratory frequencies smaller than a second preset threshold value to obtain an average result RR6. Wherein the first preset threshold is 10 times/min, and RR6 is used to indicate that the breathing rate is too slow. RR1, RR2, RR3, RR4, RR5, and RR6 are respiratory rate signature sequences.
S36: and obtaining a sleep characteristic sequence according to a plurality of cardiac intervals, a plurality of respiratory intervals and a preset heart-lung coupling method.
In an embodiment of the present application, the sleep signature sequence includes a time to sleep, a wake time, a total sleep duration, and an effective sleep duration (i.e., total sleep duration minus sleep wake duration). And calculating the sleep stages of the user under different time scales by adopting a preset cardiopulmonary coupling method by means of the cardiac interval and the respiratory interval so as to obtain the sleep characteristic sequence under the corresponding time scale.
In an alternative embodiment, referring to fig. 4, the step S40 includes steps S41 to S42, which are as follows:
s41: traversing each second time scale, and calculating mathematical expectations, normalized variances, maximum values and mode of the heart rate feature sequence, the heart rate variability feature sequence, the respiratory rate feature sequence and the sleep feature sequence at the current second time scale, wherein the mathematical expectations, the normalized variances, the maximum values and the mode are respectively used as a reference feature sequence, a fluctuation feature sequence, a most significant feature sequence and a mode feature sequence corresponding to the current second time scale;
s42: and repeatedly traversing each second time scale until determining the reference feature sequence, the fluctuation feature sequence, the most significant feature sequence and the mode feature sequence corresponding to the heart rate feature sequence, the heart rate variability feature sequence, the respiratory rate feature sequence and the sleep feature sequence at all the second time scales.
In the embodiment of the present application, for the cases where the second time scale is 3 consecutive days, 5 consecutive days, 7 consecutive days, 10 consecutive days, 14 consecutive days, and 30 consecutive days, respectively, the corresponding reference feature sequence, fluctuation feature sequence, maximum feature sequence, and mode feature sequence are calculated.
The second time scale is 3 consecutive days, and the reference feature sequence corresponding to HR1 in the heart rate feature sequence is taken as an example for explanation. At a preset first time scale, HR1 is known for each day, and the mathematical expectation of HR1 for 3 consecutive days is taken as the value for each day in the baseline signature sequence. For day 1 values in the baseline signature sequence, only day 1 HR1 is known, and there are no day 2 and no day 3 HR1, therefore day 1 values in the baseline signature sequence are 0. For day 2 values in the baseline signature sequence, only HR1 was known for day 1 and day 2, and there was no HR1 for day 3, so the day 2 value in the baseline signature sequence was 0. For the values of day 3 in the reference signature sequence, since HR1 is known for day 1, day 2 and day 3, the mathematical expectations of HR1 for day 1, day 2 and day 3 are taken as the value X1 for day 3 in the reference signature sequence, the mathematical expectations of HR1 for day 2, day 3 and day 4 are taken as the value X2 for day 4 in the reference signature sequence, and the mathematical expectations of HR1 for day 98, day 99 and day 100 are taken as the value X98 for day 100 in the reference signature sequence, thereby obtaining the reference signature sequence of HR1 for the heart rate signature sequence of consecutive 3 days.
For the solving processes of the heart rate feature sequence, the heart rate variability feature sequence, the respiratory rate feature sequence, and the reference feature sequence, the fluctuation feature sequence, the most significant feature sequence, and the mode feature sequence corresponding to the sleep feature sequence at different second time scales, the solving processes are consistent with the solving process of the above example, and are not repeated here.
In an optional embodiment, before the step S50, a step S501 is included, specifically as follows:
s501: and inputting second characteristic sample data under a plurality of different second time scales as input, and inputting breathing state sample data under a plurality of different second time scales as output to the breathing state detection model to be trained to obtain the trained breathing state detection model.
In the embodiment of the present application, the respiratory state detection model to be trained may be a deep learning network model or a machine learning model. Specifically, the machine learning model may be one of a Decision Tree model (DT), a Logistic Regression model (LR), a Support Vector Machine (SVM), a Random Forest model (RF), and an AdaBoost model. And acquiring second characteristic sample data and corresponding respiratory state sample data under a plurality of different second time scales, inputting the second characteristic sample data under the plurality of different second time scales as input, inputting the respiratory state sample data under the plurality of different second time scales as output to the respiratory state detection model to be trained, updating the network weight parameters of the respiratory state detection model to be trained, and acquiring the trained respiratory state detection model.
Example 2
The following is an example of the apparatus of the present application that can be used to perform the method of example 1 of the present application. For details which are not disclosed in the embodiments of the apparatus of the present application, please refer to the contents of the method in embodiment 1 of the present application.
Please refer to fig. 5, which shows a schematic structural diagram of a respiratory state detection apparatus provided in an embodiment of the present application. The respiratory state detection device 6 provided by the embodiment of the application comprises:
the sign signal acquiring module 61 is configured to acquire a sign signal of a user at a preset first time scale;
a signal extraction module 62, configured to extract a ballistocardiogram signal and a respiration signal from the sign signal;
a first feature sequence determining module 63, configured to determine a first feature sequence corresponding to the ballistocardiogram signal and the respiration signal; the first characteristic sequence comprises a heart rate characteristic sequence, a heart rate variability characteristic sequence, a respiratory rate characteristic sequence and a sleep characteristic sequence;
a second feature sequence determination module 64 for determining a second feature sequence of the heart rate feature sequence, the heart rate variability feature sequence, the respiratory rate feature sequence and the sleep feature sequence at a number of different second time scales; the second feature sequence comprises a reference feature sequence, a fluctuation feature sequence, a maximum feature sequence and a mode feature sequence;
a detection result obtaining module 65, configured to input the second feature sequences at a plurality of different second time scales to the trained respiratory state detection model, so as to obtain a respiratory state detection result.
Optionally, the signal extraction module includes:
the signal filtering unit is used for filtering the physical sign signals by adopting a filter to obtain filtered physical sign signals;
the first integral sign signal unit is used for carrying out corrosion-first and expansion-second opening operation on the filtered sign signals to obtain first sign signals;
the second body sign signal unit is used for performing expansion-first and corrosion-later closed operation on the filtered physical sign signals to obtain second physical sign signals;
a respiratory signal obtaining unit, configured to average the first body sign signal and the second body sign signal, and use an average result as a respiratory signal;
and the ballistocardiogram signal obtaining unit is used for subtracting the filtered sign signal from the respiratory signal to obtain a ballistocardiogram signal.
Optionally, the first feature sequence determining module includes:
the cardiac interval acquisition unit is used for carrying out time sequence bidirectional detection on the ballistocardiogram signal according to a preset ballistocardiogram signal template, determining the coordinates of peak points of the ballistocardiogram signal and acquiring a plurality of cardiac intervals according to the coordinates of the peak points;
the heart rate characteristic sequence obtaining unit is used for obtaining a plurality of heart rates according to the plurality of cardiac intervals; acquiring heart rate data parameters corresponding to a plurality of heart rates, and taking the heart rate data parameters as a heart rate characteristic sequence; the heart rate data parameters comprise the mode of the heart rate, the mean value of the heart rate, the maximum value of the heart rate and the minimum value of the heart rate;
the heart rate variability characteristic sequence obtaining unit is used for calculating standard deviations of the plurality of cardiac intervals and obtaining a plurality of power spectrums according to the standard deviations and a preset autoregressive spectrum estimation method; comparing the plurality of power spectrums with a preset frequency spectrum range, and determining the frequency spectrum range of each power spectrum; obtaining heart rate deceleration force according to the plurality of cardiac intervals; obtaining a heart rate variability characteristic sequence according to the average value of all power spectrums in the same frequency spectrum range and the heart rate deceleration force;
the breathing interval determining unit is used for determining a plurality of breathing intervals of the breathing signal according to a preset zero crossing point detection method;
the respiratory frequency characteristic sequence obtaining unit is used for obtaining a plurality of respiratory frequencies according to the plurality of respiratory intervals; acquiring breathing data parameters corresponding to a plurality of breathing frequencies, and taking the breathing data parameters as a breathing frequency characteristic sequence; the respiratory data parameters comprise a mode of respiratory frequency, a mean value of respiratory frequency, a maximum value of respiratory frequency, a minimum value of respiratory frequency, a mean value of respiratory frequency which is greater than or equal to a first preset threshold value and a mean value of respiratory frequency which is less than a second preset threshold value;
and the sleep characteristic sequence obtaining unit is used for obtaining the sleep characteristic sequence according to the plurality of cardiac intervals, the plurality of respiratory intervals and a preset heart-lung coupling method.
Optionally, the second feature sequence determining module includes:
and the feature sequence calculating unit is used for traversing each second time scale, calculating mathematical expectations, normalized variances, maximum values and modes of the heart rate feature sequence, the heart rate variability feature sequence, the respiratory frequency feature sequence and the sleep feature sequence under the current second time scale, and respectively using the mathematical expectations, the normalized variances, the maximum values and the modes as a reference feature sequence, a fluctuation feature sequence, a maximum value feature sequence and a mode feature sequence corresponding to the current second time scale.
And the characteristic sequence determining unit is used for repeatedly traversing each second time scale until determining the heart rate characteristic sequence, the heart rate variability characteristic sequence, the respiratory rate characteristic sequence and the reference characteristic sequence, the fluctuation characteristic sequence, the most significant characteristic sequence and the mode characteristic sequence corresponding to the sleep characteristic sequence in all second time scales.
By applying the embodiment of the application, the sign signals of the user under the preset first time scale are obtained; extracting a ballistocardiogram signal and a respiration signal from the sign signal; determining a first feature sequence corresponding to the ballistocardiogram signal and the respiration signal; the first feature sequence comprises a heart rate feature sequence, a heart rate variability feature sequence, a respiratory rate feature sequence and a sleep feature sequence; determining a second sequence of features of the heart rate feature sequence, the heart rate variability feature sequence, the respiratory rate feature sequence, and the sleep feature sequence at a number of different second time scales; the second feature sequence comprises a reference feature sequence, a fluctuation feature sequence, a maximum feature sequence and a mode feature sequence; and inputting the second characteristic sequences under a plurality of different second time scales into a trained respiratory state detection model to obtain a respiratory state detection result. This application acquires user's sign signal through undisturbed sensor acquisition equipment, can monitor for a long time to the cost is reduced has improved the degree of accuracy that respiratory state detected.
Example 3
The following is an embodiment of the apparatus of the present application, which may be used to perform the method of embodiment 1 of the present application. For details which are not disclosed in the device example of the present application, reference is made to the content of the method in example 1 of the present application.
Referring to fig. 6, the present application further provides an electronic device 300, where the electronic device may be embodied as a computer, a mobile phone, a tablet computer, an interactive tablet, and the like, and in an exemplary embodiment of the present application, the electronic device 300 is an interactive tablet, and the interactive tablet may include: at least one processor 301, at least one memory 302, at least one display, at least one network interface 303, a user interface 304, and at least one communication bus 305.
The user interface 304 is mainly used for providing an input interface for a user to obtain data input by the user. Optionally, the user interface may also include a standard wired interface, a wireless interface.
The network interface 303 may optionally include a standard wired interface or a wireless interface (e.g., WI-FI interface).
Wherein a communication bus 305 is used to enable the connection communication between these components.
The Memory 302 may include a Random Access Memory (RAM) or a Read-Only Memory (Read-Only Memory). Optionally, the memory includes a non-transitory computer-readable medium. The memory may be used to store an instruction, a program, code, a set of codes, or a set of instructions. The memory may include a stored program area and a stored data area, wherein the stored program area may store instructions for implementing an operating system, instructions for at least one function (such as a touch function, a sound playing function, an image playing function, etc.), instructions for implementing the various method embodiments described above, and the like; the storage data area may store data and the like referred to in the above respective method embodiments. The memory may alternatively be at least one memory device located remotely from the aforementioned processor. The memory, which is a type of computer storage medium, may include an operating system, a network communication module, a user interface module, and an operating application.
The processor may be configured to invoke an application program of the video resolution adjustment method stored in the memory, and specifically execute the method steps in embodiment 1 shown above, and the specific execution process may refer to the specific description shown in embodiment 1, which is not described herein again.
Example 4
The present application further provides a computer-readable storage medium, on which a computer program is stored, where the instructions are suitable for being loaded by a processor and executing the method steps of embodiment 1, and specific execution processes may refer to specific descriptions shown in the embodiments and are not described herein again. The device where the storage medium is located can be an electronic device such as a personal computer, a notebook computer, a smart phone and a tablet computer.
For the apparatus embodiment, since it basically corresponds to the method embodiment, reference may be made to the partial description of the method embodiment for relevant points. The above-described device embodiments are merely illustrative, and components illustrated as separate components may or may not be physically separate, and components shown as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the modules can be selected according to actual needs to achieve the purpose of the scheme of the application. One of ordinary skill in the art can understand and implement without inventive effort.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks. These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
In a typical configuration, a computing device includes one or more processors (CPUs), input/output interfaces, network interfaces, and memory.
The memory may include forms of volatile memory in a computer readable medium, random Access Memory (RAM) and/or non-volatile memory, such as Read Only Memory (ROM) or flash memory (flash RAM). The memory is an example of a computer-readable medium.
Computer-readable media, including both non-transitory and non-transitory, removable and non-removable media, may implement information storage by any method or technology. The information may be computer readable instructions, data structures, modules of a program, or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), static Random Access Memory (SRAM), dynamic Random Access Memory (DRAM), other types of Random Access Memory (RAM), read Only Memory (ROM), electrically Erasable Programmable Read Only Memory (EEPROM), flash memory or other memory technology, compact disc read only memory (CD-ROM), digital Versatile Discs (DVD) or other optical storage, magnetic cassettes, magnetic tape magnetic disk storage or other magnetic storage devices, or any other non-transmission medium that can be used to store information that can be accessed by a computing device. As defined herein, a computer readable medium does not include a transitory computer readable medium such as a modulated data signal and a carrier wave.
It should also be noted that the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrases "comprising a," "8230," "8230," or "comprising" does not exclude the presence of additional identical elements in the process, method, article, or apparatus comprising the element.
The above are merely examples of the present application and are not intended to limit the present application. Various modifications and changes may occur to those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application should be included in the scope of the claims of the present application.
Claims (10)
1. A method of respiratory state detection, the method comprising the steps of:
acquiring a sign signal of a user at a preset first time scale;
extracting a ballistocardiogram signal and a respiration signal from the sign signal;
determining a first feature sequence corresponding to the ballistocardiogram signal and the respiration signal; the first feature sequence comprises a heart rate feature sequence, a heart rate variability feature sequence, a respiratory rate feature sequence and a sleep feature sequence;
determining a second sequence of features of the heart rate feature sequence, the heart rate variability feature sequence, the respiratory rate feature sequence, and the sleep feature sequence at a number of different second time scales; the second feature sequence comprises a reference feature sequence, a fluctuation feature sequence, a most-valued feature sequence and a mode feature sequence;
and inputting the second characteristic sequences under a plurality of different second time scales into a trained respiratory state detection model to obtain a respiratory state detection result.
2. The respiratory state detection method according to claim 1, wherein:
the step of extracting ballistocardiogram signals and respiration signals from the sign signals comprises:
filtering the sign signals by using a filter to obtain filtered sign signals;
carrying out corrosion-first and expansion-second opening operation on the filtered sign signals to obtain first sign signals;
performing closed operation of expansion and corrosion on the filtered sign signal to obtain a second sign signal;
averaging the first body sign signal and the second body sign signal, and taking an average result as a respiration signal;
and subtracting the filtered sign signal from the respiration signal to obtain a ballistocardiogram signal.
3. The respiratory state detection method according to claim 1, wherein:
the determining of the first feature sequence corresponding to the ballistocardiogram signal and the respiration signal; the steps of the first signature sequence including a heart rate signature sequence, a heart rate variability signature sequence, a respiratory rate signature sequence, and a sleep signature sequence include:
performing time sequence bidirectional detection on the ballistocardiogram signal according to a preset ballistocardiogram signal template, determining a peak point coordinate of the ballistocardiogram signal, and obtaining a plurality of cardiac intervals according to the peak point coordinate;
obtaining a plurality of heart rates according to the plurality of cardiac intervals; acquiring heart rate data parameters corresponding to a plurality of heart rates, and taking the heart rate data parameters as a heart rate characteristic sequence; the heart rate data parameters comprise the mode of the heart rate, the mean value of the heart rate, the maximum value of the heart rate and the minimum value of the heart rate;
calculating standard deviations of the plurality of cardiac intervals, and obtaining a plurality of power spectrums according to the standard deviations and a preset autoregressive spectrum estimation method; comparing the plurality of power spectrums with a preset frequency spectrum range, and determining the frequency spectrum range of each power spectrum; obtaining heart rate deceleration force according to the plurality of cardiac intervals; obtaining a heart rate variability characteristic sequence according to the average value of all power spectrums in the same frequency spectrum range and the heart rate deceleration force;
determining a plurality of breathing intervals of the breathing signal according to a preset zero crossing point detection method;
obtaining a plurality of respiratory frequencies according to the plurality of respiratory intervals; acquiring breathing data parameters corresponding to a plurality of breathing frequencies, and taking the breathing data parameters as a breathing frequency characteristic sequence; the respiratory data parameters comprise a mode of respiratory frequency, a mean value of respiratory frequency, a maximum value of respiratory frequency, a minimum value of respiratory frequency, a mean value of respiratory frequency which is greater than or equal to a first preset threshold value and a mean value of respiratory frequency which is less than a second preset threshold value;
and obtaining a sleep characteristic sequence according to a plurality of cardiac intervals, a plurality of respiratory intervals and a preset heart-lung coupling method.
4. The respiratory state detection method of claim 1, wherein:
said determining a second signature sequence of said heart rate signature sequence, said heart rate variability signature sequence, said respiratory rate signature sequence, and said sleep signature sequence at a number of different second time scales; the second feature sequence comprises a reference feature sequence, a fluctuation feature sequence, a maximum feature sequence and a mode feature sequence, and the method comprises the following steps:
traversing each second time scale, and calculating mathematical expectations, normalized variances, maximum values and mode of the heart rate feature sequence, the heart rate variability feature sequence, the respiratory rate feature sequence and the sleep feature sequence at the current second time scale, wherein the mathematical expectations, the normalized variances, the maximum values and the mode are respectively used as a reference feature sequence, a fluctuation feature sequence, a most significant feature sequence and a mode feature sequence corresponding to the current second time scale;
and repeatedly traversing each second time scale until determining the reference characteristic sequence, the fluctuation characteristic sequence, the most significant characteristic sequence and the mode characteristic sequence corresponding to the heart rate characteristic sequence, the heart rate variability characteristic sequence, the respiratory rate characteristic sequence and the sleep characteristic sequence at all the second time scales.
5. The respiratory state detection method of claim 1, wherein:
before the step of inputting the second feature sequences under a plurality of different second time scales into the trained respiratory state detection model to obtain the respiratory state detection result, the method includes:
and inputting second characteristic sample data under a plurality of different second time scales as input, and inputting the breathing state sample data under a plurality of different second time scales as output to the breathing state detection model to be trained to obtain the trained breathing state detection model.
6. A respiratory state detection device, comprising:
the sign signal acquisition module is used for acquiring a sign signal of a user at a preset first time scale;
the signal extraction module is used for extracting a ballistocardiogram signal and a respiratory signal from the sign signal;
a first signature sequence determination module, configured to determine a first signature sequence corresponding to the ballistocardiogram signal and the respiration signal; the first feature sequence comprises a heart rate feature sequence, a heart rate variability feature sequence, a respiratory rate feature sequence and a sleep feature sequence;
a second feature sequence determination module for determining a second feature sequence of the heart rate feature sequence, the heart rate variability feature sequence, the respiratory rate feature sequence, and the sleep feature sequence at a number of different second time scales; the second feature sequence comprises a reference feature sequence, a fluctuation feature sequence, a maximum feature sequence and a mode feature sequence;
and the detection result obtaining module is used for inputting the second characteristic sequences under a plurality of different second time scales into the trained respiratory state detection model to obtain a respiratory state detection result.
7. The respiratory state detection device of claim 6, wherein the signal extraction module comprises:
the signal filtering unit is used for filtering the sign signals by adopting a filter to obtain filtered sign signals;
the first integral sign signal unit is used for carrying out corrosion-first and expansion-second opening operation on the filtered sign signals to obtain first sign signals;
the second body sign signal unit is used for performing closed operation of expansion and corrosion on the filtered sign signals to obtain second sign signals;
a respiratory signal obtaining unit, configured to average the first volume sign signal and the second volume sign signal, and use an average result as a respiratory signal;
and the ballistocardiogram signal obtaining unit is used for subtracting the filtered sign signal from the respiratory signal to obtain a ballistocardiogram signal.
8. The respiratory state detection device of claim 6, wherein the first signature sequence determination module comprises:
the cardiac interval acquisition unit is used for carrying out time sequence bidirectional detection on the ballistocardiogram signal according to a preset ballistocardiogram signal template, determining the coordinates of peak points of the ballistocardiogram signal and acquiring a plurality of cardiac intervals according to the coordinates of the peak points;
the heart rate characteristic sequence obtaining unit is used for obtaining a plurality of heart rates according to the plurality of cardiac intervals; acquiring heart rate data parameters corresponding to a plurality of heart rates, and taking the heart rate data parameters as a heart rate characteristic sequence; the heart rate data parameters comprise mode of heart rate, mean value of heart rate, maximum value of heart rate and minimum value of heart rate;
the heart rate variability characteristic sequence obtaining unit is used for calculating the standard deviation of the plurality of cardiac intervals and obtaining a plurality of power spectrums according to the standard deviation and a preset autoregressive spectrum estimation method; comparing the plurality of power spectrums with a preset frequency spectrum range, and determining the frequency spectrum range of each power spectrum; obtaining heart rate deceleration force according to the plurality of cardiac intervals; obtaining a heart rate variability characteristic sequence according to the average value of all power spectrums in the same frequency spectrum range and the heart rate deceleration force;
the breathing interval determining unit is used for determining a plurality of breathing intervals of the breathing signal according to a preset zero crossing point detection method;
the respiratory frequency characteristic sequence obtaining unit is used for obtaining a plurality of respiratory frequencies according to the plurality of respiratory intervals; acquiring a plurality of respiratory data parameters corresponding to respiratory frequency, and taking the respiratory data parameters as a respiratory frequency characteristic sequence; the respiratory data parameters comprise a mode of respiratory frequency, a mean value of respiratory frequency, a maximum value of respiratory frequency, a minimum value of respiratory frequency, a mean value of respiratory frequency which is greater than or equal to a first preset threshold value and a mean value of respiratory frequency which is less than a second preset threshold value;
and the sleep characteristic sequence obtaining unit is used for obtaining the sleep characteristic sequence according to the plurality of cardiac intervals, the plurality of respiratory intervals and a preset heart-lung coupling method.
9. A computer device, comprising: processor, memory and computer program stored in the memory and executable on the processor, characterized in that the steps of the method according to any of claims 1 to 5 are implemented when the processor executes the computer program.
10. A computer-readable storage medium, in which a computer program is stored which, when being executed by a processor, carries out the steps of the method according to any one of claims 1 to 5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211004088.5A CN115486833B (en) | 2022-08-22 | 2022-08-22 | Respiratory state detection method, respiratory state detection device, computer equipment and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211004088.5A CN115486833B (en) | 2022-08-22 | 2022-08-22 | Respiratory state detection method, respiratory state detection device, computer equipment and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115486833A true CN115486833A (en) | 2022-12-20 |
| CN115486833B CN115486833B (en) | 2023-06-06 |
Family
ID=84466129
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211004088.5A Active CN115486833B (en) | 2022-08-22 | 2022-08-22 | Respiratory state detection method, respiratory state detection device, computer equipment and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115486833B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116649919A (en) * | 2023-07-28 | 2023-08-29 | 华南师范大学 | Sleep stage prediction method, device and equipment based on physiological signals |
| CN117951630A (en) * | 2024-03-26 | 2024-04-30 | 大连云间来客科技有限公司 | Abnormal breathing pattern recognition method for patient with severe respiratory disease |
Citations (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20050267362A1 (en) * | 2004-05-17 | 2005-12-01 | Beth Israel Deaconess Medical Center | Assessment of sleep quality and sleep disordered breathing based on cardiopulmonary coupling |
| US20070073169A1 (en) * | 2005-09-16 | 2007-03-29 | Averina Viktoria A | System and method for generating a trend parameter based on respiration rate distribution |
| CN106037671A (en) * | 2016-07-11 | 2016-10-26 | 西北工业大学 | Method and system for apnea event detection based on BCG signal |
| CN108056769A (en) * | 2017-11-14 | 2018-05-22 | 深圳市大耳马科技有限公司 | A kind of vital sign parameter signals analysis and processing method, device and vital sign monitoring device |
| WO2019075520A1 (en) * | 2017-10-18 | 2019-04-25 | Medibio Limited | Breathing state indicator |
| CN109674474A (en) * | 2018-11-30 | 2019-04-26 | 深圳和而泰智能控制股份有限公司 | Sleep apnea recognition methods, equipment and computer-readable medium |
| EP3492945A1 (en) * | 2017-12-01 | 2019-06-05 | Origin Wireless, Inc. | Method, apparatus, and system for periodic motion detection and monitoring |
| CN111035367A (en) * | 2019-12-31 | 2020-04-21 | 华南师范大学 | Signal detection method and system for judging sleep apnea |
| US20200253547A1 (en) * | 2017-10-31 | 2020-08-13 | Apple Inc. | Monitoring System for Assessing Control of a Disease State |
| CN111803038A (en) * | 2020-09-03 | 2020-10-23 | 平安国际智慧城市科技股份有限公司 | Method, system, terminal and storage medium for monitoring sleep apnea event |
| CN112155560A (en) * | 2020-10-15 | 2021-01-01 | 国微集团(深圳)有限公司 | Apnea detection method and system based on real-time cardiac shock signal |
| US10945659B1 (en) * | 2015-03-16 | 2021-03-16 | Dp Technologies, Inc. | Dual sleep monitor |
| CN114010185A (en) * | 2021-11-17 | 2022-02-08 | 上海跃扬医疗科技有限公司 | Respiratory disorder monitoring method, device, equipment and medium |
| CN114041780A (en) * | 2021-12-08 | 2022-02-15 | 中国科学院计算技术研究所 | Method for monitoring respiration based on data acquired by inertial sensor |
| CN114376564A (en) * | 2021-12-29 | 2022-04-22 | 华南理工大学 | Sleep staging method, system, device and medium based on cardiac shock signal |
| JP2022063926A (en) * | 2020-10-13 | 2022-04-25 | ヘルスセンシング株式会社 | Sleep state estimation system |
| CN114449943A (en) * | 2019-07-25 | 2022-05-06 | 启迪医疗仪器公司 | Sleep detection for Sleep Disordered Breathing (SDB) care |
| CN114699040A (en) * | 2022-02-21 | 2022-07-05 | 华南师范大学 | Waking-up detection method, device, equipment and storage medium based on physiological signals |
| CN114711725A (en) * | 2022-04-15 | 2022-07-08 | 华南师范大学 | Sleep apnea detection method and device based on double attention mechanism |
-
2022
- 2022-08-22 CN CN202211004088.5A patent/CN115486833B/en active Active
Patent Citations (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20050267362A1 (en) * | 2004-05-17 | 2005-12-01 | Beth Israel Deaconess Medical Center | Assessment of sleep quality and sleep disordered breathing based on cardiopulmonary coupling |
| US20070073169A1 (en) * | 2005-09-16 | 2007-03-29 | Averina Viktoria A | System and method for generating a trend parameter based on respiration rate distribution |
| US10945659B1 (en) * | 2015-03-16 | 2021-03-16 | Dp Technologies, Inc. | Dual sleep monitor |
| CN106037671A (en) * | 2016-07-11 | 2016-10-26 | 西北工业大学 | Method and system for apnea event detection based on BCG signal |
| WO2019075520A1 (en) * | 2017-10-18 | 2019-04-25 | Medibio Limited | Breathing state indicator |
| US20200253547A1 (en) * | 2017-10-31 | 2020-08-13 | Apple Inc. | Monitoring System for Assessing Control of a Disease State |
| CN108056769A (en) * | 2017-11-14 | 2018-05-22 | 深圳市大耳马科技有限公司 | A kind of vital sign parameter signals analysis and processing method, device and vital sign monitoring device |
| EP3492945A1 (en) * | 2017-12-01 | 2019-06-05 | Origin Wireless, Inc. | Method, apparatus, and system for periodic motion detection and monitoring |
| CN109674474A (en) * | 2018-11-30 | 2019-04-26 | 深圳和而泰智能控制股份有限公司 | Sleep apnea recognition methods, equipment and computer-readable medium |
| CN114449943A (en) * | 2019-07-25 | 2022-05-06 | 启迪医疗仪器公司 | Sleep detection for Sleep Disordered Breathing (SDB) care |
| CN111035367A (en) * | 2019-12-31 | 2020-04-21 | 华南师范大学 | Signal detection method and system for judging sleep apnea |
| CN111803038A (en) * | 2020-09-03 | 2020-10-23 | 平安国际智慧城市科技股份有限公司 | Method, system, terminal and storage medium for monitoring sleep apnea event |
| JP2022063926A (en) * | 2020-10-13 | 2022-04-25 | ヘルスセンシング株式会社 | Sleep state estimation system |
| CN112155560A (en) * | 2020-10-15 | 2021-01-01 | 国微集团(深圳)有限公司 | Apnea detection method and system based on real-time cardiac shock signal |
| CN114010185A (en) * | 2021-11-17 | 2022-02-08 | 上海跃扬医疗科技有限公司 | Respiratory disorder monitoring method, device, equipment and medium |
| CN114041780A (en) * | 2021-12-08 | 2022-02-15 | 中国科学院计算技术研究所 | Method for monitoring respiration based on data acquired by inertial sensor |
| CN114376564A (en) * | 2021-12-29 | 2022-04-22 | 华南理工大学 | Sleep staging method, system, device and medium based on cardiac shock signal |
| CN114699040A (en) * | 2022-02-21 | 2022-07-05 | 华南师范大学 | Waking-up detection method, device, equipment and storage medium based on physiological signals |
| CN114711725A (en) * | 2022-04-15 | 2022-07-08 | 华南师范大学 | Sleep apnea detection method and device based on double attention mechanism |
Non-Patent Citations (2)
| Title |
|---|
| 聂雨亭: "基于BCG和呼吸信号的睡眠分期方法研究" * |
| 陈培敏: "一种非接触式基于心冲击图的睡眠监测装置" * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116649919A (en) * | 2023-07-28 | 2023-08-29 | 华南师范大学 | Sleep stage prediction method, device and equipment based on physiological signals |
| CN116649919B (en) * | 2023-07-28 | 2023-11-14 | 华南师范大学 | Sleep stage prediction method, device and equipment based on physiological signals |
| CN117951630A (en) * | 2024-03-26 | 2024-04-30 | 大连云间来客科技有限公司 | Abnormal breathing pattern recognition method for patient with severe respiratory disease |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115486833B (en) | 2023-06-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN115486833B (en) | Respiratory state detection method, respiratory state detection device, computer equipment and storage medium | |
| US10646168B2 (en) | Drowsiness onset detection | |
| CN105813558B (en) | System and method for determining sleep stages based on sleep cycle | |
| US20170007137A1 (en) | Method of estimating blood pressure based on image | |
| CN106470601B (en) | System and method for determining and displaying sleep restoration levels | |
| US9986951B1 (en) | Low-distortion ECG denoising | |
| JP6515670B2 (en) | Sleep depth estimation device, sleep depth estimation method, and program | |
| CN105705079A (en) | Device, method and system for processing a physiological signal | |
| CN104720746A (en) | Sleeping stage determination method and system | |
| JP2012212362A (en) | Information processing apparatus, information processing method, and program | |
| CN110236491A (en) | A kind of sleep stage monitoring method | |
| CN106681961B (en) | Signal feature extraction method and device | |
| JP2018064900A (en) | Drowsiness determination device and drowsiness determination method, and driving support apparatus | |
| JP6468637B2 (en) | Sleepiness estimation apparatus and sleepiness estimation program | |
| JP6519344B2 (en) | Heartbeat interval specifying program, heart beat interval specifying device, and heart beat interval specifying method | |
| CN115067930B (en) | Breathing state early warning method and device, computer equipment and storage medium | |
| CN107080526B (en) | Biological signal processing method and apparatus | |
| JP6518570B2 (en) | Sleep stage estimation device, method and program | |
| CN114027825B (en) | Respiratory signal acquisition method and device and computer equipment | |
| CN115381438A (en) | Method and device for reconstructing vital sign signals, computer equipment and storage medium | |
| CN113616217B (en) | Method and device for generating baseline drift curve | |
| JP6329753B2 (en) | Information processing program, information processing apparatus, information processing system, and sound determination method | |
| CN112494008A (en) | Respiration rate measuring method and device based on PPG signal | |
| JP2017124014A (en) | Drowsiness detection program, drowsiness detection method, and drowsiness detection device | |
| US20230397890A1 (en) | Fatigue level estimation apparatus, fatigue level estimation method, and computer-readable recording medium |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |