CN115568844A - Apnea detection method and device, electronic equipment and storage medium - Google Patents
Apnea detection method and device, electronic equipment and storage mediumDownload PDFInfo
- Publication number
- CN115568844A CN115568844ACN202110761417.XACN202110761417ACN115568844ACN 115568844 ACN115568844 ACN 115568844ACN 202110761417 ACN202110761417 ACN 202110761417ACN 115568844 ACN115568844 ACN 115568844A
- Authority
- CN
- China
- Prior art keywords
- peak
- apnea
- amplitude
- time
- frequency
- 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—Measuring 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
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Surgery (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Veterinary Medicine (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Physiology (AREA)
- Molecular Biology (AREA)
- Physics & Mathematics (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Pulmonology (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Signal Processing (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
Description
Translated fromChinese技术领域technical field
本发明涉及呼吸检测领域，具体而言，涉及一种呼吸暂停检测方法、装置、电子设备和存储介质。The present invention relates to the field of breathing detection, in particular to a method, device, electronic equipment and storage medium for detecting apnea.
背景技术Background technique
呼吸暂停是指自主呼吸停止，常为暂时性或自限性。由于呼吸暂停通常发生在睡眠过程中，不易察觉，如果长期的呼吸暂停现象未被发现，得不到有效的治疗，就会出现一系列的疾病，因此，检测呼吸暂停对人体健康至关重要，呼吸暂停事件也是睡眠呼吸监测的重要指标。Apnea is the cessation of spontaneous breathing, usually temporary or self-limited. Since apnea usually occurs during sleep, it is not easy to detect. If the long-term apnea is not found and treated effectively, a series of diseases will appear. Therefore, the detection of apnea is very important to human health. Apnea events are also an important indicator of sleep apnea monitoring.
临床上常采用多导睡眠监测仪对睡眠状态进行监测，需要结合专业医师进行判别，即呼吸暂停的监测和判别需要专门的穿戴设备记录和人工判读才能完成，不仅操作复杂，而且会影响用户睡眠。Clinically, polysomnography is often used to monitor the sleep state, which needs to be judged by professional doctors, that is, the monitoring and judgment of apnea requires special wearable equipment records and manual interpretation to complete, which is not only complicated to operate, but also affects the user's sleep .
发明内容Contents of the invention
有鉴于此，本发明的目的在于提供一种呼吸暂停检测方法、装置、电子设备和存储介质，以解决现有技术中检测呼吸暂停事件的准确性不高、操作复杂的问题。In view of this, the purpose of the present invention is to provide an apnea detection method, device, electronic equipment and storage medium to solve the problems of low accuracy and complicated operation in detecting apnea events in the prior art.
为了实现上述目的，本发明实施例采用的技术方案如下：In order to achieve the above object, the technical solution adopted in the embodiment of the present invention is as follows:
第一方面，本发明提供一种呼吸暂停检测方法，所述方法包括：In a first aspect, the present invention provides a method for detecting apnea, the method comprising:
提取采集的心冲击信号的时域特征；Extracting the time-domain features of the collected heart shock signal;
在所述心冲击信号的时域特征符合第一呼吸暂停条件的情况下，提取所述心冲击信号的频域特征；If the time-domain feature of the shock signal meets the first apnea condition, extract the frequency-domain feature of the shock signal;
在所述频域特征符合第二呼吸暂停条件的情况下，确定存在呼吸暂停事件。If the frequency domain feature meets the second apnea condition, it is determined that an apnea event exists.
在可选的实施方式中，所述时域特征包括所述心冲击信号中的每个波峰的幅值，所述在所述心冲击信号的时域特征符合第一呼吸暂停条件的情况下，提取所述心冲击信号的频域特征的步骤之前，所述方法还包括：In an optional implementation manner, the time-domain feature includes the amplitude of each peak in the cardiac shock signal, and when the time-domain feature of the cardiac shock signal meets the first apnea condition, Before the step of extracting the frequency domain features of the shock signal, the method also includes:
根据所述心冲击信号的时域特征确定多个连续的目标波峰；其中，所述多个连续的目标波峰中，第一个目标波峰的前一个波峰的幅值与所述第一个目标波峰的幅值之间的比值大于预设比值，最后一个目标波峰的后一个波峰的幅值与所述最后一个目标波峰的幅值之间的比值大于所述预设比值，并且任意相邻的两个目标波峰的幅值中，较大幅值与较小幅值之间的比值小于或等于所述预设比值；A plurality of continuous target peaks are determined according to the time-domain characteristics of the shock signal; wherein, among the plurality of continuous target peaks, the amplitude of the previous peak of the first target peak is the same as that of the first target peak The ratio between the amplitudes of the last target peak is greater than the preset ratio, the ratio between the amplitude of the last target peak and the amplitude of the last target peak is greater than the preset ratio, and any adjacent two Among the amplitudes of the target peaks, the ratio between the larger amplitude and the smaller amplitude is less than or equal to the preset ratio;
若所述多个连续的目标波峰的持续时间大于或等于预设时间，则判定所述心冲击信号的时域特征符合第一呼吸暂停条件。If the duration of the plurality of consecutive target peaks is greater than or equal to a preset time, it is determined that the time-domain characteristic of the cardiac shock signal meets the first apnea condition.
在可选的实施方式中，所述根据所述心冲击信号的时域特征确定多个连续的目标波峰，包括：In an optional implementation manner, the determining a plurality of continuous target peaks according to the time-domain characteristics of the shock signal includes:
对所述心冲击信号进行寻峰，确定相邻的第一波峰和第二波峰，所述第一波峰对应的时刻小于所述第二波峰对应的时刻；Peak-finding is performed on the cardiac shock signal, and the adjacent first peak and second peak are determined, and the time corresponding to the first peak is smaller 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, then at least one continuous third peak is searched 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, the amplitude of each third peak and the adjacent previous peak Among the values, the ratio between the larger magnitude and the smaller magnitude is less than or equal to the preset ratio;
将所述第二波峰以及所有的第三波峰均确定为目标波峰，得到多个连续的目标波峰。The second peak and all the third peaks are determined as target peaks to obtain a plurality of continuous target peaks.
在可选的实施方式中，所述持续时间为所述第一个目标波峰所处的时刻与所述最后一个目标波峰所处的时刻之间的差值。In an optional implementation manner, the duration is the difference between the moment of the first target peak and the moment of the last target peak.
在可选的实施方式中，所述提取所述心冲击信号的频域特征，包括：In an optional implementation manner, the extracting the frequency-domain features of the shock signal includes:
确定所述心冲击信号的时域特征符合第一呼吸暂停条件的时间段；determining a time period during which the time-domain characteristic of the shock signal meets a first apnea condition;
计算所述心冲击信号在所述时间段内的频谱；所述频域特征包括所述频谱中的每个波峰的幅值。calculating the frequency spectrum of the cardiac shock signal within the time period; the frequency domain features include the amplitude of each peak in the frequency spectrum.
在可选的实施方式中，所述频谱包括心跳频段的频谱和呼吸频段的频谱，所述在所述频域特征符合第二呼吸暂停条件的情况下，确定存在呼吸暂停事件的步骤之前，所述方法还包括：In an optional implementation manner, the frequency spectrum includes a frequency spectrum of a heartbeat frequency band and a frequency spectrum of a respiratory frequency band, and before the step of determining that there is an apnea event when the frequency domain feature meets the second apnea condition, the The method also includes:
计算所述心跳频段对应的所有波峰的幅值的第一平均值，以及所述呼吸频段对应的所有波峰的幅值的第二平均值；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;
若所述心跳频段内存在至少一个波峰的幅值大于所述第一平均值的第一预设倍数，且所述呼吸频段内不存在波峰的幅值大于所述第二平均值的第二预设倍数，则判定所述频域特征符合第二呼吸暂停条件。If there is at least one peak in the heartbeat frequency band whose amplitude is greater than a first preset multiple of the first average value, and there is no peak in the respiratory frequency band whose amplitude is greater than a second preset value of the second average value If the multiple is set, it is determined that the frequency domain feature meets the second apnea condition.
在可选的实施方式中，所述计算所述心冲击信号在所述时间段内的频谱，包括：In an optional implementation manner, the calculating the frequency spectrum of the shock signal within the time period includes:
对所述时间段内的心冲击信号进行傅里叶变换，得到所述心冲击信号在所述时间段内的频谱。Performing Fourier transform on the cardiac shock signal within the time period to obtain a frequency spectrum of the cardiac shock signal within the time period.
第二方面，本发明提供一种呼吸暂停检测装置，所述装置包括：In a second aspect, the present invention provides an apnea detection device, the device comprising:
时域特征提取模块，用于提取采集的心冲击信号的时域特征；A time-domain feature extraction module is used to extract the time-domain features of the collected heart shock signal;
频域特征提取模块，用于在所述心冲击信号的时域特征符合第一呼吸暂停条件的情况下，提取所述心冲击信号的频域特征；A frequency domain feature extraction module, configured to extract the frequency domain feature of the cardiac shock signal when the time domain feature of the cardiac shock signal meets the first apnea condition;
呼吸检测模块，用于在所述频域特征符合第二呼吸暂停条件的情况下，确定存在呼吸暂停事件。A breathing detection module, configured to determine that there is an apnea event when the frequency domain feature meets the second apnea condition.
第三方面，本发明提供一种电子设备，包括处理器、存储器及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述计算机程序被所述处理器执行时实现如前述实施方式任一项所述的呼吸暂停检测方法的步骤。In a third aspect, the present invention provides an electronic device, including a processor, a memory, and a computer program stored on the memory and operable on the processor. When the computer program is executed by the processor, the following The steps of the apnea detection method described in any one of the preceding embodiments.
第四方面，本发明提供一种计算机可读存储介质，所述计算机可读存储介质上存储计算机程序，所述计算机程序被处理器执行时实现如前述实施方式任一项所述的呼吸暂停检测方法的步骤。In a fourth aspect, the present invention provides a computer-readable storage medium, on which a computer program is stored, and when the computer program is executed by a processor, the detection of apnea as described in any one of the preceding embodiments is realized. method steps.
本发明实施例提供的呼吸暂停检测方法、装置、电子设备和存储介质，通过提取采集的心冲击信号的时域特征，在心冲击信号的时域特征符合第一呼吸暂停条件的情况下，提取该心冲击信号的频域特征，在该频域特征符合第二呼吸暂停条件的情况下，确定存在呼吸暂停事件。如此，实现了将心冲击信号的时域特征和频域特征相结合，在时域和频域中均判定呼吸暂停的情况下确定存在呼吸暂停事件，有效减少了误判的可能，提高了呼吸暂停事件的检测准确性；整个检测过程无需人为判别，心冲击信号的获取也比较方便，提高了呼吸暂停事件检测的操作简便性。The apnea detection method, device, electronic device, and storage medium provided by the embodiments of the present invention extract the time-domain features of the collected cardiac shock signals, and extract the The frequency domain feature of the cardiac shock signal. If the frequency domain feature meets the second apnea condition, it is determined that there is an apnea event. In this way, the combination of time-domain features and frequency-domain features of cardiac shock signals is realized, and the existence of apnea events is determined when apnea is judged in both the time-domain and frequency-domain, which effectively reduces the possibility of misjudgment and improves respiratory rate. The detection accuracy of the pause event; the whole detection process does not need to be judged manually, and the acquisition of the cardiac shock signal is also relatively convenient, which improves the operation simplicity of the detection of the pause event.
为使本发明的上述目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附附图，作详细说明如下。In order to make the above-mentioned objects, features and advantages of the present invention more comprehensible, preferred embodiments will be described in detail below together with the accompanying drawings.
附图说明Description of drawings
为了更清楚地说明本发明实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本发明的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。In order to illustrate the technical solutions of the embodiments of the present invention more clearly, the accompanying drawings used in the embodiments will be briefly introduced below. It should be understood that the following drawings only show some embodiments of the present invention, and thus It should be regarded as a limitation on the scope, and those skilled in the art can also obtain other related drawings based on these drawings without creative work.
图1示出了本发明实施例提供的呼吸暂停检测方法的一种流程示意图；FIG. 1 shows a schematic flow chart of a method for detecting apnea provided by an embodiment of the present invention;
图2示出了正常呼吸节律的呼吸波形与呼吸暂停的呼吸波形的一种对比示意图；Fig. 2 shows a comparison schematic diagram of the respiratory waveform of normal respiratory rhythm and the respiratory waveform of apnea;
图3示出了本发明实施例提供的呼吸暂停检测方法的另一种流程示意图；Fig. 3 shows another schematic flowchart of the apnea detection method provided by the embodiment of the present invention;
图4示出了本发明实施例提供的呼吸暂停检测方法的又一种流程示意图；Fig. 4 shows another schematic flowchart of the apnea detection method provided by the embodiment of the present invention;
图5示出了本发明实施例提供的呼吸暂停检测方法的又一种流程示意图；Fig. 5 shows another schematic flowchart of the apnea detection method provided by the embodiment of the present invention;
图6示出了本发明实施例提供的呼吸暂停检测方法的又一种流程示意图；Fig. 6 shows another schematic flowchart of the apnea detection method provided by the embodiment of the present invention;
图7示出了本发明实施例提供的呼吸暂停检测装置的一种功能模块图；Fig. 7 shows a functional block diagram of the apnea detection device provided by the embodiment of the present invention;
图8示出了本发明实施例提供的呼吸暂停检测装置的另一种功能模块图；Fig. 8 shows another functional block diagram of the apnea detection device provided by the embodiment of the present invention;
图9示出了本发明实施例提供的电子设备的一种硬件结构示意图。FIG. 9 shows a schematic diagram of a hardware structure of an electronic device provided by an embodiment of the present invention.
图标：100-电子设备；700-呼吸暂停检测装置；111-处理器；112-存储介质；113-存储器；114-输入输出接口；115-有线或无线网络接口；116-电源；1121-操作系统；1122-数据；1123-应用程序；710-时域特征提取模块；720-频域特征提取模块；730-呼吸检测模块；740-第一判定模块；750-第二判定模块。Icons: 100-electronic equipment; 700-apnea detection device; 111-processor; 112-storage medium; 113-memory; 114-input and output interface; 115-wired or wireless network interface; 116-power supply; 1121-operating system 1122-data; 1123-application program; 710-time domain feature extraction module; 720-frequency domain feature extraction module; 730-breathing detection module; 740-first judgment module; 750-second judgment module.
具体实施方式detailed description
下面将结合本发明实施例中附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。The following will clearly and completely describe the technical solutions in the embodiments of the present invention with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only some, not all, embodiments of the present invention. The components of the embodiments of the invention generally described and illustrated in the figures herein may be arranged and designed in a variety of different configurations.
因此，以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围，而是仅仅表示本发明的选定实施例。基于本发明的实施例，本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本发明保护的范围。Accordingly, the following detailed description of the embodiments of the invention provided in the accompanying drawings is not intended to limit the scope of the claimed invention, but merely represents selected embodiments of the invention. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without making creative efforts belong to the protection scope of the present invention.
需要说明的是，术语“第一”和“第二”等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。It should be noted that relative terms such as the terms "first" and "second" are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any relationship between these entities or operations. There is no such actual relationship or order between them. Furthermore, the term "comprises", "comprises" or any other variation thereof is intended to cover a non-exclusive inclusion such that a process, method, article, or apparatus comprising a set of elements includes not only those elements, but also includes elements not expressly listed. other elements of or also include elements inherent in such a process, method, article, or device. Without further limitations, an element defined by the phrase "comprising a ..." does not exclude the presence of additional identical elements in the process, method, article or apparatus comprising said element.
在一个实施例中，如图1所示，提供了一种呼吸暂停检测方法，该呼吸暂停检测方法可以应用于压电薄膜传感器、云端服务器等电子设备中，既可以是离线检测，也可以是在线检测。该呼吸暂停检测方法包括以下步骤：In one embodiment, as shown in FIG. 1 , a method for detecting apnea is provided. The method for detecting apnea can be applied to electronic devices such as piezoelectric film sensors and cloud servers. It can be either offline detection or online test. The apnea detection method comprises the following steps:
步骤S101，提取采集的心冲击信号的时域特征。Step S101, extracting the time-domain features of the collected cardiac shock signal.
在本实施例中，该心冲击信号是由压电薄膜传感器采集得到的心冲击时域信号，包括呼吸时域信号和心跳时域信号。将压电薄膜放置于用户的胸腔下，当用户呼吸时，随着胸腔的收缩会对压电薄膜传感器施加不同的压力，通过压电薄膜传感器采集到的心冲击信号，可以反映用户的呼吸状况。该心冲击信号的获得过程简单方便，不需要专业的医疗器械和专业医师操作，也不会对用户的睡眠造成影响。In this embodiment, the cardiac shock signal is a cardiac shock time-domain signal collected by the piezoelectric film sensor, including a respiratory time-domain signal and a heartbeat time-domain signal. Place the piezoelectric film under the user's chest cavity. When the user breathes, different pressures will be applied to the piezoelectric film sensor with the contraction of the chest cavity. The heart shock signal collected by the piezoelectric film sensor can reflect the user's breathing condition. . The process of obtaining the heart shock signal is simple and convenient, does not require professional medical equipment and operation by professional doctors, and will not affect the sleep of the user.
可选地，当本实施例提供的呼吸暂停检测方法是应用在压电薄膜传感器上时，可以由该压电薄膜传感器采集心冲击信号，提取心冲击信号的时域特征并进行呼吸暂停事件的检测；当本实施例提供的呼吸暂停检测方法是应用在云端服务器上时，可以由该压电薄膜传感器采集心冲击信号，并将采集的心冲击信号发送给云端服务器，由云端服务器基于接收到的心冲击信号提取时域特征并进行呼吸暂停事件的检测。Optionally, when the apnea detection method provided in this embodiment is applied to a piezoelectric film sensor, the piezoelectric film sensor can collect cardiac shock signals, extract the time-domain characteristics of the cardiac shock signals, and perform detection of apnea events. detection; when the apnea detection method provided by this embodiment is applied on the cloud server, the piezoelectric film sensor can collect the heart shock signal, and send the collected heart shock signal to the cloud server, and the cloud server can receive the heart shock signal based on the The time-domain features are extracted from the cardiac shock signal and the apnea event is detected.
步骤S102，在心冲击信号的时域特征符合第一呼吸暂停条件的情况下，提取心冲击信号的频域特征。Step S102, if the time-domain feature of the cardiac shock signal meets the first apnea condition, extract the frequency-domain feature of the cardiac shock signal.
在本实施例中，可以预先收集出现呼吸暂停事件的呼吸波形，并基于出现呼吸暂停事件的呼吸波形设置第一呼吸暂停条件。电子设备在提取到心冲击信号的时域特征后，可判断该心冲击信号的时域特征是否符合预先设置的第一呼吸暂停条件，即在时域上判断是否出现呼吸暂停。若确定心冲击信号的时域特征符合第一呼吸暂停条件，则提取出心冲击信号的频域特征，以便根据该心冲击信号的频域特征在频域上判断是否出现呼吸暂停。In this embodiment, respiratory waveforms where apnea events occur may be collected in advance, and the first apnea condition is set based on the respiratory waveforms where apnea events occur. After the electronic device extracts the time-domain feature of the cardiac shock signal, it can judge whether the time-domain feature of the cardiac shock signal meets the preset first apnea condition, that is, judge whether apnea occurs in the time domain. If it is determined that the time-domain feature of the shock signal meets the first apnea condition, the frequency-domain feature of the shock signal is extracted, so as to determine whether apnea occurs in the frequency domain according to the frequency-domain feature of the shock signal.
步骤S103，在频域特征符合第二呼吸暂停条件的情况下，确定存在呼吸暂停事件。Step S103, if the frequency domain feature meets the second apnea condition, determine that there is an apnea event.
在本实施例中，可以预先根据出现呼吸暂停事件的呼吸波形在频域上的表现，设置第二呼吸暂停条件。电子设备在提取出心冲击信号的频域特征后，可判断该频域特征是否符合预先设置的第二呼吸暂停条件，若确定该频域特征符合预先设置的第二呼吸暂停条件，则判定该心冲击信号在时域和频域中均表现为呼吸暂停，进而确定存在呼吸暂停事件，即判定该用户出现呼吸暂停。In this embodiment, the second apnea condition may be set in advance according to the performance of the respiratory waveform in the frequency domain where the apnea event occurs. After the electronic device extracts the frequency domain feature of the cardiac shock signal, it can determine whether the frequency domain feature meets the preset second apnea condition, and if it is determined that the frequency domain feature meets the preset second apnea condition, then determine the The cardiac shock signal appears as apnea in both the time domain and the frequency domain, and then it is determined that there is an apnea event, that is, it is determined that the user has apnea.
可见，本发明实施例提供的呼吸暂停检测方法，通过提取采集的心冲击信号的时域特征，在心冲击信号的时域特征符合第一呼吸暂停条件的情况下，提取该心冲击信号的频域特征，在该频域特征符合第二呼吸暂停条件的情况下，确定存在呼吸暂停事件。如此，实现了将心冲击信号的时域特征和频域特征相结合，在时域和频域中均判定呼吸暂停的情况下确定存在呼吸暂停事件，有效减少了误判的可能，提高了呼吸暂停事件的检测准确性；整个检测过程无需人为判别，心冲击信号的获取也比较方便，提高了呼吸暂停事件检测的操作简便性。It can be seen that the apnea detection method provided by the embodiment of the present invention extracts the time-domain features of the collected cardiac shock signal, and extracts the frequency-domain feature of the cardiac shock signal when the time-domain feature of the cardiac shock signal meets the first apnea condition. feature, if the frequency domain feature meets the second apnea condition, it is determined that there is an apnea event. In this way, the combination of time-domain features and frequency-domain features of cardiac shock signals is realized, and the existence of apnea events is determined when apnea is judged in both the time-domain and frequency-domain, which effectively reduces the possibility of misjudgment and improves respiratory rate. The detection accuracy of the pause event; the whole detection process does not need to be judged manually, and the acquisition of the cardiac shock signal is also relatively convenient, which improves the operation simplicity of the detection of the pause event.
可选地，如图2所示，为正常呼吸节律的呼吸波形与呼吸暂停的呼吸波形的一种对比示意图。正常呼吸节律的呼吸波形会跟随胸腔的收缩产生有规律的波峰与波谷，而有呼吸暂停的呼吸波形不会产生有规律的波峰与波谷，因此，本实施例中的心冲击信号的时域特征可以包括心冲击信号中的每个波峰的幅值，在上述步骤S102之前，可以基于心冲击信号中的每个波峰的幅值判定是否符合第一呼吸暂停条件。请参照图3，在步骤S102之前，该呼吸暂停检测方法还可以包括：Optionally, as shown in FIG. 2 , it is a schematic diagram of a comparison between the respiratory waveform of a normal respiratory rhythm and the respiratory waveform of an apnea. The respiratory waveform of normal respiratory rhythm will follow the contraction of the chest cavity to produce regular peaks and troughs, while the respiratory waveform with apnea will not produce regular peaks and troughs. Therefore, the time domain characteristics of the cardiac shock signal in this embodiment The amplitude of each peak in the cardiac shock signal may be included, and before the above step S102, it may be determined based on the amplitude of each peak in the cardiac shock signal whether the first apnea condition is met. Please refer to FIG. 3, before step S102, the apnea detection method may also include:
步骤S301，根据心冲击信号的时域特征确定多个连续的目标波峰；其中，多个连续的目标波峰中，第一个目标波峰的前一个波峰的幅值与第一个目标波峰的幅值之间的比值大于预设比值，最后一个目标波峰的后一个波峰的幅值与最后一个目标波峰的幅值之间的比值大于预设比值，并且任意相邻的两个目标波峰的幅值中，较大幅值与较小幅值之间的比值小于或等于预设比值。Step S301, determine a plurality of continuous target peaks according to the time-domain characteristics of the cardiac shock signal; wherein, among the multiple continuous target peaks, the amplitude of the previous peak of the first target peak is the same as the amplitude of the first target peak The ratio between is greater than the preset ratio, the ratio between the amplitude of the last peak of the last target peak and the amplitude of the last target peak is greater than the preset ratio, and the amplitude of any adjacent two target peaks , the ratio between the larger magnitude and the smaller magnitude is less than or equal to the preset ratio.
在本实施例中，由于在该多个连续的目标波峰中，第一个目标波峰的前一个波峰的幅值与第一个目标波峰的幅值之间的比值大于预设比值，最后一个目标波峰的后一个波峰的幅值与最后一个目标波峰的幅值之间的比值大于预设比值，同时，任意相邻的两个目标波峰的幅值中，较大幅值与较小幅值之间的比值小于或等于预设比值，表明第一个目标波峰的前一个波峰的幅值远大于第一目标波峰的幅值，最后一个目标波峰的后一个波峰的幅值远大于最后一个目标波峰的幅值，并且该多个连续的目标波峰的幅值相差不大。在此情形下，该心冲击信号在时域上表现为有连续的多个幅值较小的波峰，则可以计算该多个幅值较小的波峰(目标波峰)的持续时间，并与预设时间进行比较，根据比较结果判断该心冲击信号的时域特征是否符合第一呼吸暂停条件。应当理解的是，本实施例中所提到的多个连续的目标波峰的幅值相差不大，并不是限定这多个连续的目标波峰的幅值都是基本一致的，而是要说明多个连续的目标波峰的幅值之间的差距，相对于第一个目标波峰的前一个波峰与第一个目标波峰的幅值之间的差距，以及最后一个目标波峰的后一个波峰与最后一个目标波峰的幅值之间的差距来说，是要更小的。In this embodiment, because among the multiple consecutive target peaks, 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 last target The ratio between the amplitude of the last peak of the peak and the amplitude of the last target peak is greater than the preset ratio, and at the same time, among the amplitudes of any adjacent two target peaks, there The ratio is less than or equal to the preset ratio, indicating that the amplitude of the previous peak of the first target peak is much greater than the amplitude of the first target peak, and the amplitude of the subsequent peak of the last target peak is much greater than that of the last target peak. amplitude, and the amplitudes of the multiple consecutive target peaks have little difference. In this case, the cardiac shock signal appears to have multiple continuous peaks with smaller amplitudes in the time domain, then the duration of the multiple smaller amplitude peaks (target peaks) can be calculated and compared with the predicted Time is set for comparison, and it is judged according to the comparison result whether the time-domain feature of the cardiac shock signal meets the first apnea condition. It should be understood that the amplitudes of the multiple consecutive target peaks mentioned in this embodiment are not much different. The difference between the amplitudes of successive target peaks, the difference between the previous peak and the first target peak relative to the first target peak, and the difference between the last target peak and the last target peak In terms of the gap between the amplitudes of the target peaks, it should be smaller.
需要说明的是，该预设比值的具体取值可以根据实际情况设定，本实施例中对此不作限制。例如，该预设比值可以为80、75、90等。It should be noted that, the specific value of the preset ratio can be set according to actual conditions, which is not limited in this embodiment. For example, the preset ratio can be 80, 75, 90 and so on.
步骤S302，若多个连续的目标波峰的持续时间大于或等于预设时间，则判定心冲击信号的时域特征符合第一呼吸暂停条件。Step S302, if the duration of multiple consecutive target peaks is greater than or equal to the preset time, it is determined that the time-domain characteristic of the shock signal meets the first apnea condition.
在本实施例中，该持续时间可以为该多个连续的目标波峰中，第一个目标波峰所处的时刻与最后一个目标波峰所处的时刻之间的差值。当多个连续的目标波峰的持续时间大于或等于预设时间，此时心冲击信号可对应于图2中呼吸暂停的呼吸波形，在时域上表现为呼吸暂停，故判定该心冲击信号的时域特征符合第一呼吸暂停条件。In this embodiment, the duration may be the difference between the moment of the first target peak and the moment of the last target peak among the multiple consecutive target peaks. When the duration of multiple continuous target peaks is greater than or equal to the preset time, the cardiac shock signal can correspond to the respiratory waveform of apnea in Fig. Temporal features correspond to the first apnea condition.
其中，该预设时间可以根据实际情况设定。例如，医学上对于呼吸暂停的定义是指睡眠过程中口鼻气流完全停止10秒以上(含10秒)，则可以设置该预设时间为10秒，当心冲击信号中存在多个连续的目标波峰并且该多个连续的目标波峰的持续时间大于或者等于10秒，则判定该心冲击信号的时域特征符合第一呼吸暂停条件。Wherein, the preset time can be set according to actual conditions. For example, the medical definition of apnea refers to the complete cessation of airflow through the mouth and nose for more than 10 seconds (including 10 seconds) during sleep. You can set the preset time to 10 seconds, and beware of multiple continuous target peaks in the shock signal And if the duration of the multiple consecutive target peaks is greater than or equal to 10 seconds, then it is determined that the time-domain feature of the cardiac shock signal meets the first apnea condition.
本发明实施例提供的呼吸暂停检测方法，在提取出心冲击信号的每个波峰的幅值后，可根据心冲击信号中的每个波峰的幅值，确定多个连续的目标波峰，在该多个连续的目标波峰的持续时间大于或等于预设时间的情况下，判定该心冲击信号的时域特征符合第一呼吸暂停条件。如此，实现了通过多个幅值较小的波峰的持续时间判定心冲击信号在时域上是否出现呼吸暂停。In the apnea detection method provided by the embodiment of the present invention, after extracting the amplitude of each peak of the cardiac shock signal, a plurality of continuous target peaks can be determined according to the amplitude of each peak in the cardiac shock signal. When the duration of multiple consecutive target peaks is greater than or equal to the preset time, it is determined that the time-domain feature of the cardiac shock signal meets the first apnea condition. In this way, the determination of whether apnea occurs in the time domain of the shock signal is realized through the duration of multiple peaks with smaller amplitudes.
下面，给出确定多个连续的目标波峰的一种具体实施方式。请参照图4，上述的步骤S301可以包括如下子步骤：A specific implementation manner of determining multiple continuous target peaks is given below. Please refer to FIG. 4, the above-mentioned step S301 may include the following sub-steps:
子步骤S3011，对心冲击信号进行寻峰，确定相邻的第一波峰和第二波峰，第一波峰对应的时刻小于第二波峰对应的时刻。Sub-step S3011, peak-seeking the cardiac shock signal, and determining adjacent first and second peaks, the time corresponding to the first peak is smaller than the time corresponding to the second peak.
在本实施例中，该心冲击信号中每个时刻对应的幅值可以表示为x(1)，x(2)，…，x(n)，若x(m-1)<x(m)且x(m+1)<x(m)，则表示x(m)为一个波峰的幅值，m是该波峰所处的时刻。通过在心冲击信号中搜索各波峰的幅值，可以找到连续的两个波峰，即相邻的第一波峰和第二波峰。例如，m和m+p分别是相邻的第一波峰和第二波峰对应的时刻，且m＜m+p，x(m)和x(m+p)分别是第一波峰和第二波峰的幅值。In this embodiment, the amplitude corresponding to each moment in the cardiac shock signal can be expressed as x(1), x(2), ..., x(n), if x(m-1)<x(m) And x(m+1)<x(m), it means that x(m) is the amplitude of a peak, and m is the moment of the peak. By searching the amplitude of each peak in the cardiac shock signal, two consecutive peaks can be found, that is, adjacent first and second peaks. For example, m and m+p are the moments corresponding to the adjacent first and second peaks respectively, and m<m+p, x(m) and x(m+p) are the first and second peaks respectively the magnitude of .
在本实施例中，电子设备在对心冲击信号进行寻峰的过程中，可以按照时间的先后顺序依次寻找连续的两个波峰，本实施例中将相邻两个波峰中的前一个波峰称为第一波峰，后一个波峰称为第二波峰。当找到相邻的第一波峰和第二波峰后，对第一波峰和第二波峰的幅值进行比较，如果第一波峰的幅值和第二波峰的幅值之间的比值大于预设比值，则执行子步骤S3012；如果第一波峰的幅值和第二波峰的幅值之间的比值不大于预设比值，则重复寻峰的动作，即将第二波峰作为新的第一波峰，并寻找与新的第一波峰相邻的下一个波峰，直到找到的相邻两个波峰中的前一个波峰的幅值与后一个波峰的幅值之间的比值大于预设比值。例如，假设第一波峰和第二波峰的幅值分别为x(m)和x(m+p)，预设比值为80，当x(m)和x(m+p)满足x(m)≤80x(m+p)，则重复寻峰的动作；当x(m)＞80x(m+p)，则执行子步骤S3012。In this embodiment, in the process of finding the peak of the cardiac shock signal, the electronic device can search for two consecutive peaks sequentially according to the order of time. In this embodiment, the previous peak of the two adjacent peaks is called is the first peak, and the latter peak is called the second peak. After finding the adjacent first peak and the second peak, compare the amplitudes of the first peak and 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 , then perform sub-step S3012; if the ratio between the amplitude of the first peak and the amplitude of the second peak is not greater than the preset ratio, then repeat the action of finding the peak, that is, the second peak is used as the new first peak, and The next peak adjacent to the new first peak is searched until the ratio between the amplitude of the previous peak and the amplitude of the next peak among the found two adjacent peaks is greater 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), repeat the peak-seeking action; when x(m)>80x(m+p), execute sub-step S3012.
子步骤S3012，若第一波峰的幅值与第二波峰的幅值之间的比值大于预设比值，则在第二波峰对应的时刻之后寻找连续的至少一个第三波峰，直到第三波峰的后一个波峰的幅值与第三波峰的幅值之间的比值大于预设比值；其中，每个第三波峰与相邻的前一个波峰的幅值中，较大幅值与较小幅值之间的比值小于或等于预设比值。Sub-step S3012, if the ratio between the amplitude of the first peak and the amplitude of the second peak is greater than the preset ratio, then search for at least one continuous third peak after the time corresponding to the second peak until the third peak The ratio between the amplitude of the latter peak and the amplitude of the third peak is greater than the preset ratio; wherein, among the amplitudes of each third peak and the adjacent previous peak, the difference between the larger amplitude and the smaller amplitude The ratio between is less than or equal to the preset ratio.
在本实施例中，若找到的相邻两个波峰中，第一波峰的幅值与第二波峰的幅值之间的比值大于预设比值，则确定第二波峰所处的时刻，并在第二波峰对应的时刻之后继续寻找连续的至少一个第三波峰，将找到的第三波峰的幅值与该第三波峰相邻的前一个波峰的幅值进行比较，如果该第三波峰与相邻的前一个波峰的幅值中，较大幅值与较小幅值之间的比值小于或等于预设比值，则在该第三波峰对应的时刻之后继续寻找连续的第三波峰，当找到的第三波峰的后一个波峰的幅值与该第三波峰的幅值之间的比值大于预设比值，则停止寻峰。In this embodiment, if the ratio between the amplitude of the first peak and the amplitude of the second peak is greater than the preset ratio among the two adjacent peaks found, the moment of the second peak is determined, and at After the moment corresponding to the second peak, continue to search for at least one continuous third peak, compare the amplitude of the found third peak with the amplitude of the previous peak adjacent to the third peak, if the third peak is the same as If the ratio between the larger amplitude and the smaller amplitude is less than or equal to the preset ratio among the amplitudes of the previous adjacent peak, then continue to search for a continuous third peak after the moment corresponding to the third peak, when the found If the ratio between the amplitude of the third peak after the third peak and the amplitude of the third peak is greater than the preset ratio, the peak finding is stopped.
例如，在第二波峰对应的时刻m+p之后，寻找连续的一个或多个第三波峰，每个第三波峰与相邻的前一个波峰的幅值中，较大幅值与较小幅值之间的比值小于或等于80。假设当前找到的第三波峰所处的时刻为m+q，幅值为x(m+q)，与该第三波峰相邻的后一个波峰所处的时刻为m+h，幅值为x(m+h)，若x(m+q)和x(m+h)满足x(m+h)＞80x(m+q)，则停止寻峰，执行子步骤S3013。For example, after the moment m+p corresponding to the second peak, one or more continuous third peaks are searched, and the amplitude of each third peak and the adjacent previous peak has a larger amplitude and a smaller amplitude The ratio between them is less than or equal to 80. Assume that the moment of the currently found third peak is m+q, and the amplitude is x(m+q), and the moment of the next peak adjacent to the third peak 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), then stop peak-seeking and execute sub-step S3013.
子步骤S3013，将第二波峰以及所有的第三波峰均确定为目标波峰，得到多个连续的目标波峰。Sub-step S3013, determining the second peak and all the third peaks as target peaks to obtain multiple consecutive target peaks.
在本实施例中，该第二波峰即为多个连续的目标波峰中的第一个目标波峰，最后一个第三波峰即为多个连续的目标波峰中的最后一个目标波峰。由于第一波峰的幅值与第二波峰的幅值之间的比值大于预设比值，而对于与第二波峰连续的所有第三波峰，每个第三波峰与相邻的前一个波峰的幅值中，较大幅值与较小幅值之间的比值小于或等于预设比值，且最后一个第三波峰的后一个波峰的幅值与该最后一个第三波峰的幅值之间的比值大于预设比值，故通过将该第二波峰和所有的第三波峰确定为目标波峰，可在第一波峰和最后一个第三波峰的后一个波峰这两个幅值较大的波峰之间，确定出连续的多个幅值较小的波峰，从而得到多个连续的目标波峰。In this embodiment, the second peak is the first target peak among the multiple consecutive target peaks, and the last third peak is the last target peak among the multiple consecutive target peaks. Since 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 the third peaks continuous with the second peak, the amplitude of each third peak to the adjacent previous peak Among the values, the ratio between the larger amplitude and the smaller amplitude is less than or equal to the preset ratio, and the ratio between the amplitude of the last third peak and the amplitude of the last third peak is greater than preset ratio, so by determining the second peak and all the third peaks as target peaks, it can be determined between the first peak and the last peak of the last third peak, which are two larger peaks. Multiple continuous peaks with smaller amplitudes are generated, so as to obtain multiple continuous target peaks.
可选地，请参照图5，上述的步骤S102可以包括如下子步骤：Optionally, please refer to FIG. 5, the above-mentioned step S102 may include the following sub-steps:
子步骤S1021，确定心冲击信号的时域特征符合第一呼吸暂停条件的时间段。Sub-step S1021, determine the time period in which the time-domain characteristics of the shock signal meet the first apnea condition.
也即是，在判定心冲击信号的时域特征符合第一呼吸暂停条件的情况下，需要确定是哪一个时间段内的心冲击信号的时域特征符合第一呼吸暂停条件。That is, when it is determined that the time-domain characteristic of the shock signal meets the first apnea condition, it needs to be determined in which time period the time-domain characteristic of the shock signal meets the first apnea condition.
在一种实施方式中，可以根据上述的多个连续的目标波峰所处的时刻确定该时间段。例如，当该多个连续的目标波峰中的第一个目标波峰所处的时刻为m+p，该多个连续的目标波峰中的最后一个目标波峰所处的时刻为m+q，则可确定心冲击信号的时域特征符合第一呼吸暂停条件的时间段为(m+p)～(m+q)，该多个连续的目标波峰的持续时间为(m+q)-(m+p)。In an implementation manner, the time period may be determined according to the times at which the above-mentioned multiple continuous target peaks are located. For example, when the time of the first target peak in the multiple continuous target peaks is m+p, and the time of the last target peak in the multiple continuous target peaks is m+q, then it can be It is determined that the time period during which the time-domain characteristics of the cardiac shock signal meet the first apnea condition is (m+p)-(m+q), and the duration of the multiple continuous target peaks is (m+q)-(m+ p).
子步骤S1022，计算心冲击信号在时间段内的频谱；频域特征包括频谱中的每个波峰的幅值。Sub-step S1022, calculating the frequency spectrum of the cardiac shock signal within a time period; the frequency domain feature includes the amplitude of each peak in the frequency spectrum.
在本实施例中，电子设备通过计算心冲击信号在该时间段内的频谱，并提取出频谱中每个波峰的幅值，可以得到该心冲击信号的频域特征。In this embodiment, the electronic device can obtain the frequency domain characteristics of the cardiac shock signal by calculating the frequency spectrum of the cardiac shock signal within the time period and extracting the amplitude of each peak in the frequency spectrum.
在一种实施方式中，可以对时间段内的心冲击信号进行傅里叶变换，得到心冲击信号在时间段内的频谱。例如，当确定心冲击信号的时域特征符合第一呼吸暂停条件的时间段为(m+p)～(m+q)，则可以对这个时间段内的心冲击信号进行连续傅里叶变换或者离散傅里叶变换，得到心冲击信号在这个时间段内的频谱。In an implementation manner, Fourier transform may be performed on the cardiac shock signal within the time period to obtain the frequency spectrum of the cardiac shock signal within the time period. For example, when it is determined that the time domain characteristic of the cardiac shock signal conforms to the first apnea condition, the time period is (m+p)~(m+q), then the continuous Fourier transform can be performed on the cardiac shock signal within this time period Or discrete Fourier transform to obtain the frequency spectrum of the cardiac shock signal within this time period.
在本实施例中，该心冲击信号在该时间段内的频谱可以包括心跳频段的频谱和呼吸频段的频谱，其中，心跳频段可以为40Hz～120Hz，呼吸频段可以为6Hz～30Hz。考虑到在正常的呼吸节律下，心冲击信号的频谱中，不仅心跳频段的频谱会存在幅值较大的波峰，而且在呼吸频段的频谱也会存在幅值较大的波峰；但在呼吸暂停的心冲击信号的频谱中，只有心跳频段的频谱中存在幅值较大的波峰，呼吸频段的频谱中并不存在幅值很大的波峰。因此，本实施例中可以基于心冲击信号在该时间段内的频谱判定是否符合第二呼吸暂停条件。请参照图6，在步骤S103之前，该呼吸暂停检测方法还可以包括：In this embodiment, the frequency spectrum of the cardiac shock signal within the time period may include a spectrum of a heartbeat frequency range and a spectrum of a respiratory frequency range, wherein the heartbeat frequency range may be 40 Hz-120 Hz, and the respiratory frequency range may be 6 Hz-30 Hz. Considering that in the normal respiratory rhythm, in the frequency spectrum of the cardiac shock signal, not only the spectrum of the heartbeat frequency band will have a relatively large peak, but also the spectrum of the respiratory frequency band will also have a relatively large peak; but in the apnea In the frequency spectrum of the cardiac shock signal, only the frequency spectrum of the heartbeat frequency band has a peak with a large amplitude, and the frequency spectrum of the respiratory frequency band does not have a peak with a large amplitude. Therefore, in this embodiment, it may be determined whether the second apnea condition is met based on the frequency spectrum of the shock signal within the time period. Please refer to FIG. 6, before step S103, the apnea detection method may also include:
步骤S601，计算心跳频段对应的所有波峰的幅值的第一平均值，以及呼吸频段对应的所有波峰的幅值的第二平均值。Step S601, 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.
在本实施例中，电子设备在计算出心冲击信号在上述时间段内的频谱后，通过提取出心跳频段的频谱中的所有波峰的幅值，并进行平均值计算，可以得到心跳频段对应的所有波峰的幅值的平均值(即第一平均值)；通过提取出呼吸频段的频谱中的所有波峰的幅值，并进行平均值计算，可以得到呼吸频段对应的所有波峰的幅值的平均值(即第二平均值)。In this embodiment, after the electronic device calculates the frequency spectrum of the cardiac shock signal within the above time period, by extracting the amplitudes of all the peaks in the frequency spectrum of the heartbeat frequency band and performing average calculation, the corresponding heartbeat frequency band can be obtained. The average value of the amplitudes of all peaks (i.e. the first average value); by extracting the amplitudes of all peaks in the frequency spectrum of the respiratory frequency band and performing average calculation, the average value of the amplitudes of all peaks corresponding to the respiratory frequency band can be obtained value (i.e. the second mean value).
步骤S602，若心跳频段内存在至少一个波峰的幅值大于第一平均值的第一预设倍数，且呼吸频段内不存在波峰的幅值大于第二平均值的第二预设倍数，则判定频域特征符合第二呼吸暂停条件。Step S602, if there is at least one peak in the heartbeat frequency band whose amplitude is greater than the first preset multiple of the first average value, and there is no peak in the respiratory frequency band whose amplitude is greater than the second preset multiple of the second average value, then determine The frequency domain features meet the second apnea condition.
在本实施例中，第一预设倍数和第二预设倍数的具体取值可以根据实际情况设置，第一预设倍数和第二预设倍数可以相同，也可以不同。In this embodiment, the specific values of the first preset multiple and the second preset multiple can be set according to actual conditions, and the first preset multiple and the second preset multiple can be the same or different.
例如，第一预设倍数和第二预设倍数可以均为2，则当心跳频段内存在幅值大于第一平均值的2倍的波峰，呼吸频段内不存在幅值大于第二平均值的2倍的波峰，表明心冲击信号在该时间段内的频谱中，心跳频段有幅值较大的波峰，而呼吸频段没有幅值较大的波峰，该心冲击信号在频域上表现为呼吸暂停，故判定该心冲击信号的频域特征符合第二呼吸暂停条件。For example, the first preset multiple and the second preset multiple can both be 2, then when there is a peak whose amplitude is greater than twice the first average value in the heartbeat frequency band, there is no wave peak whose amplitude is greater than the second average value in the respiratory frequency band. 2 times the peak, indicating that in the frequency spectrum of the heart shock signal in this time period, the heartbeat frequency band has a peak with a larger amplitude, while the respiratory frequency band has no peak with a larger amplitude. Therefore, it is determined that the frequency domain feature of the shock signal meets the second apnea condition.
本发明实施例提供的呼吸暂停检测方法，通过计算心跳频段对应的所有波峰的幅值的第一平均值以及呼吸频段对应的所有波峰的幅值的第二平均值，基于第一平均值判定心跳频段内是否存在幅值较大的波峰，基于第二平均值判定呼吸频段内是否存在幅值较大的波峰，当心跳频段内存在幅值较大的波峰并且呼吸频段内不存在幅值较大的波峰时，判定心冲击信号的频域特征符合第二呼吸暂停条件，实现了对该时间段内的心冲击信号是否在频域上表现为呼吸暂停的准确判定，排除了心跳停止情况下对呼吸暂停判定的干扰。The apnea detection method provided by the embodiment of the present invention calculates 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 respiratory frequency band, and determines the heartbeat based on the first average value Whether there is a peak with a large amplitude in the frequency band, determine whether there is a peak with a large amplitude in the respiratory frequency band based on the second average value, when there is a peak with a large amplitude in the heartbeat frequency band and there is no peak with a large amplitude in the respiratory frequency band When the peak of the peak, it is determined that the frequency domain characteristics of the cardiac shock signal meet the second apnea condition, which realizes the accurate determination of whether the cardiac shock signal in the time period is represented as apnea in the frequency domain, and excludes the occurrence of apnea in the case of cardiac arrest. Interference in apnea determination.
为了执行上述实施例及各个可能的方式中的相应步骤，下面给出一种呼吸暂停检测装置的实现方式。请参阅图7，为本发明实施例提供的呼吸暂停检测装置700的一种功能模块图。需要说明的是，本实施例所提供的呼吸暂停检测装置700，其基本原理及产生的技术效果和上述实施例相同，为简要描述，本实施例部分未提及之处，可参考上述的实施例中相应内容。该呼吸暂停检测装置700包括时域特征提取模块710、频域特征提取模块720和呼吸检测模块730。In order to execute the corresponding steps in the foregoing embodiments and various possible manners, an implementation manner of an apnea detection device is given below. Please refer to FIG. 7 , which is a functional block diagram of an
该时域特征提取模块710，用于提取采集的心冲击信号的时域特征。The time-domain
可以理解，该时域特征提取模块710可以执行上述步骤S101。It can be understood that the time-domain
该频域特征提取模块720，用于在心冲击信号的时域特征符合第一呼吸暂停条件的情况下，提取心冲击信号的频域特征。The frequency domain
可以理解，该频域特征提取模块720可以执行上述步骤S102。It can be understood that the frequency domain
该呼吸检测模块730，用于在频域特征符合第二呼吸暂停条件的情况下，确定存在呼吸暂停事件。The
可以理解，该呼吸检测模块730可以执行上述步骤S103。It can be understood that the
可选地，该频域特征提取模块720可以用于确定心冲击信号的时域特征符合第一呼吸暂停条件的时间段，计算心冲击信号在时间段内的频谱；频域特征包括频谱中的每个波峰的幅值。Optionally, the frequency domain
其中，该频域特征提取模块720具体用于对时间段内的心冲击信号进行傅里叶变换，得到心冲击信号在时间段内的频谱。Wherein, the frequency domain
可以理解，该频域特征提取模块720可以执行上述步骤S1021～S1022。It can be understood that the frequency domain
可选地，请参照图8，该时域特征包括心冲击信号中的每个波峰的幅值，该呼吸暂停检测装置700还可以包括第一判定模块740和第二判定模块750。Optionally, please refer to FIG. 8 , the time-domain feature includes the amplitude of each peak in the cardiac shock signal, and the
该第一判定模块740，用于根据心冲击信号的时域特征确定多个连续的目标波峰；其中，多个连续的目标波峰中，第一个目标波峰的前一个波峰的幅值与第一个目标波峰的幅值之间的比值大于预设比值，最后一个目标波峰的后一个波峰的幅值与最后一个目标波峰的幅值之间的比值大于预设比值，并且任意相邻的两个目标波峰的幅值中，较大幅值与较小幅值之间的比值小于或等于预设比值；若多个连续的目标波峰的持续时间大于或等于预设时间，则判定心冲击信号的时域特征符合第一呼吸暂停条件。The
其中，该持续时间为第一个目标波峰所处的时刻与最后一个目标波峰所处的时刻之间的差值。Wherein, the duration is the difference between the moment of the first target peak and the moment of the last target peak.
在本实施例中，该第一判定模块740具体用于对心冲击信号进行寻峰，确定相邻的第一波峰和第二波峰，第一波峰对应的时刻小于第二波峰对应的时刻；若第一波峰的幅值与第二波峰的幅值之间的比值大于预设比值，则在第二波峰对应的时刻之后寻找连续的至少一个第三波峰，直到第三波峰的后一个波峰的幅值与第三波峰的幅值之间的比值大于预设比值；其中，每个第三波峰与相邻的前一个波峰的幅值中，较大幅值与较小幅值之间的比值小于或等于预设比值；将第二波峰以及所有的第三波峰均确定为目标波峰，得到多个连续的目标波峰。In this embodiment, the
可以理解，该第一判定模块740可以执行上述步骤S301～S302、子步骤S3011～S3013。It can be understood that the
在本实施例中，该频谱包括心跳频段的频谱和呼吸频段的频谱，该第二判定模块750，用于计算心跳频段对应的所有波峰的幅值的第一平均值，以及呼吸频段对应的所有波峰的幅值的第二平均值；若心跳频段内存在至少一个波峰的幅值大于第一平均值的第一预设倍数，且呼吸频段内不存在波峰的幅值大于第二平均值的第二预设倍数，则判定频域特征符合第二呼吸暂停条件。In this embodiment, the frequency spectrum includes the frequency spectrum of the heartbeat frequency band and the frequency spectrum of the respiratory frequency band. The second average value of the peak amplitude; if there is at least one peak amplitude greater than the first preset multiple of the first average value in the heartbeat frequency band, and there is no peak amplitude greater than the first preset multiple of the second average value in the respiratory frequency band Two preset multiples, it is determined that the frequency domain feature meets the second apnea condition.
可以理解，该第二判定模块750可以执行上述步骤S601～S602。It can be understood that the
本发明实施例提供的呼吸暂停检测装置700，通过时域特征提取模块710提取采集的心冲击信号的时域特征，频域特征提取模块720在心冲击信号的时域特征符合第一呼吸暂停条件的情况下，提取心冲击信号的频域特征，呼吸检测模块730在频域特征符合第二呼吸暂停条件的情况下，确定存在呼吸暂停事件。如此，实现了将心冲击信号的时域特征和频域特征相结合，在时域和频域中均判定呼吸暂停的情况下确定存在呼吸暂停事件，有效减少了误判的可能，提高了呼吸暂停事件的检测准确性；整个检测过程无需人为判别，心冲击信号的获取也比较方便，提高了呼吸暂停事件检测的操作简便性。The
本发明实施例提供的电子设备可以包括处理器和存储器，该存储器中存储有至少一条指令、至少一段程序、代码集或指令集，该至少一条指令、该至少一段程序、该代码集或指令集由该处理器加载并执行以实现如上述方法实施例所提供的呼吸暂停检测方法。The electronic equipment provided by the embodiments of the present invention may include a processor and a memory, at least one instruction, at least one program, code set or instruction set are stored in the memory, the at least one instruction, the at least one program, the code set or the instruction set It is loaded and executed by the processor to realize the apnea detection method provided by the above method embodiment.
该存储器可用于存储软件程序以及模块，处理器通过运行存储在存储器的软件程序以及模块，从而执行各种功能应用以及数据处理。存储器可主要包括存储程序区和存储数据区，其中，存储程序区可存储操作系统、功能所需的应用程序等；存储数据区可存储根据所述设备的使用所创建的数据等。此外，存储器可以包括高速随机存取存储器，还可以包括非易失性存储器，例如至少一个磁盘存储器件、闪存器件、或其他易失性固态存储器件。相应地，存储器还可以包括存储器控制器，以提供处理器对存储器的访问。The memory can be used to store software programs and modules, and the processor executes various functional applications and data processing by running the software programs and modules stored in the memory. The memory may mainly include a program storage area and a data storage area, wherein the program storage area may store operating systems, application programs required by functions, etc.; the data storage area may store data created according to the use of the device, etc. In addition, 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 devices. Correspondingly, the memory may also include a memory controller to provide processor access to the memory.
图9是本发明实施例提供的电子设备100的一种硬件结构示意图。如图9所示，该电子设备100可因配置或性能不同而产生比较大的差异，可以包括一个或一个以上处理器(ProcessingUnits，CPU)111(处理器111可以包括但不限于微处理器MCU或可编程逻辑器件FPGA等的处理装置)、用于存储数据的存储器113，一个或一个以上存储应用程序1123或数据1122的存储介质112(例如一个或一个以上海量存储设备)。其中，存储器113和存储介质112可以是短暂存储或持久存储。存储在存储介质112的程序可以包括一个或一个以上模块，每个模块可以包括对电子设备100中的一系列指令操作。更进一步地，处理器111可以设置为与存储介质112通信，在电子设备100上执行存储介质112中的一系列指令操作。电子设备100还可以包括一个或一个以上电源116，一个或一个以上有线或无线网络接口115，一个或一个以上输入输出接口114，和/或，一个或一个以上操作系统1121，例如WindowsServerTM，MacOSXTM，UnixTM,LinuxTM，FreeBSDTM等。FIG. 9 is a schematic diagram of a hardware structure of an
输入输出接口114可以用于经由一个网络接收或者发送数据。上述的网络具体实例可包括电子设备100的通信供应商提供的无线网络。在一个实例中，输入输出接口114包括一个网络适配器(NetworkInterfaceController，NIC)，其可通过基站与其他网络设备相连从而可与互联网进行通讯。在一个实例中，输入输出接口114可以为射频(RadioFrequency，RF)模块，其用于通过无线方式与互联网进行通讯。The input-
本领域普通技术人员可以理解，图9所示的结构仅为示意，其并不对上述电子设备100的结构造成限定。例如，电子设备100还可包括比图9中所示更多或者更少的组件，或者具有与图9所示不同的配置。Those of ordinary skill in the art can understand that the structure shown in FIG. 9 is only a schematic diagram, and it does not limit the structure of the above-mentioned
本发明实施例提供的电子设备，通过提取采集的心冲击信号的时域特征，在心冲击信号的时域特征符合第一呼吸暂停条件的情况下，提取该心冲击信号的频域特征，在该频域特征符合第二呼吸暂停条件的情况下，确定存在呼吸暂停事件。如此，实现了将心冲击信号的时域特征和频域特征相结合，在时域和频域中均判定呼吸暂停的情况下确定存在呼吸暂停事件，有效减少了误判的可能，提高了呼吸暂停事件的检测准确性；整个检测过程无需人为判别，心冲击信号的获取也比较方便，提高了呼吸暂停事件检测的操作简便性。The electronic device provided by the embodiment of the present invention extracts the time-domain feature of the collected cardiac shock signal, and extracts the frequency-domain feature of the cardiac shock signal when the time-domain feature of the cardiac shock signal meets the first apnea condition. If the frequency domain feature meets the second apnea condition, it is determined that there is an apnea event. In this way, the combination of time-domain features and frequency-domain features of cardiac shock signals is realized, and the existence of apnea events is determined when apnea is judged in both the time-domain and frequency-domain, which effectively reduces the possibility of misjudgment and improves respiratory rate. The detection accuracy of the pause event; the whole detection process does not need to be judged manually, and the acquisition of the cardiac shock signal is also relatively convenient, which improves the operation simplicity of the detection of the pause event.
本发明实施例还提供一种计算机可读存储介质，该计算机可读存储介质上存储有计算机程序，该计算机程序被处理器执行时实现上述呼吸暂停检测方法实施例的各个过程，且能达到相同的技术效果，为避免重复，这里不再赘述。其中，所述的计算机可读存储介质，如只读存储器(Read-OnlyMemory，简称ROM)、随机存取存储器(RandomAccessMemory，简称RAM)、磁碟或者光盘等。An embodiment of the present invention also provides a computer-readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, the various processes of the above-mentioned embodiment of the apnea detection method can be realized, and the same can be achieved. To avoid repetition, the technical effects will not be repeated here. Wherein, the computer-readable storage medium is, for example, a read-only memory (Read-Only Memory, ROM for short), a random access memory (Random Access Memory, RAM for short), a magnetic disk or an optical disk, and the like.
在本申请所提供的几个实施例中，应该理解到，所揭露的装置和方法，也可以通过其它的方式实现。以上所描述的装置实施例仅仅是示意性的，例如，附图中的流程图和框图显示了根据本发明的多个实施例的装置、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上，流程图或框图中的每个方框可以代表一个模块、程序段或代码的一部分，所述模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意，在有些作为替换的实现方式中，方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如，两个连续的方框实际上可以基本并行地执行，它们有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意的是，框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合，可以用执行规定的功能或动作的专用的基于硬件的系统来实现，或者可以用专用硬件与计算机指令的组合来实现。In the several embodiments provided in this application, it should be understood that the disclosed devices and methods may also be implemented in other ways. The device embodiments described above are only illustrative. For example, the flowcharts and block diagrams in the accompanying drawings show the architecture, functions and possible implementations of devices, methods and computer program products according to multiple embodiments of the present invention. operate. In this regard, each block in a flowchart or block diagram may represent a module, program segment, or part of code that includes one or more Executable instructions. 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 in succession may, in fact, be executed substantially concurrently, or they may sometimes be executed in the reverse order, depending upon the functionality involved. It should also be noted that each block of the block diagrams and/or flowchart illustrations, and combinations of blocks in the block diagrams and/or flowchart illustrations, can be implemented by a dedicated hardware-based system that performs the specified function or action , or may be implemented by a combination of dedicated hardware and computer instructions.
另外，在本发明各个实施例中的各功能模块可以集成在一起形成一个独立的部分，也可以是各个模块单独存在，也可以两个或两个以上模块集成形成一个独立的部分。In addition, each functional module in each embodiment of the present invention can be integrated together to form an independent part, or each module can exist independently, or two or more modules can be integrated to form an independent part.
所述功能如果以软件功能模块的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。If the functions are implemented in the form of software function modules and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, the essence of the technical solution of the present invention or the part that contributes to the prior art or the part of the technical solution can be embodied in the form of a software product, and the computer software product is stored in a storage medium, including Several instructions are used to make a computer device (which may be a personal computer, a server, or a network device, etc.) execute all or part of the steps of the methods described in various embodiments of the present invention. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), magnetic disk or optical disk and other media that can store program codes. .
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。The above descriptions are only preferred embodiments of the present invention, and are not intended to limit the present invention. For those skilled in the art, the present invention may have various modifications and changes. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principles of the present invention shall be included within the protection scope of the present invention.
Claims (10)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110761417.XACN115568844B (en) | 2021-07-06 | 2021-07-06 | Apnea detection method, device, electronic device and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110761417.XACN115568844B (en) | 2021-07-06 | 2021-07-06 | Apnea detection method, device, electronic device and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115568844Atrue CN115568844A (en) | 2023-01-06 |
| CN115568844B CN115568844B (en) | 2025-08-08 |
Family
ID=84579299
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110761417.XAActiveCN115568844B (en) | 2021-07-06 | 2021-07-06 | Apnea detection method, device, electronic device and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115568844B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116172527A (en)* | 2023-03-15 | 2023-05-30 | 康原(江苏)科技有限公司 | A method, system and device for intelligent health monitoring |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4365636A (en)* | 1981-06-19 | 1982-12-28 | Medicon, Inc. | Method of monitoring patient respiration and predicting apnea therefrom |
| US20070118054A1 (en)* | 2005-11-01 | 2007-05-24 | Earlysense Ltd. | Methods and systems for monitoring patients for clinical episodes |
| WO2018134823A1 (en)* | 2017-01-20 | 2018-07-26 | Oridion Medical 1987 Ltd. | Dual sensor |
| CN109620231A (en)* | 2018-12-25 | 2019-04-16 | 湖南明康中锦医疗科技发展有限公司 | Flow limitation determination method, device, computer equipment and storage medium |
| CN111134637A (en)* | 2019-12-06 | 2020-05-12 | 深圳绿米联创科技有限公司 | Signal processing method, signal processing device, electronic equipment and storage medium |
- 2021
- 2021-07-06CNCN202110761417.XApatent/CN115568844B/enactiveActive
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4365636A (en)* | 1981-06-19 | 1982-12-28 | Medicon, Inc. | Method of monitoring patient respiration and predicting apnea therefrom |
| US20070118054A1 (en)* | 2005-11-01 | 2007-05-24 | Earlysense Ltd. | Methods and systems for monitoring patients for clinical episodes |
| WO2018134823A1 (en)* | 2017-01-20 | 2018-07-26 | Oridion Medical 1987 Ltd. | Dual sensor |
| CN109620231A (en)* | 2018-12-25 | 2019-04-16 | 湖南明康中锦医疗科技发展有限公司 | Flow limitation determination method, device, computer equipment and storage medium |
| CN111134637A (en)* | 2019-12-06 | 2020-05-12 | 深圳绿米联创科技有限公司 | Signal processing method, signal processing device, electronic equipment and storage medium |
Non-Patent Citations (3)
| Title |
|---|
| A NEW ACOUSTIC METHOD OF DIFFERENTIATING PALATAL FROM NON-PALATAL SNORING OSBORNE, JE ET AL: "A new acoustic method of differentiating palatal from non-palatal snoring", 《BLACKWELL PUBLISHING》, 30 April 1999 (1999-04-30)* |
| 王海涛;郑慧君;曹征涛;杨军;俞梦孙;: "考虑个体特征的非脑电睡眠分期", 中国生物医学工程学报, no. 02, 20 April 2010 (2010-04-20)* |
| 许益彬: "基于非接触式睡眠监测系统的呼吸暂停综合征检测研究", 《硕士电子期刊出版信息》, 15 September 2019 (2019-09-15)* |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116172527A (en)* | 2023-03-15 | 2023-05-30 | 康原(江苏)科技有限公司 | A method, system and device for intelligent health monitoring |
| CN116172527B (en)* | 2023-03-15 | 2023-09-22 | 康原(江苏)科技有限公司 | An intelligent health monitoring method, system and device |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115568844B (en) | 2025-08-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CA2805324C (en) | Feature characterization for breathing monitor | |
| US9814438B2 (en) | Methods and apparatus for performing dynamic respiratory classification and tracking | |
| CN110809429A (en) | System and method for filtering noise and analyzing venous waveform signals | |
| CN106687033B (en) | Heartbeat detection method and heartbeat detection device | |
| JP6687645B2 (en) | Biological signal processing method and biological signal processing apparatus | |
| CN108992053B (en) | A real-time, untethered method for detecting heart rate and heartbeat interval | |
| WO2020072434A1 (en) | Methods and systems for monitoring sleep apnea | |
| KR20200101385A (en) | Systems and methods for filtering medical device noise artifacts from venous waveform signals | |
| WO2017118127A1 (en) | Heartbeat signal processing method, device and system | |
| CN106709414B (en) | Physiological signal measuring system and method thereof | |
| CN114615926A (en) | Snore detection device and method based on sound analysis | |
| CN115568844A (en) | Apnea detection method and device, electronic equipment and storage medium | |
| JP7466548B2 (en) | Apparatus and system for generating a respiratory alarm | |
| JP6315633B2 (en) | Heart rate detection method and heart rate detection device | |
| CN107811610B (en) | Respiration rate detection method and device, electronic equipment and storage medium | |
| CN114222524B (en) | Systems for processing respiratory rate | |
| US11375908B2 (en) | Blood pressure detection signal sampling and compensation method and apparatus, and blood pressure signal collection system | |
| US20180333064A1 (en) | Respiration estimation method and apparatus | |
| EP3695776A1 (en) | A method and system for generating a respiration instability signal | |
| CN113616217B (en) | Method and device for generating baseline drift curve | |
| US11311243B2 (en) | Blood pressure data processing apparatus, blood pressure data processing method, and blood pressure data processing program | |
| CN115005803A (en) | Breath detection method and device | |
| CN104605886B (en) | Stridulate sound detection device and method | |
| CN115517632B (en) | Resonance breathing frequency measurement method, interactive prompt generation method, device and equipment | |
| JP2020092782A (en) | Respiratory cycle measuring device and respiratory cycle measuring program |
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 |