JP5834011B2 - 生体中の少なくとも周期的現象を表す成分を含む信号を処理する方法およびシステム - Google Patents

Description

本発明は、生体中の周期的現象を表す成分を少なくとも含む信号を処理する方法およびシステムに関する。本発明はコンピュータ・プログラムにも関する。

US5,431,170は、ヘッドバンドを装着するとセンサー・ユニットが額にあてがわれるようにヘッドバンドに取り付けられたセンサー・ユニットを有する脈拍数計を開示している。センサー・ユニットは、二つの発光器と、額の身体組織内での反射後にそれぞれの発光器からの光を受光する光センサーとを有する。一つの発光器は赤外光を放出し、この発光器からの光を受光するセンサーによって与えられる電気信号は、額における血管を通る血流の脈動に従って変動する。他方の発光器は異なる、好ましくは実質的により短い波長の光、たとえば黄色の光を放出し、そのセンサーからの信号は流れの脈動に起因する変動とは比較的独立になるようにする。しかしながら、いずれの信号も、動きに起因する額の振動に従って変動する。ノイズ打ち消し回路がそれぞれのセンサーからの出力信号の周波数解析を実行してもよい：二つの信号における対応する周波数または周波数範囲におけるピークがノイズと見なされ、赤外センサーからの信号における残りの周波数が脈拍数を表す。

この既知のシステムの問題は、比較的邪魔になり、一つの特定の身体部分、つまりセンサー・ユニットが取り付けられている身体部分の動きしか補償できないということである。たとえば腕の動きは血流に影響するであろうが、前記の、より短い波長で動作する発光器および受光器の組み合わせによっては拾われない。一方、そのような動きは、呼吸や心拍に起因する血流変動を検出するのを一層難しくする。

US5,431,170

Viola,P. and Jones, M.J., "Robust real-time object detection", Proc. of IEEE workshop on statistical and computational theories of vision, 13 July 2001 De Haan et al., "True-motion estimation with 3-D recursive search block matching", IEEE transactions on circuits and systems for video technology, 3(5), October 1993, pp.368-379 B. Widrow and S.D. Stearns, Adaptive Signal Processing, Prentice-Hall, Inc., Englewood Cliffs, New Jersey 07632, 1985

生体中の少なくとも周期的現象を表す成分を含む比較的クリーンな信号を、該周期的現象を定量化する解析のために提供できるようにするとともに、幅広い型の動きについて補償できる、上述した型の方法、システムおよびコンピュータ・プログラムを提供することが望ましい。

この目的のため、本発明のある側面によれば、生体中の少なくとも周期的現象を表す成分を含む信号を処理する方法であって：
前記生体を示す画像のシーケンスを取得する段階と；
リモート光体積変動記録法を使って前記画像のシーケンスから生体中の少なくとも周期的現象を表す成分を有する少なくとも一つの第一の信号を抽出する段階と；
前記画像のシーケンスのビデオ・ベースの動き解析を実行することによって、前記生体の動きの少なくとも周期的成分を表す別個のデータを取得する段階と；
前記データを使って、少なくとも、前記生体の動きの前記周期的成分に対応する前記第一の信号の成分を抑制する段階とを含む方法が提供される。

前記少なくとも一つの第一の信号と同期して、前記生体を示す画像のシーケンスを取得して、前記画像のシーケンスのビデオ・ベースの動き解析を実行することによって、幅広い型の動きが、邪魔にならず、比較的効率的な仕方で定量化できる。生体の大きな部分が前記画像のシーケンスにおいて表現される限り、局所化された動きでさえ検出できる（たとえば、腕または脚のみの動き）。これは、多数のセンサーではなく、単一のビデオ・カメラで行うことができる。身体に取り付けられるセンサーも回避される。他の型のリモート・センシングに比べ、多くの異なる方向における動きが、特に画像平面に平行な動きのほかカメラに向かう動きおよびカメラから離れる動きが検出されることができる。前記生体の動きの少なくとも周期的成分を表す別個のデータを取得することは、画像または該画像中の関心領域（ROI: Region Of Interest）の整列に関わる画像解析を使って、前記第一の信号が前記画像のシーケンスから抽出される場合にも有益となる。そのような整列は、動きに起因して間接的に導入される画像アーチファクト、たとえば影に対して常に堅牢ではない。よって、そのような状況にあっても、前記生体の動きの少なくとも周期的成分を表す別個のデータを取得し、前記データを使って、少なくとも、前記生体の動きの前記周期的成分に対応する前記第一の信号の成分を抑制することは、前記第一の信号の品質を改善するであろう。

本発明の好ましい実施形態は従属請求項において定義される。特許請求されるシステムおよび特許請求されるコンピュータ・プログラムは特許請求される方法および従属請求項の記載と同様および／または同一の好ましい実施形態をもつことは理解しておかねばならない。

本発明によれば、前記少なくとも一つの第一の信号を取得する段階は、少なくとも一つの第一の信号を前記画像のシーケンスから抽出することを含む。

これは、一連のデジタル画像を撮影するカメラの形の一つのセンサーのみを使って実装できるので、比較的効率的である。その場合、前記第一の信号を動きの少なくとも周期的成分を表す前記別個のデータと同期させる問題は自明になる。よって、前記データを使って、少なくとも、前記生体の動きの前記周期的成分に対応する前記第一の信号の成分を抑制する段階も実装がより簡単になる。

ある変形では、前記少なくとも一つの第一の信号を前記画像のシーケンスから抽出することは、それぞれの画像からのピクセル値の組み合わせに対応する変数の値を得ることを含む。

ピクセル値を組み合わせることにより、色および／または輝度の非常に小さな変動が検出できる。一般に合計演算、可能性としては重み付けされた和を含むであろう組み合わせる段階によってランダム・ノイズは除去される。こうして、この実施形態では、リモート光体積変動記録法（photoplethysmography）が提供できる。

あるさらなる変形では、前記少なくとも一つの第一の信号をシーケンスから抽出することは、前記シーケンスにわたる前記変数の変動のみを決定することを含む。

この変形は、周期的な生物学的現象を特徴付けるときには、輝度または色の変動のみが関心対象であるとの認識に基づく。

本発明によれば、ビデオ・カメラを使った邪魔にならない生命徴候（vital sign）（たとえば心拍数／呼吸数）モニタリング、あるいはリモート光体積変動記録法が可能にされる。ある実施形態では、アルゴリズムが、血液体積および血液酸素供給とともに変動する、ある人の平均肌色を追跡する。しかしながら、この方法は、動き、肌色の局所的変動および照明レベル／スペクトル変化に敏感であり、そのため商業的に有望な用途にはそれほど好適ではない。動きへの敏感さは、本邪魔にならない生命徴候モニタリングを、個人が自由に動きたい用途には特に不適にする。そのような場合のいい例が、フィットネス用途における心拍モニタリングである。

