一种HRV测量方法以及装置制造方法及图纸

本发明专利技术提供了一种HRV测量方法以及装置，其中，该方法包括：获取原始RR间期数据集，所述原始RR间期数据集中包括多个RR间期数据；对所述原始RR间期数据集做一次滤波处理，获取窦性心律所对应的RR间期数据集；对所述RR间期数据集进行快速傅立叶变换和/或小波变换，并根据快速傅立叶变换和/或小波变换的结果，对所述RR间期数据集进行去误差处理，获取标准RR间期数据集；根据所述标准RR间期数据集计算HRV时域指标。该HRV计算方法在保证计算精度的前提下，缩减了测量需要用到的时间。

HRV measuring method and device

The invention provides a HRV measuring method and device, wherein, the method comprises: obtaining the original RR interval data sets, including a plurality of RR interval data of the original data set RR interval; a filtering of the original data set RR interval, obtain sinus rhythm of RR interval data should be set; fast Fu Liye transform and / or wavelet transform on the RR interval data set, and according to the fast Fu Liye transform and / or the result of wavelet transform, to error processing on the RR interval RR interval data set, obtain the standard data set; according to the standard of RR interval data sets to calculate HRV time domain index. The HRV calculation method reduces the time needed for measurement under the premise of ensuring the accuracy of calculation.

