JP6379899B2 - 情報処理装置、脈波計測プログラムおよび脈波計測方法 - Google Patents

Description

本発明は、情報処理装置等に関する。

脈波伝搬速度（ＰＷＶ：Pulse Wave Velocity）は、動脈における脈波の伝搬速度を示すものであり、動脈硬化等の心血管疾患の発見に役立つことが知られている。この脈波伝搬速度を測定する方法として、例えば、血管パターンにおける心拍同期信号の検出位置を２個所決定し、両検出位置間の距離と心拍同期信号の発生時間に基づいて脈波伝搬速度を測定する装置が知られている。かかる装置では、生体を拘束（静止）したうえで、心拍同期信号の検出位置を２個所決定して脈波伝搬速度を測定する。

また、脈波において、計測領域の輝度の時間変化によって当該計測領域の脈波信号を検出する技術も知られている。また、電極を生体に装着して、脈波伝搬速度を測定する技術も知られている。例えば、心電図法が挙げられる。

特開２０００−１８３８５４号公報 特開２００１−１５８４５号公報 特開２０１０−１２４２６６号公報 米国特許出願公開第２００４／０１７５１８０号明細書

しかしながら、上記の従来技術では、非接触で脈波伝搬時差を精度よく計測できないという問題がある。

例えば、上述した脈波伝搬速度は式（１）によって測定される。式（１）において、伝搬距離は血管長を示すものである。脈波伝搬時差は、式（２）で定義されるものであり、式（２）において、ｔ_１，ｔ_２は、心臓から各点への脈波伝搬時間を示す。

ＰＷＶ＝伝搬距離／脈波伝搬時差・・・（１）

脈波伝搬時差＝ｔ_１−ｔ_２・・・（２）

式（１）で算出する脈波伝搬速度において、計測量である伝搬距離および脈波伝搬時差が大きな値となるほど、計測機の分解能や雑音の影響が小さくなるため、精度が高くなる。また、脈波伝搬速度の変化に着目する場合には、伝搬距離については血管の長さであり短期的な変化はなく、固定値として扱えるため、脈波伝搬時差の計測量が、脈波伝搬速度の精度に影響を及ぼす。

ここで、脈波伝搬速度の変化に着目して、心血管疾患の発見するために、日常的、かつ簡便に脈波伝搬速度を求めるためには、非接触で利用者を拘束しないことが好ましい。しかし、非接触で精度よく脈波伝搬時差を計測する技術は確立されていない。

１つの側面では、非接触で脈波伝搬時差を精度よく計測できる情報処理装置、脈波計測プログラムおよび脈波計測方法を提供することを目的とする。

第１の案では、情報処理装置は、分割部と、検出部と、特定部とを有する。分割部は、生体の表面が撮影された画像を複数の領域に分割する。検出部は、複数の領域毎に、領域に含まれる画素の時間変化に伴う輝度変化に基づいて脈波信号を検出する。特定部は、各領域の脈波信号をそれぞれ比較して各領域の脈波信号の遅延量をそれぞれ測定し、測定した各遅延量のうち最大の遅延量を特定する。

本発明の１実施態様によれば、非接触で脈波伝搬時差を精度よく計測することができる。

図１は、本実施例１に係る情報処理装置の構成を示す機能ブロック図である。 図２は、評価領域の一例を示す図である。 図３は、分割領域の一例を示す図である。 図４は、脈波伝搬距離情報のデータ構造の一例を示す図である。 図５は、脈波情報のデータ構造の一例を示す図である。 図６は、本実施例１に係る特定部の処理を説明するための図である。 図７は、分割領域の脈波信号の一例を示す図である。 図８は、本実施例１に係る情報処理装置の初期学習の処理手順を示すフローチャート（１）である。 図９は、本実施例１に係る情報処理装置の初期学習の処理手順を示すフローチャート（２）である。 図１０は、本実施例１に係る情報処理装置が長期計測を行う場合の処理手順を示すフローチャートである。 図１１は、本実施例１に係る出力表示処理の処理手順を示すフローチャートである。 図１２は、本実施例２に係る情報処理装置の構成を示す機能ブロック図である。 図１３は、本実施例２に係る特定部の処理を説明するための図である。 図１４は、領域集合の脈波信号の一例を示す図である。 図１５は、本実施例２に係る情報処理装置の初期学習の処理手順を示すフローチャート（１）である。 図１６は、本実施例２に係る情報処理装置の初期学習の処理手順を示すフローチャート（２）である。 図１７は、本実施例２に係る情報処理装置が長期計測を行う場合の処理手順を示すフローチャートである。 図１８は、本実施例３に係る情報処理装置の構成を示す機能ブロック図である。 図１９は、除去部の処理を説明するための図である。 図２０は、本実施例３に係る情報処理装置が長期計測を行う場合の処理手順を示すフローチャートである。 図２１は、脈波計測プログラムを実行するコンピュータの一例を示す図である。

以下に、本願の開示する情報処理装置、脈波計測プログラムおよび脈波計測方法の実施例を図面に基づいて詳細に説明する。なお、この実施例によりこの発明が限定されるものではない。

本実施例１に係る情報処理装置の構成の一例について説明する。図１は、本実施例１に係る情報処理装置の構成を示す機能ブロック図である。図１に示すように、この情報処理装置１００は、撮像機能部１１０と、入力部１２０と、表示部１３０と、通信部１３５と、記憶部１４０と、制御部１５０とを有する。

撮像機能部１１０は、ＣＣＤ（Charge Coupled Device）やＣＭＯＳ（Complementary Metal Oxide Semiconductor）などの撮像素子を用いた撮像装置である。例えば、撮像機能部１１０は、Ｒ（red）、Ｇ（green）、Ｂ（blue）など３種以上の受光素子を搭載することができる。撮像機能部１１０は、ユーザの顔の画像を撮像し、撮像した画像情報を、制御部１５０に出力する。撮像する部位は、光学的に脈波信号を観測でき、輪郭や皺等、特徴が判別出来る部位であれば、顔以外であっても良い。以下の説明において、適宜、ユーザの顔の画像情報を、生体画像情報と表記する。

入力部１２０は、各種の情報を情報処理装置１００に入力するための入力装置である。例えば、入力部１２０は、キーボードやマウス、タッチパネル、各種入力ボタンに対応する。

表示部１３０は、各種の情報を表示する表示装置である。例えば、表示部１３０は、ＰＷＶの算出結果などを出力する。表示部１３０は、モニタやディスプレイ、タッチパネル等に対応する。

通信部１３５は、ネットワークを介して他の外部装置とデータ通信を実行する処理部である。外部装置の図示を省略する。通信部１３５は、通信装置に対応する。

記憶部１４０は、画像バッファ１４１、生体特徴情報１４２、評価領域情報１４３、脈波伝搬距離情報１４４、脈波情報１４５を有する。記憶部１４０は、例えば、ＲＡＭ（Random Access Memory）、ＲＯＭ（Read Only Memory）、フラッシュメモリ（Flash Memory）などの半導体メモリ素子などの記憶装置に対応する。

画像バッファ１４１は、撮像機能部１１０から取得する各フレームの生体画像情報を蓄積する記憶領域である。例えば、画像バッファ１４１は、時系列に並んだ各フレームの生体画像情報を有する。

生体特徴情報１４２は、ユーザのどの部分を特徴物とするのかを定義する情報である。本実施例１では一例として、生体特徴情報１４２には、ユーザの両目と鼻とを特徴物とする旨の情報が含まれているものとする。また、生体特徴情報１４２は、特徴物と、かかる特徴物の特徴を示す部分テンプレートの情報を対応付けて記憶してもよい。

