JP2011072674A - Blood vessel state estimating apparatus, blood vessel state estimating method, and blood vessel state estimating program - Google Patents

Blood vessel state estimating apparatus, blood vessel state estimating method, and blood vessel state estimating program Download PDF

Info

Publication number
JP2011072674A
JP2011072674A JP2009228838A JP2009228838A JP2011072674A JP 2011072674 A JP2011072674 A JP 2011072674A JP 2009228838 A JP2009228838 A JP 2009228838A JP 2009228838 A JP2009228838 A JP 2009228838A JP 2011072674 A JP2011072674 A JP 2011072674A
Authority
JP
Japan
Prior art keywords
blood vessel
difference
pulse wave
vessel state
optical sensor
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
JP2009228838A
Other languages
Japanese (ja)
Other versions
JP5381586B2 (en
Inventor
Takayuki Yamaji
隆行 山地
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.)
Fujitsu Ltd
Original Assignee
Fujitsu 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 Fujitsu Ltd filed Critical Fujitsu Ltd
Priority to JP2009228838A priority Critical patent/JP5381586B2/en
Publication of JP2011072674A publication Critical patent/JP2011072674A/en
Application granted granted Critical
Publication of JP5381586B2 publication Critical patent/JP5381586B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

<P>PROBLEM TO BE SOLVED: To provide a blood vessel state estimating apparatus, a blood vessel state estimating method and a blood vessel state estimating program by which the state of a blood vessel can accurately and easily be estimated. <P>SOLUTION: This blood vessel state estimating apparatus includes two optical sensors. Also, the blood vessel state estimating apparatus detects the pulse wave of a human body conforming to the variation amount of light which is detected by synchronizing with the other optical sensor, and converts the pulse wave into an electric signal for each optical sensor. Then, the blood vessel state estimating apparatus executes an n-th differential (n is a natural number) to the electric signal for each optical sensor, and identifies characteristic points included in an n-th differential waveform which is the waveform obtained as the result of the n-th differential. Then, the blood vessel state estimating apparatus calculates a difference for the appearance periods of the characteristic points which have been identified by a characteristic point identifying section for each optical sensor. Then, the blood vessel state estimating apparatus estimates the state of the blood vessel of the human body conforming to the difference which has been calculated for each optical sensor, and outputs an estimation result. <P>COPYRIGHT: (C)2011,JPO&INPIT

Description

本発明は、血管状態推定装置、血管状態推定方法及び血管状態推定プログラムに関する。   The present invention relates to a blood vessel state estimation device, a blood vessel state estimation method, and a blood vessel state estimation program.

従来、血管の伸展性や弾力性が加齢により変化することを踏まえ、血管の状態を推定する推定装置がある。例えば、推定装置は、加速度脈波や血圧脈波を用いることで、血管の状態が何歳に相当するかを推定する。   Conventionally, there is an estimation device that estimates the state of a blood vessel based on the fact that the extensibility and elasticity of the blood vessel change with aging. For example, the estimation device estimates how old the blood vessel state corresponds to by using an acceleration pulse wave or a blood pressure pulse wave.

例えば、加速度脈波を用いて血管の状態を推定する場合について更に説明する。例えば、推定装置は、利用者の指尖容積脈波を検出し、検出した指尖容積脈波を2回微分することで、図16に示すような加速度脈波を生成する。なお、図16は、加速度脈波の一例について説明するための図である。図16の横軸は時間軸を示し、図16の縦軸は、指尖容積脈波を2回微分することで得られた値を示す。また、図16のa〜eは、加速度脈波に含まれる特徴点各々を示す。図16に示す例では、推定装置は、基線を基準とした場合における特徴点a、b、dの高さを算出し、特徴点b、dの高さをそれぞれ特徴点aの高さで除算する。ここで、除算することで得られた値である「b/a」や「d/a」と年齢との間には相関関係があることを踏まえ、推定装置は、「b/a」や「d/a」を用いて血管の状態を推定する。例えば、推定装置は、利用者の血管の状態が20代であると推定する。   For example, the case where the state of the blood vessel is estimated using the acceleration pulse wave will be further described. For example, the estimation device detects the fingertip volume pulse wave of the user, and generates the acceleration pulse wave as shown in FIG. 16 by differentiating the detected fingertip volume pulse wave twice. In addition, FIG. 16 is a figure for demonstrating an example of an acceleration pulse wave. The horizontal axis in FIG. 16 represents the time axis, and the vertical axis in FIG. 16 represents the value obtained by differentiating the fingertip volume pulse wave twice. Moreover, ae of FIG. 16 shows each feature point included in the acceleration pulse wave. In the example illustrated in FIG. 16, the estimation device calculates the heights of the feature points a, b, and d with respect to the baseline, and divides the heights of the feature points b and d by the height of the feature point a, respectively. To do. Here, based on the fact that there is a correlation between “b / a” or “d / a”, which is a value obtained by division, and age, the estimation device uses “b / a” or “ The state of the blood vessel is estimated using “d / a”. For example, the estimation device estimates that the blood vessel state of the user is in his twenties.

また、血圧脈波を用いて血管の状態を推定する場合について更に説明する。例えば、まず、利用者の両腕と両足首に予めカフを装着しておく。そして、推定装置は、血圧脈波が心臓から両腕に伝搬するまでの時間と、血圧脈波が心臓から両足首に伝搬するまでの時間との差を測定することで、血圧脈波が伝搬する速度である脈波伝搬速度を測定する。そして、脈波伝搬速度と年齢との間には相関関係があることを踏まえ、推定装置は、測定した脈波伝搬速度を用いて血管の状態を測定する。   Further, the case where the blood vessel state is estimated using the blood pressure pulse wave will be further described. For example, first, cuffs are previously attached to the user's arms and ankles. The estimation device then measures the difference between the time it takes for the blood pressure pulse wave to propagate from the heart to both arms and the time it takes for the blood pressure pulse wave to propagate from the heart to both ankles. Measure the pulse wave velocity, which is the velocity to perform. Then, based on the fact that there is a correlation between the pulse wave propagation speed and the age, the estimation apparatus measures the state of the blood vessel using the measured pulse wave propagation speed.

特開2006−102250号公報JP 2006-102250 A 特開2000−51166号公報JP 2000-511166 A

しかしながら、加速度脈波を用いた手法では、精度が悪いという課題があった。すなわち、脈波は様々な要因によって変化し、例えば、血圧値によって変化する。また、血圧値もまた、様々な要因によって変化し、例えば、利用者の呼吸によって変化する。このため、加速度脈波を用いる推定装置では、測定するごとに脈波が変化する結果、精度が悪かった。例えば、加速度脈波を用いる推定装置では、処理を行うごとに得られる年齢の幅が大きかった。このように、加速度脈波を用いた手法では、精度が悪いという課題があった。   However, the technique using the acceleration pulse wave has a problem that accuracy is poor. That is, the pulse wave varies depending on various factors, for example, varies depending on the blood pressure value. In addition, the blood pressure value also changes due to various factors, for example, changes due to the user's breathing. For this reason, in the estimation apparatus using the acceleration pulse wave, the accuracy is poor as a result of the pulse wave changing every measurement. For example, in the estimation apparatus that uses acceleration pulse waves, the range of ages obtained each time processing is performed is large. As described above, the technique using the acceleration pulse wave has a problem that accuracy is poor.

また、例えば、血圧脈波を用いた手法では、簡易に測定できないという課題があった。すなわち、推定装置では、利用者がカフを予め装着しなければならなかった。また、カフを装着すると、利用者の足首や両腕が締め付けられることになり、利用者への負荷が大きかった。また、カフを用いる場合には、装置を小型化することが難しかった。このように、血圧脈波を用いた手法では、簡易に測定できないという課題があった。   Further, for example, the technique using blood pressure pulse waves has a problem that it cannot be easily measured. That is, in the estimation apparatus, the user has to wear the cuff in advance. When the cuff is worn, the user's ankle and both arms are tightened, and the load on the user is large. In addition, when using a cuff, it is difficult to reduce the size of the apparatus. As described above, the technique using the blood pressure pulse wave has a problem that it cannot be easily measured.

開示の技術は、上記に鑑みてなされたものであって、血管の状態を精度良く簡易に推定可能な血管状態推定装置、血管状態推定方法及び血管状態推定プログラムを提供することを目的とする。   The disclosed technology has been made in view of the above, and an object of the present invention is to provide a blood vessel state estimation device, a blood vessel state estimation method, and a blood vessel state estimation program capable of accurately and simply estimating a blood vessel state.

開示する血管状態推定装置の一つの態様によれば、人体へ光を照射する発光素子と、前記発光素子から放射された後人体を透過した光の変化量、又は、前記発光素子から放射された後人体により反射された光の変化量を検出する受光素子とを有する光センサを備える。また、血管状態推定装置は、少なくとも2つの光センサを有する。また、血管状態推定装置は、前記光センサごとに、他の光センサと同期をとって検出された光の変化量に基づいて前記人体の脈波を検出し、該脈波を電気信号に変換する脈波検出部を備える。また、血管状態推定装置は、前記光センサごとに、前記脈波検出部によって変換された電気信号に対してn次微分(nは自然数)を実行し、該n次微分の結果得られた波形であるn次微分波形に含まれる特徴点を識別する特徴点識別部を備える。また、血管状態推定装置は、前記光センサごとに、前記特徴点識別部によって識別された特徴点の出現時間の差分を算出する差分算出部を備える。また、血管状態推定装置は、前記差分算出部によって前記光センサごとに算出された差分に基づいて前記人体の血管の状態を推定し、推定結果を出力する血管状態推定部を備える。   According to one aspect of the disclosed blood vessel state estimation device, a light emitting element that irradiates light to a human body, and a change amount of light that is emitted from the light emitting element and then passes through the human body, or is emitted from the light emitting element An optical sensor having a light receiving element for detecting a change amount of light reflected by the rear human body is provided. Further, the blood vessel state estimation device has at least two optical sensors. In addition, the blood vessel state estimation device detects a pulse wave of the human body based on a change amount of light detected in synchronization with another optical sensor for each optical sensor, and converts the pulse wave into an electrical signal. A pulse wave detection unit. In addition, the blood vessel state estimation device performs n-order differentiation (n is a natural number) on the electrical signal converted by the pulse wave detection unit for each of the optical sensors, and a waveform obtained as a result of the n-order differentiation. A feature point identifying unit for identifying feature points included in the nth-order differential waveform. In addition, the blood vessel state estimation device includes a difference calculation unit that calculates a difference between appearance times of feature points identified by the feature point identification unit for each of the optical sensors. The blood vessel state estimation device further includes a blood vessel state estimation unit that estimates a blood vessel state of the human body based on the difference calculated for each of the optical sensors by the difference calculation unit and outputs an estimation result.

開示する血管状態推定装置の一つの態様によれば、血管の状態を精度良く簡易に推定可能であるという効果を奏する。   According to one aspect of the disclosed blood vessel state estimation device, there is an effect that the blood vessel state can be easily and accurately estimated.

図1は、実施例1に係る血管状態推定装置の構成の一例について説明するブロック図である。FIG. 1 is a block diagram illustrating an example of the configuration of the blood vessel state estimation device according to the first embodiment. 図2は、実施例2に係る血管状態推定装置の構成の一例について説明するブロック図である。FIG. 2 is a block diagram illustrating an example of the configuration of the blood vessel state estimation device according to the second embodiment. 図3は、実施例2における光センサの使用態様の一例について説明する図である。FIG. 3 is a diagram illustrating an example of usage of the optical sensor in the second embodiment. 図4は、実施例2における光センサの構造の一例について説明する図である。FIG. 4 is a diagram illustrating an example of the structure of the optical sensor according to the second embodiment. 図5は、実施例2における搭載基板の詳細な一例について説明する図である。FIG. 5 is a diagram illustrating a detailed example of the mounting substrate in the second embodiment. 図6は、実施例2における光センサによって変換された電気信号が示す波形の一例について説明する図である。FIG. 6 is a diagram illustrating an example of a waveform indicated by an electrical signal converted by the optical sensor according to the second embodiment. 図7は、実施例2における閾値テーブルに記憶された情報の一例について説明する図である。FIG. 7 is a diagram illustrating an example of information stored in the threshold value table according to the second embodiment. 図8は、実施例2における脈波間隔算出部による脈飛びがあった箇所を識別する処理の一例について説明する図である。FIG. 8 is a diagram illustrating an example of a process for identifying a place where a pulse jump has occurred by the pulse wave interval calculation unit according to the second embodiment. 図9Aは、実施例2における2次微分波形の一例について説明する図である。FIG. 9A is a diagram illustrating an example of a second-order differential waveform in the second embodiment. 図9Bは、実施例2における2次微分波形の一例について説明する図である。FIG. 9B is a diagram illustrating an example of a second-order differential waveform in the second embodiment. 図10は、実施例2における特徴点識別部による特徴点識別処理について更に説明する図である。FIG. 10 is a diagram for further explaining the feature point identifying process by the feature point identifying unit in the second embodiment. 図11は、実施例2における差分算出部による出現時間の差分算出処理の一例について説明する図である。FIG. 11 is a diagram illustrating an example of a difference calculation process of appearance time by a difference calculation unit according to the second embodiment. 図12Aは、センサ間差分の分散と年齢との関係について説明するための図である。FIG. 12A is a diagram for explaining the relationship between the variance of the difference between sensors and the age. 図12Bは、センサ間差分の分散と年齢との関係について説明するための図である。FIG. 12B is a diagram for explaining the relationship between the variance of the inter-sensor difference and the age. 図12Cは、センサ間差分の分散と年齢との関係について説明するための図である。FIG. 12C is a diagram for explaining the relationship between the variance of differences between sensors and the age. 図13は、センサ間差分の分散と年齢との関係について説明するための図である。FIG. 13 is a diagram for explaining the relationship between the variance of differences between sensors and the age. 図14は、実施例2に係る血管状態推定装置による処理の流れの一例について説明するフローチャートである。FIG. 14 is a flowchart illustrating an example of a process flow performed by the blood vessel state estimation device according to the second embodiment. 図15は、実施例2に係る血管状態推定プログラムを実行するコンピュータの一例について説明する図である。FIG. 15 is a schematic diagram illustrating an example of a computer that executes a blood vessel state estimation program according to the second embodiment. 図16は、加速度脈波の一例について説明するための図である。FIG. 16 is a diagram for explaining an example of the acceleration pulse wave.

以下に、開示する血管状態推定装置、血管状態推定方法及び血管状態推定プログラムの実施例を図面に基づいて詳細に説明する。なお、本実施例により開示する発明が限定されるものではない。各実施例は、処理内容を矛盾させない範囲で適宜組み合わせることが可能である。   Hereinafter, embodiments of the disclosed vascular state estimation device, vascular state estimation method, and vascular state estimation program will be described in detail with reference to the drawings. Note that the invention disclosed by this embodiment is not limited. Each embodiment can be appropriately combined within a range in which processing contents do not contradict each other.