よって、好ましい実施形態によれば、平均肌色信号における、動きに起因する望ましくない成分を抑制することが提案される。これは、動きが周期的であり、動きの周波数が関心のある周波数範囲、たとえば40拍毎分（BPM: beats per minute）から210BPMの間にはいる場合に特に重要である。この場合、心拍信号を動き信号から区別することは難しく、この望ましくない成分の除去は安定した心拍数推定のためには非常に重要である。しかしながら、本方法は周期的な動きのみに制約されない。

したがって、本方法のある実施形態は、画像において表現されている少なくとも一つの被験体について少なくとも一つの動きベクトルを計算し、動きベクトル値の反復性（repetitiveness）を解析することを含む。

この実施形態は、ある種の型の、たとえばある大きさの、周期的動きのみに関する情報を単離することを許容する。それはまた、画像のシーケンスにおいて表現されている複数の生体の間の区別をするためにも有用である。

ある変形は、画像において表現されているある被験体の複数の部分のそれぞれについて動きベクトル値を計算することを含む。

この変形は、生体中の少なくとも周期的現象を表す成分を有する複数の信号が取得される場合に有用である。それら各信号について適切な訂正信号が形成できるからである。さらに、前記複数の部分の動きがある量以上に異なる場合には、それらの信号のそれぞれを、生体中の周期的現象を表す一つの信号に組み合わせる前に別個に訂正することが可能である。他方、差がより小さい場合には、生体中の少なくとも周期的現象を表す成分を有する複数の信号が組み合わされることができ、該組み合わせが、前記生体の動きの前記周期的成分に対応する成分を抑制することによって訂正されることができる。これは本方法の比較的効率的な実施形態である。

ある実施形態では、動きの少なくとも周期的成分を表す別個のデータを取得する段階は、前記生体の少なくとも一部の空間的な変位を定量化することを含む。

この実施形態は、ある種の型の周期的動き、たとえば大規模な変位のみを表す前記第一の信号の成分を抑制することが可能であるという効果をもつ。

提案される方法のさらなる実施形態は、適応エコー打ち消し（adaptive echo cancellation）に基づく。時間において関心のある皮膚領域の動きを記述する動きベクトルが非常にクリーンな信号であることが見出されている。さらに、動きベクトルはいかなる心拍信号も含まず、よって適応的フィルタのための参照信号として使うのに非常に好適である。一般に、動きアーチファクトを著しく軽減するために非常に少数の係数しか要求されないことがわかる。信号から動きアーチファクトを除去する直接的な結果は、ずっとクリーンな心拍信号を得ることができるということであり、それは、ピーク・ツー・ピーク検出によって時間領域で代替的な方法で心拍数の推定を可能にする。クリーンな信号の追加的な利点は、HRV（heart rate variability［心拍数変動度］）が推定できるということである。

したがって、本方法のある実施形態は、画像において表現されている少なくとも一つの被験体について少なくとも一つの動きベクトルを計算する段階と、前記少なくとも一つの動きベクトルをフィルタ処理する段階と、前記第一の信号から前記フィルタ処理された動きベクトルを減算して、前記第一の信号の、前記生体の動きの前記周期的成分に対応する成分を抑制し、擾乱が軽減された第一の信号を得る段階とを含む。好ましくは、前記フィルタ処理は、前記フィルタ処理の係数が適応的に更新される適応的フィルタ処理を含む。

好ましくは、規格化（normalization）、特に最小平均平方（least-mean-square）（NLMS）アルゴリズムが適用される。NLMSアルゴリズムは、LMSアルゴリズムの特別な場合である。この規格化は動き信号（参照信号）に、よって間接的には係数更新に適用される。さらに、動きベクトルのエネルギーが近似的にわかるシステムについては、必ずしもNLMSアルゴリズムを使う必要はなく、係数を更新する間にこのエネルギーを考慮に入れることが可能である。

前記フィルタ処理の係数は好ましくは、擾乱が軽減された第一の信号（すなわち第一の信号からフィルタ処理された動きを引いたもの）のエネルギーが最小化されるよう、更新される。

ある好ましい実施形態では、特に第一の方向（たとえば水平方向またはx方向）および前記第一の方向に垂直な第二の方向（たとえば垂直方向またはy方向）の、少なくとも二つの動きベクトルが別個に計算され、別個にフィルタ処理される。

本発明のもう一つの側面によれば、生体中の少なくとも周期的現象を表す成分を含む信号を処理するシステムであって：
前記生体を示す画像のシーケンスを取得するインターフェースと；
リモート光体積変動記録法を使って前記画像のシーケンスから生体中の少なくとも周期的現象を表す成分を有する少なくとも一つの第一の信号を抽出する少なくとも一つのコンポーネントと；
前記画像のシーケンスのビデオ・ベースの動き解析を実行する画像データ処理システムを含む、前記生体の動きの少なくとも周期的成分を表す別個のデータを取得する少なくとも一つのコンポーネントと；
前記データを使って、少なくとも、前記生体の動きの前記周期的成分に対応する前記第一の信号の成分を抑制するよう構成された少なくとも一つのデータ処理コンポーネントとを含むシステムが提供される。

ある実施形態では、本システムは、本発明に基づく方法を実行するよう構成される。

本発明のもう一つの側面では、機械可読媒体に組み込まれたときに、情報処理機能をもつシステムに本発明に基づく方法を実行させることのできる一組の命令を含むコンピュータ・プログラムが提供される。

本発明は、付属の図面を参照しつつさらに詳細に説明される。
生体中の周期的現象を表す成分を少なくとも含む時間変動する信号を形成するシステムの概略的なブロック図である。 前記時間変動する信号を形成する方法の第一の一般的な実施形態を示すフローチャートである。 時間変動する信号を形成する方法の第二の実施形態を示すフローチャートである。 画像のシーケンスの限られた区間をカバーする諸画像から得られるピクセル値のサブセットの諸グループを示す図である。 ピクセル値のサブセットの選択されたグループからのある信号セクションの抽出を示す図である。 生体中の周期的現象を表す成分を少なくとも含む信号を処理する代替的なシステムを概略的に示す図である。 動き誘起された擾乱打ち消しなしの、ルームランナー上を走る個人から得られた信号のスペクトログラムを示す図である。 心拍信号における動き誘起された汚染を軽減するための動き補償ユニットのある実施形態のブロック図である。 動き誘起された擾乱打ち消し後の、ルームランナー上を走る個人から得られた信号のスペクトログラムを示す図である。