評価領域情報１４３は、例えば、評価領域の情報と、分割領域の情報を有する。図２は、評価領域の一例を示す図である。例えば、特徴点を両目１０ａ，１０ｂと鼻１０ｃとした場合の評価領域２０を示している。例えば、評価領域２０は、両目１０ａ，１０ｂを結ぶ線分を一辺とし、その辺の対辺の中点を鼻１０ｃとする平行四辺形に対応する。後述する評価領域設定部１５３は、評価領域情報１４３の評価領域２０に従って、生体画像情報上に、評価領域を設定する。

図３は、分割領域の一例を示す図である。図３において、両目１０ａ，１０ｂ間の辺とその対辺をＭ等分、残りの２辺をＮ等分し、Ｍ×Ｎ個の等面積の平行四辺形に分割した各領域が、分割領域となる。例えば、図３に示す例では、評価領域２０がＭ×Ｎ等分にされ、（１、１）〜（Ｍ、Ｎ）の分割領域が存在する。後述する分割部１５４は、評価領域情報１４３の分割領域に従って、生体画像情報上の評価領域を、複数の分割領域に分割する。

脈波伝搬距離情報１４４は、ユーザのある２点の座標と、脈波伝搬距離とが対応付けられた情報である。図４は、脈波伝搬距離情報のデータ構造の一例を示す図である。図４に示すように、この脈波伝搬距離情報は、始点座標と、終点座標と、脈波伝搬距離とを対応付ける。始点座標、終点座標は、生体画像情報上の分割領域を特定する座標番号を示し、脈脈伝搬距離は、心臓から始点座標までと、心臓から終点座標までの脈波伝搬距離の差である。例えば、図４の脈波伝搬距離情報１４４の１行目において、始点座標（ｘｓ１、ｙｓ１）から終点座標（ｘｅ１、ｙｅ１）までの脈波伝搬距離の差は、Ｌ１であることが示されている。

脈波情報１４５は、各分割領域から抽出された脈波信号の情報を示す。図５は、脈波情報のデータ構造の一例を示す図である。図５に示すように、この脈波情報１４５は、分割領域識別情報と、脈波信号とを対応付ける。分割領域識別情報は、分割領域を一意に識別する情報である。脈波信号は、分割領域から特定される脈波信号である。

図１の説明に戻る。制御部１５０は、受付部１５１と、特徴点座標決定部１５２と、評価領域設定部１５３と、分割部１５４と、脈波信号検出部１５５と、特定部１５６と、ＰＷＶ算出部１５７とを有する。制御部１５０は、例えば、ＡＳＩＣ（Application Specific Integrated Circuit）や、ＦＰＧＡ（Field Programmable Gate Array）などの集積装置に対応する。また、制御部１５０は、例えば、ＣＰＵ（Central Processing Unit）やＭＰＵ（Micro Processing Unit）等の電子回路に対応する。

受付部１５１は、撮像機能部１１０、入力部１２０、通信部１３５から受け付けた情報を、記憶部１４０に格納する処理部である。受付部１５１は、撮像機能部１１０から生体画像情報を取得し、取得した生体画像情報を、画像バッファ１４１に登録する。

受付部１５１は、入力部１２０から、生体特徴情報を取得し、取得した生体特徴情報を、生体特徴情報１４２に登録する。受付部１５１は、入力部１２０から、評価領域情報を取得し、取得した評価領域情報を、評価領域情報１４３に登録する。受付部１５１は、入力部１２０から、脈波伝搬距離情報を取得し、取得した脈波伝搬距離情報を、脈波伝搬距離情報１４４に登録する。なお、受付部１５１は、生体特徴情報、評価領域情報、脈波伝搬距離情報を通信部１３５から取得してもよい。

特徴点座標決定部１５２は、生体特徴情報１４２を基にして、画像バッファ１４１に蓄積された生体画像情報から特徴点の座標を決定する処理部である。特徴点座標決定部１５２は、生体画像情報１４２に含まれる特徴物の名称に対応する特徴物を、生体画像情報から認識する。特徴点座標決定部１５２は、認識した特徴物の中心に位置する特徴点の座標を決定する。特徴点座標決定部１５２は、生体画像情報と、決定した特徴点の座標とを対応付けて、評価領域設定部１５３に出力する。特徴点座標検出部１５２は、画像バッファ１４１に格納された生体画像情報を、時系列順に順次取得し、上記処理を繰り返し実行する。

評価領域設定部１５３は、特徴点座標決定部１５２から受け付けた情報を基にして、生体画像情報上に評価領域を設定する処理部である。評価領域設定部１５３は、評価領域情報１４３を参照し、どのように評価領域を生体画像情報上に設定するのかを判定する。例えば、評価領域設定部１５３は、図２で説明したように、両目を結ぶ線分を一辺とし、その辺の対辺の中心を鼻とする平行四辺形を、評価領域に設定する。評価領域設定部１５３は、生体画像情報と、この生体画像情報上に設定した評価領域の情報を、分割部１５４に出力する。評価領域設定部１５３は、特徴点座標決定部１５２から情報を受け付ける度に、上記処理を繰り返し実行する。

分割部１５４は、評価領域設定部１５３から受け付けた情報を基にして、生体画像情報上の評価領域を分割し、分割領域を決定する処理部である。分割部１５４は、評価領域情報１４３を参照し、どのように評価領域を複数の分割領域に分割するのかを判定する。例えば、分割部１５４は、図３で説明したように、両目間の辺とその対辺をＭ等分し、残りの２辺をＮ等分することで、（１、１）〜（Ｍ、Ｎ）の分割領域を生体画像情報上に設定する。分割部１５４は、生体画像情報と、この生体画像情報上に設定した分割領域の情報を、脈波信号検出部１５５に出力する。分割部１５４は、評価領域設定部１５３から情報を受け付ける度に、上記処理を繰り返し実行する。

脈波信号検出部１５５は、検出部の一例である。脈波信号検出部１５５は、分割部１５４から受け付けた情報を基にして、各分割領域に対応する脈波信号を検出する処理部である。脈波信号検出部１５５は、分割領域と、脈波信号とを対応付けて、脈波情報１４５に登録する。以下において、脈波信号検出部１５５の処理の一例について説明する。

脈波信号検出部１５５は、生体画像情報上の全分割領域に含まれる画素の輝度を、分割領域毎に平均することで、Ｍ×Ｎ個の輝度を算出する。脈波信号検出部１５５は、Ｍ×Ｎ個の輝度をそれぞれ該当する分割領域に対応付けて脈波情報１４５に記憶し、脈波信号（ｍ、ｎ）とする。脈波信号検出部１５５は、分割部１５４から情報を受け付ける度に、上記処理を繰り返し実行する。分割領域毎に、抽出した輝度をつなげていくことで、時間と輝度との関係を示す脈波信号が形成されていく。

特定部１５６は、脈波情報１４５を基にして、各分割領域の脈波信号をそれぞれ比較して、各分割領域の遅延量をそれぞれ測定し、測定した各遅延量のうち、最大の遅延量を特定する。以下において、特定部１５６の処理の一例について説明する。

図６は、本実施例１に係る特定部の処理を説明するための図である。図６に示す例では説明の便宜上、各分割領域を分割領域１〜Ｎとして説明する。図７は、分割領域の脈波信号の一例を示す図である。図７において各縦軸は脈波信号の大きさに対応する軸であり、横軸は時刻に対応する軸である。図７では一例として、分割領域２の脈波信号を、脈波信号３０ａとする。分割領域６の脈波信号を、脈波信号３０ｂとする。分割領域Ｎの脈波信号を、脈波信号３０ｃとする。

例えば、特定部１５６は、各脈波信号の強さが極大となる時刻をそれぞれ特定し、特定した極大となる時刻の差を、遅延量として特定する。なお、特定部１５６は、同一の脈波信号の組について、複数の遅延量が存在する場合には、最小の遅延量を、かかる脈波信号の組の遅延量として算出する。また、一定期間の脈波信号において、一方の脈波信号に時差をつけて相関係数を求め、相関係数が最大となる時差を遅延量とする手段などでも良い。