図1を用いて、実施例1に係る血管状態推定装置100の構成の一例について説明する。図1は、実施例1に係る血管状態推定装置の構成の一例について説明するブロック図である。血管状態推定装置100は、図1に示す例では、光センサ101と、光センサ102と、脈波検出部103と、特徴点識別部104と、差分算出部105と、血管状態推定部106とを有する。   An example of the configuration of the blood vessel state estimation apparatus 100 according to the first embodiment will be described with reference to FIG. FIG. 1 is a block diagram illustrating an example of the configuration of the blood vessel state estimation device according to the first embodiment. In the example illustrated in FIG. 1, the blood vessel state estimation device 100 includes an optical sensor 101, an optical sensor 102, a pulse wave detection unit 103, a feature point identification unit 104, a difference calculation unit 105, and a blood vessel state estimation unit 106. Have

光センサ101や光センサ102は、人体へ光を照射する発光素子と、発光素子から放射された後人体を透過した光の変化量、又は、発光素子から放射された後人体により反射された光の変化量を検出する受光素子とを有する。   The optical sensor 101 and the optical sensor 102 are a light emitting element that irradiates light to the human body, a change amount of light that is transmitted from the light emitting element and then transmitted through the human body, or light that is reflected from the rear human body that is emitted from the light emitting element. And a light receiving element for detecting a change amount of the light receiving element.

脈波検出部103は、光センサごとに、他の光センサと同期をとって検出された光の変化量に基づいて人体の脈波を検出し、脈波を電気信号に変換する。また、特徴点識別部104は、光センサごとに、脈波検出部103によって変換された電気信号に対してn次微分(nは自然数)を実行し、n次微分の結果得られた波形であるn次微分波形に含まれる特徴点を識別する。   For each optical sensor, the pulse wave detection unit 103 detects the pulse wave of the human body based on the change amount of light detected in synchronization with the other optical sensors, and converts the pulse wave into an electrical signal. In addition, the feature point identification unit 104 performs n-order differentiation (n is a natural number) on the electrical signal converted by the pulse wave detection unit 103 for each optical sensor, and uses the waveform obtained as a result of the n-order differentiation. A feature point included in a certain nth-order differential waveform is identified.

差分算出部105は、光センサごとに、特徴点識別部104によって識別された特徴点の出現時間の差分を算出する。また、血管状態推定部106は、差分算出部105によって光センサごとに算出された差分に基づいて人体の血管の状態を推定し、推定結果を出力する。   The difference calculation unit 105 calculates the difference in appearance time of the feature points identified by the feature point identification unit 104 for each optical sensor. The blood vessel state estimation unit 106 estimates the blood vessel state of the human body based on the difference calculated for each optical sensor by the difference calculation unit 105, and outputs an estimation result.

このように、実施例1によれば、血管の状態を精度良く推定可能である。すなわち、脈波は様々な要因によって変化し、例えば、血圧値によって変化する。また、血圧値もまた、様々な要因によって変化し、例えば、利用者の呼吸によって変化する。この結果、脈波の時間軸を基準とした場合における特徴点の高さを用いて血管の状態を推定すると、精度が悪かった。しかしながら、実施例1によれば、特徴点各々の出現時間は血圧値などによる変動を受けにくいことを踏まえ、特徴点の出現時間の差分に基づいて血管の状態を推定するので、精度良く推定することが可能である。   Thus, according to Example 1, the state of the blood vessel can be estimated with high accuracy. That is, the pulse wave varies depending on various factors, for example, varies depending on the blood pressure value. In addition, the blood pressure value also changes due to various factors, for example, changes due to the user's breathing. As a result, when the state of the blood vessel was estimated using the height of the feature point when the time axis of the pulse wave was used as a reference, the accuracy was poor. However, according to the first embodiment, based on the fact that the appearance time of each feature point is less susceptible to fluctuation due to the blood pressure value or the like, the state of the blood vessel is estimated based on the difference in the appearance time of the feature point. It is possible.

また、実施例1によれば、カフを用いる手法と比較して、血管の状態を簡易に推定可能である。   Further, according to the first embodiment, it is possible to easily estimate the state of the blood vessel as compared with the method using the cuff.

[血管状態推定装置の構成]
次に、実施例2に係る血管状態推定装置200について説明する。まず、図2を用いて、実施例2に係る血管状態推定装置200の構成の一例について説明する。図2は、実施例2に係る血管状態推定装置の構成の一例について説明するブロック図である。血管状態推定装置200は、図2に示す例では、光センサ201と、光センサ202と、出力部203と、記憶部300と、制御部400とを有する。
[Configuration of blood vessel state estimation device]
Next, the blood vessel state estimation device 200 according to the second embodiment will be described. First, an example of the configuration of the blood vessel state estimation apparatus 200 according to the second embodiment will be described with reference to FIG. FIG. 2 is a block diagram illustrating an example of the configuration of the blood vessel state estimation device according to the second embodiment. In the example illustrated in FIG. 2, the blood vessel state estimation device 200 includes an optical sensor 201, an optical sensor 202, an output unit 203, a storage unit 300, and a control unit 400.

光センサ201と光センサ202とは、制御部400と接続される。また、光センサ201や光センサ202は、人体へ光を照射する発光素子と、発光素子から放射された後人体を透過した光の変化量、又は、発光素子から放射された後人体により反射された光の変化量を検出する受光素子とを有する。また、光センサ201と光センサ202とは、それぞれ異なる太さの血管がある人体の部分に光を照射する。   The optical sensor 201 and the optical sensor 202 are connected to the control unit 400. In addition, the optical sensor 201 and the optical sensor 202 are reflected by a light emitting element that irradiates light to the human body and a change amount of light that has been emitted from the light emitting element and transmitted through the human body, or by the rear human body that is emitted from the light emitting element. And a light receiving element for detecting the amount of change in light. In addition, the optical sensor 201 and the optical sensor 202 irradiate light to a portion of a human body having blood vessels with different thicknesses.

ここで、発光素子は、例えば、LED(Light Emitting Diode)が該当する。また、受光素子は、例えば、PD(Photo Diode)が該当する。発光素子として用いられるLEDは、例えば、760nm付近の波長の光を照射する。なお、以下では、発光素子として用いられるLEDが760nm付近の波長の光を照射する場合を例に説明するが、本発明はこれに限定されるものではない。例えば、発光素子として用いられるLEDは、760nmよりも長波長側の光を照射しても良く、760nmよりも短波長側の光を照射しても良く、利用者が任意の波長を選択して良い。   Here, the light emitting element corresponds to, for example, an LED (Light Emitting Diode). The light receiving element corresponds to, for example, PD (Photo Diode). The LED used as the light emitting element irradiates light having a wavelength near 760 nm, for example. In the following description, an example in which an LED used as a light emitting element emits light having a wavelength near 760 nm will be described, but the present invention is not limited to this. For example, an LED used as a light emitting element may irradiate light having a wavelength longer than 760 nm, or irradiate light having a wavelength shorter than 760 nm, and the user selects an arbitrary wavelength. good.

図3を用いて、光センサ201や光センサ202の使用態様の一例について説明する。図3は、実施例2における光センサの使用態様の一例について説明する図である。図3では、血管状態推定装置200が移動体端末である場合を例に示した。図3に示す例では、光センサ201は、クリップ状の形状を有し、利用者によって耳朶に装着されることで耳と接触する。また、図3に示す例では、光センサ201は、ケーブルを介して血管状態推定装置200の筐体と接続される。また、図3に示す例では、光センサ202は、血管状態推定装置200の筐体に設けられ、利用者によって指が押しつけられることで指と接触する。   An example of usage of the optical sensor 201 and the optical sensor 202 will be described with reference to FIG. FIG. 3 is a diagram illustrating an example of usage of the optical sensor in the second embodiment. FIG. 3 shows an example in which the blood vessel state estimation device 200 is a mobile terminal. In the example shown in FIG. 3, the optical sensor 201 has a clip-like shape, and comes into contact with the ear when worn by the user on the earlobe. In the example illustrated in FIG. 3, the optical sensor 201 is connected to the housing of the blood vessel state estimation device 200 via a cable. In the example illustrated in FIG. 3, the optical sensor 202 is provided in the housing of the blood vessel state estimation device 200 and comes into contact with the finger when the finger is pressed by the user.

なお、以下では、図3に示すように、光センサ201は耳と接触し、光センサ202が指と接触する場合を例に説明するが、本発明はこれに限定されるものではない。例えば、光センサ201や光センサ202が接触する箇所は、利用者が任意に選択して良い。   In the following, as shown in FIG. 3, a case where the optical sensor 201 is in contact with the ear and the optical sensor 202 is in contact with the finger will be described as an example. However, the present invention is not limited to this. For example, the user may arbitrarily select a location where the optical sensor 201 or the optical sensor 202 contacts.

また、以下では、図3に示すように、光センサ201は、クリップ状の形状を有して耳に装着される場合を例に説明するが、本発明はこれに限定されるものではない。例えば、光センサ201は、クリップ状ではなく板状であっても良く、利用者が任意の形状を選択して良い。   In the following description, as shown in FIG. 3, the optical sensor 201 has a clip-like shape and is attached to an ear. However, the present invention is not limited to this. For example, the optical sensor 201 may have a plate shape instead of a clip shape, and the user may select an arbitrary shape.

また、以下では、図3に示すように、光センサ201は、ケーブルを介して血管状態推定装置200の筐体と接続される場合を例に説明するが、本発明はこれに限定されるものではない。例えば、光センサ201は、無線を介して血管状態推定装置200の筐体と接続されても良く、利用者が任意の接続形態を選択して良い。また、同様に、光センサ201は、血管状態推定装置100の筐体に設けられても良い。   In the following, as illustrated in FIG. 3, the optical sensor 201 will be described as an example in which the optical sensor 201 is connected to the housing of the blood vessel state estimation device 200 via a cable, but the present invention is not limited thereto. is not. For example, the optical sensor 201 may be connected to the housing of the blood vessel state estimation device 200 via wireless, and the user may select any connection form. Similarly, the optical sensor 201 may be provided in the housing of the blood vessel state estimation device 100.

また、以下では、図3に示すように、光センサ202は、血管状態推定装置200の筐体に設けられる場合を例に説明するが、本発明はこれに限定されるものではない。例えば、光センサ202は、ケーブルや無線を介して血管状態推定装置200の筐体と接続されても良く、利用者が任意の接続形態を選択して良い。   In the following, as illustrated in FIG. 3, the optical sensor 202 is described as an example provided in the housing of the blood vessel state estimation device 200, but the present invention is not limited to this. For example, the optical sensor 202 may be connected to the housing of the blood vessel state estimation device 200 via a cable or wirelessly, and the user may select an arbitrary connection form.

光センサ201の構造の一例について、図4を用いて説明する。図4は、実施例2における光センサの構造の一例について説明する図である。図4には、光センサ201について、光センサ201を形成する各部を分解した上で示した。図4に示すように、光センサ201は、カバーケース504に、開閉可能なふたとしてクリップ状に可動する開口制御ラバー501が取り付けられる。その上で、開口制御ラバー501のカバーケース504側の面には、可視光をカットする可視光カットフィルタ502と、LED(Light Emitting Diode)やPD(Photo Diode)が搭載される搭載基板503とが順に積層される。   An example of the structure of the optical sensor 201 will be described with reference to FIG. FIG. 4 is a diagram illustrating an example of the structure of the optical sensor according to the second embodiment. FIG. 4 shows the optical sensor 201 after disassembling each part forming the optical sensor 201. As shown in FIG. 4, in the optical sensor 201, an opening control rubber 501 that is movable in a clip shape as a lid that can be opened and closed is attached to a cover case 504. In addition, on the surface of the opening control rubber 501 on the cover case 504 side, a visible light cut filter 502 for cutting visible light, and a mounting substrate 503 on which an LED (Light Emitting Diode) and a PD (Photo Diode) are mounted. Are sequentially stacked.

ここで、可視光カットフィルタ502は、開口制御ラバー501を透過した可視光線を遮断することで、搭載基板503に搭載されたPDが、開口制御ラバー501を透過した可視光線を検知しないようにする。搭載基板503に搭載されたLEDは、発光素子に対応し、搭載基板503に搭載されたPDは、受光素子に対応する。なお、光センサ202は、例えば、搭載基板503が血管状態推定装置200の筐体に設けられることで実現される。   Here, the visible light cut filter 502 blocks the visible light transmitted through the aperture control rubber 501 so that the PD mounted on the mounting substrate 503 does not detect the visible light transmitted through the aperture control rubber 501. . The LED mounted on the mounting substrate 503 corresponds to a light emitting element, and the PD mounted on the mounting substrate 503 corresponds to a light receiving element. Note that the optical sensor 202 is realized, for example, by providing the mounting substrate 503 in the housing of the blood vessel state estimation device 200.

また、図5を用いて、搭載基板503の詳細な一例について説明する。図5は、実施例2における搭載基板の詳細な一例について説明する図である。図5に示す例では、光センサ201の搭載基板503である場合を例に示した。つまり、図5に示す例では、搭載基板503が耳507と接触する場合を例に示した。   A detailed example of the mounting substrate 503 will be described with reference to FIG. FIG. 5 is a diagram illustrating a detailed example of the mounting substrate in the second embodiment. In the example illustrated in FIG. 5, the case where the optical sensor 201 is the mounting substrate 503 is illustrated as an example. That is, in the example illustrated in FIG. 5, the case where the mounting substrate 503 is in contact with the ear 507 is illustrated as an example.

図5に示すように、搭載基板503の平面上には、LEDなどの発光素子505と、PDなどの受光素子506とが設けられる。また、搭載基板503は、耳507と接触する。ここで、図5に示すように、受光素子506は、発光素子505から照射された光のうち、耳507によって反射された光を受光し、受光した光の変化量を検知する。あるいは、受光素子506は、発光素子505から照射された光のうち、耳507を透過した光を受光し、受光した光の変化量を検知する。   As shown in FIG. 5, a light emitting element 505 such as an LED and a light receiving element 506 such as a PD are provided on the plane of the mounting substrate 503. Further, the mounting substrate 503 is in contact with the ear 507. Here, as shown in FIG. 5, the light receiving element 506 receives the light reflected by the ear 507 among the light emitted from the light emitting element 505, and detects the amount of change in the received light. Alternatively, the light receiving element 506 receives light transmitted through the ear 507 among the light emitted from the light emitting element 505, and detects the amount of change in the received light.

光センサ201は、光センサ202と同期をとって検出された光の変化量に基づいて人体の脈波を検出し、脈波を電気信号に変換する。例えば、光センサ201は、耳にて脈波を検出し、電気信号に変換する。また、光センサ202は、光センサ201と同期をとって検出された光の変化量に基づいて人体の脈波を検出し、脈波を電気信号に変換する。例えば、光センサ202は、指にて脈波を検出し、電気信号に変換する。   The optical sensor 201 detects a human body pulse wave based on the amount of change in light detected in synchronization with the optical sensor 202, and converts the pulse wave into an electrical signal. For example, the optical sensor 201 detects a pulse wave with an ear and converts it into an electrical signal. The optical sensor 202 detects a pulse wave of the human body based on the change amount of light detected in synchronization with the optical sensor 201, and converts the pulse wave into an electrical signal. For example, the optical sensor 202 detects a pulse wave with a finger and converts it into an electrical signal.

光センサ201や光センサ202が脈波を検出する処理について、更に説明する。まず、発光素子が光を照射する箇所の血流量と、受光素子が受光する光の量との間には、相関関係にある。具体的には、発光素子が光を照射する箇所の血流量が増加すればするほど、発光素子が光を照射する箇所の血流によって吸収される光の量が増加し、受光素子が受光する光の量が減少する。次に、心臓の鼓動により発生した脈波が到達すると、血管が広がり、発光素子が光を照射する箇所の血流量が増加する。つまり、光センサ201や光センサ202は、受光素子によって受光される光の変化量を検出することで、発光素子が光を照射する箇所の血流量の変化を検知する。また、この結果、光センサ201や光センサ202は、血流量の変化を検知することで、脈波を検出する。   Processing for detecting the pulse wave by the optical sensor 201 and the optical sensor 202 will be further described. First, there is a correlation between the blood flow volume at the location where the light emitting element irradiates light and the amount of light received by the light receiving element. Specifically, as the blood flow rate at the location where the light emitting element irradiates light increases, the amount of light absorbed by the blood flow at the location where the light emitting element illuminates increases, and the light receiving element receives light. The amount of light decreases. Next, when the pulse wave generated by the heartbeat arrives, the blood vessels widen, and the blood flow volume at the location where the light emitting element emits light increases. That is, the optical sensor 201 or the optical sensor 202 detects a change in blood flow at a location where the light emitting element emits light by detecting the amount of change in the light received by the light receiving element. As a result, the optical sensor 201 and the optical sensor 202 detect a pulse wave by detecting a change in blood flow.

図6を用いて、光センサ201や光センサ202によって変換された電気信号が示す波形の一例について説明する。なお、図6は、光センサによって変換された電気信号が示す波形の一例について説明する図である。ここで、図6の(1)は、光センサ201によって変換された電気信号が示す波形の一例を示す。言い換えると、図6の(1)は、耳と接触した光センサについての波形である。図6の(2)は、光センサ202によって変換された電気信号が示す波形の一例を示す。言い換えると、図6の(2)は、指と接触した光センサについての波形である。図6の(1)や(2)の横軸は、時間軸を示し、図6の(1)や(2)の縦軸は、受光素子によって受光された光の量を示す。つまり、縦軸は、光センサが光を照射する箇所における血流量を示す。また、この結果、図6の(1)や(2)において、血流量があがっている箇所は、心臓の鼓動により発生した脈波が到達したことを示す。   An example of a waveform indicated by an electrical signal converted by the optical sensor 201 or the optical sensor 202 will be described with reference to FIG. FIG. 6 is a diagram illustrating an example of a waveform indicated by an electrical signal converted by the optical sensor. Here, (1) of FIG. 6 shows an example of a waveform indicated by the electrical signal converted by the optical sensor 201. FIG. In other words, (1) in FIG. 6 is the waveform for the optical sensor in contact with the ear. FIG. 6 (2) shows an example of a waveform indicated by the electrical signal converted by the optical sensor 202. FIG. In other words, (2) of FIG. 6 is a waveform about the optical sensor in contact with the finger. The horizontal axes of (1) and (2) in FIG. 6 indicate the time axis, and the vertical axes of (1) and (2) in FIG. 6 indicate the amount of light received by the light receiving element. That is, the vertical axis indicates the blood flow volume at the location where the light sensor irradiates light. As a result, in (1) and (2) of FIG. 6, the portion where the blood flow rate is increased indicates that the pulse wave generated by the heartbeat has arrived.

図2の説明に戻る。出力部203は、制御部400と接続される。出力部203は、制御部400から情報を受け付け、受け付けた情報を出力する。出力部203は、例えば、モニタやディスプレイ、タッチパネル、又は、スピーカなどが該当する。なお、出力部203によって出力される情報の詳細については、ここでは説明を省略し、関係する各部について説明する際に併せて説明する。   Returning to the description of FIG. The output unit 203 is connected to the control unit 400. The output unit 203 receives information from the control unit 400 and outputs the received information. The output unit 203 corresponds to, for example, a monitor, a display, a touch panel, or a speaker. The details of the information output by the output unit 203 will be omitted here, and will be described together with the description of each related unit.

記憶部300は、制御部400と接続される。また、記憶部300は、制御部400による各種処理に用いるデータを記憶する。記憶部300は、例えば、RAM(Random Access Memory)やROM(Read Only Memory)、フラッシュメモリ(Flash Memory)などの半導体メモリ素子、又は、ハードディスクや光ディスクなどの記憶装置である。記憶部300は、図2に示す例では、閾値テーブル301を有する。   Storage unit 300 is connected to control unit 400. The storage unit 300 stores data used for various processes performed by the control unit 400. The storage unit 300 is, for example, a semiconductor memory device such as a random access memory (RAM), a read only memory (ROM), or a flash memory, or a storage device such as a hard disk or an optical disk. The storage unit 300 includes a threshold value table 301 in the example illustrated in FIG.

閾値テーブル301は、後述するように、血管状態推定部406によって用いられる閾値を記憶する。例えば、図7を用いて、閾値テーブル301に記憶された情報の一例について説明する。図7は、実施例2における閾値テーブルに記憶された情報の一例について説明する図である。図7に示す例では、血管の状態が相当する年代を示す「年代」を記憶する場合を示した。また、図7の「閾値」に含まれる「s」は、後述するように、血管状態推定部406によって算出される値である。図7の「s」の詳細については、血管状態推定部406について説明する際に併せて説明するため、ここでは説明を省略する。   The threshold value table 301 stores threshold values used by the blood vessel state estimation unit 406, as will be described later. For example, an example of information stored in the threshold value table 301 will be described with reference to FIG. FIG. 7 is a diagram illustrating an example of information stored in the threshold value table according to the second embodiment. In the example shown in FIG. 7, the case where “age” indicating the age corresponding to the state of the blood vessel is stored is shown. Further, “s” included in the “threshold value” in FIG. 7 is a value calculated by the blood vessel state estimation unit 406 as described later. Since the details of “s” in FIG. 7 will be described together with the description of the blood vessel state estimation unit 406, description thereof is omitted here.

図7に示すように、閾値テーブル301は、年代に対応付けて閾値を記憶する。例えば、図7に示す例では、閾値テーブル301は、年代「20代」に対応付けて閾値「0<s<2」を記憶し、年代「30代」に対応付けて閾値「2<s<4」を記憶する。また、図7に示す例では、閾値テーブル301は、年代「40代」や「50代」についても閾値を記憶する。すなわち、例えば、閾値テーブル301は、血管状態推定部406によって「s=1」であると算出された場合には、血管の状態が「20代」に相当することを記憶する。   As shown in FIG. 7, the threshold value table 301 stores threshold values in association with ages. For example, in the example illustrated in FIG. 7, the threshold table 301 stores a threshold “0 <s <2” associated with the age “20s” and a threshold “2 <s <2” associated with the age “30s”. 4 "is stored. Further, in the example illustrated in FIG. 7, the threshold value table 301 stores threshold values for the ages “40s” and “50s”. That is, for example, when the blood vessel state estimation unit 406 calculates “s = 1”, the threshold value table 301 stores that the blood vessel state corresponds to “20's”.

なお、以下では、図7に示すように、閾値テーブル301は、年代ごとに閾値を記憶する場合を例に説明するが、本発明はこれに限定されるものではない。例えば、閾値テーブル301は、年齢ごとに閾値を記憶しても良く、利用者が任意に設定した範囲ごとに閾値を記憶して良い。   In the following, as illustrated in FIG. 7, the threshold value table 301 is described as an example in which threshold values are stored for each age, but the present invention is not limited to this. For example, the threshold value table 301 may store a threshold value for each age, or may store a threshold value for each range arbitrarily set by the user.

制御部400は、光センサ201、光センサ202、出力部203及び記憶部300と接続される。制御部400は、各種の処理手順などを規定したプログラムを記憶する内部メモリを有し、種々の処理を制御する。制御部400は、例えば、ASIC(Application Specific Integrated Circuit)、FPGA(Field Programmable Gate Array)、CPU(Central Processing Unit)、MPU(Micro Processing Unit)などの電子回路である。制御部400は、図2に示す例では、脈波間隔算出部401と、脈波間隔算出部402と、n次微分処理部403と、特徴点識別部404と、差分算出部405と、血管状態推定部406とを有する。   The control unit 400 is connected to the optical sensor 201, the optical sensor 202, the output unit 203, and the storage unit 300. The control unit 400 includes an internal memory that stores a program that defines various processing procedures and the like, and controls various processes. The control unit 400 is an electronic circuit such as an application specific integrated circuit (ASIC), a field programmable gate array (FPGA), a central processing unit (CPU), or a micro processing unit (MPU). In the example illustrated in FIG. 2, the control unit 400 includes a pulse wave interval calculation unit 401, a pulse wave interval calculation unit 402, an nth-order differentiation processing unit 403, a feature point identification unit 404, a difference calculation unit 405, and a blood vessel. A state estimation unit 406.

なお、図2に示す例では、脈波間隔算出部401と脈波間隔算出部402とを別の制御部として記載したが、本発明はこれに限定されるものではなく、統合して1つの制御部としても良い。また、脈波間隔算出部401と脈波間隔算出部402とは、それぞれ、光センサ201や光センサ202と統合しても良い。   In the example shown in FIG. 2, the pulse wave interval calculation unit 401 and the pulse wave interval calculation unit 402 are described as separate control units. However, the present invention is not limited to this, and is integrated into one unit. It is good also as a control part. Moreover, the pulse wave interval calculation unit 401 and the pulse wave interval calculation unit 402 may be integrated with the optical sensor 201 and the optical sensor 202, respectively.

脈波間隔算出部401は、光センサ201によって変換された電気信号が示す波形から、連続する2つの脈波の間の脈波間隔を算出する。具体的には、脈波間隔算出部401は、電気信号が示す波形のうち最大点を識別し、隣接する最大点間の間隔を算出することで脈波間隔を算出する。また、脈波間隔算出部402は、脈波間隔算出部401と同様に、光センサ202によって変換された電気信号が示す波形から、脈波間隔を算出する。   The pulse wave interval calculation unit 401 calculates a pulse wave interval between two consecutive pulse waves from the waveform indicated by the electrical signal converted by the optical sensor 201. Specifically, the pulse wave interval calculation unit 401 identifies the maximum point in the waveform indicated by the electrical signal, and calculates the pulse wave interval by calculating the interval between adjacent maximum points. Similarly to the pulse wave interval calculation unit 401, the pulse wave interval calculation unit 402 calculates the pulse wave interval from the waveform indicated by the electrical signal converted by the optical sensor 202.

また、脈波間隔算出部401や脈波間隔算出部402は、それぞれ、脈波間隔を用いて脈飛びがあった箇所を識別する。脈飛びがあった箇所を識別する処理の一例について、図8を用いて説明する。図8は、実施例2における脈波間隔算出部による脈飛びがあった箇所を識別する処理の一例について説明する図である。図8の(1)は、光センサ201についての波形を示し、図6の(1)に対応する。つまり、図8の(1)は、耳と接触した光センサについての波形である。図8の(2)は、光センサ202についての波形を示し、図6の(2)に対応する。つまり、図8の(2)は、指と接触した光センサについての波形である。図8の(1)や(2)の横軸は、時間軸を示し、図8の(1)や(2)の縦軸は、血流量を示す。また、図8の(1)の横軸と図8の(2)の横軸とは、それぞれ対応するものとして説明する。言い換えると、横軸が同じ位置である場合には、図8の(1)と(2)とが同じ時間におけるデータを示すものとして説明する。また、図8の「脈波n」(nは自然数)は、波形に含まれる脈波それぞれを示し、図8の「最大点」は、脈波それぞれに含まれる最大点を示す。また、図8の「脈波間隔」は、脈波間隔算出部401や脈波間隔算出部402によって算出される脈波間隔を示す。   Also, the pulse wave interval calculation unit 401 and the pulse wave interval calculation unit 402 each identify a place where a pulse jump has occurred using the pulse wave interval. An example of the process for identifying the place where the pulse jumps will be described with reference to FIG. FIG. 8 is a diagram illustrating an example of a process for identifying a place where a pulse jump has occurred by the pulse wave interval calculation unit according to the second embodiment. (1) in FIG. 8 shows a waveform for the optical sensor 201 and corresponds to (1) in FIG. That is, (1) of FIG. 8 is a waveform about the optical sensor in contact with the ear. (2) in FIG. 8 shows a waveform for the optical sensor 202 and corresponds to (2) in FIG. That is, (2) of FIG. 8 is a waveform about the optical sensor in contact with the finger. The horizontal axes of (1) and (2) in FIG. 8 indicate the time axis, and the vertical axis of (1) and (2) in FIG. 8 indicates the blood flow. Further, description will be made assuming that the horizontal axis of (1) in FIG. 8 and the horizontal axis of (2) in FIG. 8 correspond to each other. In other words, when the horizontal axis is the same position, description will be made assuming that (1) and (2) in FIG. 8 indicate data at the same time. Further, “pulse wave n” (n is a natural number) in FIG. 8 indicates each pulse wave included in the waveform, and “maximum point” in FIG. 8 indicates the maximum point included in each pulse wave. 8 indicates a pulse wave interval calculated by the pulse wave interval calculation unit 401 or the pulse wave interval calculation unit 402.

ここで、図8の(1)に示す例では、脈波1〜6間における連続する2つの脈波間の脈波間隔には、他の脈波間隔のn倍(nは自然数)となる脈波間隔はない。言い換えると、図8の(1)に示す例では、脈波間隔は、どれも同じような値になる。この場合、脈波間隔算出部401は、脈飛びがないと判定する。一方、図8の(2)に示す例では、脈波1〜2間の脈波間隔は、脈波2〜5間における連続する2つの脈波間の脈波間隔と比較して2倍の値になる。この場合、脈波間隔算出部401は、脈波1〜2間において脈飛びがあったと識別する。なお、脈波間隔算出部401による処理結果は、後述するように、特徴点識別部404によって用いられる。   Here, in the example shown in (1) of FIG. 8, the pulse wave interval between two continuous pulse waves between pulse waves 1 to 6 is a pulse that is n times (n is a natural number) the other pulse wave intervals. There is no wave interval. In other words, in the example shown in (1) of FIG. 8, all pulse wave intervals have the same value. In this case, the pulse wave interval calculation unit 401 determines that there is no pulse skip. On the other hand, in the example shown in (2) of FIG. 8, the pulse wave interval between the pulse waves 1 and 2 is twice the pulse wave interval between two continuous pulse waves between the pulse waves 2 to 5. become. In this case, the pulse wave interval calculation unit 401 identifies that there is a pulse jump between the pulse waves 1 and 2. Note that the processing result by the pulse wave interval calculation unit 401 is used by the feature point identification unit 404 as described later.

