CN115568844A - Apnea detection method and device, electronic equipment and storage medium - Google Patents
Apnea detection method and device, electronic equipment and storage medium Download PDFInfo
- Publication number
- CN115568844A CN115568844A CN202110761417.XA CN202110761417A CN115568844A CN 115568844 A CN115568844 A CN 115568844A CN 202110761417 A CN202110761417 A CN 202110761417A CN 115568844 A CN115568844 A CN 115568844A
- Authority
- CN
- China
- Prior art keywords
- peak
- apnea
- amplitude
- time
- condition
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/08—Detecting, measuring or recording devices for evaluating the respiratory organs
- A61B5/0826—Detecting or evaluating apnoea events
-
- 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/7242—Details of waveform analysis using integration
Abstract
The embodiment of the invention provides an apnea detection method and device, electronic equipment and a storage medium, and relates to the field of respiration detection. According to the method, the time domain characteristics of the acquired cardiac shock signals are extracted, the frequency domain characteristics of the cardiac shock signals are extracted under the condition that the time domain characteristics of the cardiac shock signals accord with a first apnea condition, and the apnea event is determined to exist under the condition that the frequency domain characteristics accord with a second apnea condition. Therefore, the time domain characteristic and the frequency domain characteristic of the cardiac shock signal are combined, the apnea event is determined to exist under the condition that the apnea is judged in both the time domain and the frequency domain, the possibility of misjudgment is effectively reduced, and the detection accuracy of the apnea event is improved; the whole detection process does not need to be artificially distinguished, the acquisition of the cardiac shock signal is more convenient, and the operation simplicity of the detection of the apnea event is improved.
Description
Technical Field
The invention relates to the field of breath detection, in particular to an apnea detection method, an apnea detection device, electronic equipment and a storage medium.
Background
Apnea refers to spontaneous cessation of breathing, often temporary or self-limiting. The apnea usually occurs in the sleep process and is not easy to detect, and if a long-term apnea phenomenon is not found, effective treatment cannot be obtained, a series of diseases can occur, so that the detection of the apnea is important for human health, and the apnea event is also an important index for sleep apnea monitoring.
The sleep state is monitored by the polysomnography instrument commonly used clinically, and the sleep state needs to be judged by a professional doctor, namely, the monitoring and the judgment of apnea can be completed only by special wearable equipment recording and manual interpretation, so that the operation is complex, and the sleep of a user can be influenced.
Disclosure of Invention
In view of the above, an object of the present invention is to provide an apnea detecting method, an apnea detecting device, an electronic device and a storage medium, so as to solve the problems of low accuracy and complex operation in detecting an apnea event in the prior art.
In order to achieve the above purpose, the embodiment of the present invention adopts the following technical solutions:
in a first aspect, the present invention provides a method of apnea detection, the method comprising:
extracting time domain characteristics of the acquired heart impact signals;
under the condition that the time domain features of the cardiac shock signals accord with a first apnea condition, extracting the frequency domain features of the cardiac shock signals;
determining that an apnea event is present if the frequency domain feature meets a second apnea condition.
In an alternative embodiment, the time-domain feature includes a magnitude of each peak in the ballistocardiogram signal, and the method further includes, before the step of extracting the frequency-domain feature of the ballistocardiogram signal in the case that the time-domain feature of the ballistocardiogram signal meets the first apnea condition, the step of:
determining a plurality of continuous target wave crests according to the time domain characteristics of the ballistocardiogram signal; wherein, in the plurality of continuous target peaks, the ratio of the amplitude of the former peak of the first target peak to the amplitude of the first target peak is greater than a preset ratio, the ratio of the amplitude of the latter peak of the last target peak to the amplitude of the last target peak is greater than the preset ratio, and the ratio of the larger amplitude to the smaller amplitude in the amplitudes of any two adjacent target peaks is smaller than or equal to the preset ratio;
and if the duration of the continuous target wave crests is greater than or equal to the preset time, judging that the time domain characteristics of the cardiac shock signals meet a first apnea condition.
In an alternative embodiment, the determining a plurality of continuous target peaks according to the time-domain feature of the ballistocardiogram signal includes:
searching peaks of the cardiac shock signal, and determining a first peak and a second peak which are adjacent to each other, wherein the time corresponding to the first peak is less than the time corresponding to the second peak;
if the ratio between the amplitude of the first peak and the amplitude of the second peak is greater than the preset ratio, searching for at least one continuous third peak after the moment corresponding to the second peak until the ratio between the amplitude of the next peak of the third peak and the amplitude of the third peak is greater than the preset ratio; wherein, in the amplitude of each third peak and the adjacent previous peak, the ratio of the larger amplitude to the smaller amplitude is less than or equal to the preset ratio;
and determining the second wave crest and all the third wave crests as target wave crests to obtain a plurality of continuous target wave crests.
In an alternative embodiment, the duration is the difference between the time at which the first target peak is located and the time at which the last target peak is located.
In an optional embodiment, the extracting the frequency-domain feature of the ballistocardiograph signal includes:
determining a time period during which a time domain feature of the ballistocardiographic signal meets a first apnea condition;
calculating a frequency spectrum of the ballistocardiogram signal over the time period; the frequency domain features include the magnitude of each peak in the frequency spectrum.
In an alternative embodiment, the frequency spectrum includes a frequency spectrum of a heartbeat frequency band and a frequency spectrum of a respiratory frequency band, and the method further includes, before the step of determining that an apnea event exists in the case that the frequency domain feature meets the second apnea condition:
calculating a first average value of the amplitudes of all peaks corresponding to the heartbeat frequency band and a second average value of the amplitudes of all peaks corresponding to the respiratory frequency band;
and if the amplitude of at least one peak in the heartbeat frequency band is greater than a first preset multiple of the first average value, and the amplitude of the peak in the respiratory frequency band is not greater than a second preset multiple of the second average value, judging that the frequency domain characteristic accords with a second apnea condition.
In an alternative embodiment, the calculating the frequency spectrum of the ballistocardiographic signal in the time period includes:
and carrying out Fourier transform on the ballistocardiogram signal in the time period to obtain a frequency spectrum of the ballistocardiogram signal in the time period.
In a second aspect, the present invention provides an apnea detection apparatus, said apparatus comprising:
the time domain characteristic extraction module is used for extracting the time domain characteristics of the acquired heart attack signals;
the frequency domain feature extraction module is used for extracting the frequency domain feature of the cardiac shock signal under the condition that the time domain feature of the cardiac shock signal meets a first apnea condition;
a breath detection module to determine that an apnea event is present if the frequency domain feature meets a second apnea condition.
In a third aspect, the invention provides an electronic device comprising a processor, a memory and a computer program stored on the memory and executable on the processor, the computer program, when executed by the processor, implementing the steps of the apnea detection method according to any of the preceding embodiments.
In a fourth aspect, the invention provides a computer-readable storage medium having stored thereon a computer program which, when executed by a processor, carries out the steps of the apnea detection method of any one of the preceding embodiments.
According to the apnea detection method, the apnea detection device, the electronic equipment and the storage medium provided by the embodiment of the invention, the time domain characteristics of the acquired cardiac shock signals are extracted, the frequency domain characteristics of the cardiac shock signals are extracted under the condition that the time domain characteristics of the cardiac shock signals meet a first apnea condition, and the apnea event is determined to exist under the condition that the frequency domain characteristics meet a second apnea condition. Therefore, the time domain characteristic and the frequency domain characteristic of the cardiac shock signal are combined, the apnea event is determined to exist under the condition that the apnea is judged in both the time domain and the frequency domain, the possibility of misjudgment is effectively reduced, and the detection accuracy of the apnea event is improved; the whole detection process does not need to be artificially distinguished, the acquisition of the cardiac shock signal is more convenient, and the operation simplicity of the detection of the apnea event is improved.
In order to make the aforementioned and other objects, features and advantages of the present invention comprehensible, preferred embodiments accompanied with figures are described in detail below.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings required in the embodiments will be briefly described below, it should be understood that the following drawings only illustrate some embodiments of the present invention and therefore should not be considered as limiting the scope, and those skilled in the art can also obtain other related drawings based on the drawings without inventive efforts.
Fig. 1 is a schematic flow chart of an apnea detection method provided by an embodiment of the present invention;
FIG. 2 is a schematic diagram showing a comparison of a respiratory waveform of a normal respiratory rhythm with a respiratory waveform of an apnea;
fig. 3 is a schematic flow chart of an apnea detecting method according to an embodiment of the present invention;
fig. 4 is a schematic flow chart of an apnea detecting method according to an embodiment of the present invention;
fig. 5 shows a schematic flow chart of an apnea detection method provided by the embodiment of the invention;
fig. 6 is a schematic flow chart of an apnea detecting method according to an embodiment of the present invention;
FIG. 7 illustrates a functional block diagram of an apnea detection apparatus provided by embodiments of the present invention;
FIG. 8 is another functional block diagram of an apnea detection apparatus provided in accordance with an embodiment of the present invention;
fig. 9 is a schematic diagram illustrating a hardware structure of an electronic device according to an embodiment of the present invention.
An icon: 100-an electronic device; 700-apnea detecting means; 111-a processor; 112-a storage medium; 113-a memory; 114-input-output interface; 115-wired or wireless network interface; 116-a power supply; 1121 — operating system; 1122-data; 1123-applications; 710-a time domain feature extraction module; 720-frequency domain feature extraction module; 730-a breath detection module; 740 — a first decision module; 750-a second determination module.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. The components of embodiments of the present invention generally described and illustrated in the figures herein may be arranged and designed in a wide variety of different configurations.
Thus, the following detailed description of the embodiments of the present invention, presented in the figures, is not intended to limit the scope of the invention, as claimed, but is merely representative of selected embodiments of the invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments of the present invention without making any creative effort, shall fall within the protection scope of the present invention.
It is noted that relational terms such as "first" and "second," and the like, may be used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions. Also, 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 like elements in a process, method, article, or apparatus that comprises the element.
In one embodiment, as shown in fig. 1, an apnea detecting method is provided, and the apnea detecting method may be applied to electronic devices such as a piezoelectric film sensor and a cloud server, and may be used for offline detection or online detection. The apnea detection method comprises the following steps:
and S101, extracting the time domain characteristics of the acquired heart attack signals.
In this embodiment, the cardiac shock signal is a cardiac shock time-domain signal acquired by the piezoelectric film sensor, and includes a respiration time-domain signal and a heartbeat time-domain signal. The piezoelectric film is placed under the chest of a user, when the user breathes, different pressures can be applied to the piezoelectric film sensor along with the contraction of the chest, and the breathing condition of the user can be reflected through the cardiac shock signals collected by the piezoelectric film sensor. The acquisition process of the heart impact signal is simple and convenient, professional medical instruments and professional doctors are not needed for operation, and the sleep of a user is not affected.
Optionally, when the apnea detecting method provided by this embodiment is applied to a piezoelectric film sensor, the piezoelectric film sensor may collect a cardiac shock signal, extract a time-domain feature of the cardiac shock signal, and detect an apnea event; when the apnea detection method provided by the embodiment is applied to a cloud server, the piezoelectric film sensor can be used for collecting a cardiac shock signal, the collected cardiac shock signal is sent to the cloud server, and the cloud server extracts time domain features based on the received cardiac shock signal and detects an apnea event.
And S102, extracting the frequency domain feature of the ballistocardiogram signal under the condition that the time domain feature of the ballistocardiogram signal meets the first apnea condition.
In this embodiment, the respiratory waveform at which the apnea event occurs may be collected in advance, and the first apnea condition may be set based on the respiratory waveform at which the apnea event occurs. After the time domain characteristics of the ballistocardiogram signal are extracted, the electronic equipment can judge whether the time domain characteristics of the ballistocardiogram signal meet a preset first apnea condition, namely, whether apnea occurs is judged in a time domain. And if the time domain characteristic of the cardiac shock signal is determined to accord with the first apnea condition, extracting the frequency domain characteristic of the cardiac shock signal so as to judge whether apnea occurs on a frequency domain according to the frequency domain characteristic of the cardiac shock signal.
And step S103, determining that an apnea event exists under the condition that the frequency domain features accord with a second apnea condition.
In this embodiment, the second apnea condition may be set in advance according to the behavior of the respiratory waveform in the frequency domain in which the apnea event occurs. After the frequency domain feature of the cardiac shock signal is extracted, the electronic device can judge whether the frequency domain feature meets a preset second apnea condition, if the frequency domain feature is determined to meet the preset second apnea condition, the cardiac shock signal is judged to be apnea in both a time domain and a frequency domain, and then an apnea event is determined to exist, namely the user is judged to have apnea.
Therefore, in the apnea detection method provided by the embodiment of the present invention, by extracting the time domain feature of the acquired cardiac shock signal, the frequency domain feature of the cardiac shock signal is extracted when the time domain feature of the cardiac shock signal meets the first apnea condition, and the apnea event is determined to exist when the frequency domain feature meets the second apnea condition. Therefore, the time domain characteristic and the frequency domain characteristic of the cardiac shock signal are combined, the apnea event is determined to exist under the condition that the apnea is judged in both the time domain and the frequency domain, the possibility of misjudgment is effectively reduced, and the detection accuracy of the apnea event is improved; the whole detection process does not need to be artificially distinguished, the acquisition of the cardiac shock signal is more convenient, and the operation simplicity of the detection of the apnea event is improved.
Alternatively, as shown in fig. 2, a schematic diagram of a comparison of a respiratory waveform of a normal respiratory rhythm and a respiratory waveform of an apnea. The respiration waveform of the normal respiration rhythm may generate regular peaks and troughs following the contraction of the thoracic cavity, and the respiration waveform of the apnea may not generate regular peaks and troughs, so the time domain feature of the ballistocardiogram signal in this embodiment may include the amplitude of each peak in the ballistocardiogram signal, and before the step S102, it may be determined whether the first apnea condition is met based on the amplitude of each peak in the ballistocardiogram signal. Referring to fig. 3, before step S102, the apnea detecting method may further include:
step S301, determining a plurality of continuous target wave crests according to the time domain characteristics of the heart attack signals; in the plurality of continuous target wave crests, the ratio of the amplitude of the former wave crest of the first target wave crest to the amplitude of the first target wave crest is larger than a preset ratio, the ratio of the amplitude of the latter wave crest of the last target wave crest to the amplitude of the last target wave crest is larger than the preset ratio, and in the amplitudes of any two adjacent target wave crests, the ratio of the larger amplitude to the smaller amplitude is smaller than or equal to the preset ratio.
In this embodiment, because the ratio between the amplitude of the previous peak of the first target peak and the amplitude of the first target peak is greater than the preset ratio, the ratio between the amplitude of the next peak of the last target peak and the amplitude of the last target peak is greater than the preset ratio, and meanwhile, the ratio between the larger amplitude and the smaller amplitude of any two adjacent target peaks is smaller than or equal to the preset ratio, it indicates that the amplitude of the previous peak of the first target peak is far greater than the amplitude of the first target peak, the amplitude of the next peak of the last target peak is far greater than the amplitude of the last target peak, and the amplitudes of the consecutive target peaks are not different from each other. In this case, if the ballistocardiogram signal shows a plurality of continuous peaks with smaller amplitude in the time domain, the duration of the plurality of peaks with smaller amplitude (target peaks) can be calculated and compared with the preset time, and it is determined whether the time domain feature of the ballistocardiogram signal meets the first apnea condition according to the comparison result. It should be understood that, in the embodiment, the amplitudes of the plurality of consecutive target peaks are not different greatly, and the amplitudes of the plurality of consecutive target peaks are not limited to be substantially identical, but the difference between the amplitudes of the plurality of consecutive target peaks is smaller than the difference between the amplitude of the first target peak and the amplitude of the first target peak, and the difference between the amplitude of the last target peak and the amplitude of the last target peak.
It should be noted that the specific value of the preset ratio can be set according to an actual situation, and this is not limited in this embodiment. For example, the predetermined ratio may be 80, 75, 90, etc.
Step S302, if the duration of the multiple continuous target peaks is greater than or equal to the preset time, determining that the time-domain feature of the ballistocardiogram signal meets the first apnea condition.
In this embodiment, the duration may be a difference between a time at which a first target peak and a time at which a last target peak of the plurality of consecutive target peaks is located. When the duration of the plurality of consecutive target peaks is greater than or equal to the preset time, the ballistocardiogram signal may correspond to the respiration waveform of the apnea in fig. 2, and represents the apnea in the time domain, so that it is determined that the temporal characteristic of the ballistocardiogram signal meets the first apnea condition.
Wherein, the preset time can be set according to the actual situation. For example, if the medical definition of apnea refers to complete stop of oronasal airflow for more than 10 seconds (including 10 seconds) during sleep, the preset time may be set to 10 seconds, and when there are multiple continuous target peaks in the ballistocardiogram signal and the duration of the multiple continuous target peaks is greater than or equal to 10 seconds, the temporal feature of the ballistocardiogram signal is determined to meet the first apnea condition.
According to the apnea detection method provided by the embodiment of the invention, after the amplitude of each peak of the cardiac shock signal is extracted, a plurality of continuous target peaks can be determined according to the amplitude of each peak in the cardiac shock signal, and under the condition that the duration time of the plurality of continuous target peaks is greater than or equal to the preset time, the time domain characteristic of the cardiac shock signal is judged to meet a first apnea condition. Therefore, whether the apnea occurs in the cardiac shock signal in the time domain or not is judged according to the duration of the plurality of peaks with smaller amplitude.
In the following, an embodiment for determining a plurality of consecutive target peaks is given. Referring to fig. 4, the step S301 may include the following sub-steps:
and a substep S3011 of performing peak searching on the cardiac shock signal, and determining a first peak and a second peak which are adjacent to each other, wherein the time corresponding to the first peak is less than the time corresponding to the second peak.
In this embodiment, the amplitude corresponding to each time in the ballistocardiograph signal can be represented as x (1), x (2), \ 8230;, x (n), if x (m-1) < x (m) and x (m + 1) < x (m), then x (m) is the amplitude of a peak, and m is the time at which the peak is located. By searching the amplitude of each peak in the ballistocardiogram signal, two continuous peaks, namely a first peak and a second peak which are adjacent to each other, can be found. For example, m and m + p are the time instants at which the adjacent first and second peaks correspond, respectively, and m < m + p, and x (m) and x (m + p) are the amplitudes of the first and second peaks, respectively.
In this embodiment, in the process of searching peaks for the heart shock signal, the electronic device may sequentially search two consecutive peaks according to the time sequence, in this embodiment, a former peak of the two adjacent peaks is referred to as a first peak, and a latter peak of the two adjacent peaks is referred to as a second peak. After finding the adjacent first peak and second peak, comparing the amplitudes of the first peak and the second peak, and if the ratio between the amplitude of the first peak and the amplitude of the second peak is greater than a preset ratio, performing substep S3012; if the ratio between the amplitude of the first peak and the amplitude of the second peak is not larger than the preset ratio, repeating the peak searching action, namely taking the second peak as a new first peak, and searching the next peak adjacent to the new first peak until the ratio between the amplitude of the former peak and the amplitude of the latter peak in the two adjacent peaks found is larger than the preset ratio. For example, assuming that the amplitudes of the first peak and the second peak are x (m) and x (m + p), respectively, and the preset ratio is 80, when x (m) and x (m + p) satisfy x (m) ≦ 80x (m + p), the peak searching action is repeated; when x (m) > 80x (m + p), sub-step S3012 is performed.
In the substep S3012, if the ratio between the amplitude of the first peak and the amplitude of the second peak is greater than the preset ratio, searching for at least one continuous third peak after the time corresponding to the second peak until the ratio between the amplitude of the next peak of the third peak and the amplitude of the third peak is greater than the preset ratio; and in the amplitude values of each third peak and the adjacent previous peak, the ratio of the larger amplitude value to the smaller amplitude value is smaller than or equal to a preset ratio.
In this embodiment, if the ratio between the amplitude of the first peak and the amplitude of the second peak in the two adjacent found peaks is greater than the preset ratio, the time at which the second peak is located is determined, at least one continuous third peak is continuously found after the time corresponding to the second peak, the amplitude of the found third peak is compared with the amplitude of the previous peak adjacent to the third peak, if the ratio between the larger amplitude and the smaller amplitude in the amplitudes of the third peak and the previous peak adjacent to the third peak is less than or equal to the preset ratio, the continuous third peak is continuously found after the time corresponding to the third peak, and when the ratio between the amplitude of the next peak of the found third peak and the amplitude of the third peak is greater than the preset ratio, the peak finding is stopped.
For example, after the time m + p corresponding to the second peak, one or more consecutive third peaks are found, and the ratio of the larger amplitude to the smaller amplitude in the amplitudes of each third peak and the adjacent previous peak is less than or equal to 80. Assuming that the time at which the currently found third peak is located is m + q, the amplitude is x (m + q), the time at which the next peak adjacent to the third peak is located is m + h, and the amplitude is x (m + h), if x (m + q) and x (m + h) satisfy x (m + h) > 80x (m + q), the peak searching is stopped, and substep S3013 is executed.
And a substep S3013, determining the second peak and all the third peaks as target peaks, and obtaining a plurality of continuous target peaks.
In this embodiment, the second peak is a first target peak of the plurality of consecutive target peaks, and the last third peak is a last target peak of the plurality of consecutive target peaks. Because the ratio between the amplitude of the first peak and the amplitude of the second peak is greater than the preset ratio, and for all third peaks continuous to the second peak, in the amplitudes of each third peak and the adjacent previous peak, the ratio between the larger amplitude and the smaller amplitude is smaller than or equal to the preset ratio, and the ratio between the amplitude of the next peak of the last third peak and the amplitude of the last third peak is greater than the preset ratio, the second peak and all third peaks are determined as target peaks, and a plurality of continuous peaks with smaller amplitudes can be determined between the two peaks with larger amplitudes, namely the next peak of the first peak and the last third peak, so as to obtain a plurality of continuous target peaks.
Optionally, referring to fig. 5, the step S102 may include the following sub-steps:
and a substep S1021, determining a time period during which the temporal characteristic of the ballistocardiogram signal meets the first apnea condition.
That is, in the case where it is determined that the temporal characteristic of the cardiac shock signal meets the first apnea condition, it is necessary to determine which temporal characteristic of the cardiac shock signal in which time period meets the first apnea condition.
In one embodiment, the time period may be determined according to the time instants at which the plurality of consecutive target peaks are located. For example, when a first target peak of the plurality of consecutive target peaks is at a time m + p and a last target peak of the plurality of consecutive target peaks is at a time m + q, a time period during which the temporal feature of the ballistocardiogram signal meets the first apnea condition may be determined as (m + p) - (m + q) and a duration of the plurality of consecutive target peaks is (m + q) - (m + p).
Substep S1022, calculating a frequency spectrum of the ballistocardiogram signal in the time period; the frequency domain features include the magnitude of each peak in the frequency spectrum.
In this embodiment, the electronic device may obtain the frequency domain characteristic of the ballistocardiogram signal by calculating the frequency spectrum of the ballistocardiogram signal in the time period and extracting the amplitude of each peak in the frequency spectrum.
In one embodiment, the fourier transform may be performed on the ballistocardiograph signal over the time period to obtain a frequency spectrum of the ballistocardiograph signal over the time period. For example, when the time period during which the time domain feature of the ballistocardiograph signal meets the first apnea condition is determined to be (m + p) to (m + q), the ballistocardiograph signal in the time period can be subjected to continuous fourier transform or discrete fourier transform, and the frequency spectrum of the ballistocardiograph signal in the time period can be obtained.
In this embodiment, the frequency spectrum of the ballistocardiogram signal in the time period may include a frequency spectrum of a heartbeat frequency band and a frequency spectrum of a respiratory frequency band, where the heartbeat frequency band may be 40Hz to 120Hz, and the respiratory frequency band may be 6Hz to 30Hz. Considering that under a normal respiratory rhythm, in the frequency spectrum of the heart attack signal, not only the frequency spectrum of the heartbeat frequency band has a peak with a larger amplitude value, but also the frequency spectrum of the respiratory frequency band has a peak with a larger amplitude value; however, in the frequency spectrum of the apnea heart attack signal, only the frequency spectrum of the heartbeat frequency band has a peak with a large amplitude, and the frequency spectrum of the respiration frequency band does not have a peak with a large amplitude. Therefore, whether the second apnea condition is met can be determined in the present embodiment based on the frequency spectrum of the ballistocardiographic signal during the time period. Referring to fig. 6, before step S103, the apnea detecting method may further include:
step S601, calculating a first average value of amplitudes of all peaks corresponding to the heartbeat frequency band and a second average value of amplitudes of all peaks corresponding to the respiratory frequency band.
In this embodiment, after calculating the frequency spectrum of the cardiac shock signal in the time period, the electronic device may obtain an average value (i.e., a first average value) of the amplitudes of all peaks corresponding to the heartbeat frequency band by extracting the amplitudes of all peaks in the frequency spectrum of the heartbeat frequency band and performing average value calculation; by extracting the amplitudes of all peaks in the frequency spectrum of the respiratory frequency band and performing average value calculation, the average value (i.e., the second average value) of the amplitudes of all peaks corresponding to the respiratory frequency band can be obtained.
Step S602, if at least one peak exists in the heartbeat frequency band, and the amplitude of the peak is greater than a first preset multiple of the first average value, and a second preset multiple of the peak is not present in the respiratory frequency band, it is determined that the frequency domain characteristic meets a second apnea condition.
In this embodiment, the specific values of the first preset multiple and the second preset multiple may be set according to actual conditions, and the first preset multiple and the second preset multiple may be the same or different.
For example, the first preset multiple and the second preset multiple may both be 2, and when a peak with an amplitude greater than 2 times of the first average value exists in the heartbeat frequency band and a peak with an amplitude greater than 2 times of the second average value does not exist in the respiration frequency band, it indicates that in the frequency spectrum of the cardiac shock signal in the time period, the heartbeat frequency band has a peak with a larger amplitude, and the respiration frequency band does not have a peak with a larger amplitude, the cardiac shock signal shows apnea in the frequency domain, so that it is determined that the frequency domain characteristic of the cardiac shock signal meets the second apnea condition.
According to the apnea detection method provided by the embodiment of the invention, whether a peak with a larger amplitude exists in the heartbeat frequency band is judged based on the first average value by calculating the first average value of the amplitudes of all peaks corresponding to the heartbeat frequency band and the second average value of the amplitudes of all peaks corresponding to the respiration frequency band, whether a peak with a larger amplitude exists in the respiration frequency band is judged based on the second average value, when the peak with the larger amplitude exists in the heartbeat frequency band and the peak with the larger amplitude does not exist in the respiration frequency band, the frequency domain characteristic of the cardiac shock signal is judged to accord with the second apnea condition, the accurate judgment that whether the cardiac shock signal in the time period shows apnea in the frequency domain is realized, and the apnea judgment interference under the condition of heartbeat stopping is eliminated.
In order to perform the corresponding steps in the above embodiments and in each possible manner, an implementation of the apnea detecting device is given below. Referring to fig. 7, a functional block diagram of an apnea detecting apparatus 700 according to an embodiment of the present invention is shown. It should be noted that the basic principle and the resulting technical effects of the apnea detecting device 700 provided by the present embodiment are the same as those of the above embodiments, and for the sake of brief description, reference may be made to the corresponding contents in the above embodiments for the parts that are not mentioned in the present embodiment. The apnea detection apparatus 700 includes a time domain feature extraction module 710, a frequency domain feature extraction module 720, and a breath detection module 730.
The time domain feature extraction module 710 is configured to extract a time domain feature of the acquired ballistocardiogram signal.
It is understood that the time-domain feature extraction module 710 may perform the above step S101.
The frequency domain feature extraction module 720 is configured to extract the frequency domain feature of the ballistocardiogram signal under the condition that the time domain feature of the ballistocardiogram signal meets the first apnea condition.
It is understood that the frequency domain feature extraction module 720 can perform the step S102.
The breath detection module 730 is configured to determine that an apnea event exists if the frequency-domain feature satisfies a second apnea condition.
It is understood that the breath detection module 730 may perform the step S103.
Optionally, the frequency-domain feature extraction module 720 may be configured to determine a time period during which the time-domain feature of the ballistocardiogram signal meets the first apnea condition, and calculate a frequency spectrum of the ballistocardiogram signal within the time period; the frequency domain features include the magnitude of each peak in the frequency spectrum.
The frequency domain feature extraction module 720 is specifically configured to perform fourier transform on the ballistocardiograph signal in the time period to obtain a frequency spectrum of the ballistocardiograph signal in the time period.
It is understood that the frequency domain feature extraction module 720 may perform the steps S1021 to S1022.
Optionally, referring to fig. 8, the time-domain feature includes the amplitude of each peak in the ballistocardiograph signal, and the apnea detecting apparatus 700 may further include a first determining module 740 and a second determining module 750.
The first determining module 740 is configured to determine a plurality of consecutive target peaks according to a time domain feature of the ballistocardiogram signal; the ratio of the amplitude of the front peak of the first target peak to the amplitude of the first target peak in the plurality of continuous target peaks is greater than a preset ratio, the ratio of the amplitude of the rear peak of the last target peak to the amplitude of the last target peak is greater than the preset ratio, and the ratio of the larger amplitude to the smaller amplitude in any two adjacent target peaks is less than or equal to the preset ratio; and if the duration time of the continuous target wave crests is greater than or equal to the preset time, judging that the time domain characteristic of the cardiac shock signal meets the first apnea condition.
Wherein the duration is a difference between a time at which the first target peak is located and a time at which the last target peak is located.
In this embodiment, the first determining module 740 is specifically configured to perform peak searching on the cardiac shock signal, and determine a first peak and a second peak that are adjacent to each other, where a time corresponding to the first peak is less than a time corresponding to the second peak; if the ratio of the amplitude of the first peak to the amplitude of the second peak is greater than the preset ratio, searching for at least one continuous third peak after the moment corresponding to the second peak until the ratio of the amplitude of the next peak of the third peak to the amplitude of the third peak is greater than the preset ratio; in the amplitude values of each third peak and the adjacent previous peak, the ratio of the larger amplitude value to the smaller amplitude value is smaller than or equal to a preset ratio; and determining the second peak and all the third peaks as target peaks to obtain a plurality of continuous target peaks.
It is understood that the first determining module 740 may perform the above steps S301 to S302 and the sub steps S3011 to S3013.
In this embodiment, the frequency spectrum includes a frequency spectrum of a heartbeat frequency band and a frequency spectrum of a respiratory frequency band, and the second determining module 750 is configured to calculate a first average value of amplitudes of all peaks corresponding to the heartbeat frequency band and a second average value of amplitudes of all peaks corresponding to the respiratory frequency band; and if the amplitude of at least one peak in the heartbeat frequency band is greater than a first preset multiple of the first average value and the amplitude of the peak in the respiratory frequency band is not greater than a second preset multiple of the second average value, judging that the frequency domain characteristic accords with a second apnea condition.
It is understood that the second determining module 750 can perform the above steps S601 to S602.
In the apnea detection apparatus 700 provided in the embodiment of the present invention, the time domain feature extraction module 710 extracts the time domain feature of the acquired cardiac shock signal, the frequency domain feature extraction module 720 extracts the frequency domain feature of the cardiac shock signal when the time domain feature of the cardiac shock signal meets a first apnea condition, and the apnea detection module 730 determines that an apnea event exists when the frequency domain feature meets a second apnea condition. Therefore, the time domain characteristic and the frequency domain characteristic of the cardiac shock signal are combined, the apnea event is determined to exist under the condition that the apnea is judged in both the time domain and the frequency domain, the possibility of misjudgment is effectively reduced, and the detection accuracy of the apnea event is improved; the whole detection process does not need to be artificially distinguished, the acquisition of the cardiac shock signal is more convenient, and the operation simplicity of the detection of the apnea event is improved.
An electronic device provided by embodiments of the present invention may include a processor and a memory, where at least one instruction, at least one program, a set of codes, or a set of instructions is stored in the memory, and the at least one instruction, the at least one program, the set of codes, or the set of instructions is loaded and executed by the processor to implement the apnea detection method provided by the above method embodiments.
The memory may be used to store software programs and modules, and the processor may execute various functional applications and data processing by operating the software programs and modules stored in the memory. The memory can mainly comprise a program storage area and a data storage area, wherein the program storage area can store an operating system, application programs needed by functions and the like; the storage data area may store data created according to use of the apparatus, and the like. Further, the memory may include high speed random access memory, and may also include non-volatile memory, such as at least one magnetic disk storage device, flash memory device, or other volatile solid state storage device. Accordingly, the memory may also include a memory controller to provide the processor access to the memory.
Fig. 9 is a schematic diagram of a hardware structure of the electronic device 100 according to an embodiment of the present invention. As shown in fig. 9, the electronic device 100 may have a relatively large difference due to different configurations or performances, and may include one or more processors (CPUs) 111 (the processors 111 may include but are not limited to processing devices such as a microprocessor MCU or a programmable logic device FPGA), a memory 113 for storing data, and one or more storage media 112 (e.g., one or more mass storage devices) for storing application programs 1123 or data 1122. The memory 113 and the storage medium 112 may be, among other things, transient storage or persistent storage. The program stored on the storage medium 112 may include one or more modules, each of which may include a sequence of instructions operating on the electronic device 100. Still further, the processor 111 may be configured to communicate with the storage medium 112 to execute a series of instruction operations in the storage medium 112 on the electronic device 100. Electronic apparatus 100 may also include one or more power supplies 116, one or more wired or wireless network interfaces 115, one or more input-output interfaces 114, and/or one or more operating systems 1121, such as windowssserver, macOSXTM, unix, linux, freeBSDTM, and the like.
The input output interface 114 may be used to receive or transmit data via a network. Specific examples of the network described above may include a wireless network provided by a communication provider of the electronic device 100. In one example, the input/output interface 114 includes a network adapter (NIC) that can be connected to other network devices through a base station to communicate with the internet. In one example, the input/output interface 114 can be a Radio Frequency (RF) module, which is used to communicate with the internet in a wireless manner.
It will be understood by those skilled in the art that the structure shown in fig. 9 is only an illustration and is not intended to limit the structure of the electronic device 100. For example, electronic device 100 may also include more or fewer components than shown in FIG. 9, or have a different configuration than shown in FIG. 9.
According to the electronic device provided by the embodiment of the invention, the time domain characteristic of the acquired cardiac shock signal is extracted, the frequency domain characteristic of the cardiac shock signal is extracted under the condition that the time domain characteristic of the cardiac shock signal meets a first apnea condition, and an apnea event is determined to exist under the condition that the frequency domain characteristic meets a second apnea condition. Therefore, the time domain characteristic and the frequency domain characteristic of the cardiac shock signal are combined, the apnea event is determined to exist under the condition that the apnea is judged in both the time domain and the frequency domain, the possibility of misjudgment is effectively reduced, and the detection accuracy of the apnea event is improved; the whole detection process does not need to be artificially distinguished, the acquisition of the cardiac shock signal is more convenient, and the operation simplicity of the detection of the apnea event is improved.
An embodiment of the present invention further provides a computer-readable storage medium, where a computer program is stored on the computer-readable storage medium, and when the computer program is executed by a processor, the computer program implements each process of the above-mentioned apnea detecting method embodiment, and can achieve the same technical effect, and in order to avoid repetition, the computer program is not described herein again. The computer-readable storage medium may be a Read-only memory (ROM), a Random Access Memory (RAM), a magnetic disk or an optical disk.
In the embodiments provided in the present application, it should be understood that the disclosed apparatus and method can be implemented in other ways. The apparatus embodiments described above are merely illustrative and, for example, the flowchart and block diagrams in the figures illustrate the architecture, functionality, and operation of possible implementations of apparatus, methods and computer program products according to various embodiments of the present invention. In this regard, each block in the flowchart or block diagrams may represent a module, segment, or portion of code, which comprises one or more executable instructions for implementing the specified logical function(s). It should also be noted that, in some alternative implementations, the functions noted in the block may occur out of the order noted in the figures. For example, two blocks shown in succession may, in fact, be executed substantially concurrently, or the blocks may sometimes be executed in the reverse order, depending upon the functionality involved. It will also be noted that each block of the block diagrams and/or flowchart illustration, and combinations of blocks in the block diagrams and/or flowchart illustration, can be implemented by special purpose hardware-based systems which perform the specified functions or acts, or combinations of special purpose hardware and computer instructions.
In addition, the functional modules in the embodiments of the present invention may be integrated together to form an independent part, or each module may exist alone, or two or more modules may be integrated to form an independent part.
The functions, if implemented in the form of software functional modules and sold or used as a stand-alone product, may be stored in a computer readable storage medium. Based on such understanding, the technical solution of the present invention or a part thereof which substantially contributes to the prior art may be embodied in the form of a software product, which is stored in a storage medium and includes several instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk, or an optical disk, and various media capable of storing program codes.
The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (10)
1. A method of apnea detection, said method comprising:
extracting time domain characteristics of the collected heart impact signals;
under the condition that the time domain features of the cardiac shock signals accord with a first apnea condition, extracting the frequency domain features of the cardiac shock signals;
determining that an apnea event is present if the frequency domain feature meets a second apnea condition.
2. The method of claim 1, wherein the time-domain feature comprises a magnitude of each peak in the ballistocardiogram signal, and wherein the method further comprises, before the step of extracting the frequency-domain feature of the ballistocardiogram signal if the time-domain feature of the ballistocardiogram signal meets the first apnea condition:
determining a plurality of continuous target peaks according to the time domain characteristics of the heart attack signal; wherein, in the plurality of continuous target peaks, the ratio of the amplitude of the former peak of the first target peak to the amplitude of the first target peak is greater than a preset ratio, the ratio of the amplitude of the latter peak of the last target peak to the amplitude of the last target peak is greater than the preset ratio, and the ratio of the larger amplitude to the smaller amplitude in the amplitudes of any two adjacent target peaks is smaller than or equal to the preset ratio;
and if the duration of the continuous target wave crests is greater than or equal to the preset time, judging that the time domain characteristics of the cardiac shock signals meet a first apnea condition.
3. The method of claim 2, wherein determining a plurality of consecutive target peaks according to the time-domain feature of the ballistocardiogram signal comprises:
searching peaks of the cardiac shock signal, and determining a first peak and a second peak which are adjacent to each other, wherein the time corresponding to the first peak is less than the time corresponding to the second peak;
if the ratio of the amplitude of the first peak to the amplitude of the second peak is greater than the preset ratio, searching for at least one continuous third peak after the moment corresponding to the second peak until the ratio of the amplitude of the next peak to the amplitude of the third peak is greater than the preset ratio; wherein, in the amplitude of each third peak and the adjacent previous peak, the ratio of the larger amplitude to the smaller amplitude is less than or equal to the preset ratio;
and determining the second wave crest and all the third wave crests as target wave crests to obtain a plurality of continuous target wave crests.
4. The method of claim 2, wherein the duration is a difference between a time at which the first target peak is located and a time at which the last target peak is located.
5. The method according to any one of claims 1-4, wherein the extracting the frequency domain features of the ballistocardiogram signal comprises:
determining a time period during which a time domain feature of the ballistocardiographic signal meets a first apnea condition;
calculating a frequency spectrum of the ballistocardiogram signal over the time period; the frequency domain features include the magnitude of each peak in the frequency spectrum.
6. The method of claim 5, wherein the spectrum comprises a spectrum of heartbeat frequency bands and a spectrum of respiratory frequency bands, and wherein the method further comprises, prior to the step of determining that an apneic event is present if the frequency domain feature meets a second apneic condition:
calculating a first average value of the amplitudes of all peaks corresponding to the heartbeat frequency band and a second average value of the amplitudes of all peaks corresponding to the respiratory frequency band;
and if the amplitude of at least one peak in the heartbeat frequency band is greater than a first preset multiple of the first average value and the amplitude of the peak in the respiratory frequency band is not greater than a second preset multiple of the second average value, judging that the frequency domain characteristic meets a second apnea condition.
7. The method of claim 5, wherein the calculating the frequency spectrum of the ballistocardiographic signal over the time period comprises:
and carrying out Fourier transform on the heart attack signal in the time period to obtain the frequency spectrum of the heart attack signal in the time period.
8. An apnea detection apparatus, characterized in that said apparatus comprises:
the time domain characteristic extraction module is used for extracting the time domain characteristics of the acquired heart attack signals;
the frequency domain feature extraction module is used for extracting the frequency domain feature of the cardiac shock signal under the condition that the time domain feature of the cardiac shock signal meets a first apnea condition;
a breath detection module to determine that an apnea event is present if the frequency domain feature meets a second apnea condition.
9. An electronic device, comprising a processor, a memory and a computer program stored on the memory and executable on the processor, the computer program, when executed by the processor, implementing the steps of the apnea detection method of any of claims 1 to 7.
10. A computer-readable storage medium, on which a computer program is stored which, when being executed by a processor, carries out the steps of the apnea detection method of any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110761417.XA CN115568844A (en) | 2021-07-06 | 2021-07-06 | Apnea detection method and device, electronic equipment and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110761417.XA CN115568844A (en) | 2021-07-06 | 2021-07-06 | Apnea detection method and device, electronic equipment and storage medium |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115568844A true CN115568844A (en) | 2023-01-06 |
Family
ID=84579299
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110761417.XA Pending CN115568844A (en) | 2021-07-06 | 2021-07-06 | Apnea detection method and device, electronic equipment and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115568844A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116172527A (en) * | 2023-03-15 | 2023-05-30 | 康原(江苏)科技有限公司 | Intelligent health monitoring method, system and device |
-
2021
- 2021-07-06 CN CN202110761417.XA patent/CN115568844A/en active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116172527A (en) * | 2023-03-15 | 2023-05-30 | 康原(江苏)科技有限公司 | Intelligent health monitoring method, system and device |
| CN116172527B (en) * | 2023-03-15 | 2023-09-22 | 康原(江苏)科技有限公司 | Intelligent health monitoring method, system and device |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10470719B2 (en) | Machine learnt model to detect REM sleep periods using a spectral analysis of heart rate and motion | |
| EP3410924B1 (en) | Machine learnt model to detect rem sleep periods using a spectral analysis of heart rate and motion | |
| KR102267105B1 (en) | Apparatus and method for inspecting sleep disorder based on deep-learning | |
| US11389118B2 (en) | System and method for breathing pattern extraction from PPG signals | |
| CN110742585A (en) | Sleep staging method based on BCG (BCG-broadcast) signals | |
| CN102551699B (en) | The system that non-invasive cardiac output is determined | |
| KR20200024855A (en) | Systems and Methods for Filtering Noise and Analyzing Vein Waveform Signals | |
| CN108324271B (en) | Electrocardiosignal identification method and system and electrocardiosignal monitoring equipment | |
| CN108992053B (en) | Method for real-time non-binding detection of heart rate and heartbeat interval | |
| WO2017053058A1 (en) | Devices, systems, and methods for determining heart rate of a subject from noisy electrocardiogram data | |
| US20170281021A1 (en) | Heartbeat detection method and heartbeat detection device | |
| CN109872820B (en) | Method, device, equipment and storage medium for measuring blood pressure without cuff | |
| CN115568844A (en) | Apnea detection method and device, electronic equipment and storage medium | |
| CN107811610B (en) | Respiration rate detection method and device, electronic equipment and storage medium | |
| JP6842306B2 (en) | Biological signal processing method, its device, and computer-readable recording medium | |
| Mishra et al. | Performance Evaluation of Various Window Techniques for Noise Cancellation from ECG Signal | |
| JP6315633B2 (en) | Heart rate detection method and heart rate detection device | |
| CN113616217B (en) | Method and device for generating baseline drift curve | |
| US20190183362A1 (en) | Systems and methods for filtering medical device noise artifacts from venous waveform signals | |
| Yan et al. | A resource-efficient, robust QRS detector using data compression and time-sharing architecture | |
| JP6922790B2 (en) | Fatigue estimation device and program | |
| WO2018042261A1 (en) | Method and system for monitoring of mental effort | |
| JP6538620B2 (en) | Breathing estimation method and apparatus | |
| CN115517632B (en) | Resonance respiratory frequency measuring method, interactive prompt generating method, device and equipment | |
| CN109567749B (en) | Sleep data processing method and device, computer equipment and storage 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 |