CN103961080B - 生物体信息检测装置以及生物体信息检测方法 - Google Patents
生物体信息检测装置以及生物体信息检测方法 Download PDFInfo
- Publication number
- CN103961080B CN103961080B CN201410044180.3A CN201410044180A CN103961080B CN 103961080 B CN103961080 B CN 103961080B CN 201410044180 A CN201410044180 A CN 201410044180A CN 103961080 B CN103961080 B CN 103961080B
- Authority
- CN
- China
- Prior art keywords
- signal
- parameter
- observation
- acceleration
- value
- 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
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 28
- 230000001133 acceleration Effects 0.000 claims abstract description 163
- 230000009471 action Effects 0.000 claims abstract description 128
- 210000003462 vein Anatomy 0.000 claims abstract description 49
- 210000001367 artery Anatomy 0.000 claims abstract description 48
- JEIPFZHSYJVQDO-UHFFFAOYSA-N ferric oxide Chemical compound O=[Fe]O[Fe]=O JEIPFZHSYJVQDO-UHFFFAOYSA-N 0.000 claims abstract description 36
- 238000004458 analytical method Methods 0.000 claims abstract description 13
- 239000000203 mixture Substances 0.000 claims abstract description 13
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 8
- 238000012360 testing method Methods 0.000 claims abstract description 8
- 238000003786 synthesis reaction Methods 0.000 claims abstract description 7
- 230000003068 static effect Effects 0.000 claims description 55
- 230000017531 blood circulation Effects 0.000 claims description 12
- 230000000284 resting effect Effects 0.000 claims description 11
- 238000000034 method Methods 0.000 description 77
- 238000012545 processing Methods 0.000 description 50
- 230000008569 process Effects 0.000 description 48
- 238000009532 heart rate measurement Methods 0.000 description 41
- 210000000707 wrist Anatomy 0.000 description 24
- 230000015654 memory Effects 0.000 description 20
- 238000010586 diagram Methods 0.000 description 11
- 238000005259 measurement Methods 0.000 description 9
- 230000008859 change Effects 0.000 description 8
- 238000010606 normalization Methods 0.000 description 8
- 239000008280 blood Substances 0.000 description 7
- 210000004369 blood Anatomy 0.000 description 7
- 210000004204 blood vessel Anatomy 0.000 description 7
- 230000005622 photoelectricity Effects 0.000 description 6
- 230000001932 seasonal effect Effects 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 230000036541 health Effects 0.000 description 4
- 238000005286 illumination Methods 0.000 description 4
- 238000004321 preservation Methods 0.000 description 4
- 210000000709 aorta Anatomy 0.000 description 3
- 238000013500 data storage Methods 0.000 description 3
- 238000009434 installation Methods 0.000 description 3
- 230000005012 migration Effects 0.000 description 3
- 238000013508 migration Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000004020 luminiscence type Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 102000001554 Hemoglobins Human genes 0.000 description 1
- 108010054147 Hemoglobins Proteins 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 210000003423 ankle Anatomy 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000003862 health status Effects 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000007493 shaping process Methods 0.000 description 1
- 238000007920 subcutaneous administration Methods 0.000 description 1
- 238000002834 transmittance Methods 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, 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/024—Detecting, measuring or recording pulse rate or heart rate
- A61B5/02438—Detecting, measuring or recording pulse rate or heart rate with portable devices, e.g. worn by the patient
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, 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/024—Detecting, measuring or recording pulse rate or heart rate
- A61B5/02416—Detecting, measuring or recording pulse rate or heart rate using photoplethysmograph signals, e.g. generated by infrared radiation
- A61B5/02427—Details of sensor
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/68—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
- A61B5/6801—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
- A61B5/6802—Sensor mounted on worn items
- A61B5/681—Wristwatch-type devices
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
- A61B5/7207—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise induced by motion artifacts
- A61B5/721—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise induced by motion artifacts using a separate sensor to detect motion or using motion information derived from signals other than the physiological signal to be measured
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2562/00—Details of sensors; Constructional details of sensor housings or probes; Accessories for sensors
- A61B2562/04—Arrangements of multiple sensors of the same type
- A61B2562/046—Arrangements of multiple sensors of the same type in a matrix array
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Cardiology (AREA)
- Public Health (AREA)
- Surgery (AREA)
- Veterinary Medicine (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Animal Behavior & Ethology (AREA)
- Physiology (AREA)
- Signal Processing (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Psychiatry (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
Abstract
本发明提供一种生物体信息检测装置以及生物体信息检测方法。检测利用者的脉拍的生物体信息检测装置,具备:检测部,其输出根据所述利用者的至少1个观测部位的脉波而检测出的观测信号;加速度测量部,其输出伴随所述利用者的动作而测量的多个轴方向的多个加速度信号。并且,基于根据多个参数合成所述多个加速度信号而得的合成加速度信号与所述观测信号的比较,推定在所述观测信号中包含的、与基于所述利用者的动作的加速度成分对应的所述各参数的特定的值,并且根据从所述观测信号中去除与所推定的所述各参数的所述特定的值对应的特定的合成加速度信号所得的差分,计算所述利用者的脉拍数。由此,可以比较准确地测量动作中的利用者的脉拍数。
Description
本申请要求了申请日为2013年2月6日、申请号为特愿2013-021247的日本专利申请的优先权,其内容通过引用被合并于此。
技术领域
本发明涉及生物体信息检测装置以及生物体信息检测方法、生物体信息检测程序,特别涉及具备在运动时安装在人体上来测定脉拍的脉拍测定功能的生物体信息检测装置以及生物体信息检测方法。
背景技术
近年来,由于健康意向的提高,日常进行跑步或散步、骑车等运动来维持、促进健康状态的人不断增加。这样的人群中,为了掌握自己的健康状态或运动状态,测量或记录各种生物体信息。
作为用于掌握人体状态的生物体信息,有各种生理指标。作为该生理指标之一,例如公知1分钟的心脏的跳动数、即心拍数。
作为心拍数的测量方法,一般已知心电图方式。在该心电图方式中需要在胸部安装多个电极。在安装了该电极的状态下,有时会限制日常生活或运动时的行动,电极的安装也复杂,因此有时给测定设备的使用者(用户)带来很大负担。
因此,如今作为可以更简便地进行测量的生理指标,代替心拍而测量脉拍的方法被广泛使用。
作为脉拍的测量方法,例如已知光电脉波法(或光学式脉波检测法)。该光电脉波法的原理大致是利用血液中的血红蛋白的光吸收特性,检测与脉波对应的观测信号。即,脉波是由于心拍而引起的动脉内的压力变化作为波动而传递到末梢动脉而得。并且,通过使红外线等的光透过皮肤照射到末梢动脉中的血液,测量由该血液散射的光的反射光的强度的时间变化来作为观测信号,可以检测表示末梢动脉的血流的波动的流量变化的脉波。根据这样的光电脉波法,可以从手指或耳朵、手腕等取得脉波,并据此简单地求出脉拍。
但是,血流伴随日常生活时或运动时的身体的动作(体动)而变化。因此,在光电脉波法中具有在很大程度上受到该体动引起的血流变化(体动噪音)的影响,在光电脉波法的观测信号中混入该体动噪音的问题。
与此相对,作为从混合了脉波信号和体动噪音的观测信号中去除体动噪音的信号成分,得到脉波信号的方法之一,例如在日本特开2003-102694号公报中记载了将体动噪音视为通过加速度计取得的加速度信号,将观测信号和加速度信号的差分信号设为脉波信号的方法。
在上述那样的取得脉波信号的方法中,由于将加速度信号视为与体动噪音等同来进行信号处理,因此可以通过比较简单的信号处理取得脉波信号。
但是,根据本申请的发明人的验证,在人体的动作中取得的加速度信号和体动噪音未必相同,例如已知加速度信号和体动噪音振幅不同,或者从加速度产生的时刻到其影响显现在观测信号中为止存在时间差(时滞)。而且,该时间差不固定,已知时间差根据人体的动作状态的变化或脉拍的测量位置的变化而变化。
在上述那样的方法中未完全考虑到这样的问题。因此,在上述那样的方法中无法恰当地去除在观测信号中包含的由体动导致的噪音成分,无法准确地测量人体的动作中的脉拍数。
发明内容
本发明具有能够提供可以恰当地推定根据利用者的脉波检测出的观测信号中的、伴随利用者的动作的体动噪音成分,检测出准确的脉波的生物体信息检测装置以及生物体信息检测方法、生物体信息检测程序的优点。
本发明提供一种生物体信息检测装置,其中,具备:
检测部,其输出根据利用者的至少1个观测部位的脉波而检测出的观测信号;
加速度测量部,其输出伴随所述利用者的动作而测量的、与互不相同的多个轴方向的各个轴方向对应的多个加速度信号;
参数推定部,其基于根据多个参数合成所述多个加速度信号而得的合成加速度信号与所述观测信号的比较,推定在所述观测信号中包含的、与基于所述利用者的动作的加速度成分对应的所述各参数的特定的值;以及
脉拍数计算部,其根据从所述观测信号中去除与所述各参数的所述特定的值对应的特定的合成加速度信号所得的差分,计算所述利用者的脉拍数。
另外,本发明还提供一种生物体信息检测方法,其中,
取得基于利用者的至少1个观测部位的脉波的观测信号,并且取得伴随该利用者的动作的、与互不相同的多个轴方向的各个轴方向对应的多个加速度信号,
基于根据多个参数合成所述多个加速度信号而得的合成加速度信号与所述观测信号的比较,推定在所述观测信号中包含的、与基于所述利用者的动作的加速度成分对应的所述各参数的特定的值,
根据从所述观测信号中去除与所述各参数的所述特定的值对应的特定的合成加速度信号所得的差分,计算所述利用者的脉拍数。
附图说明
图1A、图1B是表示本发明的生物体信息检测装置的安装例以及外观结构例的概要图。
图2A、图2B、图2C是表示本发明的生物体信息检测装置的测量面的结构例的概要图。
图3是表示第1实施方式的生物体信息检测装置的一个结构例的框图。
图4是表示在第1实施方式的生物体信息检测装置的生物体信息检测方法中执行的静止时脉波测量动作的流程图。
图5A、图5B、图5C、图5D是表示在第1实施方式的静止时脉波测量动作中执行的脉波的多点观测的一例的概念图。
图6是表示通过第1实施方式的静止时脉波测量动作取得的脉波信号的一例的波形图。
图7是表示在第1实施方式的生物体信息检测装置的生物体信息检测方法中执行的动作时脉波测量动作的流程图。
图8是用于说明在第1实施方式的动作时脉波测量动作中计算的极值间隔的概念图。
图9A、图9B、图9C是表示通过第1实施方式的动作时脉波测量动作取得的各信号的一例的波形图。
图10是表示在第1实施方式的动作时脉波测量动作中执行的时滞、旋转角度推定处理的流程图。
图11是用于说明在第1实施方式的动作时脉波测量动作中定义的3轴方向的概念图。
图12是表示通过第1实施方式的时滞、旋转角度推定处理计算的归一化相互相关系数的一例的图。
图13是表示通过第1实施方式的时滞、旋转角度推定处理取得的旋转角度和极大值的迁移的一例的图。
图14是表示在第1实施方式的动作时脉波测量动作中执行的振幅推定处理的流程图。
图15A、图15B是表示第2实施方式的生物体信息检测装置的测量面的结构例的概要图。
图16是表示在第2实施方式的生物体信息检测方法中执行的动作时脉波测量动作的流程图。
具体实施方式
以下,关于本发明的生物体信息检测装置以及生物体信息检测方法、生物体信息检测程序,表示出实施方式来详细说明。
<第1实施方式>
(生物体信息检测装置)
图1A、图1B是表示本发明的生物体信息检测装置的安装例以及外观结构的概要图。
在此,图1A是表示将本发明的生物体信息检测装置安装在人体上的状态的概要图,图1B是表示本发明的生物体信息检测装置的正面以及侧面的概要结构图。
图2A、图2B、图2C是表示本发明的生物体信息检测装置的测量区域的结构例的概要图。
本发明的生物体信息检测装置100例如如图1A所示,具有安装在用户(利用者)US的手腕等上的腕表型(或腕带型)的外观形状。
生物体信息检测装置100例如如图1B所示,大体上具备:具有测量用户US的脉拍并且向用户US提供预定的信息的功能的设备本体101;通过缠绕在用户US的手腕USh上,用于将上述设备本体101安装在手腕USh上并使其贴紧的带部102。
在设备本体101的与手腕USh接触的一侧的面(图1B的右图的、沿II-II线方向视图的右侧面)的预定区域中设置了测量区域MS。
在该测量区域MS中,例如如图2A~2C所示,以预定的图案二维排列了1至多个发光元件E1~E9和1至多个受光元件R1~R4。
在测量区域MS中,发光元件和受光元件例如如图2A、2B、2C所示那样配置。
在图2A所示的配置中,在1个发光元件E1的周围以包围方式配置多个(4个)受光元件R1~R4。即,发光元件和受光元件以1对多的关系排列。
在图2B所示的配置中,在1个受光元件R1的周围以包围方式配置多个(4个)发光元件E1~E4。即,发光元件和受光元件以多对1的关系排列。
在图2C所示的配置中,在多个(4个)受光元件R1~R4的各自的周围以包围方式配置多个发光元件E1~E9。即,发光元件和受光元件以多对多的关系排列。
这样,在本实施方式中,具有将1至多个发光元件以及1至多个受光元件中的至少某一方配置了多个的结构。
此外,测量区域MS中排列的发光元件和受光元件的个数或配置方式不限于图2A~图2C所示的图案。也可以将任意个数的发光元件或受光元件按照交错状、格子状、圆弧状等任意图案交替排列。
图3是表示本实施方式的生物体信息检测装置的一个结构例的框图。
生物体信息检测装置100具体来说,例如如图3所示,大体上具备:发光部(检测部)10、发光控制部15、受光部(检测部)20、加速度测量部30、信号放大部40、滤波器部50、存储器部60、静止时脉波振幅记录部(存储部)65、信号处理部(观测信号选择部、参数推定部、脉拍数计算部)70、显示部80和操作部90。
发光部10具有上述1至多个发光元件E1~E9,如图2A~图2C所示,在设备本体101的与手腕USh接触的一侧的面的测量区域MS中以预定的图案排列。
发光元件E1~E9例如可以应用发光二极管(LED:LightEmittingDiode)等。发光元件E1~E9按照后述的发光控制部15的驱动控制,以预定的发光强度(或发光量)发出可见光,并将发出的光照射到手腕USh的皮肤表面(体表面)SF。
在此,在使用可见光的反射式的脉波检测法中,由于可见光在体内的透射性低,因此具有难以受到来自存在于体内深部的静脉或动脉的血流的反射光的影响,难以受到在各个血管中发生的血流路长度导致的拍动的传递时滞的影响的优点。
此外,作为从发光元件发出的可见光,例如可以良好地应用波长525nm左右的绿色可见光。
发光控制部15按照来自后述的信号处理部70的控制,使构成发光部10的1至多个发光元件E1~E9以预定的点亮模式(即预定的顺序且预定的发光强度)分别发光。
受光部20具有上述的1至多个受光元件R1~R4,如图2A~图2C所示,在设备本体101的测量区域MS中以预定的图案排列。
受光元件R1~R4例如可以应用光电晶体管或照度传感器等。受光元件R1~R4接受从上述的1至多个发光元件E1~E9分别发光并照射皮肤表面SF的观测脉波的观测部位Pm后,通过观测部位Pm附近的血管中的血液而散射的光作为反射光,由此输出与受光量对应的输出信号(观测信号)。
加速度测量部30具有3轴加速度传感器。3轴加速度传感器将在用户US的动作中施加给生物体信息检测装置100的移动速度的变化的比例(加速度)作为加速度信号而输出。
从该加速度测量部30输出的加速度信号如后所述,作为与由x轴、y轴、z轴构成的相互垂直的3轴方向的各自对应的3个加速度信号而被输出。
信号放大部40将通过受光部20取得的观测信号以及通过加速度测量部30测量的加速度信号放大到适合于后述的信号处理部70中的信号处理的预定的信号水平。
滤波器部50使通过信号放大部40放大后的上述观测信号以及加速度信号中的预定频带的信号成分通过,提供给信号处理部70。
存储器部60例如具有数据保存用存储器(以下记作“数据存储器”)、程序保存用存储器(以下记作“程序存储器”)、作业数据保存用存储器(以下记作“作业用存储器”)。
数据存储器具有闪速存储器等非易失性存储器,在用户US的动作时或运动时,将上述通过受光部20取得的观测信号、通过加速度测量部30测量的加速度信号与时间数据关联起来保存(记录)在预定的存储区域。
程序存储器具有ROM(只读存储器),保存用于实现生物体信息检测装置100的各结构(发光部10或受光部20、加速度测量部30、后述的显示部80或操作部90等)的预定的功能的控制程序、用于实现基于上述观测信号或加速度信号计算脉拍数的功能的算法程序。
作业用存储器具有RAM(随机存取存储器),临时保存在执行上述控制程序以及算法程序时使用或生成的各种数据。
此外,数据存储器的一部分或全部例如可以具有作为存储卡等可移动存储介质的形态,以能够相对于生物体信息检测装置100的设备本体101装卸的方式构成。
静止时脉波振幅记录部65,在用户US不动作的静止时或安静时,将上述通过受光部20取得的观测信号的信号波(脉波)的振幅与时间数据关联起来,保存(记录)在预定的存储区域中。
信号处理部70是CPU(中央运算装置)或MPU(微处理器单元),按照在上述存储器60中保存的控制程序进行处理。由此,信号处理部70控制存储器部60中的各种数据的保存或读出动作、显示部80中的各种信息的显示动作、操作部90中的输入操作的检测动作等。
信号处理部70按照在上述存储器部60中保存的算法程序进行处理,由此如后述的生物体信息检测方法所示,执行根据通过受光部20取得的观测信号、通过加速度测量部30测量的加速度信号计算脉拍数的动作等。
此外,信号处理部70中执行的控制程序或算法程序可以预先被安装在信号处理部70的内部。
显示部80例如具有能够彩色或单色显示的液晶显示面板或有机EL显示面板等显示装置,至少显示通过信号处理部70计算出的脉拍数。
此外,显示部80可以除了脉拍数以外或代替脉拍数,用文字或数字信息、图像信息等显示脉波(脉的波形数据)或移动速度、步数、当前时刻等。
在此,例如在脉拍的波形数据(脉波数据)中包含与血流关联的各种信息。
即,可以将脉拍数据作为用于判定例如健康或身体状况(血管的堵塞或血管年龄、紧张状态的判定等)、运动状态等的重要参数来应用,将与它们对应的判定结果用特定的文字或数字信息、图像信息、发光图案等显示在显示部80中。
此外,在本实施方式中,作为向用户US提供或通知各种信息的输出接口,仅表示出显示部80,但是不限于此。例如,在显示部80以外,可以具备发出特定的音色或声音消息的蜂鸣器或扬声器等音响部、以特定的振动模式振动的振动部等其他接口。
操作部90具有按钮开关或滑动开关、键盘、在显示部80的前面配置或一体化形成的触摸面板等。操作部90用于进行生物体信息检测装置100中的电源的接通、断开动作、脉波或加速度的测量动作、显示部80中的显示动作等各种动作的选择或执行、设定值等的输入操作。
(生物体信息检测方法)
接着,说明上述生物体信息检测装置中的生物体信息检测方法。
具有上述那样的结构的生物体信息检测装置中的生物体信息检测方法,概要地执行用于取得静止时的脉波的观测信号的静止时脉波测量动作、根据动作时取得的脉波的观测信号和加速度信号计算脉拍数的动作时脉波测量动作。
(静止时脉波测量动作)
图4是表示在本实施方式的生物体信息检测装置的生物体信息检测方法中执行的静止时脉波测量动作的流程图。
图5A、图5B、图5C、图5D是表示本实施方式的静止时脉波测量动作中执行的脉波的多点观测的一例的概念图。
图6是表示通过本实施方式的静止时脉波测量动作取得的脉波信号的一例的波形图。
在静止时脉波测量动作中,如图4所示,首先以一定时间取得用户US不进行运动等动作的静止状态或安静状态下的脉波的观测信号和加速度信号(步骤S101)。
具体来说,信号处理部70在显示部80中显示测量静止时的脉波的意思的文字信息或图像信息等,催促用户US保持静止状态或安静状态。
接着,信号处理部70指定发光部10的特定的发光元件和受光部20的特定的受光元件的组合,通过发光控制部15使指定的发光元件以预定的发光强度发光。由此,在用户US的皮肤表面SF的观测脉波的区域(观测部位Pm)照射从发光元件出射的光。
所照射的光的一部分由观测部位Pm附近的血管的血液散射,从皮肤表面SF作为发射光而出射。
该反射光由上述指定的受光元件接受。然后,与受光元件的受光量对应的输出信号经由信号放大部40以及滤波器部50作为观测信号被输出到信号处理部70。
在此,以发光元件和受光元件的排列为图2C所示的发光元件E1~E9和受光元件R1~R4的排列图案时为例,详细说明基于特定的发光元件和特定的受光元件的组合的观测信号的取得动作的一例。
首先,信号处理部70,例如如图5A所示,指定发光元件E1和受光元件R1、发光元件E3和受光元件R2、发光元件E7和受光元件R3、发光元件E9和受光元件R4的组合。
接着,通过发光控制部15使各发光元件E1、E3、E7、E9以预定的发光强度发光,使光照射皮肤表面SF的各观测部位Pm11、Pm32、Pm73、Pm94,由各受光元件R1、R2、R3、R4接受其反射光。
由此,取得皮肤表面SF的各观测部位Pm11、Pm32、Pm73、Pm94的静止时的脉波的观测信号。
在此,各观测部位Pm11、Pm32、Pm73、Pm94的观测信号的取得动作,例如按照观测部位Pm11、Pm32、Pm73、Pm94的顺序按时序执行。此外,观测信号的取得动作也可以在各观测部位Pm11、Pm32、Pm73、Pm94同时并行地执行。
然后,信号处理部70例如如图5B所示,指定发光元件E5和各受光元件R1~R4的组合。
接着,通过发光控制部15使发光元件E5以预定的发光强度发光,使光照射皮肤表面SF的各观测部位Pm51、Pm52、Pm53、Pm54,由各受光元件R1~R4接受其反射光。
由此,取得皮肤表面SF的各观测部位Pm51~Pm54的静止时的脉波的观测信号。
在此,与图5A所示的情况相同,各观测部位Pm51~Pm54的观测信号的取得动作分别针对观测部位Pm51~Pm54按时序来执行。此外,也可以在各观测部位Pm51~Pm54同时并行执行。
以下,同样地通过信号处理部70例如如图5C、图5D所示,分别指定发光元件E2和各受光元件R1、R2的组合、发光元件E8和各受光元件R3、R4的组合、以及发光元件E4和各受光元件R1、R3的组合、发光元件E6和各受光元件R2、R4的组合,取得各观测部位Pm21、Pm22、Pm83、Pm84以及Pm41、Pm62、Pm43、Pm64的静止时的脉波的观测信号。
通过这样的一系列动作(多点观测),取得在测量区域MS内排列的、相邻的发光元件和受光元件之间的各观测部位的脉波的观测信号。
另外,这样的脉波的观测信号的取得动作,以包含几个至十几个左右表示脉波的波形的任意时间、例如几秒至10秒左右的时间持续执行。
另一方面,与脉波的观测信号的取得动作并行地,信号处理部70控制加速度测量部30来测量用户US的3轴方向的加速度。
在此,在上述脉波的观测信号的取得动作的期间中持续执行3轴加速度的测量动作。
通过加速度测量部30测量到的3轴加速度经由信号放大部40以及滤波器部50作为加速度信号被输出到信号处理部70。
如此取得的脉波的观测信号和加速度信号,在根据时间数据相互关联后被保存在存储器部60的预定的存储区域。
接着,判断在上述步骤S101中取得的3轴方向的各加速度信号的振幅是否在预定的阈值以下(步骤S102)。
具体来说,信号处理部70从存储器部60读出在脉波的观测信号的取得动作中取得的加速度信号,针对该3轴方向的各加速度信号判断信号波形的极大值与极小值的差分、即振幅的最大值(最大振幅)是否在与用户US不进行运动等动作的静止状态或安静状态时对应的预定的阈值以下。即,根据该阈值判断用户US是否处于静止状态或安静状态。
在此,用于判断用户US的静止状态或安静状态的阈值,可以设定为跑步等动作时的振幅的例如5%左右。
该阈值例如可以根据过去的用户US的动作中的加速度信号来设定,也可以根据从不特定多数的样本取得的一般动作中的加速度信号来设定,也可以在用户US为静止状态或安静状态时任意地设定。
此外,如后所述,本申请的发明人发现z轴方向的加速度信号对脉波信号几乎没有影响。因此,在上述各加速度信号的振幅是否在预定阈值以下的判断中,可以不对z轴方向的加速度信号的振幅进行判断。
在上述步骤S102中判定为所取得的3轴方向的各加速度信号的振幅在阈值以下时,判断出用户US处于静止状态或安静状态,将此时的各观测部位的脉波的观测信号的振幅的平均值作为静止时的观测信号的振幅来记录(步骤S103)。
具体来说,信号处理部70在判定为所取得的3轴方向的各加速度信号的振幅都在上述阈值以下时,读出根据时间数据与该加速度信号关联起来保存在存储器部60中的、各观测部位的脉波的观测信号,计算这些信号波形的极大值与极小值的差分、即振幅的平均值。
然后,信号处理部70将计算出的平均值作为静止时的观测信号的振幅,保存(记录)在静止时脉波振幅记录部65中,结束静止时脉波测量动作。
另一方面,在上述步骤S102中判定为所取得的3轴的各加速度信号的振幅比阈值大时,判断出用户US不处于静止状态或安静状态,进行催促用户US静止的错误显示(步骤S104)。
具体来说,信号处理部70在判定为所取得的3轴方向的各加速度信号的振幅的某个比上述阈值大时,判断出用户US不处于静止状态或安静状态,在显示部80中显示停止运动等动作,要求静止的意思的文字信息或图像信息等,催促用户US保持静止状态或安静状态。
接着,进行消除或舍弃在上述步骤S101中取得并保存在存储器部60中的脉波的观测信号以及加速度信号的复位动作(步骤S105)。
然后,返回步骤S101,再次执行上述一系列处理(步骤S101~S105)。
可以规定由通过这样的静止时脉波测量动作取得的静止时的观测信号、即各观测部位的脉波的观测信号的振幅的平均值组成的信号波形,是实质上不包含由于用户US的动作引起的体动噪音或者大体可以忽视体动噪音的状态的脉波信号。该脉波信号例如成为图6所示那样的波形。
此外,在图6中表示了在10秒期间观测静止时的脉波的情况下的信号波形的一例。
另外,在图6中,纵轴是对通过受光部20(受光元件)取得的观测信号进行A/D变换后的数字值。
(动作时脉波测量动作)
图7是表示在本实施方式的生物体信息检测装置的生物体信息检测方法中执行的动作时脉波测量动作的流程图。
图8是用于说明在本实施方式的动作时脉波测量动作中计算的极值间隔的概念图。
图9A、图9B、图9C是表示通过本实施方式的动作时脉波测量动作取得的各信号的一例的波形图。
在动作时脉波测量动作中,如图7所示,首先以一定时间取得用户US进行运动等动作的动作状态下的脉波的观测信号和加速度信号(步骤S201)。
具体来说,与上述静止时脉波测量动作同样,在用户US进行运动等动作的动作过程中,信号处理部70以一定时间取得在测量区域MS内排列的、相邻的发光元件和受光元件之间的各观测部位的脉波的观测信号。
在此,脉波的观测信号的取得动作与上述静止时脉波测量动作的步骤S101同样,在包含几个至十几个左右表示脉波的波形的任意时间内执行即可。该时间可以设定为与静止时脉波测量动作相同的时间(例如几秒至10秒左右),也可以设定为与之不同的时间。
另一方面,在脉波的观测信号的取得动作的期间中,信号处理部70持续取得用户US的动作引起的3轴方向的加速度信号。
所取得的脉波的观测信号和加速度信号,在根据时间数据相互关联后被保存在存储器部60的预定存储区域。
接着,与上述静止时脉波测量动作同样,判断在上述步骤S201中取得的3轴方向的各加速度信号的振幅的最大值(最大振幅)是否在与静止状态对应的预定阈值以下(步骤S202)。
接着,信号处理部70在判定为上述加速度信号的振幅在阈值以下时,判断出用户US处于静止状态。然后,从上述步骤S201中取得的各观测部位的脉波的观测信号中,将振幅最大的观测信号视为(判定为)最良好地测量了脉波的脉波信号,选择该观测信号(步骤S203)。
在此,在该步骤S203中选择的观测信号,由于是静止状态,因此判断为几乎不受体动噪音(加速度成分)的影响,可以视为与在上述静止时脉波测量动作中被保存在静止时脉波振幅记录部65中的观测信号(参照图6)具有同等或近似的信号波形。
接着,信号处理部70针对每一个波形从所选择的观测信号中搜索极值,计算其极值间隔(步骤S204)。
具体来说,在步骤S203中选择的观测信号例如具有图8所示的信号波形的情况下,信号处理部70计算在观测信号中包含的各波形中的振幅例如为极小值Pmin的时间Ta、Tb相互的差分时间,作为上述极值间隔。
此外,步骤S204中的极值间隔的计算动作,可以对所选择的观测信号中包含的波形中的任意时间的波形(即代表波形)来执行。或者,可以是将针对观测信号的一定时间内包含的多个波形所计算出的多个极值间隔进行平均化的结果(平均值)、或从多个极值间隔的分布中提取中央值的结果。
接着,根据在上述步骤S204中计算出的极值间隔,计算每单位时间(例如1分钟)的脉拍数(步骤S205)。
具体来说,信号处理部70在根据所选择的观测信号计算出的极值间隔的时间单位是秒时,通过用极值间隔分割(除)60,换算成1分钟的脉拍数。
接着,信号处理部70在显示部80中通过数值信息或图像信息等显示计算出的脉拍数,提供或通知给用户US(步骤S206)。
接着,在继续后续脉拍数的测量时,返回步骤S201。另一方面,在不继续测量时(结束时),结束动作时脉波测量动作(步骤S207)。
此外,在上述步骤S203中,表示了从在各观测部位取得的多个脉波的观测信号中选择一个振幅最大的观测信号的方法。但是,本发明不限定于此。
在本发明的生物体信息检测方法中,例如可以针对多个观测信号的各个观测信号计算极值间隔并换算成脉拍数,最终针对这些多个脉拍数取平均值或中央值等,提供给用户US。
另一方面,在上述步骤S202中判定为所取得的3轴方向的各加速度信号的振幅比阈值大的情况下,判断出该观测信号受到了体动噪音(加速度成分)的影响。
在这种情况下执行以下所示的降低体动噪音的影响的处理。
具体来说,首先,信号处理部70从上述步骤S201中取得的多个脉波的观测信号中,将与静止时的观测信号的振幅最近似的观测信号视为(判定为)体动噪音的影响最少的脉波信号,选择该观测信号(步骤S208)。
通过这样的观测信号的选择处理,可以减小脉波信号由于体动噪音而几乎完全消失的风险。换句话说,这可以避免脉波信号由于体动噪音被完全消除,无法判别的状态。
在此,在该步骤S208中选择的观测信号,例如具有图9A中实线所示的混合了脉波成分和体动噪音成分的信号波形。
在图9A中,虚线是不包含体动噪音或者大致可以忽视体动噪音的状态的脉波信号(例如通过上述静止时脉波测量动作取得的观测信号;以下记作“参照脉波信号”)。
在此,在图9A所示的观测信号(实线)的情况下,由于体动噪音的影响,其相位从参照脉波信号的相位偏移。
接着,信号处理部70执行推定在上述步骤S208中从用户US的动作发生的时刻到该动作所引起的加速度的影响显现在所选择的加速度信号的合成波形以及脉波的观测信号中为止的时间差(时滞)、与观测部位的主要血流方向和加速度信号的轴方向的角度差相对应的旋转角度的时滞、旋转角度推定处理(步骤S300)。
图10是表示在本实施方式的动作时脉波测量动作中执行的时滞、旋转角度推定处理的流程图。
图11是用于说明在本实施方式的动作时脉波测量动作中定义的3轴方向的概念图。
图12是表示通过本实施方式的时滞、旋转角度推定处理计算出的归一化相互相关系数的一例的图。
图13是表示通过本实施方式的时滞、旋转角度推定处理取得的旋转角度和极大值的迁移的一例的图。
首先,说明加速度信号的合成。
考虑到上述时滞的合成加速度信号,可以使用如下的数学式(1)来计算。
【数学式1】
A(t)=c1×Ax(t-d1)+c2×Ay(t-d2)+c3×Az(t-d3)…(1)
在此,A(t)是合成3轴方向的加速度信号所得的合成加速度信号,与观测信号中包含的加速度成分相对应。
Ax、Ay、Az分别是x轴方向、y轴方向、z轴方向的加速度信号,t表示时刻。
c1、c2、c3分别是与加速度信号Ax、Ay、Az相乘的比例系数,是设定合成加速度信号A(t)的振幅的系数。
d1、d2、d3分别表示直到在脉波的观测信号中显现加速度信号Ax、Ay、Az的影响为止的时间差(时滞)。
在此,关于x轴、y轴、z轴例如如图11所示,将手腕USh的长轴方向(腕的延伸方向;附图左右方向)定义为x轴方向,将与该x轴方向垂直的手腕USh的短轴方向(腕的宽度方向;附图左上右下方向)定义为y轴方向,将与x、y轴方向垂直的手腕USh的表里方向(附图上下方向)定义为z轴方向。
即,x轴和y轴被规定为沿着手腕USh的皮肤表面SF的方向。
合成在图11中定义的x、y、z的3轴方向的加速度信号所得的合成加速度信号A(t),在原理上可以使用上述数学式(1)计算。
但是,本申请的发明人根据各种验证的结果发现:(a)z轴方向的加速度信号Az对脉波信号几乎没有影响;(b)时滞d1、d2、d3几乎不依赖于轴方向,得到大致同等的值;(c)合成加速度信号A(t)通过使x轴方向的加速度信号Ax和y轴方向的加速度信号Ay旋转来计算,得到与真(本来)的合成加速度信号大致同等的值。
在此,可以通过轴的旋转来规定x轴方向和y轴方向的加速度信号的系数之比的理由,考虑是由于存在于观测部位的皮下(皮肤表面SF的下层)的多个动脉或毛细血管中的主要血流方向(即图11所示的血管VS的延长方向)因每个观测部位而不同。
即,图11所示的旋转角度θ对应于观测部位的主要血流方向和加速度信号的轴方向(图11中为x轴方向)的角度差。
根据这样的验证的结果,上述数学式(1)所示的3轴方向的合成加速度信号A(t)可以使用如下的数学式(2)来计算。
【数学式2】
A(t)=c×cosθ×Ax(t-d)-c×sinθ×Ay(t-d)…(2)
并且,步骤S300中执行的时滞、旋转角度推定处理,在上述数学式(2)中将与加速度信号Ax、Ay相乘的比例系数c固定为c=1的状态下推定时滞d的值、以及加速度信号Ax、Ay的旋转角度θ的值。
在时滞、旋转角度推定处理中,如图10所示,首先,信号处理部70设定相对于x轴方向(x轴)和y轴方向(y轴)的加速度信号的旋转角度θ(步骤S301)。
该旋转角度θ,每当重复后述的一系列处理(步骤S301~S305)时,在-90°(=-π/2)~+90°(=π/2)的范围内每次以预定角度依次被更新(增加或减少)。由此,搜索最佳的旋转角度θ。
在此,为了简化说明,作为初始值的一例,表示将旋转角度θ设定为0°,从0°到+90°以预定间隔依次增加角度的情况。
信号处理部70将不产生时滞的状态(即时滞d=0)设定为初始状态。
接着,信号处理部70根据设定为初始值的旋转角度θ(=0°)和比例系数c=1、时滞d=0,使用上述数学式(2)合成x轴方向的加速度信号Ax(t)和y轴方向的加速度信号Ay(t)(步骤S302)。
在此,在该步骤S302中生成的合成加速度信号A(t),例如成为图9B中实线所示的信号波形。在图9B中,虚线是上述参照脉波信号。
接着,信号处理部70根据在上述步骤S208中选择的观测信号和在步骤S302中生成的合成加速度信号A(t),计算对应于时滞d的归一化相互相关系数(步骤S303)。
在该步骤S303中计算出的归一化相互相关系数,例如如图12所示。
在此,存在以下因果关系,即用户US的动作中产生某加速度的结果,其经过某程度的时间后影响到脉波的观测信号。该时间是时滞。
信号处理部70,在图12所示的归一化相互相关系数中,在该因果关系成立的方向(图12的情况下,时滞为从0到正的方向),以预定的间隔依次更新时滞的值。由此,搜索相关系数最初达到极大值Dmax的时滞的值。
然后,提取出相关系数最初达到极大值Dmax时的时滞d,将其保存在存储器部60的预定存储区域中(步骤S304)。
此外,在图12中用粗线表示达到极大值Dmax的位置。
在此,在图12所示的归一化相互相关系数中,提取时滞d的范围,例如可以是使时滞d的值依次增加,在相关系数达到极大值Dmax的时刻结束处理的范围。或者,可以将时滞d的值规定为不达到特定时间例如1秒以上,在该时间之前计算归一化相互相关系数,此后,在该时间范围内求出相关系数的极大值Dmax。
接着,判断在上述步骤S304中计算出的此次的相关系数的极大值是否比前次计算出的极大值小(步骤S305)。
具体来说,信号处理部70从存储器部60读出此次和前次的相关系数的极大值。然后,在判定此次的极大值比前次的极大值小(即,前次的极大值比此次的极大值大)时,将前次的极大值的位置的时滞d和前次的旋转角度θ保存(记录)在存储器部60的预定存储区域中(步骤S306),结束时滞、旋转角度推定处理。
另一方面,在上述步骤S305中信号处理部70判定此次的相关系数的极大值为前次的极大值以上时,将此次的极大值的位置的时滞d和此次的旋转角度θ保存在存储器部60的预定存储区域中。
然后,返回步骤S301,再设定旋转角度θ。然后,再次执行上述的一系列搜索处理(步骤S301~S305)。
此外,在上述步骤S305中,在初次的判断处理的情况下,由于不存在前次的极大值,因此,在这种情况下无条件地返回步骤S301。然后,再设定旋转角度θ后,再次执行上述的一系列搜索处理(步骤S301~S305)。
此外,在上述的时滞、旋转角度推定处理中应用在步骤S305中所示的判定处理,是由于根据本申请发明人的验证,相关系数的极大值相对于旋转角度θ的变化得到了呈现单峰性的结果。但是,本发明不限于此。
例如,针对-90°~+90°的范围内的全部旋转角度θ计算归一化相互相关系数。并且,可以应用从该计算结果中选择相关系数的极大值达到最大的位置(图中为Pmax)的时滞d和此时的旋转角度θ,将其保存在存储器部60中的方法。
应用了这种方法的情况下的旋转角度θ和相关系数的极大值的关系(迁移),例如成为图13所示的形状。
在本实施方式中,作为旋转角度θ的设定范围,说明了在-90°~+90°的180°的范围内更新的情况。但是,本发明不限于此。旋转角度θ的设定范围至少具有180°的范围即可,例如可以将360°(全周)作为设定范围。
接着,信号处理部70根据在上述时滞、旋转角度推定处理中推定的时滞d以及旋转角度θ,推定设定合成加速度信号的振幅的比例系数。并且,执行生成去除了用户US的动作引起的体动噪音的影响的、与真的脉波信号近似的信号的振幅推定处理(步骤S400)。
即,在振幅推定处理中,执行在上述数学式(2)中,作为设定合成加速度信号A(t)的振幅的系数,推定与x、y各方向的加速度信号Ax、Ay相乘的比例系数c的处理。
图14是表示在本实施方式的动作时脉波测量动作中执行的振幅推定处理的流程图。
在振幅推定处理中,如图14所示,首先,信号处理部70设定在上述的数学式(2)中与x、y各方向的加速度信号Ax、Ay相乘的比例系数c(步骤S401)。
在此,首先,在上述数学式(2)中应用在上述的时滞、旋转角度推定处理中推定的时滞d、加速度信号的旋转角度θ的值,计算将比例系数设定为1时的合成加速度信号A(t)的振幅。
然后,比较该合成加速度信号A(t)的振幅和在上述步骤S208中选择出的脉波的观测信号的振幅。然后,计算合成加速度信号A(t)的振幅与观测信号的振幅相等的比例系数c的值,作为比例系数c的初始值,将比例系数c设定为该初始值。
接着,信号处理部70根据所设定的比例系数c、上述时滞、旋转角度推定处理中推定的时滞d、加速度信号的旋转角度θ的值,使用上述数学式(2)生成合成加速度信号A(t)(步骤S402)。然后,信号处理部70取得在上述步骤S208中选择的脉波的观测信号和在上述步骤S402中生成的合成加速度信号A(t)的差,生成差分的信号,作为差分信号(步骤S403)。
接着,信号处理部70计算所生成的差分信号的振幅(步骤S404)。
然后,判断在上述静止时脉波测量动作中取得的静止时的脉波的观测信号的振幅与上述差分信号的振幅的差的绝对值是否比预定的阈值小(步骤S405)。
信号处理部70在判断出上述绝对值比阈值小时,将此时的比例系数c的值保存(记录)在存储器部60的预定存储区域中(步骤S406),结束振幅推定处理。
另一方面,在上述步骤S405中,信号处理部70在判断出上述绝对值在阈值以上的情况下,将比例系数c的值再设定为别的值,然后更新(步骤S407)。
然后,返回步骤S402,再次执行上述一系列处理(步骤S402~S405)。
在此,再设定的比例系数c的值,以预定间隔依次增加或减少。
在此,在生成差分信号的处理(步骤S403)中生成的差分信号,例如具有图9C中实线所示那样的信号波形。图9C中虚线是上述参照脉波信号。
在此,图9C表示通过执行上述的时滞、旋转角度推定处理以及振幅推定处理的一系列处理,作为上述差分信号,得到与参照脉波信号的相位大体一致的信号波形的情况。
此外,在上述振幅推定处理中应用了以静止时的观测信号的振幅为基准,判断是否重复一系列处理的方法。但是,本发明不限于此。
例如,在认为与脉波信号的振幅相比,体动噪音的振幅足够大的情况下,可以应用依次更新比例系数c的值来搜索差分信号的振幅的最小值,并据此判断是否重复一系列处理的方法。
然后,在上述时滞、旋转角度推定处理(步骤S300)以及振幅推定处理(步骤S400)的结束后,如图7所示,信号处理部70将生成的差分信号视为脉波信号,计算极值间隔(步骤S209)。
此外,步骤S209中的极值间隔的计算动作与上述步骤S204同样,可以对所生成的差分信号中包含的波形中、任意时间的波形来执行。可以是将针对多个波形计算出的多个极值间隔进行平均化所得的结果(平均值)、或从多个极值间隔的分布中提取中央值的结果。
接着,根据在上述步骤S209中计算出的极值间隔,计算1分钟的脉拍数(步骤S205)。
接着,通过在显示部80中显示计算出的脉拍数,提供或通知给用户US(步骤S206)。
如上所述,在本实施方式中,在从通过光电脉波法在用户US的运动时取得的脉波的观测信号中去除由该运动引起的体动噪音成分来计算脉拍数的方法中,当运动时取得的3轴方向的加速度信号的振幅超过了预定的阈值时,使用从脉波的观测信号中去除沿着体表的特定方向(包含x轴以及y轴的x-y平面中的旋转角度θ)的加速度成分后的信号(差分信号)来计算脉拍数。
在此,在本实施方式中,当求出要从脉波的观测信号中去除的加速度成分时,应用推定与脉波的观测信号的时间差(时滞d)、决定振幅的大小的系数(比例系数c)、各方向的加速度的旋转角度θ这3个参数的方法。
此外,各参数的值根据用户US进行怎样的动作,即根据身体的各部位的动作状态而变化。例如,在慢走的状态、快走的状态、不抖动手腕的状态、大幅度抖动手腕的状态等状态下,加速度对脉波的观测信号的影响方式不同,上述各参数的值变化。因此,每次进行脉波的测量时,需要进行该各参数的值的推定。
在本实施方式中,具有在测量区域中排列多个发光元件和受光元件的结构,进行在测量区域内的不同的多个观测部位测量脉波的多点观测,在运动时取得的3轴方向的加速度信号的振幅超过预定阈值时,重新选择体动噪音的影响少的观测部位的脉波的观测信号来计算脉拍数。
在此,观测信号的选择,应用根据静止时在各观测部位取得的观测信号的振幅,选择体动噪音的影响最小的观测信号的方法。
这样,在本实施方式中,通过从运动中的脉波的观测信号中去除根据新推定的参数计算出的加速度信号(加速度成分),可以取得与真(本来)的脉波信号同相位的信号(差分信号)。并且,根据该差分信号可以比较准确地测量运动时的瞬时脉拍。
<第2实施方式>
接着,说明本发明的生物体信息检测装置的第2实施方式。
在此,适当参照上述附图来说明与第1实施方式同等的结构和动作。
在上述的第1实施方式中,说明了如图2所示具有在生物体信息检测装置100的测量区域MS中配置的发光元件以及受光元件中的至少某一方被配置了多个的结构,从通过多点观测而取得的多个脉波的观测信号中选择最佳的观测信号的情况。
在第2实施方式中,具有在测量区域MS中分别仅配置1个发光元件以及受光元件的结构,具有从1个地方的观测部位仅取得一个脉波的观测信号的方法(一点观测)。
图15A、图15B是表示第2实施方式的生物体信息检测装置的测量面的结构例的概要图,图15A是表示发光元件和受光元件的配置例的概要图,图15B是表示脉波的观测部位的概念图。
第2实施方式的生物体信息检测装置,在上述第1实施方式所示的结构(参照图1)中,如图15A所示,具有在设备本体101的测量区域MS中配置了1个发光元件E1和1个受光元件R1的结构。
即,在本实施方式中,发光元件和受光元件以1对1的关系排列。
并且,在图3所示的生物体信息检测装置的结构中,通过发光控制部15例如如图15B所示,使发光元件E1以预定的发光强度发光,使光照射到皮肤表面SF的观测部位Pm11,通过皮肤表面SF附近的血管中的血液被散射的光作为反射光由受光元件R1接受。由此,取得皮肤表面SF的观测部位Pm11的脉波的观测信号。
接着,说明本实施方式的生物体信息检测方法。
在此,适当参照上述附图来说明与第1实施方式同等的动作或处理。
本实施方式的生物体信息检测方法与上述第1实施方式同样地执行静止时脉波测量动作、动作时脉波测量动作。
首先,在本实施方式的静止时脉波测量动作中,在第1实施方式所示的图4的流程图中催促用户US静止。然后,使发光元件E1发光,通过受光元件R1接受其反射光。由此,以一定时间取得静止状态下的脉波的观测信号和加速度信号(步骤S101)。
在此,在本实施方式中具有发光元件和受光元件以1对1的关系排列的结构,因此,通过该步骤S101从一个地方的观测部位Pm11仅取得一个脉波的观测信号(一点观测)。
然后,在步骤S102中判定为该观测信号的取得时所测量的加速度信号的振幅在预定阈值以下的情况下,判断出用户US处于静止状态或安静状态,在步骤S101中取得的脉波的观测信号的振幅的平均值作为静止时的观测信号的振幅,被保存在静止时脉波振幅记录部65中(步骤S103)。
另一方面,在步骤S102中判定为观测信号的取得时所测量的加速度信号的振幅比预定阈值大的情况下,判断出用户US不处于静止状态或安静状态,与上述第1实施方式同样,在进行步骤S104、S105的动作后,再次执行上述的一系列处理(步骤S101~S105)。
图16是表示在本实施方式的生物体信息检测方法中执行的动作时脉波测量动作的流程图。
在本实施方式的动作时脉波测量动作中,相对于第1实施方式中表示的图7的流程图,省略了步骤S203以及S208的观测信号的选择处理。
即,在本实施方式的动作时脉波测量动作中,如图16的流程图所示,以一定时间取得用户US的运动状态下的脉波的观测信号和加速度信号(步骤S211)。在该步骤S211中也通过一点观测仅取得一个脉波的观测信号。
然后,在步骤S212中判定为该观测信号的取得时所测量的加速度信号的振幅在预定阈值以下时,在步骤S211中取得的脉波的观测信号被视为良好地测量了脉波的脉波信号。
然后,计算该观测信号的极值间隔(步骤S213)。
接着,根据在上述步骤S213中计算出的极值间隔,计算1分钟的脉拍数(步骤S214)。
然后,通过在显示部80中显示计算出的脉拍数,提供或通知给用户US(步骤S215)。
另一方面,在上述步骤S212中判定为所取得的加速度信号的振幅比阈值大的情况下,执行与上述第1实施方式同样的处理。即,执行推定在合成加速度信号A(t)以及脉波的观测信号中显现加速度的影响之前的时滞、以及与观测部位的主要血流方向和加速度信号的轴方向的角度差相对应的旋转角度的时滞、旋转角度推定处理(步骤S300)、以及生成去除了体动噪音的影响的与真的脉波信号近似的信号来作为差分信号的振幅推定处理(步骤S400)。
然后,将通过上述步骤S300以及S400的一系列处理生成的差分信号视为脉波信号,计算极值间隔(步骤S217)。
然后,根据在上述步骤S217中计算出的极值间隔,计算1分钟的脉拍数(步骤S214)。
然后,通过在显示部80中显示计算出的脉拍数,提供或通知给用户US(步骤S215)。
然后,在继续后续脉拍数的测量的情况下,再次执行上述一系列处理(步骤S211~S217)。
这样,在本实施方式中,也与上述第1实施方式同样地,从运动中的脉波的观测信号中去除根据新推定的与脉波的观测信号的时间差(时滞d)、决定振幅的大小的系数(比例系数c)、各轴方向的加速度信号的旋转角度θ这3个参数计算出的加速度信号(加速度成分),由此可以取得与真(本来)的脉波信号同相位的信号(差分信号)。并且,根据该差分信号可以比较准确地测量运动时的瞬时脉拍。
在此,在本实施方式中,具有在测量区域中分别排列1个发光元件以及受光元件的结构,应用了在一个地方的观测部位测量脉波的一点观测,因此,可以根据单一的脉波的观测信号通过简易的处理测量运动时的瞬时脉拍。
此外,在上述各实施方式中,说明了生物体信息检测装置100具有腕表型的形状,将具备测量区域MS的设备本体101以贴紧用户US的手腕USh的手背侧的方式安装的情况。但是,本发明不限于此。例如也可以贴紧安装在手心侧。
此外,如上述实施方式所示,在安装在手腕的手背侧的情况下,与安装在手心侧的情况相比,难以受到手腕的筋的***等导致的安装状态(测量区域与皮肤表面的贴紧状态)的变化的影响。因此,可以良好地取得脉波的观测信号。
在上述的各实施方式中,说明了生物体信息检测装置100具有腕表型的形状,安装在用户US的手腕USh的情况,但是,本发明不限于此。
即,本发明中,在测量区域MS中以预定的图案排列发光元件以及受光元件的生物体信息检测装置,被贴紧安装在可以观测人体的动作中的脉波的部位即可。
例如,可以具有在上述的手腕或上臂等的腕部、除去指尖的指部、耳朵、脚腕等观测部位,通过带等缠绕或者夹入来安装的形状。
以上,说明了本发明的若干实施方式,但是本发明不限于上述实施方式,也包含在请求专利保护的范围内记载的发明及其等同的范围。
Claims (13)
1.一种生物体信息检测装置,其特征在于,
具备:
检测部,其输出根据利用者的至少1个观测部位的脉波而检测出的观测信号;
加速度测量部,其输出伴随所述利用者的动作而测量的、与互不相同的多个轴方向的各个轴方向对应的多个加速度信号;
参数推定部,其基于根据多个参数合成所述多个加速度信号而得的合成加速度信号与所述观测信号的比较,推定在所述观测信号中包含的、与基于所述利用者的动作的加速度成分对应的所述各参数的特定的值;以及
脉拍数计算部,其根据从所述观测信号中去除与所述各参数的所述特定的值对应的特定的合成加速度信号所得的差分,计算所述利用者的脉拍数,
所述参数推定部,根据将所述各参数的值设定为互不相同的多个值时的、互不相同的多个所述合成加速度信号与所述观测信号的相互相关系数的值,推定所述各参数的所述特定的值。
2.根据权利要求1所述的生物体信息检测装置,其特征在于,
所述参数推定部,取得与所述多个合成加速度信号的各个相对应的所述相互相关系数的值,将该相互相关系数的值达到极大时的所述各参数的值推定为所述各参数的所述特定的值。
3.根据权利要求1所述的生物体信息检测装置,其特征在于,
所述多个参数包括:与从所述利用者的动作发生时刻到该动作的影响产生在所述观测信号中为止的时间差对应的第1参数;与所述多个加速度信号的各自的、所述各轴方向和所述观测部位的主要血流方向的角度差对应的第2参数,
所述参数推定部,推定所述第1参数的特定的值和所述第2参数的特定的值作为所述特定的值。
4.根据权利要求3所述的生物体信息检测装置,其特征在于,
具有存储部,其将在所述利用者不进行动作的静止状态下,通过所述检测部检测出的所述观测信号的振幅的值作为静止时振幅来存储,
所述多个参数还包括作为设定所述合成加速度信号的振幅的比例系数的第3参数,
所述参数推定部,基于所述静止时振幅与所述观测信号和所述合成加速度信号的所述差分的信号的振幅的比较,推定所述第3参数的特定的值作为所述特定的值。
5.根据权利要求4所述的生物体信息检测装置,其特征在于,
所述加速度测量部,
将沿着所述利用者的所述观测部位的体表面的、互相垂直的方向的x轴以及y轴、与所述x轴以及y轴垂直的方向的z轴这3轴设为所述多个轴方向,
取得所述x轴方向的第1加速度信号、所述y轴方向的第2加速度信号、所述z轴方向的第3加速度信号作为所述多个加速度信号,
所述参数推定部,至少使用所述第1加速度信号和所述第2加速度信号,来推定所述各参数的所述特定的值。
6.根据权利要求4所述的生物体信息检测装置,其特征在于,
所述检测部,基于所述利用者的互不相同的多个所述观测部位的各个脉波,输出多个所述观测信号,
所述生物体信息检测装置具备从所述多个观测信号中选择振幅与所述静止时振幅最近似的特定的观测信号的观测信号选择部,
所述参数推定部,根据所述特定的观测信号推定所述多个参数的特定的值。
7.根据权利要求6所述的生物体信息检测装置,其特征在于,
所述观测信号选择部,将所述多个加速度信号的各自的振幅与预定的阈值进行比较,当所述各加速度信号的振幅比所述阈值大时,从所述多个观测信号中进行所述特定的观测信号的选择。
8.根据权利要求6所述的生物体信息检测装置,其特征在于,
所述检测部具有:对所述多个观测部位的各个观测部位照射光的发光部;接受从该发光部照射并由所述多个观测部位的各个观测部位反射的光,输出所述多个观测信号的受光部,
所述发光部具备出射光的1个到多个发光元件,
所述受光部具备接受光的1个到多个受光元件,
该发光部以及该受光部至少将所述发光元件或所述受光元件的某一方具备多个。
9.根据权利要求1所述的生物体信息检测装置,其特征在于,
具备显示通过所述脉拍数计算部计算出的所述脉拍数的显示部。
10.一种生物体信息检测方法,其特征在于,
取得基于利用者的至少1个观测部位的脉波的观测信号,并且取得伴随该利用者的动作的、与互不相同的多个轴方向的各个轴方向对应的多个加速度信号,
基于根据多个参数合成所述多个加速度信号而得的合成加速度信号与所述观测信号的比较,推定在所述观测信号中包含的、与基于所述利用者的动作的加速度成分对应的所述各参数的特定的值,
根据从所述观测信号中去除与所述各参数的所述特定的值对应的特定的合成加速度信号所得的差分,计算所述利用者的脉拍数,
根据将所述各参数的值设定为互不相同的多个值时的、互不相同的多个所述合成加速度信号与所述观测信号的相互相关系数的值,进行所述各参数的所述特定的值的推定。
11.根据权利要求10所述的生物体信息检测方法,其特征在于,
在所述各参数的所述特定的值的推定中,取得与所述多个合成加速度信号的各个对应的所述相互相关系数的值,将该相互相关系数的值达到极大时的所述各参数的值推定为所述各参数的所述特定的值。
12.根据权利要求10所述的生物体信息检测方法,其特征在于,
所述多个参数包括:与从所述利用者的动作发生时刻到该动作的影响产生在所述观测信号中为止的时间差对应的第1参数;与所述各加速度信号的所述各轴方向和所述观测部位的主要血流方向的角度差对应的第2参数,
推定所述第1参数的特定的值和所述第2参数的特定的值作为所述多个参数的特定的值。
13.根据权利要求12所述的生物体信息检测方法,其特征在于,
所述多个参数还包括作为设定所述合成加速度信号的振幅的比例系数的第3参数,
将在所述利用者不进行动作的静止状态下检测出的所述观测信号的振幅的值作为静止时振幅来存储,
基于所述静止时振幅与所述观测信号和所述合成加速度信号的所述差分的信号的振幅的比较,推定所述第3参数的特定的值作为所述多个参数的特定的值。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2013021247A JP5979604B2 (ja) | 2013-02-06 | 2013-02-06 | 生体情報検出装置及び生体情報検出方法、生体情報検出プログラム |
JP2013-021247 | 2013-02-06 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103961080A CN103961080A (zh) | 2014-08-06 |
CN103961080B true CN103961080B (zh) | 2016-01-20 |
Family
ID=51231446
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410044180.3A Active CN103961080B (zh) | 2013-02-06 | 2014-01-30 | 生物体信息检测装置以及生物体信息检测方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US20140221848A1 (zh) |
JP (1) | JP5979604B2 (zh) |
CN (1) | CN103961080B (zh) |
Families Citing this family (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20170164878A1 (en) * | 2012-06-14 | 2017-06-15 | Medibotics Llc | Wearable Technology for Non-Invasive Glucose Monitoring |
WO2015109239A1 (en) | 2014-01-17 | 2015-07-23 | Medtronic, Inc. | Movement disorder symptom control |
US20150374245A1 (en) * | 2014-06-26 | 2015-12-31 | Heptagon Micro Optics Pte. Ltd. | Monitoring arterial blood flow |
US10552746B2 (en) | 2014-09-25 | 2020-02-04 | International Business Machines Corporation | Identification of time lagged indicators for events with a window period |
KR102302844B1 (ko) * | 2014-10-23 | 2021-09-16 | 삼성전자주식회사 | 정맥 패턴을 이용한 사용자 인증 방법 및 장치 |
WO2016106768A1 (zh) * | 2014-12-30 | 2016-07-07 | 北京阿纽山医药科技有限公司 | 一种智能健康腕带 |
CN104814726B (zh) * | 2015-05-28 | 2016-08-24 | 京东方科技集团股份有限公司 | 一种智能手环 |
US20170164847A1 (en) * | 2015-12-15 | 2017-06-15 | Texas Instruments Incorporated | Reducing Motion Induced Artifacts in Photoplethysmography (PPG) Signals |
JP6645190B2 (ja) | 2016-01-04 | 2020-02-14 | オムロンヘルスケア株式会社 | 診断支援装置、診断支援方法、及び、診断支援プログラム |
CN105832533B (zh) * | 2016-06-06 | 2018-06-08 | 厦门精图信息技术有限公司 | 一种智能的穿戴式心脏保护器 |
US10960202B2 (en) * | 2016-09-27 | 2021-03-30 | Medtronic, Inc. | Adaptive deep brain stimulation using movement desynchronization |
JP6723132B2 (ja) * | 2016-09-29 | 2020-07-15 | ルネサスエレクトロニクス株式会社 | 脈拍測定装置、光量制御方法、及びプログラム |
EP3530181B1 (en) | 2016-11-02 | 2023-10-04 | Huawei Technologies Co., Ltd. | Smart wearable device |
CN106388786A (zh) * | 2016-12-14 | 2017-02-15 | 广东技术师范学院 | 一种基于智能指环采集数据的处理方法 |
JP2018099409A (ja) * | 2016-12-21 | 2018-06-28 | セイコーエプソン株式会社 | 測定装置および測定方法 |
EP3417770A1 (en) * | 2017-06-23 | 2018-12-26 | Koninklijke Philips N.V. | Device, system and method for detection of pulse and/or pulse-related information of a patient |
US10617311B2 (en) * | 2017-08-09 | 2020-04-14 | Samsung Electronics Co., Ltd. | System and method for real-time heartbeat events detection using low-power motion sensor |
CN108175388B (zh) * | 2017-12-01 | 2021-01-26 | 中国联合网络通信集团有限公司 | 基于可穿戴设备的行为监测方法及装置 |
TWI676464B (zh) * | 2017-12-15 | 2019-11-11 | 劉佳星 | 電容式加速度計裝置及其感測方法 |
EP3570260A1 (en) * | 2018-05-15 | 2019-11-20 | Koninklijke Philips N.V. | Wrist-worn medical alert device for conveying an emergency message to a caregiver |
KR102640331B1 (ko) * | 2018-10-19 | 2024-02-26 | 삼성전자주식회사 | 생체정보 추정 장치 및 방법과, 생체정보 추정 지원 장치 |
CN113397497A (zh) * | 2021-06-10 | 2021-09-17 | 维沃移动通信有限公司 | 信号处理方法、装置和电子设备 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6361501B1 (en) * | 1997-08-26 | 2002-03-26 | Seiko Epson Corporation | Pulse wave diagnosing device |
CN1954771A (zh) * | 2003-03-19 | 2007-05-02 | 精工爱普生株式会社 | 脉搏计及其控制方法、以及手表型信息装置 |
CN101803925A (zh) * | 2010-03-31 | 2010-08-18 | 上海交通大学 | 运动状态下的血氧饱和度监测装置 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS6384517A (ja) * | 1986-09-30 | 1988-04-15 | アイシン精機株式会社 | 車上心拍計 |
GB9011887D0 (en) * | 1990-05-26 | 1990-07-18 | Le Fit Ltd | Pulse responsive device |
JP2001008908A (ja) * | 1999-06-28 | 2001-01-16 | Omron Corp | 電子血圧計 |
JP2003102694A (ja) * | 2001-09-28 | 2003-04-08 | Kenji Sunakawa | 心拍数測定装置 |
JP3726832B2 (ja) * | 2003-03-19 | 2005-12-14 | セイコーエプソン株式会社 | 脈拍計、腕時計型情報機器、制御プログラムおよび記録媒体 |
JP4345459B2 (ja) * | 2003-12-01 | 2009-10-14 | 株式会社デンソー | 生体状態検出装置 |
JP5060186B2 (ja) * | 2007-07-05 | 2012-10-31 | 株式会社東芝 | 脈波処理装置及び方法 |
CN102551685B (zh) * | 2010-12-30 | 2015-04-01 | 世意法(北京)半导体研发有限责任公司 | 对象监视器 |
-
2013
- 2013-02-06 JP JP2013021247A patent/JP5979604B2/ja active Active
-
2014
- 2014-01-30 CN CN201410044180.3A patent/CN103961080B/zh active Active
- 2014-02-04 US US14/172,577 patent/US20140221848A1/en not_active Abandoned
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6361501B1 (en) * | 1997-08-26 | 2002-03-26 | Seiko Epson Corporation | Pulse wave diagnosing device |
CN1954771A (zh) * | 2003-03-19 | 2007-05-02 | 精工爱普生株式会社 | 脉搏计及其控制方法、以及手表型信息装置 |
CN101803925A (zh) * | 2010-03-31 | 2010-08-18 | 上海交通大学 | 运动状态下的血氧饱和度监测装置 |
Also Published As
Publication number | Publication date |
---|---|
JP2014150869A (ja) | 2014-08-25 |
CN103961080A (zh) | 2014-08-06 |
JP5979604B2 (ja) | 2016-08-24 |
US20140221848A1 (en) | 2014-08-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103961080B (zh) | 生物体信息检测装置以及生物体信息检测方法 | |
EP3456251B1 (en) | Bio-information measuring apparatus and bio-information measuring method | |
CN105263403B (zh) | 生物体信息处理装置以及生物体信息处理方法 | |
US10722177B2 (en) | Cardiovascular monitoring using combined measurements | |
EP3692895B1 (en) | Apparatus and method for estimating bio-information | |
JP6201469B2 (ja) | 生体情報処理装置、生体情報処理方法 | |
JP6599883B2 (ja) | 心拍数モニターシステム、心拍数モニタリング方法及びコンピュータプログラム | |
US20150057508A1 (en) | Physiological information measuring apparatus | |
CN106560155A (zh) | 用于测量生物信息的设备和方法 | |
CN107920766A (zh) | 睡眠现象监测 | |
JP2007007075A (ja) | 血圧測定装置 | |
JP2005160640A (ja) | 生体状態検出装置 | |
WO2017162686A1 (en) | Fetal movement monitoring system and method | |
CN109640798A (zh) | 血压测定装置、血压测定装置的控制方法和程序 | |
KR20190029889A (ko) | 생체정보 추정 장치 및 방법 | |
EP3878352A1 (en) | Apparatus and method for estimating bio-information | |
JP6197926B2 (ja) | 生体情報検出装置及び生体情報検出方法、生体情報検出プログラム | |
JP2018068556A (ja) | 胸腔内圧推定装置 | |
JP6066451B2 (ja) | 生体情報検出装置及び生体情報検出方法、生体情報検出プログラム | |
CN112515629A (zh) | 用于估计生物信息的设备 | |
CN110267587A (zh) | 光学感测装置和相应的光学感测方法 | |
EP3854301A1 (en) | Apparatus and method for estimating bio-information | |
JP2017225839A (ja) | 生体情報検出装置及び生体情報検出方法、生体情報検出プログラム | |
JP5970712B2 (ja) | 生体情報測定回路、装置、プログラム及び方法 | |
CN109924955B (zh) | 脑血管动力学参数的确定方法、装置、终端及存储介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |