发明内容
(一)要解决的技术问题
鉴于现有技术的上述缺点、不足,本发明提供一种基于毫米波雷达的睡眠质量评估方法、设备、系统及介质,其解决了现有方案难以针对用户睡眠质量作出客观、长周期持续监测以及兼顾舒适度的技术问题。
(二)技术方案
为了达到上述目的,本发明采用的主要技术方案包括:
第一方面,本发明实施例提供一种基于毫米波雷达的睡眠质量评估方法,包括:
通过毫米波雷达监测设备采集被测者的生命体征数据;
依据获取的所述生命体征数据,通过构建生命体征数据曲线执行被测者的体动频次分析、呼吸混乱程度分析以及在床时长分析的一种或多种操作;
基于得到的所述体动频次、所述呼吸混乱程度以及所述在床时长之中的一种或多种信息,进行被测者的睡眠质量评估。
可选地,所述生命体征数据包括:体动功率、呼吸波形以及在/离床标志,且所述体动功率、所述呼吸波形以及所述在/离床标志在毫米波雷达监测设备侧的更新频次均为1秒钟;
其中,
所述体动频率为毫米波雷达监测设备输出的一个变量,指的是感兴趣距离范围内雷达回波功率之和,感兴趣距离范围指的是雷达视线范围内人体所处的距离区间;
所述呼吸波形为毫米波雷达监测设备输出的一个变量,指的是毫米波雷达设备提取的反映人体呼吸运动的信号波形,该呼吸波形代表了人体胸腔相对于毫米波雷达监测设备的距离变化情况;
所述在/离床标志为毫米波雷达监测设备输出的一个变量,指的是毫米波雷达设备对监测对象是否在床上做出的判断,该值为0时表示监测对象在床上,该值为1时表示监测对象不在床上。
可选地,依据获取的所述生命体征数据,通过构建生命体征数据曲线执行被测者的体动频次分析、呼吸混乱程度分析以及在床时长分析的一种或多种操作包括:
执行被测者的体动频次分析操作包括:
对获取的体动功率数据进行处理,得到T 1分钟内体动频次曲线;
采用阈值法对所述T 1分钟内体动频次曲线上的各个取值点进行二值划分,将大于阈值N 1的点赋值为1,否则赋值为0,得到体动频次1值初步区段与体动频次0值初步区段;
将所述体动频次0值初步区段中持续时间小于T 2分钟的点重新赋值为1,得到体动频次重新赋值区段与体动频次0值最终区段,并将体动频次0值最终区段视为体动稀少区段;
基于所述体动频次1值初步区段与所述体动频次重新赋值区段,得体动频繁区段;
执行被测者的呼吸混乱程度分析操作包括:
对获取的所述呼吸波形数据进行处理,得到T 3分钟内呼吸混乱程度曲线;
采用阈值法对所述T 3分钟内呼吸混乱程度曲线上的各个取值点进行二值划分,将低于阈值τ 1的点赋值为1,否则赋值为0,得到呼吸混乱程度1值初步区段与呼吸混乱程度0值初步区段;
将所述呼吸混乱程度1值初步区段中持续时间小于T 4分钟的点重新赋值为0,得到呼吸混乱程度重新赋值区段与呼吸混乱程度1值最终区段,并将呼吸混乱程度1值最终区段视为呼吸平稳区段;
基于所述呼吸混乱程度0值初步区段与所述呼吸混乱程度重新赋值区段,得到呼吸非平稳区段;
执行被测者的在床时长分析操作包括:依据获取的在/离床标志数据,得到在床时长;
其中,T 1=2~10min、T 2=2~10min、T 3=2~5min、T 4=2~5min、N 1=3~15次、τ 1=0.1~0.3。
可选地,对获取的体动功率数据进行处理,得到T 1分钟内体动频次曲线包括:
依据经10Hz采样得到的所述体动功率对被测者的体动事件进行滑窗检测,采用有序统计CFAR检测器,并设置参考单元对应时间长度为36秒、保护单元对应时间长度为20秒以及滑窗步进为1秒,经处理得到体动事件检测二值结果;
基于所述体动事件检测二值结果,采用T 1分钟窗长和1秒钟滑窗步进计算T 1分钟内体动次数,得到T 1分钟内体动频次曲线。
可选地,对经10Hz采样的所述呼吸波形数据进行处理,得到T 3分钟内呼吸混乱程度曲线包括:
采用20秒窗长和1秒步进对经10Hz采样得到的呼吸波形进行处理,得到细粒度的呼吸混乱程度曲线;
在每个处理窗口内,对当前20秒时长的呼吸波形数据进行快速傅里叶变换,得到当前20秒时长呼吸波形数据的傅里叶谱;
在傅里叶谱上0.1Hz~0.7Hz频率范围内寻找功率最强的谱峰点,根据下式计算谱峰功率占比:
式中,r为谱峰功率占比,p(f)为傅里叶谱,f为频率,f peak为最强的谱峰点所对应的频率点;
依据所述谱峰功率占比通过下式得到呼吸混乱程度值:
ChaosIndex=2/(1+exp(r/(1-r)))
采用T 3分钟窗长和1秒步进对所述细粒度的呼吸混乱程度曲线进行平滑处理,得到T 3分钟内呼吸混乱程度曲线。
可选地,基于得到的所述体动频次、所述呼吸混乱程度以及所述在床时长之中的一种或多种信息,进行被测者的睡眠质量评估包括:
依据得到的所述体动稀少区段得到体动稀少时长;
依据得到的所述呼吸平稳区段得到呼吸平稳时长;
依据监测时间段内用户首次进入在床状态到其进入睡眠状态的时间点的时长得到入睡时长;其中,所述进入在床状态由在/离床标志确定,所述进入睡眠状态的时间点指的是被测者在床状态内首次体动稀少时长超过10分钟的起始时刻;
基于所述在床时长、所述入睡时长、所述体动稀少时长以及所述呼吸平稳时长,进行被测者的睡眠质量评估。
可选地,被测者的睡眠质量评估是通过如下睡眠质量评估模型实现:
其中,SleepQuality表示被测者的睡眠质量评估值,T Bed为在床时长、T Sopite 为入睡时长、T Quiet 为体动稀少时长、T Steady 为呼吸平稳时长。
第二方面,本发明实施例提供一种基于毫米波雷达的睡眠质量评估设备,包括:
数据获取模块,用于通过毫米波雷达监测设备采集被测者的生命体征数据;
数据处理模块,用于依据获取的所述生命体征数据,通过构建生命体征数据曲线执行被测者的体动频次分析、呼吸混乱程度分析以及在床时长分析的一种或多种操作;
睡眠评估模块,用于基于得到的所述体动频次、所述呼吸混乱程度以及所述在床时长之中的一种或多种信息,进行被测者的睡眠质量评估。
第三方面,本发明实施例提供一种基于毫米波雷达的睡眠质量评估系统,包括:
设置于距离被测者的预设范围内的毫米波雷达监测设备;
与所述毫米波雷达监测设备通信连接的至少一个数据分析设备;
与所述至少一个数据分析设备通信连接的手持移动端;以及
与所述至少一个数据分析设备通信连接的存储器;
其中,所述存储器存储有可被所述至少一个数据分析设备执行的指令,所述指令被所述至少一个数据分析设备执行,以使所述至少一个数据分析设备能够执行如上所述的基于毫米波雷达的睡眠质量评估方法。
第四方面,本发明实施例提供一种计算机可读介质,其上存储有计算机可执行指令,所述可执行指令被处理器执行时实现如上所述的基于毫米波雷达的睡眠质量评估方法。
(三)有益效果
本发明的有益效果是:本发明基于毫米波雷达监测到的人体生命体征数据对被测者睡眠质量进行评估,监测手段具有非接触的优势,无需用户佩戴,因而不会干扰用户的正常睡眠。本发明基于被测者真实生命体征数据,依据被测者夜间睡眠时的在/离床状态、体动情况、呼吸周期平稳度这三个维度对被测者睡眠质量进行评估,能够客观反映被测者的睡眠质量。另外本发明借助于毫米波雷达非接触探测的优势,实施过程让用户无感知,这能够更客观真实地反映其睡眠质量,十分适合长期居家监测,便于形成连续监测数据,具备较大的应用前景与实际价值。
具体实施方式
为了更好地解释本发明,以便于理解,下面结合附图,通过具体实施方式,对本发明作详细描述。
如图1所示,本发明实施例提出的一种基于毫米波雷达的睡眠质量评估方法,其包括:首先,通过毫米波雷达监测设备采集被测者的生命体征数据;其次,依据获取的生命体征数据,通过构建生命体征数据曲线执行被测者的体动频次分析、呼吸混乱程度分析以及在床时长分析的一种或多种操作;最后,基于得到的体动频次、呼吸混乱程度以及在床时长之中的一种或多种信息,进行被测者的睡眠质量评估。
本发明基于毫米波雷达监测到的人体生命体征数据对被测者睡眠质量进行评估,监测手段具有非接触的优势,无需用户佩戴,因而不会干扰用户的正常睡眠。本发明基于被测者真实生命体征数据,依据被测者夜间睡眠时的在/离床状态、体动情况、呼吸周期平稳度这三个维度对被测者睡眠质量进行评估,能够客观反映被测者的睡眠质量。另外本发明借助于毫米波雷达非接触探测的优势,实施过程让用户无感知,这能够更客观真实地反映其睡眠质量,十分适合长期居家监测,便于形成连续监测数据,具备较大的应用前景与实际价值。
为了更好地理解上述技术方案,下面将参照附图更详细地描述本发明的示例性实施例。虽然附图中显示了本发明的示例性实施例,然而应当理解,可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了能够更清楚、透彻地理解本发明,并且能够将本发明的范围完整的传达给本领域的技术人员。
具体地,本发明提供的一种基于毫米波雷达的睡眠质量评估方法,包括:
S1、参考图2,本发明通过预先设于距离被测人员的预设范围内的毫米波雷达监测设备采集被测者睡眠过程中的生命体征数据。用户按照要求安装好毫米波雷达监测设备,用户准备入睡前通过手机微信小程序或APP开启睡眠监测功能,睡眠结束后用户通过手机微信小程序或APP关闭睡眠监测功能,在睡眠监测期间,毫米波雷达监测设备将自动采集被测者的生命体征数据。
本发明中使用的毫米波雷达监测设备采集的生命体征数据包括:图3示出的在/离床状态、图4示出的体动功率、图5示出的呼吸波形。体动频率为毫米波雷达监测设备输出的一个变量,指的是感兴趣距离范围内雷达回波功率之和,感兴趣距离范围指的是雷达视线范围内人体所处的距离区间;呼吸波形为毫米波雷达监测设备输出的一个变量,指的是毫米波雷达设备利用特定的信号处理方法提取的反映人体呼吸运动的信号波形,该呼吸波形代表了人体胸腔相对于毫米波雷达监测设备的距离变化情况;在/离床标志为毫米波雷达监测设备输出的一个变量,指的是毫米波雷达设备通过一定的处理逻辑对监测对象是否在床上做出的判断,该值为0时表示监测对象在床上,该值为1时表示监测对象不在床上。
且体动功率、在/离床标志以及呼吸波形的在毫米波雷达监测设备侧的更新频次均为1秒钟。更新频次均为1秒钟指的是每隔1秒钟毫米波雷达设备输出一次这些变量,更新频次需要与这些变量的采样频率相区分,体动功率的采样频率是10Hz,呼吸波形的采样频率是10Hz,在/离床标志的采样频率是1Hz。
S2、针对获取的生命体征数据,执行被测者的体动频次分析、呼吸混乱程度分析以及在床时长确定的一种或多种操作。
进一步地,
如图6所示,体动频次分析包括:
S21a、对获取的体动功率数据进行处理,得到T 1分钟内体动频次曲线。
更进一步地,如图7所示,步骤S21a包括:
S21a-1、利用获取的体动功率对被测者的体动事件进行滑窗检测,采用有序统计CFAR检测器,并设置参考单元对应时间长度为36秒、保护单元对应时间长度为20秒以及滑窗步进为1秒,经处理得到体动事件检测二值结果。
S21a-2、针对体动事件检测二值结果采用5分钟窗长和1秒钟滑窗步进计算5分钟内体动次数,得到5分钟内体动频次曲线。
S22a、采用阈值法对T 1分钟内体动频次曲线上的各个取值点进行二值划分,将大于阈值N 1的点赋值为1,否则赋值为0,得到体动频次1值初步区段与体动频次0值初步区段。
S23a、将体动频次0值初步区段中持续时间小于T 2分钟的点重新赋值为1,得到体动频次重新赋值区段与体动频次0值最终区段,并将体动频次0值最终区段视为体动稀少区段。
S24a、基于体动频次1值初步区段与体动频次重新赋值区段,得体动频繁区段。
进一步地,
如图8所示,呼吸混乱程度分析包括:
S21b、对获取的呼吸波形数据进行处理,得到T 3分钟内呼吸混乱程度曲线。
更进一步地,如图9所示,步骤S21b包括:
S21b-1、采用20秒窗长和1秒步进对呼吸波形进行处理,得到细粒度的呼吸混乱程度曲线。
S21b-2、在每个处理窗口内,对当前20秒时长的呼吸波形数据进行快速傅里叶变换,得到当前20秒时长呼吸波形数据的傅里叶谱。
S21b-3、在傅里叶谱上0.1Hz~0.7Hz频率范围内寻找功率最强的谱峰点,根据下式计算谱峰功率占比。
式中,r为谱峰功率占比,p(f)为傅里叶谱,f为频率,f peak为最强的谱峰点所对应的频率点。
S21b-4、依据谱峰功率占比通过下式得到呼吸混乱程度值。
ChaosIndex=2/(1+exp(r/(1-r)))。
S21b-5、采用5分钟窗长和1秒步进对细粒度的呼吸混乱程度曲线进行平滑处理,得到5分钟内呼吸混乱程度曲线。
S22b、采用阈值法对T 3分钟内呼吸混乱程度曲线上的各个取值点进行二值划分,将低于阈值τ 1的点赋值为1,否则赋值为0,得到呼吸混乱程度1值初步区段与呼吸混乱程度0值初步区段。
S23b、将呼吸混乱程度1值初步区段中持续时间小于分钟的点重新赋值为0,得到呼吸混乱程度重新赋值区段与呼吸混乱程度1值最终区段,并将呼吸混乱程度1值最终区段视为呼吸平稳区段。
S24b、基于呼吸混乱程度0值初步区段与呼吸混乱程度重新赋值区段,得到呼吸非平稳区段。
其中,T 1=2~10min、T 2=2~10min、T 3=2~5min、T 4=2~5min、N 1=3~15次、τ 1=0.1~0.3。
以及,在床时长分析包括:依据预设监测时间段内的在/离床标志,得到在床时长。
S3、依据体动频次、呼吸混乱程度以及在床时长之中的一种或多种信息,进行被测者的睡眠质量评估。
如图10所示,步骤S3包括:
S31、依据得到的体动频次分析结果得到体动稀少时长。
S32、依据得到的呼吸混乱程度分析结果得到呼吸平稳时长。
S33、依据得到在床时长并结合监测时间段内用户首次进入在床状态到其进入睡眠状态的时间长度得到入睡时长;若中间出现离床,离床时长不计入入睡时长。其中,进入睡眠状态的时间点指的是被测者在床状态内首次体动稀少时长超过10分钟的起始时刻,若中间出现离床,离床时长不计入入睡时长。
S34、依据在床时长、入睡时长、体动稀少时长以及呼吸平稳时长,通过如下睡眠质量评估模型进行被测者的睡眠质量评估:
其中,SleepQuality表示被测者的睡眠质量评估值,T Bed为在床时长、T Sopite 为入睡时长、T Quiet 为体动稀少时长、T Steady 为呼吸平稳时长。
另外地,本发明还提供一种基于毫米波雷达的睡眠质量评估系统,包括:
数据获取模块,用于通过毫米波雷达监测设备采集被测者的生命体征数据。
数据处理模块,用于依据获取的生命体征数据,通过构建生命体征数据曲线执行被测者的体动频次分析、呼吸混乱程度分析以及在床时长分析的一种或多种操作;
睡眠评估模块,用于基于得到的体动频次、呼吸混乱程度以及在床时长之中的一种或多种信息,进行被测者的睡眠质量评估。
由于本发明上述实施例所描述的系统/装置,为实施本发明上述实施例的方法所采用的系统/装置,故而基于本发明上述实施例所描述的方法,本领域所属技术人员能够了解该系统/装置的具体结构及变形,因而在此不再赘述。凡是本发明上述实施例的方法所采用的系统/装置都属于本发明所欲保护的范围。
同时,本发明还提供一种基于毫米波雷达的睡眠质量评估设备,包括:设置于距离被测者的预设范围内的毫米波雷达监测设备;
与毫米波雷达监测设备通信连接的至少一个数据分析设备;
与至少一个数据分析设备通信连接的手持移动端;以及
与至少一个数据分析设备通信连接的存储器;
其中,存储器存储有可被至少一个数据分析设备执行的指令,指令被至少一个数据分析设备执行,以使至少一个数据分析设备能够执行如上述的一种基于毫米波雷达的睡眠质量评估方法。
以及,本发明提供一种计算机可读介质,其上存储有计算机可执行指令,可执行指令被处理器执行时实现如上所述的一种基于毫米波雷达的睡眠质量评估方法。
综上所述,本发明提供一种基于毫米波雷达的睡眠质量评估方法、系统、设备及介质,如图11所示,本发明包括如下具体步骤:
(1)使用毫米波雷达监测设备采集被测者整晚的睡眠数据。
(2)对体动功率数据进行加工,得到5分钟内体动频次曲线。首先利用体动功率数据对体动事件进行滑窗检测,采用有序统计CFAR检测器,参考单元对应时间长度为36秒,保护单元对应时间长度为20秒,滑窗步进为1秒,经处理后得到体动事件检测二值结果,如图4所示。在体动事件检测二值结果数据上采用5分钟窗长和1秒钟滑窗步进计算5分钟内体动次数,得到5分钟内体动频次曲线。
(3)体动频繁/稀少段划分。首先采用阈值法对5分钟内体动频次进行二值划分,将5分钟内体动频次大于阈值的点赋值为1,否则赋值为0。然后将连续0值段中持续时间小于5分钟的段内的点重新赋值为1。至此,1代表体动频繁,0代表体动稀少。
(4)对呼吸波形数据进行加工,得到5分钟内呼吸混乱程度曲线。首先采用20秒窗长和1秒步进对呼吸波形数据进行处理,得到细粒度的呼吸混乱程度曲线;在每个处理窗口内,对当前20秒时长的呼吸波形数据进行快速傅里叶变换,得到当前20秒时长呼吸波形数据的傅里叶谱,在0.1Hz~0.7Hz频率范围内寻找功率最强的谱峰点机器对应频率点,并计算计算谱峰功率占比,利用谱峰功率占比通过得到呼吸混乱程度值;然后对细粒度的呼吸混乱程度曲线进行平滑处理得到5分钟内呼吸混乱程度曲线,采用5分钟窗长和1秒步进对细粒度的呼吸混乱程度曲线进行处理,在每个处理窗口内对窗口内5分钟时长的呼吸混乱程度值求平均值,作为当前窗口输出值。
(5)呼吸平稳/非平稳段划分。首先采用阈值法对5分钟内呼吸混乱程度曲线进行二值划分,将5分钟内呼吸混乱程度低于阈值的点赋值为1,否则赋值为0。然后将连续1值段中持续时间小于5分钟的段内的点重新赋值为0。至此,1代表呼吸平稳,0代表呼吸非平稳。
(6)在床时长提取。本发明中在床时长指的是监测时间段内用户在床上的时间长度,根据监测时间段内的在/离床标志,计算得到在床时长。
(7)入睡时长提取。本发明中入睡时长指的是监测时间段内用户首次进入在床状态到其进入睡眠状态的时间长度,如果中间出现离床,离床时长不计入入睡时长。其中进入睡眠状态时间点指的是被测者在床段内首次体动稀少段持续时间超过10分钟的起始时刻。
(8)体动稀少时长提取。利用步骤(3)中得到的体动频繁/稀少段划分结果计算体动稀少段的总时长,该值为体动稀少时长。
(9)呼吸平稳时长提取。利用步骤S5中得到的呼吸平稳/非平稳划分结果计算呼吸平稳段的总时长,该值为呼吸平稳时长。
(10)睡眠质量评估。本发明将被测者睡眠质量在0~100之间进行量化,采用睡眠质量评估数学模型结合步骤(1)—(9)得到的在床时长、入睡时长、体动稀少时长、呼吸平稳时长这些信息进行被测者的睡眠质量评估。
基于上述描述,可知鉴于毫米波雷达具有非接触探测的优势,无需用户佩戴即可探测用户的生命体征,本发明利用毫米波雷达采集的被测者睡眠过程中的生命体征数据对被测者的睡眠质量进行评估,在实现无感知探测的前提下,还能得到一个详细且精准的评估数据,具备较大的应用前景。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例,或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。
应当注意的是,在权利要求中,不应将位于括号之间的任何附图标记理解成对权利要求的限制。词语“包含”不排除存在未列在权利要求中的部件或步骤。位于部件之前的词语“一”或“一个”不排除存在多个这样的部件。本发明可以借助于包括有若干不同部件的硬件以及借助于适当编程的计算机来实现。在列举了若干装置的权利要求中,这些装置中的若干个可以是通过同一个硬件来具体体现。词语第一、第二、第三等的使用,仅是为了表述方便,而不表示任何顺序。可将这些词语理解为部件名称的一部分。
此外,需要说明的是,在本说明书的描述中,术语“一个实施例”、“一些实施例”、“实施例”、“示例”、“具体示例”或“一些示例”等的描述,是指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外,在不相互矛盾的情况下,本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
尽管已描述了本发明的优选实施例,但本领域的技术人员在得知了基本创造性概念后,则可对这些实施例作出另外的变更和修改。所以,权利要求应该解释为包括优选实施例以及落入本发明范围的所有变更和修改。
显然,本领域的技术人员可以对本发明进行各种修改和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也应该包含这些修改和变型在内。