n次微分処理部403は、光センサごとに、電気信号が示す波形に対してn次微分(nは自然数)を実行する。具体的には、n次微分処理部403は、光センサ201についての波形に対してn次微分(nは自然数)を実行し、また、光センサ202についての波形に対してn次微分(nは自然数)を実行する。以下では、特に言及しない限り、n次微分処理部403は、2次微分を実行する場合を例に説明するが、本発明はこれに限定されるものではない。例えば、n次微分処理部403は、1次微分や3次微分を実行しても良く、利用者が任意の微分回数を設定して良い。   The n-th order differential processing unit 403 performs n-th order differentiation (n is a natural number) for the waveform indicated by the electrical signal for each optical sensor. Specifically, the nth-order differentiation processing unit 403 performs nth-order differentiation (n is a natural number) on the waveform for the optical sensor 201, and nth-order differentiation (n for the waveform for the optical sensor 202). Is a natural number). Hereinafter, the n-th order differentiation processing unit 403 will be described as an example where the second-order differentiation is performed unless otherwise specified, but the present invention is not limited to this. For example, the nth-order differentiation processing unit 403 may execute first-order differentiation or third-order differentiation, and the user may set an arbitrary number of differentiations.

例えば、n次微分処理部403は、2次微分を実行した場合には、図9Aや図9Bに示すような2次微分波形が得られる。ここで、図9Aと図9Bとは、実施例2における2次微分波形の一例について説明する図である。図9Aの(1)は、光センサ201についての波形を示し、図6の(1)や図8の(1)に対応する。つまり、図9Aの(1)は、耳と接触した光センサについての波形である。図9Aの(2)は、図9Aの(1)に対して2次微分することで得られた2次微分波形である。つまり、図9Aの(2)は、耳についての2次微分波形である。図9Bの(1)は、光センサ202についての波形を示し、図6の(2)や図8の(2)に対応する。つまり、図9Bの(1)は、指と接触した光センサについての波形である。図9Bの(2)は、図9Bの(1)に対して2次微分することで得られた2次微分波形である。つまり、図9Bの(2)は、指についての2次微分波形である。   For example, the n-th order differential processing unit 403 obtains a second-order differential waveform as shown in FIGS. 9A and 9B when the second-order differentiation is executed. Here, FIG. 9A and FIG. 9B are diagrams for explaining an example of the secondary differential waveform in the second embodiment. 9A shows the waveform of the optical sensor 201, and corresponds to (1) in FIG. 6 and (1) in FIG. That is, (1) of FIG. 9A is a waveform about the optical sensor in contact with the ear. (2) in FIG. 9A is a second-order differential waveform obtained by second-order differentiation with respect to (1) in FIG. 9A. That is, (2) in FIG. 9A is a second-order differential waveform for the ear. (1) in FIG. 9B shows a waveform for the optical sensor 202, and corresponds to (2) in FIG. 6 and (2) in FIG. That is, (1) of FIG. 9B is a waveform about the optical sensor in contact with the finger. (2) in FIG. 9B is a second-order differential waveform obtained by second-order differentiation with respect to (1) in FIG. 9B. That is, (2) in FIG. 9B is a second-order differential waveform for the finger.

また、図9Aや図9Bの横軸は、時間軸を示す。また、図9Aの(1)や図9Bの(1)の縦軸は、血流量を示す。図9の(1)の横軸と図9の(2)の横軸とは、それぞれ対応するものとして説明する。また、図9Aの(2)や図9Bの(2)の縦軸は、2次微分により得られた値を示す。   Moreover, the horizontal axis of FIG. 9A and FIG. 9B shows a time axis. Moreover, the vertical axis | shaft of (1) of FIG. 9A and (1) of FIG. 9B shows a blood flow rate. The description will be made assuming that the horizontal axis of (1) in FIG. 9 and the horizontal axis of (2) in FIG. 9 correspond to each other. Moreover, the vertical axis | shaft of (2) of FIG. 9A and (2) of FIG. 9B shows the value obtained by the secondary differentiation.

特徴点識別部404は、光センサごとに、2次微分波形に含まれる特徴点を識別する。具体的には、特徴点識別部404は、光センサ201に対応する2次微分波形に含まれる脈波ごとに特徴点を複数識別する。また、特徴点識別部404は、光センサ202に対応する2次微分波形に含まれる脈波ごとに特徴点を複数識別する。   The feature point identification unit 404 identifies the feature points included in the secondary differential waveform for each optical sensor. Specifically, the feature point identifying unit 404 identifies a plurality of feature points for each pulse wave included in the secondary differential waveform corresponding to the optical sensor 201. The feature point identifying unit 404 identifies a plurality of feature points for each pulse wave included in the secondary differential waveform corresponding to the optical sensor 202.

また、特徴点識別部404は、脈波間隔算出部401による処理結果に基づいて、両方の2次微分波形に脈波が含まれるタイミングを識別する。そして、特徴点識別部404は、識別したタイミングに含まれる脈波ごとに、少なくとも4つの特徴点を識別する。つまり、特徴点識別部404は、2次微分波形に含まれる脈波ごとに少なくとも4つの特徴点を識別する。   Further, the feature point identifying unit 404 identifies the timing at which the pulse wave is included in both secondary differential waveforms based on the processing result by the pulse wave interval calculating unit 401. The feature point identification unit 404 identifies at least four feature points for each pulse wave included in the identified timing. That is, the feature point identifying unit 404 identifies at least four feature points for each pulse wave included in the secondary differential waveform.

図10を用いて、特徴点識別部404による特徴点識別処理について更に説明する。図10は、実施例2における特徴点識別部による特徴点識別処理について更に説明する図である。なお、図10の(1)は、耳についての2次微分波形であり、図9Aの(2)に対応する。また、図10の(2)は、指についての2次微分波形であり、図9Bの(2)に対応する。図10の(1)に示すように、耳についての2次微分波形には脈飛びがなく、図10の(2)に示すように、指についての2次微分波形には脈波1と脈波2との間に脈飛びがある場合を例に説明する。   The feature point identifying process by the feature point identifying unit 404 will be further described with reference to FIG. FIG. 10 is a diagram for further explaining the feature point identifying process by the feature point identifying unit in the second embodiment. In addition, (1) of FIG. 10 is a secondary differential waveform about the ear, and corresponds to (2) of FIG. 9A. Further, (2) in FIG. 10 is a second-order differential waveform for the finger, and corresponds to (2) in FIG. 9B. As shown in (1) of FIG. 10, the secondary differential waveform for the ear has no pulse skipping, and as shown in (2) of FIG. 10, the secondary differential waveform for the finger has pulse wave 1 and pulse. A case where there is a pulse jump between the waves 2 will be described as an example.

また、図10の(1)や(2)の横軸は、時間軸を示し、図10の(1)や(2)の縦軸は、2次微分処理により得られた値を示す。図10の(1)の横軸と図10の(2)の横軸とは、それぞれ対応するものとして説明する。図10の「脈波n」(nは自然数)は、波形に含まれる脈波それぞれを示し、図10の「特徴点a〜e」は、特徴点識別部404によって脈波ごとに識別された特徴点を示す。また、図10の「脈飛びあり」に示すように、脈飛びがあった箇所を示した。なお、特徴点「a」〜「e」は、それぞれ、1番特徴点〜5番特徴点とも称する。なお、1番特徴点は、脈波に含まれる特徴点のうち出現時間が1番目になる特徴点を示す。また、同様に、2番特徴点〜5番特徴点は、それぞれ、脈波に含まれる特徴点のうち出現時間が2番目〜5番目になる特徴点を示す。   Also, the horizontal axes of (1) and (2) in FIG. 10 indicate the time axis, and the vertical axes of (1) and (2) in FIG. 10 indicate values obtained by the secondary differentiation process. The description will be made assuming that the horizontal axis of (1) in FIG. 10 and the horizontal axis of (2) in FIG. 10 correspond to each other. “Pulse wave n” (n is a natural number) in FIG. 10 indicates each pulse wave included in the waveform, and “feature points a to e” in FIG. 10 are identified for each pulse wave by the feature point identifying unit 404. Indicates feature points. In addition, as shown in “with pulse skip” in FIG. Note that the feature points “a” to “e” are also referred to as the first feature point to the fifth feature point, respectively. The first feature point is a feature point having the first appearance time among the feature points included in the pulse wave. Similarly, the second feature point to the fifth feature point indicate feature points whose appearance times are the second to fifth among the feature points included in the pulse wave, respectively.

この場合、特徴点識別部404は、図10の(1)の脈拍1、3〜6と同じタイミングにおいて、図10の(2)に示す2次微分波形にも脈波が含まれると識別する。この結果、特徴点識別部404は、図10の(1)の脈拍1、3〜6について特徴点を識別する。また、特徴点識別部404は、脈拍1、3〜6と同じタイミングにある脈波について特徴点を識別する。例えば、特徴点識別部404は、図10の(1)の脈拍1〜5について特徴点を識別する。なお、図10に示すように、図10の(1)の脈拍1、3〜6は、それぞれ、図10の(2)の脈拍1〜5と同じタイミングになる。   In this case, the feature point identification unit 404 identifies that the pulse wave is also included in the secondary differential waveform shown in (2) of FIG. 10 at the same timing as the pulses 1, 3 to 6 in (1) of FIG. . As a result, the feature point identifying unit 404 identifies feature points for the pulses 1, 3 to 6 in (1) of FIG. Further, the feature point identifying unit 404 identifies feature points for the pulse wave at the same timing as the pulses 1 and 3 to 6. For example, the feature point identifying unit 404 identifies feature points for the pulses 1 to 5 in (1) of FIG. As shown in FIG. 10, pulses 1 and 3 to 6 in (1) in FIG. 10 have the same timing as pulses 1 to 5 in (2) in FIG.

また、図10に示すように、特徴点識別部404は、脈拍ごとに少なくても4つの特徴点を識別する。例えば、特徴点識別部404は、特徴点a〜eあるいは特徴点a〜dを脈拍ごとに識別する。なお、図10の(1)や(2)に示すように、特徴点a〜dについては、ほとんどの2次微分波形において識別可能である。図10の(1)に示すように、特徴点eについては、識別できないことが多々あることが知られている。   As shown in FIG. 10, the feature point identifying unit 404 identifies at least four feature points for each pulse. For example, the feature point identifying unit 404 identifies feature points a to e or feature points a to d for each pulse. As shown in (1) and (2) of FIG. 10, the feature points a to d can be identified in most secondary differential waveforms. As shown in (1) of FIG. 10, it is known that there are many cases where the feature point e cannot be identified.

なお、図10に示す例では、図10の(1)に示すように、耳についての2次微分波形において、特徴点識別部404が特徴点eを識別できなかった場合を例に示した。また、図10の(2)に示すように、指についての2次微分波形において、特徴点識別部404が特徴点eを識別できた場合を例に示した。しかしながら、本発明はこれに限定されるものではなく、耳についての2次微分波形において特徴点eを識別できた場合であっても良く、指についての2次微分波形において特徴点eを識別できなかった場合であっても良い。   In the example illustrated in FIG. 10, as illustrated in FIG. 10 (1), the case where the feature point identification unit 404 cannot identify the feature point e in the second-order differential waveform for the ear is illustrated as an example. Further, as shown in (2) of FIG. 10, the case where the feature point identifying unit 404 can identify the feature point e in the secondary differential waveform for the finger is shown as an example. However, the present invention is not limited to this, and the feature point e may be identified in the secondary differential waveform for the ear, and the feature point e can be identified in the secondary differential waveform for the finger. There may be no case.

差分算出部405は、光センサごとに、特徴点識別部404によって識別された特徴点の出現時間の差分を算出する。具体的には、差分算出部405は、光センサ201に対応する2次微分波形の特徴点について差分を算出する。また、差分算出部405は、光センサ202に対応する2次微分波形の特徴点について差分を算出する。ここで、出現時間の差分は、時間が該当する。例えば、出現時間の差分は、1秒や5m秒などの時間が該当する。   The difference calculation unit 405 calculates the difference between the appearance times of the feature points identified by the feature point identification unit 404 for each optical sensor. Specifically, the difference calculation unit 405 calculates a difference for the feature point of the secondary differential waveform corresponding to the optical sensor 201. Further, the difference calculation unit 405 calculates a difference for the feature point of the secondary differential waveform corresponding to the optical sensor 202. Here, the difference in appearance time corresponds to time. For example, the difference in appearance time corresponds to a time such as 1 second or 5 milliseconds.

例えば、差分算出部405は、他方の光センサと比較して太い血管に光を照射する光センサについて、脈波ごとに、3番目の特徴点と4番目の特徴点との間の第1の差分を算出する。また、差分算出部405は、他方の光センサと比較して細い血管に光を照射する光センサについて、脈波ごとに、1番目の特徴点と2番目の特徴点との間の第2の差分を算出する。   For example, for the optical sensor that irradiates light on a blood vessel that is thicker than the other optical sensor, the difference calculation unit 405 includes a first feature point between the third feature point and the fourth feature point for each pulse wave. Calculate the difference. In addition, the difference calculation unit 405 has a second difference between the first feature point and the second feature point for each pulse wave with respect to the optical sensor that irradiates light to a thin blood vessel as compared with the other optical sensor. Calculate the difference.

図11を用いて更に説明する。図11は、実施例2における差分算出部による出現時間の差分算出処理の一例について説明する図である。図11の(1)は、耳についての2次微分波形であり、図9Aの(2)や図10の(1)に対応する。また、図11の(2)は、指についての2次微分波形であり、図9Bの(2)や図10の(2)に対応する。また、図11の(1)や(2)の横軸は、時間軸を示し、図11の(1)や(2)の縦軸は、2次微分処理により得られた値を示す。図11の(1)の横軸と図11の(2)の横軸とは、それぞれ対応するものとして説明する。図11の「脈波n」(nは自然数)は、波形に含まれる脈波それぞれを示す。図11の「特徴点a〜d」は、特徴点識別部404によって識別された特徴点を示す。   This will be further described with reference to FIG. FIG. 11 is a diagram illustrating an example of a difference calculation process of appearance time by a difference calculation unit according to the second embodiment. (1) in FIG. 11 is a second-order differential waveform for the ear, and corresponds to (2) in FIG. 9A and (1) in FIG. Further, (2) in FIG. 11 is a second-order differential waveform for the finger, and corresponds to (2) in FIG. 9B and (2) in FIG. Also, the horizontal axes of (1) and (2) in FIG. 11 indicate the time axis, and the vertical axes of (1) and (2) in FIG. 11 indicate values obtained by the secondary differentiation process. The horizontal axis of (1) in FIG. 11 and the horizontal axis of (2) in FIG. 11 are assumed to correspond to each other. “Pulse wave n” (n is a natural number) in FIG. 11 indicates each pulse wave included in the waveform. “Feature points a to d” in FIG. 11 indicate feature points identified by the feature point identifying unit 404.

この場合、差分算出部405は、図11の(1)の「耳差分1」〜「耳差分5」に示すように、耳についての2次微分波形に含まれる脈拍1、3〜6について、特徴点aが出現する時間と、特徴点bが出現する時間との間の差分の絶対値を算出する。また、差分算出部405は、図11の(2)の「指差分1」〜「指差分5」に示すように、指についての2次微分波形に含まれる脈拍1〜5について、特徴点cが出現する時間と、特徴点dが出現する時間との間の差分の絶対値を算出する。   In this case, as shown in “ear difference 1” to “ear difference 5” in (1) of FIG. 11, the difference calculation unit 405 performs the following operations on the pulses 1 and 3 to 6 included in the secondary differential waveform for the ear. The absolute value of the difference between the time when the feature point a appears and the time when the feature point b appears is calculated. Moreover, the difference calculation part 405 is the feature point c about the pulses 1-5 contained in the secondary differential waveform about a finger | toe, as shown to "finger difference 1"-"finger difference 5" of (2) of FIG. The absolute value of the difference between the time when appears and the time when the feature point d appears is calculated.

