(57)【要約】【課題】血流速度の計測において体動に起因するノイズを低減することが可能な生体情報処理装置を提供すること。【解決手段】拍動間隔に基づいて第１の周波数帯域を選択し、血流速度信号の中から前記第１の周波数帯域の信号を通過させて、拍動成分信号を抽出する第１のフィルタ部を備える、生体情報処理装置を提供する。【選択図】図３
2020006014.tif 000002
【特許請求の範囲】
【請求項１】
拍動間隔に基づいて第１の周波数帯域を選択し、血流速度信号の中から前記第１の周波数帯域の信号を通過させて、拍動成分信号を抽出する第１のフィルタ部を備える、生体情報処理装置。
【請求項２】
前記拍動成分信号のピーク位置を検出し、ピーク時間と該ピーク時間における拍動成分強度とを出力するピーク検出部を備え、
前記ピーク検出部は、
体動発生時以外においては、前記拍動間隔に基づいて前記拍動成分信号のピーク位置を検出して、拍動成分強度を履歴情報として保持し、
体動発生時においては、前記拍動間隔及び前記履歴情報に基づいて前記拍動成分信号のピーク位置を検出した後、前記履歴情報に基づいて前記ピーク位置の検出が誤りであるか否かを判定し、前記ピーク位置の検出が誤りである場合は、前記拍動間隔を用いてピーク位置を補正し、前記履歴情報を用いて拍動成分強度を補正する、請求項１に記載の生体情報処理装置。
【請求項３】
前記第１のフィルタ部は、複数のバンドパスフィルタで構成される第１のフィルタバンクを備え、前記拍動間隔に基づいて前記第１のフィルタバンクから選択されたバンドパスフィルタを用いて前記拍動成分信号を抽出する、請求項１に記載の生体情報処理装置。
【請求項４】
前記ピーク検出部は、体動発生時以外において、前記拍動間隔に基づいて第２の周波数帯域を選択し、前記拍動成分信号の中から前記第２の周波数帯域の信号を通過させて、前記拍動成分信号のピーク位置を検出する第２のフィルタ部を備える、請求項２に記載の生体情報処理装置。
【請求項５】
前記第２のフィルタ部は、複数のバンドパスフィルタで構成される第２のフィルタバンクを備え、体動発生時以外において、前記拍動間隔に基づいて前記第２のフィルタバンクから選択されたバンドパスフィルタを用いて前記拍動成分信号のピーク位置を検出する、請求項４に記載の生体情報処理装置。
【請求項６】
前記ピーク検出部は、体動発生時において、前記ピーク位置の検出が誤りである場合は、前記拍動間隔により前記ピーク位置の確率分布の尤度を計算し、前記尤度が最大となるピーク位置を補正後のピーク位置とする、請求項２に記載の生体情報処理装置。
【請求項７】
拍動間隔に基づいて第１の周波数帯域を選択し、血流速度信号の中から前記第１の周波数帯域の信号を通過させて、拍動成分信号を抽出する第１のフィルタ工程を含む、生体情報処理方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本技術は、生体情報処理装置及び生体情報処理方法に関する。より詳しくは、血流速度を計測するための生体情報処理装置及び生体情報処理方法に関する。
【背景技術】
【０００２】
近年、ヘルスケアやウェルネスのブームにともない、日常生活において人の生理反応を計測し、健康状態や心理状態をセンシングする技術が注目されている。人の心理状態が変化すると、脳から自律神経系を介して信号が伝達され、呼吸、皮膚温度、発汗、心臓、血管活動などの各機能に変化を引き起こす。心理状態の変化に伴って心室の拡張と収縮のタイミングや心拍出量が変化することから、指先や手首、耳朶、額などにおいて心室の収縮期と拡張期に応じた血流速度の変化を検出することにより、心理状態をセンシングできることが知られている。このように、血流速度は人の心理状態を表す１つの指標として注目されている。
【０００３】
人の生理反応を計測する際に、計測している部位が動くと、体動に起因するノイズ（体動ノイズ）が観測信号中に発生してＳ／Ｎ比が低下することから、計測の精度を上げるべく、体動ノイズを低減する技術が研究されている。
【０００４】
例えば、特許文献１には、正確な脈拍数を得ることを課題とする技術が記載されている。特許文献２には、精度の高い心拍測定を可能とすることを課題とする技術が記載されている。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開平１１−２７６４４８号公報
【特許文献２】国際公開第２０１８／０５５９６９号
【発明の概要】
【発明が解決しようとする課題】
【０００６】
上記特許文献１及び２の技術は、脈拍数や心拍数の計測に関する技術である。しかしながら、心理状態を表す指標として近年注目されている血流速度については、体動ノイズを低減する技術が確立しておらず、計測精度の向上が求められている。
【０００７】
そこで、本技術は、血流速度の計測において体動に起因するノイズを低減することが可能な生体情報処理装置を提供することを主目的とする。
【課題を解決するための手段】
【０００８】
すなわち、本技術は、
拍動間隔に基づいて第１の周波数帯域を選択し、血流速度信号の中から前記第１の周波数帯域の信号を通過させて、拍動成分信号を抽出する第１のフィルタ部を備える、生体情報処理装置を提供する。
前記生体情報処理装置は、
前記拍動成分信号のピーク位置を検出し、ピーク時間と該ピーク時間における拍動成分強度とを出力するピーク検出部を備えてもよく、
前記ピーク検出部は、
体動発生時以外においては、前記拍動間隔に基づいて前記拍動成分信号のピーク位置を検出して、拍動成分強度を履歴情報として保持し、
体動発生時においては、前記拍動間隔及び前記履歴情報に基づいて前記拍動成分信号のピーク位置を検出した後、前記履歴情報に基づいて前記ピーク位置の検出が誤りであるか否かを判定し、前記ピーク位置の検出が誤りである場合は、前記拍動間隔を用いてピーク位置を補正し、前記履歴情報を用いて拍動成分強度を補正してもよい。
前記第１のフィルタ部は、複数のバンドパスフィルタで構成される第１のフィルタバンクを備えてもよく、前記拍動間隔に基づいて前記第１のフィルタバンクから選択されたバンドパスフィルタを用いて前記拍動成分信号を抽出してもよい。
前記ピーク検出部は、体動発生時以外において、前記拍動間隔に基づいて第２の周波数帯域を選択し、前記拍動成分信号の中から前記第２の周波数帯域の信号を通過させて、前記拍動成分信号のピーク位置を検出する第２のフィルタ部を備えてもよい。
前記第２のフィルタ部は、複数のバンドパスフィルタで構成される第２のフィルタバンクを備えてもよく、体動発生時以外において、前記拍動間隔に基づいて前記第２のフィルタバンクから選択されたバンドパスフィルタを用いて前記拍動成分信号のピーク位置を検出してもよい。
前記ピーク検出部は、体動発生時において、前記ピーク位置の検出が誤りである場合は、前記拍動間隔により前記ピーク位置の確率分布の尤度を計算し、前記尤度が最大となるピーク位置を補正後のピーク位置としてもよい。
また、本技術は、拍動間隔に基づいて第１の周波数帯域を選択し、血流速度信号の中から前記第１の周波数帯域の信号を通過させて、拍動成分信号を抽出する第１のフィルタ工程を含む、生体情報処理方法を提供する。
【発明の効果】
【０００９】
本技術によれば、血流速度の計測において体動に起因するノイズを低減することが可能な生体情報処理装置を提供することができる。なお、本技術の効果は、ここに記載された効果に必ずしも限定されるものではなく、本明細書中に記載されたいずれかの効果であってもよい。
【図面の簡単な説明】
【００１０】
【図１】血流速度信号の例を示す模式図である。
【図２】第１実施形態に係る生体情報処理装置の構成例を示す図である。
【図３】第１実施形態に係る生体情報処理装置の機能構成の一例を示すブロック図である。
【図４】第１のフィルタバンクに記憶されたデータの一例を示す表である。
【図５】バンドパスフィルタの一例を示すグラフである。
【図６】第１のフィルタ工程の一例を示すフローチャートである。
【図７】ピーク検出工程の一例を示すフローチャートである。
【図８】体動発生時におけるピーク位置の補正の一例を説明するための図である。
【図９】第２のフィルタバンクに記憶されたデータの一例を示す表である。
【図１０】第２のフィルタ工程の一例を示すフローチャートである。
【図１１】第２のフィルタ工程の一例を説明するための図である。
【図１２】第１実施形態の第１変形例に係る生体情報処理装置の機能構成の一例を示すブロック図である。
【図１３】第１実施形態の第２変形例に係る生体情報処理装置の機能構成の一例を示すブロック図である。
【発明を実施するための形態】
【００１１】
以下、本技術を実施するための好適な形態について図面を参照しながら説明する。なお、以下に説明する実施形態は、本技術の代表的な実施形態を示したものであり、これにより本技術の範囲が狭く解釈されることはない。なお、説明は以下の順序で行う。
１．安静時及び体動発生時の血流速度信号についての説明
２．第１実施形態
（１）生体情報処理装置１の構成
（２）生体情報処理装置１の動作
３．第１実施形態の第１変形例
４．第１実施形態の第２変形例
【００１２】
＜１．安静時及び体動発生時の血流速度信号についての説明＞
図１は、血流速度信号の例を示す模式図である。図１Ａは安静時における血流速度信号であり、図１Ｂは体動発生時における血流速度信号である。図１Ｂ中のピークＰ及びピークＱは、体動に起因する偽のピークを示している。
【００１３】
例えば、人の心理状態の変化を心室の収縮期と拡張期に応じた血流速度の変化によりセンシングする場合、収縮期と拡張期の血流を測定するために血流速度信号のピーク位置をより正確に検出する必要がある。図１Ａに示すように測定対象者の体動が発生していない時、即ち測定対象者が安静にしている時には、血流速度信号中のピークは心臓の拍動に伴って観察されるため、ピーク位置の検出に誤りが生じる可能性は低い。一方で、図１Ｂに示すように、測定対象者の体動に起因して体動ノイズが発生し、拍動に伴わない偽のピーク（ピークＰ及びピークＱ）が観察される時には、ピーク位置の誤検出が起きやすくなる。このため、ピーク位置の誤検出を回避し、血流速度の変化をより正確に検出するためには、体動ノイズを低減する必要がある。
【００１４】
以上のとおり、血流速度の計測において血流速度信号のピーク位置をより正確に検出すべく体動ノイズを低減する技術が求められている。本技術は、血流速度の計測において体動ノイズを低減する技術を提供する。
【００１５】
＜２．第１実施形態＞
（１）生体情報処理装置１の構成
図２を参照して、第１実施形態に係る生体情報処理装置１の構成について説明する。
【００１６】
図２は、第１実施形態に係る生体情報処理装置１の構成例を示す図である。生体情報処理装置１は、血流センサ部１０と制御部２０とを備える。
【００１７】
生体情報処理装置１の実施形態としては、例えば、ヘッドバンド型、ネックバンド型、ベルト型などの種々のウェアラブル装置；スマートフォンやタブレット端末などの任意の携帯情報端末（PDA：Personal Digital Assistant）；医療機器、ゲーム機器、家電機器などの任意の電子機器；などが挙げられる。生体情報処理装置１は、物理的に１つの装置などにより構成されてもよく、物理的に複数の装置などにより構成されてもよい。例えば、血流センサ１０の機能を有するウェアラブル装置と、当該ウェアラブル装置と通信可能に接続し制御部２０の機能を有する携帯情報端末と、により生体情報処理装置１が構成されてもよい。
【００１８】
血流センサ部１０は、血流を測定することができればその構成は特に限定されない。血流センサ部１０として、例えば、レーザードップラー血流計（Laser Doppler Flowmetry：ＬＤＦ）を採用することができる。レーザードップラー血流計は、人の皮膚表面にレーザー光を照射して毛細血管内血流を非侵襲に且つ連続的に計測することが可能である上に、小型である。このため、レーザードップラー血流計は、測定対象者が日常生活において簡便に装着できるウェアラブル装置に適している。
【００１９】
制御部２０は、例えばＣＰＵやメモリ（ＲＡＭ、ＲＯＭ）などのコンピュータに必要なハードウェアから構成される。ＣＰＵがＲＯＭなどに記録されているプログラムやデータをＲＡＭ上に呼び出して処理を実行することで、生体情報処理装置１の全体の制御や機能が実現される。制御部２０として、例えばＦＰＧＡ（Field Programmable Gate Array）などのＰＬＤ（Programmable Logic Device）、その他ＡＳＩＣ（Application Specific Integrated Circuit）などのデバイスが用いられてもよい。
【００２０】
図３は、第１実施形態に係る生体情報処理装置１の機能構成の一例を示すブロック図である。心拍変動検出部２１、血流速度計算部２２、第１のフィルタ部２３、ピーク検出部２４、血流速度ベースライン復元部２５及び包絡線検出部２６は、制御部２０により実現される。
【００２１】
（２）生体情報処理装置１の動作
次に、図３を参照して、第１実施形態に係る生体情報処理装置１の動作について説明する。
【００２２】
心拍変動検出部２１は、測定対象者の心拍変動を検出する心拍変動検出工程を実行し、拍動間隔を出力する。図３において、「ＰＰＧ」は光電容積脈波方式（Photoplethysmography：ＰＰＧ方式）の脈波センサを示しており、「ＥＣＧ」は心電図（Electrocardiogram）を意味している。本実施形態において、心拍変動検出部２１は、ＰＰＧ方式の脈波センサ又は心電図から得られる情報に基づいて心拍変動を検出し、拍動間隔を算出する。
【００２３】
血流速度計算部２２は、測定対象者の血流速度を計算する血流速度計算工程を実行し、血流速度信号を出力する。図２において、「ＬＤＦ」はレーザードップラー血流計を意味しており、血流センサ部１０（図２）がレーザードップラー血流計である場合を示している。人の皮膚表面にレーザー光を照射し、反射光として観測されるビート信号から、下記数式（Ｉ）で血流速度（Perfusion）を算出することができる。
【００２４】
【数１】
2020006014.tif 000003
【００２５】
上記数式（Ｉ）中、ωは周波数、Ｐ（ω）は周波数スペクトラム密度関数を表す。なお、血流速度を算出する手法は上述した手法に限定されるものではない。
【００２６】
第１のフィルタ部２３は、血流速度計算部２２から出力された血流速度信号の中から拍動成分信号を抽出する第１のフィルタ工程を実行する。拍動成分信号は、血流速度信号に含まれる拍動に伴う変動成分（Alternating Current成分：ＡＣ成分）である。
【００２７】
体動発生時以外（安静時）において、血流速度信号中のピークは拍動に伴って検出され、ピークとピークの間の時間は拍動間隔分だけ離れる。これに対して、体動発生時には、血流速度信号中に拍動間隔との相関性が低い偽のピークが発生する。第１のフィルタ部２３は、拍動間隔の情報を用いて血流速度信号の中から拍動成分信号を抽出する第１のフィルタ工程を実行することにより、体動ノイズである拍動に伴わない偽のピークを除去する。つまり、本実施形態の生体情報処理装置１は、第１のフィルタ工程を実行することにより、血流速度信号中の体動ノイズを低減することができる。
【００２８】
上記第１のフィルタ工程では、第１のフィルタ部２３が、心拍変動検出部２１から出力された拍動間隔に基づいて第１の周波数帯域を選択し、上記血流速度信号の中から第１の周波数帯域の信号を通過させるフィルタ処理を行う。このフィルタ処理により、第１のフィルタ部２３は拍動成分信号を抽出する。
【００２９】
第１のフィルタ部２３における上記フィルタ処理として、バンドパスフィルタ処理を採用することができる。バンドパスフィルタ処理を実行するため、第１のフィルタ部２３は、複数のバンドパスフィルタで構成される第１のフィルタバンク２３ａを備えうる。第１のフィルタ部２３は、心拍変動検出部２１から出力された拍動間隔に基づいて、第１のフィルタバンク２３ａからバンドパスフィルタを選択し、当該バンドパスフィルタを用いることにより上記拍動成分信号を抽出することができる。第１のフィルタバンク２３ａは、複数のバンドパスフィルタを、各フィルタに対応するフィルタ係数として保持している。
【００３０】
図４は、第１のフィルタバンク２３ａに記憶されたデータの一例を示す表である。第１のフィルタバンク２３ａには、各バンドパスフィルタのフィルタＩＤ（１、２、・・・Ｎ）、拍動間隔の範囲、及びフィルタ係数がそれぞれ対応付けられて記憶されている。図４に示す例では、４０ｂｐｍから２００ｂｐｍまでの拍動間隔に対応するＮ種類のフィルタ係数が示されている。
【００３１】
例えば、フィルタＩＤが２である項目には、拍動間隔の範囲として、拍動間隔（ｈｒ＿ｔｒｅｎｄ）が５５ｂｐｍ以上６５ｂｐｍ以未満の場合についての条件式（５５ｂｐｍ≦ｈｒ＿ｔｒｅｎｄ＜６５ｂｐｍ）が記憶されている。また、フィルタ係数として、ＢＰ１ｃｏｅｆ＿ｖｅｃｔｏｒ２のデータが記憶されている。
【００３２】
図５は、バンドパスフィルタの一例を示すグラフである。図５Ａは、フィルタＩＤが２である場合のフィルタ係数（ＢＰ１ｃｏｅｆ＿ｖｅｃｔｏｒ２）のグラフである。図５Ａの横軸はタップ数であり、縦軸はフィルタ係数である。図５Ａに示したフィルタ係数を用いてバンドパスフィルタが構成される。
【００３３】
図５Ｂは、フィルタＩＤが２のフィルタ係数を使ったバンドパスフィルタの周波数特性である。図５Ｂの横軸は周波数であり、縦軸はバンドパスフィルタを通過する信号の強度（透過率）を表している。図５ＢにおいてＢＰ１ｃｏｅｆ＿ｖｅｃｔｏｒ２の周波数特性を確認すると、拍動間隔が６０ｂｐｍの周辺の帯域のみを通過させるフィルタ特性を有することが理解できる。なお、バンドパスフィルタの周波数特性は、血流速度信号から拍動成分信号を抽出することを目的に適宜設定されればよい。
【００３４】
このように、第１のフィルタバンク２３ａに記憶された複数のフィルタ係数から、各拍動間隔に対応するバンドパスフィルタを選択することが可能である。これにより、拍動間隔の範囲ごとに適正なフィルタ処理を行うことが可能であり、血流速度信号から拍動成分信号を抽出することが可能となる。
【００３５】
図６は、第１のフィルタ部２３が実行する第１のフィルタ工程の一例を示すフローチャートである。まず、第１のフィルタ部２３に拍動間隔が入力される（ステップＳ１０１）。第１のフィルタ部２３は、第１のフィルタバンク２３ａから拍動間隔の該当範囲を選択し（ステップＳ１０２）、対応するフィルタＩＤを決定する（ステップＳ１０３）。図５に示す例で説明すると、ステップＳ１０１において入力された拍動間隔が６０ｂｐｍである場合には、５５ｂｐｍ以上６５ｂｐｍ未満の範囲に含まれると判定され、フィルタＩＤとしてＩＤ＝２が決定される。
【００３６】
第１のフィルタ部２３は、フィルタＩＤからフィルタ係数を決定する（ステップＳ１０４）。図５に示す例で説明すると、ＩＤ＝２と決定された場合のフィルタ係数はＢＰ１ｃｏｅｆ＿ｖｅｃｔｏｒ２となる。このように、入力された拍動間隔を含む周波数帯域を通過させるバンドパスフィルタが決定される。決定されたフィルタ係数の通過帯域は、本実施形態において、上述した第１の周波数帯域に相当する。
【００３７】
第１のフィルタ部２３は、決定されたフィルタ係数からバンドパスフィルタを選択し、当該バンドパスフィルタを用いてバンドパスフィルタ処理を行う（ステップＳ１０５）。バンドパスフィルタ処理によって、拍動成分以外の帯域が除去されて、血流速度信号の中から拍動成分信号が抽出される。このようにして、第１のフィルタ部２３は、血流速度信号から体動に起因するノイズ（体動ノイズ）を除去することができる。
【００３８】
図３に戻り、第１のフィルタ部２３について更に説明する。第１のフィルタ部２３は、第１のフィルタ工程において、拍動成分信号を抽出する処理で除去された拍動成分以外の帯域のうち拍動成分信号よりも低域側の信号を血流速度ベースライン成分信号として保持し、出力する。血流速度ベースライン成分信号は、血流速度信号に含まれる拍動に伴わない基線成分（Direct Current成分:ＤＣ成分）である。血流速度ベースライン成分信号は、後述する血流速度ベースライン復元部２５に入力される。
【００３９】
ピーク検出部２４は、第１のフィルタ部２３から出力された拍動成分信号のピーク位置を検出するピーク検出工程を実行し、ピーク時間、及び当該ピーク時間における拍動成分強度を出力する。ピーク検出部２４は、体動発生時か否かに応じて、拍動成分強度を履歴情報として保持するか又は拍動成分強度を補正する。これにより、体動ノイズをより低減させて、ピーク位置検出の精度を高めることができる。
【００４０】
図７は、ピーク検出部２４が実行するピーク検出工程の一例を示すフローチャートである。ピーク検出部２４における処理に関して、まず、図７を参照して全体の流れを説明する。
【００４１】
ピーク検出部２４は、体動発生時であるか体動発生時以外であるかを判断する（ステップＳ２０１）。体動発生時であるか否かは、体動を検知する公知の技術を採用又は応用することにより判断すればよい。例えば、加速度センサやジャイロセンサの出力結果に基づいて判断してもよいし、血流センサ部１０から出力される血流信号の波形に乱れが生じたか否かで判断してもよい。また、血流センサ部１０から出力される血流信号のデータを蓄積して人工知能（Artificial Intelligence：ＡＩ）に学習させ、人工知能により体動発生時か否かを決定してもよい。
【００４２】
体動発生時以外であると判断した場合（ステップＳ２０１：ＮＯ）、拍動間隔に基づいて拍動成分信号のピーク位置を検出し（ステップＳ２０２）、拍動成分強度を履歴情報として保持する（ステップＳ２０３）。ステップＳ２０２において、拍動間隔は心拍変動検出部２１から出力された拍動間隔であり、拍動成分信号は第１のフィルタ部２３から出力された拍動成分信号である。
【００４３】
体動発生時以外の場合は、体動ノイズが少なく拍動成分信号の信号品質が高いため、拍動に伴うピーク成分が十分に保持されている。このため、ステップＳ２０２及びＳ２０３においては、ピーク位置の検出を拍動間隔に基づいて行って体動の影響が少なく信頼度の高い拍動成分強度を履歴情報として保持する。
【００４４】
一方、体動発生時であると判断した場合（ステップＳ２０１：ＹＥＳ）、拍動間隔及び履歴情報に基づいて拍動成分信号のピーク位置を検出する（ステップＳ２０４）。ステップＳ２０４において、拍動間隔は心拍変動検出部２１から出力された拍動間隔であり、履歴情報は上記ステップＳ２０３で保持した拍動成分強度を含む。当該履歴情報は、ステップＳ２０３で保持した拍動成分強度以外の拍動成分強度情報を含んでもよい。例えば、ピーク検出部２４は、ピーク検出工程を実行する前に取得した拍動成分強度を、履歴情報として予め保持してもよい。例えば、ピーク検出部２４は、体動の影響を受けないように測定対象者が安静な状態で計測した血流速度の結果に基づく拍動成分強度を、履歴情報として予め保持してもよい。
【００４５】
ピーク検出部２４は、履歴情報に基づいてピーク位置の検出が誤りであるか否かを判定する（ステップＳ２０５）。例えば、検出されたピークの拍動成分強度と履歴情報に含まれる拍動成分強度との差分を算出し、当該差分が閾値以上であれば検出されたピーク位置は誤りであると判断することができる。
【００４６】
ピーク位置の検出が誤りではない場合は（ステップＳ２０５：ＮＯ）、処理を終了する。ピーク位置の検出が誤りである場合は（ステップＳ２０５：ＹＥＳ）、拍動間隔を用いてピーク位置を補正し（ステップＳ２０６）、履歴情報を用いて拍動成分強度を補正する（ステップＳ２０７）。ステップＳ２０６の拍動間隔とステップＳ２０７の履歴情報は、上記ステップＳ２０４で説明した拍動間隔及び履歴情報と同様である。なお、ステップＳ２０６とＳ２０７の順序は、どちらが先でもよく、同時でもよい。
【００４７】
体動発生時の場合、第１のフィルタ部２３によって体動ノイズを低減してもなお、拍動成分信号の信号品質を更に高めたい場合がある。ピーク検出部２４は、信頼性の高い拍動成分強度が履歴として保持されている履歴情報を拍動間隔と組み合わせて用いることで、ピーク位置検出の精度を高めている（ステップＳ２０４）。また、ピーク検出部２４は、誤検出発見時にピーク位置と拍動成分強度を補正することで、信号品質を向上させている（ステップＳ２０６、２０７）。このようにして、ピーク検出部２４は体動ノイズの更なる低減を実現する。
【００４８】
拍動成分強度を補正する処理（ステップＳ２０７）では、例えば、履歴情報として保持されている拍動成分強度を、補正後の拍動成分強度として出力すればよい。
【００４９】
次に、ピーク検出部２４における処理に関して、図８を参照して、体動発生時に拍動間隔を用いてピーク位置を補正するステップ（図７中のステップＳ２０６）の一例を説明する。
【００５０】
図８は、体動発生時におけるピーク位置の補正の一例を説明するための図である。図８中、血流速度信号上に付した×印は体動ノイズに起因する偽のピークであり、黒色の丸印は真のピークである。
【００５１】
ピーク検出部２４は、体動発生時において、ピーク位置の検出が誤りであると判定した場合は、拍動間隔によりピーク位置の確率分布の尤度を計算し、当該尤度が最大となるピーク位置を補正後のピーク位置とすることができる。このようにして、ピーク検出部２４は、ピーク位置を補正することができる。
【００５２】
図８に示す例では、検出された真のピークＲのピーク位置を起点として、拍動間隔の分だけ離れた次のピークＳが真のピークである確率が最も高くなる。偽のピークＺを真のピークとして誤って検出した場合、拍動間隔によってピーク位置の確率分布の尤度を計算し、尤度が最大となるピークＳのピーク位置を用いて、偽のピークＺのピーク位置を補正する。具体的には、ピークＲの次のピークの位置を、ピークＺの位置ではなくピークＳの位置に補正する。
【００５３】
次に、ピーク検出部２４における処理に関して、図３及び図９を参照して、体動発生時以外においてピーク位置を検出するステップ（図７中のステップＳ２０２）について、更に説明する。
【００５４】
体動発生時以外では、体動ノイズが少なく拍動成分信号の信号品質は高い。しかしながら、拍動成分信号中に偽のピークが発生する可能性を完全に排除することができない場合、偽のピークが発生してピークの誤検出につながりうる。このため、ピーク検出部２４は、以下にて説明する第２のフィルタ工程を実行することにより拍動間隔に応じたピーク位置の検出を行い、更なる体動ノイズ低減を実現する。
【００５５】
図３に示すように、ピーク検出部２４は、第２のフィルタ部２４１を備えることがきる。第２のフィルタ部２４１は、体動発生時以外において、拍動間隔に基づいて第２の周波数帯域を選択し、拍動成分信号の中から第２の周波数帯域の信号を通過させるフィルタ処理を行って拍動成分信号のピーク位置を検出する第２のフィルタ工程を実行する。この拍動成分信号は、第１のフィルタ部２３から出力された拍動成分信号である。
【００５６】
第２のフィルタ部２４１におけるフィルタ処理として、バンドパスフィルタ処理を採用することができる。バンドパスフィルタ処理を実行するため、第２のフィルタ部２４１は、複数のバンドパスフィルタで構成される第２のフィルタバンク２４１ａを備えうる。第２のフィルタ部２４１は、体動発生時以外において、心拍変動検出部２１から出力された拍動間隔に基づいて、第２のフィルタバンク２４１ａからバンドパスフィルタを選択し、当該バンドパスフィルタを用いることにより上記拍動成分信号のピーク位置を検出することができる。第２のフィルタバンク２４１ａは、複数のバンドパスフィルタを、各フィルタに対応するフィルタ係数として保持している。
【００５７】
図９は、第２のフィルタバンク２４１ａに記憶されたデータの一例を示す表である。第２のフィルタバンク２４１ａには、各バンドパスフィルタのフィルタＩＤ（１、２、・・・Ｎ）、拍動間隔の範囲、及びフィルタ係数がそれぞれ対応付けられて記憶されている。図９に示す例では、４０ｂｐｍから２００ｂｐｍまでの拍動間隔に対応するＮ種類のフィルタ係数が示されている。
【００５８】
例えば、フィルタＩＤが２である項目には、拍動間隔の範囲として、拍動間隔（ｈｒ＿ｔｒｅｎｄ）が５５ｂｐｍ以上６５ｂｐｍ以未満の場合についての条件式（５５ｂｐｍ≦ｈｒ＿ｔｒｅｎｄ＜６５ｂｐｍ）が記憶されている。また、フィルタ係数として、ＢＰ２ｃｏｅｆ＿ｖｅｃｔｏｒ２のデータが記憶されている。
【００５９】
図１０は、第２のフィルタ部２４１が実行する第２のフィルタ工程の一例を示すフローチャートである。まず、第２のフィルタ部２４１に拍動間隔が入力される（ステップＳ３０１）。第２のフィルタ部２４１は、第２のフィルタバンク２４１ａから拍動間隔の該当範囲を選択し（ステップＳ３０２）、対応するフィルタＩＤを決定する（ステップＳ３０３）。図９に示す例で説明すると、ステップＳ３０１において入力された拍動間隔が６０ｂｐｍである場合には、５５ｂｐｍ以上６５ｂｐｍ未満の範囲に含まれると判定され、フィルタＩＤとしてＩＤ＝２が決定される。
【００６０】
第２のフィルタ部２４１は、フィルタＩＤからフィルタ係数を決定する（ステップＳ３０４）。図９に示す例で説明すると、ＩＤ＝２と決定された場合のフィルタ係数はＢＰ２ｃｏｅｆ＿ｖｅｃｔｏｒ２となる。このように、入力された拍動間隔を含む周波数帯域を通過させるバンドパスフィルタが決定される。決定されたフィルタ係数の通過帯域は、本実施形態において、上述した第２の周波数帯域に相当する。
【００６１】
第２のフィルタ部２４１は、決定されたフィルタ係数からバンドパスフィルタを選択し、当該バンドパスフィルタを用いてバンドパスフィルタ処理を行う（ステップＳ３０５）。バンドパスフィルタ処理によって、拍動成分信号のピーク位置を検出する。
【００６２】
図１１は、第２のフィルタ工程の一例を説明するための図である。図１１Ａ中、拍動成分信号上に付した白色の丸印は現在のピーク位置、黒色の丸印は真のピーク、×印は偽のピークである。例えば、図１１Ｂに示すように第２のフィルタバンク２４１ａ（図９）からフィルタ係数をＢＰ２ｃｏｅｆ＿ｖｅｃｔｏｒ２と決定し、図１１Ａに示す拍動成分信号に重畳する。その結果、図１１Ｃに示すように拍動間隔に応じた真のピークの強度のみを強調し、偽のピークの強度を低減することができる。このようにして、第２のフィルタ部２４１は、体動ノイズを低減して拍動成分信号のピーク位置を精度よく検出することができる。
【００６３】
図３に戻り、次に、血流速度ベースライン復元部２５について説明する。血流速度ベースライン復元部２５は、第１のフィルタ部２３から出力された血流速度ベースライン成分信号と、ピーク検出部２４から出力されたピーク時間及び拍動成分強度と、を加算する、血流速度ベースライン復元工程を実行する。これにより、血流速度ベースライン復元部２５は血流強度を復元し、ピーク時間及び血流強度を出力する。
【００６４】
第１のフィルタ部２３から出力された任意のピーク時間（ｔｐｋｓ）における血流速度ベースライン成分を「ｔｐｋｓ，ＤＣｐｋｓ」と表記し、ピーク検出部２４から出力されたピーク時間（ｔｐｋｓ）における拍動成分強度を「ｔｐｋｓ，ＡＣｐｋｓ」と表記する。この場合、血流速度ベースライン復元部２５から出力されるピーク時間（ｔｐｋｓ）における血流強度は「ｔｐｋｓ，ＡＣｐｋｓ＋ＤＣｐｋｓ」である。
【００６５】
包絡線検出部２６は、血流速度ベースライン復元部２５から出力されたピーク時間及び血流強度（ｔｐｋｓ，ＡＣｐｋｓ＋ＤＣｐｋｓ）のデータ列に対し、各ピーク間の一次元データ内挿により血流速度信号の包絡線を計算する包絡線検出工程を実行する。一次元データ内挿では、公知であるＣｕｂｉｃ補間などのピーク間を滑らかに補間するアルゴリズムを使用し、包絡線の曲線を計算する。
【００６６】
以上詳述した動作により、本実施形態の生体情報処理装置１は、体動ノイズが低減された血流速度信号を出力する。
【００６７】
本技術の生体情報処理装置は、血流速度の計測における体動ノイズを低減することが可能であるため、血流速度を計測するあらゆる場面で効果を発揮するが、測定対象者が日常生活環境下で装着できるウェアラブル装置に適している。
【００６８】
例えば、ウォーキングやジョギング、ランニングなどの運動時における主な体動ノイズの要因は腕の周期的な動きであるため、腕の動きを加速度センサなどで計測して参照信号とすることで、血流速度の情報に混入する体動ノイズを低減することは可能である。しかしながら、日常生活において血流速度を常時計測する場合、体動ノイズは指や手首などの非周期的な動きに起因するものが多く、腕の周期的な動きに起因するものは少ない。このため、上述のように加速度センサなどを用いても低減できる体動ノイズの範囲は狭く、日常生活の環境下で血流速度の変化を常時、精度よく検出することは困難であった。
【００６９】
これに対して、本技術の生体情報処理装置は、周期的か否かといった体動ノイズの性質にかかわらず、日常生活で発生する体動ノイズを効果的に低減することが可能である。また、本技術の生体情報処理装置は、ＬＤＦの技術を用いることなどにより、小型化が可能である上、非侵襲的且つ連続的な計測が可能である。このため、本技術の生体情報処理装置は、ウェアラブル装置に適しており、具体的にはウェアラブル血流計に適している。本技術を用いたウェアラブル血流計は、日常生活で発生する動作時においても頑健に血流速度の変化を検出することができる。
【００７０】
血流速度の変化は、人の健康状態や心理状態を表す指標として知られている。このため、本技術の生体情報処理装置は、健康状態及び／又は心理状態をセンシングするウェアラブル装置、アプリケーション、システムなどに用いることも可能である。
【００７１】
＜３．第１実施形態の第１変形例＞
図１２を参照し、第１実施形態の第１変形例に係る生体情報処理装置１Ａの動作について、上記第１実施形態と異なる点を説明する。
【００７２】
図１２は、第１実施形態の第１変形例に係る生体情報処理装置１Ａの機能構成の一例を示すブロック図である。本変形例の生体情報処理装置１Ａは、上記第１実施形態の心拍変動検出部２１に代えて拍動間隔検出処理部２１Ａを備える。拍動間隔検出処理部２１Ａは、制御部２０（図２）により実現される。
【００７３】
拍動間隔検出処理部２１Ａは、測定対象者の拍動間隔を検出する拍動間隔検出工程を実行し、拍動間隔を出力する。拍動間隔検出処理部２１Ａは、例えば、ＬＤＦから出力されるレーザー光の反射光の信号を用いて拍動間隔を算出し、出力する。
【００７４】
本変形例の生体情報処理装置１Ａは、上記第１実施形態で例示したＰＰＧやＥＣＧといった他の装置を必要としないため、構成をより単純化することが可能である。
【００７５】
＜４．第１実施形態の第２変形例＞
図１３を参照し、第１実施形態の第２変形例に係る生体情報処理装置１Ｂの動作について、上記第１実施形態と異なる点を説明する。
【００７６】
図１３は、第１実施形態の第２変形例に係る生体情報処理装置１Ｂの機能構成の一例を示すブロック図である。本変形例の生体情報処理装置１Ｂは、上記第１実施形態の心拍変動検出部２１に代えて拍動間隔データ保持部２１Ｂを備える。拍動間隔データ保持部２１Ｂは、制御部２０（図２）により実現される。
【００７７】
拍動間隔データ保持部２１Ｂは、ピーク検出部２４における過去の検出結果に基づいて拍動間隔を算出し保持する拍動間隔データ保持処理を実行する。ピーク検出部２４における過去の検出結果とは、例えば、ピーク検出部２４により過去に検出されたピーク時間、拍動成分強度などである。第１のフィルタ部２３及びピーク検出部２４は、拍動間隔データ保持部２１Ｂが保持する過去の拍動間隔データに基づいて処理を行う。
【００７８】
図１３に示す例では、ピーク検出部２４は、ピーク時間（ｔｐｋｓ）における拍動成分強度（ｔｐｋｓ，ＡＣｐｋｓ）を出力するにあたって、拍動間隔データ保持部２１Ｂが保持している過去（時間：ｔ−１）の拍動間隔のデータを用いる。
【００７９】
本変形例の生体情報処理装置１Ｂは、上記第１実施形態で例示したＰＰＧやＥＣＧといった他の装置を必要としないため、構成をより単純化することが可能である。
【００８０】
なお、本技術は以下のような構成も採ることができる。
〔１〕拍動間隔に基づいて第１の周波数帯域を選択し、血流速度信号の中から前記第１の周波数帯域の信号を通過させて、拍動成分信号を抽出する第１のフィルタ部を備える、生体情報処理装置。
〔２〕前記拍動成分信号のピーク位置を検出し、ピーク時間と該ピーク時間における拍動成分強度とを出力するピーク検出部を備え、
前記ピーク検出部は、
体動発生時以外においては、前記拍動間隔に基づいて前記拍動成分信号のピーク位置を検出して、拍動成分強度を履歴情報として保持し、
体動発生時においては、前記拍動間隔及び前記履歴情報に基づいて前記拍動成分信号のピーク位置を検出した後、前記履歴情報に基づいて前記ピーク位置の検出が誤りであるか否かを判定し、前記ピーク位置の検出が誤りである場合は、前記拍動間隔を用いてピーク位置を補正し、前記履歴情報を用いて拍動成分強度を補正する、〔１〕に記載の生体情報処理装置。
〔３〕前記第１のフィルタ部は、複数のバンドパスフィルタで構成される第１のフィルタバンクを備え、前記拍動間隔に基づいて前記第１のフィルタバンクから選択されたバンドパスフィルタを用いて前記拍動成分信号を抽出する、〔１〕又は〔２〕に記載の生体情報処理装置。
〔４〕前記ピーク検出部は、体動発生時以外において、前記拍動間隔に基づいて第２の周波数帯域を選択し、前記拍動成分信号の中から前記第２の周波数帯域の信号を通過させて、前記拍動成分信号のピーク位置を検出する第２のフィルタ部を備える、〔２〕に記載の生体情報処理装置。
〔５〕前記第２のフィルタ部は、複数のバンドパスフィルタで構成される第２のフィルタバンクを備え、体動発生時以外において、前記拍動間隔に基づいて前記第２のフィルタバンクから選択されたバンドパスフィルタを用いて前記拍動成分信号のピーク位置を検出する、〔４〕に記載の生体情報処理装置。
〔６〕前記ピーク検出部は、体動発生時において、前記ピーク位置の検出が誤りである場合は、前記拍動間隔により前記ピーク位置の確率分布の尤度を計算し、前記尤度が最大となるピーク位置を補正後のピーク位置とする、〔２〕に記載の生体情報処理装置。
〔７〕拍動間隔に基づいて第１の周波数帯域を選択し、血流速度信号の中から前記第１の周波数帯域の信号を通過させて、拍動成分信号を抽出する第１のフィルタ工程を含む、生体情報処理方法。
【符号の説明】
【００８１】
１，１Ａ，１Ｂ  生体情報処理装置
１０  血流センサ部
２０  制御部
２１  心拍変動検出部
２１Ａ  拍動間隔検出処理部
２１Ｂ  拍動間隔データ保持部
２２  血流速度計算部
２３  第１のフィルタ部
２３ａ  第１のフィルタバンク
２４  ピーク検出部
２４１  第２のフィルタ部
２４１ａ  第２のフィルタバンク
２５  血流速度ベースライン復元部
２６  包絡線検出部
【図１】
2020006014.tif 000004
【図２】
2020006014.tif 000005
【図３】
2020006014.tif 000006
【図４】
2020006014.tif 000007
【図５】
2020006014.tif 000008
【図６】
2020006014.tif 000009
【図７】
2020006014.tif 000010
【図８】
2020006014.tif 000011
【図９】
2020006014.tif 000012
【図１０】
2020006014.tif 000013
【図１１】
2020006014.tif 000014
【図１２】
2020006014.tif 000015
【図１３】
2020006014.tif 000016