リモート光体積変動記録法〔フォトプレチスモグラフィー〕（PPG: photoplethysmography）を実行するのに好適なシステムのいくつかの実施形態をここで述べる。光体積変動記録法は、皮膚反射率変動を使ってある種の周期的な生理現象を特徴付ける方法である。人間の皮膚は、少なくとも二つの層をもつオブジェクトとしてモデル化できる。二つの層の一つは表皮（薄い表面の層）、他方は真皮（表皮の下にあるより厚い層）である。入射光線の約5%は表皮において反射され、これはあらゆる波長および皮膚の色について成り立つ。残りの光は、体反射（body reflectance）として知られる現象（二色性反射モデル（Dichromatic Reflection Model）において記述される）においてそれら二つの皮膚の層内で散乱され、吸収される。表皮は光フィルタのように振る舞い、主として光を吸収する。光はその波長および表皮におけるメラニン濃度に依存して透過される。真皮では、光は散乱および吸収の両方を受ける。吸収は血液組成、すなわち血液の内容およびヘモグロビン、ビリルビンおよびベータカロテンのようなその構成成分に依存し、よって、吸収は血流変動に敏感である。真皮の光学的属性は一般にすべての人種について同じである。真皮は、成人の全血管網の約10%をなす、血管の密なネットワークを含む。これらの血管が、体内の血流に従って収縮する。その結果、これらの血管は真皮の構造を変え、それが皮膚層の反射率に影響する。結果として、皮膚反射率変動から心拍数が決定できる。

記載されるシステムおよび方法は光体積変動記録法に限定されるものではなく、関心対象の周期現象を表す信号が捕捉され、動きについて補正される他の状況においても使用できる。光体積変動記録法の場合、本稿に記載される方法は、関心対象の周期現象（光体積変動記録法の場合、脈動する血流）に無関係な動きアーチファクトや照明の変化に対して堅牢である。

本システムはいくつものコンポーネントを有する。そのうち若干のみが図１に示されている。本システムは、データ処理システム１を含む。これは汎用コンピュータまたは特定用途のための特製コンピューティング装置であってもよく、データ処理ユニット２およびメイン・メモリ３を有する。データ処理システム１はさらに、本稿で概説される方法を実行するためのソフトウェアを含むソフトウェアを記憶する大容量記憶装置４を含む。データ処理システム１は、デジタル画像のシーケンスを取得するために、ビデオ・カメラ６へのインターフェース５を有する。データ処理システム１は、出力を知覚できる形に表現するための出力装置８、たとえば視覚表示ユニットへのインターフェース７も有する。もちろん、データ処理システム１はさらに、ユーザー入力装置、ネットワーク・インターフェースなどといった、図１に別個には示さない一般的な機能要素を含む。

ある実施形態（図示せず）では、データ処理システム１はさらに、少なくとも一つの光源へのインターフェースを含む。該光源は、電磁放射のスペクトルの可視部分または近赤外部分の特定の範囲に調整された光源であることができる。

データ処理システム１は、ビデオ・カメラ６によって捕捉された画像のシーケンスの各画像からのピクセル値に基づくある値の少なくとも変動を表す第一の信号を生成するよう構成される。データ処理システム１はさらに、捕捉された画像の別個の解析によって決定される動きを抑制または除去するために、前記第一の信号を処理するよう構成される。

第一の実施形態（図２）では、データ処理システム１によって実行される方法は、画像のシーケンスを取得することを含む（ステップ９）。それらの画像の少なくとも一つはセグメント化され（ステップ１０）、一つまたは複数の関心領域が選択される（ステップ１１）。

セグメント化および選択ステップ１０、１１は、ある種の型の表面を表す、特にある基準を満たす領域を同定するために実行される。ある実施形態では、画像において表現されている人物の顔の一部または全部を表す関心領域が同定される。この実施形態を実装するための好適なアルゴリズムは非特許文献１に記載されている。ある種の形、色および／またはパターンをもつ画像セグメントを認識するための他の好適なアルゴリズムが知られており、このアルゴリズムの代わりに、またはこのアルゴリズムと組み合わせて使用できる。

本記載は、一つの関心領域が選択されるという想定に基づいて進められる。この関心領域が画像シーケンスを通じて追跡される。これを行うための好適なアルゴリズムが非特許文献２に記載されている。

シーケンス中の最初の一対の画像については、データ処理システム１は、関心領域内から、少なくとも一つのピクセル値の諸サブセットを選択することに進む。特に、第一の画像において関心領域から複数のサブセットが選択される（ステップ１２）。選択された各サブセットについて、サブシーケンス中の隣接する画像（一般には現在画像の次の画像だが、直前の画像であってもよい）において、対応するサブセットが見出される（ステップ１３）。

一般に、各サブセットは複数のピクセル値をもち、さらなる画像中で対応するサブセットを見出すステップ１３は何らかの基準に従って類似したサブセットを探すことを含む。これはすなわち、次の画像における関心領域内の最良一致〔ベストマッチ〕するサブセットが見出されるということである。第一の画像からの選択された各サブセットおよび次の画像からの最良一致するサブセットはサブセットのグループをなす。第一の画像からのサブセットはみな互いに異なるので、諸グループの任意の対について、第一の画像および次の画像両方からのサブセットは異なっている。

次のフレーム中の対応するサブセットを見出すことは、第一の画像からの選択された諸サブセットの集合に対応する、諸サブセットの組み合わされた集合を、次の画像において探すグローバル探索を含むことができる。これは低コストの解決策であるが、比較的不正確である可能性が高い。あるいはまた、第一の画像から選択されたピクセル値の各サブセットについて網羅的探索が実行されることができる。これは比較的正確な方法だが、計算的には比較的高価である。もう一つの実施形態では、類似するサブセットは、そのサブセットについての予測される動きベクトルを取得し、予測された動きベクトルを修正して諸候補サブセットを選択し、少なくとも一つの類似性基準と突き合わせて評価することによって見出される。修正は更新ベクトルの形であり、これはランダムであってもよい。そのようなアプローチは、予測探索を表し、これは、現在画像において諸サブセットにわたって数回逐次反復されることができる。初期ステップ１２において選択された各サブセットについて個々に動きベクトルが見出されるという事実のため、あるグループ内のサブセットどうしは比較的よく一致する。予測された動きベクトルの使用は、比較的低い計算コストでの実装を許容する。予測された動きベクトルは、たとえば、現在フレームにおいて選択された諸サブセットの第一のものについて網羅的な探索を実行し、次いで、結果として得られる動きベクトルを他の諸サブセットについての予測された動きベクトルとして使うことによって得ることができる。

次に（ステップ１４）、各グループが解析される。少なくとも一つの所定の基準を満たさないグループは破棄される（ステップ１５）。ある実施形態では、前記少なくとも一つの基準は、最良一致するグループ（単数または複数）が選択されるという意味で相対的な基準である。しかしながら、一般には、前記少なくとも一つの基準は絶対的な基準であろう。形成されるグループの数は非常に多数になり、前記少なくとも一つの基準を満たすものが常に若干はあるだろうからである。