なお、以下では、耳についての2次微分波形と指についての2次微分波形とについて、それぞれ、特徴点c−d間の出現時間の差分と特徴点a−b間の出現時間の差分とを算出する場合について説明するが、本発明はこれに限定されるものではない。例えば、差分算出部405は、耳についての2次微分波形について、特徴点a−b間の出現時間の差分や特徴点b−c間の出現時間の差分を算出しても良い。また、例えば、差分算出部405は、指についての2次微分波形について、特徴点b−c間の出現時間の差分や特徴点c−d間の出現時間の差分を算出しても良い。また、差分算出部405は、耳についての2次微分波形と指についての2次微分波形との両方において、同じ特徴点間の差分を算出しても良い。   In the following, regarding the secondary differential waveform for the ear and the secondary differential waveform for the finger, the difference in the appearance time between the feature points cd and the difference in the appearance time between the feature points ab, respectively. Although the case where it calculates is demonstrated, this invention is not limited to this. For example, the difference calculation unit 405 may calculate a difference in appearance time between feature points ab and a difference in appearance time between feature points bc for the second-order differential waveform for the ear. Further, for example, the difference calculation unit 405 may calculate the difference in appearance time between feature points bc and the difference in appearance time between feature points cd for the secondary differential waveform for the finger. Further, the difference calculation unit 405 may calculate a difference between the same feature points in both the second-order differential waveform for the ear and the second-order differential waveform for the finger.

血管状態推定部406は、差分算出部405によって光センサごとに算出された差分に基づいて人体の血管の状態を推定し、推定結果を出力する。具体的には、血管状態推定部406は、光センサごとに算出された差分間の差分であるセンサ間差分を、差分算出部405が算出した複数の差分ごとに算出する。例えば、血管状態推定部406は、第1の差分と第2の差分との間の差分を脈波ごとに算出することでセンサ間差分各々を算出する。図10に示す例を用いて説明すると、血管状態推定部406は、耳差分1〜5と、指差分1〜5との間の差分をそれぞれ算出することで、センサ間差分1〜5を算出する。   The blood vessel state estimation unit 406 estimates the blood vessel state of the human body based on the difference calculated for each optical sensor by the difference calculation unit 405, and outputs an estimation result. Specifically, the blood vessel state estimation unit 406 calculates an inter-sensor difference, which is a difference between differences calculated for each optical sensor, for each of a plurality of differences calculated by the difference calculation unit 405. For example, the blood vessel state estimation unit 406 calculates each inter-sensor difference by calculating a difference between the first difference and the second difference for each pulse wave. If it demonstrates using the example shown in FIG. 10, the blood-vessel state estimation part 406 will calculate the differences 1-5 between sensors by calculating the difference between the ear differences 1-5 and the finger differences 1-5, respectively. To do.

なお、血管状態推定部406は、5つのセンサ間差分を算出する場合を例に説明したが、本発明はこれに限定されるものではない。例えば、4つ以下のセンサ間差分を算出しても良く、6つ以上のセンサ間差分を算出しても良く、利用者が任意に設定して良い。   In addition, although the blood vessel state estimation unit 406 has been described as an example of calculating the difference between the five sensors, the present invention is not limited to this. For example, four or less inter-sensor differences may be calculated, six or more inter-sensor differences may be calculated, and the user may arbitrarily set them.

また、血管状態推定部406は、センサ間差分各々が広く分散した場合には狭く分散した場合よりも血管年齢が高いと推定し、センサ間差分各々が狭く分散した場合には広く分散した場合よりも血管年齢が低いと推定する。なお、血管年齢は、血管の状態が相当する年齢を示す。つまり、血管年齢が「20代」であると推定された場合には、血管の状態が、20代の人の血管の状態と同等であることを示す。   In addition, the blood vessel state estimation unit 406 estimates that the age of the blood vessel is higher than the case where the difference between the sensors is widely dispersed than when the difference is narrowly distributed, and the case where the difference between the sensors is narrowly dispersed than the case where the differences are widely dispersed. Also assume that the vascular age is low. The blood vessel age indicates an age corresponding to a blood vessel state. That is, when it is estimated that the blood vessel age is “20s”, it indicates that the blood vessel state is equivalent to the blood vessel state of a person in their 20s.

図12A〜図12Cを用いて、センサ間差分の分散と年齢との関係について説明する。なお、図12A〜図12Cは、センサ間差分の分散と年齢との関係について説明するための図である。図12Aは、30歳代の利用者のデータから算出された指差分と耳差分との関係を示した図である。図12Bは、40歳代の利用者のデータから算出された指差分と耳差分との関係を示した図である。図12Cは、50歳代の利用者のデータから算出された指差分と耳差分との関係を示した図である。図12A〜図12Cにおいて、横軸は、脈拍数を示し、縦軸は、指差分又は耳差分の絶対値を示す。   The relationship between the variance of the difference between sensors and the age will be described with reference to FIGS. 12A to 12C. 12A to 12C are diagrams for explaining the relationship between the variance of the difference between sensors and the age. FIG. 12A is a diagram illustrating a relationship between finger differences and ear differences calculated from data of users in their thirties. FIG. 12B is a diagram showing a relationship between finger differences and ear differences calculated from data of users in their 40s. FIG. 12C is a diagram illustrating a relationship between finger differences and ear differences calculated from data of users in their 50s. 12A to 12C, the horizontal axis indicates the pulse rate, and the vertical axis indicates the absolute value of the finger difference or the ear difference.

図12Aに示すように、30歳の利用者のデータでは、指差分と耳差分とが全体的に同じように変化した。つまり、図12Aに示す例では、指差分が増加した場合には耳差分も増加し、指差分が減少した場合には耳差分も減少した。一方、図12Bに示すように、43歳の利用者のデータでは、指差分と耳差分とが同じように変化したり、同じように変化しなかったりした。つまり、図12Bの601に示す例では、指差分が増加した場合には耳差分も増加し、指差分が減少した場合には耳差分も減少した。言い換えると、図12Bの601に示す例では、指差分と耳差分とが同じように変化した。これに対して、図12Bの602に示す例では、指差分が増加した場合には耳差分が減少し、指差分が減少した場合には耳差分は増加した。言い換えると、図12Bの602に示す例では、指差分と耳差分とがばらばらに変化した。また、図12Cに示すように、50歳の利用者のデータでは、指差分と耳差分とが全体的にばらばらに変化した。   As shown in FIG. 12A, in the data of a 30-year-old user, the finger difference and the ear difference changed in the same manner as a whole. That is, in the example shown in FIG. 12A, when the finger difference increases, the ear difference also increases, and when the finger difference decreases, the ear difference also decreases. On the other hand, as shown in FIG. 12B, in the data of a 43-year-old user, the finger difference and the ear difference changed in the same way or did not change in the same way. That is, in the example indicated by 601 in FIG. 12B, the ear difference increases when the finger difference increases, and the ear difference decreases when the finger difference decreases. In other words, in the example shown at 601 in FIG. 12B, the finger difference and the ear difference changed in the same way. On the other hand, in the example shown by 602 in FIG. 12B, the ear difference decreases when the finger difference increases, and the ear difference increases when the finger difference decreases. In other words, in the example shown by 602 in FIG. 12B, the finger difference and the ear difference changed in a discontinuous manner. Further, as shown in FIG. 12C, the finger difference and the ear difference changed in a whole manner in the data of a 50-year-old user.

つまり、年齢が低い場合には、センサ間差分がどれも同じような値になる結果、複数のセンサ間差分の分散の程度は、年齢が高い場合と比較して低くなると考えられる。一方、年齢が高い場合には、センサ間差分がばらばらな値になる結果、複数のセンサ間差分の分散の程度は、年齢が低い場合と比較して高くなると考えられる。   That is, when the age is low, all the differences between the sensors have the same value, and as a result, the degree of dispersion of the differences among the plurality of sensors is considered to be lower than when the age is high. On the other hand, when the age is high, the difference between the sensors becomes a different value. As a result, the degree of dispersion of the differences between the plurality of sensors is considered to be higher than that when the age is low.

図13を用いて、センサ間差分の分散と年齢との関係について更に説明する。図13は、センサ間差分の分散と年齢との関係について説明するための図である。図13の横軸は、センサ間差分の値を示し、図13の縦軸は、対応するセンサ間差分の値が得られる確率を示す。つまり、縦軸の値が高ければ高いほど、対応するセンサ間差分の値が得られる確率が高いことを示す。言い換えると、血管状態推定部406が複数のセンサ間差分を算出した場合に、同じ値のセンサ間差分が算出される確率が高いことを示す。   The relationship between the variance of the difference between sensors and the age will be further described with reference to FIG. FIG. 13 is a diagram for explaining the relationship between the variance of differences between sensors and the age. The horizontal axis of FIG. 13 shows the value of the inter-sensor difference, and the vertical axis of FIG. 13 shows the probability that the corresponding inter-sensor difference value is obtained. That is, the higher the value on the vertical axis, the higher the probability that a corresponding inter-sensor difference value is obtained. In other words, when the blood vessel state estimation unit 406 calculates a plurality of inter-sensor differences, it indicates that there is a high probability that the same inter-sensor difference is calculated.

図13に示すように、20代の場合には、他の年代と比較してピークの縦軸の値が最も高く、他の年代と比較してピークの横軸の幅が最も小さくなる。また、ピークの縦軸の値は、20代から50代へと加齢に従って値が小さくなり、ピークの横軸の幅は、20代から50代へと加齢に従って値が大きくなる。   As shown in FIG. 13, in the twenties, the value of the vertical axis of the peak is the highest compared to other ages, and the width of the horizontal axis of the peak is the smallest compared to other ages. Further, the value of the vertical axis of the peak decreases with age from the 20s to the 50s, and the width of the horizontal axis of the peak increases with age from the 20s to the 50s.

このように、センサ間差分の分散と年齢との間には、年齢が高くなるに従って、センサ間差分の分散の程度が大きくなるという関係がある。このことを踏まえ、血管状態推定部406は、血管の状態の推定し、推定結果を出力する。具体的には、利用者によって、センサ間差分の標準偏差の範囲が年代ごとに予め算出され、閾値テーブル301に格納される。その上で、血管状態推定部406は、センサ間差分各々の標準偏差を算出し、算出した標準偏差の値がいずれの閾値に該当するかを判定し、該当すると判定した閾値に対応する年齢を推定結果とする。なお、図7の「s」は、センサ間差分の標準偏差を示す。   As described above, there is a relationship between the variance of the difference between sensors and the age that the degree of variance of the difference between sensors increases as the age increases. Based on this, the blood vessel state estimation unit 406 estimates the blood vessel state and outputs an estimation result. Specifically, the standard deviation range of the inter-sensor difference is calculated in advance for each age by the user and stored in the threshold value table 301. Then, the blood vessel state estimation unit 406 calculates a standard deviation of each inter-sensor difference, determines which threshold value the calculated standard deviation value corresponds to, and determines an age corresponding to the threshold value determined to be applicable. Estimated result. Note that “s” in FIG. 7 indicates the standard deviation of the inter-sensor difference.

なお、例えば、血管状態推定部406は、下記の(数1)や(数2)を用いることで、標準偏差「s」を算出する。   For example, the blood vessel state estimation unit 406 calculates the standard deviation “s” by using the following (Equation 1) and (Equation 2).

Figure 2011072674
Figure 2011072674

Figure 2011072674
Figure 2011072674

なお、(数1)は、「n」個算出したセンサ間差分の平均値を算出する式を示す。(数1)の「x」は、センサ間差分の値を示す。(数2)に示すように、標準偏差「s」は、(数1)にて算出した平均値を用いて算出される。   Note that (Equation 1) represents an equation for calculating an average value of “n” calculated inter-sensor differences. “X” in (Expression 1) indicates the value of the difference between sensors. As shown in (Expression 2), the standard deviation “s” is calculated using the average value calculated in (Expression 1).

つまり、血管状態推定部406は、センサ間差分各々について標準偏差を算出し、算出したセンサ間差分に対応する年代を閾値テーブル301から取得する。そして、取得した年代を推定結果として利用者に出力する。例えば、図7に示す例を用いて具体的に説明すると、血管状態推定部406は、算出した標準偏差が「1」である場合には、「20代」であると推定する。そして、血管状態推定部406は、推定結果となる「20代」を出力部203を介して利用者に出力する。   That is, the blood vessel state estimation unit 406 calculates a standard deviation for each inter-sensor difference, and acquires the age corresponding to the calculated inter-sensor difference from the threshold table 301. And the acquired age is output to a user as an estimation result. For example, specifically described with reference to the example shown in FIG. 7, the blood vessel state estimation unit 406 estimates that it is “20's” when the calculated standard deviation is “1”. The blood vessel state estimation unit 406 then outputs “20's” as an estimation result to the user via the output unit 203.

なお、血管状態推定装置200は、既知のパーソナルコンピュータ、ワークステーション、サーバ、携帯電話、PHS(Personal Handyphone System)端末、移動体通信端末又はPDA(Personal Digital Assistant)を利用して実現しても良い。   The blood vessel state estimation apparatus 200 may be realized using a known personal computer, workstation, server, mobile phone, PHS (Personal Handyphone System) terminal, mobile communication terminal, or PDA (Personal Digital Assistant). .

例えば、サーバに、図2に示した閾値テーブル301を搭載する。また、サーバに、脈波間隔算出部401と、脈波間隔算出部402と、n次微分処理部403と、特徴点識別部404と、差分算出部405と、血管状態推定部406との各機能を搭載することによって実現しても良い。この場合、例えば、サーバは、利用者が利用した光センサからデータを受信すると、血管の状態を推定し、推定結果を利用者が利用する情報端末に出力する。   For example, the threshold table 301 shown in FIG. 2 is installed in the server. Further, the server includes a pulse wave interval calculation unit 401, a pulse wave interval calculation unit 402, an nth-order differentiation processing unit 403, a feature point identification unit 404, a difference calculation unit 405, and a blood vessel state estimation unit 406. You may implement | achieve by mounting a function. In this case, for example, when the server receives data from the optical sensor used by the user, the server estimates the state of the blood vessel and outputs the estimation result to the information terminal used by the user.

また、例えば、図2に示した各部を複数の装置に分散して搭載することで、血管の状態を推定しても良い。例えば、利用者が利用する情報端末(例えば、携帯電話)に、図2に示した各部のうち、光センサ201と、光センサ202と、脈波間隔算出部401と、脈波間隔算出部402と、n次微分処理部403と、特徴点識別部404と、差分算出部405とを搭載する。また、サーバに、閾値テーブル301と、血管状態推定部406とを搭載する。   Further, for example, the state of the blood vessel may be estimated by distributing and mounting each unit illustrated in FIG. 2 in a plurality of apparatuses. For example, an information terminal (for example, a mobile phone) used by a user may include an optical sensor 201, an optical sensor 202, a pulse wave interval calculation unit 401, and a pulse wave interval calculation unit 402 among the units illustrated in FIG. And an nth-order differentiation processing unit 403, a feature point identification unit 404, and a difference calculation unit 405. In addition, a threshold table 301 and a blood vessel state estimation unit 406 are mounted on the server.