特定部１５６は、図７に示す例において、分割領域２の脈波信号３０ａと、分割領域６の脈波信号３０ｂとの遅延量を、０．５ｓとして算出する。特定部１５６は、分割領域６の脈波信号３０ｂと、分割領域Ｎの脈波信号３０ｃとの遅延量を、０．３５ｓとして算出する。特定部１５６は、分割領域２の脈波信号３０ｂと、分割領域Ｎの脈波信号３０ｃとの遅延量を、０．１５ｓとして算出する。特定部１５６は、その他の分割領域についても同様にして、遅延量を算出する。

特定部１５６は、各分割領域の組の遅延量のうち、最大の遅延量と、係る遅延量に対応する分割領域の組を特定する。例えば、図７に示す例では、分割領域２の脈波信号３０ａと、分割領域６の脈波信号３０ｂとの遅延量「０．５ｓ」が、他の分割領域の組の遅延量と比較した最大の遅延量となっている。特定部１５６は、特定した最大の遅延量と、係る遅延量に対応する分割領域の組の情報を、ＰＷＶ算出部１５７に出力する。

ここで、より具体的には、上述した特徴点座標決定部１５２、評価領域設定部１５３、分割部１５４、脈波信号検出部１５５、特定部１５６は、初期学習動作と、長期計測動作の２段階で動作を行う。

初期学習動作において、各処理部１５２〜１５６は、所定の回数、上記処理を繰り返し実行し、特定部１５６は、遅延量が最大となる分割領域の組を、順次記憶する。特定部１５６は、遅延量が最大となる分割領域の組のうち、出現頻度が最大となる分割領域の組を、「評価対象の分割領域の組」として特定する。特定部１５６は、特定した評価対象の分割領域の組の情報を、評価領域情報１４３に登録する。例えば、特定部１５６は、（１、１）〜（Ｍ、Ｎ）の分割領域のうち、分割領域（ｍ１、ｎ１）、（ｍ２、ｎ２）を、評価対象の分割領域とする。この場合には、かかる分割領域（ｍ１、ｎ１）、（ｍ２、ｎ２）が、評価対象の分割領域の組である旨を、評価領域情報１４３に登録する。

続いて、特定部１５６は、長期計測動作において、評価領域情報１４３に登録された評価対象の分割領域の組（ｍ１、ｎ１）、（ｍ２、ｎ２）に着目し、評価対象の分割領域の組の脈波信号の遅延量を特定する。特定部１５６は、特定した脈波信号の遅延量と、評価対象の分割領域の情報を、ＰＷＶ算出部１５７に出力する。

ＰＷＶ算出部１５７は、決定部の一例である。ＰＷＶ算出部１５７は、特定部１５６から受け付けた情報と、脈波伝搬距離情報１４４とを基にして、脈波伝搬速度（ＰＷＶ）を算出する処理部である。具体的に、ＰＷＶ算出部１５７は、式（３）を利用する。式（３）において、脈波伝搬距離は、最大の遅延量に対応する分割領域の組の各座標と、脈波伝搬距離情報１４４の始点座標、終点座標の組とを比較して、特定される距離である。分割領域の座標は、例えば、分割座標の中心座標とする。ＰＷＶ算出部１５７は、最大の遅延量に対応する分割領域の組の各座標と、脈波伝搬距離情報１４４の座標の組とを比較し、最も類似する始点座標、終点座標の組に対応するレコードの脈波伝搬距離を特定する。

ＰＷＶ＝脈波伝搬距離／最大の遅延量・・・（３）

ＰＷＶ算出部１５７は、今回算出したＰＷＶの値と、前回算出したＰＷＶの値とを比較し、差分の絶対値が、閾値よりも大きいか否かを判定する。ＰＷＶ算出部１５７は、差分の絶対値が、閾値よりも大きい場合には、表示部１３０に、今回算出したＰＷＶの値を表示させると共に、警告表示を行う。ＰＷＶ算出部１５７は、差分の絶対値が、閾値以下である場合には、今回算出したＰＷＶの値を、表示部１３０に表示させる。

次に、本実施例１に係る情報処理装置１００の処理手順について説明する。図８および図９は、本実施例１に係る情報処理装置の初期学習の処理手順を示すフローチャートである。図８に示すように、情報処理装置１００の受付部１５１は、生体特徴情報を取得し、記憶部１４０に格納する（ステップＳ１０１）。受付部１５１は、評価領域情報を取得し、記憶部１４０に格納する（ステップＳ１０２）。受付部１５１は、脈波伝搬距離情報を取得し、記憶部１４０に格納する（ステップＳ１０３）。受信部１５１は、タイマを設定する（ステップＳ１０４）。ステップＳ１０４において、例えば、受付部１５１は、タイマ長Ｔ１、タイマｔ１＝０、タイマ長Ｔ２、タイマｔ２＝０をそれぞれ設定する。

情報処理装置１００の撮像機能部１１０は、生体画像を撮影し、受付部１５１に出力する（ステップＳ１０５）。ステップＳ１０５において、受付部１５１は、生体画像情報を、画像バッファ１４１に蓄積する。情報処理装置１００の特徴点座標決定部１５２は、特徴点の座標を決定する（ステップＳ１０６）。情報処理装置１００の評価領域設定部１５３は、評価領域を決定する（ステップＳ１０７）。情報処理装置１００の分割部１５４は、分割領域を決定し（ステップＳ１０８）、図９のステップＳ１０９に移行する。

図９について説明する。情報処理装置１００の脈波信号検出部１５５は、全分割領域に対応する輝度を算出する（ステップＳ１０９）。脈波信号検出部１５５は、脈波信号を検出し、脈波信号を脈波情報１４５に格納する（ステップＳ１１０）。脈波信号検出部１５５は、ｔ１＝ｔ１＋１により、ｔ１の値を更新する（ステップＳ１１１）。

脈波信号検出部１５５は、ｔ１の値とＴ１の値とが等しいか否かを判定する（ステップＳ１１２）。脈波信号検出部１５５は、ｔ１の値とＴ１の値とが等しくない場合に（ステップＳ１１２，Ｎｏ）、図８のステップＳ１０５に移行する。

一方、ｔ１の値とＴ１の値とが等しい場合には（ステップＳ１１２，Ｙｅｓ）、情報処理装置１００の特定部１５６は、各分割領域の脈波信号を比較し、遅延量が最大となる分割領域の組を特定する（ステップＳ１１３）。特定部１５６は、特定した分割領域の組を記憶する（ステップＳ１１４）。

特定部１５６は、ｔ２＝ｔ２＋１により、ｔ２の値を更新する（ステップＳ１１５）。特定部１５６は、ｔ２の値とＴ２の値とが等しいか否かを判定する（ステップＳ１１６）。特定部１５６は、ｔ２の値とＴ２の値とが等しくない場合には（ステップＳ１１６，Ｎｏ）、図８のステップＳ１０５に移行する。

一方、特定部１５６は、ｔ２の値とＴ２の値とが等しい場合には（ステップＳ１１６，Ｙｅｓ）、ステップＳ１１７に移行する。この場合には、Ｔ２個の分割領域の組が、特定部１５６に記憶されていることになる。特定部１５６は、Ｔ２個の分割領域の組に基づいて、評価領域情報１４３を更新する（ステップＳ１１７）。例えば、ステップＳ１１７において、特定部１５６は、Ｔ２個の分割領域の組のうち、出現頻度が最大となる分割領域の組を、「評価対象の分割領域の組」として、評価領域情報１４３に設定する。