ある個別的実施形態では、前記少なくとも一つの基準は、ピクセル値の各サブセットが最小の空間的一様性をもつべきであるという基準を含む。もう一つの実施形態では、前記少なくとも一つの基準は少なくとも一つの類似性基準を含む。前記類似性基準は、たとえば色および／またはテクスチャーの類似性に関係することができる。ある実施形態では、画像から画像にかけて、接続されたピクセルはある距離より多くは動いてはならないという意味で、動き基準が使われる。これらの基準のそれぞれは、関心対象の現象（光体積変動記録法の場合、周期的な生理現象）に関係ない動きおよび／または照明変化に対する堅牢さを補償するよう設計される。

次いで、サブセットの選択されたグループのみのそれぞれについて、諸信号セグメントが構築される（ステップ１６）。ある実施形態では、あるサブセット内からのピクセル値が、たとえば重み付けられた平均を取ることによって、あるいは平均値を見出すことによって、単一の値に組み合わされる。これは、グループ内の各サブセットについて行われるので、画像毎に一つのそのような値がある。その効果は、ノイズに対する堅牢性が比較的高いということである。ただし、信号値を抽出する他の方法も可能である。ある実施形態では、画像のシーケンスは、画像フレームの二つ以上のシーケンスを含み、各シーケンスが異なるチャネルに対応し、電磁放射のスペクトルの限られた範囲において捕捉された電磁放射の強度を表すピクセル値からなる。たとえば、画像のシーケンスは、複数の異なる色チャネル（たとえば赤、緑および青）における画像フレームのシーケンスを含むことができる。この代替では、組み合わせは、完全な画像をなす各画像フレームからのピクセル値の重み付けされた和を取ることに関わることができる。組み合わせはまた、減算演算に関わることもできる。たとえば、緑の光に対応するチャネルにおけるピクセル値から、すべてのチャネルを横断した平均輝度を減算する。

次いで、あるグループのそれぞれのサブセットに関連付けられた値の間の差が取られる。各グループが二つの画像からのサブセットのみを含む実施形態における信号セグメントは、一つの値のみを含む。選択ステップ１５は、信号セグメントが構築され終わったあとまで遅らされることができることを注意しておく。これは特に、値の組み合わせまたは画像間の該組み合わせの差に基づく選択基準を使うためである。こうして、強い信号値を与えるグループのみが選択されることができる。

ステップ１２〜１６は、次の各画像対について繰り返される。それにより、画像のシーケンスによって定義される時間ベース上の種々の時間区間を表す種々の信号セグメントが得られる。

これらが最終ステップにおいて融合され（ステップ１７）、画像のシーケンスによってカバーされる時間区間に対応する時間区間をカバーする信号が得られる。信号セグメントの融合のためには、さまざまな方法を使用できる。融合は、たとえば、信号間の最高の相関を見出すこと、または外れ値（outliers）の照明（illumination）に基づくことができる。任意的に、諸信号セクションは、融合ステップ１７の結果である時間変動する信号に影響する度合いを決定する異なる重みをもって、融合プロセスに参加することができる。たとえば、重み付け値は、類似したサブセットを見出すステップ１３において、グループ中のサブセットどうしが一致すると決定される信頼度に依存することができる。もう一つの実施形態では、グループ中のサブセット間の空間的距離が少なくとも部分的に重み付けを決定する。距離が大きければ、重みはより低くなる。

図２の方法は、隣り合うサブセットからのサブセットの対のグループの例を使って説明したが、図２の方法は、シーケンスの相続く画像からの三つ以上のサブセットのグループを使うこともできる。それらのグループによって定義される諸区間は重なり合うことができる。たとえば、あるグループは第一から第三の画像のサブセットを含むことができ、別のグループは第二から第四の画像を含むことができる。

さらにもう一つの代替的な実施形態では、ある種のグループは、隣接し合う画像からのサブセットから形成されるのではなく、ある種の画像をスキップすることによって形成される。たとえば、ある実装では、第一の画像から諸サブセットが選択され、第二の画像における最も類似する諸サブセットが見出され、それにより、第一および第二の画像からのサブセットの諸対が得られる。特定の対がある基準を満たさない場合、その対は破棄され、第三の画像からの類似したサブセットが決定される。次いで、サブセットの対を形成するために補間を使うことができ、その対からたとえば信号セクションが抽出される。こうして、この実施形態では、少なくとも一つの基準の集合を満たさないピクセルのサブセットの諸対を破棄するステップ１５に続いて、画像のシーケンスにおいて現在画像からさらに隔たった画像に対して、さらなる諸画像において類似する諸サブセットを見出すステップ１３が繰り返される。

図２は、動きに対する堅牢さを高めるために追加的な補正に含まれるステップを示している。並行ステップ１８において、画像のシーケンスにおいて表される少なくとも一つの被験体の少なくとも一部の動きの指標が決定される。画像のシーケンスにおいて表される被験体の少なくともある種の型の少なくとも周期的な動きを表す別個の信号を得るためである。次いでこの情報は、融合プロセスの結果として得られた信号の、ある成分を抑制する（ステップ１９）ために使われる。関心領域に対応する被験体のみの動きが解析されることは理解されるであろう。

動き解析は、画像のシーケンスにわたる関心領域の動きを表す動きベクトルのシーケンスを決定することによって達成できる。次いで、それらの動きベクトルの反復性を解析することによって、生体の動きの周期的成分を表すデータが決定される。ある実施形態では、振幅のある範囲内の動きのみが考慮に入れられる。

動きベクトルを使うことに対する代替として、ピクセル値またはピクセル値の集合における照明変化が、動きを推定するために使用できる。しかしながら、これはそれほど信頼できず、動きの周波数もあまりよく決定できない。

ある代替的な実施形態では、部分的に重なり合っていてもよい複数の関心領域がある。各関心領域は、画像のシーケンスにおいて表されている単一の被験体の諸部分を表す。各関心領域について、ピクセル値に基づくある値の少なくとも変動を表す信号が、図２に示される仕方で決定される。同様に、各関心領域の動きが別個に解析される。次いで、さまざまな関心領域の動きの相違がある値未満であるかどうかの判定がなされる。そうであれば、ピクセル値に基づくある値の少なくとも変動を表す信号は、さらに一つの信号に融合される。それらの関心領域に関連付けられた諸動き信号に基づく単一の動き信号も形成される。この単一の動き信号は、ピクセル値に基づくある値の少なくとも変動を表す信号をさらに融合することによって得られる前記信号の成分を抑制するために使われる。他方、種々の関心領域に関連付けられた動き信号の相違がある量より大きければ、ピクセル値に基づくある値の少なくとも変動を表す信号は、関連付けられた関心領域に関連付けられた動き信号を使って個々に補正される。次いで、ピクセル値に基づくある値の少なくとも変動を表す補正された諸信号がさらに一つの信号に融合される。