この場合、利用者が利用する情報端末は、差分算出部405によって算出された差分のデータをサーバに送信する。そして、サーバは、差分のデータを受信すると、受信した差分のデータに基づいて血管の状態を推定する。その後、サーバは、推定処理が終了すると、利用者が利用する情報端末に推定結果を送信する。また、例えば、サーバは、推定結果をデータベースに記憶しておき、利用者からサーバに対するアクセスがあると、データベースから推定結果を読み出して出力する。なお、利用者が利用する情報端末がサーバに送信する情報は、差分算出部405によって算出された差分である場合を例に説明したが、本発明はこれに限定されるものではない。例えば、利用者が利用する情報端末は、センサ間差分を算出した上でサーバに送信しても良く、脈波間隔算出部401〜差分算出部405によって実行される処理のうちどの段階のデータを送信しても良い。   In this case, the information terminal used by the user transmits the difference data calculated by the difference calculation unit 405 to the server. When the server receives the difference data, the server estimates the state of the blood vessel based on the received difference data. Thereafter, when the estimation process ends, the server transmits the estimation result to the information terminal used by the user. In addition, for example, the server stores the estimation result in a database, and reads and outputs the estimation result from the database when the user accesses the server. In addition, although the case where the information that the information terminal used by the user transmits to the server is a difference calculated by the difference calculation unit 405, the present invention is not limited to this. For example, the information terminal used by the user may calculate the difference between the sensors and transmit it to the server, and which stage data among the processes executed by the pulse wave interval calculation unit 401 to the difference calculation unit 405 may be used. You may send it.

[血管状態推定装置による処理]
次に、図14を用いて、実施例2に係る血管状態推定装置200による処理の流れの一例について説明する。図14は、実施例2に係る血管状態推定装置による処理の流れの一例について説明するフローチャートである。
[Processing by blood vessel state estimation device]
Next, an example of the flow of processing performed by the blood vessel state estimation apparatus 200 according to the second embodiment will be described with reference to FIG. FIG. 14 is a flowchart illustrating an example of a process flow performed by the blood vessel state estimation device according to the second embodiment.

図14に示すように、光センサは、他の光センサと同期をとって検出された光の変化量に基づいて人体の脈波を検出すると(ステップS101肯定)、脈波を電気信号に変換する(ステップS102)。つまり、光センサ201は、光センサ202と同期をとって検出された光の変化量に基づいて人体の脈波を検出し、脈波を電気信号に変換する。また、光センサ202は、光センサ201と同期をとって検出された光の変化量に基づいて人体の脈波を検出し、脈波を電気信号に変換する。   As shown in FIG. 14, when the optical sensor detects a pulse wave of a human body based on the amount of change in light detected in synchronization with another optical sensor (Yes in step S101), the optical sensor converts the pulse wave into an electrical signal. (Step S102). That is, the optical sensor 201 detects the pulse wave of the human body based on the change amount of the light detected in synchronization with the optical sensor 202, and converts the pulse wave into an electrical signal. The optical sensor 202 detects a pulse wave of the human body based on the change amount of light detected in synchronization with the optical sensor 201, and converts the pulse wave into an electrical signal.

そして、脈波間隔算出部401や脈波間隔算出部402は、それぞれ、光センサによって変換された電気信号が示す波形ごとに脈波間隔を算出する(ステップS103)。つまり、脈波間隔算出部401は、電気信号が示す波形のうち最大点を識別し、隣接する最大点間の間隔を算出することで脈波間隔を算出する。そして、脈波間隔算出部401や脈波間隔算出部402は、それぞれ、脈波間隔を用いて脈飛びがあった箇所を識別する(ステップS104)。   Then, the pulse wave interval calculation unit 401 and the pulse wave interval calculation unit 402 respectively calculate the pulse wave interval for each waveform indicated by the electrical signal converted by the optical sensor (step S103). That is, the pulse wave interval calculation unit 401 identifies the maximum point in the waveform indicated by the electrical signal, and calculates the pulse wave interval by calculating the interval between adjacent maximum points. Then, the pulse wave interval calculation unit 401 and the pulse wave interval calculation unit 402 each identify a place where a pulse jump has occurred using the pulse wave interval (step S104).

そして、n次微分処理部403は、光センサごとに、電気信号が示す波形に対してn次微分(nは自然数)を実行する(ステップS105)。例えば、n次微分処理部403は、光センサ201によって変換された電気信号に対して2次微分を実行し、また、光センサ202によって変換された電気信号に対して2次微分を実行する。   Then, the nth order differential processing unit 403 performs nth order differentiation (n is a natural number) on the waveform indicated by the electrical signal for each optical sensor (step S105). For example, the nth-order differentiation processing unit 403 performs second-order differentiation on the electric signal converted by the optical sensor 201 and performs second-order differentiation on the electric signal converted by the optical sensor 202.

そして、特徴点識別部404は、光センサごとに、n次微分の結果得られた波形であるn次微分波形に含まれる特徴点を識別する(ステップS106)。例えば、図10の(1)に示すように、特徴点識別部404は、耳についての2次微分波形に含まれる脈波ごとに特徴点を複数識別する。また、図10の(2)に示すように、特徴点識別部404は、指についての2次微分波形に含まれる脈波ごとに特徴点を複数識別する。   Then, the feature point identifying unit 404 identifies the feature points included in the nth-order differential waveform that is a waveform obtained as a result of the nth-order differentiation for each optical sensor (step S106). For example, as shown in (1) of FIG. 10, the feature point identifying unit 404 identifies a plurality of feature points for each pulse wave included in the secondary differential waveform for the ear. As shown in (2) of FIG. 10, the feature point identification unit 404 identifies a plurality of feature points for each pulse wave included in the secondary differential waveform for the finger.

そして、差分算出部405は、光センサごとに、特徴点識別部404によって識別された特徴点の出現時間の差分を算出する(ステップS107)。例えば、図11の(1)に示すように、差分算出部405は、耳についての2次微分波形に含まれる脈拍について、耳差分1〜4を算出する。また、図11の(2)に示すように、差分算出部405は、指についての2次微分波形に含まれる脈拍について、指差分1〜4を算出する。   Then, the difference calculation unit 405 calculates the difference in appearance time of the feature points identified by the feature point identification unit 404 for each optical sensor (step S107). For example, as illustrated in (1) of FIG. 11, the difference calculation unit 405 calculates ear differences 1 to 4 for the pulse included in the secondary differential waveform for the ear. Moreover, as shown in (2) of FIG. 11, the difference calculation part 405 calculates the finger differences 1-4 about the pulse contained in the secondary differential waveform about a finger.

そして、血管状態推定部406は、センサ間差分を複数算出する(ステップS108)。つまり、図10に示す例を用いて説明すると、血管状態推定部406は、耳差分1〜4と、指差分1〜4との間の差分をそれぞれ算出することで、センサ間差分1〜4を算出する。   Then, the blood vessel state estimation unit 406 calculates a plurality of inter-sensor differences (step S108). In other words, using the example illustrated in FIG. 10, the blood vessel state estimation unit 406 calculates the difference between the ear differences 1 to 4 and the finger differences 1 to 4, thereby calculating the inter-sensor differences 1 to 4. Is calculated.

そして、血管状態推定部406は、血管の状態を推定する(ステップS109)。つまり、血管状態推定部406は、センサ間差分各々が広く分散した場合には狭く分散した場合よりも血管年齢が高いと推定し、センサ間差分各々が狭く分散した場合には広く分散した場合よりも血管年齢が低いと推定する。   Then, the blood vessel state estimation unit 406 estimates the blood vessel state (step S109). In other words, the blood vessel state estimation unit 406 estimates that the age of the blood vessel is higher than when the inter-sensor differences are widely dispersed, and the case where the inter-sensor differences are widely dispersed, than when the inter-sensor differences are widely dispersed. Also assume that the vascular age is low.

そして、血管状態推定部406は、推定結果を出力する(ステップS110)。例えば、血管状態推定部406は、血管状態推定装置200を利用する利用者に対して、出力部203を介して推定結果を出力する。例えば、血管状態推定部406は、「20代」であると出力する。   Then, the blood vessel state estimation unit 406 outputs an estimation result (step S110). For example, the blood vessel state estimation unit 406 outputs an estimation result via the output unit 203 to a user who uses the blood vessel state estimation device 200. For example, the blood vessel state estimation unit 406 outputs that it is “20's”.

なお、図14に示した処理手順は、上述した順番に限定されるものではなく、処理内容を矛盾させない範囲で適宜変更しても良い。例えば、血管状態推定装置200は、ステップS109において血管の状態を推定すると、推定結果を記憶部300に格納して処理を終了しても良い。この場合、例えば、血管状態推定装置200は、利用者から推定結果の要求を別途受け付けると、記憶部300に格納された推定結果を読み出して出力する。また、例えば、血管状態推定装置200は、ステップS104とS105やS106とを並列して実行しても良い。この場合、血管状態推定装置200は、脈飛びの箇所を識別する処理を実行するとともに、2次微分処理を実行して特徴点を識別する。そして、血管状態推定装置200は、脈飛びの箇所を識別した処理結果と、識別した特徴点とを用いて、センタ間差分を算出する。   Note that the processing procedure shown in FIG. 14 is not limited to the above-described order, and may be appropriately changed within a range that does not contradict the processing contents. For example, when the blood vessel state estimation device 200 estimates the blood vessel state in step S109, the blood vessel state estimation device 200 may store the estimation result in the storage unit 300 and terminate the process. In this case, for example, when the blood vessel state estimation apparatus 200 separately receives a request for an estimation result from the user, the blood vessel state estimation apparatus 200 reads and outputs the estimation result stored in the storage unit 300. Further, for example, the blood vessel state estimation device 200 may execute steps S104 and S105 or S106 in parallel. In this case, the blood vessel state estimation device 200 executes a process of identifying a pulse skipping portion and executes a secondary differentiation process to identify a feature point. Then, the blood vessel state estimation apparatus 200 calculates the inter-center difference using the processing result that identifies the skipped portion and the identified feature point.

[実施例2の効果]
上述したように、実施例2によれば、血管状態推定装置200は、少なくとも2つの光センサを有する。また、血管状態推定装置200は、光センサごとに、他の光センサと同期をとって検出された光の変化量に基づいて人体の脈波を検出し、脈波を電気信号に変換する。そして、血管状態推定装置200は、光センサごとに、電気信号に対してn次微分(nは自然数)を実行し、n次微分の結果得られた波形であるn次微分波形に含まれる特徴点を識別する。そして、血管状態推定装置200は、光センサごとに、特徴点の出現時間の差分を算出し、光センサごとに算出された差分に基づいて人体の血管の状態を推定する。そして、血管状態推定装置200は、推定結果を出力する。この結果、実施例2によれば、血管の状態を精度良く簡易に推定可能である。
[Effect of Example 2]
As described above, according to the second embodiment, the blood vessel state estimation device 200 includes at least two optical sensors. Moreover, the blood vessel state estimation apparatus 200 detects the pulse wave of a human body based on the variation | change_quantity of the light detected synchronizing with another optical sensor for every optical sensor, and converts a pulse wave into an electrical signal. The blood vessel state estimation apparatus 200 performs an n-order differentiation (n is a natural number) on the electrical signal for each optical sensor, and is included in the n-order differentiation waveform that is a waveform obtained as a result of the n-order differentiation. Identify points. Then, the blood vessel state estimation device 200 calculates the difference in appearance time of feature points for each optical sensor, and estimates the state of the human blood vessel based on the difference calculated for each optical sensor. Then, the blood vessel state estimation device 200 outputs an estimation result. As a result, according to the second embodiment, it is possible to easily and accurately estimate the state of the blood vessel.

すなわち、実施例2によれば、特徴点各々の出現時間は血圧値などによる変動を受けにくいことを踏まえ、特徴点の出現時間の差分に基づいて血管の状態を推定するので、精度良く推定することが可能である。また、カフを用いる手法と比較して、実施例2によれば、血管の状態を簡易に推定可能である。   That is, according to the second embodiment, since the appearance time of each feature point is not easily affected by fluctuations due to the blood pressure value or the like, the state of the blood vessel is estimated based on the difference in the appearance time of the feature point. It is possible. Further, according to the second embodiment, the state of the blood vessel can be easily estimated as compared with the method using the cuff.

また、実施例2によれば、血管状態推定装置200は、光センサごとに差分を複数算出する。そして、血管状態推定装置200は、光センサごとに算出された差分間の差分であるセンサ間差分を、差分算出部が算出した複数の差分ごとに算出する。そして、血管状態推定装置200は、センサ間差分各々が広く分散した場合には狭く分散した場合よりも血管年齢が高いと推定し、センサ間差分各々が狭く分散した場合には広く分散した場合よりも血管年齢が低いと推定する。この結果、実施例2によれば、センタ間差分の分散について算出するだけで、血管の状態を精度良く簡易に推定可能である。   Further, according to the second embodiment, the blood vessel state estimation device 200 calculates a plurality of differences for each optical sensor. Then, the blood vessel state estimation device 200 calculates an inter-sensor difference, which is a difference between differences calculated for each optical sensor, for each of a plurality of differences calculated by the difference calculation unit. The blood vessel state estimation apparatus 200 estimates that the age of the blood vessel is higher than when the inter-sensor differences are widely dispersed, and more than when the inter-sensor differences are widely dispersed. Also assume that the vascular age is low. As a result, according to the second embodiment, it is possible to easily and accurately estimate the state of the blood vessel simply by calculating the variance of the difference between the centers.

また、実施例2によれば、少なくとも2つある光センサは、それぞれ異なる太さの血管がある人体の部分に光を照射するので、血管の状態を精度良く推定可能である。つまり、血管の太さが違えば、加齢による血管の伸展性や弾力性に対する影響の大きさが異なることを踏まえ、それぞれ異なる太さの血管において光の変化量を検知して血管の状態を推定するので、血管の状態を精度良く推定可能である。   Further, according to the second embodiment, at least two photosensors irradiate light to a part of a human body having blood vessels of different thicknesses, so that the state of the blood vessels can be estimated with high accuracy. In other words, if the blood vessel thickness is different, the magnitude of the effect of aging on blood vessel extensibility and elasticity is different. Since the estimation is performed, the state of the blood vessel can be estimated with high accuracy.

また、実施例2によれば、血管状態推定装置200は、光センサごとに変換された電気信号に対して2次微分を実行し、2次微分波形に含まれる脈波ごとに少なくとも4つの特徴点を識別する。そして、血管状態推定装置200は、他方の光センサと比較して太い血管に光を照射する光センサについて、脈波ごとに、3番目の特徴点と4番目の特徴点との間の第1の差分を算出する。また、血管状態推定装置200は、他方の光センサと比較して細い血管に光を照射する光センサについて、脈波ごとに、1番目の特徴点と2番目の特徴点との間の第2の差分を算出する。そして、血管状態推定装置200は、第1の差分と第2の差分との間の差分を脈波ごとに算出することでセンサ間差分各々を算出する。この結果、血管の状態を精度良く推定可能である。   In addition, according to the second embodiment, the blood vessel state estimation device 200 performs secondary differentiation on the electrical signal converted for each optical sensor, and at least four features for each pulse wave included in the secondary differential waveform. Identify points. Then, the blood vessel state estimation device 200 uses the first feature point between the third feature point and the fourth feature point for each pulse wave with respect to the light sensor that irradiates light to a blood vessel that is thicker than the other light sensor. The difference is calculated. In addition, the blood vessel state estimation device 200 uses a second between the first feature point and the second feature point for each pulse wave with respect to an optical sensor that irradiates light to a thin blood vessel as compared with the other optical sensor. The difference is calculated. And the vascular state estimation apparatus 200 calculates each difference between sensors by calculating the difference between a 1st difference and a 2nd difference for every pulse wave. As a result, the state of the blood vessel can be estimated with high accuracy.