次に、本実施例１に係る情報処理装置１００が長期計測を行う場合の処理手順について説明する。図１０は、本実施例１に係る情報処理装置が長期計測を行う場合の処理手順を示すフローチャートである。図１０に示すように、情報処理装置１００の受付部１５１は、タイマを設定する（ステップＳ１５１）。ステップＳ１５１において、例えば、受付部１５１は、タイマ長Ｔ１、タイマｔ１＝０をそれぞれ設定する。

情報処理装置１００の撮像機能部１１０は、生体画像を撮影し、受付部１５１に出力する（ステップＳ１５２）。ステップＳ１５２において、受付部１５１は、生体画像情報を、画像バッファ１４１に蓄積する。情報処理装置１００の特徴点座標決定部１５２は、特徴点の座標を決定する（ステップＳ１５３）。情報処理装置１００の評価領域設定部１５３は、評価領域を決定する（ステップＳ１５４）。情報処理装置１００の分割部１５４は、評価領域を分割し、情報処理装置１００の脈波信号検出部１５５は、評価領域情報１４３を基にして、評価対象の分割領域の組を決定する（ステップＳ１５５）。

脈波信号検出部１５５は、評価対象の分割領域の組の輝度を算出する（ステップＳ１５６）。脈波信号検出部１５５は、評価対象の分割領域の組の脈波信号を特定し、特定した脈波信号を脈波情報１４５に格納する（ステップＳ１５７）。脈波信号検出部１５５は、ｔ１＝ｔ１＋１により、ｔ１の値を更新する（ステップＳ１５８）。

脈波信号検出部１５５は、ｔ１の値とＴ１の値とが等しいか否かを判定する（ステップＳ１５９）。脈波信号検出部１５５は、ｔ１の値とＴ１の値とが等しくない場合には（ステップＳ１５９，Ｎｏ）、ステップＳ１５１に移行する。

一方、ｔ１の値とＴ１の値とが等しい場合には（ステップＳ１５９，Ｙｅｓ）、情報処理装置１００の特定部１５６は、評価対象の分割領域の脈波信号を比較して、遅延量を算出する（ステップＳ１６０）。情報処理装置１００のＰＷＶ算出部１５７は、遅延量と脈波伝搬距離情報１４４とを基にして、ＰＷＶを算出する（ステップＳ１６１）。ＰＷＶ算出部１５７は、出力表示処理を行う（ステップＳ１６２）。

続いて、図１０のステップＳ１６２に示した出力表示処理の一例について説明する。図１１は、本実施例１に係る出力表示処理の処理手順を示すフローチャートである。図１１に示すように、情報処理装置１００のＰＷＶ算出部１５７は、１回目の処理であるか否かを判定する（ステップＳ１７１）。ＰＷＶ算出部１５７は、１回目の処理でない場合には（ステップＳ１７１，Ｎｏ）、ステップＳ１７３に移行する。

一方、ＰＷＶ算出部１５７は、１回目の処理である場合には（ステップＳ１７１，Ｙｅｓ）、ＰＷＶの記憶値Ｐｂを初期化する（ステップＳ１７２）。ステップＳ１７２において、ＰＷＶ計算部１５７は、「Ｐｂ」の値を、今回、ＰＷＶ算出部１５７が算出したＰＷＶの値「Ｐ」に設定する。

ＰＷＶ算出部１５７は、｜Ｐ−Ｐｂ｜の値が閾値ＴＨよりも大きいか否かを判定する（ステップＳ１７３）。ＰＷＶ算出部１５７は、｜Ｐ−Ｐｂ｜の値が閾値ＴＨよりも大きくない場合には（ステップＳ１７３，Ｎｏ）、ステップＳ１７５に移行する。

一方、ＰＷＶ算出部１５７は、｜Ｐ−Ｐｂ｜の値が閾値ＴＨよりも大きい場合には（ステップＳ１７３，Ｙｅｓ）、表示部１３０に警報表示させる（ステップＳ１７４）。ＰＷＶ算出部１５７は、表示部１３０にＰＷＶの現状値Ｐを表示させる（ステップＳ１７５）。ＰＷＶ算出部１５７は、ＰＷＶの記憶値Ｐｂを現状値Ｐによって更新する（ステップＳ１７６）。なお、図１１に示す処理は、スタートからエンドまで繰り返し実行される。

次に、本実施例１に係る情報処理装置１００の効果について説明する。情報処理装置１００は、生体の表面が撮影された生体画像情報を複数の分割領域に分割し、各分割領域に含まれる画素の時間変化に伴う輝度変化に基づいて脈波信号を検出する。そして、情報処理装置１００は、各分割領域の脈波信号をそれぞれ比較して各分割領域の脈波信号の遅延量をそれぞれ測定し、測定した各遅延量のうち最大の遅延量を特定する。このように、情報処理装置１００は、最大の遅延量を特定することで、ＰＷＶを計測する場合の精度を向上することができる。例えば、ＰＷＶは、式（１）（式（３））によって算出され、この式（１）の分母および分子共に、より大きな値を用いると、精度が向上するためである。

次に、本実施例２に係る情報処理装置２００の構成の一例について説明する。図１２は、本実施例２に係る情報処理装置の構成を示す機能ブロック図である。図１２に示すように、この情報処理装置２００は、撮像機能部２１０と、入力部２２０と、表示部２３０と、通信部２３５と、記憶部２４０と、制御部２５０とを有する。図１２において、撮像機能部２１０、入力部２２０、表示部２３０、通信部２３５に関する説明は、実施例１に示した撮像機能部１１０、入力部１１０、表示部１３０、通信部１３５と同様であるため、説明を省略する。

記憶部２４０は、画像バッファ２４１、生体特徴情報２４２、評価領域情報２４３、脈波伝搬距離情報２４４、脈波情報２４５を有する。記憶部２４０は、例えば、ＲＡＭ、ＲＯＭ、フラッシュメモリなどの半導体メモリ素子などの記憶装置に対応する。画像バッファ２４１、生体特徴情報２４２、評価領域情報２４３、脈波伝搬距離情報２４４、脈波情報２４５に関する説明は、図１に示した画像バッファ１４１、生体特徴情報１４２、評価領域情報１４３、脈波伝搬距離情報１４４、脈波情報１４５に関する説明と同様である。

制御部２５０は、受付部２５１、特徴点座標決定部２５２、評価領域設定部２５３、分割部２５４、脈波信号検出部２５５、特定部２５６、ＰＷＶ算出部２５７を有する。制御部２５０は、例えば、ＡＳＩＣや、ＦＰＧＡなどの集積装置に対応する。また、制御部２５０は、例えば、ＣＰＵやＭＰＵ等の電子回路に対応する。このうち、受付部２５１、特徴点座標決定部２５２、評価領域設定部２５３、分割部２５４、脈波信号検出部２５５に関する説明は、図１に示した受付部１５１、特徴点座標決定部１５２、評価領域設定部１５３、分割部１５４、脈波信号検出部１５５に関する説明と同様である。

特定部２５６は、脈波情報２４５を基にして、各分割領域の脈波信号をそれぞれ比較して、位相の差が所定値未満となる各分割領域を統合する。以下の説明において、統合した後の各分割領域を、領域集合と表記する。この領域集合には、複数の分割領域が含まれる場合もあるし、単一の分割領域しか含まれない場合もある。

特定部２５６は、領域集合の各脈波信号を平均化することで、領域集合に対応する脈波信号を決定する。特定部２５６は、各領域集合の脈波信号を比較して、各領域集合の遅延量をそれぞれ測定し、測定した各遅延量のうち、最大の遅延量を特定する。以下において、特定部２５６の処理の一例について説明する。

図１３は、本実施例２に係る特定部の処理を説明するための図である。図１３に示す例では説明の便宜上、各分割領域を分割領域１〜Ｎとして説明する。例えば、分割領域２，３，５の脈波信号の位相の差が所定値未満であるとすると、特定部２５６は、分割領域２，３，５を統合して、領域集合４０ａを生成する。特定部２５６は、分割領域２，３，５の脈波信号を平均化したものを、領域集合４０ａの脈波信号とする。