動き抑制は、時間領域または周波数領域において実行できる。ある実施形態では、動き信号と、動きの周期成分に対応する周期成分が抑制されるべき信号との両方のウェーブレット解析が実行される。

融合ステップ１７の結果である信号、動き抑制ステップ１９の結果である信号およびさまざまな区間について得られた（ステップ１６）諸信号セグメントの任意のものが、さらなる処理を受けることができ、一般にはそうされる。そうしたさらなる処理は、関心対象のバイオメトリック信号に対応する情報を抽出するステップ、たとえばアルファ・トレンド平均フィルタ（alpha trend mean filter）を使ったフィルタ処理を含んでいてもよい。さらなる可能なフィルタ処理操作は帯域通過フィルタ処理である。

結果として得られる信号は、たとえば心拍数または呼吸数に対応する、たとえば関心対象の周期的現象の周波数の値を得るために解析されることができる。撮像装置（たとえばCTスキャナ）のゲーティングのために使われることもできる。

画像のシーケンスからのピクセル値に基づくある値の少なくとも変動を表す時間変動する信号を形成する代替的な方法（図３）は、ピクセル値のサブセットのグループをなす二つ以上の画像における類似したサブセットを見出すステップ１３を省く。その代わり、位置において対応する諸サブセットの諸グループが形成される。

この方法は、ある数の画像のバッファを利用する。本方法の各反復工程で、該バッファには最も最近の画像が追加され（ステップ２０）、該バッファからより最近でない画像が除去される。こうして、バッファは画像のシーケンスのサブシーケンスを含む。これらは、画像のセット２２の各画像上にグリッド（図４）をかぶせることによって、空間時間的体積を構成するために使われる（ステップ２１）。セット２２中の画像は、少なくともビデオ・カメラ６によって与えられる画像に基づく。それらの画像は、動き、カメラ動きなどを補償するためにフィルタ処理および／または画像処理操作（図示せず）した結果であってもよい。

グリッドは、幅Δx、高さΔyの空間的開口を画定する。各空間時間的体積は、セット２２中の各画像の諸開口のうち所与の開口中の諸ピクセル値からなる。ある空間的開口内の諸ピクセル値に対応するピクセル値の諸サブセットのシーケンス２３の例が図５に示されている。セット２２内の画像の数が時間区間Δtを決定する。

各空間時間的体積内のピクセル値の諸サブセットは、少なくとも一つの基準からなる集合と突き合わせて評価される（ステップ２４）。前記少なくとも一つの基準を満たすもののみが保持され、生体中の周期現象を表す成分を少なくとも含む信号を取得するために使われる。

一般に、図２の実施形態の選択ステップ１４との関連で述べたのと同じ種類の諸選択基準が、図３の実施形態において使用できる。グリッドが画像全体の上にかぶせられるところでは、図３の方法の選択ステップ２４は、少なくとも、ピクセル値の諸サブセットがある型の表面、たとえば人間の露出した皮膚を表すかどうかを評価するステップをも含むであろう。さらなる選択基準は、ピクセル値の各サブセットの一様性や、諸サブセットのシーケンス２３の時間的一様性に関係する。最も一様な空間時間的諸体積だけ、あるいは何らかの絶対的な一様基準を満たすものだけが保持される。空間領域および時間領域における一様性の計算は、ダイナミックレンジ、分散、エントロピーなどの計算を含むいくつもある方法のいずれを使って行ってもよい。

空間時間的諸体積を構成するステップ２１は、デフォルトの空間的開口および時間区間Δtを使うことができる。ある実施形態では、開口の異なるサイズが、図３の方法の各反復工程で、すなわち、処理される全画像シーケンスの各区間について、使用される。追加的または代替的に、区間の長さ、すなわちセット２２中の画像の枚数が変えられる。開口サイズ（幅Δxまたは高さΔy）は解析される領域の一様性にリンクされることができる。その際、より一様な領域はより大きな開口につながるであろう。他の好適な影響する因子は、ビデオ・カメラ６のセンサーの空間的分解能やそうしたセンサーの他の特徴を含む。他の実施形態では、グリッドは、異なる開口が異なる寸法をもつという意味で不規則である。適切な寸法は、オブジェクト・セグメンテーション、色セグメンテーションおよび背景‐前景セグメンテーションのうちの少なくとも一つを画像のセット２２の少なくとも一つの画像に対して実行することによって決定されることもできる。セグメンテーションが比較的小さな諸セグメントを与える場合には、グリッドの諸開口も小さいべきである。画像のセット２２によってカバーされる区間Δtの適切なサイズは、一般に、各画像中に存在する動きの量に依存する。このように、選択ステップ２４に関係する解析の一部はこれらの実施形態では、空間時間的体積を構成するために実行されていることになる。区間Δtを決定できる他の因子は、ビデオ・カメラ６のフレーム・レート、画像において表現されている被験体の動きの強度および照明の変化を含む。

空間時間的体積のあるものが選択されたのち、選択された各空間時間的体積について、その空間時間的体積のピクセル値に基づくある値の時間的変動が抽出される（ステップ２５）。このステップ２５において、あるサブセット内からのピクセル値は、たとえば重み付き平均を取ることにより、あるいは平均値を見出すことにより、単一の値に組み合わされる。これは、シーケンス２３（図５）中の各サブセットについて行われ、そのため、画像毎にそのような一つの値がある。ピクセル値の組み合わせの効果は、ノイズに対する堅牢性が比較的高いということである。ただし、信号値を抽出する他の方法も可能である。ある実施形態では、各画像は二つ以上の画像フレームを有し、各フレームが異なるチャネルに対応し、電磁放射のスペクトルの限られた範囲において捕捉された電磁放射の強度を表すピクセル値からなる。たとえば、画像のシーケンスは、種々の色チャネル（たとえば赤、緑および青）における画像フレームのシーケンスを有することができる。この代替では、組み合わせは、完全な画像をなす各画像フレームからのピクセル値の重み付けされた和を取ることに関わってもよい。組み合わせはまた、減算演算に関わることもできる。たとえば、緑の光に対応するチャネルにおけるピクセル値から、平均輝度を減算する。

図５に示されるように、ピクセル値の一つのサブセットについて計算された値と、シーケンス２３中のすべてのサブセットについての値の時間平均との間の差が計算される。一般に、ある時間期間にわたる色値の空間平均輝度の時間的ゆらぎを抽出するために、他のいかなる方法を使ってもよい。これらの変動は、振幅変化および周波数値のある範囲内にあるべきである。