すなわち、脈波には、心臓の鼓動により発生する駆動波と、駆動波が血管の分岐や壁にぶつかって反射することにより発生する反射波とがある。このため、電気信号が示す波形や2次微分波形に含まれる脈波は、駆動波と反射波との両方の成分が含まれる。また、血管の太さが同じであれば、血管の抹消になればなるほど、反射波の強さが強くなり、脈波のピークにおいて反射波が占める成分が大きくなると考えられる。また、血管の太さが細ければ細いほど、反射波の強さが強くなり、脈波のピークにおいて反射波が占める成分が大きくなると考えられる。また、血管は、老化すると、一度伸びた(拡張した)後の戻りが遅くなる結果、反射波に対する反応が悪くなると考えられる。   That is, the pulse wave includes a drive wave generated by the heartbeat and a reflected wave generated when the drive wave collides with a blood vessel branch or wall and reflects. For this reason, the pulse wave included in the waveform indicated by the electric signal or the second-order differential waveform includes both components of the drive wave and the reflected wave. Further, if the blood vessels have the same thickness, the more the blood vessels are erased, the stronger the intensity of the reflected wave, and the larger the component occupied by the reflected wave at the peak of the pulse wave. Further, it is considered that the thinner the blood vessel, the stronger the reflected wave, and the larger the component occupied by the reflected wave at the peak of the pulse wave. In addition, when the blood vessel ages, it is considered that the response to the reflected wave becomes worse as a result of the slow return after being stretched (expanded).

このことを踏まえ、血管状態推定装置200は、血管の太さに応じて、駆動波の後にやってくると考えられる反射波の影響が相対的に大きい考えられる箇所の特徴点を用いて差分を算出する。耳と指とで光センサを使用する場合を例に説明すると、耳の場合には、指の血管と比較して血管が細い結果、指の場合と比較して反射波の影響が強くなると考えられ、脈波形のピーク付近に反射波の影響が現れると考えられる。このことを踏まえ、血管状態推定装置200は、耳の場合には、2次微分波形の特徴点のうち前の方にある特徴点a−bを用いて差分を算出する。また、指の場合には、耳の血管と比較して血管が太い結果、耳の場合と比較して反射波の影響が弱いと考えられ、脈波形のピークよりも後ろに反射波の影響が現れると考えられる。このことを踏まえ、血管状態推定装置200は、後ろの方にある特徴点c−dを用いて差分を算出する。このように、血管状態推定装置200は、反射波の影響が強くでると考えられる箇所を用いて血管の状態を推定するので、加齢による影響を精度良く検知でき、血管の状態を精度良く推定可能である。   Based on this, the blood vessel state estimation device 200 calculates the difference using the feature points of the places where the influence of the reflected wave that is considered to come after the driving wave is relatively large according to the thickness of the blood vessel. . In the case of using an optical sensor between the ear and the finger, for example, in the case of the ear, the blood vessel is thinner than the finger blood vessel, and the influence of the reflected wave is considered to be stronger than that of the finger. Therefore, the influence of the reflected wave appears near the peak of the pulse waveform. In consideration of this, in the case of the ear, the blood vessel state estimation device 200 calculates a difference using the feature point ab located earlier in the feature points of the secondary differential waveform. In the case of the finger, the blood vessel is thicker than the ear blood vessel. As a result, the influence of the reflected wave is considered to be weaker than that of the ear, and the influence of the reflected wave is behind the peak of the pulse waveform. It is thought to appear. Based on this, the blood vessel state estimation device 200 calculates a difference using the feature point cd at the back. In this way, since the blood vessel state estimation apparatus 200 estimates the blood vessel state using a place where the influence of the reflected wave is considered to be strong, the blood vessel state can be accurately detected and the blood vessel state can be accurately estimated. Is possible.

さて、これまで本発明の実施例1及び2について説明したが、本発明は上記した実施例1及び2以外にも、その他の実施例にて実施されても良い。そこで、以下では、その他の実施例について説明する。   Although the first and second embodiments of the present invention have been described so far, the present invention may be implemented in other embodiments in addition to the first and second embodiments described above. Therefore, other embodiments will be described below.

例えば、上述の実施例では、血管状態推定装置が、推定結果として年代を出力する場合を例に説明したが、本発明はこれに限定されるものではない。例えば、血管状態推定装置は、血管の状態が良好か良好ではないかを出力しても良く、血管年齢が利用者の年齢よりも若いか老いているかを出力しても良く、利用者が任意に設定して良い。   For example, in the above-described embodiment, the case where the blood vessel state estimation device outputs the age as the estimation result has been described as an example, but the present invention is not limited to this. For example, the blood vessel state estimation device may output whether the blood vessel state is good or not, may output whether the blood vessel age is younger or older than the user's age, and the user can arbitrarily May be set to

[システム構成]
また、本実施例において説明した各処理のうち、自動的に行われるものとして説明した処理の全部又は一部を手動的に行うこともでき、あるいは、手動的に行われるものとして説明した処理の全部又は一部を公知の方法で自動的に行うこともできる。例えば、上述した実施例では、血管状態推定装置200が、血管状態を推定し、推定結果を出力する場合について説明したが、本発明はこれに限定されるものではない。例えば、血管状態推定装置200は、血管状態を推定することなく、センサ間差分についての標準偏差の値と閾値テーブルとを併せて出力するにとどめても良い。
[System configuration]
Also, among the processes described in this embodiment, all or part of the processes described as being performed automatically can be performed manually, or the processes described as being performed manually can be performed. All or a part can be automatically performed by a known method. For example, in the above-described embodiment, the case where the blood vessel state estimation device 200 estimates the blood vessel state and outputs the estimation result has been described, but the present invention is not limited to this. For example, the blood vessel state estimation apparatus 200 may output only the standard deviation value and the threshold value table for the inter-sensor difference without estimating the blood vessel state.

この他、上記文書中や図面中で示した処理手順、制御手順、具体的名称、各種のデータやパラメータを含む情報(図1〜図14)については、特記する場合を除いて任意に変更することができる。例えば、図3に示した閾値テーブル301のデータは、利用者が任意の値に変更しても良い。   In addition, the processing procedures, control procedures, specific names, and information including various data and parameters (FIGS. 1 to 14) shown in the documents and drawings are arbitrarily changed unless otherwise specified. be able to. For example, the data in the threshold table 301 shown in FIG. 3 may be changed by the user to an arbitrary value.

また、図示した各装置の各構成要素は機能概念的なものであり、必ずしも物理的に図示の如く構成されていることを要しない。すなわち、各装置の分散・統合の具体的形態は図示のものに限られず、その全部又は一部を、各種の負荷や使用状況などに応じて、任意の単位で機能的又は物理的に分散・統合して構成することができる。例えば、光センサ201と光センサ202と出力部203とを血管状態推定装置200の外部装置としてネットワーク経由で接続するようにしても良い。また、血管状態推定部406と閾値テーブル301とを別の装置が有し、ネットワーク経由で接続されて協働することで、上述した血管状態推定装置200の機能を実現するようにしても良い。   Further, each component of each illustrated apparatus is functionally conceptual, and does not necessarily need to be physically configured as illustrated. That is, the specific form of distribution / integration of each device is not limited to the one shown in the figure, and all or a part of the distribution / integration may be functionally or physically distributed in arbitrary units according to various loads or usage conditions. Can be integrated and configured. For example, the optical sensor 201, the optical sensor 202, and the output unit 203 may be connected as an external device of the blood vessel state estimation device 200 via a network. Further, the blood vessel state estimation unit 406 and the threshold value table 301 may be provided in different devices, and may be connected via a network to cooperate to realize the above-described function of the blood vessel state estimation device 200.

[コンピュータ]
また、上記の実施例で説明した各種の処理は、予め用意されたプログラムをパーソナルコンピュータやワークステーションなどのコンピュータで実行することによって実現することができる。そこで、以下では、図15を用いて、上記の実施例と同様の機能を有する血管状態推定プログラムを実行するコンピュータの一例について説明する。なお、図15は、実施例2に係る血管状態推定プログラムを実行するコンピュータの一例について説明する図である。
[Computer]
The various processes described in the above embodiments can be realized by executing a prepared program on a computer such as a personal computer or a workstation. In the following, an example of a computer that executes a blood vessel state estimation program having the same function as that of the above-described embodiment will be described with reference to FIG. FIG. 15 is a schematic diagram illustrating an example of a computer that executes a blood vessel state estimation program according to the second embodiment.

図15に示すように、実施例2におけるコンピュータ3000は、操作部3001、マイク3002、スピーカ3003、ディスプレイ3004、光センサ3005、通信部3006、CPU3010、ROM3011、HDD(Hard Disk Drive)3012、RAM(Random Access Memory)3013をバス3009などで接続して構成されている。   As illustrated in FIG. 15, the computer 3000 according to the second embodiment includes an operation unit 3001, a microphone 3002, a speaker 3003, a display 3004, an optical sensor 3005, a communication unit 3006, a CPU 3010, a ROM 3011, a HDD (Hard Disk Drive) 3012, a RAM ( Random Access Memory) 3013 is connected by a bus 3009 or the like.

ROM3011には、上記の実施例2で示した脈波間隔算出部401〜血管状態推定部406の各部と同様の機能を発揮する制御プログラムが格納される。すなわち、図15に示す例では、ROM3011には、脈波間隔算出プログラム3011aと、n次微分処理プログラム3011bと、特徴点識別プログラム3011cとが予め記憶される。また、ROM3011には、差分算出プログラム3011dと、血管状態推定プログラム3011eとが予め記憶される。なお、これらのプログラム3011a〜3011eについては、図2に示した血管状態推定装置200の各構成要素と同様、適宜統合又は分離しても良い。   The ROM 3011 stores a control program that exhibits the same functions as those of the pulse wave interval calculation unit 401 to the blood vessel state estimation unit 406 shown in the second embodiment. In other words, in the example illustrated in FIG. 15, the pulse wave interval calculation program 3011a, the nth-order differentiation processing program 3011b, and the feature point identification program 3011c are stored in advance in the ROM 3011. Further, the ROM 3011 stores a difference calculation program 3011d and a blood vessel state estimation program 3011e in advance. In addition, about these programs 3011a-3011e, you may integrate or isolate | separate suitably like each component of the blood-vessel state estimation apparatus 200 shown in FIG.

ここで、脈波間隔算出プログラム3011aは、脈波間隔算出部402と同様の機能を発揮する制御プログラムである。n次微分処理プログラム3011bは、n次微分処理部403と同様の機能を発揮する制御プログラムである。特徴点識別プログラム3011cは、特徴点識別部404と同様の機能を発揮する制御プログラムである。   Here, the pulse wave interval calculation program 3011 a is a control program that exhibits the same function as the pulse wave interval calculation unit 402. The n-th order differential processing program 3011b is a control program that exhibits the same function as the n-th order differential processing unit 403. The feature point identification program 3011c is a control program that exhibits the same function as the feature point identification unit 404.

そして、CPU3010は、プログラム3011a〜3011bをROM3011から読み出して実行することで、脈波間隔算出プロセス3010aと、n次微分処理プロセス3010bとして機能する。また、CPU3010は、プログラム3011c〜3011eをROM3011から読み出して実行することで、n次微分処理プロセス3011cとして機能する。また、CPU3010は、プログラム3011d〜3011eをROM3011から読み出して実行することで、差分算出プロセス3010dと、血管状態推定プロセス3010eとして機能する。   The CPU 3010 reads out the programs 3011a to 3011b from the ROM 3011 and executes them, thereby functioning as a pulse wave interval calculation process 3010a and an nth-order differentiation process 3010b. In addition, the CPU 3010 functions as an nth-order differentiation process 3011c by reading and executing the programs 3011c to 3011e from the ROM 3011. The CPU 3010 functions as a difference calculation process 3010d and a blood vessel state estimation process 3010e by reading the programs 3011d to 3011e from the ROM 3011 and executing them.

なお、各プロセス3010a〜3010eは、図2に示した、脈波間隔算出部401と、脈波間隔算出部402と、n次微分処理部403と、特徴点識別部404と、差分算出部405と、血管状態推定部406とにそれぞれ対応する。   Each of the processes 3010a to 3010e includes a pulse wave interval calculation unit 401, a pulse wave interval calculation unit 402, an nth-order differential processing unit 403, a feature point identification unit 404, and a difference calculation unit 405 illustrated in FIG. And the blood vessel state estimation unit 406, respectively.

そして、HDD3012には、閾値テーブル3012aが設けられている。なお、閾値テーブル3012aは、図2に示した、閾値テーブル301に対応する。   The HDD 3012 is provided with a threshold table 3012a. Note that the threshold value table 3012a corresponds to the threshold value table 301 illustrated in FIG.

そして、CPU3010は、閾値テーブル3012aを読み出してRAM3013に格納する。そして、CPU3010は、RAM3013に格納された閾値データ3013aと、光の変化量データ3013bと、電気信号の波計データ3013cと、2次微分波計データ3013dと、差分データ3013eとを用いて、血管状態推定プログラムを実行する。   Then, the CPU 3010 reads the threshold table 3012a and stores it in the RAM 3013. Then, the CPU 3010 uses the threshold data 3013a, the light change amount data 3013b, the electric wave meter data 3013c, the second derivative wave meter data 3013d, and the difference data 3013e stored in the RAM 3013, using the blood vessel. Run the state estimation program.

[その他]
なお、本実施例で説明した血管状態推定プログラムは、インターネットなどのネットワークを介して配布することができる。また、血管状態推定プログラムは、ハードディスク、フレキシブルディスク(FD)、CD−ROM、MO、DVDなどのコンピュータで読み取り可能な記録媒体に記録され、コンピュータによって記録媒体から読み出されることによって実行することもできる。
[Others]
The blood vessel state estimation program described in this embodiment can be distributed via a network such as the Internet. The blood vessel state estimation program can also be executed by being recorded on a computer-readable recording medium such as a hard disk, a flexible disk (FD), a CD-ROM, an MO, or a DVD, and being read from the recording medium by the computer. .

