CN114176521B - Sleep staging method and device based on radar information and terminal - Google Patents

Sleep staging method and device based on radar information and terminal Download PDF

Info

Publication number
CN114176521B
CN114176521B CN202111498745.1A CN202111498745A CN114176521B CN 114176521 B CN114176521 B CN 114176521B CN 202111498745 A CN202111498745 A CN 202111498745A CN 114176521 B CN114176521 B CN 114176521B
Authority
CN
China
Prior art keywords
sleep
heartbeat
index
amplitude
moment
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.)
Active
Application number
CN202111498745.1A
Other languages
Chinese (zh)
Other versions
CN114176521A (en
Inventor
程毅
彭诚诚
胡承帅
赵洛伟
司孟昌
胡倩婷
何文彦
刘子华
秦屹
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Whst Co Ltd
Original Assignee
Whst Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Whst Co Ltd filed Critical Whst Co Ltd
Priority to CN202111498745.1A priority Critical patent/CN114176521B/en
Publication of CN114176521A publication Critical patent/CN114176521A/en
Application granted granted Critical
Publication of CN114176521B publication Critical patent/CN114176521B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/48Other medical applications
    • A61B5/4806Sleep evaluation
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/0205Simultaneously evaluating both cardiovascular conditions and different types of body conditions, e.g. heart and respiratory condition
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/024Detecting, measuring or recording pulse rate or heart rate
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/0507Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves  using microwaves or terahertz waves
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/08Detecting, measuring or recording devices for evaluating the respiratory organs
    • A61B5/0816Measuring devices for examining respiratory frequency
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/48Other medical applications
    • A61B5/4806Sleep evaluation
    • A61B5/4812Detecting sleep stages or cycles
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Physics & Mathematics (AREA)
  • Public Health (AREA)
  • Pathology (AREA)
  • Physiology (AREA)
  • Cardiology (AREA)
  • Pulmonology (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Psychiatry (AREA)
  • Signal Processing (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)

Abstract

The invention provides a sleep staging method, a sleep staging device, a sleep staging terminal and a sleep staging storage medium based on radar information. The method comprises the following steps: acquiring respiratory heartbeat information in a sleep period monitored by a vital sign monitoring radar; calculating a sleep staging index at each moment in a sleep period based on the breathing heartbeat information; calculating a sleep stage characteristic value of each moment in a sleep period based on the sleep stage index; and comparing the sleep stage characteristic value with a preset stage threshold value, and determining the sleep stage corresponding to each moment in the sleep period according to the comparison result. The invention carries out sleep staging based on radar information, acquires information under the conditions of no perception and no influence on normal sleep of a user, has simple and accurate sleep staging processing process and has higher reliability.

Description

Sleep staging method and device based on radar information and terminal
Technical Field
The invention relates to the technical field of health monitoring, in particular to a sleep staging method and device based on radar information, a terminal and a storage medium.
Background
Currently, in some health monitoring systems, it is necessary to comprehensively analyze the health condition of a human body and output a health report, where sleep monitoring for sleep condition is an important part of the health report, and what is more important in sleep monitoring is how to accurately perform sleep staging. Generally, the sleep process can be divided into a wake-up period, a light sleep period and a deep sleep period.
In the prior art, sleep staging is generally performed by using physiological state signals such as electroencephalogram signals and electrooculogram signals, and for example, patent document No. CN113303770A is disclosed. However, physiological state signals such as electroencephalogram signals and electrooculogram signals are often monitored by some contact-type devices, for example, patch-type or wearable contact-type monitoring devices, which are easy to press human tissues for a long time and damage human bodies, and when some sensitive people use patch-type or wearable sleep monitoring instruments, the phenomenon of insomnia caused by tension is easy to occur, and it is difficult to effectively monitor the real sleep condition of the detected people.
For example, patent document No. CN113456030A discloses a sleep staging method based on heart rate monitoring data, however, this method needs to extract 8 time domain characteristic indexes and 4 frequency domain characteristic indexes from heart rate data, and needs to construct a bidirectional gate cycle network for sleep staging, which is complex in processing process, and reliability of staging results is not high enough due to single sleep staging based on heart rate.
Disclosure of Invention
The invention provides a sleep staging method, a sleep staging device, a sleep staging terminal and a storage medium based on radar information, and aims to solve the problems that the sleep staging method in the prior art is complex in processing process and low in reliability.
In a first aspect, the present invention provides a sleep staging method based on radar information, which is applied to a health monitoring system, wherein the health monitoring system includes a vital sign monitoring radar, and the method includes:
acquiring respiratory heartbeat information in a sleep period monitored by a vital sign monitoring radar;
calculating a sleep staging index of each moment in a sleep period based on the breathing heartbeat information;
calculating a sleep stage characteristic value of each moment in a sleep period based on the sleep stage index;
and comparing the sleep stage characteristic value with a preset stage threshold value, and determining the sleep stage corresponding to each moment in the sleep period according to the comparison result.
In a possible implementation manner, the respiratory heartbeat information includes: the respiratory frequency, heartbeat frequency, respiratory amplitude and heartbeat amplitude of each moment in the sleep period;
the sleep staging index comprises: a respiratory frequency stability index, a heartbeat frequency stability index, a respiratory amplitude variation index and a heartbeat amplitude variation index;
correspondingly, the calculating the sleep staging index at each moment in the sleep session based on the respiratory heartbeat information includes:
calculating the average value of the respiratory frequency of the set time length according to the respiratory frequency of each time in the previous set time length of the target time length, and calculating the respiratory frequency stability index of the target time length according to the respiratory frequency of each time in the previous set time length of the target time length and the average value of the respiratory frequency of the set time length;
calculating the heartbeat frequency average value of the set time length according to the heartbeat frequency of each time in the previous set time length of the target time length, and calculating the heartbeat frequency stability index of the target time length according to the heartbeat frequency of each time in the previous set time length of the target time length and the heartbeat frequency average value of the set time length;
calculating the average value of the breathing amplitude of the set time length according to the breathing amplitude of each time in the previous set time length of the target time, and calculating the breathing amplitude change index of the target time according to the breathing amplitude of each time in the previous set time length of the target time and the average value of the breathing amplitude of the set time length;
calculating the average value of the heartbeat amplitude of the set time according to the heartbeat amplitude of each time in the previous set time of the target time, and calculating the heartbeat amplitude change index of the target time according to the heartbeat amplitude of each time in the previous set time of the target time and the average value of the heartbeat amplitude of the set time.
In one possible implementation, the method of calculating the respiratory rate stability index at a target time comprises:
Figure 658138DEST_PATH_IMAGE001
wherein the content of the first and second substances,Twhich represents the time of the object and the time of day,Lwhich indicates the set time period of time,
Figure 542656DEST_PATH_IMAGE002
an index representing the stability of the breathing frequency at the target moment,
Figure DEST_PATH_IMAGE003
to representtThe breathing frequency at the moment of time is,
Figure 791234DEST_PATH_IMAGE004
represents the mean value of the breathing frequency at the target moment; wherein the content of the first and second substances,
Figure DEST_PATH_IMAGE005
the method for calculating the heart beat frequency stability index at the target moment comprises the following steps:
Figure 858547DEST_PATH_IMAGE006
wherein the content of the first and second substances,
Figure DEST_PATH_IMAGE007
an index representing the stability of the frequency of the heartbeat at the target time,
Figure 141761DEST_PATH_IMAGE008
to representtThe frequency of the heart beats at a moment in time,
Figure DEST_PATH_IMAGE009
representing a mean value of the heartbeat frequency at the target time; wherein the content of the first and second substances,
Figure 698645DEST_PATH_IMAGE010
the method for calculating the respiratory amplitude variation index at the target moment comprises the following steps:
Figure DEST_PATH_IMAGE011
wherein the content of the first and second substances,
Figure 434519DEST_PATH_IMAGE012
an index representing the change in the amplitude of breathing at the target moment,
Figure DEST_PATH_IMAGE013
representing the amplitude of the breathing at the target moment,
Figure 102261DEST_PATH_IMAGE014
indicating the front of the target timeLMean value of breathing amplitude over a time period;
the method for calculating the heart beat amplitude change index at the target moment comprises the following steps:
Figure DEST_PATH_IMAGE015
wherein the content of the first and second substances,
Figure 944709DEST_PATH_IMAGE016
an index of change in the amplitude of the heartbeat at the target time,
Figure DEST_PATH_IMAGE017
representing the amplitude of the heartbeat at the target time,
Figure 406914DEST_PATH_IMAGE018
indicating the front of the target timeLMean value of the amplitude of the heartbeat over the duration.
In a possible implementation manner, the calculating a sleep stage characteristic value of each time in a sleep period based on the sleep stage index includes:
normalizing the sleep staging indexes;
carrying out weighted summation on the normalized sleep stage indexes based on the preset weight of each sleep stage index to obtain a sleep stage characteristic value of each moment in a sleep period;
the weights of the respiratory frequency stability index and the respiratory amplitude change index are both larger than the heartbeat frequency stability index and the heartbeat amplitude change index.
In a possible implementation manner, the normalizing the sleep staging index includes;
acquiring a maximum value of a respiratory frequency stability index, a maximum value of a heartbeat frequency stability index, a maximum value of a respiratory amplitude change index and a maximum value of a heartbeat amplitude change index in a sleep period;
taking the ratio of the respiratory frequency stability index at each moment to the maximum value of the respiratory frequency stability index as the normalized respiratory frequency stability index at the moment;
taking the ratio of the heartbeat frequency stability index at each moment to the maximum value of the heartbeat frequency stability index as the normalized heartbeat frequency stability index at the moment;
taking the ratio of the respiration amplitude change index of each moment to the maximum value of the respiration amplitude change index as the normalized respiration amplitude change index of the moment;
and taking the ratio of the heartbeat amplitude change index of each moment to the maximum value of the heartbeat amplitude change index as the normalized heartbeat amplitude change index of the moment.
In a possible implementation manner, the staging threshold includes a wakefulness threshold and a shallow sleep threshold, the comparing the sleep staging characteristic value with a preset staging threshold, and determining the sleep period corresponding to each time in the sleep interval according to the comparison result includes:
if the sleep stage characteristic value of the target time is greater than or equal to a preset waking threshold, the sleep stage corresponding to the time is a waking stage;
if the sleep stage characteristic value of the target time is smaller than the waking threshold and is greater than or equal to a preset light sleep threshold, the sleep stage corresponding to the time is a light sleep stage;
and if the sleep stage characteristic value of the target time is smaller than the light sleep threshold, the sleep stage corresponding to the time is a deep sleep stage.
In a possible implementation manner, before the obtaining of the respiratory heartbeat information in the sleep period monitored by the vital sign monitoring radar, the method further includes:
acquiring in-bed information and rising information monitored by a vital sign monitoring radar;
determining the getting-on time and the getting-up time according to the in-bed information and the getting-up information;
and determining the sleep time interval according to the time of getting up and the time of getting up.
In a second aspect, the present invention provides a sleep staging device based on radar information, which is applied to a health monitoring system, wherein the health monitoring system includes a vital sign monitoring radar, and the sleep staging device includes:
the first acquisition unit is used for acquiring the respiration heartbeat information in the sleep period monitored by the vital sign monitoring radar;
the staging index calculating unit is used for calculating the sleep staging index of each moment in the sleep period based on the respiratory heartbeat information acquired by the first acquiring unit;
the staging characteristic value calculating unit is used for calculating the sleep staging characteristic value of each moment in the sleep period based on the sleep staging index calculated by the staging index calculating unit;
and the sleep period determining unit is used for comparing the sleep stage characteristic value calculated by the stage characteristic value calculating unit with a preset stage threshold value and determining the sleep period corresponding to each moment in the sleep period according to the comparison result.
In a possible implementation manner, the respiratory heartbeat information includes: the respiratory frequency, the heartbeat frequency, the respiratory amplitude and the heartbeat amplitude of each moment in the sleep period; the sleep staging index comprises: a respiratory frequency stability index, a heartbeat frequency stability index, a respiratory amplitude variation index and a heartbeat amplitude variation index;
correspondingly, the staging index calculating unit is specifically configured to calculate a respiratory frequency average value of the set time according to the respiratory frequency of each time in the previous set time of the target time, and calculate a respiratory frequency stability index of the target time according to the respiratory frequency of each time in the previous set time of the target time and the respiratory frequency average value of the set time;
calculating the heartbeat frequency average value of the set time length according to the heartbeat frequency of each time in the previous set time length of the target time length, and calculating the heartbeat frequency stability index of the target time length according to the heartbeat frequency of each time in the previous set time length of the target time length and the heartbeat frequency average value of the set time length;
calculating the average value of the breathing amplitude of the set time length according to the breathing amplitude of each time in the previous set time length of the target time, and calculating the breathing amplitude change index of the target time according to the breathing amplitude of each time in the previous set time length of the target time and the average value of the breathing amplitude of the set time length;
calculating the average value of the heartbeat amplitude of the set time according to the heartbeat amplitude of each time in the previous set time of the target time, and calculating the heartbeat amplitude change index of the target time according to the heartbeat amplitude of each time in the previous set time of the target time and the average value of the heartbeat amplitude of the set time.
In a possible implementation manner, the method for calculating the respiratory rate stability index at the target moment by the staging index calculation unit specifically includes:
Figure 426822DEST_PATH_IMAGE001
wherein the content of the first and second substances,Twhich represents the time of the object and the time of the object,Lit indicates that the set time period is,
Figure 632676DEST_PATH_IMAGE002
an index representing the stability of the breathing frequency at the target moment,
Figure 359323DEST_PATH_IMAGE003
to representtThe breathing frequency at the moment of time is,
Figure 54747DEST_PATH_IMAGE004
represents the mean value of the breathing frequency at the target moment; wherein the content of the first and second substances,
Figure 765214DEST_PATH_IMAGE005
the method for calculating the heart beat frequency stability index at the target moment comprises the following steps:
Figure 243600DEST_PATH_IMAGE006
wherein the content of the first and second substances,
Figure 152650DEST_PATH_IMAGE007
an index representing the stability of the frequency of the heartbeat at the target time,
Figure 691079DEST_PATH_IMAGE008
to representtThe frequency of the heart beats at a moment in time,
Figure 685580DEST_PATH_IMAGE009
representing a mean value of the heartbeat frequency at the target time; wherein the content of the first and second substances,
Figure 466192DEST_PATH_IMAGE010
the method for calculating the respiratory amplitude variation index at the target moment comprises the following steps:
Figure 229748DEST_PATH_IMAGE011
wherein the content of the first and second substances,
Figure 204657DEST_PATH_IMAGE012
an index representing the change in the amplitude of breathing at the target moment,
Figure 624137DEST_PATH_IMAGE013
representing the amplitude of breathing at a target moment,
Figure 772222DEST_PATH_IMAGE014
Indicating the front of the target timeLMean value of breathing amplitude over a time period;
the method for calculating the heart beat amplitude change index at the target moment comprises the following steps:
Figure 327968DEST_PATH_IMAGE015
wherein the content of the first and second substances,
Figure 739358DEST_PATH_IMAGE016
an index of change in the amplitude of the heartbeat at the target time,
Figure 442872DEST_PATH_IMAGE017
representing the amplitude of the heartbeat at the target time,
Figure 66751DEST_PATH_IMAGE018
indicating the front of the target timeLMean value of the amplitude of the heartbeat over the duration.
In one possible implementation manner, the sleep staging apparatus further includes:
a normalization processing unit, configured to perform normalization processing on the sleep stage index;
correspondingly, the stage characteristic value calculating unit is specifically configured to perform weighted summation on the normalized sleep stage indexes based on preset weights of the sleep stage indexes to obtain a sleep stage characteristic value at each moment in a sleep period;
the weights of the respiratory frequency stability index and the respiratory amplitude change index are both larger than the heartbeat frequency stability index and the heartbeat amplitude change index.
In one possible implementation manner, the sleep staging apparatus further includes:
the second acquisition unit is used for acquiring the maximum value of the respiratory frequency stability index, the maximum value of the heartbeat frequency stability index, the maximum value of the respiratory amplitude change index and the maximum value of the heartbeat amplitude change index in a sleep period;
correspondingly, the staging characteristic value calculating unit is specifically configured to use a ratio of the respiratory frequency stability index at each time to a maximum value of the respiratory frequency stability index acquired by the second acquiring unit as the normalized respiratory frequency stability index at the time;
taking the ratio of the heartbeat frequency stability index at each moment to the maximum value of the heartbeat frequency stability index acquired by the second acquisition unit as the normalized heartbeat frequency stability index at the moment;
the ratio of the respiration amplitude change index at each moment to the maximum value of the respiration amplitude change index acquired by the second acquisition unit is used as the normalized respiration amplitude change index at the moment;
and taking the ratio of the heartbeat amplitude change index at each moment to the maximum value of the heartbeat amplitude change index acquired by the second acquisition unit as the normalized heartbeat amplitude change index at the moment.
In a possible implementation manner, the staging threshold includes a wakefulness threshold and a light sleep threshold, and the sleep stage determining unit is specifically configured to, if the sleep stage feature value of the target time is greater than or equal to a preset wakefulness threshold, determine that the sleep stage corresponding to the target time is a wakefulness stage;
if the sleep stage characteristic value of the target time is smaller than the waking threshold and is greater than or equal to a preset light sleep threshold, the sleep stage corresponding to the time is a light sleep stage;
if the sleep stage characteristic value of the target moment is smaller than the light sleep threshold value, the sleep stage corresponding to the moment is a deep sleep stage.
In one possible implementation manner, the sleep staging apparatus further includes:
a third obtaining unit, configured to obtain in-bed information and rising information monitored by the vital sign monitoring radar before obtaining the respiratory heartbeat information during the sleep period monitored by the vital sign monitoring radar;
the starting time determining unit is used for determining the getting-on time and the getting-up time according to the in-bed information and the getting-up information acquired by the third acquiring unit;
and the sleep period determining unit is used for determining the sleep period according to the getting-up time and the getting-up time determined by the starting time determining unit.
In a third aspect, the present invention provides a terminal, comprising a memory, a processor and a computer program stored in the memory and executable on the processor, wherein the processor, when executing the computer program, implements the steps of the radar information based sleep staging method according to the first aspect or any one of the possible implementations of the first aspect.
In a fourth aspect, the present invention provides a computer-readable storage medium storing a computer program which, when executed by a processor, implements the steps of the radar information based sleep staging method according to the first aspect as such or according to any one of the possible implementations of the first aspect.
The invention provides a sleep stage method, a device, a terminal and a storage medium based on radar information. On one hand, the monitoring information for sleep staging is acquired based on radar, and belongs to non-contact monitoring, so that a user has no perception and normal sleep of the user is not influenced; on the other hand, compared with the method that the sleep state of a human body can be more comprehensively reflected by relying on heart rate information based on the respiratory heartbeat information, the reliability of sleep staging is higher, in addition, the sleep staging index of each moment can be calculated, the sleep staging characteristic value of the moment can be determined, the sleep stage can be measured in seconds (each moment), and the sleep stage can be distinguished, so that the processing process is simple, convenient and accurate.
Drawings
In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the drawings needed to be used in the embodiments or the prior art descriptions will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without inventive exercise.
FIG. 1 is a flowchart of an implementation of a sleep staging method based on radar information according to an embodiment of the present invention;
FIG. 2 is a schematic structural diagram of a sleep staging apparatus based on radar information according to an embodiment of the present invention;
fig. 3 is a schematic diagram of a terminal according to an embodiment of the present invention.
Detailed Description
In the following description, for purposes of explanation and not limitation, specific details are set forth, such as particular system structures, techniques, etc. in order to provide a thorough understanding of the embodiments of the invention. It will be apparent, however, to one skilled in the art that the present invention may be practiced in other embodiments that depart from these specific details. In other instances, detailed descriptions of well-known systems, devices, circuits, and methods are omitted so as not to obscure the description of the present invention with unnecessary detail.
In order to make the objects, technical solutions and advantages of the present invention more apparent, the following description is made by way of specific embodiments with reference to the accompanying drawings.
The sleep staging method based on radar information provided by the embodiment of the invention can be applied to a health monitoring system, and the health monitoring system comprises a vital sign monitoring radar capable of monitoring human vital signs, such as the vital sign monitoring radar with the model number of Senstein IRS 60-5.
The monitored human vital sign signals may include respiration information, heartbeat information, person in bed information, person rising information, and the like. The embodiment of the invention carries out sleep staging based on the radar information monitored by the vital sign monitoring radar, a user does not sense the sleep staging in the monitoring process, the normal sleep of the user is not influenced, and the obtained radar information is accurate and reliable.
Referring to fig. 1, it shows a flowchart of an implementation of a sleep staging method based on radar information according to an embodiment of the present invention. As shown in fig. 1, a sleep staging method based on radar information may include:
in step 101, respiration and heartbeat information during a sleep period monitored by a vital sign monitoring radar is acquired.
In the embodiment of the invention, the respiration and heartbeat information in the sleep period monitored by the vital sign monitoring radar is firstly acquired, wherein the sleep period refers to the period from getting up to getting up of a user. The respiratory heartbeat information may include monitored heartbeat information and respiratory information of the user.
Specifically, the respiratory heartbeat information may include: respiratory rate, heartbeat rate, respiratory amplitude, and heartbeat amplitude at each moment in the sleep session. Wherein each time may be measured in seconds.
The sleep time interval can be manually set, and can also be determined according to the in-bed information and the rising information of the personnel monitored by the vital sign monitoring radar.
In an optional embodiment, before the step 101, the method may further include:
acquiring in-bed information and rising information monitored by a vital sign monitoring radar;
determining the getting-on time and the getting-up time according to the in-bed information and the getting-up information;
and determining a sleep period according to the time of getting up and the time of getting up.
The sleep period is determined according to the in-bed information and the rising information monitored by the vital sign monitoring radar.
Specifically, the radar can analyze the time of getting on and off the bed of the human body according to the in-bed mark (for example, 1 in the bed and 0 out of the bed) at each moment. The specific method for obtaining the time to bed is to find all the times tonbed (K) with the bed mark 1, where K =1,2 … K1. K1 is the number at the time when the bed flag is 1. Then starting from K =2 and traversing backwards to K = K1, calculating tone (K) -tone (K-1), if a preset threshold thresh1 is exceeded, ending the traversal, and recording the bed time as tone (K). The specific method of obtaining the time to bed exit is to find all times toffbed (K) at which the bed flag is 0, where K =1,2 … K2. K2 is the number at the time when the bed flag is 0. Then starting a forward traversal from K = K2 to K =2, computing tOffBed (K) -tOffBed (K-1), and if a preset threshold thresh2 is exceeded, ending the traversal, noting the bed exit time getUpTime as tofbd (K).
In step 102, a sleep staging index for each moment in a sleep session is calculated based on the respiratory heartbeat information.
In the embodiment of the present invention, the sleep staging index is a measure for defining characteristics of each sleep stage, and may include: a respiratory frequency stability index, a heartbeat frequency stability index, a respiratory amplitude variation index and a heartbeat amplitude variation index; these indices may be calculated from the respiratory heartbeat information monitored by the vital signs monitoring radar.
In the embodiment of the invention, the respiratory frequency stability index is used for representing the respiratory stability from the perspective of the speed of body movement (body chest micromotion caused by respiration), and the larger the respiratory frequency stability index is, the more unstable the respiration of a person at the current moment is, the higher the possibility that the person belongs to the shallow sleep period at the current moment is; the smaller the respiratory rate stability index, the more stable the respiration of the person at the current moment, the higher the probability that the person is in deep sleep at the current moment.
In the embodiment of the invention, the heartbeat frequency stability index is used for representing the stability of the heartbeat from the angle of the speed of body movement (body chest micromotion caused by the heartbeat), and the larger the heartbeat frequency stability index is, the more unstable the heartbeat of a person at the current moment is, the higher the possibility that the current moment of the person belongs to a shallow sleep period is; the smaller the heartbeat frequency stability index is, the more stable the heartbeat of the person at the current moment is, and the higher the possibility that the current moment of the person belongs to the deep sleep period is.
In the embodiment of the invention, the respiratory amplitude variation index and the heartbeat amplitude variation index can also define whether the patient is in the waking period or not besides the function of distinguishing the deep sleep period from the light sleep period, and the amplitude variation reflected in the radar waveform is larger due to the human body movement caused by other factors in the waking period.
In the embodiment of the invention, because the radar has more accurate measurement on the amplitude change, and the existence of the respiration amplitude change index and the heartbeat amplitude change index can better distinguish the waking period from the sleeping period, the sleep staging algorithm of the embodiment combines the respiration amplitude change index, the heartbeat amplitude change index, the respiration frequency stability index and the heartbeat frequency stability index to perform staging estimation, and has higher reliability.
In one embodiment, a method of calculating a respiratory rate stability index may include: calculating the average value of the respiratory frequency of the set time according to the respiratory frequency of each time in the previous set time of the target time, and calculating the respiratory frequency stability index of the target time according to the respiratory frequency of each time in the previous set time of the target time and the average value of the respiratory frequency of the set time.
Specifically, the calculation formula can be expressed as:
Figure 804900DEST_PATH_IMAGE001
wherein the content of the first and second substances,Twhich represents the time of the object and the time of the object,Lit indicates that the set time period is,
Figure 121612DEST_PATH_IMAGE002
an index representing the stability of the breathing frequency at the target moment,
Figure 312422DEST_PATH_IMAGE003
to representtThe breathing frequency at the moment of time is,
Figure 739992DEST_PATH_IMAGE004
represents a mean value of the breathing frequency at the target time; wherein the content of the first and second substances,
Figure 67068DEST_PATH_IMAGE005
in one embodiment, the method of calculating a heart beat frequency stability index may include: calculating the heartbeat frequency average value of the set time length according to the heartbeat frequency of each time in the previous set time length of the target time length, and calculating the heartbeat frequency stability index of the target time length according to the heartbeat frequency of each time in the previous set time length of the target time length and the heartbeat frequency average value of the set time length.
Specifically, the calculation formula can be expressed as:
Figure 321725DEST_PATH_IMAGE006
wherein the content of the first and second substances,
Figure 937515DEST_PATH_IMAGE007
an index representing the stability of the frequency of the heartbeat at the target time,
Figure 231093DEST_PATH_IMAGE008
to representtThe frequency of the heart beats at a moment in time,
Figure 615938DEST_PATH_IMAGE009
representing a mean value of the heartbeat frequency at the target time; wherein the content of the first and second substances,
Figure 71190DEST_PATH_IMAGE010
in one embodiment, a method of calculating a respiratory amplitude variation index may include: calculating the average value of the breathing amplitude of the set time length according to the breathing amplitude of each time in the previous set time length of the target time, and calculating the breathing amplitude change index of the target time according to the breathing amplitude of each time in the previous set time length of the target time and the average value of the breathing amplitude of the set time length.
Specifically, the calculation formula can be expressed as:
Figure 174275DEST_PATH_IMAGE011
wherein the content of the first and second substances,
Figure 5965DEST_PATH_IMAGE012
an index representing the change in the amplitude of breathing at the target moment,
Figure 510895DEST_PATH_IMAGE013
which represents the amplitude of the breathing at the target moment,
Figure 137049DEST_PATH_IMAGE014
indicating the front of the target timeLMean value of breathing amplitude over time.
In one embodiment, the method of calculating a heart beat amplitude variation index may include: calculating the average value of the heartbeat amplitude of the set time according to the heartbeat amplitude of each time in the previous set time of the target time, and calculating the heartbeat amplitude change index of the target time according to the heartbeat amplitude of each time in the previous set time of the target time and the average value of the heartbeat amplitude of the set time.
Specifically, the calculation formula can be expressed as:
Figure 727430DEST_PATH_IMAGE015
wherein the content of the first and second substances,
Figure 566073DEST_PATH_IMAGE016
an index of change in the amplitude of the heartbeat at the target time,
Figure 456669DEST_PATH_IMAGE017
representing the amplitude of the heartbeat at the target time,
Figure 519303DEST_PATH_IMAGE018
indicating the front of the target timeLMean value of the amplitude of the heartbeat over the duration.
In step 103, a sleep stage characteristic value at each moment in a sleep session is calculated based on the sleep stage index.
In the embodiment of the invention, the sleep stage characteristic value of each moment in the sleep period can be calculated according to the sleep stage index of each moment, so that the sleep stage is defined by the sleep stage characteristic value. The sleep stage characteristic value is used for representing the sleep state of the human body in a digital form according to the monitoring result of the radar, the probability that the human body is in a deep sleep stage is higher when the sleep stage characteristic value is smaller, and the probability that the human body is in a light sleep stage or even a waking stage is higher when the sleep stage characteristic value is larger. The sleep staging is carried out based on the sleep staging characteristic value, so that the algorithm provided by the embodiment of the invention is simpler, more direct and faster in processing speed.
The sleep stage index comprises four indexes of a respiratory frequency stability index, a heartbeat frequency stability index, a respiratory amplitude change index and a heartbeat amplitude change index, and the influence weights of the indexes on the sleep stage characteristic value are different, so that normalization processing and comprehensive calculation are required.
In one embodiment, the step 103 may include:
normalizing the sleep staging index;
in this embodiment, the normalization process is intended to simplify the process, and the four indices are normalized into a unified scalar quantity to facilitate the weight assignment and the weight verification.
Optionally, the step of normalizing the sleep stage index may include;
acquiring the maximum value of a respiratory frequency stability index, the maximum value of a heartbeat frequency stability index, the maximum value of a respiratory amplitude change index and the maximum value of a heartbeat amplitude change index in a sleep period;
different persons have large differences with respect to four sleep stage indexes, and even the same person has the same differences in the four sleep stage indexes in different sleep periods. Therefore, the embodiment of the invention adopts the maximum value of the respiratory frequency stability index, the maximum value of the heartbeat frequency stability index, the maximum value of the respiratory amplitude change index and the maximum value of the heartbeat amplitude change index in each sleep period as the benchmarks to carry out normalization processing on the four sleep stage indexes. Namely, taking the ratio of the respiratory frequency stability index at each moment to the maximum value of the respiratory frequency stability index as the normalized respiratory frequency stability index at the moment; taking the ratio of the heartbeat frequency stability index at each moment to the maximum value of the heartbeat frequency stability index as the normalized heartbeat frequency stability index at the moment; taking the ratio of the respiration amplitude change index of each moment to the maximum value of the respiration amplitude change index as the normalized respiration amplitude change index of the moment; and taking the ratio of the heartbeat amplitude change index at each moment to the maximum value of the heartbeat amplitude change index as the normalized heartbeat amplitude change index at the moment.
After the normalization processing, the normalized sleep stage indexes may be subjected to weighted summation based on preset weights of the sleep stage indexes, so as to obtain a sleep stage characteristic value at each moment in the sleep period.
It should be noted that, because the influence of respiration on body movement is larger than the influence of heartbeat on body movement, the respiratory frequency stability index and the respiratory amplitude variation index calculated based on the radar information are more accurate than the heartbeat frequency stability index and the heartbeat amplitude variation index, and therefore, in practical application, the weights of the respiratory frequency stability index and the respiratory amplitude variation index are both larger than the heartbeat frequency stability index and the heartbeat amplitude variation index. Of course, the final weight of the four indexes can be finally determined by testing and verifying a large amount of data, and a most accurate weight distribution mode is finally determined.
In the embodiment of the present invention, the four sleep stage indexes have different influence weights on the sleep stage characteristic value, and the four normalized sleep stage indexes may be subjected to weighted summation to finally obtain the sleep stage characteristic value at each time in the sleep period.
In step 104, the sleep stage characteristic value is compared with a preset stage threshold, and a sleep stage corresponding to each time in the sleep period is determined according to the comparison result.
In the embodiment of the present invention, the sleep stage characteristic value corresponding to each time may be compared with a preset stage threshold, and finally, the sleep stage corresponding to each time in the sleep period is determined.
In one embodiment, the staging threshold includes an awake threshold and a light sleep threshold, and the step 103 may specifically include:
if the sleep stage characteristic value of the target time is greater than or equal to a preset waking threshold, the sleep stage corresponding to the time is a waking stage;
if the sleep stage characteristic value of the target time is smaller than the waking threshold and is greater than or equal to a preset light sleep threshold, the sleep stage corresponding to the time is a light sleep stage;
and if the sleep stage characteristic value of the target moment is smaller than the light sleep threshold, the sleep stage corresponding to the moment is a deep sleep stage.
In the embodiment of the present invention, in order to correctly distinguish the waking period, the light sleep period and the deep sleep period, two defined thresholds, namely, the waking threshold and the light sleep threshold, may be set; in practical application, a wakefulness threshold is taken as a boundary, a sleep stage characteristic value (representing the digital characteristic of a human body sleep state) is greater than or equal to the wakefulness threshold, which indicates that a human body currently has larger action frequency and action amplitude and belongs to a wakefulness stage; the sleep stage characteristic value (the digitalized characteristic representing the sleep state of the human body) is smaller than the wakefulness threshold value, which indicates that the human body currently has smaller action frequency and action amplitude and currently belongs to the sleep stage (including a light sleep stage and a deep sleep stage). Taking a shallow sleep threshold as a boundary, wherein the sleep stage characteristic value is smaller than the waking threshold but larger than or equal to the shallow sleep threshold, and indicates that the human body currently has smaller action frequency and action amplitude and currently belongs to the shallow sleep stage; the characteristic value of the sleep stage is smaller than the light sleep threshold value, which indicates that the human body has extremely small action frequency and action amplitude at present and belongs to the deep sleep stage at present.
In the embodiment of the invention, the wakefulness threshold and the light sleep threshold can be obtained according to a large amount of empirical data statistics. By carrying out division confirmation on the waking period, the light sleep period and the deep sleep period at each moment, the waking period duration, the light sleep period duration and the deep sleep period duration corresponding to the whole sleep period can be accurately and reliably obtained finally.
In the embodiment of the invention, the waking period duration, the light sleep period duration and the deep sleep period duration can be quickly and reliably obtained by using the sleep staging method, and the waking period duration, the light sleep period duration and the deep sleep period duration can be used for generating a sleep quality analysis report so as to reflect the health condition of a human body from the sleep perspective and provide corresponding solutions.
In the sleep quality analysis report, in addition to the sleep quality information such as the waking period duration, the light sleep period duration, and the deep sleep period duration, a respiratory frequency variation curve, a heartbeat frequency variation curve, a respiratory frequency distribution histogram, a heartbeat frequency distribution histogram, the getting-up time, and the getting-up time may be generated based on radar information, so as to more comprehensively reflect the sleep condition of the human body.
As can be seen from the above, in the embodiment of the present invention, the respiration heartbeat information in the sleep period monitored by the vital sign monitoring radar is obtained, the sleep stage index of each time in the sleep period is calculated based on the respiration heartbeat information, the sleep stage characteristic value of each time in the sleep period is calculated based on the sleep stage index, the sleep stage characteristic value is compared with the preset stage threshold, and the sleep stage corresponding to each time in the sleep period is determined according to the comparison result. On one hand, the monitoring information for sleep staging is acquired based on radar, and belongs to non-contact monitoring, so that a user has no perception and normal sleep of the user is not influenced; on the other hand, compared with the method that the sleep state of a human body can be more comprehensively reflected by relying on heart rate information based on the respiratory heartbeat information, the reliability of sleep staging is higher, in addition, the sleep staging index of each moment can be calculated, the sleep staging characteristic value of the moment can be determined, the sleep stage can be measured in seconds (each moment), and the sleep stage can be distinguished, so that the processing process is simple, convenient and accurate.
It should be understood that, the sequence numbers of the steps in the foregoing embodiments do not imply an execution sequence, and the execution sequence of each process should be determined by its function and inherent logic, and should not constitute any limitation to the implementation process of the embodiments of the present invention.
The following are embodiments of the apparatus of the invention, reference being made to the corresponding method embodiments described above for details which are not described in detail therein.
Fig. 2 is a schematic structural diagram of a sleep staging apparatus based on radar information according to an embodiment of the present invention, which only shows the relevant parts related to the embodiment of the present invention for convenience of illustration, and the details are as follows:
as shown in fig. 2, the sleep staging apparatus 2 based on radar information may include: a first acquisition unit 21, a staging index calculation unit 22, a staging feature value calculation unit 23, and a sleep session determination unit 24.
The first acquisition unit 21 is configured to acquire respiratory heartbeat information within a sleep period monitored by the vital sign monitoring radar;
a staging index calculation unit 22, configured to calculate a sleep staging index at each time in the sleep session based on the respiratory heartbeat information acquired by the first acquisition unit 21;
a staging feature value calculation unit 23 for calculating a sleep staging feature value at each time in the sleep session based on the sleep staging index calculated by the staging index calculation unit 22;
the sleep period determining unit 24 is configured to compare the sleep stage feature value calculated by the stage feature value calculating unit 23 with a preset stage threshold, and determine a sleep period corresponding to each moment in the sleep period according to a comparison result.
In a possible implementation manner, the respiratory heartbeat information includes: the respiratory frequency, the heartbeat frequency, the respiratory amplitude and the heartbeat amplitude of each moment in the sleep period; the sleep staging index includes: a respiratory frequency stability index, a heartbeat frequency stability index, a respiratory amplitude variation index and a heartbeat amplitude variation index;
correspondingly, the staging index calculating unit 22 is specifically configured to calculate a respiratory frequency average value of the set time according to the respiratory frequency of each time in the previous set time of the target time, and calculate a respiratory frequency stability index of the target time according to the respiratory frequency of each time in the previous set time of the target time and the respiratory frequency average value of the set time;
calculating the heartbeat frequency average value of the set time length according to the heartbeat frequency of each time in the previous set time length of the target time length, and calculating the heartbeat frequency stability index of the target time length according to the heartbeat frequency of each time in the previous set time length of the target time length and the heartbeat frequency average value of the set time length;
calculating the average value of the breathing amplitude of the set time length according to the breathing amplitude of each time in the previous set time length of the target time, and calculating the heartbeat amplitude change index of the target time according to the breathing amplitude of each time in the previous set time length of the target time and the average value of the breathing amplitude of the set time length;
calculating the average value of the heartbeat amplitude of the set time according to the heartbeat amplitude of each time in the previous set time of the target time, and calculating the heartbeat amplitude change index of the target time according to the heartbeat amplitude of each time in the previous set time of the target time and the average value of the heartbeat amplitude of the set time.
In a possible implementation manner, in particular, the method for calculating the respiratory rate stability index at the target time by the staging index calculation unit 22 includes:
Figure 596980DEST_PATH_IMAGE001
wherein the content of the first and second substances,Twhich represents the time of the object and the time of the object,Lit indicates that the set time period is,
Figure 472270DEST_PATH_IMAGE002
an index representing the stability of the breathing frequency at the target moment,
Figure 748531DEST_PATH_IMAGE003
to representtThe breathing frequency of the moment in time,
Figure 919749DEST_PATH_IMAGE004
represents the mean value of the breathing frequency at the target moment; wherein the content of the first and second substances,
Figure 281460DEST_PATH_IMAGE005
the method for calculating the heart beat frequency stability index at the target moment comprises the following steps:
Figure 461906DEST_PATH_IMAGE006
wherein the content of the first and second substances,
Figure 530356DEST_PATH_IMAGE007
an index representing the stability of the frequency of the heartbeat at the target time,
Figure 934792DEST_PATH_IMAGE008
to representtThe frequency of the heart beats at a moment in time,
Figure 987062DEST_PATH_IMAGE009
representing a mean value of the heartbeat frequency at the target time; wherein the content of the first and second substances,
Figure 440040DEST_PATH_IMAGE010
the method for calculating the respiratory amplitude variation index at the target moment comprises the following steps:
Figure 690893DEST_PATH_IMAGE011
wherein, the first and the second end of the pipe are connected with each other,
Figure 469493DEST_PATH_IMAGE012
an index representing the change in the amplitude of breathing at the target moment,
Figure 743480DEST_PATH_IMAGE013
which represents the amplitude of the breathing at the target moment,
Figure 62465DEST_PATH_IMAGE014
indicating the front of the target timeLMean value of breathing amplitude over a time period;
the method for calculating the heart beat amplitude change index at the target moment comprises the following steps:
Figure 606973DEST_PATH_IMAGE015
wherein the content of the first and second substances,
Figure 353212DEST_PATH_IMAGE016
an index of change in the amplitude of the heartbeat at the target time,
Figure 848915DEST_PATH_IMAGE017
the amplitude of the heartbeat at the target time is indicated,
Figure 971592DEST_PATH_IMAGE018
indicating the front of the target timeLMean value of the amplitude of the heartbeat over the duration.
In a possible implementation manner, the sleep staging apparatus 2 further includes:
a normalization processing unit, configured to perform normalization processing on the sleep staging index;
correspondingly, the stage characteristic value calculating unit 23 is specifically configured to perform weighted summation on the normalized sleep stage indexes based on preset weights of the sleep stage indexes to obtain a sleep stage characteristic value at each moment in a sleep period;
the weights of the respiratory frequency stability index and the respiratory amplitude change index are both larger than the heartbeat frequency stability index and the heartbeat amplitude change index.
In a possible implementation manner, the sleep staging apparatus 2 further includes:
the second acquisition unit is used for acquiring the maximum value of the respiratory frequency stability index, the maximum value of the heartbeat frequency stability index, the maximum value of the respiratory amplitude change index and the maximum value of the heartbeat amplitude change index in a sleep period;
correspondingly, the staging characteristic value calculating unit 23 is specifically configured to use a ratio of the respiratory frequency stability index at each time to a maximum value of the respiratory frequency stability index acquired by the second acquiring unit as the normalized respiratory frequency stability index at the time;
taking the ratio of the heartbeat frequency stability index at each moment to the maximum value of the heartbeat frequency stability index acquired by the second acquisition unit as the normalized heartbeat frequency stability index at the moment;
the ratio of the respiration amplitude change index at each moment to the maximum value of the respiration amplitude change index acquired by the second acquisition unit is used as the normalized respiration amplitude change index at the moment;
and taking the ratio of the heartbeat amplitude change index at each moment to the maximum value of the heartbeat amplitude change index acquired by the second acquisition unit as the normalized heartbeat amplitude change index at the moment.
In a possible implementation manner, the staging threshold includes an awake threshold and a light sleep threshold, and the sleep stage determining unit 24 is specifically configured to, if the sleep stage feature value of the target time is greater than or equal to a preset awake threshold, determine that the sleep stage corresponding to the target time is an awake stage;
if the sleep stage characteristic value of the target time is smaller than the waking threshold and is greater than or equal to a preset light sleep threshold, the sleep stage corresponding to the time is a light sleep stage;
and if the sleep stage characteristic value of the target time is smaller than the light sleep threshold, the sleep stage corresponding to the time is a deep sleep stage.
In a possible implementation manner, the sleep staging apparatus 2 further includes:
a third obtaining unit, configured to obtain in-bed information and rising information monitored by the vital sign monitoring radar before obtaining the respiratory heartbeat information during the sleep period monitored by the vital sign monitoring radar;
the starting time determining unit is used for determining the getting-on time and the getting-up time according to the in-bed information and the getting-up information acquired by the third acquiring unit;
and the sleep period determining unit is used for determining the sleep period according to the getting-up time and the getting-up time determined by the starting time determining unit.
As can be seen from the above, in the embodiment of the present invention, the respiration heartbeat information in the sleep period monitored by the vital sign monitoring radar is obtained, the sleep stage index of each time in the sleep period is calculated based on the respiration heartbeat information, the sleep stage characteristic value of each time in the sleep period is calculated based on the sleep stage index, the sleep stage characteristic value is compared with the preset stage threshold, and the sleep stage corresponding to each time in the sleep period is determined according to the comparison result. On one hand, the monitoring information for sleep staging is acquired based on radar, and belongs to non-contact monitoring, so that a user has no perception and normal sleep of the user is not influenced; on the other hand, compared with the method that the sleep state of a human body can be more comprehensively reflected by relying on heart rate information based on the respiratory heartbeat information, the reliability of sleep staging is higher, in addition, the sleep staging index of each moment can be calculated, the sleep staging characteristic value of the moment can be determined, the sleep stage can be measured in seconds (each moment), and the sleep stage can be distinguished, so that the processing process is simple, convenient and accurate.
Fig. 3 is a schematic diagram of a terminal according to an embodiment of the present invention. As shown in fig. 3, the terminal 3 of this embodiment includes: a processor 30, a memory 31 and a computer program 32 stored in said memory 31 and executable on said processor 30. The processor 30, when executing the computer program 32, implements the steps in the various radar-information-based sleep staging method embodiments described above, such as steps 101-104 shown in fig. 1. Alternatively, the processor 30, when executing the computer program 32, implements the functions of the modules/units in the above-mentioned device embodiments, such as the functions of the units 21 to 24 shown in fig. 2.
Illustratively, the computer program 32 may be divided into one or more units, which are stored in the memory 31 and executed by the processor 30 to accomplish the present invention. The one or more units may be a series of computer program instruction segments capable of performing specific functions, which are used to describe the execution of the computer program 32 in the terminal 3. For example, the computer program 32 may be divided into the units 21 to 24 shown in fig. 2.
The terminal 3 may be a desktop computer, a notebook, a palm computer, a cloud server, or other computing devices. The terminal 3 may include, but is not limited to, a processor 30, a memory 31. It will be appreciated by those skilled in the art that fig. 3 is only an example of a terminal 3 and does not constitute a limitation of the terminal 3 and may comprise more or less components than those shown, or some components may be combined, or different components, e.g. the terminal may further comprise input output devices, network access devices, buses, etc.
The Processor 30 may be a Central Processing Unit (CPU), other general purpose Processor, a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), a Field Programmable Gate Array (FPGA) or other Programmable logic device, discrete Gate or transistor logic device, discrete hardware component, or the like. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like.
The memory 31 may be an internal storage unit of the terminal 3, such as a hard disk or a memory of the terminal 3. The memory 31 may also be an external storage device of the terminal 3, such as a plug-in hard disk, a Smart Media Card (SMC), a Secure Digital (SD) Card, a Flash memory Card (Flash Card) and the like provided on the terminal 3. Further, the memory 31 may also comprise both an internal memory unit and an external memory device of said terminal 3. The memory 31 is used for storing the computer program and other programs and data required by the terminal. The memory 31 may also be used to temporarily store data that has been output or is to be output.
It will be apparent to those skilled in the art that, for convenience and brevity of description, only the above-mentioned division of the functional units and modules is illustrated, and in practical applications, the above-mentioned function distribution may be performed by different functional units and modules according to needs, that is, the internal structure of the apparatus is divided into different functional units or modules to perform all or part of the above-mentioned functions. Each functional unit and module in the embodiments may be integrated in one processing unit, or each unit may exist alone physically, or two or more units are integrated in one unit, and the integrated unit may be implemented in a form of hardware, or in a form of software functional unit. In addition, specific names of the functional units and modules are only used for distinguishing one functional unit from another, and are not used for limiting the protection scope of the present application. The specific working processes of the units and modules in the system may refer to the corresponding processes in the foregoing method embodiments, and are not described herein again.
In the above embodiments, the descriptions of the respective embodiments have respective emphasis, and reference may be made to the related descriptions of other embodiments for parts that are not described or illustrated in a certain embodiment.
Those of ordinary skill in the art will appreciate that the various illustrative elements and algorithm steps described in connection with the embodiments disclosed herein may be implemented as electronic hardware or combinations of computer software and electronic hardware. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the implementation. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present invention.
In the embodiments provided in the present invention, it should be understood that the disclosed apparatus/terminal and method may be implemented in other ways. For example, the above-described apparatus/terminal embodiments are merely illustrative, and for example, the division of the modules or units is only one logical division, and there may be other divisions when actually implemented, for example, a plurality of units or components may be combined or may be integrated into another system, or some features may be omitted, or not executed. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be an indirect coupling or communication connection through some interfaces, devices or units, and may be in an electrical, mechanical or other form.
The units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of the embodiment.
In addition, functional units in the embodiments of the present invention may be integrated into one processing unit, or each unit may exist alone physically, or two or more units are integrated into one unit. The integrated unit can be realized in a form of hardware, and can also be realized in a form of a software functional unit.
The integrated modules/units, if implemented in the form of software functional units and sold or used as separate products, may be stored in a computer readable storage medium. Based on such understanding, all or part of the flow of the method according to the above embodiments may be implemented by a computer program, which may be stored in a computer readable storage medium and used by a processor to implement the steps of the embodiments of the sleep staging method based on radar information. Wherein the computer program comprises computer program code, which may be in the form of source code, object code, an executable file or some intermediate form, etc. The computer-readable medium may include: any entity or device capable of carrying the computer program code, recording medium, usb disk, removable hard disk, magnetic disk, optical disk, computer Memory, Read-Only Memory (ROM), Random Access Memory (RAM), electrical carrier wave signals, telecommunications signals, software distribution medium, and the like. It should be noted that the computer readable medium may contain other components which may be suitably increased or decreased as required by legislation and patent practice in jurisdictions, for example, in some jurisdictions, computer readable media which may not include electrical carrier signals and telecommunications signals in accordance with legislation and patent practice.
The above-mentioned embodiments are only used for illustrating the technical solutions of the present invention, and not for limiting the same; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; such modifications and substitutions do not depart from the spirit and scope of the embodiments of the present invention, and they should be construed as being included therein.

Claims (10)

1. A sleep staging method based on radar information is applied to a health monitoring system, the health monitoring system comprises a vital sign monitoring radar, and the method comprises the following steps:
acquiring respiratory heartbeat information in a sleep period monitored by a vital sign monitoring radar;
calculating a sleep staging index at each moment in a sleep period based on the breathing heartbeat information;
calculating a sleep stage characteristic value of each moment in a sleep period based on the sleep stage index;
comparing the sleep stage characteristic value with a preset stage threshold value, and determining a sleep stage corresponding to each moment in a sleep period according to a comparison result;
the respiratory heartbeat information includes: the respiratory amplitude and the heartbeat amplitude at each moment in the sleep period; the sleep staging index includes: a respiratory amplitude variation index and a heartbeat amplitude variation index;
correspondingly, the calculating the sleep staging index at each moment in the sleep period based on the respiratory heartbeat information includes:
calculating the average value of the breathing amplitude of the set time length according to the breathing amplitude of each time in the previous set time length of the target time, and calculating the breathing amplitude change index of the target time according to the breathing amplitude of each time in the previous set time length of the target time and the average value of the breathing amplitude of the set time length;
calculating the heartbeat amplitude mean value of the set time length according to the heartbeat amplitude of each time in the previous set time length of the target time length, and calculating the heartbeat amplitude change index of the target time length according to the heartbeat amplitude of each time in the previous set time length of the target time length and the heartbeat amplitude mean value of the set time length;
the method for calculating the respiratory amplitude variation index at the target moment comprises the following steps:
Figure FDA0003634266130000021
wherein, breath amplitude variation index at the target time is represented by breath amplitude variation (T), breath amplitude (T) at the target time is represented by breath amplitude at the target time,
Figure FDA0003634266130000022
represents the mean value of the breathing amplitude within the preceding L duration of the target moment;
the method for calculating the heart beat amplitude change index at the target moment comprises the following steps:
Figure FDA0003634266130000023
wherein, heartbeat amplitude change index (T) represents heartbeat amplitude change index at the target moment, heartbeat amplitude (T) represents heartbeat amplitude at the target moment,
Figure FDA0003634266130000024
representing the average value of the heart beat amplitudes within the preceding L time period of the target time.
2. The radar information-based sleep staging method of claim 1, wherein the respiratory heartbeat information further includes: respiratory and heartbeat frequencies at each time during a sleep session;
the sleep staging index further comprises: a respiratory frequency stability index and a heartbeat frequency stability index;
correspondingly, the calculating the sleep stage index of each moment in the sleep period based on the breathing heartbeat information further comprises:
calculating the average value of the respiratory frequency of the set time length according to the respiratory frequency of each time in the previous set time length of the target time length, and calculating the respiratory frequency stability index of the target time length according to the respiratory frequency of each time in the previous set time length of the target time length and the average value of the respiratory frequency of the set time length;
calculating the heartbeat frequency average value of the set time length according to the heartbeat frequency of each time in the previous set time length of the target time length, and calculating the heartbeat frequency stability index of the target time length according to the heartbeat frequency of each time in the previous set time length of the target time length and the heartbeat frequency average value of the set time length.
3. The radar-information-based sleep staging method according to claim 2, wherein the method of calculating the respiratory rate stability index at the target time includes:
Figure FDA0003634266130000031
wherein T represents a target time, L represents a set time period, breathfreqstability (T) represents a respiratory frequency stability index at the target time, breathfreq (T) represents a respiratory frequency at time T, and breathfreqmean (T) represents a respiratory frequency mean value for a set time period immediately preceding the target time; wherein the content of the first and second substances,
Figure FDA0003634266130000032
the method for calculating the heart beat frequency stability index at the target moment comprises the following steps:
Figure FDA0003634266130000033
wherein, heartFreqStability (T) represents the heartbeat frequency stability index at the target time, heartFreq (T) represents the heartbeat frequency at the t time, and heartFreqmean (T) represents the heartbeat frequency mean value of the previous set time length at the target time; wherein the content of the first and second substances,
Figure FDA0003634266130000034
4. the radar-information-based sleep staging method of claim 2, wherein the calculating a sleep staging characteristic value for each time instant within a sleep session based on the sleep staging index comprises:
normalizing the sleep staging index;
carrying out weighted summation on the normalized sleep stage indexes based on the preset weight of each sleep stage index to obtain a sleep stage characteristic value of each moment in a sleep period;
the weights of the respiratory frequency stability index and the respiratory amplitude change index are both larger than the heartbeat frequency stability index and the heartbeat amplitude change index.
5. The radar-information-based sleep staging method of claim 4, wherein the normalizing the sleep staging index comprises;
acquiring a maximum value of a respiratory frequency stability index, a maximum value of a heartbeat frequency stability index, a maximum value of a respiratory amplitude change index and a maximum value of a heartbeat amplitude change index in a sleep period;
taking the ratio of the respiratory frequency stability index at each moment to the maximum value of the respiratory frequency stability index as the normalized respiratory frequency stability index at the moment;
taking the ratio of the heartbeat frequency stability index at each moment to the maximum value of the heartbeat frequency stability index as the normalized heartbeat frequency stability index at the moment;
taking the ratio of the respiration amplitude change index of each moment to the maximum value of the respiration amplitude change index as the normalized respiration amplitude change index of the moment;
and taking the ratio of the heartbeat amplitude change index at each moment to the maximum value of the heartbeat amplitude change index as the normalized heartbeat amplitude change index at the moment.
6. A sleep staging method based on radar information according to any one of claims 1 to 5, wherein the staging threshold comprises an awake threshold and a shallow sleep threshold;
the step of comparing the sleep stage characteristic value with a preset stage threshold value and determining the sleep stage corresponding to each moment in the sleep period according to the comparison result comprises:
if the sleep stage characteristic value of the target time is greater than or equal to the waking threshold, the sleep stage corresponding to the time is the waking stage;
if the sleep stage characteristic value of the target moment is smaller than the waking threshold and is greater than or equal to a preset light sleep threshold, the sleep stage corresponding to the moment is a light sleep stage;
and if the sleep stage characteristic value of the target moment is smaller than the light sleep threshold, the sleep stage corresponding to the moment is a deep sleep stage.
7. The radar-information-based sleep staging method of claim 6, further comprising, prior to the obtaining respiratory heartbeat information for a sleep session monitored by a vital signs monitoring radar:
acquiring in-bed information and rising information monitored by a vital sign monitoring radar;
determining the getting-on time and the getting-up time according to the in-bed information and the getting-up information;
and determining a sleep period according to the time of getting up and the time of getting up.
8. A sleep staging device based on radar information, applied to a health monitoring system including a vital sign monitoring radar, the sleep staging device comprising:
the first acquisition unit is used for acquiring the respiration heartbeat information in the sleep period monitored by the vital sign monitoring radar;
the staging index calculating unit is used for calculating the sleep staging index of each moment in the sleep period based on the respiratory heartbeat information;
the staging characteristic value calculating unit is used for calculating a sleep staging characteristic value of each moment in a sleep period based on the sleep staging index;
the sleep stage determining unit is used for comparing the sleep stage characteristic value with a preset stage threshold value and determining a sleep stage corresponding to each moment in a sleep period according to a comparison result;
the respiration heartbeat information comprises respiration amplitude and heartbeat amplitude; the sleep staging index comprises a respiration amplitude variation index and a heartbeat amplitude variation index;
correspondingly, the staging index calculation unit is specifically configured to calculate a respiratory amplitude mean value of the set time length according to the respiratory amplitude of each time in the previous set time length of the target time length, and calculate a respiratory amplitude variation index of the target time length according to the respiratory amplitude of each time in the previous set time length of the target time length and the respiratory amplitude mean value of the set time length; calculating the heartbeat amplitude mean value of the set time length according to the heartbeat amplitude of each time in the previous set time length of the target time length, and calculating the heartbeat amplitude change index of the target time length according to the heartbeat amplitude of each time in the previous set time length of the target time length and the heartbeat amplitude mean value of the set time length;
the method for calculating the breathing amplitude change index of the target moment by the staging index calculation unit comprises the following steps:
Figure FDA0003634266130000061
wherein, breath amplitude variation index at the target time is represented by breath amplitude variation (T), breath amplitude (T) at the target time is represented by breath amplitude at the target time,
Figure FDA0003634266130000062
represents a mean value of the breathing amplitudes within the preceding L duration of the target time;
the method for calculating the heartbeat amplitude change index at the target moment by the staging index calculation unit comprises the following steps:
Figure FDA0003634266130000063
wherein, heartbeat amplitude change index (T) represents heartbeat amplitude change index at the target moment, heartbeat amplitude (T) represents heartbeat amplitude at the target moment,
Figure FDA0003634266130000071
representing the average value of the heart beat amplitudes within the preceding L time period of the target time.
9. A terminal comprising a memory, a processor and a computer program stored in the memory and executable on the processor, characterized in that the processor when executing the computer program implements the steps of the radar information based sleep staging method according to any of the preceding claims 1 to 7.
10. A computer-readable storage medium, in which a computer program is stored, which, when being executed by a processor, carries out the steps of the radar-information based sleep staging method according to any one of the preceding claims 1 to 7.
CN202111498745.1A 2021-12-09 2021-12-09 Sleep staging method and device based on radar information and terminal Active CN114176521B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111498745.1A CN114176521B (en) 2021-12-09 2021-12-09 Sleep staging method and device based on radar information and terminal

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111498745.1A CN114176521B (en) 2021-12-09 2021-12-09 Sleep staging method and device based on radar information and terminal

Publications (2)

Publication Number Publication Date
CN114176521A CN114176521A (en) 2022-03-15
CN114176521B true CN114176521B (en) 2022-07-01

Family

ID=80542900

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111498745.1A Active CN114176521B (en) 2021-12-09 2021-12-09 Sleep staging method and device based on radar information and terminal

Country Status (1)

Country Link
CN (1) CN114176521B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115778352B (en) * 2023-02-02 2023-05-09 北京华声医学科技开发公司 Sleep quality assessment method, equipment, system and medium based on millimeter wave radar
CN115868941B (en) * 2023-03-03 2023-05-05 深圳市魔样科技有限公司 Information management method for intelligent ring
CN115968080A (en) * 2023-03-17 2023-04-14 森思泰克河北科技有限公司 Night lamp control method and device based on radar, radar and storage medium

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108416367A (en) * 2018-02-08 2018-08-17 南京理工大学 Sleep stage method based on multi-sensor data decision level fusion
CN109480787A (en) * 2018-12-29 2019-03-19 中国科学院合肥物质科学研究院 A kind of contactless sleep monitor equipment and sleep stage method based on ULTRA-WIDEBAND RADAR
CN113180596A (en) * 2021-04-07 2021-07-30 中山大学 Non-contact sleep analysis method and device and storage medium

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112205964A (en) * 2019-07-11 2021-01-12 京东方科技集团股份有限公司 Sleep intervention device and sleep intervention management system

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108416367A (en) * 2018-02-08 2018-08-17 南京理工大学 Sleep stage method based on multi-sensor data decision level fusion
CN109480787A (en) * 2018-12-29 2019-03-19 中国科学院合肥物质科学研究院 A kind of contactless sleep monitor equipment and sleep stage method based on ULTRA-WIDEBAND RADAR
CN113180596A (en) * 2021-04-07 2021-07-30 中山大学 Non-contact sleep analysis method and device and storage medium

Also Published As

Publication number Publication date
CN114176521A (en) 2022-03-15

Similar Documents

Publication Publication Date Title
CN114176521B (en) Sleep staging method and device based on radar information and terminal
CN108388912B (en) Sleep staging method based on multi-sensor feature optimization algorithm
RU2689177C1 (en) Sleep control device
US8155733B2 (en) Lie detection method and system
CN109276241B (en) Pressure identification method and equipment
US20110257556A1 (en) Systems and methods for quantitatively characterizing slow wave activities and states in sleep
WO2016029166A1 (en) Automated diagnosis based at least in part on pulse waveforms
Phillips et al. WristO2: Reliable peripheral oxygen saturation readings from wrist-worn pulse oximeters
Chen et al. Algorithm for heart rate extraction in a novel wearable acoustic sensor
US11617545B2 (en) Methods and systems for adaptable presentation of sensor data
Wang et al. PulseDB: A large, cleaned dataset based on MIMIC-III and VitalDB for benchmarking cuff-less blood pressure estimation methods
JP2008173160A (en) Method of analyzing fluctuation of heart rate and method of determining health state using the same
CN113520343A (en) Sleep risk prediction method and device and terminal equipment
Mamun et al. Swallowing accelerometry signal feature variations with sensor displacement
CN107811610B (en) Respiration rate detection method and device, electronic equipment and storage medium
US20200330040A1 (en) Method and system for detecting noise in vital sign signal
WO2018104970A1 (en) Pulse detection, measurement and analysis based health management system, method and apparatus
Ma et al. A new measure to characterize multifractality of sleep electroencephalogram
CN106504235B (en) Rhythm of the heart method based on image procossing
CN112598033B (en) Physiological signal processing method, device, equipment and storage medium
CN107280650B (en) Method and device for acquiring characteristic parameters of living body
Lu et al. Pulse waveform analysis for pregnancy diagnosis based on machine learning
Pinheiro et al. Multi-sensor wearable health device framework for real-time monitoring of elderly patients using a mobile application and high-resolution parameter estimation
CN115023180A (en) Method and system for determining cardiovascular parameters
Mouney et al. A new blood pressure estimation approach using ppg sensors: Subject specific evaluation over a long-term period

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