図示した実施形態では、空間時間的諸体積を形成するためにスライディング窓が使われる。こうして、各反復工程において、画像の現在セット２２について得られる諸信号セグメントが、互いと、また先行する諸反復工程において得られた重なる諸信号セグメントと融合される（ステップ２６）。融合は、たとえば信号間の最高の相関を見出すことに基づく、または照明外れ値（illuminating outliers）に基づくことができる。これらの信号セグメントは、異なる重み付けパラメータをもって、融合処理に参加できる。これらは、たとえば関連する空間時間的体積の時間区間Δtの長さに依存することができる（より大きな区間Δtをもつ体積により大きな重みが割り当てられる）。

融合された信号は、最も重要な情報を抽出するために使われる。最も重要な情報とは、心拍信号モニタリングの場合には心拍数である。融合された信号から、該融合された信号の、関心のある部分を抽出するために、アルファ・トレンド平均フィルタを使うことができる。

図のように、図２に関係して述べた周期的動きを表す成分を抑制するのと同様の方法が使われる。すなわち、処理される画像のシーケンスにおいて表される生体の動きの周期的成分を少なくとも表す別個のデータを取得する並行ステップ２７は、画像のシーケンスのビデオ・ベースの解析を実行することによって得られる。図示した実施形態では、一つの動き信号が得られて、その周期的動きに対応する融合された信号の成分を少なくとも抑制する（ステップ２８）ために使われる状況が示されている。

ある代替的な実施形態では、諸画像の別個の諸部分について、別個の動き信号が取得され、それらの部分に対応する諸空間時間的体積について得られた諸信号セクションの成分を抑制するために使われることができる。ある個別的実施形態では、別個の空間時間的体積から抽出された信号セクションの少なくとも周期的成分を抑制するために別個の動き信号を取得するかどうかは、動き信号間の相違の度合いに依存する。

動き抑制をある型の動きに限定するために、解析ステップ２７は一般に、たとえば動きベクトルを計算することによって、動きを定量化する。この場合、動き信号は、動きベクトルの反復性を決定することによって得ることができる。反復周波数に対応する範囲内で、融合された信号の、あるいは抽出された信号セクションの諸成分を抑制するためである。

図６は、上記の諸方法の基礎になる動き抑制の原理を実装する代替的なシステムを示している。本システムのうち図１のシステムに対応する部分は、同様の参照符号を与えられている。図６のデータ処理システム６０１は、生体中の周期的現象を表す信号を捕捉するための少なくとも一つのセンサー６１０へのインターフェース６０９をさらに有することによって異なっている。

ある実施形態では、センサー６１０は、生体の心拍を表す成分を少なくとも含む信号を捕捉するよう構成される。ある変形では、センサー６１０は、生体と接触した際にそのような信号を捕捉するよう適応される。たとえば、センサー６１０は、一つまたは複数の心電図信号を捕捉するために人の身体の近くに保持されるよう構成された電極を有することができる。

データ処理システム６０１は、ビデオ・カメラ６０６から画像のシーケンスを同時に受信する。センサー６１０によってモニタリングされる生体の身体の少なくとも一部の動きの少なくとも周期的成分を表す信号を決定するために、ビデオ・ベースの解析が実行される。これは、必ずしも、センサー６１０があてがわれるのと同じ身体部分である必要はない。たとえば、センサー６１０がルームサイクル〔トレーニング用自転車〕のハンドルに統合されていて、ビデオ・カメラ６０６はルームサイクルに乗った人の上体の画像を取り込むことができる。センサー６１０からの信号は一般に、運動する人の上体の揺れによって影響されるので、センサー６１０からの信号の対応する成分を抑制するために、この揺れる動きを表す信号を決定することは有用である。これは、出力装置６０８を介してユーザーに補正された信号に基づく情報を提供するデータ処理システム１によって実行される。

ある実施形態では、データ処理システムは、ビデオ・ベースの解析の一部として、動き、少なくとも一つの時間変動する動きベクトルを計算する。次いで動きベクトル値の反復性が解析されて、動きの少なくとも周期成分を表す信号が得られる。この信号は、センサー６１０からの信号を補正するために使われる。

ある種の実施形態では、ある型の動きの少なくとも周期成分を表す信号が得られる。たとえば、特定の平面内の動きのみが判別されることができる。代替的にまたは追加的に、ある大きさまたはある周波数範囲内の動きのみが判別されることができる。

センサー６１０からの信号は、信号中の周期的な動きに対応するブロッキング周波数によって補正されることができる。もう一つの実施形態では、時間‐周波数（たとえばウェーブレット）解析が動き信号およびセンサー６１０からの信号の両方に対して実行される。これは、周期成分が比較的小さな一連の時間的窓について別個に解析されるという意味で、非反復性の動きの推定および補償を許容する。そのような時間的窓の大きさは、動きの期待される強度およびビデオ・カメラ６０６と被験体の間の距離を含む種々のパラメータに依存することができる。

次に、本発明のあるさらなる実施形態について述べる。肌色の変動を測定する方法が、顔の画像を撮るカメラに基づき、手動で選択された皮膚領域（たとえば頬）にわたるピクセル平均を計算する。この平均信号の周期的変動を見ることによって、心拍数が実験室場面において抽出できる。

前記カメラによって得られる平均信号は、心拍信号に対応する非常に小さな信号を含み、個人の動きによって容易に汚染されてしまう。ルームランナーで走っている場合のように個人の動きが周期的であり、その周波数が関心のある周波数範囲内、たとえば40BPMから210BPMまでの間にはいる場合、動きの周波数を心拍の周波数から区別することは難しい。一例として、ルームランナーで走っている個人から得られた信号のスペクトログラムが図７にプロットされている。

図７では、横軸は分単位での時間に対応し、縦軸は毎秒拍数（周波数）に対応する。150BPMのまわりに水平方向の暗い線が明瞭に見て取れる。この線が暗いほど、この周波数での信号のエネルギーが大きい。この暗い線が動きに対応する。さらに、150BPM、300BPMおよび若干薄いが450BPM周辺により低い強度で三本の線が見える。これらの線は心拍信号の基本周波数および二つの高調波に対応する。動き信号のエネルギーが心拍信号に比べて大きいので、明らかに、動きが信号を有意に汚染している。よって、心拍数を正しく推定するためには、動きの影響を減らす必要がある。さらに、信号から動きを除去することが可能であれば、ずっとクリーンな心拍信号を得ることができ、それによりピーク・ツー・ピーク検出によって時間領域で心拍数を推定することが可能になる。クリーンな信号の追加的な利点は、HRV（心拍数変動度）が推定できるということである。

