JP7083185B2 - Biological information calculation system - Google Patents
Biological information calculation system Download PDFInfo
- Publication number
- JP7083185B2 JP7083185B2 JP2020206629A JP2020206629A JP7083185B2 JP 7083185 B2 JP7083185 B2 JP 7083185B2 JP 2020206629 A JP2020206629 A JP 2020206629A JP 2020206629 A JP2020206629 A JP 2020206629A JP 7083185 B2 JP7083185 B2 JP 7083185B2
- Authority
- JP
- Japan
- Prior art keywords
- unit
- data
- pulse wave
- extraction
- wave signal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Landscapes
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Description
本発明は、生体情報演算システムに関する。 The present invention relates to a biometric information calculation system.
従来、ユーザから血液を採取する侵襲測定方法では、ユーザから生体情報を直接測定する場合がある。例えば生体情報として血糖値を測定する際には、穿刺によって血液を採取し、酵素電極法などの血糖値センサで測定する。このような生体情報の測定のための侵襲測定方法は、血液採取においてユーザの心身への負担が生じる。 Conventionally, in the invasive measurement method of collecting blood from a user, biological information may be directly measured from the user. For example, when measuring a blood glucose level as biological information, blood is collected by puncture and measured by a blood glucose level sensor such as an enzyme electrode method. Such an invasive measurement method for measuring biological information causes a burden on the user's mind and body in blood collection.
このため、ユーザの心身への負担を軽減するために非侵襲生体情報測定法の技術が必要とされている。この非侵襲生体情報測定法によれば、例えばひずみセンサ、ジャイロセンサ、光電容積脈波計によって測定された脈波の波形信号から、予め定めた脈波と生体情報の相関関係に基づき、ユーザの生体情報を求める技術が注目されている。 Therefore, a technique for non-invasive biometric information measurement is required in order to reduce the burden on the user's mind and body. According to this non-invasive biometric information measurement method, for example, from a waveform signal of a pulse wave measured by a strain sensor, a gyro sensor, or a photoelectric volume pulse wave meter, a user can use a predetermined correlation between the pulse wave and the biometric information. Technology for obtaining biometric information is attracting attention.
また、脈波には血糖値や血圧などの様々な生体情報と相関関係があることを利用し、脈波の波形信号から非侵襲的方法で生体情報を取得する生体情報推定装置及びその方法が提案されている(例えば特許文献1参照)。 In addition, a biometric information estimation device and its method that acquire biometric information from the waveform signal of the pulse wave by a non-invasive method by utilizing the fact that the pulse wave has a correlation with various biometric information such as blood pressure level and blood pressure. It has been proposed (see, for example, Patent Document 1).
特許文献1では、生体信号から複数の特徴値を抽出して、特徴値に基づいてスケール因子を決定し、スケール因子と特徴値に基づいて正確かつ安定した生体情報を推定する生体情報推定装置及びその方法が開示されている。
In
ここで、脈波の波形信号は、抽出条件や処理の方法に応じて、取得できる生体情報の数値にばらつきが生じる場合がある。例えば抽出条件として微分を行った脈波を処理した場合と、抽出条件として微分を行わなかった脈波を処理した場合とを比較したとき、取得できる生体情報の数値にばらつきが生じることが懸念される。すなわち、脈波の波形信号を利用して十分な精度がある生体情報を取得するためには、複数種類の抽出条件等を経て同時に出力された複数の生体情報から、脈波の波形信号を多角的に評価する必要がある。 Here, in the waveform signal of the pulse wave, the numerical value of the biometric information that can be acquired may vary depending on the extraction conditions and the processing method. For example, when comparing the case where the pulse wave differentiated as the extraction condition is processed and the case where the pulse wave without differentiation is processed as the extraction condition, there is a concern that the numerical value of the biometric information that can be acquired may vary. To. That is, in order to acquire biometric information with sufficient accuracy by using the pulse wave waveform signal, the pulse wave waveform signal is polygonized from a plurality of biometric information output simultaneously through a plurality of types of extraction conditions and the like. It is necessary to evaluate it.
一方、特許文献1では、1つの脈波の波形信号に対して、複数種類の抽出条件等を経て、複数の生体情報を同時に出力することが開示されていない。このため、特許文献1の開示技術では抽出条件等によって取得できる生体情報にばらつきが生じるので、十分な精度が得られないという問題点があった。
On the other hand,
そこで本発明は、上述した問題点に鑑みて案出されたものであり、その目的とするところは、多角的な評価によって高精度な評価結果を得ることができる生体情報演算システムを提供することにある。 Therefore, the present invention has been devised in view of the above-mentioned problems, and an object of the present invention is to provide a biometric information calculation system capable of obtaining highly accurate evaluation results by multifaceted evaluation. It is in.
第1発明に係る生体情報演算システムは、速度脈波を脈波信号として取得する取得手段と、第1抽出条件を参照し、前記取得手段によって取得された脈波信号に基づく第1評価データを抽出する第1抽出手段と、前記第1抽出条件に紐づく第1処理条件を参照し、前記第1抽出手段によって抽出された第1評価データに対する第1血糖値を、第1データとして取得する第1データ取得手段と、第2抽出条件を参照し、前記取得手段によって取得された脈波信号に基づく第2評価データを抽出する第2抽出手段と、前記第2抽出条件に紐づく第2処理条件を参照し、前記第2抽出手段によって抽出された第2評価データに対する、前記第1データと異なる種類の生体情報からなる第2データを取得する第2データ取得手段と、を備え、前記第1評価データと、前記第2評価データとは、同一の前記脈波信号に基づいて抽出されることを特徴とする。 The biological information calculation system according to the first invention refers to an acquisition means for acquiring a velocity pulse wave as a pulse wave signal and a first extraction condition, and obtains first evaluation data based on the pulse wave signal acquired by the acquisition means. With reference to the first extraction means to be extracted and the first processing condition associated with the first extraction condition, the first blood glucose level for the first evaluation data extracted by the first extraction means is acquired as the first data . The first data acquisition means, the second extraction means for extracting the second evaluation data based on the pulse wave signal acquired by the acquisition means with reference to the second extraction condition, and the second extraction means associated with the second extraction condition. The second data acquisition means for acquiring the second data consisting of biological information of a type different from the first data with respect to the second evaluation data extracted by the second extraction means with reference to the processing conditions is provided. The first evaluation data and the second evaluation data are characterized in that they are extracted based on the same pulse wave signal.
第2発明に係る生体情報演算システムは、第1発明において、前記第1抽出手段は、前記取得手段によって取得された脈波信号を微分しない前記第1抽出条件を参照し、前記第2抽出手段は、前記取得手段によって取得された脈波信号を微分する前記第2抽出条件を参照することを特徴とする。 In the first invention, the biological information calculation system according to the second invention refers to the first extraction condition in which the first extraction means does not differentiate the pulse wave signal acquired by the acquisition means, and the second extraction means. Is characterized by referring to the second extraction condition for differentiating the pulse wave signal acquired by the acquisition means.
第3発明に係る生体情報演算システムは、第1発明又は第2発明において、前記第1データ取得手段は、予め取得された第1評価データと生体情報の実測値との相関関係に基づいて構築された第1検量モデルを用いて、前記第1抽出手段によって抽出された第1評価データに対する第1データを取得する第1回帰分析手段を有し、前記第1回帰分析手段は、前記第1検量モデルを用いて、前記第1抽出手段によって抽出された第1評価データに対する第1血糖値を、前記第1データとして取得することを特徴とする。 The biological information calculation system according to the third invention is constructed in the first invention or the second invention by the first data acquisition means based on the correlation between the first evaluation data acquired in advance and the measured value of the biological information. The first regression analysis means has a first regression analysis means for acquiring the first data for the first evaluation data extracted by the first extraction means by using the first calibration model, and the first regression analysis means is said to be the first. Using a calibration model, the first blood glucose level for the first evaluation data extracted by the first extraction means is acquired as the first data.
第1発明~第3発明によれば、本発明の生体情報演算システムは、第1抽出条件を参照した第1抽出手段によって、脈波信号に基づく第1評価データを抽出し、さらに第1抽出条件に紐づく第1処理条件を参照した第1データ取得手段によって、第1評価データに対する第1データを取得する。また、本発明の生体情報演算システムは、第2抽出条件を参照した第2抽出手段によって、脈波信号に基づく第2評価データを抽出し、さらに第2抽出条件に紐づく第2処理条件を参照した第2データ取得手段によって、第2評価データに対する第2データを取得する。これらによって、入力された1つの脈波信号から、異なる抽出方法及び処理方法を施した複数の生体情報を同時に取得できる。この複数の生体情報によって、一つの脈波信号で生体情報の多角的な評価によって高精度な評価結果を得ることが可能になる。 According to the first to third inventions, the biometric information calculation system of the present invention extracts the first evaluation data based on the pulse wave signal by the first extraction means with reference to the first extraction condition, and further first extracts. The first data for the first evaluation data is acquired by the first data acquisition means referring to the first processing condition associated with the condition. Further, the biometric information calculation system of the present invention extracts the second evaluation data based on the pulse wave signal by the second extraction means with reference to the second extraction condition, and further sets the second processing condition associated with the second extraction condition. The second data for the second evaluation data is acquired by the referenced second data acquisition means. As a result, it is possible to simultaneously acquire a plurality of biometric information to which different extraction methods and processing methods are applied from one input pulse wave signal. With this plurality of biometric information, it is possible to obtain highly accurate evaluation results by multifaceted evaluation of biometric information with one pulse wave signal.
特に、第2発明によれば、本発明の生体情報演算システムは、第1抽出手段で脈波信号を微分することにより、加速度脈波を第1評価データとして抽出する。また、本発明の生体情報演算システムは、第2抽出手段で脈波信号を微分しないことにより、速度脈波を第2評価データとして抽出する。これらによって、一つの脈波信号から加速度脈波と速度脈波を取得できる。加速度脈波と速度脈波から得られる特徴点が異なるため、これらのデータから生体情報を取得することにより、生体情報のより多角的な評価によって高精度な評価結果を得ることが可能になる。 In particular, according to the second invention, the biometric information calculation system of the present invention extracts the acceleration pulse wave as the first evaluation data by differentiating the pulse wave signal by the first extraction means. Further, the biometric information calculation system of the present invention extracts the velocity pulse wave as the second evaluation data by not differentiating the pulse wave signal by the second extraction means. With these, the acceleration pulse wave and the velocity pulse wave can be acquired from one pulse wave signal. Since the feature points obtained from the acceleration pulse wave and the velocity pulse wave are different, by acquiring the biological information from these data, it becomes possible to obtain a highly accurate evaluation result by more multifaceted evaluation of the biological information.
特に、第3発明によれば、第1回帰分析手段は、第1検量モデルを用いて、第1抽出手段によって抽出された第1評価データに対する第1血糖値を、前記第1データとして取得する。これによって、誤検出の少ない速度脈波を第1評価データとして用いることができ、血糖値を算出することができるため、さらに高精度な評価結果を得ることが可能になる。 In particular, according to the third invention, the first regression analysis means uses the first calibration model to acquire the first blood glucose level for the first evaluation data extracted by the first extraction means as the first data. .. As a result, the velocity pulse wave with less erroneous detection can be used as the first evaluation data, and the blood glucose level can be calculated, so that more accurate evaluation results can be obtained.
以下、本発明を適用した生体情報演算システム1について、図面を参照しながら詳細に説明をする。
Hereinafter, the biological
本発明を適用した生体情報演算システム1は、例えば図1に示すような構成により具現化される。この生体情報演算システム1は、電子機器2と、この電子機器2に対して公衆通信網3を介して接続されるサーバ4、センサ5とを備えている。
The biometric
公衆通信網3は、電子機器2と、サーバ4と、センサ5とを通信回線を介して接続されるインターネット網等である。公衆通信網3は、生体情報演算システム1を一定の狭いエリア内で運用する場合には、LAN(Local Area Network)で構成されてもよい。また、公衆通信網3は、いわゆる光ファイバ通信網で構成されてもよい。また、この公衆通信網3は、有線通信網に限定されるものではなく、無線通信網で実現されてもよい。
The
サーバ4は、公衆通信網3を介して送られてきた情報が蓄積されたデータベースで構築されている。また、このサーバ4は、電子機器2からの要求に基づいて、この蓄積した情報を公衆通信網3を介して電子機器2へと送信する。
The
図2は、センサ5の具体的な構成例を示している。センサ5は、取得部50と、通信I/F51(通信I/F27)と、メモリ52と、命令部53とがそれぞれ内部バス54で接続されている。センサ5は、例えば図3のようにリストバンド55に埋め込まれ、リストバンド55は、手首に装着される。それによって、センサ5は、取得部50が手首の脈に接触するので、脈波信号を容易に測定できる。また、センサ5は、衣服に埋め込まれていてもよい。また、センサ5は、ユーザとして、人間に取り付けられることを想定して説明しているが、人間に限らず、ペットや牛、豚などの家畜、又は養殖中の魚類にも使用してもよい。
FIG. 2 shows a specific configuration example of the
脈波信号は、ユーザの速度脈波に対応するデジタル信号又はアナログ信号を示し、例えば縦軸を速度脈波の振幅、横軸を時間とした平面上のグラフに波形としてプロットでき、振幅の大きさの時間変化を可視化した信号として扱うことができる。 The pulse wave signal indicates a digital signal or an analog signal corresponding to the user's velocity pulse wave. For example, the vertical axis can be plotted as a waveform on a plane graph with the amplitude of the velocity pulse wave and the horizontal axis with time, and the amplitude can be large. It can be treated as a signal that visualizes the change over time.
取得部50は、ユーザから脈波信号を測定するための少なくとも一つの測定器を備える。例えば、取得部50は、ファイバブラッググレーティング(FBG)センサなどのひずみセンサと、ジャイロセンサと、脈波信号測定のための1つ以上の電極と、光電容積脈波(PPG)センサと、圧力センサと、光源及び検出器を含んだ光検出モジュールとの内の少なくとも1つを含む測定器を備えている。取得部50は、測定器を通じてユーザと直接インターフェースされて脈波信号を取得することができる。取得部50は、取得した脈波信号を通信I/F51、又はメモリ52へと送信する。
The
通信I/F51は、公衆通信網3を介して取得部50から送信された脈波信号を電子機器2、又はサーバ4に送信する。また、通信I/F51は、公衆通信網3と接続するための回線制御回路や、電子機器2及びサーバ4との間でデータ通信を行うための信号変換回路等が実装されている。通信I/F51は、内部バス54からの各種命令に変換処理を施して、これを公衆通信網3側へ送出するとともに、公衆通信網3からのデータを受信した場合には、これに所定の変換処理を施して内部バス54へ送信する。
The communication I /
メモリ52は、取得部50から送信された脈波信号を保存する。メモリ52は、公衆通信網3を介して接続される他の末端装置から命令を受けることにより、保存した脈波信号を通信I/F51へ送信する。
The
命令部53は、取得部50へ脈波信号の取得の命令を入力するための操作ボタンやキーボード等で構成される。この命令部53は、脈波信号の取得の命令を受け付けた場合に、これを取得部50に通知する。この通知を受けた取得部50は、脈波信号を取得する。
The
電子機器2は、例えば、パーソナルコンピュータ(PC)、携帯電話、スマートフォン、タブレット型端末、ウェアラブル端末等であり、少なくともユーザの操作に基づいて、公衆通信網3を介して通信可能なデバイスである。なお、この電子機器2は、センサ5を内蔵し、センサ5に対して公衆通信網3を介した通信を行わない機器であってもよい。以下の例では、この電子機器2に、PCを適用する場合を例にとり説明をする。
The
図4は、電子機器2の具体的な構成例を示している。この電子機器2は、ROM(Read Only Memory)20と、RAM(Random Access Memory)21と、CPU(Central Processing Unit)22と、操作部23と、出力I/F16と、記憶部24と、データ入出力部25と、通信I/F27とが内部バス26にそれぞれ接続されている。また、出力I/F16には、表示部28が接続されている。
FIG. 4 shows a specific configuration example of the
また、電子機器2は、CPU22が、RAM21を作業領域として、記憶部24等に記憶された生体情報演算プログラムを実行することにより、脈波信号から生体情報を取得する。
Further, in the
ROM20は、電子機器2全体のハードウェア資源を制御するためのプログラムが格納されている。
The
RAM21は、電子機器2全体のハードウェア資源を制御するときの各種命令を一時的に記憶するデータの蓄積や展開等に使用する作業領域である。
The
CPU22は、内部バス26を介して制御信号を送信することにより、電子機器2内に実装された各構成要素を制御するためのいわゆる中央演算ユニットである。また、このCPU22は、操作部23を介したユーザの操作に応じて各種制御用の指令を、内部バス26を介して伝達する。
The
操作部23は、各種制御用の指令を入力するためのマウスやキーボード、タッチパネル、及び操作ボタン等のデバイスで構成されている。操作部23は、ユーザが実際に取得したい生体情報に関する情報が入力される他、生体情報演算プログラムを実行するための実行命令がユーザから入力される。操作部23は、上述した実行命令がユーザにより入力された場合には、これをCPU22に通知する。この通知を受けたCPU22は、生体情報演算プログラムを記憶部24から読み出して実行する。
The
出力I/F16は、各種情報の表示を制御する。出力I/F16はCPU22による制御に基づいて、表示画像を作り出すグラフィックコントローラにより構成されている。この出力I/F16に接続される表示部28は、例えば、液晶ディスプレイ(LCD)等によって実現される。
The output I /
記憶部24は、ハードディスク等に代表される記録媒体であり、実行すべき生体情報取得プログラムが格納されている。記憶部24は、ハードディスクで構成される場合において、CPU22による制御に基づき、各アドレスに対して所定の情報が書き込まれるとともに、必要に応じてこれが読み出される。記憶部24はCPU22により読み出されて生体情報取得プログラムを出力されることになる。
The
通信I/F27は、公衆通信網3と接続するための回線制御回路や、他の端末装置との間でデータ通信を行うための信号変換回路等が実装されている。通信I/F27は、内部バス26からの各種命令に変換処理を施して、これを公衆通信網3側へ送出するとともに、公衆通信網3からのデータを受信した場合には、これに所定の変換処理を施して、内部バス26、又はCPU22へ送信する。
The communication I /
データ入出力部25は、外部から電子機器2内へデータを入力し、又は電子機器2において生成されたデータを外部へ出力する。データ入出力部25は、USB(Universal Serial Bus)メモリや、記録媒体との間でデータを入力するためのインターフェースとして構成されている。
The data input /
本発明を適用した生体情報演算システム1は、電子機器2内にインストールされた生体情報演算プログラムを介して実行していくこととなる。即ち、ユーザは、電子機器2、又はセンサ5を操作し、電子機器2にインストールされている生体情報演算プログラムを通じて、脈波信号から生体情報を取得する。
The biometric
図5は、電子機器2の生体情報演算プログラムを実現するためのシーケンスを示している。電子機器2は、脈波信号取得部60と、脈波信号取得部60に接続された第1抽出部61及び第2抽出部63と、第1抽出部61に接続された第1データ取得部62と、第2抽出部63に接続された第2データ取得部64と、第1データ取得部62及び第2データ取得部64とに接続された最適血糖値算出部65と、を備えている。
FIG. 5 shows a sequence for realizing the biometric information calculation program of the
脈波信号取得部60は、公衆通信網3を介してセンサ5、サーバ4及び他の電子機器から送信された脈波信号を取得する。脈波信号取得部60は、取得した脈波信号を第1抽出部61及び第2抽出部63へ出力する。
The pulse wave
第1抽出部61は、第1抽出条件を参照し、脈波信号取得部60から入力された脈波信号に基づく第1評価データを抽出する。第1抽出部61は、抽出した第1評価データを第1データ取得部62へ出力する。
The
第1データ取得部62は、第1抽出条件に紐づく第1処理条件を参照し、第1抽出部61から入力された第1評価データを処理し、第1データを取得する。第1データ取得部62は、取得した第1データを最適血糖値算出部65へ出力する。
The first
第2抽出部63は、第2抽出条件を参照し、脈波信号取得部60から入力された脈波信号に基づく第2評価データを抽出する。第2抽出部63は、抽出した第2評価データを第2データ取得部64へ出力する。
The
第2データ取得部64は、第2抽出条件に紐づく第2処理条件を参照し、第2抽出部63から入力された第2評価データを処理し、第2データを取得する。第2データ取得部は、取得した第2データを最適血糖値算出部65へ出力する。
The second
最適血糖値算出部65は、第1データ取得部62から入力された第1データと第2データ取得部64から入力された第2データとに基づき、最適値となる生体情報を算出する。
The optimum blood glucose
また、最適血糖値算出部65は、電子機器2に必ずしも備えられているわけではなく、第1データおよび第2データを生体情報として出力してもよい。
Further, the optimum blood glucose
図6は、第1抽出部61及び第1データ取得部62の具体的な構成例を示している。第1抽出部61は、フィルタ処理部610と、フィルタ処理部610に接続された微分部611と、フィルタ処理部610に接続されたピーク位置算出部614と、微分部611に接続された分割部612と、分割部612に接続された規格化部613と、ピーク位置算出部614に接続されたピーク間隔平均算出部615と、ピーク位置算出部614に接続されたピーク間隔プロット部616と、ピーク位置算出部614に接続されたフーリエ変換部617と、フーリエ変換部617に接続された最大周波数検出部618と、を備えている。
FIG. 6 shows a specific configuration example of the
フィルタ処理部610は、取得した脈波信号にフィルタリング処理を施す。フィルタ処理部610は、フィルタリングに例えば0.5~5Hzのバンドパスフィルタを用いるが、この限りではない。また、フィルタ処理部610は、第1抽出条件を参照し、取得した脈波信号から第1評価データを抽出する抽出方法を決定する。フィルタ処理部610は、フィルタリング処理された脈波信号を決定された抽出方法に基づき、微分部611と、ピーク位置算出部614と、フーリエ変換部617のいずれか少なくとも一つに出力する。
The
第1抽出部61は、一つの第1評価データを取得するために第1抽出部61に備えられた抽出方法の必ずしもすべてを使用するわけではなく、フィルタ処理部610によって決められた少なくとも一つの抽出方法で脈波信号から第1評価データを抽出する。
The
微分部611は、フィルタ処理部610から入力された脈波信号を微分する。微分部611は、フィルタ処理部610によって、微分処理が必要と判断された場合、入力された脈波信号に微分処理を行う。微分部611は、処理した脈波信号を分割部612に出力する。
The
分割部612は、微分部611から入力された複数の波形信号のそれぞれを、整数周期分の分割波形データに分割する。本実施形態では分割部612は、整数周期は1周期としているが、複数周期にしてもよい。分割部612は、分割波形データを規格化部613へ出力する。
The
規格化部613は、分割部612から入力された複数の分割波形信号の時間幅を統一するために規格化をし、複数の分割波形信号の平均となる平均波形信号を取得し、平均波形信号の振幅の最大値を1、最小値を0とした規格化を行う。規格化部613は、規格化された平均波形信号を回帰分析部620へ出力する。
The
規格化部613は、分割波形信号の時間幅を統一するために分割部612から入力された複数の分割波形信号を一定の時間幅又は一定のサンプリング数でトリミングを行ってもよい。規格化部613は、時間幅を統一するための処理方法がフィルタ処理部610によって決められる。
The
規格化部613は、平均波形信号を取得するとき、複数の分割波形信号を必要とするが、必要な分割波形信号の数はフィルタ処理部610によって決められる。
The
ピーク位置算出部614は、フィルタ処理部610から入力された脈波信号のピーク位置と、隣り合うピーク同士の距離であるピーク間隔を算出する。ピーク位置算出部614は、抽出方法に基づき、算出したピーク間隔をピーク間隔平均算出部615と、ピーク間隔プロット部616と、フーリエ変換部617の少なくともいずれか一つへと出力する。
The peak
ピーク間隔平均算出部615は、ピーク位置算出部614から入力されたピーク間隔の平均を算出し、ピーク間隔の平均値を測定器のサンプリングレートで割り、秒数に変換する。ピーク間隔平均算出部615は、秒数に変換した平均値を脈拍処理部621へ出力する。
The peak interval
ピーク間隔プロット部616は、ピーク位置算出部614から入力されたピーク間隔を横軸とし、上述したピーク間隔の隣のピーク間隔を縦軸としたグラフをプロットし、脈波のピーク間隔のポアンカレプロットを取得する。ピーク間隔プロット部616は、脈波のピーク間隔のポアンカレプロットをストレスプロット処理部622へ出力する。
The peak
フーリエ変換部617は、ピーク位置算出部614からピーク間隔を入力された場合は、ピーク間隔を時系列データに変換したデータをフーリエ変換する。また、フーリエ変換部617は、フィルタ処理部610で処理された脈波信号を入力として、フーリエ変換を行ってもよい。フーリエ変換部617は、抽出方法に基づき、フーリエ変換した信号を最大周波数検出部618と、ストレスフーリエ処理部623の少なくともいずれか一つに出力する。
When the peak interval is input from the peak
最大周波数検出部618は、フーリエ変換部617から入力された信号を0.15~0.35Hzの間で最大値を示す周波数である最大周波数を検出する。最大周波数検出部618は、検出した最大周波数を呼吸数処理部624へと出力する。
The maximum
第1データ取得部62は、規格化部613に接続された回帰分析部620と、ピーク間隔平均算出部615に接続された脈拍処理部621と、ピーク間隔プロット部616に接続されたストレスプロット処理部622と、フーリエ変換部617に接続されたストレスフーリエ処理部623と、最大周波数検出部618に接続された呼吸数処理部624と、を備える。
The first
第1データ取得部62は、例えば処理の方法として、生体情報の実測値と予め取得された脈波信号との相関関係に基づいて構築された検量モデルを用いて、第1評価データから第1データを取得する。
The first
回帰分析部620は、規格化部613から入力された規格化された平均波形信号から、検量モデルに基づいて、第1データとして例えば血糖値、血圧、血中の酸素飽和度、血中二酸化炭素濃度などを取得する。
The
回帰分析部620は、あらかじめ構築しておいた汎用的に利用できる検量モデル用いて、第1データ取得部62に記憶させて第1データを取得することもできる。
The
脈拍処理部621は、ピーク間隔平均算出部615から入力されたピーク間隔の平均値をサンプリングレートで割り、秒数に変換する。脈拍処理部621は、60秒を算出した秒数で割り、一分当たりの脈拍数(bpm)を算出し、第1データとして脈拍数を取得する。
The
ストレスプロット処理部622は、ピーク間隔プロット部616から入力されたポアンカレプロットから、第1データとして例えばストレス度を取得する。ストレスプロット処理部622は、ポアンカレプロットの分散値を計算し、分散値の大きさからストレス度を推定する。
The stress
ストレスフーリエ処理部623は、フーリエ変換部617から入力された信号の積分比から、第1データとして例えばストレス度を取得する。具体的な方法として、ストレスフーリエ処理部623は、例えばフーリエ変換されたピーク間隔の時系列データから自己回帰モデル(autoregressive model)を用いて、PSD(power spectral density)を計算し、パワースペクトルの0.5Hzから0.15Hzまでの領域を高周波LF(Low Frequency)とし、0.15Hzから0.40Hzまでの領域を低周波HF(Hi Frequency)の強度をそれぞれ合計した積分値の比によってストレス度を決定する。
The stress
呼吸数処理部624は、最大周波数検出部618から入力された最大値を示す周波数に60秒をかけて一分あたりの呼吸数(bpm)に換算し、第1データとして呼吸数を取得する。
The respiratory
図7は、第2抽出部63及び第2データ取得部64の具体的な構成例を示している。第2抽出部63は、フィルタ処理部630と、フィルタ処理部630に接続された微分部631と、フィルタ処理部630に接続されたピーク位置算出部634と、微分部631に接続された分割部632と、分割部632に接続された規格化部633と、ピーク位置算出部634に接続されたピーク間隔平均算出部635と、ピーク位置算出部634に接続されたピーク間隔プロット部636と、ピーク位置算出部634に接続されたフーリエ変換部637と、フーリエ変換部637に接続された最大周波数検出部638と、を備える。第2抽出部63は、第1抽出部61と同じもので構成されているが、第1抽出部61において、第1抽出条件を参照していたものが、第2抽出部63の場合は、第1抽出条件ではなく、第2抽出条件を参照している。
FIG. 7 shows a specific configuration example of the
フィルタ処理部630は、取得した脈波信号にフィルタリング処理を施す。フィルタ処理部630は、フィルタリングに例えば0.5~5Hzのバンドパスフィルタを用いるが、この限りではない。また、フィルタ処理部630は、第2抽出条件を参照し、取得した脈波信号から第1評価データを抽出する方法を決定する。フィルタ処理部630は、フィルタリング処理された脈波信号を決定された抽出方法に基づき、微分部631と、ピーク位置算出部634と、フーリエ変換部637のいずれか少なくとも一つに出力する。
The
第2抽出部63は、一つの第2評価データを取得するために第2抽出部63に備えられた抽出方法を必ずしもすべて使用するわけではなく、フィルタ処理部630によって決められた少なくとも一つの抽出方法で脈波信号から第2評価データを抽出する。
The
微分部631は、フィルタ処理部630から入力された脈波信号を微分する。微分部631は、フィルタ処理部630によって、微分処理が必要と判断された場合、入力された脈波信号に微分処理を行う。微分部631は、微分処理した脈波信号を分割部632に出力する。
The
分割部632は、微分部631から入力された複数の波形信号のそれぞれを、整数周期分の分割波形データに分割する。本実施形態では分割部632は、整数周期は1周期としているが、複数周期にしてもよい。分割部632は、分割波形データを規格化部633へ出力する。
The
規格化部633は、分割部632から入力された複数の分割波形信号の時間幅を統一するために規格化をし、複数の分割波形信号の平均となる平均波形信号を取得し、平均波形信号の振幅の最大値を1、最小値を0とした規格化を行う。規格化部633は、規格化された平均波形信号を回帰分析部640へ出力する。
The
規格化部633は、分割波形信号の時間幅を統一するために分割部632から入力された複数の分割波形信号を一定の時間幅又は一定のサンプリング数でトリミングを行ってもよい。規格化部633は、時間幅を統一するための処理方法がフィルタ処理部630によって決められる。
The
規格化部633は、平均波形信号を取得するとき、複数の分割波形信号を必要とするが、必要な分割波形信号の数はフィルタ処理部630によって決められる。
The
ピーク位置算出部634は、フィルタ処理部630から入力された脈波信号のピーク位置と、隣り合うピーク同士の距離であるピーク間隔を算出する。ピーク位置算出部634は、抽出方法に基づき、算出したピーク間隔をピーク間隔平均算出部635と、ピーク間隔プロット部636と、フーリエ変換部637の少なくともいずれか一つへと出力する。
The peak
ピーク間隔平均算出部635は、ピーク位置算出部634から入力されたピーク間隔の平均を算出し、ピーク間隔の平均値を測定器のサンプリングレートで割り、秒数に変換する。ピーク間隔平均算出部635は、秒数に変換した平均値を脈拍処理部641へ出力する。
The peak interval
ピーク間隔プロット部636は、ピーク位置算出部634から入力されたピーク間隔を横軸とし、上述したピーク間隔の隣のピーク間隔を縦軸としたグラフをプロットし、脈波のピーク間隔のポアンカレプロットを取得する。ピーク間隔プロット部636は、脈波のピーク間隔のポアンカレプロットをストレスプロット処理部642へ出力する。
The peak
フーリエ変換部637は、ピーク位置算出部634からピーク間隔を入力された場合は、ピーク間隔を時系列データに変換したデータをフーリエ変換する。また、フーリエ変換部637は、フィルタ処理部630で処理された脈波信号を入力として、フーリエ変換を行ってもよい。フーリエ変換部637は、抽出方法に基づき、フーリエ変換した信号を最大周波数検出部638と、ストレスフーリエ処理部643の少なくともいずれか一つに出力する。
When the peak interval is input from the peak
最大周波数検出部638は、フーリエ変換部637から入力された信号を0.15~0.35Hzの間で最大値を示す周波数である最大周波数を検出する。最大周波数検出部638は、検出した最大周波数を呼吸数処理部644へと出力する。
The maximum
第2データ取得部64は、規格化部633に接続された回帰分析部640と、ピーク間隔平均算出部635に接続された脈拍処理部641と、ピーク間隔プロット部636でプロットしたグラフからストレス度を取得するストレスプロット処理部642と、フーリエ変換部637に接続されたストレスフーリエ処理部643と、最大周波数検出部638に接続された呼吸数処理部644と、を備える。
The second
第2データ取得部64は、第2抽出条件に紐づく第2処理条件を参照し、第2抽出部63によって抽出された第2評価データを処理し、第2データを取得する。
The second
第2データ取得部64は、例えば処理の方法として、生体情報の実測値と予め取得された脈波信号との相関関係に基づいて構築された検量モデルを用いて構築し、第2評価データから第2データを取得する。
The second
回帰分析部640は、規格化部633から入力された規格化された平均波形信号から、検量モデルに基づいて、第2データとして例えば血糖値、血圧、血中の酸素飽和度、血中二酸化炭素濃度などを取得する。
The
回帰分析部640は、あらかじめ構築しておいた汎用的に利用できる検量モデルを汎用データとして、第2データ取得部64に記憶させて第2データを取得することもできる。
The
脈拍処理部641は、ピーク間隔平均算出部635から入力されたピーク間隔の平均値をサンプリングレートで割り、秒数に変換する。脈拍処理部641は、60秒を算出した秒数で割り、一分当たりの脈拍数(bpm)を算出し、第2データとして脈拍数を取得する。
The
ストレスプロット処理部642は、ピーク間隔プロット部636から入力されたポアンカレプロットから、第2データとして例えばストレス度を取得する。ストレスプロット処理部642は、ポアンカレプロットの分散値を計算し、分散値の大きさからストレス度を推定する。
The stress
ストレスフーリエ処理部643は、フーリエ変換部637から入力された信号の積分比から、第2データとして例えばストレス度を取得する。
The stress
呼吸数処理部644は、最大周波数検出部638から入力された最大値を示す周波数に60秒をかけて一分あたりの呼吸数(bpm)に換算し、第2データとして呼吸数を取得する。
The respiratory
生体情報演算システム1によって取得される第1データ及び第2データは、例えば血圧、血糖値、血中の酸素飽和度、血中二酸化炭素濃度、脈拍数、呼吸数、ストレス度、血管年齢、糖尿病の程度などのうち少なくとも一つを含める。
The first data and the second data acquired by the biological
次に、本実施形態における生体情報演算システム1の一例について説明する。図8は、本実施形態における生体情報演算システム1の一例を示すフローチャートである。
Next, an example of the biometric
まず脈波信号取得ステップS10において、生体情報演算システム1は、取得部50が脈波信号を取得し、内部バス54を介して、通信I/F51へ脈波信号を出力する。またこの時、取得部50が取得した脈波信号の代わりに、生体情報演算システム1は、メモリ52に記録されていた脈波信号を、内部バス54を介して、通信I/F51へ脈波信号を出力してもよい。
First, in the pulse wave signal acquisition step S10, the biometric
脈波信号の具体的な取得方法として、FBGセンサを用いて測定する場合を例として説明する。 As a specific acquisition method of the pulse wave signal, a case of measuring using an FBG sensor will be described as an example.
FBGセンサは、1本の光ファイバ内に所定間隔をあけて回折格子構造を形成したである。本実施形態では、FBGセンサとして、センサ部分の長さ10mm、波長分解能±0.1pm、波長範囲1550±0.5nmのものを使用し、ファイバの直径は145μm、コア径10.5μmとした。しかし、これらの設定は、FBGセンサの一例であり、本発明においてはこの限りではない。このセンサを手首に取り付け、FBGセンサを皮膚に接触させた状態で測定をする。 The FBG sensor has a diffraction grating structure formed in one optical fiber at predetermined intervals. In the present embodiment, the FBG sensor has a sensor portion length of 10 mm, a wavelength resolution of ± 0.1 pm, and a wavelength range of 1550 ± 0.5 nm, and has a fiber diameter of 145 μm and a core diameter of 10.5 μm. However, these settings are examples of FBG sensors, and are not limited to this in the present invention. This sensor is attached to the wrist, and the measurement is performed with the FBG sensor in contact with the skin.
光ファイバに用いる光源として、波長範囲1525~1570nmのASE(Amplified Spontaneous Emission)光源を使用してもよい。光源からの出射光は、サーキュレータを介してFBGセンサに入射させる。FBGセンサからの反射光は、サーキュレータを介してマッハツェンダー干渉計に導き、マッハツェンダー干渉計からの出力光を、光検出器によって検知する。マッハツェンダー干渉計は、ビームスプリッタにより光路差のある2つの光路に分離し、再びビームスプリッタにより一つに重ね合わせて干渉光を作り出すためのものである。光路差をつけるため、本実施形態では一方の光ファイバの長さを長くしている。コヒーレント光は、光路差に応じて干渉縞が生じるから、干渉縞のパターンを測定することによって、FBGセンサに生じた歪の変化、すなわち脈波を検知することができる。 As the light source used for the optical fiber, an ASE (Amplified Spontaneous Emission) light source having a wavelength range of 1525 to 1570 nm may be used. The light emitted from the light source is incident on the FBG sensor via the circulator. The reflected light from the FBG sensor is guided to the Mach-Zehnder interferometer via the circulator, and the output light from the Mach-Zehnder interferometer is detected by the optical detector. The Mach-Zehnder interferometer is for splitting into two optical paths having optical path differences by a beam splitter and superimposing them on one again by a beam splitter to generate interferometric light. In this embodiment, the length of one of the optical fibers is increased in order to make an optical path difference. Since coherent light produces interference fringes according to the optical path difference, it is possible to detect a change in strain generated in the FBG sensor, that is, a pulse wave, by measuring the pattern of the interference fringes.
FBGセンサの歪み量を検出して、脈波の波形を検出するシステムである光ファイバセンサシステムは、FBGセンサに入射させる光源の他に、広い帯域のASE光源とする手段、サーキュレータ、マッハツェンダー干渉計、ビームスプリッタといった光学系や、光検出器が備える受光センサや、波長シフト量を解析する解析手段を含む。光ファイバセンサシステムは、使用するFBGセンサの特性に応じて光源や帯域光を選択して使用することができ、検波方法等の解析手段についても種々の方法を採用することができる。また、本発明は、光ファイバセンサシステムの機能や方式が限定されるものではない。 The optical fiber sensor system, which is a system that detects the amount of distortion of the FBG sensor and detects the waveform of the pulse wave, has a means for using a wide band ASE light source, a circulator, and Mach Zender interference in addition to the light source incident on the FBG sensor. It includes an optical system such as a meter and a beam splitter, a light receiving sensor provided in the optical detector, and an analysis means for analyzing the wavelength shift amount. The optical fiber sensor system can be used by selecting a light source or band light according to the characteristics of the FBG sensor to be used, and various methods can be adopted as analysis means such as a detection method. Further, the present invention is not limited to the functions and methods of the optical fiber sensor system.
これらの方法によって、取得部50が脈波信号を取得し、内部バス54を介して、通信I/F51へ脈波信号を出力する。
By these methods, the
次に取得部50から脈波信号を入力された通信I/F51は、公衆通信網3を介して、脈波信号を脈波信号取得部60へ送信する。また、この時、取得部50が取得した脈波信号の代わりに、サーバ4に保存された脈波信号を脈波信号取得部60へ送信してもよい。
Next, the communication I /
次に公衆通信網3を介して、脈波信号を送信された脈波信号取得部60は、脈波信号を第1抽出部61及び第2抽出部63へ出力する。
Next, the pulse wave
次に第1抽出条件決定ステップS11において、第1抽出部61は、脈波信号取得部60から入力された脈波信号をフィルタ処理部610へ入力する。その後、フィルタ処理部610は入力された脈波信号に施す抽出方法を、第1抽出条件を参照して、決定する。
Next, in the first extraction condition determination step S11, the
フィルタ処理部610は、第1抽出条件の中から、取得した脈波信号の状態、取得したい生体情報、外部要因によって、第1抽出部61が脈波信号に施す抽出方法を決定する。外部要因は、例えばユーザの年齢や、性別や、病歴や、生活習慣や、投薬情報や、動脈硬化の程度や、健康状態や、遺伝情報などのユーザの情報、温度や、湿度などの環境情報などのうち少なくとも一つを含める。例えば、取得したい生体情報として、血糖値を取得したい場合、フィルタ処理部610は、脈波信号がフィルタ処理部610で処理を施された後、微分部611へと出力されるように第1抽出部61に命令をする。第1抽出条件は、脈波信号から第1評価データを抽出するための抽出方法の一覧を含めるデータ群である。第1抽出条件のデータ群は複数の抽出方法を含めてもよいし、予め決められた単一の抽出方法を含めていてもよい。
The
次に第1抽出ステップS12において、第1抽出部61は、脈波信号取得部60から入力された脈波信号から第1評価データを抽出する。第1評価データとは、第1データ取得部62で生体情報を取得するために第1抽出部61によって、脈波信号から抽出された波形データである。第1評価データは、例えば脈波信号が、フィルタリング処理、微分処理、規格化処理、及び平均化処理の少なくとも何れかによって1周期分等の波形に処理された波形データである。本実施形態では、例として血糖値を取得するための生体情報演算システム1を説明する。
Next, in the first extraction step S12, the
まずフィルタ処理部610は、脈波信号取得部60から入力された脈波信号にフィルタリング処理を施した後、例えば血糖値を生体情報として取得する場合、第1抽出条件決定ステップS11で決定した抽出方法に基づいて、脈波信号を微分部611へと出力する。
First, the
次に微分部611は、フィルタ処理部610から入力された脈波信号に微分をするか、もしくはしないかを、第1抽出条件決定ステップS11で決定した抽出方法に基づいて判断し、処理を施した後、脈波信号を分割部612へと出力する。
Next, the
第1抽出条件決定ステップS11で決定した抽出方法に基づいて、微分部611が脈波信号を微分するかしないかを判断する理由は、脈波信号を微分するかしないかで得られる第1評価データの特徴に差が生じ、取得したい第1データによって、それに適した第1評価データを得るためである。また、「脈波信号を微分する」とは、脈波信号を加速度脈波として抽出することであり、「脈波信号を微分しない」とは、脈波信号を速度脈波として抽出することである。
The reason why the
次に分割部612は、微分部611から入力された複数の波形信号のそれぞれを、平均化するために1周期分の分割波形データに分割する。その後、分割部612は、分割波形データを規格化部613へ出力する。
Next, the
規格化部613は、分割部612から入力された複数の分割波形信号の時間幅を統一するために横軸の規格化をし、複数の分割波形信号の平均となる平均波形信号を取得し、平均波形信号の最大値を1、最小値を0とした縦軸の規格化を行う。この時、第1抽出部61は平均波形信号を第1評価データとして取得する。その後、規格化部613は、規格化された平均波形信号を第1データ取得部62へ出力する。
The
規格化部613で時間幅を統一するために横軸の規格化処理を行う理由は、脈波の終端側で差が大きく表れることから、この部分を削除し、脈波の本体部分を解析対象とするためである。また平均波形信号の最大値を1、最小値を0とした縦軸の規格化を行う理由は、FBGセンサを測定部位に取り付ける際の押しつけ圧力のばらつきや、測定時にFBGセンサが位置ずれすることによる測定データのばらつきを平均化するし、測定時のばらつきに起因するノイズを抑え、脈波信号と生体情報の実測値の相関関係の精度を向上させるためである。
The reason why the normalization process on the horizontal axis is performed in the
次に生体情報演算システム1は、第1データ取得ステップS13へ移行し、第1抽出条件に紐づく第1処理条件を参照し、第1抽出部61から入力された第1評価データを第1データ取得部62で処理し、第1データを取得する。
Next, the biometric
第1処理条件は、第1抽出部61から入力された第1評価データに、第1データ取得部62で施す処理の第1抽出条件に紐づけられた方法を含めるデータ群である。第1データ取得部62は第1処理条件の中から処理方法を決定する。第1処理条件のデータ群は複数の処理方法を含めてもよいし、予め決められた単一の処理方法を含めていてもよい。
The first processing condition is a data group including the method associated with the first extraction condition of the processing performed by the first
例えば、上述した抽出方法で取得した平均波形信号を第1評価データとして、第1データ取得部62で処理する場合、第1データ取得部62は、平均波形信号を回帰分析部620へ出力するように決定し、さらに回帰分析部620で平均波形信号に施す処理方法を決定する。
For example, when the average waveform signal acquired by the above-mentioned extraction method is processed by the first
第1抽出部61から平均波形信号を入力された回帰分析部620は、実測値と脈波信号の相関関係を示した検量モデルを用いて、平均波形信号から例えば血糖値を取得し、第1データとして出力する。
The
検量モデルは、予め測定された平均波形脈波を説明変数とし、生体情報の実測値を目的変数として、例えば回帰分析により解析し、その解析結果に基づいて構築されたものである。検量モデルは、あらかじめ構築しておいた汎用的に利用できる検量モデルを記憶部24に記憶させて、第1データを測定することもできる。検量モデルの構築は、例えば定期的にキャリブレーションする場合や、ユーザが変わったときに構築し直しするといった場合に必要になることがある。
The calibration model is constructed based on the analysis result, for example, by performing regression analysis using the measured average waveform pulse wave measured in advance as the explanatory variable and the measured value of the biological information as the objective variable. As the calibration model, a calibration model that has been constructed in advance and can be used for general purposes can be stored in the
検量モデルを構築するための回帰分析の方法として、例えばPLS回帰分析、クラス毎に主成分分析を行って主成分モデルとして検量モデルを得るSIMCA(Soft Independent Modeling of Class Analogy)法を利用した回帰分析、AIを用いたニューラルネットワークなどによる回帰分析方法などの少なくともひとつの方法を含める。 As a method of regression analysis for constructing a calibration model, for example, PLS regression analysis, regression analysis using the SIMCA (Soft Independent Modeling of Class Analogy) method for obtaining a calibration model as a principal component model by performing principal component analysis for each class. , Includes at least one method, such as a regression analysis method using a neural network using AI.
また、例えば血中二酸化炭素濃度などの異常値の値を観測しにくく、異常値のデータが集めにくい生体情報を推定する際に、回帰分析部620は第1抽出部61から入力された平均波形信号と、検量モデルとの乖離度から、生体情報の異常値を推定してもよい。
Further, when estimating biological information in which it is difficult to observe an abnormal value such as a blood carbon dioxide concentration and it is difficult to collect abnormal value data, the
回帰分析部620は複数の検量モデルを有し、入力された平均波形データに対して、どの検量モデルを用いるかは、第1データ取得部62において第1処理条件を参照して決定した処理方法によって決定される。例えばフィルタ処理部610で、取得したい生体情報として、血糖値が選択された場合、第1データ取得部62は、第1抽出部61から入力された平均波形信号を回帰分析部620へ出力し、回帰分析部620で、予め測定された脈波の波形データを説明変数とし、血糖値の実測値を目的変数として、回帰分析により解析し、その解析結果に基づいて構築された検量モデルを用いて、入力された平均波形データから第1データを取得するように決定する。
The
また、例えばフィルタ処理部610が、外部要因として、ユーザの年齢に基づいて、信号の抽出条件を決定した場合、第1データ取得部62は、第1抽出条件に紐付けられた第1処理条件を参照することで、第1抽出部61から入力された平均波形信号を回帰分析部620へ出力し、回帰分析部620で、予め測定された脈波の波形データを説明変数とし、ユーザの年齢に近いユーザの血糖値の実測値を目的変数として、回帰分析により解析し、その解析結果に基づいて構築された検量モデルを用いて、入力された平均波形データから第1データを取得するように決定する。これらによって、抽出方法に適した処理方法の決定が可能となる。
Further, for example, when the
次に第2抽出部63は、第2抽出条件決定ステップS14において、脈波信号取得部60から入力された脈波信号をフィルタ処理部630へ入力する。その後、フィルタ処理部630は入力された脈波信号に施す抽出方法を、第2抽出条件を参照して、決定する。
Next, in the second extraction condition determination step S14, the
第2抽出条件は、脈波信号から第2評価データを抽出するための抽出方法の一覧を含めるデータ群である。第2抽出条件のデータ群は複数の抽出方法を含めてもよいし、予め決められた単一の抽出方法を含めていてもよい。また、第2抽出条件は、第1抽出条件と同じものであってもよい。 The second extraction condition is a data group including a list of extraction methods for extracting the second evaluation data from the pulse wave signal. The data group of the second extraction condition may include a plurality of extraction methods, or may include a single predetermined extraction method. Further, the second extraction condition may be the same as the first extraction condition.
フィルタ処理部630は、第2抽出条件の中から、取得した脈波信号の状態、取得したい生体情報、外部要因によって、第2抽出部63が脈波信号に施す抽出方法を決定する。外部要因は、例えばユーザの年齢や、性別や、病歴や、生活習慣や、健康状態や、遺伝情報などのユーザの情報、温度や、湿度などの環境情報などのうち少なくとも一つを含める。また、フィルタ処理部630は、第1データ取得部62で取得した第1データの内容によって、第2抽出部63が脈波信号に施す抽出方法を決定してもよい。例えば、第1データ取得部62で取得した第1データの血糖値の精度が低い場合、フィルタ処理部630は、第1抽出部61で行わなかった抽出方法を決定する。また、例えば第1データとして血糖値を取得した場合、フィルタ処理部630は、血圧を第2データとして取得するために、規格化部633で、分割波形信号の時間幅を統一するために一定のサンプリング数でトリミングを行う抽出条件を決定してもよい。これによって、第1データの変化や特徴に合わせた第2データの取得が可能となり、高精度な評価結果をえることができる。また、フィルタ処理部630は、第1データ取得部62で取得した第1データの内容によって、第2データ取得部64が第2評価データに施す処理方法を決定してもよい。具体的な方法としては、後述する第2データ取得部64が、第1データ取得部62で取得した第1データの内容によって、第2評価データに施す処理方法を決定する方法と同様なものでもよい。
The
また、フィルタ処理部630は、第1データ取得部62で取得した第1評価データの抽出方法によって、第2抽出部63が脈波信号に施す抽出方法を決定してもよい。
例えば、第1評価データを取得するときに、微分部611で微分を行った場合、フィルタ処理部630は、微分部631で微分を行わないように決定する。これらによって、第1データの変化に伴い、最適な第2データを取得することが可能となる。
Further, the
For example, when the
次に第2抽出部63は、第2抽出ステップS15において、脈波信号取得部60から入力された脈波信号から第2評価データを抽出する。第2評価データとは、第2データ取得部64で処理をするために第2抽出部63によって、脈波信号から抽出された波形データである。第2評価データは、例えば脈波信号が、フィルタリング処理、微分処理、規格化処理、及び平均化処理の少なくとも何れかによって1周期分等の波形に処理された波形データである。
Next, in the second extraction step S15, the
まずフィルタ処理部630は、脈波信号取得部60から入力された脈波信号にフィルタリング処理を施した後、例えば血糖値を生体情報として取得する場合、第2抽出条件決定ステップS14で決定した抽出方法に基づいて、脈波信号を微分部631へと出力する。
First, the
次に微分部631は、フィルタ処理部630から入力された脈波信号に微分をするか、もしくはしないかを、第2抽出条件決定ステップS14で決定した抽出方法に基づいて判断し、処理を施した後、脈波信号を分割部632へと出力する。
Next, the
次に分割部632は、微分部631から入力された複数の波形信号のそれぞれを、1周期分の分割波形データに分割する。その後、分割部632は、分割波形データを規格化部633へ出力する。
Next, the
規格化部633は、分割部632から入力された複数の分割波形信号の時間幅を統一するために規格化をし、複数の分割波形信号の平均となる平均波形信号を取得し、平均波形信号の最大値を1、最小値を0とした規格化を行う。この時、第2抽出部63は平均波形信号を第2評価データとして取得する。その後、規格化部633は、規格化された平均波形信号を第2データ取得部64へ出力する。
The
次に生体情報演算システム1は、第2データ取得ステップS16へ移行し、第2抽出条件に紐づく第2処理条件を参照し、第2抽出部63から入力された第2評価データを第2データ取得部64で処理し、第2データを取得する。
Next, the biometric
第2処理条件は、第2抽出部63から入力された第2評価データに、第2データ取得部64で施す処理の第2抽出条件に紐づけられた方法を含めるデータ群である。第2データ取得部64は第2処理条件の中から処理方法を決定する。第2処理条件のデータ群は複数の処理方法を含めてもよいし、予め決められた単一の処理方法を含めていてもよい。また、第2処理条件は、第1処理条件と同じものであってもよい。
The second processing condition is a data group including the method associated with the second extraction condition of the processing performed by the second
また、第2データ取得部64は、第1データの処理方法に応じて、第2データ取得部64で施す第2評価データの処理方法を決定してもよい。例えば、若年層のユーザから測定した実測値と脈波信号の相関関係を示した検量モデルを用いて、平均波形信号から第1データを取得する処理方法を施した場合、第2データ取得部64は、上記の処理方法に応じて、よりユーザに適した検量モデルを用いて、平均波形信号から第2データを取得してもよい。これによって、異なる処理の方法で取得した複数の生体情報を取得することが可能となり、より多角的で高精度な評価が可能となる。
Further, the second
また、第2データ取得部64は、第2処理条件を参照し、第1データの内容に対応する処理方法を決定してもよい。例えば、第1データによって、ユーザが低血糖の傾向があることが判明した場合、第2データ取得部64は、第2データ取得部64で施す第2評価データの処理方法として、低血糖のユーザの血糖値の実測値と脈波信号の相関関係を示した検量モデルを用いて、血糖値を取得する処理方法を決定してもよい。取得できる血糖値の値には大きな誤差がある場合が想定されるため、第1データの内容から低血糖帯、通常血糖帯、高血糖帯、超高血糖帯に分類し、その血糖帯に合わせた検量モデルを用いることで第2データとして取得する血糖値の精度を大幅に向上することが可能となる。また、第1データの内容からユーザが糖尿病かどうかを判断し、その結果に合わせて第2データ取得部64で施す第2評価データの処理方法を決定してもよい。これらによって、第1データの変化に追従して、第2評価データの処理を決定することで、より高精度な評価が可能になる。
Further, the second
処理方法は、例えば回帰分析部620で、平均波形信号から第1データを取得する処理で用いられる検量モデルの内容によっても分類することができる。例えば、若年層のユーザから測定した実測値と脈波信号の相関関係を示した検量モデルを用いて、平均波形信号から第1データを取得する処理方法と別の年代層のユーザから測定した実測値と脈波信号の相関関係を示した検量モデルを用いて、平均波形信号から第1データを取得する処理方法はそれぞれ異なる処理方法である。
The processing method can also be classified according to the content of the calibration model used in the processing of acquiring the first data from the average waveform signal, for example, in the
また、第2データ取得部64は、第1評価データの分類パターンを参照して、第1評価データを分類した結果を第1データとし、上述した第1データを参照して、第2データ取得部64で施す第2評価データの処理方法を決定してもよい。これによって、例えば分類しやすい加速度脈波を用いて、信号の特徴を分類したうえで、誤検出を抑制できる速度脈波を用いた高精度な評価が可能となる。これによって、第1データの変化に追従して、よりユーザに適した第2評価データの処理を決定することが可能となる。
Further, the second
分類パターンは、信号の波形の特徴に応じて、信号を2つ以上のグループに分類するための分類表のことである。 The classification pattern is a classification table for classifying signals into two or more groups according to the characteristics of the waveform of the signal.
図9は、加速度脈波の分類パターンの一例である。図9は、加速度脈波にはa~eの変曲点が存在し、加速度脈波における最大のピークをa点とし、a点から順に各変曲点をb点、c点、d点、e点とし、a点を1とし、最小値であるb点もしくはd点を0とした規格化を行い、各変曲点の値と、その差の大小関係により分類する方法で、脈波を7パターンに分類した図である。まず、変曲点の値がb<dの場合は、パターンAまたはBに分類する。b<dでさらにc≧0.5であればA、そうでなければBに分類する。次に変曲点の値がb≒dの場合、パターンCまたはDに分類する。b≒dでさらにc≒0の場合はパターンD、そうでなければパターンCに分類する。最後に、b>dの場合は、パターンE、F、Gのいずれかに分類できる。b>dでさらにb<cであればパターンEに、b≒cであればパターンF,b>cであればパターンGに分類する。 FIG. 9 is an example of the classification pattern of the acceleration pulse wave. In FIG. 9, inflection points a to e exist in the acceleration pulse wave, the maximum peak in the acceleration pulse wave is point a, and each inflection point is point b, point c, point d in order from point a. Normalization is performed with point e, point a as 1, and the minimum value b or d as 0, and the pulse wave is classified according to the value of each inflection point and the magnitude relationship of the difference. It is a figure classified into 7 patterns. First, when the value of the inflection point is b <d, it is classified into pattern A or B. If b <d and c ≧ 0.5, it is further classified as A, and if not, it is classified as B. Next, when the value of the inflection point is b≈d, it is classified into the pattern C or D. If b≈d and c≈0, it is classified into pattern D, and if not, it is classified into pattern C. Finally, if b> d, it can be classified into any of patterns E, F, and G. If b> d and b <c, it is classified into pattern E, if b≈c, it is classified into pattern F, and if b> c, it is classified into pattern G.
第1データ取得部62は、微分部611で脈波信号を微分した加速度脈波に基づいた第1評価データが入力された場合、第1評価データが、例えば図9のどのパターンに当てはまるかを判断し、分類パターンを決定する。例えば、入力された第1評価データの変曲点bが変曲点dより小さく、さらに変曲点c≧0.5であれば、パターンAを第1評価データの分類パターンとする。
When the first evaluation data based on the acceleration pulse wave obtained by differentiating the pulse wave signal by the
加速度脈波の分類ごとに、適した検量モデルが異なるので、第2データ取得部64は、第1評価データの分類パターンを参照して、第1評価データを分類した結果を第1データとし、上述した第1データを参照して、第2データ取得部64で施す第2評価データの処理の方法を決定することで、ユーザごとに適した処理方法を実現することができる。例えば、第1評価データがパターンAに分類された結果を第1データとするとき、第2データ取得部64は、回帰分析部640で、予め測定された脈波のパターンAの波形データを説明変数とし、生体情報の実測値を目的変数として、回帰分析により解析し、その解析結果に基づいて構築された検量モデルを用いて、入力された平均波形データから第2データを取得するように決定する。
Since the suitable calibration model differs depending on the classification of the acceleration pulse wave, the second
図10は、速度脈波の分類パターンの一例である。図10は、速度脈波における最大のピークを1とし、最小値を0とした規格化を行い、最大のピーク以外にピークが見られたパターンをグループ1とし、最大のピーク以外にピークが見られなかったパターンをグループ2とした分類パターンである。これを上述したように用いることで、第1データ取得部62は、微分部611で脈波信号を微分しない速度脈波に基づいた第1評価データが入力された場合でも分類パターンを決定することができる。
FIG. 10 is an example of a classification pattern of velocity pulse waves. In FIG. 10, normalization is performed with the maximum peak in the velocity pulse wave being 1 and the minimum value being 0, and the pattern in which the peak is seen other than the maximum peak is
次に、最適血糖値算出ステップS17において、第1データ取得部62で取得した第1データと、第2データ取得部64で取得した第2データを最適血糖値算出部65へ入力し、第1データと第2データから最適な生体情報を取得する。例えば、第1データとして取得した血糖値を第1血糖値とし、第2データとして取得した血糖値を第2血糖値とし、第1血糖値と第2血糖値から最適血糖値を算出する。最適値となる生体情報の算出方法については、例えば第1データと第2データにそれぞれの測定精度に伴って重み付けを行い、重み付けに基づき、複数の第1データおよび複数の第2データから最適値となる生体情報を算出してもよい。また、他の例として、他センサ等によって取得された血糖値を参照値とし、参照値と第1データ及び第2データのプロットグラフを作成し二つのプロットグラフのうち、エラーグリッド上で良好な値を示すデータを出力する方法と、第1データと第2データをそれぞれ複数取得し、バラつきの小さいデータを出力する方法と、予め設定された許容範囲に、第1データと第2データが範囲内に含まれるかを評価し、含まれるデータを出力する方法などがあげられる。
Next, in the optimum blood glucose level calculation step S17, the first data acquired by the first
1 生体情報演算システム
2 電子機器
3 公衆通信網
4 サーバ
5 センサ
16 出力I/F
20 ROM
21 RAM
22 CPU
23 操作部
24 記憶部
25 データ入出力部
26 内部バス
27 通信I/F
28 表示部
50 取得部
51 通信I/F
52 メモリ
53 命令部
54 内部バス
55 リストバンド
60 脈波信号取得部
61 第1抽出部
62 第1データ取得部
63 第2抽出部
64 第2データ取得部
65 最適血糖値算出部
610 フィルタ処理部
611 微分部
612 分割部
613 規格化部
614 ピーク位置算出部
615 ピーク間隔平均算出部
616 ピーク間隔プロット部
617 フーリエ変換部
618 最大周波数検出部
620 回帰分析部
621 脈拍処理部
622 ストレスプロット処理部
623 ストレスフーリエ処理部
624 呼吸数処理部
630 フィルタ処理部
631 微分部
632 分割部
633 規格化部
634 ピーク位置算出部
635 ピーク間隔平均算出部
636 ピーク間隔プロット部
637 フーリエ変換部
638 最大周波数検出部
640 回帰分析部
641 脈拍処理部
642 ストレスプロット処理部
643 ストレスフーリエ処理部
644 呼吸数処理部
S10 脈波信号取得ステップ
S11 第1抽出条件決定ステップ
S12 第1抽出ステップ
S13 第1データ取得ステップ
S14 第2抽出条件決定ステップ
S15 第2抽出ステップ
S16 第2データ取得ステップ
S17 最適血糖値算出ステップ
1 Biometric
20 ROM
21 RAM
22 CPU
23
28
52
Claims (3)
第1抽出条件を参照し、前記取得手段によって取得された脈波信号に基づく第1評価データを抽出する第1抽出手段と、
前記第1抽出条件に紐づく第1処理条件を参照し、前記第1抽出手段によって抽出された第1評価データに対する第1血糖値を、第1データとして取得する第1データ取得手段と、
第2抽出条件を参照し、前記取得手段によって取得された脈波信号に基づく第2評価データを抽出する第2抽出手段と、
前記第2抽出条件に紐づく第2処理条件を参照し、前記第2抽出手段によって抽出された第2評価データに対する、前記第1データと異なる種類の生体情報からなる第2データを取得する第2データ取得手段と、
を備え、
前記第1評価データと、前記第2評価データとは、同一の前記脈波信号に基づいて抽出されること
を特徴とする生体情報演算システム。 An acquisition means for acquiring a velocity pulse wave as a pulse wave signal,
With reference to the first extraction condition, the first extraction means for extracting the first evaluation data based on the pulse wave signal acquired by the acquisition means, and the first extraction means.
A first data acquisition means for acquiring the first blood glucose level for the first evaluation data extracted by the first extraction means as the first data with reference to the first processing condition associated with the first extraction condition.
With reference to the second extraction condition, the second extraction means for extracting the second evaluation data based on the pulse wave signal acquired by the acquisition means, and the second extraction means.
With reference to the second processing condition associated with the second extraction condition, the second data consisting of biological information of a type different from the first data is acquired with respect to the second evaluation data extracted by the second extraction means. 2 Data acquisition means and
Equipped with
A biological information calculation system characterized in that the first evaluation data and the second evaluation data are extracted based on the same pulse wave signal.
前記第2抽出手段は、前記取得手段によって取得された脈波信号を微分する前記第2抽出条件を参照すること
を特徴とする請求項1に記載の生体情報演算システム。 The first extraction means refers to the first extraction condition that does not differentiate the pulse wave signal acquired by the acquisition means.
The biometric information calculation system according to claim 1, wherein the second extraction means refers to the second extraction condition for differentiating the pulse wave signal acquired by the acquisition means.
前記第1回帰分析手段は、前記第1検量モデルを用いて、前記第1抽出手段によって抽出された第1評価データに対する第1血糖値を、前記第1データとして取得すること
を特徴とする請求項1又は2に記載の生体情報演算システム。 The first data acquisition means is the first extracted by the first extraction means using the first calibration model constructed based on the correlation between the first evaluation data acquired in advance and the measured value of the biological information. 1 Has a first regression analysis means for acquiring the first data for the evaluation data,
The first regression analysis means is characterized in that the first blood glucose level for the first evaluation data extracted by the first extraction means is acquired as the first data by using the first calibration model. Item 2. The biometric information calculation system according to Item 1 or 2.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2020206629A JP7083185B2 (en) | 2020-09-03 | 2020-12-14 | Biological information calculation system |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2020148048A JP6851665B1 (en) | 2020-09-03 | 2020-09-03 | Biological information calculation system |
JP2020206629A JP7083185B2 (en) | 2020-09-03 | 2020-12-14 | Biological information calculation system |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2020148048A Division JP6851665B1 (en) | 2020-09-03 | 2020-09-03 | Biological information calculation system |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2022042920A JP2022042920A (en) | 2022-03-15 |
JP7083185B2 true JP7083185B2 (en) | 2022-06-10 |
Family
ID=87654632
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2020206629A Active JP7083185B2 (en) | 2020-09-03 | 2020-12-14 | Biological information calculation system |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP7083185B2 (en) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010187993A (en) | 2009-02-19 | 2010-09-02 | Denso Corp | Blood pressure measuring device and program as well as recording medium |
JP2011217784A (en) | 2010-04-05 | 2011-11-04 | Seiko Epson Corp | Measuring device and program |
WO2016147795A1 (en) | 2015-03-13 | 2016-09-22 | 国立大学法人信州大学 | Non-invasive blood glucose level measurement method and non-invasive blood glucose level measurement device |
WO2016167202A1 (en) | 2015-04-17 | 2016-10-20 | 太陽誘電株式会社 | Vibration waveform sensor and waveform analysis device |
JP2017023472A (en) | 2015-07-23 | 2017-02-02 | 国立大学法人 和歌山大学 | Bio-information estimation device, bio-information estimation method, and computer program |
WO2017187710A1 (en) | 2016-04-28 | 2017-11-02 | 太陽誘電株式会社 | Vibration waveform sensor and pulse wave detector |
JP2019154774A (en) | 2018-03-13 | 2019-09-19 | 株式会社東芝 | Biometric measurement device and biometric measurement method |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10295657A (en) * | 1997-04-24 | 1998-11-10 | Matsushita Electric Ind Co Ltd | Blood pressure measuring device |
-
2020
- 2020-12-14 JP JP2020206629A patent/JP7083185B2/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010187993A (en) | 2009-02-19 | 2010-09-02 | Denso Corp | Blood pressure measuring device and program as well as recording medium |
JP2011217784A (en) | 2010-04-05 | 2011-11-04 | Seiko Epson Corp | Measuring device and program |
WO2016147795A1 (en) | 2015-03-13 | 2016-09-22 | 国立大学法人信州大学 | Non-invasive blood glucose level measurement method and non-invasive blood glucose level measurement device |
WO2016167202A1 (en) | 2015-04-17 | 2016-10-20 | 太陽誘電株式会社 | Vibration waveform sensor and waveform analysis device |
JP2017023472A (en) | 2015-07-23 | 2017-02-02 | 国立大学法人 和歌山大学 | Bio-information estimation device, bio-information estimation method, and computer program |
WO2017187710A1 (en) | 2016-04-28 | 2017-11-02 | 太陽誘電株式会社 | Vibration waveform sensor and pulse wave detector |
JP2019154774A (en) | 2018-03-13 | 2019-09-19 | 株式会社東芝 | Biometric measurement device and biometric measurement method |
Also Published As
Publication number | Publication date |
---|---|
JP2022042920A (en) | 2022-03-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108113683B (en) | Apparatus and method for estimating biological components | |
JP6544751B2 (en) | Non-invasive blood sugar level measuring method and non-invasive blood sugar level measuring device | |
JP6525138B2 (en) | Blood pressure measuring device | |
US20060224074A1 (en) | Heartbeat measuring apparatus | |
WO2005110211A1 (en) | Bioinformation detector | |
JP6899609B1 (en) | Biometric information calculation system and server | |
JP7083185B2 (en) | Biological information calculation system | |
US20190298190A1 (en) | Pulse detection, measurement and analysis based health management system, method and apparatus | |
JP6851665B1 (en) | Biological information calculation system | |
JP6851664B1 (en) | Biological information calculation system | |
KR102547612B1 (en) | A method for generating heart rate variability inforamtion related to an external object using a plurality filters and an electronic device thereof | |
TWI640297B (en) | Non-invasive blood glucose measuring device, method, and system with identification function | |
WO2022050333A1 (en) | Biometric information computing system, server, and data structure | |
WO2022050334A1 (en) | Biological information calculation system | |
JP6845520B1 (en) | Biological information calculation system | |
JP7083194B1 (en) | Biological information calculation system | |
CN117045217B (en) | Cuff-free blood pressure measurement method and related equipment thereof | |
JP7455191B2 (en) | Biological abnormality detection device, biological abnormality detection method, and program | |
JP7083195B1 (en) | Biological information calculation system | |
JP2010200911A (en) | Apparatus for measuring blood component | |
Almeida et al. | Validation of a waveform delineator device for cardiac studies: Repeatability and data mining analysis | |
JP2023153709A (en) | Biological information operation system, server, and biological information operation method | |
CN116530942A (en) | Monitoring device based on spectrum sensor | |
CN116157062A (en) | Biological information calculation system | |
KR20230062852A (en) | Biometric information operation system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20201214 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20211214 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20220209 |
|
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: 20220426 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20220524 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 7083185 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |