CN106264475A - 单一光电传感器睡眠呼吸多生理参数监测方法及装置 - Google Patents

单一光电传感器睡眠呼吸多生理参数监测方法及装置 Download PDF

Info

Publication number
CN106264475A
CN106264475A CN201610893587.2A CN201610893587A CN106264475A CN 106264475 A CN106264475 A CN 106264475A CN 201610893587 A CN201610893587 A CN 201610893587A CN 106264475 A CN106264475 A CN 106264475A
Authority
CN
China
Prior art keywords
pulse frequency
photoplethysmographic
module
detection
information
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201610893587.2A
Other languages
English (en)
Other versions
CN106264475B (zh
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.)
Shaanxi Huachen Petroleum Technology Co ltd
Suzhou Sanlang Medical Technology Co ltd
Original Assignee
Xian Jiaotong University
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 Xian Jiaotong University filed Critical Xian Jiaotong University
Priority to CN201610893587.2A priority Critical patent/CN106264475B/zh
Publication of CN106264475A publication Critical patent/CN106264475A/zh
Application granted granted Critical
Publication of CN106264475B publication Critical patent/CN106264475B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • A61B5/02405Determining heart rate variability
    • 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
    • A61B5/02416Detecting, measuring or recording pulse rate or heart rate using photoplethysmograph signals, e.g. generated by infrared radiation
    • 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
    • A61B5/02416Detecting, measuring or recording pulse rate or heart rate using photoplethysmograph signals, e.g. generated by infrared radiation
    • A61B5/02427Details of sensor
    • A61B5/02433Details of sensor for infrared radiation
    • 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/026Measuring blood flow
    • A61B5/0261Measuring blood flow using optical means, e.g. infrared light
    • 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/026Measuring blood flow
    • A61B5/0295Measuring blood flow using plethysmography, i.e. measuring the variations in the volume of a body part as modified by the circulation of blood therethrough, e.g. impedance plethysmography
    • 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
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/103Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/11Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/145Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
    • A61B5/1455Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters
    • A61B5/14551Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters for measuring blood gases
    • 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/4818Sleep apnoea
    • 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

Landscapes

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

Abstract

本发明单一光电传感器睡眠呼吸多生理参数监测方法及装置,能够在仅使用单一光电传感器检测的基础上提取得到五项生理参数。所述方法包括,步骤1,采集光电容积脉搏波;步骤2,对采集到的光电容积脉搏波进行预处理;步骤3,根据检测的峰值点、谷值点以及基线变化信息判断当前测量时刻是否处于运动状态;当处于非运行状态,将识别的非运动状态数据视为正常光电容积脉搏波;步骤4,当检测到当前时刻处于运动状态,则进行体动次数计数;步骤5,根据正常光电容积脉搏波,检测脉率信息、血氧饱和度信息、呼吸信息和脉率变异性信息;步骤6,根据步骤4和5提取脉率、血氧饱和度、呼吸、脉率变异性以及体动次数信息,进行数据传输、显示和记录。

Description