【技术实现步骤摘要】
一种HRV测量方法以及装置
本专利技术涉及计算机应用
，具体而言，涉及一种HRV测量方法以及装置。
技术介绍
HRV(Heartratevariability，心率变异性)是指心率节奏快慢随时间所发生的变化。经临床实践证明，HRV可作为心源性猝死危险性的一个独立预测指标，在评价心血管疾病的预后和预测急性心血管事件方面具有重要的意义。现有的HRV测量方法包括短程测试和长程测试。短程测试是指在短时间内(例如5分钟)通过专有设备进行测量，该方法虽然具备测量时间短，使用方便的有点，但是测量出来的结果数据波动大，重复性差，且结果误差和较大，往往只有研究领域才会使用该方法对HRV进行测量。而长程测试虽然精确度比较高，但是一般的测试时间需要24小时左右，被测量的患者需要全天佩戴动态心电图监护仪，用户的很多行动都会受到限制。因此，一种耗时短、测量精确的心率变异性测量方法成为目前亟待解决的问题。
技术实现思路
有鉴于此，本专利技术实施例的目的在于提供一种心率变异性的测量方法以及装置，能够在较短的时间内，以更高的准确性测量到HRV。第一方面，本专利技术实施例提供了一种HRV测量方法，包括：获取原始RR间期数据集，所述原始RR间期数据集中包括多个RR间期数据；对所述原始RR间期数据集做一次滤波处理，获取窦性心律所对应的RR间期数据集；对所述RR间期数据集进行快速傅立叶变换和/或小波变换，并根据快速傅立叶变换和/或小波变换的结果，对所述RR间期数据集进行去误差处理，获取标准RR间期数据集；根据所述标准RR间期数据集计算HRV时域指标。结合第一方面，本专利技术实施例提供了第一方面的第一种可能的实施方式，其中：所述对所述原始RR间期数据集做一次滤波处理，获取窦性心律所对应的RR间期数据集，具体包括：计算所述原始RR间期数据集中每相邻的两个RR间期数据之间的差值；判断所述差值的绝对值是否小于等于预设的第一阈值；如果否，则将该相邻的两个所述RR间期数据从所述原始RR间期数据集中剔除，形成所述RR间期数据集。结合第一方面，本专利技术实施例提供了第一方面的第二种可能的实施方式，其中：所述对所述RR间期数据集进行快速傅立叶变换和/或小波变换，并根据快速傅立叶变换和/或小波变换的结果，对所述RR间期数据集进行去误差处理，获取标准RR间期数据集，具体包括：对所述RR间期数据集进行快速傅立叶变换，分别获取HRV的四个标准频率段ULF、VLF、LF以及HF的能量：EULF、EVLF、ELF以及EHF；判断EULF+EVLF是否小于预设的阈第二值，和/或，判断(EULF+EVLF)/(EULF+EVLF+ELF+EHF)*100％是否小于等于预设的第三阈值；如果是，则将所述RR间期数据集作为标准RR间期数据集；如果否，则对所述RR间期数据集进行小波变换，获取多个阶层的波形；按照从高阶至低阶，或者从低阶至高阶的顺序，执行下述处理过程：对当前阶层的波形进行快速傅立叶变换，获取当前阶层中HRV的四个标准频率段：ULF、VLF、LF以及HF的能量EULF、EVLF、ELF以及EHF判断ULF+VLF是否小于预设的阈第二值，和/或，判断(ULF+VLF)/(ULF+VLF+LF+HF)*100％是否小于等于预设的第三阈值；如果是，则将当前阶层所对应的RR间期数据集作为标准间期数据集。结合第一方面，本专利技术实施例提供了第一方面的第三种可能的实施方式，其中：所述根据所述标准RR间期数据集计算HRV时域指标具体包括：根据所述标准RR间期数据集，以RR间期作为纵坐标，以i为横坐标，绘制HRV信号曲线；根据下述公式计算HRV时域技术指标：NNVGR为全部正常NN间期的平均值；SDNN为全部正常NN间期的标准差；RNSSD为全部相邻NN间期之差的均方根值；SDSD为全部相邻NN间期之差的标准差；NN50，为全部相邻NN间期之差大于50ms的心搏数；PNN50为全部相邻NN间期之差大于50ms的心搏数的百分比；在上述公式中，N为正常心搏总数，RRi是所述标准RR间期数据集中的第i个RR间期数据，RRi+1是第i+1个RR间期数据。RR是N个心搏的RR间期数据的平均值；NN间期为窦性心搏间期；RRi′＝RRi-RRi-1；结合第一方面，本专利技术实施例提供了第一方面的第四种可能的实施方式，其中：所述对所述RR间期数据集进行快速傅立叶变换和小波变换之前，还包括：对所述RR间期数据集进行可用性检测。结合第一方面，本专利技术实施例提供了第一方面的第五种可能的实施方式，其中：所述对所述RR间期数据集进行可用性检测具体包括：按照时间轴顺序将所述RR间期数据集中的RR间期数据进行分组，形成多个RR间期数据组；每个所述RR间期数据组中包含预设数量的RR间期数据；计算每个所述RR间期数据组中所有RR间期数据的平均值；根据每个所述RR间期数据组中所有RR间期数据的平均值，计算该平均值所对应的心率值；判断所有心率值中最大值与最小值之差是否大于预设的数值，如果是，则数据不可用；和/或，按照时间轴顺序将所述RR间期数据集中的RR间期数据进行分组，形成第一分组以及第二分组；其中，所述第一分组中的RR间期数据的数量和所述第二分组中RR间期数据的数量不同；计算第一分组以及第二分组中RR间期数据差值的平均值；根据所述第一分组以及第二分组中RR间期数据差值的平均值计算检测准确值；判断所述准确标准值是否在预设范围之内，如果否，则数据不可用。结合第一方面，本专利技术实施例提供了第一方面的第六种可能的实施方式，其中：所述对所述RR间期数据集进行快速傅立叶变换和/或小波变换之前，还包括：对一次滤波所获得的窦性心律所对应的RR间期数据集做二次滤波处理；其中，所述二次滤波处理按照如下方式进行：计算所述RR间期数据集中所有RR间期数据的标准差；使用所述RR间期数据的标准差，消除RR间期数据集中的峰值数据。结合第一方面，本专利技术实施例提供了第一方面的第七种可能的实施方式，其中：所述使用所述RR间期数据的标准差，消除RR间期数据集中的峰值数据，具体包括：根据下述公式计算RR间期数据取值区间的两个端点RR下端以及RR上端：RR下端＝RR平-3×RR标；RR上端＝RR平+3×RR标；其中，RR平为RR间期数据集中所有RR间期数据的平均值；所述RR标为RR间期数据的标准差；判断RR间期数据集中的每一个RR间期数据是否落入该RR间期数据取值区间以内，如果否，则将该RR间期数据从所述RR间期数据集中剔除。第二方面，本专利技术实施例还提供一种HRV测量装置，包括：数据获取模块，用于获取原始RR间期数据集，所述原始RR间期数据集中包括多个RR间期数据；一次滤波模块，用于对所述原始RR间期数据集做一次滤波处理，获取窦性心律所对应的RR间期数据集；误差处理模块，用于对所述RR间期数据集进行快速傅立叶变换和/或小波变换，并根据快速傅立叶变换和/或小波变换的结果，对所述RR间期数据集进行去误差处理，获取标准RR间期数据集；HRV计算模块，用于根据所述标准RR间期数据集计算HRV时域指标。结合第二方面，本专利技术实施例提供了第二方面的第一种可能的实施方式，其中：还包括：数据可用性检测模块，用于对所述RR间期数据集进行可用性检测；和/或，还包括：二次滤波模块，用于对一次滤波所获本文档来自技高网...