時間において関心のある皮膚領域の動きを記述する（平均）動きベクトルが非常にクリーンな信号であり、精密に動きを記述することが見出されている。さらに、動きベクトルはいかなる心拍信号も含まず、よって（たとえば非特許文献３に記載されるような）適応フィルタなどのための、平均肌色信号における動きに誘起された歪みを低減させるための参照信号として使うのに非常に好適である（適応エコー打ち消しの用語では、これは遠方端信号（far-end signals）と呼ばれる）。

各ピクセルの個々の動きベクトルに目を向けると、これらの動きベクトルはみな同じではないことが見て取れる（たとえば、頭の回転がある場合、並進についてはこれらのベクトルはより似通ったものになる）。好ましくは、単にすべてのピクセル（またはあるグループのピクセル）の平均動きベクトルが取られ、参照信号として使われる。オブジェクト、たとえば顔を諸領域に分割し、ローカルに補償するためにそれらの領域の平均動きベクトルを使うことも可能である。

図８は、心拍信号における（一般に、図２または図３を参照して上記した処理によって得られる信号、特にステップ１７またはステップ２６の出力における）動きに誘起された汚染を減らすために、上記のシステムにおいて、好ましくは処理ユニット７または６０７の一部として使うための動き補償ユニット８００のある実施形態のブロック図を描いている。このように、前記動き補償ユニット８００は、たとえば、図２または図３を参照して例解した前記諸方法のステップ１９またはステップ２８を実装するよう適応される。

動き補償ユニット８００は、入力x動き信号８０９および入力y動き信号８１０として与えられる、関心のある皮膚領域の（平均）x動きおよびy動きをフィルタ処理する二つの適応フィルタ８０１、８０２を有する。これらの動き信号８０９、８１０はたとえば、提案される方法の上記の諸実施形態のステップ１８またはステップ２７によって与えられてもよいし、および／または図６に示されるセンサー６１０によって与えられてもよい。（x方向およびy方向の動きについての）二つのフィルタという数および二つの動き信号の使用は必須ではないことを注意しておくべきである。一般に、所望される方向（たとえば動きが最も強いと想定される方向）における一つのフィルタおよび一つの動き信号、あるいは（所望される諸方向における）より多くの動き信号を使うこともできる。

この実施形態では、適応フィルタ８０１、８０２の係数は、規格化ユニット８０３、８０４内の規格化最小平均二乗（NLMS: normalized least-mean-square）アルゴリズムによって更新されるが、他の適応アルゴリズムが使われることもできる。フィルタ処理されたx動き信号およびy動き信号８０５、８０６は、減算ユニット８０７において入力信号８０８から減算される。入力信号８０８は、動きに誘起された信号により汚染された、カメラによって得られた心拍信号に対応する。規格化ユニット８０３、８０４は、出力信号８１１（理想的には、動きに誘起された汚染が除去された心拍信号）のエネルギーが最小化されるよう、係数を適応させる。x動き信号およびy動き信号８０９、８１０は心拍信号を含まないので、この最小に到達するのは、出力信号８１１が動きと相関したいかなる信号も含まないときである。よって、本実施形態によれば、適応フィルタ８０１、８０２は好ましくは単に、動きに起因する歪みが減らされるよう、参照信号の振幅および位相を変える。

図９では、関心対象の皮膚領域の平均動きベクトルが上記適応フィルタによって心拍信号中の動きに誘起された信号を減らすために使われるときの、ルームランナーで走っている個人から得られた結果のスペクトログラムがプロットされている。図７に描かれたプロットに比べ、動きに誘起された汚染が著しく除去されていることが明らかに見て取れる。この場合、両適応フィルタ８０１、８０２について四つの係数のみが使われる。

図７ないし図９を参照して上述した実施形態では、x動きおよびy動きが完全に無相関であり、二つの規格化ユニット８０３、８０４が独立して適応することが想定されていた。実際上、これは必ずしも成り立たず、一般には、二つの入力動き信号８０９、８１０の間の相関を取り入れるために、適応フィルタ８０１、８０２としての複数入力単一出力の適応フィルタが使われる。しかしながら、相関が小さい場合、上記の実施形態は上例が示すように非常によく機能する。たとえば、本願によれば、心拍信号を信頼できる形で抽出することが可能であり、該心拍信号は動きに誘起された歪みの削減によってさらに改善される。
上述の実施形態が本発明を限定するのではなく例解しているのであり、付属の請求項の範囲から外れることなく当業者は多くの代替的な実施形態を設計できるであろうことを注意しておくべきである。請求項において、括弧内に参照符号があったとしても、その請求項を限定するものと解釈してはならない。「有する」「含む」の語は、請求項に挙げられている以外の要素やステップの存在を排除するものではない。要素の単数形の表現はそのような要素の複数の存在を排除するものではない。ある種の施策が互いに異なる従属請求項に記載されているというだけの事実がそれらの施策の組み合わせが有利に使用できないことを示すものではない。

コンピュータ・プログラムは、他のハードウェアと一緒にまたは他のハードウェアの一部として供給される光記憶媒体または半導体媒体のような好適な非一時的媒体上で記憶／頒布されてもよいが、インターネットまたは他の有線もしく無線遠隔通信システムを介するなど他の形で頒布されてもよい。