单一光电传感器睡眠呼吸多生理参数监测方法及装置
技术领域
本发明涉及生物医学工程领域,具体为单一光电传感器睡眠呼吸多生理参数监测方法及装置。
背景技术
随着现代社会节奏的不断加快,睡眠质量越来越严重地影响到个人生理和精神各方面的健康状况。大量文献表明,约有1/3的高血压病、1/5的心脏病是由睡眠呼吸障碍引起的,其它睡眠呼吸障碍并发症包括糖尿病、肾病、肥胖、癫痫、性功能减退、植物神经功能紊乱、红细胞增多症、严重的还有心肌梗塞、脑血栓、脑卒等。因此,对睡眠呼吸障碍的监测对临床监护和病理生理研究有着十分重要的意义。
对睡眠呼吸障碍的常规监测方法是采用多导睡眠监护仪对测量对象进行睡眠监护,以数量众多的不同传感器同时获取测量对象血氧、脉搏、血压、心电、脑电、肌电、眼电、口鼻气流、气道压力、胸腹呼吸、鼾声等多项生理参数以进行综合分析。通常每项生理参数的检测都需要使用一种与之相配合的传感器。然而睡眠呼吸障碍监测领域一直存在一对难以调和的矛盾:一方面,睡眠是涉及众多生理变化的复杂行为,只有检测尽可能多的生理参数才能全面、客观地反映和检测睡眠呼吸障碍;另一方面:为不影响、干扰正常的睡眠,希望由测量引入的导联、传感器的数目越少越好。这一矛盾导致临床上不得不折衷考虑检测参数和使用传感器的数量。为此,美国睡眠医学会(American Academy of Sleep Medicine,AASM)专门制定了用于睡眠呼吸暂停综合征监测的四级睡眠监护设备标准,其中I级为标准多导睡眠仪;II级为至少能够检测七个生理参数的综合便携式多导睡眠仪;III级为至少能够检测呼吸、心率以及血氧饱和度等四个生理参数的改良型便携式睡眠呼吸障碍检测设备;IV级为连续单双生理指标记录仪。四级睡眠监护设备按照检测参数的数目分别使用在从大型医院睡眠监护中心到中小型医院、社区诊所以及家庭应用等不同的场所。
作为IV级连续单双生理指标记录仪的典型代表,仅使用光电传感器检测光电容积脉搏波的脉搏血氧仪目前已在临床监护领域广泛应用,可用于检测脉率(心率)和血氧饱和度。
光电容积脉搏波描记法是使用光电传感器从活体组织中检测毛细血管床容积变化的一种无创检测方法。其核心原理是1851年August Beer(德)发表的光吸收的基本定律,即Lambert-Beer定律。受心脏搏动的影响,毛细血管床容积呈现出周期性变化,以Lambert-Beer定律检测出其变化周期即可相应地检测出脉率,而在能够检测出脉率的情况下,其数值与心率几乎是完全相等的,此即为使用光电传感器检测脉率(心率)的原理。1935年Matthes在Lambert-Beer定律基础上结合血液中氧合血红蛋白与脱氧血红蛋白对不同波长的光吸收系数的不同的现象,提出了使用光电传感器无创检测血氧饱和度的检测方法,并使用至今。
但是,仅使用光电传感器检测光电容积脉搏波的脉搏血氧仪,在成本低廉、使用方便的同时也因为检测参数数目过少而限制了其单独在大型医院睡眠监护中心以及中小型医院的使用。
发明内容
针对现有技术中存在的问题,本发明提供一种单一光电传感器睡眠呼吸多生理参数监测方法及装置,能够在仅使用单一光电传感器检测光电容积脉搏波的基础上,从原始光电容积波中提取得到脉率、血氧饱和度、呼吸、脉率变异性以及体动五项生理参数,以低成本实现睡眠呼吸多个生理参数的同时监测。
本发明是通过以下技术方案来实现:
单一光电传感器睡眠呼吸多生理参数监测方法,包括下述步骤:
步骤1,以模拟方式或数字方式采集采样率在500Hz以上的光电容积脉搏波;光电容积脉搏波包含呼吸信息和未经滤波处理的原始光电容积脉搏波信号,原始光电容积脉搏波信号包含以反射或者透射方式穿过被检组织的红光和红外光;
步骤2,对采集到的光电容积脉搏波进行预处理,包括基线的提取检测、峰值点与谷值点的提取与检测;
步骤3,根据检测的峰值点、谷值点以及基线变化信息判断当前测量时刻是否处于运动状态;
当峰值点、谷值点以及基线变化幅度超过所设定的阈值时认为当前时刻处于运动状态,将丢弃当前测量值重新测量,直至当前测量值脱离运动状态进入正常光电容积脉搏波范围内,将识别的非运动状态数据视为正常光电容积脉搏波;
步骤4,当检测到当前时刻处于运动状态,则进行体动次数计数;根据提取的运动状态信息,将所检测到的运动状态信息在当前运动停止后记为体动增加1次,累加至体动次数计数值;
步骤5,根据正常光电容积脉搏波,检测脉率信息、血氧饱和度信息、呼吸信息和脉率变异性信息;
步骤6,根据步骤4和5提取脉率、血氧饱和度、呼吸、脉率变异性以及体动次数信息,进行数据传输、显示和记录。
优选的,所述步骤2中,峰值点、谷值点以及基线由公式(1)、(2)、(3)提取和检测;
Peakvalue=max(x(k:k+100)),1≤k≤N-100 (1)
Vallyvalue=min(x(k:k+100)),1≤k≤N-100 (2)
Base=(Sup+Sdown)/2 (3)
其中,x为红光或红外光的光电容积脉搏波信号,N为所采数据样点的总数,k为1至N-100之间的自然数,Sup,Sdown分别为光电容积脉搏波信号的上包络和下包络。
优选的,所述步骤3中峰值点、谷值点以及基线变化阈值由实验确定,运动阈值为10~20%。
优选的,所述步骤5中,脉率由如下公式(4)提取和检测;
HR=60/RR (4)
其中,RR为脉率间期,单位为秒。
优选的,所述步骤5中,血氧饱和度由如下公式(5)或(6)提取和检测;
SpO2=α-βR (5)
SpO2=α-βR-γR2 (6)
其中, 分别为红光和红外光强度的交流分量,分别为红光和红外光强度的直流分量;α、β以及γ分别为经验值,通过定标来确定。
优选的,所述步骤5中,呼吸的提取和检测采用数字滤波法,其输入输出关系如下公式(8)所示;
y ( n ) = { Σ k = 0 N a k x ( n - k ) } - { Σ k = 1 M b k y ( n - k ) } - - - ( 8 )
其中,ak和bk为IIR数字滤波器的系数,N、M分别为数字滤波器系统传递函数的分子和分母多项式系数阶数,n为自然数;光电容积脉搏波通过数字滤波得到呼吸波波形,然后通过过零法或峰值法检测计算得到呼吸率。
优选的,所述步骤5中,脉率变异性的提取和检测采用脉率变异性的非线性分析方法,计算脉率间期彭家莱图的标准差SD1、SD2以及二者的比值SD1/SD2,具体如下公式(9)和(10);
S D 1 2 = 1 N - 1 Σ i = 1 N - 1 ( D i ′ ) 2 - - - ( 9 )
S D 2 2 = 1 N - 1 Σ i = 1 N - 1 ( D i ′ ′ ) 2 - - - ( 10 )
其中,N为RR间期总数,RRn和RRn+1分别为第n与第n+1个脉率间期,为全部脉率间期序列的间期均值,n为自然数。
单一光电传感器睡眠呼吸多生理参数监测装置,包括下述模块:
原始数据采集模块,用于以模拟方式或数字方式采集采样率在500Hz以上的光电容积脉搏波;
预处理模块,用于对原始数据采集模块采集到的光电容积脉搏波进行预处理,包括基线的提取检测、峰值点与谷值点的提取与检测;
运动识别模块,用于根据预处理模块检测的峰值点、谷值点以及基线变化信息判断当前测量时刻是否处于运动状态;当峰值点、谷值点以及基线变化幅度超过所设定的阈值时认为当前时刻处于运动状态,将丢弃当前测量值重新测量,直至当前测量值脱离运动状态进入正常光电容积脉搏波范围内;所检测到的运动信息送至体动检测模块进行体动次数计数;
正常光电容积脉搏波模块,用于将运动识别模块识别的非运动状态数据视为正常光电容积脉搏波,用于后续信号处理;
脉率检测模块、血氧饱和度检测模块、呼吸检测模块和脉率变异性检测模块,分别用于根据正常光电容积脉搏波模块提取的正常光电容积脉搏波,检测脉率信息、血氧饱和度信息、呼吸信息和脉率变异性信息,并输送至数据后处理模块;
体动检测模块,用于根据运动识别模块提取的运动信息,将所检测到的运动信息在当前运动停止后记为体动增加1次,累加至体动次数计数值,并输送至数据后处理模块;
数据后处理模块,对脉率检测模块、血氧饱和度检测模块、呼吸检测模块、脉率变异性检测模块和体动检测模块提取的脉率、血氧饱和度、呼吸、脉率变异性以及体动次数信息,与上位机进行数据传输、显示和记录。
优选的,所述原始数据采集模块包括发光驱动、红光光二极管、红外光发光二极管和光频转换器;外部微处理器控制的发光驱动用于交替产生驱动电流,分别驱动红光发光二极管以及红外光发光二极管;光频转换器用于将所发出的红光、红外光以反射或者透射方式穿过被测组织的红光和红外光的光强转化为频率信号,输出到微处理器捕捉计数。
优选的,所述数据后处理模块包括无线蓝牙数据传输模块、MicroSD卡数据存储模块和显示模块。
与现有技术相比,本发明具有以下有益的技术效果:
由于毛细血管床容积变化同时受到呼吸系统、循环系统、神经系统、以及运动系统的共同影响,本发明在仅使用单一光电传感器检测光电容积脉搏波的基础上,以数字信号处理的方法及装置从光电容积脉搏波中提取上述系统状态特征中间值的多个生理参数,则能够在一定程度上平衡睡眠呼吸生理参数监测中传感器数目与检测参数数目间的矛盾。
本发明通过采集脉率和脉率变异性替代现有技术的心率和心率变异性的采集,省去了体表心电传感器;以呼吸对脉率参量的影响监测呼吸信息,省去了呼吸传感器;以运动对光电容积脉搏波波形的影响监测体动次数,不仅充分利用了现有技术中无法消除的运动干扰,还省去了运动传感器;使用数字式数据采集模块,简化了系统设计,提高了信噪比;使用数字信号处理的方式同时从光电容积脉搏波中提取检测多项生理参数,各检测模块之间互不影响,不会因为增加检测参数的数目而影响信噪比,降低检测质量。
采用本发明所述的方法和装置能够将现有的IV级连续单双生理指标记录仪的代表、基于光电容积脉搏波的脉搏血氧仪在几乎不增加生产成本的情况下升级至能够检测呼吸、心率以及血氧饱和度等四个以上生理参数的III级改良型便携式睡眠呼吸障碍检测设备。扩大了其使用范围,可以满足中小型医院、社区诊所以及家庭等不同的场所的应用。
附图说明
图1为本发明方法的步骤流程示意图。
图2为原始光电容积脉搏波变量示意图。
图3为本发明装置的模块结构示意图。
图4为原始数据采集模块结构示意图。
图5为本发明装置上位机软件分析截屏图。
具体实施方式
下面结合具体的实施例对本发明做进一步的详细说明,所述是对本发明的解释而不是限定。
本发明为了解决现有睡眠呼吸生理参数监测中传感器数目与检测参数数目间的矛盾,针对III级改良型便携式睡眠呼吸障碍检测设备的要求,提出了一种数字信号处理方法及装置,
本发明的发明点在于:仅使用单一光电传感器,以数字信号处理的方式同时从光电容积脉搏波中提取检测脉率、血氧饱和度、呼吸、脉率变异性以及体动等多项生理参数,较好地平衡了睡眠呼吸生理参数监测中传感器数目与检测参数数目间的矛盾。
参见图1,本发明的单一光电传感器睡眠呼吸多生理参数监测方法包括下述步骤:
步骤1,采集光电容积脉搏波,以模拟方式或数字方式,所采集的光电容积脉搏波是采样率在500Hz以上;优选地,所采集的光电容积脉搏波是以光频变换的数字方式采集,采样率为1KHz,且包含呼吸信息、未经滤波处理的原始光电容积脉搏波信号,信号包含以透射方式穿过被检组织的两种波长的光,两种波长分别为660nm红光和940nm红外光;
步骤2,对步骤1采集到的光电容积脉搏波进行预处理,包括基线的提取检测、峰值点与谷值点的提取与检测;
步骤3,根据步骤2检测的峰值点、谷值点以及基线变化信息判断当前测量时刻是否处于运动状态;当峰值点、谷值点以及基线变化幅度超过所设定的阈值时认为当前时刻处于运动状态,将丢弃当前测量值重新测量,直至当前测量值脱离运动状态进入正常光电容积脉搏波范围内,将识别的非运动状态数据视为正常光电容积脉搏波;所检测到的运动信息送至步骤4进行体动次数计数;
步骤4,当检测到当前时刻处于的运动状态,则进行体动次数计数;根据提取的运动状态信息,将所检测到的运动状态信息在当前运动停止后记为体动增加1次,累加至体动次数计数值;根据步骤3识别的非运动状态数据视为正常光电容积脉搏波,用于后续信号处理;
步骤5,根据步骤3提取的正常光电容积脉搏波,检测脉率信息、血氧饱和度信息、呼吸信息和脉率变异性信息并送至步骤6;
步骤6,根据步骤4和5提取的脉率、血氧饱和度、呼吸、脉率变异性以及体动次数信息,进行数据传输、显示和记录;优选的对超出安全阈值的信息进行声、光报警;
其中,步骤2中,所述峰值点、谷值点以及基线参见图2,由公式(1)、(2)、(3)描述:
Peakvalue=max(x(k:k+100)),1≤k≤N-100 (1)
Vallyvalue=min(x(k:k+100)),1≤k≤N-100 (2)
Base=(Sup+Sdown)/2 (3)
其中,x为红光或红外光的光电容积脉搏波信号,N为所采数据样点的总数,k为1至N-100之间的自然数,以100点为窗宽,在光电容积脉搏波信号上滑动,当最大值处于窗的中间位置时,认为此点为峰值点;同理,当最小值处于窗的中间位置时,认为此点为谷值点;Sup,Sdown分别为光电容积脉搏波信号的上包络和下包络。
步骤3中,峰值点、谷值点以及基线变化阈值由实验确定,优选地,运动阈值为:10~20%,即峰值点、谷值点以及基线变化中的任一项超过前次测量值的10~20%即视为测量对象进入体动状态;
步骤5中,脉率的提取与检测如公式(4)描述;
HR=60/RR (4)
其中,RR为脉率间期,单位为秒。
步骤5中,血氧饱和度的提取与检测如公式(5),(6),(7)描述;
SpO2=α-βR (5)
SpO2=α-βR-γR2 (6)
R = I A C λ 1 / I D C λ 1 I A C λ 2 / I D C λ 2 - - - ( 7 )
其中,分别为红光、红外光强度的交流分量,分别为红光、红外光强度的直流分量。因此,在一个心搏间期内获得了红光和红外光强度的直流和交流分量时,就可以计算出R值,进而计算出该心搏周期内的血氧饱和度。
α、β以及γ分别为经验值,通过定标来确定。公式(5)为血氧饱和度与R值的一阶线性关系,在一阶关系不能满足线性度要求时,可以采用二阶公式(6)来描述.
步骤5中,呼吸的提取与检测包括但不限于数字滤波法、数学形态学方法、小波变换法、时频联合分析法等,优选地,采用数字滤波法,其输入输出关系如公式(8)描述;
y ( n ) = { Σ k = 0 N a k x ( n - k ) } - { Σ k = 1 M b k y ( n - k ) } - - - ( 8 )
其中,ak和bk为IIR数字滤波器的系数,N、M分别为数字滤波器系统传递函数的分子和分母多项式系数阶数,n为自然数。原始光电容积脉搏波通过该带通数字滤波即可得到呼吸波波形,然后通过过零法或峰值法检测即可计算得到呼吸率。
步骤5中,脉率变异性的提取与检测包括但不限于时域分析法、频域分析法以及非线性分析法;优选地,采用脉率变异性的非线性分析方法,计算脉率间期彭家莱图的标准差SD1、SD2以及二者的比值SD1/SD2,如公式(9),(10)描述;
S D 1 2 = 1 N - 1 Σ i = 1 N - 1 ( D i ′ ) 2 - - - ( 9 )
S D 2 2 = 1 N - 1 Σ i = 1 N - 1 ( D i ′ ′ ) 2 - - - ( 10 )
其中,D′i和D″i如公式(11)、(12)描述,N为RR间期总数:
D i ′ = 1 2 | RR n - RR n + 1 | - - - ( 11 )
D i ′ ′ = 1 2 | RR n + RR n + 1 - 2 R R ‾ | - - - ( 12 )
其中,RRn和RRn+1分别为第n与第n+1个脉率间期,为全部脉率间期序列的间期均值,n为自然数:
所述步骤4中,体动次数计数值初值为0,所检测到的运动信息在当前运动停止后记为体动增加1次,累加至体动次数计数值;
所述步骤5中检测脉率信息、血氧饱和度信息、呼吸信息和脉率变异性信息为并列关系,执行时不存在先后关系;
所述步骤6中,脉率、血氧饱和度等参数超出安全阈值时进行声、光报警的阈值由实验确定,其安全阈值分别为:60(beat/min)≤脉率≤100(beat/min),血氧饱和度≥90(%)。
与之相应地,单一光电传感器睡眠呼吸多生理参数监测装置,参见图3,包括下述模块:
原始数据采集模块,用于以模拟方式或数字方式采集采样率在500Hz以上的光电容积脉搏波;优选的采用光频变换的数字方式;所采集的光电容积脉搏波是采样率为1KHz,且包含呼吸信息、未经滤波处理的原始光电容积脉搏波信号;信号包含以反射或透射方式穿过被检组织的两种波长的光,两种波长分别为660nm红光和940nm红外光;
预处理模块,对原始数据采集模块采集到的光电容积脉搏波进行预处理,包括基线的提取检测、峰值点与谷值点的提取与检测;
运动识别模块,根据预处理模块检测的峰值点、谷值点以及基线变化信息判断当前测量时刻是否处于运动状态,当峰值点、谷值点以及基线变化幅度超过所设定的阈值时认为当前时刻处于运动状态,将丢弃当前测量值重新测量,直至当前测量值脱离运动状态进入正常光电容积脉搏波范围内;所检测到的运动信息送至体动检测模块进行体动次数计数;
正常光电容积脉搏波模块,将运动识别模块识别的非运动状态数据视为正常光电容积脉搏波,用于后续信号处理;
脉率检测模块,根据正常光电容积脉搏波模块提取的正常光电容积脉搏波,检测脉率信息,并送至数据后处理模块;
血氧饱和度检测模块,根据正常光电容积脉搏波模块提取的正常光电容积脉搏波,检测血氧饱和度信息,并送至数据后处理模块;
呼吸检测模块,根据正常光电容积脉搏波模块提取的正常光电容积脉搏波,检测呼吸信息,并送至数据后处理模块;
脉率变异性检测模块,根据正常光电容积脉搏波模块提取的正常光电容积脉搏波,检测脉率变异性信息,并送至数据后处理模块;
体动检测模块,根据运动识别模块提取的运动信息,将所检测到的运动信息在当前运动停止后记为体动增加1次,累加至体动次数计数值,并送至数据后处理模块;
数据后处理模块,对脉率检测模块、血氧饱和度检测模块、呼吸检测模块、脉率变异性检测模块、体动检测模块提取的脉率、血氧饱和度、呼吸、脉率变异性以及体动次数信息,与上位机进行数据传输、显示、记录,并对超出安全阈值的信息进行声、光报警等数据后处理;
上述装置原始数据采集模块中,如图4所示,由发光驱动、红光、红外光发光二极管、光频转换器构成,其中发光驱动在外部微处理器的控制下交替产生驱动电流,分别驱动红光以及红外光发光二极管,所发出的红光、红外光以反射或者透射方式穿过被测组织,被光频转换器直接将光强转化为频率信号,供微处理器捕捉计数。这里微处理器采用MSP430,光频转换器采用TSL235。
上述装置数据后处理模块中,包括无线蓝牙数据传输、MicroSD卡数据存储和显示模块,以及蜂鸣器、LED报警灯构成的报警组件。
本发明所采集和检测的数据传输到上位机后得到的分析图像如图5所示。