【技术保护点】
一种HRV测量方法，其特征在于，包括：获取原始RR间期数据集，所述原始RR间期数据集中包括多个RR间期数据；对所述原始RR间期数据集做一次滤波处理，获取窦性心律所对应的RR间期数据集；对所述RR间期数据集进行快速傅立叶变换和/或小波变换，并根据快速傅立叶变换和/或小波变换的结果，对所述RR间期数据集进行去误差处理，获取标准RR间期数据集；根据所述标准RR间期数据集计算HRV时域指标。

【技术特征摘要】
1.一种HRV测量方法，其特征在于，包括：获取原始RR间期数据集，所述原始RR间期数据集中包括多个RR间期数据；对所述原始RR间期数据集做一次滤波处理，获取窦性心律所对应的RR间期数据集；对所述RR间期数据集进行快速傅立叶变换和/或小波变换，并根据快速傅立叶变换和/或小波变换的结果，对所述RR间期数据集进行去误差处理，获取标准RR间期数据集；根据所述标准RR间期数据集计算HRV时域指标。2.根据权利要求1所述的方法，其特征在于，所述对所述原始RR间期数据集做一次滤波处理，获取窦性心律所对应的RR间期数据集，具体包括：计算所述原始RR间期数据集中每相邻的两个RR间期数据之间的差值；判断所述差值的绝对值是否小于等于预设的第一阈值；如果否，则将该相邻的两个所述RR间期数据从所述原始RR间期数据集中剔除，形成所述RR间期数据集。3.根据权利要求2所述的方法，其特征在于，所述对所述RR间期数据集进行快速傅立叶变换和/或小波变换，并根据快速傅立叶变换和/或小波变换的结果，对所述RR间期数据集进行去误差处理，获取标准RR间期数据集，具体包括：对所述RR间期数据集进行快速傅立叶变换，分别获取HRV的四个标准频率段ULF、VLF、LF以及HF的能量：EULF、EVLF、ELF以及EHF；判断EULF+EVLF是否小于预设的阈第二值，和/或，判断(EULF+EVLF)/(EULF+EVLF+ELF+EHF)*100％是否小于等于预设的第三阈值；如果是，则将所述RR间期数据集作为标准RR间期数据集；如果否，则对所述RR间期数据集进行小波变换，获取多个阶层的波形；按照从高阶至低阶，或者从低阶至高阶的顺序，执行下述处理过程：对当前阶层的波形进行快速傅立叶变换，获取当前阶层中HRV的四个标准频率段：ULF、VLF、LF以及HF的能量EULF、EVLF、ELF以及EHF判断ULF+VLF是否小于预设的阈第二值，和/或，判断(ULF+VLF)/(ULF+VLF+LF+HF)*100％是否小于等于预设的第三阈值；如果是，则将当前阶层所对应的RR间期数据集作为标准间期数据集。4.根据权利要求1-3任意一项所述的方法，其特征在于，所述根据所述标准RR间期数据集计算HRV时域指标具体包括：根据所述标准RR间期数据集，以RR间期作为纵坐标，以i为横坐标，绘制HRV信号曲线；根据下述公式计算HRV时域技术指标：NNVGR为全部正常NN间期的平均值；SDNN为全部正常NN间期的标准差；RNSSD为全部相邻NN间期之差的均方根值；SDSD为全部相邻NN间期之差的标准差；NN50，为全部相邻NN间期之差大于50ms的心搏数；PNN50为全部相邻NN间期之差大于50ms的心搏数的百分比；在上述公式中，N为正常心搏总数，RRi是所述标准RR间期数据集中的第i个RR间期数据，RRi+1是第i+1个RR间期数据。是N个心搏的RR间期数据的平均值；NN间期为窦性心搏间期；RR′i＝RRi-RRi-1；

【专利技术属性】
技术研发人员：陈远翔，马斌，
申请(专利权)人：北京心量科技有限公司，
类型：发明
国别省市：北京,11