図１３において、例えば、分割領域６，７の脈波信号の位相の差が所定値未満であるとすると、特定部２５６は、分割領域６，７を統合して、領域集合４０ｂを生成する。特定部２５６は、分割領域６，７の脈波信号を平均化したものを、領域集合４０ｂの脈波信号とする。

図１３において、例えば、分割領域Ｎ−６，Ｎの脈波信号の位相の差が所定値未満であるとすると、特定部２５６は、分割領域Ｎ−６，Ｎを統合して、領域集合４０ｃを生成する。特定部２５６は、分割領域Ｎ−６，Ｎの脈波信号を平均化したものを、領域集合４０ｃの脈波信号とする。

図１４は、領域集合の脈波信号の一例を示す図である。図１４において各縦軸は脈波信号の大きさに対応する軸であり、横軸は時刻に対応する軸である。図１４では一例として、領域集合４０ａの脈波信号を、脈波信号５０ａとする。領域集合４０ｂの脈波信号を、脈波信号５０ｂとする。領域集合４０ｃの脈波信号を、脈波信号５０ｃとする。

例えば、特定部２５６は、各領域集合において、各脈波信号の強さが極大となる時刻をそれぞれ特定し、特定した極大となる時刻の差を、遅延量として特定する。なお、特定部２５６は、同一の脈波信号の組について、複数の遅延量が存在する場合には、最小の遅延量を、かかる脈波信号の組の遅延量として算出する。

特定部２５６は、図１４に示す例において、領域集合４０ａの脈波信号５０ａと、領域集合４０ｂの脈波信号５０ｂとの遅延量を、０．５ｓとして算出する。特定部２５６は、領域集合４０ｂの脈波信号５０ｂと、領域集合４０ｃの脈波信号５０ｃとの遅延量を、０．３５ｓとして算出する。特定部２５６は、領域集合４０ａの脈波信号５０ａと、領域集合４０ｃの脈波信号５０ｃとの遅延量を、０．１５ｓとして算出する。特定部２５６は、その他の分割領域についても同様にして、遅延量を算出する。

特定部２５６は、各領域集合の組の遅延量のうち、最大の遅延量と、係る遅延量に対応する領域集合の組を特定する。例えば、図１４に示す例では、領域集合４０ａの脈波信号５０ａと、領域集合４０ｂの脈波信号５０ｂとの遅延量「０．５ｓ」が、他の領域集合の組の遅延量と比較した最大の遅延量となっている。特定部２５６は、特定した最大の遅延量と、係る遅延量に対応する領域集合の組の情報を、ＰＷＶ算出部２５７に出力する。

ここで、より具体的には、上述した特徴点座標決定部２５２、評価領域設定部２５３、分割部２５４、脈波信号検出部２５５、特定部２５６は、初期学習動作と、長期計測動作の２段階で動作を行う。

初期学習動作において、各処理部２５２〜２５６は、所定の回数、上記処理を繰り返し実行し、特定部２５６は、遅延量が最大となる領域集合の組を、順次記憶する。特定部２５６は、遅延量が最大となる領域集合の組のうち、出現頻度が最大となる領域集合の組を、「評価対象の領域集合の組」として特定する。特定部２５６は、特定した評価対象の領域集合の組の情報を、評価領域情報２４３に登録する。

続いて、特定部２５６は、長期計測動作において、評価領域情報２４３に登録された評価対象の領域集合の組に着目し、評価対象の領域集合の組の脈波信号の遅延量を特定する。特定部２５６は、特定した脈波信号の遅延量と、評価対象の領域集合の情報を、ＰＷＶ算出部２５７に出力する。

ＰＷＶ算出部２５７は、特定部２５６から受け付けた情報と、脈波伝搬距離情報２４４とを基にして、脈波伝搬速度（ＰＷＶ）を算出する処理部である。具体的に、ＰＷＶ算出部２５７は、実施例１に示した式（３）を利用する。式（３）において、脈波伝搬距離は、最大の遅延量に対応する領域集合に含まれる分割領域の組の各座標と、脈波伝搬距離情報２４４の始点座標、終点座標の組とを比較して、特定される距離である。例えば、領域集合の組が図１３に示した領域集合４０ａと、領域集合４０ｂとの組であるものとする。この場合には、ＰＷＶ計算部２５７は、分割領域２，３，５のいずれかの分割領域と、分割領域６、７のいずれかの分割領域とを特定し、特定した分割領域の座標と、脈波伝搬距離情報２４４とを基にして、脈波伝搬距離を特定する。なお、ＰＷＶ算出部２５７は、脈波伝搬距離が最大となる、分割領域の組合せを、領域集合の組から検索し、利用してもよい。ＰＷＶ算出部２５７に関するその他の処理は、図１に示したＰＷＶ算出部１５７と同様である。

次に、本実施例２に係る情報処理装置２００の処理手順について説明する。図１５および図１６は、本実施例２に係る情報処理装置の初期学習の処理手順を示すフローチャートである。図１５に示すように、情報処理装置２００の受付部２５１は、生体特徴情報を取得し、記憶部２４０に格納する（ステップＳ２０１）。受付部２５１は、評価領域情報を取得し、記憶部２４０に格納する（ステップＳ２０２）。受付部２５１は、脈波伝搬距離情報を取得し、記憶部２４０に格納する（ステップＳ２０３）。受信部２５１は、タイマを設定する（ステップＳ２０４）。ステップＳ２０４において、例えば、受付部２５１は、タイマ長Ｔ１、タイマｔ１＝０、タイマ長Ｔ２、タイマｔ２＝０をそれぞれ設定する。

情報処理装置２００の撮像機能部２１０は、生体画像を撮影し、受付部２５１に出力する（ステップＳ２０５）。ステップＳ２０５において、受付部２５１は、生体画像情報を、画像バッファ２４１に蓄積する。情報処理装置２００の特徴点座標決定部２５２は、特徴点の座標を決定する（ステップＳ２０６）。情報処理装置２００の評価領域設定部２５３は、評価領域を決定する（ステップＳ２０７）。情報処理装置２００の分割部２５４は、分割領域を決定し（ステップＳ２０８）、図１６のステップＳ２０９に移行する。

図１６について説明する。情報処理装置２００の脈波信号検出部２５５は、全分割領域に対応する輝度を算出する（ステップＳ２０９）。脈波信号検出部２５５は、脈波信号を検出し、脈波信号を脈波情報２４５に格納する（ステップＳ２１０）。脈波信号検出部２５５は、ｔ１＝ｔ１＋１により、ｔ１の値を更新する（ステップＳ２１１）。

脈波信号検出部２５５は、ｔ１の値とＴ１の値とが等しいか否かを判定する（ステップＳ２１２）。脈波信号検出部２５５は、ｔ１の値とＴ１の値とが等しくない場合に（ステップＳ２１２，Ｎｏ）、図１５のステップＳ２０５に移行する。

一方、ｔ１の値とＴ１の値とが等しい場合には（ステップＳ２１２，Ｙｅｓ）、情報処理装置２００の特定部２５６は、各分割領域の脈波信号を比較して遅延量を算出し、記憶する（ステップＳ２１３）。特定部２５６は、ｔ２＝ｔ２＋１により、ｔ２の値を更新する（ステップＳ２１４）。

特定部２５６は、ｔ２の値とＴ２の値とが等しいか否かを判定する（ステップＳ２１５）。特定部２５６は、ｔ２の値とＴ２の値とが等しくない場合には（ステップＳ２１５，Ｎｏ）、図１５のステップＳ２０５に移行する。