図２の方法において類似のサブセットを見出すために使われる類似性の指標はテクスチャーの類似性に限定されることができる。ピクセル値の厳密な類似性は要求されない。
いくつかの態様を記載しておく。
〔態様１〕
生体中の少なくとも周期的現象を表す成分を含む信号を処理する方法であって：
生体中の少なくとも周期的現象を表す成分を有する少なくとも一つの第一の信号を取得する段階と；
前記少なくとも一つの第一の信号と同期して、前記生体を示す画像のシーケンスを取得して、前記画像のシーケンスのビデオ・ベースの動き解析を実行することによって、前記生体の動きの少なくとも周期的成分を表す別個のデータを取得する段階と；
前記データを使って、少なくとも、前記生体の動きの前記周期的成分に対応する前記第一の信号の成分を抑制する段階とを含む、
方法。
〔態様２〕
前記少なくとも一つの第一の信号を取得する段階は、少なくとも一つの第一の信号を前記画像のシーケンスから抽出することを含む、態様１記載の方法。
〔態様３〕
前記少なくとも一つの第一の信号を前記画像のシーケンスから抽出することが、それぞれの画像からのピクセル値の組み合わせに対応する変数の値を得ることを含む、態様２記載の方法。
〔態様４〕
前記少なくとも一つの第一の信号を前記シーケンスから抽出することが、前記シーケンスにわたる前記変数の変動のみを決定することを含む、態様３記載の方法。
〔態様５〕
前記画像において表現されている少なくとも一つの被験体について少なくとも一つの動きベクトルを計算し、動きベクトル値の反復性を解析することを含む、態様１記載の方法。
〔態様６〕
前記画像において表現されているある被験体の複数の部分のそれぞれについて動きベクトル値を計算することを含む、態様５記載の方法。
〔態様７〕
動きの少なくとも周期的成分を表す別個のデータを取得する段階が、前記生体の少なくとも一部の空間的な変位を定量化することを含む、態様１記載の方法。
〔態様８〕
前記画像において表現されている少なくとも一つの被験体について少なくとも一つの動きベクトルを計算する段階と、
前記少なくとも一つの動きベクトルをフィルタ処理する段階と、
前記第一の信号から前記フィルタ処理された動きベクトルを減算して、前記第一の信号の、前記生体の動きの前記周期的成分に対応する成分を抑制し、擾乱が軽減された第一の信号を得る段階とを含む、
態様１記載の方法。
〔態様９〕
前記フィルタ処理が、前記フィルタ処理の係数が適応的に更新される適応的フィルタ処理を含む、態様８記載の方法。
〔態様１０〕
前記フィルタ処理の係数が、擾乱が軽減された第一の信号のエネルギーが最小化されるよう更新される、態様９記載の方法。
〔態様１１〕
特に第一の方向および前記第一の方向に垂直な第二の方向の、少なくとも二つの動きベクトルが別個に計算され、別個にフィルタ処理される、態様８記載の方法。
〔態様１２〕
生体中の少なくとも周期的現象を表す成分を含む信号を処理するシステムであって：
生体中の少なくとも周期的現象を表す成分を有する少なくとも一つの第一の信号を取得する少なくとも一つのコンポーネントと；
前記少なくとも一つの第一の信号と同期して、前記生体を示す画像のシーケンスを取得するインターフェースおよび前記画像のシーケンスのビデオ・ベースの動き解析を実行する画像データ処理システムを含む、前記生体の動きの少なくとも周期的成分を表す別個のデータを取得する少なくとも一つのコンポーネントと；
前記データを使って、少なくとも、前記生体の動きの前記周期的成分に対応する前記第一の信号の成分を抑制するよう構成された少なくとも一つのデータ処理コンポーネントとを含むシステム。
〔態様１３〕
態様１ないし１１のうちいずれか一項記載の方法を実行するよう構成された、態様８記載のシステム。
〔態様１４〕
機械可読媒体に組み込まれたときに、情報処理機能をもつシステムに態様１ないし１１のうちいずれか一項記載の方法を実行させることのできる一組の命令を含むコンピュータ・プログラム。

Claims (9)

1. 生体中の少なくとも周期的現象を表す成分を含む信号を処理する方法であって：
   前記生体を示す画像のシーケンスを取得する段階と；
   複数の関心領域のそれぞれについてリモート光体積変動記録法を使って前記画像のシーケンスから生体中の少なくとも周期的現象を表す成分を有する少なくとも一つの第一の信号を抽出する段階と；
   前記画像のシーケンスのビデオ・ベースの動き解析を実行することによって、各関心領域について前記生体の動きの少なくとも周期的成分を表す別個のデータを取得する段階であって、前記別個のデータは前記画像の間の少なくとも一つの動きベクトルを含む、段階と；
   前記少なくとも一つの動きベクトルに対して、適応的に更新される係数を用いた適応的フィルタ処理を行なう段階と；
   前記第一の信号から前記フィルタ処理された動きベクトルを減算して、前記第一の信号の、前記生体の動きの前記周期的成分に対応する成分を抑制し、擾乱が軽減された第一の信号を得る段階であって、前記抑制を、さまざまな前記関心領域の動きの相違が所定の量より大きい場合には、各関心領域について、関連する関心領域の関連する別個のデータを使って個々に行い、さまざまな前記関心領域の動きの相違が前記所定の量より小さい場合には、それらの関心領域に関連する前記別個のデータに基づいて形成される単一の別個のデータに基づいて行う段階とを含み、
   前記適応的フィルタ処理の係数が、擾乱が軽減された第一の信号のエネルギーが最小化されるよう更新される、
   方法。
2. 前記少なくとも一つの第一の信号を前記画像のシーケンスから抽出することが、それぞれの画像からのピクセル値の組み合わせに対応する変数の値を得ることを含む、請求項１記載の方法。
3. 前記少なくとも一つの第一の信号を前記シーケンスから抽出することが、前記シーケンスにわたる前記変数の変動のみを決定することを含む、請求項２記載の方法。
4. 前記画像において表現されている少なくとも一つの被験体について少なくとも一つの動きベクトルを計算し、動きベクトル値の反復性を解析することを含む、請求項１記載の方法。
5. 前記画像において表現されているある被験体の複数の部分のそれぞれについて動きベクトル値を計算することを含む、請求項４記載の方法。
6. 第一の方向および前記第一の方向に垂直な第二の方向の、少なくとも二つの動きベクトルが別個に計算され、別個にフィルタ処理される、請求項１記載の方法。
7. 生体中の少なくとも周期的現象を表す成分を含む信号を処理するシステムであって：
   前記生体を示す画像のシーケンスを取得するインターフェースと；
   複数の関心領域のそれぞれについてリモート光体積変動記録法を使って前記画像のシーケンスから生体中の少なくとも周期的現象を表す成分を有する少なくとも一つの第一の信号を抽出する少なくとも一つのコンポーネントと；
   前記画像のシーケンスのビデオ・ベースの動き解析を実行する画像データ処理システムを含む、各関心領域について前記生体の動きの少なくとも周期的成分を表す別個のデータを取得する少なくとも一つのコンポーネントであって、前記別個のデータは前記画像の間の少なくとも一つの動きベクトルを含む、コンポーネントと；
   前記少なくとも一つの動きベクトルに対して、適応的に更新される係数を用いた適応的フィルタ処理を行なうコンポーネントと；
   前記第一の信号から前記フィルタ処理された動きベクトルを減算して、前記第一の信号の、前記生体の動きの前記周期的成分に対応する成分を抑制し、擾乱が軽減された第一の信号を得るコンポーネントであって、前記抑制を、さまざまな前記関心領域の動きの相違が所定の量より大きい場合には、各関心領域について、関連する関心領域の関連する別個のデータを使って個々に行い、さまざまな前記関心領域の動きの相違が前記所定の量より小さい場合には、それらの関心領域に関連する前記別個のデータに基づいて形成される単一の別個のデータに基づいて行うコンポーネントとを含み、
   前記適応的フィルタ処理の係数が、擾乱が軽減された第一の信号のエネルギーが最小化されるよう更新される、
   システム。
8. 請求項１ないし６のうちいずれか一項記載の方法を実行するよう構成されたシステム。
9. 機械可読媒体に組み込まれたときに、情報処理機能をもつシステムに請求項１ないし６のうちいずれか一項記載の方法を実行させることのできる一組の命令を含むコンピュータ・プログラム。

Images