(57)【特許請求の範囲】
【請求項１】
生体の一部を照明する光源部と、照明された生体の一部を撮像する撮像部と、撮像画像中に設定された解析領域を解析し生体情報を得るデータ処理部と、得られた生体情報を出力する出力部と、各種操作等を行うための操作部と備え、
データ処理部は、撮像画像の特徴パターンを抽出する機能を備え、第１解析領域が設定済の第１撮像画像の特徴パターンを解析領域未設定の第２撮像画像の特徴パターンに対して相対移動させながら特徴パターンの一致の程度を表す評価値を算出し、その評価値を所定の条件に合致させるために必要な特徴パターンの相対移動量に基いて、第１解析領域に含まれる生体部位と同一の生体部位を含む第２解析領域を第２撮像画像中に設定する機能、および、第１解析領域および第２解析領域を解析することにより得られた生体情報の時系列変化を表示する表示画面を出力部に表示させる機能を備えることを特徴とする無侵襲生体計測装置。
【請求項２】
上記データ処理部の機能として、第１撮像画像の特徴パターンを第２撮像画像の特徴パターンに対して所定距離ずつ相対移動させ、その都度特徴パターンの一致の程度を表す評価値を算出する機能を有することを特徴とする請求項１記載の無侵襲生体計測装置。
【請求項３】
撮像画像の特徴パターンを抽出する機能が、撮像画像を構成する着目画素ごとにその着目画素を中心とする局所領域の平均的濃度を算出しこれら各平均的濃度を用いて各着目画素を２値化していくことにより撮像画像の特徴パターンを抽出するものであることを特徴とする請求項１記載の無侵襲生体計測装置。
【請求項４】
撮像画像の特徴パターンが、血管部分を抽出して得られた血管パターンであることを特徴とする請求項１記載の無侵襲生体計測装置。
【請求項５】
各解析領域は同一の生体に対して同一の血管部位に設定され、生体情報として血管寸法又は血液成分濃度に関する情報を計測することを特徴とする請求項１記載の無侵襲生体計測装置。
【請求項６】
生体の一部を照明する照明ステップと、
照明された生体の一部を撮像し、撮像画像を取得する撮像ステップと、
撮像ステップにおいて取得した第２撮像画像の特徴パターンを抽出する抽出ステップと、
第１解析領域が設定済である第１撮像画像の特徴パターンをメモリーから読み出す読出ステップと、
抽出ステップにおいて抽出した特徴パターンを読出ステップにおいてメモリーから読み出した特徴パターンに対して相対移動させながら特徴パターンの一致の程度を表す評価値を算出し、その評価値を所定の条件に合致させるために必要な特徴パターンの相対移動量に基いて、撮像ステップにおいて取得した第２撮像画像中に、第１解析領域に含まれる生体部位と同一の生体部位を含む第２解析領域を設定する解析領域設定ステップと、
第１解析領域および第２解析領域を解析することにより得られた生体情報の時系列変化を表示する表示画面を表示する表示ステップと、を有することを特徴とする無侵襲生体計測方法。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、生体に光を照射し生体から得られた画像情報を解析し生体の形態や生体成分濃度などの生体情報を得る無侵襲生体計測装置および無侵襲生体計測方法に関するものである。
【０００２】
【従来の技術】
採血することなく（無侵襲的に）ヘモグロビンなどの生体成分を計測しようとする装置として国際公開第ＷＯ９７／２４０６６号公報に開示されたものがある。その装置の場合、検出部に指を配置し指の透過画像を得、その画像中に解析領域（コントラストのよい血管部位）を設定し、その解析領域を解析しヘモグロビン濃度を算出している。
【０００３】
【発明が解決しようとする課題】
検出部から指を抜き出して、再度指を配置すると指の位置が先程とは異なってしまう。すると得られる撮像画像が左右あるいは上下にずれた画像となるので、同一の生体部位を計測することが難しくなる。このことは、ある特定の血管部位について血管幅や生体成分濃度の経時変化を計測しようとする場合に不都合を生じさせる。
【０００４】
例えば、特開平７−２１３７３号公報には撮像された血管像から得られた血管パターンから特徴点（端点と分岐点）をピックアップし、その特徴点の位置的な配列、血管の方向ベクトル、特徴点間の連結状態、あるいはこれらの組合せを特徴量とし、個人識別を行う装置が開示されているが、無侵襲生体計測装置における生体の位置ずれという問題に対する解決策については何ら開示されていない。
ただ、「登録時と照合時とで指の方向がずれる場合も予想されるため、比較の際には所定の許容範囲を設けてもよい。」との記載があるのみである。
【０００５】
一方、画像処理技術として、濃度むらが生じている画像からその特徴を抽出しようとして画像の２値化を行うことも知られているが、上記無侵襲生体計測装置における生体の位置ずれの問題と関連づけて説明されたものは知られていない。
【０００６】
本発明は、装置に対する生体の配置位置が異なっても同一生体に対して同一部位を計測することができる無侵襲生体計測装置を提供することを課題とする。
【０００７】
【課題を解決するための手段】
本発明の無侵襲生体計測装置は、生体の一部を照明する光源部と、照明された生体の一部を撮像する撮像部と、撮像画像中に設定された解析領域を解析し生体情報を得るデータ処理部と、得られた生体情報を出力する出力部と、各種操作等を行うための操作部と備え、データ処理部は、撮像画像の特徴パターンを抽出する機能を備え、第１解析領域が設定済の第１撮像画像の特徴パターンを解析領域未設定の第２撮像画像の特徴パターンに対して相対移動させながら特徴パターンの一致の程度を表す評価値を算出し、その評価値を所定の条件に合致させるために必要な特徴パターンの相対移動量に基いて、第１解析領域に含まれる生体部位と同一の生体部位を含む第２解析領域を第２撮像画像中に設定する機能、および、第１解析領域および第２解析領域を解析することにより得られた生体情報の時系列変化を表示する表示画面を出力部に表示させる機能を備えることを特徴とする。
また、本発明の無侵襲生体計測方法は、生体の一部を照明する照明ステップと、照明された生体の一部を撮像し、撮像画像を取得する撮像ステップと、撮像ステップにおいて取得した第２撮像画像の特徴パターンを抽出する抽出ステップと、第１解析領域が設定済である第１撮像画像の特徴パターンをメモリーから読み出す読出ステップと、抽出ステップにおいて抽出した特徴パターンを読出ステップにおいてメモリーから読み出した特徴パターンに対して相対移動させながら特徴パターンの一致の程度を表す評価値を算出し、その評価値を所定の条件に合致させるために必要な特徴パターンの相対移動量に基いて、撮像ステップにおいて取得した第２撮像画像中に、第１解析領域に含まれる生体部位と同一の生体部位を含む第２解析領域を設定する解析領域設定ステップと、第１解析領域および第２解析領域を解析することにより得られた生体情報の時系列変化を表示する表示画面を表示する表示ステップと、を有することを特徴とする。
【０００８】
わかりやすく説明する。
第１の撮像画像Ｇ１と第２の撮像画像Ｇ２があり、第１の撮像画像Ｇ１については解析領域Ａ１が設定されていて、第２の撮像画像については解析領域が設定されていないとする。各撮像画像Ｇ１、Ｇ２について特徴パターンＨ１、Ｈ２を抽出し、両特徴パターンＨ１、Ｈ２のパターンをよく一致させるために必要な相対移動量Δを求める。第１撮像画像Ｇ１における解析領域Ａ１の位置と上記相対移動量Δとから第２の撮像画像Ｇ２における解析領域Ａ２の位置を決定することができる。
【０００９】
このようにして同一生体に対して得られた異なる撮像画像に対して同一の生体部位を解析領域として設定し計測することができ、このことにより、生体情報の再現性の向上や経時変化の測定が可能となる。
【００１０】
【発明の実施の形態】
この発明の無侵襲生体計測装置において、生体の一部とは生体（好適例としては哺乳類）から分離した組織ではなく生体のありのままの組織の一部であり、計測に適する部分であればよい。ヒトの場合には指や耳朶などがその一例である。
ラットやマウスの場合には尻尾であってもよい。生体情報とは生体組織の形態（形状や大きさや数など）や生体成分の濃度等に関する情報である。具体的には血管の寸法や血液成分の濃度（例えば、ヘモグロビン、ヘマトクリット等）や血液成分の濃度比（例えば、血液の酸素化率等）などである。
【００１１】
操作や取り扱いの利便性の観点から、本装置は生体から画像情報を得る検出部（検出ユニット、検出装置）と得られた画像情報を解析する解析部（解析ユニット、解析装置）とに分離的に構成されていることが好ましい。両者はコネクタあるいはケーブル等の情報伝送手段により接続される。検出部は光源部及び撮像部を備え、解析部には生体情報を解析するデータ処理部、出力部及び操作部を備えることができる。
【００１２】
光源部には、半導体レーザ（ＬＤ）やＬＥＤやハロゲンランプなどの光源が使用でき、光源からの光を生体に直接、あるいは光ファイバなどを介して照射する。照射する光の波長としては生体組織を透過し、水の吸収が大きくない６００〜９５０ｎｍの範囲にあることが好ましい。
【００１３】
撮像部は、レンズなどの光学素子とＣＣＤなどの撮像素子により構成することができる。生体からの画像情報は撮像素子により生成される。
【００１４】
より好適な画像情報を得るために、検出部に生体の一部を光学系に対して保持する保持部材を備えることが好ましい。生体の一部が例えばヒトの手の指である場合には、保持部材として光源部と撮像部との間において指を離脱可能に圧迫することなく保持するものを使用することができる。例えば、可動片で指を両側から挟持する方式のものを用いることができる。
【００１５】
解析部としては、ＣＰＵ、ＲＯＭ、ＲＡＭおよびＩ／Ｏポートからなるマイクロコンピュータあるいは市販のパーソナルコンピュータを利用することができ、撮像画像から生体情報を解析するデータ処理部、各種データを入力したり各種操作を行ったりする操作部及び生体の画像情報の解析結果を出力する出力部を備えてもよい。出力部としてはＣＲＴや液晶ディスプレイなどの表示装置やプリンタ等の印字装置を使用でき、操作部としてはキーボードやテンキーやタッチキーなど使用することができる。
【００１６】
データ処理部は、撮像画像の特徴パターンを抽出する機能と、特徴パターン同士をよく一致させるために必要な特徴パターンの相対移動量に基いて解析領域未設定の撮像画像の解析領域を設定する機能を備えている。このことにより、同一の生体について所定時間置いて（所定日数置く場合も含む）再配置され得られた複数の撮像画像に対して同一部位に解析領域を設定し、その解析領域における生体情報を計測することができる。すなわち、同一生体に対して得られた異なる撮像画像に対して同一の部位を計測することができるので生体情報の精度や再現性を向上させることができる。また、生体情報の経時変化（経日変化を含む）を計測することも可能となる。各解析領域が同一血管部位に設定される場合には、生体情報としてその血管部位の血管寸法やその血管部位を流れる血液成分の濃度を得ることができる。
【００１７】
上記特徴パターン同士をよく一致させるに必要な相対移動量を決定するについて、一方の特徴パターンを他方の特徴パターンに対して相対移動させながらその都度特徴パターンの一致の程度を表わす評価値を算出しその評価値に基いて決定するようにすることが実用的である。
【００１８】
なお、特徴パターンの一致度が所定以上とならない場合には警告を発することが好ましい。
【００１９】
濃度むらのある撮像画像からより好ましく特徴パターンを抽出するために、撮像画像を構成する着目画素ごとにその着目画素を中心とする局所領域の平均的濃度を算出しこれら各平均的濃度を用いて各着目画素を２値化していくことにより撮像画像の特徴パターンを抽出するものであることが好ましい。
【００２０】
さらに、処理の簡略化、高速化の観点から、各局所領域の平均的濃度を算出する際、その局所領域から選出された一部の画素についてそれらの平均濃度を算出しその平均濃度をその局所領域の平均的濃度とすることが好ましい。
【００２１】
また、撮像画像の特徴パターンとして、血管部分を抽出して得られた血管パターンであっても良い。血管パターンは、撮像画像を構成する輝度に着目し、その輝度が周囲に比べて低くなっている点を結ぶことにより得ることができる。
【００２２】
【実施例】
図１は本発明の無侵襲生体計測装置の一実施例の外観図であり、マウス様の検出部１をノート型パソコン様の解析部３の載置部４に載置した状態を示している。図２は図１の無侵襲生体計測装置を機能的に表現した全体ブロック図である。
【００２３】
載置部３に載置された検出部１は、検出部１側のコネクタ７と解析部３側のコネクタ８とで電気的にまた機構的にも連結されている。計測や取り扱いの自由度を得るために検出部１を解析部３の載置部４から取り外し、コネクタ付き接続コードで検出部１と解析部３とを接続して使用することも可能である。
【００２４】
検出部１はアーム５とアーム５に対して回転移動できるハウジング６とからなり、アーム５の内側に光源部１１が内蔵され、アーム５に対向するハウジング６の内側に撮像部１２を内蔵している。アーム５とハウジング６の間にヒトの指１４を保持し、指１４の第２関節部分において指の背側に光を照射し指の腹側からその透過光像を撮像する。指は指の両側から弾性的に挟持されている。
【００２５】
光源部１１は、波長の異なる複数のＬＥＤを備えた発光素子からなる。第１のＬＥＤとして、中心波長８３０ｎｍ、半値幅４０ｎｍのＬ３９８９（浜松ホトニクス（株）製）を使用し、第２のＬＥＤとして中心波長８９０ｎｍ、半値幅５０ｎｍのＬ２５６５６（同上製）を使用している。血管幅を計測する際は第１のＬＥＤのみを点灯させ、血液成分濃度を計測する際は第１及び第２のＬＥＤを点灯させている。
【００２６】
解折部３はデータ処理部２、出力部（液晶モニター）１０及び操作部（複数のキー）９からなる。解析部３には外部記憶媒体であるフロッピーディスクの挿入口１３が設けられ測定情報等を外部記憶することができる。
【００２７】
図２は本実施例の無侵襲生体計測装置の全体ブロック図である。
解析部３のデータ処理部２の実体は、ＣＰＵ、メモリを備えるコンピュータであり、検出部１の撮像部１２で得られた指１４の撮像画像の特徴パターン生成する特徴抽出機能、特徴パターンを記憶する機能、異なる撮像画像の特徴パターン同士をよく一致させるに必要な特徴パターンの相対移動量を決定する機能、その相対移動量に基づき撮像画像に対して同一の血管部位を解析領域として設定する解析領域設定機能、解析領域内の血管部分の濃度分布を濃度プロファィルとして抽出する濃度プロファイル抽出機能、抽出された濃度プロファイルの形態的特徴を定量化する定量化機能、定量化された特徴に基づいて血管寸法や血液成分濃度などの生体情報を演算する演算機能などを備える。
【００２８】
図３、４、５のフローチャートを参照しつつ、具体的に血管幅の計測手順を説明する。図４、図５はそれぞれ解析領域設定、相対移動量算出のフローチャートである。ここでは一つの波長を使って同一人の指を経時的に複数回撮像し、同一血管部位における血管幅の経時変化を計測している。
【００２９】
まず、第１のＬＥＤ（第１波長）によって指１４を照明して撮像し、図６に示すように、指１４の皮膚表面に存在する血管（静脈）像を合む画像（ｘ軸方向６４０画素×ｙ軸方向４８０画素）を得る（ステップＳ１）。
【００３０】
本装置は被検者ごとに記憶ファイルを有している。この被検者について過去に計測が行われていなければ（つまり記憶ファイルがない、あるいは記憶ファイルがあっても該当データがない）、第１回目の測定であるとし（ステップＳ３１）、その撮像画像Ｇ１（ｘ，ｙ）中から血管が最もコントラストよく写っている領域を探索しその血管を含む領域を解析領域Ａ１（ｘ，ｙ）として設定する（ステップＳ３２）。解析領域Ａ１（ｘ，ｙ）は通常、自動設定されるが、使用者が出力部１０に表示されるモニタ像を見ながら操作部９を操作し手動的に設定してもよい。
【００３１】
次に、撮像画像Ｇ１（ｘ，ｙ）の特徴パターンを抽出し（ステップＳ３３）、特徴パターンをメモリに記憶する。解折領域Ａ１（ｘ，ｙ）の位置情報も解析領域Ａ１（ｘ，ｙ）の領域境界である四角形の各頂点Ｐ、Ｑ、Ｒ、Ｓの座標（ｐ１，ｐ２）、（ｑ１，ｑ２）、（ｒ１，ｒ２）、（ｓ１，ｓ２）としてメモリに記憶する（ステップＳ３４）。
この被検者について第１回目の計測であれば処理はこれで終了する。
【００３２】
ところで、画像の特徴パターンを抽出するため画像を２値化するのであるが、撮像画像Ｇ１（ｘ，ｙ）は光の照射むらや生体の厚みや特性の違い等により濃度むらの大きい画像となっているので、一定の閾値で２値化すると画像の特徴をよく捉えることができない。
【００３３】
そこで画像の各画素ごとに異なる閾値を設定することによりその画像が特徴的に有している二次元パターンを抽出する。すなわち、図１０に示すように着目画素（ｘ，ｙ）（黒色で示す）に対してその着目画素（ｘ，ｙ）を中心とする局所領域Ｄ（ｘ，ｙ）を想定し、その局所領域Ｄ（ｘ，ｙ）の平均的濃度Ｄｘｙを着目画素（ｘ，ｙ）の閾値とする。本実施例では、局所領域Ｄ（ｘ，ｙ）を１７画素×１７画素の領域とし、さらに、高速処理を行うため１画素置きに選んだ、着目画素を含む８１個の画素（灰色で示す）の濃度の平均値を局所領域Ｄ（ｘ，ｙ）の平均的濃度としている（選択画素は｛（ｘ−ｉ，ｙ−ｊ）｜ｉ，ｊ＝０，±２，±４，±６，±８｝）。ところで、本実施例では局所領域としては９画素×９画素から３３画素×３３画素が適していた。
【００３４】
すべての画素について各着目画素（ｘ，ｙ）の濃度と各閾値Ｄｘｙの大小比較を行うことにより、図８に示すような、撮像画像Ｇ１（ｘ，ｙ）の特徴がよく現れた特徴パターンを得ている。この特徴パターンには血管や指のしわなどの形態が現れている。この特徴パターンは６４０画素×４８０画素の全領域ではなく明るさ中心付近のｘ軸方向２４０画素×ｙ軸方向２６０画素の領域を特徴パターンＨ１（ｘ，ｙ）としてを記憶するようにしている。これは処理の高速化に寄与している。図８中にこの特徴パターンＨ１（ｘ，ｙ）の領域を大きい方の長方形の線で示し、特徴パターンＨ１（ｘ，ｙ）の領域内に設定される解析領域Ａ１（ｘ，ｙ）の領域を小さい方の長方形の線で示した。
【００３５】
ここで、明るさ中心付近は、輝度重心を中心にした指定領域とした。輝度重心を求める式は、数１、数２で示され、輝度重心は（ＸＧ，ＹＧ）と表示される。
【数１】
0004554739.tif 000002
【数２】
0004554739.tif 000003
【００３６】
さて、次に、計測がその人にとって２回目であるときには（ステップＳ３１）、図７に示すような撮像画像Ｇ２（ｘ，ｙ）が得られると（ステップＳ１）、前記と同様にして図９のごとくｘ軸方向６４０画素×ｙ軸方向４８０画素の特徴パターンＨ２（ｘ，ｙ）を得る（ステップＳ３６）。なお、図７は図６より指が少し左側にシフトしているときに得られる画像である。
【００３７】
次に、上記第１の特徴パターンＨ１（ｘ，ｙ）をメモリから読み出し（ステップＳ３７）、図１１に示すように、第２の特徴パターンＨ２（ｘ，ｙ）との比較を行うことにより両特徴パターンがよく一致する相対移動量を求める（ステップＳ３８）。具体的には第２の特徴パターンＨ２（ｘ，ｙ）に対して第１の特徴パターンＨ１（ｘ，ｙ）を移動させ特徴パターンがよく一致するときの相対移動量Δｘ、Δｙを求める。そして、第１の撮像画像Ｇ１（ｘ，ｙ）における解析領域Ａ１（ｘ，ｙ）の位置と上記相対移動量Δｘ、Δｙとから第２の撮像画像Ｇ２（ｘ、ｙ）における解析領域Ａ２（ｘ，ｙ）の位置を決定する。
【００３８】
図５を参照し、上記特徴パターンＨ１（ｘ，ｙ）、Ｈ２（ｘ，ｙ）の比較について具体的に説明する。
第１の特徴パターンＨ１（ｘ，ｙ）を第２の特徴パターンＨ２（ｘ，ｙ）に対してｘ軸方向にΔｘ、ｙ軸方向にΔｙだけ相対移動させ（ステップＳ４２）、そのときの両パターンの一致度合いを反映する評価値ＥΔｘΔｙを算出する（ステップＳ４３）。評価値ＥΔｘΔｙはそれぞれ対応する２つの画素について値（１か０）が一致しているか否かを調べ、一致している画素の総数としている。この処理を相対移動量を予め設定された範囲内で少しずつ変えながら評価値ＥΔｘΔｙの算出を行い（ステップＳ４１、Ｓ４６）、評価値ＥΔｘΔｙが最大となったｘ軸方向の移動量Δｘ及びｙ軸方向の移動量Δｙを記憶する（ステップＳ４４、４５）。なお、指は指の軸方向（ｘ軸方向）にずれて配置されやすいので、評価値ＥΔｘΔｙの算出プロセスにおいてｘ軸方向には数画素置きに相対移動させながら計算を繰り返すことが高速化の点で好ましい。
【００３９】
ここで、図２１、図２２を用いて、第１の特徴パターンＨ１（ｘ，ｙ）を第２の特徴パターンＨ２（ｘ，ｙ）に対してＸ軸方向△ｘ、Ｙ軸方向に△ｙだけ相対移動させる方法について、詳しく説明する。
【００４０】
第１の特徴パターンＨ１（ｘ，ｙ）は第２の特徴パターンＨ２（ｘ，ｙ）に対して図２１に示すように移動される。図２２のａの状態、すなわち、第１の特徴パターンＨ１（ｘ，ｙ）をＸ軸方向に△ｘ＝−１４０、Ｙ軸方向に△ｙ＝−３０移動させたときの位置をスタートとする。そして、Ｘ軸方向に数画素ずつ右側にずらして行き、Ｘ軸方向に△ｘ＝１４０、Ｙ軸方向に△ｙ＝−３０まで移動させる（ｂ）。次に、Ｘ軸方向に△ｘ＝−１４０まで戻し、Ｙ軸方向に数画素だけ下側にずらした後、Ｘ軸方向に数画素ずつ右側にずらしていく。これを繰り返し行い、Ｘ軸方向に△ｘ＝−１４０、Ｙ軸方向に△ｙ＝３０移動させ（ｃ）、Ｘ軸方向△ｘ＝１４０、Ｙ軸方向に△ｙ＝３０移動させる（ｄ）。移動はこれで完了する。このように、第１の特徴パターンＨ１（ｘ，ｙ）を第２の特徴パターンＨ２（ｘ，ｙ）に対して数画素ずつ移動させ、その都度、評価値を求め、その評価値が最大になるときの相対移動量△ｘ、△ｙを特徴パターンがよく一致するときの相対移動量Δｘ、Δｙとして求めることができる。
【００４１】
本実施例では、相対移動量Δｘ＝６３、Δｙ＝−７であり、そのときの特徴パターンの一致度は６２．１％であった。一致度とは一方の特徴パターンを他方の特徴パターンに対して所定量相対移動させたときの評価値を、両パターンが完全に一致する場合の評価値で割った割合である。
【００４２】
以上のような処理の結果、まだ解析領域が設定されていない第２の撮像画像Ｇ２（ｘ、ｙ）における解析領域Ａ２（ｘ，ｙ）の頂点Ｐ’、Ｑ’、Ｒ’、Ｓ’の座標はそれぞれ（ｐ１＋Δｘ，ｑ１＋Δｙ）、（ｐ２＋Δｘ，ｑ２＋Δｙ）、（ｐ３＋Δｘ，ｑ３＋Δｙ）（ｐ４＋Δｘ，ｑ４＋Δｙ）とすることができる（ステップＳ３９）。このようにして、第２の撮像画像Ｇ２（ｘ，ｙ）に対して第１の撮像画像Ｇ１（ｘ，ｙ）で設定したのと同一の血管部位を解析領域として設定することができる。一人の被検者の指について時系列的（例えば、２時間おき）にｎ回撮像し毎回その撮像画像の相対位置が異なってても、常に同一の血管部位を解析領域Ａ１，Ａ２，・・・，Ａｎとして設定することができる。
【００４３】
次に、図３のステップＳ３おいて、設定された解析領域において血管に垂直な方向の濃度プロファイル（図１２）を作成し、濃度プロファイルをベースラインで規格化することにより、図１３に示すように入射光量に依存しない規格化濃度プロファイルを得る（ステップＳ４）。
【００４４】
この規格化した濃度プロファイル（図１３）の形態的特徴であるピーク高さｈ１、（１／２）ｈ１における分布輻（半値輻）ｗ１を算出し（ステップＳ５）、生体情報を求め記憶する。ここでは半値幅ｗ１を血管幅として算出した（ステップＳ６）。
【００４５】
計測が完了すると、算出した血管幅からその時系列的変化を表わすグラフや表を作成して表示する（ステップＳ７、Ｓ８）。
【００４６】
以上、１波長を使って血管幅を計測する例を紹介したが、複数波長を使えばヘモグロビン濃度を計測することができる。すなわち、各波長に対して得られた規格化された濃度プロファイルの形態からヘモグロビン濃度を算出することができる。詳細は国際公開第ＷＯ９７／２４０６６号公報等を参考にすることができるので、ここでは説明は省略する。
【００４７】
図１４は出力部１０（液晶ディスプレイ）の表示例である。表示画面はグラフ表示領域５１、測定データ表示領域５２、測定情報領域５３、メッセージ領域５４、データ表示形式切替領域５５、コマンド領域５６を備える。
【００４８】
グラフ表示領域５１には、計測された生体情報であるヘモグロビン濃度ＨＧＢや血管幅が時系列のグラフとして表示されている。
【００４９】
測定データ表示領域５２には、上記グラフ上の測定時点の測定データが表示される。測定情報領域５３には、測定番号、測定者名、生年月日、性別等の測定者情報が表示される。メッセージ領域５４はシステムの状態を示すメッセージや、測定者に動作を促すメッセージが表示される。
【００５０】
データ表示形式切替領域５５にはデータ表示形式（グラフ、データ一覧表、撮像画像）を切り替えるためのアイコンが表示される。撮像画像を表示するときは解析領域も示すことが好ましい。
【００５１】
コマンド領域５６は、各種コマンドを実行するアイコンが表示され、アイコン選択によりＰＣカードへの測定データの記録、ファイルの削除設定、データの転送、印刷、メンテナンス等が実行できる。
【００５２】
図１５のグラフ領域５１に示されたグラフは、同一の被験者について、毎日の運動前に計測して得られた血管幅、ヘモグロビンＨＧＢ、酸素化率ＶＯＩと、別途入力されたタイムトライアルの成績（計測毎にデータを入力）の時系列グラフである。また、「カーソル」キーで時系列データのグラフ上の計測点を示すバー（破線で示す）を動かすことができ、該当する測定日、測定時間、測定結果、コメントを表示領域５２に表示することができる。
【００５３】
次に、前記に説明した特徴パターンとは異なる特徴パターンを使用する例について説明する。
【００５４】
まず、画像の特徴パターンとして、輝度が周囲に比べて低く谷状になっている点を結んだ折れ線（血管パターン）を抽出し、記録する。
すなわち、図１５（Ａ）に示すように着目画素（ｘ，ｙ）に対してｙ軸方向に±１０画素、つまり（ｘ，ｙ−１０）から（ｘ，ｙ＋１０）の間の範囲における輝度の平均値の９５％を下回る点を抽出して、血管候補点とする。つまり、ｙ軸方向の連続する２１個の画素集団において、その画素集団の輝度が平均輝度の９５％以下となる画素を血管候補点として選択する。そして、その画素集団をｙ軸方向に１画素ずつ、ずらしていき、ｙ軸方向の端から端まで血管候補点の抽出を繰り返す。次ぎに、ｘ軸方向に１画素分ずらし、前記操作を繰り返し、撮像画像Ｇ１（ｘ，ｙ）の全領域において血管候補点の抽出を行う。
【００５５】
次ぎに、図１５（Ｂ）に示すように±１０画素の範囲内で輝度が最も低い血管候補点、輝度最下点（ａ１、ｂ１）（黒色で示す）を抽出する。ここで、斜線部は血管候補点を示している。隣列において、輝度最下点（ａ１、ｂ１）に対応する点（ａ２，ｂ１）対して±１０画素の範囲内で輝度が最も低い輝度最下点（ａ２、ｂ２）（黒色で示す）を抽出し、輝度最下点（ａ１、ｂ１）と（ａ２、ｂ２）を連結する。撮像画像Ｇ１（ｘ，ｙ）の全領域において、前記の操作を繰り返し、血管パターンを形成する。
【００５６】
このようにして、撮像画像Ｇ１（ｘ，ｙ）は図１８に太い破線示した血管パターンとして表される。この特徴パターンは、第１の実施例と同様に、明るさ中心付近の領域を特徴パターンＨ１（ｘ，ｙ）として記憶している。図１８中にこの特徴パターンＨ１（ｘ，ｙ）の領域を大きい方の長方形の破線で示し、特徴パターンＨ１（ｘ，ｙ）の領域内に設定される解析領域Ａ１（ｘ，ｙ）の領域を小さい方の長方形の破線で示した。
２回目の計測のときに得られる撮像画像Ｇ２（ｘ，ｙ）は、前記と同様にして図１９のごとくｘ軸方向６４０画素×ｙ軸方向４８０画素の特徴パターンＨ２（ｘ，ｙ）を得る。
【００５７】
次ぎに、図１６、１７を参照し、特徴パターンＨ１（ｘ，ｙ）、Ｈ２（ｘ，ｙ）の比較について具体的に説明する。
第１の特徴パターンＨ１（ｘ，ｙ）を第２の特徴パターンＨ２（ｘ，ｙ）に対してｘ軸方向にΔｘ、ｙ軸方向にΔｙだけ相対移動させ、両パターンの一致度合いを反映する評価値ＦΔｘΔｙを算出する（ステップＳ５２、Ｓ５３）。図１６に示したように、ｙ軸方向においてＨ１（ｘ，ｙ）の血管パターン上の点と、Ｈ２（ｘ，ｙ）の血管パターン上の点との距離を調べ、距離が最短となる点間の最短距離の２乗値を求める。そして、これを全てのパターン上の点における最短距離（図１６中の矢印で示した部分）の２乗値を求め、合計することにより、評価値ＦΔｘΔｙを得る。この処理を相対移動量を予め設定された範囲内で少しずつ変えながら評価値ＦΔｘΔｙの算出を行い（ステップＳ５１、Ｓ５６）、評価値ＦΔｘΔｙが最小となったｘ軸方向の移動量Δｘ及びｙ軸方向の移動量Δｙを記憶する（ステップＳ５４、５５）。
【００５８】
このようにして、第１の特徴パターンＨ１（ｘ，ｙ）と第２の特徴パターンＨ２（ｘ，ｙ）が最もよく一致する位置が得られる（図２０）。
【００５９】
【発明の効果】
各撮像画像について特徴パターンを抽出し、両特徴パターンをよく一致させるに必要な相対移動量を用いて解析領域未設定の撮像画像に対して同一生体部位に解析領域を設定することができる。異なる撮像画像に対しても同一生体部位を解析領域として設定し計測することができるので、生体情報計測の再現性の向上や経時変化の測定が可能となる。
【００６０】
上記特徴パターン同士をよく一致させるに必要な相対移動量を決定するについて、一方の特徴パターンを他方の特徴パターンに対して相対移動させながらその都度特徴パターンの一致の程度を表わす評価値を算出しその評価値に基いて決定するようにすれば実用的である。
【００６１】
特徴パターンを得るについて、着目画素ごとにその着目画素を中心とする局所領域の平均的濃度を算出しこれら各平均的濃度を用いて各着目画素を２値化していく場合には、撮像画像Ｇ１（ｘ，ｙ）が濃度むらのある画像であっても画像の特徴をよく捉えることができる。
【００６２】
また、撮像画像の特徴パターンが、撮像画像を構成する輝度に着目し、その輝度が周囲に比べて低くなっている点を結んだ血管パターンであることから、撮像画像Ｇ１（ｘ，ｙ）の画素数が少ないので、特徴パターンの比較において高速化できる。
【図面の簡単な説明】
【図１】この発明の実施例の外観を示す斜視図である。
【図２】この発明の実施例のブロック図である。
【図３】この発明の実施例の動作を示すフローチャートである。
【図４】この発明の実施例の動作を示すフローチャートである。
【図５】この発明の実施例の動作を示すフローチャートである。
【図６】この発明の実施例により得られる画像の一例である。
【図７】この発明の実施例により得られる画像の一例である。
【図８】図６の画像の特徴パターンである。
【図９】図７の画像の特徴パターンである。
【図１０】着目画素と局所領域の説明図である。
【図１１】特徴パターンの比較の説明図である。
【図１２】濃度プロファイルを示す説明図である。
【図１３】正規化された濃度プロファイルを示す説明図である。
【図１４】この発明の表示例を示す説明図である。
【図１５】この発明の血管パターンを求める方法の説明図である。
【図１６】この発明の実施例の評価値を求める方法の説明図である。
【図１７】この発明の実施例の動作を示すフローチャートである。
【図１８】１回目に撮像した画像の特徴パターンである。
【図１９】２回目に撮像した画像の特徴パターンである。
【図２０】この発明の特徴パターンを一致させた図である。
【図２１】この発明の特徴パターンの移動方法の説明図である。
【図２２】この発明の特徴パターンの移動方法の説明図である。
【符号の説明】
１  検出部
２  データ処理部
３  解析部
９  操作部
１０  出力部
１１  光源部
１２  撮像部
１４  ヒトの指
【図１】
0004554739.tif 000004
【図２】
0004554739.tif 000005
【図３】
0004554739.tif 000006
【図４】
0004554739.tif 000007
【図５】
0004554739.tif 000008
【図６】
0004554739.tif 000009
【図７】
0004554739.tif 000010
【図８】
0004554739.tif 000011
【図９】
0004554739.tif 000012
【図１０】
0004554739.tif 000013
【図１１】
0004554739.tif 000014
【図１２】
0004554739.tif 000015
【図１３】
0004554739.tif 000016
【図１４】
0004554739.tif 000017
【図１５】
0004554739.tif 000018
【図１６】
0004554739.tif 000019
【図１７】
0004554739.tif 000020
【図１８】
0004554739.tif 000021
【図１９】
0004554739.tif 000022
【図２０】
0004554739.tif 000023
【図２１】
0004554739.tif 000024
【図２２】
0004554739.tif 000025