一方、特定部２５６は、ｔ２の値とＴ２の値とが等しい場合には（ステップＳ２１５，Ｙｅｓ）、ステップＳ２１６に移行する。この場合には、Ｔ２個の分割領域の組が、特定部２５６に記憶されていることになる。特定部２５６は、分割領域の組合せ毎に、遅延量の統計を取る（ステップＳ２１６）。ステップＳ２１６において、特定部２５６は、例えば、Ｔ２個の遅延量の平均を、統計値として算出する。

特定部２５６は、統計値が閾値α未満であるか否かを判定する（ステップＳ２１７）。特定部２５６は、統計値が閾値α未満でない場合には（ステップＳ２１７，Ｎｏ）、評価領域情報２４３を更新し（ステップＳ２１８）、処理を終了する。例えば、ステップＳ２１８において、特定部２５６は、Ｔ２個の分割領域の組のうち、遅延量が最大となる分割領域の組を、「評価対象の領域集合の組」として、評価領域情報２４３に設定する。または、特定部２５６は、分割領域が統合されている場合には、遅延量が最大となる領域集合の組を、「評価対象の領域集合の組」として、評価領域情報２４３に設定する。

一方、特定部２５６は、統計値が閾値α未満である場合には（ステップＳ２１７，Ｙｅｓ）、統合前の分割領域の数と、領域集合の数との合計が３以上であるか否かを判定する（ステップＳ２１９）。特定部２５６は、統合前の分割領域の数と、領域集合の数との合計が３以上でない場合には（ステップＳ２１９，Ｎｏ）、ステップＳ２１８に移行する。

一方、特定部２５６は、統合前の分割領域の数と、領域集合の数との合計が３以上の場合には（ステップＳ２１９，Ｙｅｓ）、分割領域を統合する（ステップＳ２２０）。ステップＳ２２０において、特定部２５６は、遅延量が最小となる分割領域の組合せを、１つに統合して、領域集合に設定する。特定部２５６は、脈波情報２４５を更新し（ステップＳ２２１）、ステップＳ２１３に移行する。ステップＳ２２１において、例えば、特定部２５６は、領域集合に含まれる各脈波信号を平均化した脈波信号を、係る領域集合の脈波信号に設定することで、脈波情報２４５を更新する。

次に、本実施例２に係る情報処理装置２００が長期計測を行う場合の処理手順について説明する。図１７は、本実施例２に係る情報処理装置が長期計測を行う場合の処理手順を示すフローチャートである。図１７に示すように、情報処理装置２００の受付部２５１は、タイマを設定する（ステップＳ２５１）。ステップＳ２５１において、例えば、受付部２５１は、タイマ長Ｔ１、タイマｔ１＝０をそれぞれ設定する。

情報処理装置２００の撮像機能部２１０は、生体画像を撮影し、受付部２５１に出力する（ステップＳ２５２）。ステップＳ２５２において、受付部２５１は、生体画像情報を、画像バッファ２４１に蓄積する。情報処理装置２００の特徴点座標決定部２５２は、特徴点の座標を決定する（ステップＳ２５３）。情報処理装置２００の評価領域設定部２５３は、評価領域を決定する（ステップＳ２５４）。情報処理装置２００の分割部２５４は、評価領域を分割し、情報処理装置２００の脈波信号検出部２５５は、評価領域情報２４３を基にして、評価対象の領域集合の組を決定する（ステップＳ２５５）。

脈波信号検出部２５５は、評価対象の領域集合の組の輝度を算出する（ステップＳ２５６）。脈波信号検出部２５５は、評価対象の領域集合の組の脈波信号を特定し、特定した脈波信号を脈波情報２４５に格納する（ステップＳ２５７）。脈波信号検出部２５５は、ｔ１＝ｔ１＋１により、ｔ１の値を更新する（ステップＳ２５８）。

脈波信号検出部２５５は、ｔ１の値とＴ１の値とが等しいか否かを判定する（ステップＳ２５９）。脈波信号検出部２５５は、ｔ１の値とＴ１の値とが等しくない場合には（ステップＳ２５９，Ｎｏ）、ステップＳ２５１に移行する。

一方、ｔ１の値とＴ１の値とが等しい場合には（ステップＳ２５９，Ｙｅｓ）、情報処理装置２００の特定部２５６は、評価対象の領域集合の脈波信号を比較して、遅延量を算出する（ステップＳ２６０）。情報処理装置２００のＰＷＶ算出部２５７は、遅延量と脈波伝搬距離情報２４４とを基にして、ＰＷＶを算出する（ステップＳ２６１）。ＰＷＶ算出部２５７は、出力表示処理を行う（ステップＳ２６２）。出願表示処理は、図１１に示した出願表示処理と同様である。

次に、本実施例２に係る情報処理装置２００の効果について説明する。情報処理装置２００は、各分割領域の脈波信号をそれぞれ比較して遅延量が閾値未満となる分割領域の組を統合し、統合した分割領域間の遅延量、または、統合した分割領域と結合されていない分割領域との遅延量をそれぞれ測定し、測定した各遅延量のうち最大の遅延量を特定する。このように遅延量が類似する分割領域を統合することで、輝度のＳＮ比（signal-noise ratio）を改善することができ、精度よく脈波信号を検出することができる。

また、情報処理装置２００は、各分割領域の脈波信号をそれぞれ比較して遅延量が類似する分割領域を統合する処理を、統合後の領域の数が３つ以上である場合のみ実行し、統合された領域の脈波信号の遅延量を特定する。これによって、遅延量が小さい評価領域であっても、遅延量の最も大きな２つの統合された領域が残り、精度よく、ＰＷＶを算出することができる。

次に、本実施例３に係る情報処理装置の構成の一例について説明する。図１８は、本実施例３に係る情報処理装置の構成を示す機能ブロック図である。図１８に示すように、この情報処理装置３００は、撮像機能部３１０と、入力部３２０と、表示部３３０と、通信部３３５と、記憶部３４０と、制御部３５０とを有する。図１８において、撮像機能部３１０、入力部３２０、表示部３３０、通信部３３５に関する説明は、実施例１に示した撮像機能部１１０、入力部１１０、表示部１３０、通信部１３５と同様であるため、説明を省略する。

記憶部３４０は、画像バッファ３４１、生体特徴情報３４２、評価領域情報３４３、脈波伝搬距離情報３４４、脈波情報３４５を有する。記憶部３４０は、例えば、ＲＡＭ、ＲＯＭ、フラッシュメモリなどの半導体メモリ素子などの記憶装置に対応する。画像バッファ３４１、生体特徴情報３４２、評価領域情報３４３、脈波伝搬距離情報３４４、脈波情報３４５に関する説明は、図１に示した画像バッファ１４１、生体特徴情報１４２、評価領域情報１４３、脈波伝搬距離情報１４４、脈波情報１４５に関する説明と同様である。

制御部３５０は、受付部３５１、特徴点座標決定部３５２、評価領域設定部３５３、分割部３５４、脈波信号検出部３５５、特定部３５６、ＰＷＶ算出部３５７、除去部３５８を有する。制御部３５０は、例えば、ＡＳＩＣや、ＦＰＧＡなどの集積装置に対応する。また、制御部３５０は、例えば、ＣＰＵやＭＰＵ等の電子回路に対応する。このうち、受付部３５１、評価領域設定部３５３、分割部３５４、脈波信号検出部３５５に関する説明は、図１に示した受付部１５１、評価領域設定部１５３、分割部１５４、脈波信号検出部１５５に関する説明と同様である。また、特定部３５６およびＰＷＶ算出部３５７に関する説明は、図１に示した特定部１５６およびＰＷＶ算出部１５７に関する説明と同様である。