100 血管状態推定装置
101 光センサ
102 光センサ
103 脈波検出部
104 特徴点識別部
105 差分算出部
106 血管状態推定部
200 血管状態推定装置
201 光センサ
202 光センサ
203 出力部
300 記憶部
301 閾値テーブル
400 制御部
401 脈波間隔算出部
402 脈波間隔算出部
403 n次微分処理部
404 特徴点識別部
405 差分算出部
406 血管状態推定部
501 開口制御ラバー
502 可視光カットフィルタ
503 搭載基板
504 カバーケース
505 発光素子
506 受光素子
DESCRIPTION OF SYMBOLS 100 Blood vessel state estimation apparatus 101 Optical sensor 102 Optical sensor 103 Pulse wave detection part 104 Feature point identification part 105 Difference calculation part 106 Blood vessel state estimation part 200 Blood vessel state estimation apparatus 201 Optical sensor 202 Optical sensor 203 Output part 300 Storage part 301 Threshold table 400 Control Unit 401 Pulse Wave Interval Calculation Unit 402 Pulse Wave Interval Calculation Unit 403 nth Order Differentiation Processing Unit 404 Feature Point Identification Unit 405 Difference Calculation Unit 406 Blood Vessel State Estimation Unit 501 Aperture Control Rubber 502 Visible Light Cut Filter 503 Mounted Substrate 504 Cover Case 505 Light emitting element 506 Light receiving element

Claims (6)

人体へ光を照射する発光素子と、前記発光素子から放射された後人体を透過した光の変化量、又は、前記発光素子から放射された後人体により反射された光の変化量を検出する受光素子とを有する少なくとも2つの光センサごとに、他の光センサと同期をとって検出された光の変化量に基づいて前記人体の脈波を検出し、該脈波を電気信号に変換する脈波検出部と、
前記光センサごとに、前記脈波検出部によって変換された電気信号に対してn次微分(nは自然数)を実行し、該n次微分の結果得られた波形であるn次微分波形に含まれる特徴点を識別する特徴点識別部と、
前記光センサごとに、前記特徴点識別部によって識別された特徴点の出現時間の差分を算出する差分算出部と、
前記差分算出部によって前記光センサごとに算出された差分に基づいて前記人体の血管の状態を推定し、推定結果を出力する血管状態推定部と
を備えたことを特徴とする血管状態推定装置。
A light emitting element that irradiates light to the human body, and a light receiving element that detects a change amount of light transmitted from the light emitting element and then transmitted through the human body, or a change amount of light reflected from the human body after being emitted from the light emitting element. A pulse for detecting a pulse wave of the human body based on a change amount of light detected in synchronization with another optical sensor and converting the pulse wave into an electric signal for each of at least two optical sensors having an element; A wave detector;
For each of the optical sensors, n-order differentiation (n is a natural number) is performed on the electrical signal converted by the pulse wave detection unit, and the waveform is obtained as a result of the n-order differentiation. A feature point identifying unit for identifying the feature point to be
A difference calculation unit that calculates a difference in appearance time of feature points identified by the feature point identification unit for each of the optical sensors;
A blood vessel state estimation device comprising: a blood vessel state estimation unit that estimates a blood vessel state of the human body based on a difference calculated for each of the optical sensors by the difference calculation unit and outputs an estimation result.
前記差分算出部は、前記光センサごとに差分を複数算出し、
前記血管状態推定部は、前記光センサごとに算出された差分間の差分であるセンサ間差分を、前記差分算出部が算出した複数の差分ごとに算出し、該センサ間差分各々が広く分散した場合には狭く分散した場合よりも血管年齢が高いと推定し、該センサ間差分各々が狭く分散した場合には広く分散した場合よりも血管年齢が低いと推定することを特徴とする請求項1に記載の血管状態推定装置。
The difference calculation unit calculates a plurality of differences for each of the optical sensors,
The blood vessel state estimation unit calculates an inter-sensor difference, which is a difference between differences calculated for each optical sensor, for each of a plurality of differences calculated by the difference calculation unit, and each of the inter-sensor differences is widely dispersed. The blood vessel age is estimated to be higher than the case where the blood vessel age is higher than that in the case of narrow dispersion, and the blood vessel age is estimated to be lower than that in the case of wide dispersion when each of the differences between the sensors is narrowly distributed. The blood vessel state estimation device according to 1.
少なくとも2つある前記光センサは、それぞれ異なる太さの血管がある人体の部分に光を照射することを特徴とする請求項1に記載の血管状態推定装置。   2. The blood vessel state estimation apparatus according to claim 1, wherein at least two of the optical sensors irradiate light to a part of a human body having blood vessels of different thicknesses. 前記特徴点識別部は、前記光センサごとに、前記脈波検出部によって変換された電気信号に対して2次微分を実行し、2次微分波形に含まれる脈波ごとに少なくとも4つの特徴点を識別し、
前記差分算出部は、他方の光センサと比較して太い血管に光を照射する光センサについて、脈波ごとに、少なくとも4つの特徴点のうち出現時間が前から3番目になる3番特徴点と4番目になる4番特徴点との間の第1の差分を算出し、他方の光センサと比較して細い血管に光を照射する光センサについて、脈波ごとに、少なくとも4つの特徴点のうち出現時間が前から1番目になる1番特徴点と2番目になる2番特徴点との間の第2の差分を算出し、
前記血管状態推定部は、第1の差分と第2の差分との間の差分を脈波ごとに算出することでセンサ間差分各々を算出することを特徴とする請求項3に記載の血管状態推定装置。
The feature point identification unit performs second order differentiation on the electrical signal converted by the pulse wave detection unit for each of the optical sensors, and at least four feature points for each pulse wave included in the second order differential waveform Identify and
For the optical sensor that irradiates light to a blood vessel that is thicker than the other optical sensor, the difference calculating unit is the third feature point that has the third appearance time among the at least four feature points for each pulse wave. And at least four feature points for each pulse wave with respect to an optical sensor that irradiates light to a thin blood vessel compared to the other optical sensor. A second difference between the first feature point of which appearance time is the first and the second feature point of the second from the previous time,
The blood vessel state estimation unit according to claim 3, wherein the blood vessel state estimation unit calculates each difference between sensors by calculating a difference between the first difference and the second difference for each pulse wave. Estimating device.
情報処理装置により実行される方法であって、
人体へ光を照射する発光素子と、前記発光素子から放射された後人体を透過した光の変化量、又は、前記発光素子から放射された後人体により反射された光の変化量を検出する受光素子とを有する少なくとも2つの光センサごとに、他の光センサと同期をとって検出された光の変化量に基づいて前記人体の脈波を検出し、該脈波を電気信号に変換する脈波検出ステップと、
前記光センサごとに、前記脈波検出ステップによって変換された電気信号に対してn次微分(nは自然数)を実行し、該n次微分の結果得られた波形であるn次微分波形に含まれる特徴点を識別する特徴点識別ステップと、
前記光センサごとに、前記特徴点識別ステップによって識別された特徴点の出現時間の差分を算出する差分算出ステップと、
前記差分算出ステップによって前記光センサごとに算出された差分に基づいて前記人体の血管の状態を推定し、推定結果を出力する血管状態推定ステップと
を含んだことを特徴とする血管状態推定方法。
A method executed by an information processing apparatus,
A light emitting element that irradiates light to the human body, and a light receiving element that detects a change amount of light transmitted from the light emitting element and then transmitted through the human body, or a change amount of light reflected from the human body after being emitted from the light emitting element. A pulse for detecting a pulse wave of the human body based on a change amount of light detected in synchronization with another optical sensor and converting the pulse wave into an electric signal for each of at least two optical sensors having an element; A wave detection step;
For each of the optical sensors, n-order differentiation (n is a natural number) is performed on the electrical signal converted by the pulse wave detection step, and the waveform is obtained as a result of the n-order differentiation. A feature point identifying step for identifying the feature points to be
A difference calculating step for calculating a difference between appearance times of feature points identified by the feature point identifying step for each of the optical sensors;
A blood vessel state estimation method comprising: estimating a blood vessel state of the human body based on the difference calculated for each of the optical sensors in the difference calculation step, and outputting an estimation result.
情報処理装置に、
人体へ光を照射する発光素子と、前記発光素子から放射された後人体を透過した光の変化量、又は、前記発光素子から放射された後人体により反射された光の変化量を検出する受光素子とを有する少なくとも2つの光センサごとに、他の光センサと同期をとって検出された光の変化量に基づいて前記人体の脈波を検出し、該脈波を電気信号に変換する脈波検出手順と、
前記光センサごとに、前記脈波検出手順によって変換された電気信号に対してn次微分(nは自然数)を実行し、該n次微分の結果得られた波形であるn次微分波形に含まれる特徴点を識別する特徴点識別手順と、
前記光センサごとに、前記特徴点識別手順によって識別された特徴点の出現時間の差分を算出する差分算出手順と、
前記差分算出手順によって前記光センサごとに算出された差分に基づいて前記人体の血管の状態を推定し、推定結果を出力する血管状態推定手順と
を実行させることを特徴とする血管状態推定プログラム。
In the information processing device,
A light emitting element that irradiates light to the human body, and a light receiving element that detects a change amount of light transmitted from the light emitting element and then transmitted through the human body, or a change amount of light reflected from the human body after being emitted from the light emitting element. A pulse for detecting a pulse wave of the human body based on a change amount of light detected in synchronization with another optical sensor and converting the pulse wave into an electric signal for each of at least two optical sensors having an element; Wave detection procedure;
For each of the optical sensors, n-order differentiation (n is a natural number) is performed on the electrical signal converted by the pulse wave detection procedure, and the waveform is obtained as a result of the n-order differentiation. A feature point identification procedure for identifying a feature point to be
A difference calculation procedure for calculating a difference between appearance times of feature points identified by the feature point identification procedure for each of the optical sensors;
A blood vessel state estimation program for executing a blood vessel state estimation procedure for estimating a blood vessel state of the human body based on a difference calculated for each of the optical sensors by the difference calculation procedure and outputting an estimation result.
JP2009228838A 2009-09-30 2009-09-30 Blood vessel state estimation device, blood vessel state estimation method, and blood vessel state estimation program Expired - Fee Related JP5381586B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2009228838A JP5381586B2 (en) 2009-09-30 2009-09-30 Blood vessel state estimation device, blood vessel state estimation method, and blood vessel state estimation program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2009228838A JP5381586B2 (en) 2009-09-30 2009-09-30 Blood vessel state estimation device, blood vessel state estimation method, and blood vessel state estimation program

Publications (2)

Publication Number Publication Date
JP2011072674A true JP2011072674A (en) 2011-04-14
JP5381586B2 JP5381586B2 (en) 2014-01-08

Family

ID=44017298

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2009228838A Expired - Fee Related JP5381586B2 (en) 2009-09-30 2009-09-30 Blood vessel state estimation device, blood vessel state estimation method, and blood vessel state estimation program

Country Status (1)

Country Link
JP (1) JP5381586B2 (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014157714A1 (en) * 2013-03-26 2014-10-02 株式会社志成データム Vascular viscoelasticity evaluation device, vascular viscoelasticity evaluation method, and program
JP2016147052A (en) * 2015-02-13 2016-08-18 台医光電科技股▲ふん▼有限公司 Multi-site measurement accessory, multi-site measurement device and multi-site measurement system
JP2020195708A (en) * 2019-06-05 2020-12-10 富士通コネクテッドテクノロジーズ株式会社 Pulse wave analyzer, pulse wave analysis method, and pulse wave analysis program
US11660025B2 (en) 2017-07-11 2023-05-30 Nec Corporation Identification device, identification method, and recording medium with recorded identification program

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000051166A (en) * 1998-08-05 2000-02-22 Kenji Takazawa Organism information processor
JP2006102250A (en) * 2004-10-06 2006-04-20 Terumo Corp Circulation index measuring device, circulation index measuring method, control program, and computer-readable storage medium
JP2007044429A (en) * 2005-08-12 2007-02-22 Tottori Univ Pulse wave diagnostic system using self-organizing map (som), generating program of self-organizing map (som) and generating process
JP2009095511A (en) * 2007-10-18 2009-05-07 Hitachi Ltd Bio-instrumentation system

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000051166A (en) * 1998-08-05 2000-02-22 Kenji Takazawa Organism information processor
JP2006102250A (en) * 2004-10-06 2006-04-20 Terumo Corp Circulation index measuring device, circulation index measuring method, control program, and computer-readable storage medium
JP2007044429A (en) * 2005-08-12 2007-02-22 Tottori Univ Pulse wave diagnostic system using self-organizing map (som), generating program of self-organizing map (som) and generating process
JP2009095511A (en) * 2007-10-18 2009-05-07 Hitachi Ltd Bio-instrumentation system

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014157714A1 (en) * 2013-03-26 2014-10-02 株式会社志成データム Vascular viscoelasticity evaluation device, vascular viscoelasticity evaluation method, and program
JP2014188035A (en) * 2013-03-26 2014-10-06 Shisei Deetamu:Kk Vascular viscoelasticity evaluation apparatus, method, and program
US10506932B2 (en) 2013-03-26 2019-12-17 Shisei Datum Co., Ltd. Vascular viscoelasticity evaluation device, vascular viscoelasticity evaluation method, and program
JP2016147052A (en) * 2015-02-13 2016-08-18 台医光電科技股▲ふん▼有限公司 Multi-site measurement accessory, multi-site measurement device and multi-site measurement system
US11660025B2 (en) 2017-07-11 2023-05-30 Nec Corporation Identification device, identification method, and recording medium with recorded identification program
JP2020195708A (en) * 2019-06-05 2020-12-10 富士通コネクテッドテクノロジーズ株式会社 Pulse wave analyzer, pulse wave analysis method, and pulse wave analysis program

Also Published As

Publication number Publication date
JP5381586B2 (en) 2014-01-08

Similar Documents

Publication Publication Date Title
JP5742441B2 (en) Biological information processing device
CN103961080B (en) Bioinformation detecting device and Biont information detection method
US11633115B2 (en) Apparatus and method for estimating bio-information
RU2675399C2 (en) Heart rate monitoring system
US20180078215A1 (en) Feature extraction apparatus and method for biometric information detection, biometric information detection apparatus, and wearable device
JP5549092B2 (en) Blood pressure value measuring apparatus and blood pressure value measuring method
JP5381586B2 (en) Blood vessel state estimation device, blood vessel state estimation method, and blood vessel state estimation program
US11179045B2 (en) Blood pressure measurement device and blood pressure measurement method
JP2012095940A (en) Pulse rate measuring method and blood oxygen saturation measuring method
JP2014236775A (en) Organism information measuring device and organism information measurement method
US11844631B2 (en) Apparatus and method for estimating bio-information
JP6213943B1 (en) Blood pressure estimation device, blood pressure estimation method, and computer program
JP5540662B2 (en) Respiration estimation device, respiration estimation method, and respiration estimation program
JP6115329B2 (en) Biological information processing apparatus and biological information processing method
KR20190041837A (en) Apparatus and method for bio information processing
CN109106383A (en) mood sensing system and method
JP6152728B2 (en) Biological information detection device
JP6436186B2 (en) Biological information detection device
EP3854301A1 (en) Apparatus and method for estimating bio-information
JP2014008200A (en) Pulse wave signal processing device, blood pressure measuring apparatus, and pulse wave signal processing program
JP5655721B2 (en) Biological information processing apparatus and biological information processing method
JP2011212383A (en) Biological information processor
JP6036916B2 (en) Biological information processing device
US11813039B2 (en) Apparatus and method for estimating bio-information
JC Spot and continuous monitoring of heart rate by combining time and frequency domain analysis of photoplethysmographic signals at rest conditions.

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20120605

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20130829

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20130903

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20130916

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

LAPS Cancellation because of no payment of annual fees