Claims (10)

1.单一光电传感器睡眠呼吸多生理参数监测方法,其特征在于,包括下述步骤:
步骤1,以模拟方式或数字方式采集采样率在500Hz以上的光电容积脉搏波;光电容积脉搏波包含呼吸信息和未经滤波处理的原始光电容积脉搏波信号,原始光电容积脉搏波信号包含以反射或者透射方式穿过被检组织的红光和红外光;
步骤2,对采集到的光电容积脉搏波进行预处理,包括基线的提取检测、峰值点与谷值点的提取与检测;
步骤3,根据检测的峰值点、谷值点以及基线变化信息判断当前测量时刻是否处于运动状态;
当峰值点、谷值点以及基线变化幅度超过所设定的阈值时认为当前时刻处于运动状态,将丢弃当前测量值重新测量,直至当前测量值脱离运动状态进入正常光电容积脉搏波范围内,将识别的非运动状态数据视为正常光电容积脉搏波;
步骤4,当检测到当前时刻处于运动状态,则进行体动次数计数;根据提取的运动状态信息,将所检测到的运动状态信息在当前运动停止后记为体动增加1次,累加至体动次数计数值;
步骤5,根据正常光电容积脉搏波,检测脉率信息、血氧饱和度信息、呼吸信息和脉率变异性信息;
步骤6,根据步骤4和5提取脉率、血氧饱和度、呼吸、脉率变异性以及体动次数信息,进行数据传输、显示和记录。
2.如权利要求1所述的单一光电传感器睡眠呼吸多生理参数监测方法,其特征在于,所述步骤2中,峰值点、谷值点以及基线由公式(1)、(2)、(3)提取和检测;
Peakvalue=max(x(k:k+100)),1≤k≤N-100 (1)
Vallyvalue=min(x(k:k+100)),1≤k≤N-100 (2)
Base=(Sup+Sdown)/2 (3)
其中,x为红光或红外光的光电容积脉搏波信号,N为所采数据样点的总数,k为1至N-100之间的自然数,Sup,Sdown分别为光电容积脉搏波信号的上包络和下包络。
3.如权利要求1所述的单一光电传感器睡眠呼吸多生理参数监测方法,其特征在于,所述步骤3中峰值点、谷值点以及基线变化阈值由实验确定,运动阈值为10~20%。
4.如权利要求1所述的单一光电传感器睡眠呼吸多生理参数监测方法,其特征在于,所述步骤5中,脉率由如下公式(4)提取和检测;
HR=60/RR (4)
其中,RR为脉率间期,单位为秒。
5.如权利要求1所述的单一光电传感器睡眠呼吸多生理参数监测方法,其特征在于,所述步骤5中,血氧饱和度由如下公式(5)或(6)提取和检测;
SpO2=α-βR (5)
SpO2=α-βR-γR2 (6)
其中, 分别为红光和红外光强度的交流分量,分别为红光和红外光强度的直流分量;α、β以及γ分别为经验值,通过定标来确定。
6.如权利要求1所述的单一光电传感器睡眠呼吸多生理参数监测方法,其特征在于,所述步骤5中,呼吸的提取和检测采用数字滤波法,其输入输出关系如下公式(8)所示;
y ( n ) = { Σ k = 0 N a k x ( n - k ) } - { Σ k = 1 M b k y ( n - k ) } - - - ( 8 )
其中,ak和bk为IIR数字滤波器的系数,N、M分别为数字滤波器系统传递函数的分子和分母多项式系数阶数,n为自然数;光电容积脉搏波通过数字滤波得到呼吸波波形,然后通过过零法或峰值法检测计算得到呼吸率。
7.如权利要求1所述的单一光电传感器睡眠呼吸多生理参数监测方法,其特征在于,所述步骤5中,脉率变异性的提取和检测采用脉率变异性的非线性分析方法,计算脉率间期彭家莱图的标准差SD1、SD2以及二者的比值SD1/SD2,具体如下公式(9)和(10);
S D 1 2 = 1 N - 1 Σ i = 1 N - 1 ( D i ′ ) 2 - - - ( 9 )
S D 2 2 = 1 N - 1 Σ i = 1 N - 1 ( D i ′ ′ ) 2 - - - ( 10 )
其中,N为RR间期总数,RRn和RRn+1分别为第n与第n+1个脉率间期,为全部脉率间期序列的间期均值,n为自然数。
8.单一光电传感器睡眠呼吸多生理参数监测装置,其特征在于,包括下述模块:
原始数据采集模块,用于以模拟方式或数字方式采集采样率在500Hz以上的光电容积脉搏波;
预处理模块,用于对原始数据采集模块采集到的光电容积脉搏波进行预处理,包括基线的提取检测、峰值点与谷值点的提取与检测;
运动识别模块,用于根据预处理模块检测的峰值点、谷值点以及基线变化信息判断当前测量时刻是否处于运动状态;当峰值点、谷值点以及基线变化幅度超过所设定的阈值时认为当前时刻处于运动状态,将丢弃当前测量值重新测量,直至当前测量值脱离运动状态进入正常光电容积脉搏波范围内;所检测到的运动信息送至体动检测模块进行体动次数计数;
正常光电容积脉搏波模块,用于将运动识别模块识别的非运动状态数据视为正常光电容积脉搏波,用于后续信号处理;
脉率检测模块、血氧饱和度检测模块、呼吸检测模块和脉率变异性检测模块,分别用于根据正常光电容积脉搏波模块提取的正常光电容积脉搏波,检测脉率信息、血氧饱和度信息、呼吸信息和脉率变异性信息,并输送至数据后处理模块;
体动检测模块,用于根据运动识别模块提取的运动信息,将所检测到的运动信息在当前运动停止后记为体动增加1次,累加至体动次数计数值,并输送至数据后处理模块;
数据后处理模块,对脉率检测模块、血氧饱和度检测模块、呼吸检测模块、脉率变异性检测模块和体动检测模块提取的脉率、血氧饱和度、呼吸、脉率变异性以及体动次数信息,与上位机进行数据传输、显示和记录。
9.如权利要求8所述的单一光电传感器睡眠呼吸多生理参数监测装置,其特征在于,所述原始数据采集模块包括发光驱动、红光光二极管、红外光发光二极管和光频转换器;外部微处理器控制的发光驱动用于交替产生驱动电流,分别驱动红光发光二极管以及红外光发光二极管;光频转换器用于将所发出的红光、红外光以反射或者透射方式穿过被测组织的红光和红外光的光强转化为频率信号,输出到微处理器捕捉计数。
10.如权利要求8所述的单一光电传感器睡眠呼吸多生理参数监测装置,其特征在于,所述数据后处理模块包括无线蓝牙数据传输模块、MicroSD卡数据存储模块和显示模块。
CN201610893587.2A 2016-10-13 2016-10-13 单一光电传感器睡眠呼吸多生理参数监测方法及装置 Active CN106264475B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610893587.2A CN106264475B (zh) 2016-10-13 2016-10-13 单一光电传感器睡眠呼吸多生理参数监测方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610893587.2A CN106264475B (zh) 2016-10-13 2016-10-13 单一光电传感器睡眠呼吸多生理参数监测方法及装置