特徴点座標決定部３５２は、実施例１に示した特徴点座標決定部１５２と同様にして、生体特徴情報３４２を基にして、画像バッファ３４１に蓄積された生体画像情報から特徴点の座標を決定する処理部である。特徴点座標決定部３５２は、生体画像情報と、決定した特徴物と、特徴点の座標とを対応付けて、評価領域設定部３５３、除去部３５８に出力する。

除去部３５８は、時間変化に伴う特徴点の移動量の変化を示す移動量信号を検出し、移動量信号の周波数成分を特定する。除去部３５８は、脈波情報３４５に記憶された脈波信号から、移動量信号の周波数成分を除去する。以下において、除去部３５８の処理の一例について説明する。

除去部３５８は、顔の中心座標の移動量を、式（４）に基づいて算出する。顔の中心座標は、例えば、鼻の特徴点の座標に対応する。式（４）において、ｘ_ｔ、ｙ_ｔは、現フレームの中心座標を示す。ｘ_ｔ−ｔ’、ｙ_ｔ−ｔ’は、所定フレーム前の中心座標を示す。除去部３５８は、フレーム毎に上記処理を繰り返し実行することで、移動量信号を算出する。

除去部３５８は、周知技術により、移動量信号を周波数解析し、周波数毎のスペクトル密度を算出する。除去部３５８は、スペクトル密度が閾値以上となる周波数を特定する。除去部３５８は、特定した周波数を中心とする帯域除外フィルタを脈波情報３４５の脈波信号にかけることで、脈波情報３４５に記憶された脈波信号から、移動量信号の周波数成分を除去する。

図１９は、除去部の処理を説明するための図である。図１９の信号６０は、脈波信号を示す。図１９の横軸は周波数に対応し、縦軸はスペクトル密度に対応する。例えば、除去部３５８は、上記処理により、スペクトル密度が閾値以上となる周波数を２Ｈｚとして特定した場合には、係る２Ｈｚの周波数を中心とする帯域６５をカットするフィルタを脈波信号６０にかける。

次に、本実施例３に係る情報処理装置３００の処理手順について説明する。ここで、初期学習に関する処理手順は、実施例１に示したものと同様であるため、説明を省略する。図２０は、本実施例３に係る情報処理装置が長期計測を行う場合の処理手順を示すフローチャートである。図２０に示すように、情報処理装置３００の受付部３５１は、タイマを設定する（ステップＳ３０１）。ステップ３５１において、例えば、受付部３５１は、タイマ長Ｔ１、タイマｔ１＝０をそれぞれ設定する。

情報処理装置３００の撮像機能部３１０は、生体画像を撮影し、受付部３５１に出力する（ステップＳ３０２）。ステップＳ３０２において、受付部３５１は、生体画像情報を、画像バッファ３４１に蓄積する。情報処理装置３００の特徴点座標決定部３５２は、特徴点の座標を決定する（ステップＳ３０３）。情報処理装置３００の評価領域設定部３５３は、評価領域を決定する（ステップＳ３０４）。情報処理装置３００の分割部３５４は、評価領域を分割し、情報処理装置３００の脈波信号検出部３５５は、評価領域情報３４３を基にして、評価対象の分割領域の組を決定する（ステップＳ３０５）。

情報処理装置３００の除去部３５８は、時間変化に伴う特徴点の移動量の変化を基にして、移動量信号を検出する（ステップＳ３０６）。脈波信号検出部３５５は、評価対象の分割領域の組の輝度を算出する（ステップＳ３０７）。脈波信号検出部３５５は、評価対象の分割領域の組の脈波信号を特定し、特定した脈波信号を脈波情報３４５に格納する（ステップＳ３０８）。脈波信号検出部３５５は、ｔ１＝ｔ１＋１により、ｔ１の値を更新する（ステップＳ３０９）。

脈波信号検出部３５５は、ｔ１の値とＴ１の値とが等しいか否かを判定する（ステップＳ３１０）。脈波信号検出部３５５は、ｔ１の値とＴ１の値とが等しくない場合に（ステップＳ３１０，Ｎｏ）、ステップＳ３０１に移行する。

一方、ｔ１の値とＴ１の値とが等しい場合には（ステップＳ３１０，Ｙｅｓ）、除去部３５８は、体動除去処理を実行する（ステップＳ３１１）。ステップＳ３１１において、例えば、除去部３５８は、移動量信号を周波数解析し、周波数毎のスペクトル密度を算出する。除去部３５８は、スペクトル密度が閾値以上となる周波数を特定する。除去部３５８は、特定した周波数を中心とする帯域除外フィルタを脈波情報３４５の脈波信号にかけることで、脈波情報３４５に記憶された脈波信号から、移動量信号の周波数成分を除去する。

情報処理装置３００の特定部３５６は、評価対象の分割領域の脈波信号を比較して、遅延量を算出する（ステップＳ３１２）。情報処理装置３００のＰＷＶ算出部３５７は、遅延量と脈波伝搬距離情報３４４とを基にして、ＰＷＶを算出する（ステップＳ３１３）。ＰＷＶ算出部３５７は、出力表示処理を行う（ステップＳ３１４）。ステップＳ３１４における出力表示処理の処理手順は、図１１に示した処理手順と同様である。

次に、本実施例３に係る情報処理装置３００の効果について説明する。情報処理装置３００は、生体画像情報から生体の特徴点の位置を特定し、時間変化に伴う特徴点の移動量の変化を示す移動量信号の周波数成分を特定し、脈波信号の周波数成分から移動量信号の周波数成分を除去する。このため、ユーザの移動に起因する影響を脈波信号から取り除くことができ、ユーザのＰＷＶの精度を更に向上させることができる。

次に、上記実施例に示した情報処理装置１００，２００，３００と同様の機能を実現するデータ配置プログラムを実行するコンピュータの一例について説明する。図２１は、脈波計測プログラムを実行するコンピュータの一例を示す図である。

図２１に示すように、コンピュータ４００は、各種演算処理を実行するＣＰＵ４０１と、ユーザからのデータの入力を受け付ける入力装置４０２と、ディスプレイ４０３とを有する。また、コンピュータ４００は、記憶媒体からプログラム等を読取る読み取り装置４０４と、ネットワークを介して他のコンピュータとの間でデータの授受を行うインターフェース装置４０５と、カメラ４０６とを有する。また、コンピュータ４００は、各種情報を一時記憶するＲＡＭ４０７と、ハードディスク装置４０８とを有する。そして、各装置４０１〜４０８は、バス４０９に接続される。

ハードディスク装置４０８は、分割プログラム４０８ａ、検出プログラム４０８ｂ、特定プログラム４０８ｃを有する。ＣＰＵ４０１は、分割プログラム４０８ａ、検出プログラム４０８ｂ、特定プログラム４０８ｃを読み出して、ＲＡＭ４０７に展開する。分割プログラム４０８ａは、分割プロセス４０７ａとして機能する。検出プログラム４０８ｂは、検出プロセス４０７ｂとして機能する。特定プログラム４０８ｃは、特定プロセス４０７ｃとして機能する。例えば、分割プロセス４０７ａは、分割部１５４，２５４，３５４に対応する。検出プロセス４０７ｂは、脈波信号検出部１５５，２５５，３５５に対応する。特定プロセス４０７ｃは、特定部１５６，２５６，３５６に対応する。

なお、分割プログラム４０８ａ、検出プログラム４０８ｂ、特定プログラム４０８ｃについては、必ずしも最初からハードディスク装置４０８に記憶させておかなくても良い。例えば、コンピュータ４００に挿入されるフレキシブルディスク（ＦＤ）、ＣＤ−ＲＯＭ、ＤＶＤディスク、光磁気ディスク、ＩＣカードなどの「可搬用の物理媒体」に各プログラムを記憶させておく。そして、コンピュータ４００が各プログラム４０８ａ〜４０８ｃを読み出して実行するようにしてもよい。