Publications (2)

Publication Number Publication Date
CN106264475A true CN106264475A (zh) 2017-01-04
CN106264475B CN106264475B (zh) 2019-11-08

Family

ID=57717323

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610893587.2A Active CN106264475B (zh) 2016-10-13 2016-10-13 单一光电传感器睡眠呼吸多生理参数监测方法及装置

Country Status (1)

Country Link
CN (1) CN106264475B (zh)

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107874750A (zh) * 2017-11-28 2018-04-06 华南理工大学 脉率变异性和睡眠质量融合的心理压力监测方法及装置
CN108464836A (zh) * 2018-02-09 2018-08-31 重庆东渝中能实业有限公司 一种面向社区医疗的血氧饱和度检测系统及方法
WO2019000230A1 (zh) * 2017-06-27 2019-01-03 华为技术有限公司 一种睡眠监测方法、装置和可穿戴设备
CN109199355A (zh) * 2018-09-18 2019-01-15 深圳和而泰数据资源与云技术有限公司 心率信息检测方法、装置和检测设备
CN109528214A (zh) * 2018-11-07 2019-03-29 深圳市新元素医疗技术开发有限公司 一种多功能腕式血氧仪
CN110192846A (zh) * 2019-05-31 2019-09-03 出门问问信息科技有限公司 可穿戴设备
TWI674882B (zh) * 2018-10-05 2019-10-21 國立成功大學 一種呼吸速率連續評估方法
CN110463215A (zh) * 2017-01-31 2019-11-15 西门子股份公司 用于提供技术设备的测量值的方法、技术系统和用于运行技术系统的方法
CN111358454A (zh) * 2020-03-17 2020-07-03 乐普(北京)医疗器械股份有限公司 一种基于差异波占比检测信号质量的方法和装置
CN112005204A (zh) * 2018-04-27 2020-11-27 宁波吉利汽车研究开发有限公司 多媒体效果
CN114176527A (zh) * 2021-12-29 2022-03-15 深圳融昕医疗科技有限公司 体动状态的确定方法及计算机可读存储介质
CN114257683A (zh) * 2020-09-21 2022-03-29 Oppo广东移动通信有限公司 消息提示方法及装置、可穿戴设备、计算机可读存储介质
CN115137316A (zh) * 2022-06-10 2022-10-04 佳禾智能科技股份有限公司 一种多生理参数监测手表及其监测方法
CN117426758A (zh) * 2023-12-20 2024-01-23 武汉纺织大学 一种基于多传感信息融合的智慧服装系统及方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101732050A (zh) * 2009-12-04 2010-06-16 西安交通大学 一种基于光电容积波的呼吸率监测方法
CN102028457A (zh) * 2010-11-24 2011-04-27 北京麦邦光电仪器有限公司 脉率测量方法及指环式脉率测量仪
CN102631199A (zh) * 2012-04-20 2012-08-15 天津大学 一种快速的动态光谱数据提取方法
JP2014168710A (ja) * 2009-02-02 2014-09-18 Fukuda Denshi Co Ltd 在宅酸素療法管理装置、生体情報測定装置及び動作情報取得装置
CN104622444A (zh) * 2015-01-30 2015-05-20 中国科学院电子学研究所 一种多光电传感模块腕式监测系统

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2014168710A (ja) * 2009-02-02 2014-09-18 Fukuda Denshi Co Ltd 在宅酸素療法管理装置、生体情報測定装置及び動作情報取得装置
CN101732050A (zh) * 2009-12-04 2010-06-16 西安交通大学 一种基于光电容积波的呼吸率监测方法
CN102028457A (zh) * 2010-11-24 2011-04-27 北京麦邦光电仪器有限公司 脉率测量方法及指环式脉率测量仪
CN102631199A (zh) * 2012-04-20 2012-08-15 天津大学 一种快速的动态光谱数据提取方法
CN104622444A (zh) * 2015-01-30 2015-05-20 中国科学院电子学研究所 一种多光电传感模块腕式监测系统

Cited By (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11131606B2 (en) 2017-01-31 2021-09-28 Siemens Energy Global GmbH & Co. KG Method for providing measured values of a technical installation, technical system, and method for operating the technical system
CN110463215A (zh) * 2017-01-31 2019-11-15 西门子股份公司 用于提供技术设备的测量值的方法、技术系统和用于运行技术系统的方法
WO2019000230A1 (zh) * 2017-06-27 2019-01-03 华为技术有限公司 一种睡眠监测方法、装置和可穿戴设备
CN107874750A (zh) * 2017-11-28 2018-04-06 华南理工大学 脉率变异性和睡眠质量融合的心理压力监测方法及装置
CN107874750B (zh) * 2017-11-28 2020-01-10 华南理工大学 脉率变异性和睡眠质量融合的心理压力监测方法及装置
CN108464836A (zh) * 2018-02-09 2018-08-31 重庆东渝中能实业有限公司 一种面向社区医疗的血氧饱和度检测系统及方法
CN108464836B (zh) * 2018-02-09 2021-09-03 重庆东渝中能实业有限公司 一种面向社区医疗的血氧饱和度检测系统及方法
CN112005204A (zh) * 2018-04-27 2020-11-27 宁波吉利汽车研究开发有限公司 多媒体效果
CN109199355A (zh) * 2018-09-18 2019-01-15 深圳和而泰数据资源与云技术有限公司 心率信息检测方法、装置和检测设备
CN109199355B (zh) * 2018-09-18 2021-09-28 深圳和而泰数据资源与云技术有限公司 心率信息检测方法、装置和检测设备
TWI674882B (zh) * 2018-10-05 2019-10-21 國立成功大學 一種呼吸速率連續評估方法
CN109528214A (zh) * 2018-11-07 2019-03-29 深圳市新元素医疗技术开发有限公司 一种多功能腕式血氧仪
CN110192846A (zh) * 2019-05-31 2019-09-03 出门问问信息科技有限公司 可穿戴设备
CN111358454A (zh) * 2020-03-17 2020-07-03 乐普(北京)医疗器械股份有限公司 一种基于差异波占比检测信号质量的方法和装置
CN111358454B (zh) * 2020-03-17 2022-07-29 乐普(北京)医疗器械股份有限公司 一种基于差异波占比检测信号质量的方法和装置
CN114257683A (zh) * 2020-09-21 2022-03-29 Oppo广东移动通信有限公司 消息提示方法及装置、可穿戴设备、计算机可读存储介质
CN114176527A (zh) * 2021-12-29 2022-03-15 深圳融昕医疗科技有限公司 体动状态的确定方法及计算机可读存储介质
CN115137316A (zh) * 2022-06-10 2022-10-04 佳禾智能科技股份有限公司 一种多生理参数监测手表及其监测方法
CN117426758A (zh) * 2023-12-20 2024-01-23 武汉纺织大学 一种基于多传感信息融合的智慧服装系统及方法
CN117426758B (zh) * 2023-12-20 2024-04-05 武汉纺织大学 一种基于多传感信息融合的智慧服装系统及方法

Also Published As

Publication number Publication date
CN106264475B (zh) 2019-11-08

Similar Documents

Publication Publication Date Title
CN106264475A (zh) 单一光电传感器睡眠呼吸多生理参数监测方法及装置
US11229374B2 (en) Plethysmograph variability processor
US6529752B2 (en) Sleep disorder breathing event counter
JP7019611B2 (ja) 被検者の呼吸情報を決定するための方法及び装置
EP3454730B1 (en) System for determining the frequency of a periodic physiological process of a subject
US20090024044A1 (en) Data recording for patient status analysis
US8197415B2 (en) Non-invasive heart monitoring apparatus and method
CN104545888B (zh) 一种基于动态心电与呼吸波采集的睡眠呼吸暂停采集分析系统
CN109512409B (zh) 一种心血管健康评估装置及其使用方法
CN104622444B (zh) 一种多光电传感模块腕式监测系统
Murali et al. A wearable device for physical and emotional health monitoring
EP3500154B1 (en) A sensor system and method for determining a breathing type
CN204971245U (zh) 一种人体生命体征多参数监控系统
TWM542444U (zh) 穿戴式多功能脈搏波診斷與分析裝置
CN107174225A (zh) 一种臂带式血管生理参数监测装置
CN209770350U (zh) 一种心血管健康评估装置
CN1073392C (zh) 多参数监护仪及其检测方法
US20220287592A1 (en) Behavior task evaluation system and behavior task evaluation method
Leier et al. Respiration signal extraction from photoplethysmogram using pulse wave amplitude variation
Swetha et al. Analytical Study of Functional Trends in Cardio Respiratory Activity
de Gevigney et al. Advanced Photodetector for Athlete Monitoring
CN109758124A (zh) 一种手功能康复的客观评价方法及装置和系统
Jamwal et al. Pulse Signal Processing for Analysis of Cardiovascular Diseases
El Barike Belhcen Development of a non-invasive device for assessing diabetic and non-diabetic individuals based on heart rate analysis
Yang Physiological indexes collection system for multifunctional gyms based on FFT and photoelectric sensor

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20231110

Address after: West side of the second floor of Building 3, No. 8-3 Chunwang Road, Panyang Industrial Park, Huangdai Town, Xiangcheng District, Suzhou City, Jiangsu Province, 215000

Patentee after: Suzhou Sanlang Medical Technology Co.,Ltd.

Address before: 710000 workshop 2 in yard 5, West Qinling Third Road, Caotang science and technology industry base, high tech Zone, Xi'an, Shaanxi Province

Patentee before: SHAANXI HUACHEN PETROLEUM TECHNOLOGY CO.,LTD.

Effective date of registration: 20231110

Address after: 710000 workshop 2 in yard 5, West Qinling Third Road, Caotang science and technology industry base, high tech Zone, Xi'an, Shaanxi Province

Patentee after: SHAANXI HUACHEN PETROLEUM TECHNOLOGY CO.,LTD.

Address before: Beilin District Xianning West Road 710049, Shaanxi city of Xi'an province No. 28

Patentee before: XI'AN JIAOTONG University

TR01 Transfer of patent right