以上の各実施例を含む実施形態に関し、さらに以下の付記を開示する。

（付記１）生体の表面が撮影された画像を複数の領域に分割する分割部と、
複数の領域毎に、領域に含まれる画素の時間変化に伴う輝度変化に基づいて脈波信号を検出する検出部と、
各領域の脈波信号をそれぞれ比較して各領域の脈波信号の遅延量をそれぞれ測定し、測定した各遅延量のうち最大の遅延量を特定する特定部と
を有することを特徴とする情報処理装置。

（付記２）前記特定部は、各領域の脈波信号をそれぞれ比較して遅延量が閾値未満となる領域の組を統合し、統合した領域間の遅延量、または、統合した領域と統合されていない領域との遅延量をそれぞれ測定し、測定した各遅延量のうち最大の遅延量を特定する処理を更に実行することを特徴とする付記１に記載の情報処理装置。

（付記３）前記特定部は、各領域の脈波信号をそれぞれ比較して遅延量が類似する領域を統合する処理を、統合後の領域の数が３つ以上である場合のみ繰り返し実行し、統合された領域の脈波信号の遅延量を特定する処理を更に実行することを特徴とする付記１に記載の情報処理装置。

（付記４）前記画像から前記生体の特徴点の位置を特定し、時間変化に伴う前記特徴点の移動量の変化を示す信号の周波数成分を特定し、前記脈波信号の周波数成分から前記脈波信号の周波数成分を除去する除去部を更に有することを特徴とする付記１、２または３に記載の情報処理装置。

（付記５）前記特定部が特定した遅延量を基にして、脈波伝搬速度を決定する決定部を更に有することを特徴とする付記１〜４の何れか一つに記載の情報処理装置。

（付記６）コンピュータに、
生体の表面が撮影された画像を複数の領域に分割し、
複数の領域毎に、領域に含まれる画素の時間変化に伴う輝度変化に基づいて脈波信号を検出し、
各領域の脈波信号をそれぞれ比較して各領域の脈波信号の遅延量をそれぞれ測定し、測定した各遅延量のうち最大の遅延量を特定する
処理を実行させることを特徴とする脈波計測プログラム。

（付記７）前記特定する処理は、各領域の脈波信号をそれぞれ比較して遅延量が閾値未満となる領域の組を統合し、統合した領域間の遅延量、または、統合した領域と統合されていない領域との遅延量をそれぞれ測定し、測定した各遅延量のうち最大の遅延量を特定する処理を更に実行することを特徴とする付記６に記載の脈波計測プログラム。

（付記８）前記特定する処理は、各領域の脈波信号をそれぞれ比較して遅延量が類似する領域を統合する処理を、統合後の領域の数が３つ以上の場合のみ繰り返し実行し、統合された領域の脈波信号の遅延量を特定する処理を更に実行することを特徴とする付記６に記載の脈波計測プログラム。

（付記９）前記画像から前記生体の特徴点の位置を特定し、時間変化に伴う前記特徴点の移動量の変化を示す信号の周波数成分を特定し、前記脈波信号の周波数成分から前記脈波信号の周波数成分を除去する処理を更にコンピュータに実行させることを特徴とする付記６、７または８に記載の脈波計測プログラム。

（付記１０）前記遅延量を基にして、脈波伝搬速度を決定する処理を更にコンピュータに実行させることを特徴とする付記６〜９の何れか一つに記載の脈波計測プログラム。

（付記１１）コンピュータが実行する脈波計測方法であって、
生体の表面が撮影された画像を複数の領域に分割し、
複数の領域毎に、領域に含まれる画素の時間変化に伴う輝度変化に基づいて脈波信号を検出し、
各領域の脈波信号をそれぞれ比較して各領域の脈波信号の遅延量をそれぞれ測定し、測定した各遅延量のうち最大の遅延量を特定する
処理を実行することを特徴とする脈波計測方法。

（付記１２）前記特定する処理は、各領域の脈波信号をそれぞれ比較して遅延量が閾値未満となる領域の組を統合し、統合した領域間の遅延量、または、統合した領域と統合されていない領域との遅延量をそれぞれ測定し、測定した各遅延量のうち最大の遅延量を特定する処理を更に実行することを特徴とする付記１１に記載の脈波計測方法。

（付記１３）前記特定する処理は、各領域の脈波信号をそれぞれ比較して遅延量が類似する領域を統合する処理を、統合後の領域の数が３つ以上の場合のみ繰り返し実行し、統合された領域の脈波信号の遅延量を特定する処理を更に実行することを特徴とする付記１１に記載の脈波計測方法。

（付記１４）前記画像から前記生体の特徴点の位置を特定し、時間変化に伴う前記特徴点の移動量の変化を示す信号の周波数成分を特定し、前記脈波信号の周波数成分から前記脈波信号の周波数成分を除去する処理を更にコンピュータに実行させることを特徴とする付記１１、１２または１３に記載の脈波計測方法。

（付記１５）前記遅延量を基にして、脈波伝搬速度を決定する処理を更にコンピュータに実行させることを特徴とする付記１１〜１４の何れか一つに記載の脈波計測方法。

１００，２００，３００ 情報処理装置
１５４，２５４，３５４ 分割部
１５５，２５５，３５５ 脈波信号検出部
１５６，２５６，３５６ 特定部

Claims (7)

1. 生体の表面が撮影された画像を複数の領域に分割する分割部と、
   複数の領域毎に、領域に含まれる画素の時間変化に伴う輝度変化に基づいて脈波信号を検出する検出部と、
   各領域の脈波信号をそれぞれ比較して各領域の脈波信号の遅延量をそれぞれ測定し、測定した各遅延量のうち最大の遅延量を特定する特定部と
   を有することを特徴とする情報処理装置。
2. 前記特定部は、各領域の脈波信号をそれぞれ比較して遅延量が閾値未満となる領域の組を統合し、統合した領域間の遅延量、または、統合した領域と統合されていない領域との遅延量をそれぞれ測定し、測定した各遅延量のうち最大の遅延量を特定する処理を更に実行することを特徴とする請求項１に記載の情報処理装置。
3. 前記特定部は、各領域の脈波信号をそれぞれ比較して遅延量が類似する領域を統合する処理を、統合後の領域の数が３つ以上の場合のみ繰り返し実行し、統合された領域の脈波信号の遅延量を特定する処理を更に実行することを特徴とする請求項１に記載の情報処理装置。
4. 前記画像から前記生体の特徴点の位置を特定し、時間変化に伴う前記特徴点の移動量の変化を示す信号の周波数成分を特定し、前記脈波信号の周波数成分から前記特徴点の移動量の変化を示す信号の周波数成分を除去する除去部を更に有することを特徴とする請求項１、２または３に記載の情報処理装置。
5. 前記特定部が特定した遅延量を基にして、脈波伝搬速度を決定する決定部を更に有することを特徴とする請求項１〜４の何れか一つに記載の情報処理装置。
6. コンピュータに、
   生体の表面が撮影された画像を複数の領域に分割し、
   複数の領域毎に、領域に含まれる画素の時間変化に伴う輝度変化に基づいて脈波信号を検出し、
   各領域の脈波信号をそれぞれ比較して各領域の脈波信号の遅延量をそれぞれ測定し、測定した各遅延量のうち最大の遅延量を特定する
   処理を実行させることを特徴とする脈波計測プログラム。
7. コンピュータが実行する脈波計測方法であって、
   生体の表面が撮影された画像を複数の領域に分割し、
   複数の領域毎に、領域に含まれる画素の時間変化に伴う輝度変化に基づいて脈波信号を検出し、
   各領域の脈波信号をそれぞれ比較して各領域の脈波信号の遅延量をそれぞれ測定し、測定した各遅延量のうち最大の遅延量を特定する
   処理を実行することを特徴とする脈波計測方法。

Images