JP2017158718A - 評価装置、出力方法、並びにコンピュータプログラム - Google Patents

Abstract

【課題】相対的に簡易な方法で自律神経機能を評価する。【解決手段】評価装置（１０）は、被験者（９００）の血流量を取得する取得部（１２）と、血流量を変化させるための状態変化が被験者に生じた後の血流量の変化を示す時間波形が規定する所定領域の面積（Ｓ）に基づいて、被験者の生体機能を評価する評価部（１４）とを備える。【選択図】図２

Description

本発明は、例えば被験者の生体機能を評価する評価装置、被験者の生体機能の評価結果を出力する出力方法、及び、コンピュータを評価装置として機能させるコンピュータプログラムの技術分野に関する。

人（或いは、任意の生体）が急に起立した場合や起立した後に立ち続ける場合等には、一過性の血圧低下が起こる可能性がある。このような一過性の血圧低下は、健常者であれば、数秒から十数秒程度で解消される。しかしながら、自律神経機能が悪化している者であれば、自律神経機能の悪化に起因して血圧調整機能が低下している可能性が高い。このため、一過性の血圧低下が解消しにくく、血圧低下が大きくなった場合は、結果として、頭部血流の低下を招き、人はめまいを感じたり或いは失神したりしてしまいかねない。従って、人の血流状態（特に、頭部の血流状態）に基づいて自律神経機能を評価することが望ましい。

自律神経機能を評価する装置の一例は、特許文献１及び２に記載されている。具体的には、特許文献１及び２の夫々には、血流量を変化させる負荷を与える前の第１姿勢（しゃがみ姿勢）から血流量を変化させる負荷を与えた後の第２姿勢（立位姿勢）に移行する被験者の血流量の減少率及び回復時間に基づいて、被験者の自律神経機能を評価する装置が記載されている。

一方で、非特許文献１に記載されているように、起立負荷による血流量の反応型は、起立負荷によって減少した血流量が標準時間以内に回復する基本型、起立負荷によって減少した血流量の回復に標準時間以上要する回復遅延型、起立負荷によって減少した血流量が十分に回復しない不完全回復型及びその他の型という４群に分類される。血流量の反応型が不完全回復型となる場合には、血流量が十分に回復しないがゆえに、特許文献１及び２に記載された回復時間というパラメータを用いて自律神経機能を評価することが困難である。そこで、非特許文献１に記載された方法は、血流量の反応型が基本型又は回復遅延型となる場合には、被験者が起立することで血流量が減少してから回復するまでに要する反応時間に基づいて被験者の自律神経機能を評価し、血流量の反応型が不完全回復型となる場合には、起立前の血流量と起立後の血流量との比に基づいて被験者の自律神経機能を評価している。

国際公開第２０１３／０３８５５１号パンフレット 国際公開第２０１３／０３８５５２号パンフレット

畑中裕司、石川和夫、川崎富泰、窪田伸三、高木潔、丹家元陽、吉村幸男、老籾宗忠、馬場茂明、「レーザー皮膚血流量計による糖尿病性自律神経障害の定量的分析−とくに起立負荷テストについて−」、糖尿病、１９８４年、２７巻６号、ｐ６８９−６９５

しかしながら、非特許文献１に記載された手法は、被験者の血流量の反応型の違いに応じて選択される異なる複数のパラメータ（つまり、反応時間及び比）を用いて、自律神経機能を評価する必要がある。つまり、非特許文献１に記載された手法は、血流量の反応型の全ての種類に共通のパラメータを用いて自律神経機能を評価することができない。従って、血流量の反応型毎の評価精度がばらつくがゆえに、評価装置が自律神経機能を高精度に評価することができない可能性があるという技術的問題点が生ずる。

尚、このような技術的問題点は、血流量に基づいて自律神経機能を評価する装置（方法）のみならず、血流量に基づいて被験者の任意の生体機能を評価する装置（方法）においても生ずる可能性がある。

本発明が解決しようとする課題には上記のようなものが一例として挙げられる。本発明は、被験者の生体機能をより高精度に評価することが可能な評価装置、出力方法及びコンピュータプログラムを提供することを課題とする。

本発明の評価装置の第１の態様は、被験者の血流量を取得する取得部と、前記血流量を変化させるための状態変化が前記被験者に生じた後の前記血流量の変化を示す時間波形が規定する所定領域の面積に基づいて、前記被験者の生体機能を評価する評価部とを備える。

本発明の出力方法の第１の態様は、被験者の血流量を取得する取得工程と、前記血流量を変化させるための状態変化が前記被験者に生じた後の前記血流量の変化を示す時間波形が規定する所定領域の面積に基づく前記被験者の生体機能の評価結果を出力する出力工程とを備える。

本発明のコンピュータプログラムの第１の態様は、コンピュータを上述した本発明の評価装置の第１の態様として機能させる。

図１は、本実施例の評価システムの構成を示すブロック図である。 図２は、本実施例の評価装置の動作の全体の流れを示すフローチャートである。 図３（ａ）から図３（ｃ）は、夫々、血流情報が示す血流量の時間波形を示すグラフである。 図４は、評価パラメータの算出動作の流れを示すフローチャートである。 図５（ａ）から図５（ｃ）は、夫々、評価パラメータの算出動作で用いられる中間パラメータを、血流情報が示す血流量の時間波形上で示すグラフである。 図６（ａ）及び図６（ｂ）は、自律神経機能の悪化度合いと面積及び正規化面積との間の相関を示すグラフである。 図７（ａ）から図７（ｃ）は、夫々、血流情報が示す血流量の時間波形上で第１変形例の面積Ｓを示すグラフである。 図８（ａ）から図８（ｃ）は、夫々、血流情報が示す血流量の時間波形上で第２変形例の面積Ｓを示すグラフである。 図９（ａ）から図９（ｃ）は、夫々、血流情報が示す血流量の時間波形上で第３変形例の面積Ｓを示すグラフである。 図１０（ａ）から図１０（ｃ）は、夫々、血流情報が示す血流量の時間波形上で第３変形例の面積Ｓを示すグラフである。 図１１（ａ）から図１１（ｃ）は、夫々、血流情報が示す血流量の時間波形上で第３変形例の面積Ｓを示すグラフである。

以下、評価装置、出力方法、及び、コンピュータプログラムの実施形態について順に説明する。

（評価装置の実施形態）
＜１＞
本実施形態の評価装置は、被験者の血流量を取得する取得部と、前記血流量を変化させるための状態変化が前記被験者に生じた後の前記血流量の変化を示す時間波形が規定する所定領域の面積に基づいて、前記被験者の生体機能を評価する評価部とを備える。

本実施形態の評価装置によれば、被験者の血流量の反応型がどのような反応型であっても、「所定領域の面積」という共通のパラメータを用いて生体機能が評価される。従って、本実施形態の評価装置は、被験者の血流量の反応型の違いに応じて選択される異なる複数のパラメータを用いて生体機能を評価する比較例の評価装置と比較して、血流量の反応型毎の評価精度のばらつきを抑制することができる。従って、本実施形態の評価装置は、比較例の評価装置と比較して、生体機能をより高精度に評価することができる。

＜２＞
本実施形態の評価装置の他の態様では、前記評価部は、前記状態変化に起因した前記血流量の減少量を用いて前記面積を正規化することで得られる正規化結果に基づいて、前記生体機能を評価する。

この態様によれば、被験者によってばらつきが生じる血流量の減少量で面積が正規化されるので、正規化前の面積の被験者毎のばらつきに比べ、正規化後の面積の被験者毎のばらつき（誤差）を抑えられ、生体機能の悪化度合いに応じた面積の変化を捉えやすくなる。また、状態変化に起因した血流量の減少量が被験者毎に異なる場合であっても、評価装置は、生体機能を高精度に評価することができる。

＜３＞
本実施形態の評価装置の他の態様では、前記所定領域は、前記時間波形のうち前記状態変化が生じた第１時刻から前記第１時刻よりも後の第２時刻までの間の第１期間中の前記血流量を示す少なくとも一部の波形部分が規定する領域を含む。

この態様によれば、評価装置は、このような所定領域の面積に基づいて生体機能を高精度に評価することができる。

＜４＞
本実施形態の評価装置の他の態様では、前記所定領域は、前記時間波形のうち前記状態変化が生じた第１時刻よりも後の第３時刻から前記第３時刻よりも後の第２時刻までの間の第２期間中の前記血流量を示す少なくとも一部の波形部分が規定する領域を含む。

この態様によれば、評価装置は、このような所定領域の面積に基づいて生体機能を高精度に評価することができる。

＜５＞
本実施形態の評価装置の他の態様では、前記所定領域は、前記時間波形のうち前記状態変化が生じた第１時刻から前記第１時刻よりも後の第３時刻までの間の第３期間中の前記血流量を示す少なくとも一部の波形部分が規定する領域、及び、前記時間波形のうち前記第３時刻から前記第３時刻よりも後の第２時刻までの間の第２期間中の前記血流量を示す少なくとも一部の波形部分が規定する領域を含む。

この態様によれば、評価装置は、このような所定領域の面積に基づいて生体機能を高精度に評価することができる。

＜６＞
上述したように第３時刻に基づいて定まる所定領域の面積に基づいて生体機能を評価する評価装置の他の態様では、前記第３時刻は、前記状態変化に起因して前記血流量が最も減少した時刻である。

この態様によれば、評価装置は、このような第３時刻に基づいて定まる所定領域の面積に基づいて生体機能を高精度に評価することができる。

＜７＞
上述したように第２時刻に基づいて定まる所定領域の面積に基づいて生体機能を評価する評価装置の他の態様では、前記第２時刻は、前記状態変化に起因して減少した前記血流量が所定量にまで回復した回復時刻である
この態様によれば、評価装置は、回復時刻に基づいて定まる所定領域の面積に基づいて生体機能を高精度に評価することができる。

＜８＞
上述したように第２時刻に基づいて定まる所定領域の面積に基づいて生体機能を評価する評価装置の他の態様では、前記状態変化に起因して減少した前記血流量が所定量に回復する回復時刻が、所定のカットオフ時刻よりも前である場合には、前記第２時刻は前記回復時刻であり、前記回復時刻が前記カットオフ時刻よりも後である場合には、前記第２時刻は前記カットオフ時刻である。

この態様によれば、評価装置は、回復時刻又はカットオフ時刻に基づいて定まる所定領域の面積に基づいて生体機能を高精度に評価することができる。

＜９＞
上述したように第２時刻がカットオフ時刻となる評価装置の他の態様では、前記カットオフ時刻は、前記生体機能が正常でない前記被験者に前記状態変化が生じた場合における前記回復時刻に基づいて規定される。

この態様によれば、評価装置は、このようなカットオフ時刻に基づいて定まる所定領域の面積に基づいて生体機能を高精度に評価することができる。

＜１０＞
上述したように第２時刻が回復時刻となる評価装置の他の態様では、前記所定量は、前記状態変化が生ずる前の第４期間中の前記血流量の平均値である。

この態様によれば、評価装置は、状態変化に起因して減少した血流量が平均値にまで回復した回復時刻に基づいて定まる所定領域の面積に基づいて生体機能を高精度に評価することができる。

＜１１＞
上述したように状態変化に起因して減少した血流量が平均値にまで回復した時刻が回復時刻となる評価装置の他の態様では、前記少なくとも一部の波形が規定する前記領域は、前記平均値と前記少なくとも一部の波形部分とが囲む領域を含む。

この態様によれば、評価装置は、このような所定領域の面積に基づいて生体機能を高精度に評価することができる。

＜１２＞
上述したように第２時刻が回復時刻となる評価装置の他の態様では、前記所定量は、前記状態変化が生ずる前の第４期間中の前記血流量の平均値よりも小さい中間値である。

この態様によれば、評価装置は、状態変化に起因して減少した血流量が中間値にまで回復した回復時刻に基づいて定まる所定領域の面積に基づいて生体機能を高精度に評価することができる。

＜１３＞
上述したように状態変化に起因して減少した血流量が中間値にまで回復した時刻が回復時刻となる評価装置の他の態様では、前記少なくとも一部の波形が規定する前記領域は、前記中間値と前記少なくとも一部の波形部分とが囲む領域を含む。

この態様によれば、評価装置は、このような所定領域の面積に基づいて生体機能を高精度に評価することができる。

（出力方法の実施形態）
＜１４＞
本実施形態の出力方法は、被験者の血流量を取得する取得工程と、前記血流量を変化させるための状態変化が前記被験者に生じた後の前記血流量の変化を示す時間波形が規定する所定領域の面積に基づく前記被験者の生体機能の評価結果を出力する出力工程とを備える。

本実施形態の出力方法によれば、上述した本実施形態の評価装置が享受することができる各種効果と同様の効果を適切に享受することができる。尚、上述した本実施形態の評価装置における各種態様に対応して、本実施形態の出力方法も各種態様を採ることが可能である。

（コンピュータプログラムの実施形態）
＜１５＞
本実施形態のコンピュータプログラムは、コンピュータを上述した本実施形態の評価装置（但し、その各種態様を含む）として機能させる
本実施形態のコンピュータプログラムによれば、上述した本実施形態の評価装置が享受することができる各種効果と同様の効果を適切に享受することができる。尚、上述した本実施形態の評価装置における各種態様に対応して、本実施形態のコンピュータプログラムも各種態様を採ることが可能である。また、本実施形態のコンピュータプログラムは、例えばコンピュータ読取可能な記録媒体（例えば、ＣＤやＤＶＤやＢＤ等の光ディスクや、半導体メモリや、ハードディスク等の磁気ディスク）に記録されていてもよい。

本実施形態のこのような作用及び他の利得は次に説明する実施例から更に明らかにされる。

以上説明したように、本実施形態の評価装置は、取得手段と、評価手段とを備える。本実施形態の出力方法は、取得工程と、出力工程とを備える。本実施形態のコンピュータプログラムは、コンピュータを本実施形態の評価装置として機能させる。従って、被験者の生体機能がより高精度に評価される。

以下、図面を参照しながら、評価装置、出力方法及びコンピュータプログラムの実施例について説明する。尚、以下では、本発明の評価装置、出力方法及びコンピュータプログラムが、被験者９００の自律神経機能を評価する評価装置１０を備える評価システム１に適用された実施例について説明を進める。但し、本発明の評価装置、出力方法及びコンピュータプログラムは、自律神経機能とは異なる被験者９００の生体機能を評価する任意の評価装置に適用されてもよい。

（１）評価システム１の構成
図１を参照して、本実施例の評価システム１の構成について説明する。図１は、本実施例の評価システム１の構成を示すブロック図である。

図１に示すように、評価システム１は、人間等の被験者９００の自律神経機能を評価するための評価試験を行うシステムである。評価システム１は、評価装置１０と、血流計２０と、ディスプレイ３０とを備える。

評価装置１０は、制御部１１と、取得部１２と、指示部１３と、評価部１４とを備えている。尚、制御部１１、取得部１２、指示部１３及び評価部１４は、例えば専用のＩＣチップ等によって物理的に実現されてもよい。或いは、制御部１１、取得部１２、指示部１３及び評価部１４は、例えば、ＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）上で動作するソフトウェアによって論理的に実現されてもよい。

制御部１１は評価装置１０全体の動作を制御する。

取得部１２は、血流計２０が出力する被験者９００の血流情報を取得する。血流情報は、被験者９００の血流量を示す情報である。例えば、血流情報は、血流量の時間波形（つまり、血流量の推移（言い換えれば、血流量の時間的な変化）を示す血流波形）を示す情報を含んでいる。

指示部１３は、評価装置１０によって自律神経機能が評価される被験者９００に対して、所望の指示を出力する。本実施例では、指示部１３は、所望の指示を、ディスプレイ３０上に表示される画像（例えば、所望の指示を文字や図形等で示す画像）として出力する。但し、指示部１３は、所望の指示を、その他の態様で出力してもよい。例えば、指示部１３は、所望の指示を、スピーカから出力される音声（例えば、所望の指示を通知する音声等）として出力してもよい。

評価部１４は、取得部１２が取得する血流情報に基づいて、被験者９００の自律神経機能を評価する。評価部１４は、その内部に物理的に又は論理的に実現される処理ブロックとして、パラメータ演算部１４１と、評価結果出力部１４２とを備える。

パラメータ演算部１４１は、取得部１２が取得する血流情報に基づいて、自律神経機能を評価するために用いられる評価パラメータを算出する。

評価結果出力部１４２は、パラメータ演算部１４１が算出した評価パラメータに基づいて、被験者９００の自律神経機能を評価する。加えて、評価結果出力部１４２は、被験者９００の自律神経機能の評価結果を出力する。本実施例では、評価結果出力部１４２は、自律神経機能の評価結果を、ディスプレイ３０上に表示される画像（例えば、自律神経機能の評価の結果を文字や図形等で示す画像）として出力する。但し、評価結果出力部１４２は、自律神経機能の評価結果を、その他の態様（例えば、音声等）で出力してもよい。

血流計２０は、被験者９００の血流量を検出する。血流計２０は、被験者９００の頭部の血流量（或いは、頭部の近傍にある血管の血流量）を検出することが好ましい。血流計２０は、被験者９００に装着される。血流計２０は、当該検出した血流量を示す血流情報を評価装置１０（特に、取得部１２）に出力する。血流計２０としては、例えばレーザードップラーフローメトリー法を用いて血流波形を検出するレーザー血流計が一例としてあげられる。

ディスプレイ３０は、所望の画像を表示する装置である。例えば、本実施例では、ディスプレイ３０は、評価システム１の動作状態を示すための画像（言い換えれば、画面）や、指示部１３から出力される所望の指示を示す画像（画面）や、評価結果出力部１４２から出力される自律神経機能の評価の結果を示す画像（画面）等を表示する。

（２）評価装置１０の動作
続いて、図２から図６を参照しながら、本実施例の評価装置１０の動作について説明する。

（２−１）評価装置１０の動作の全体の流れ
初めに、図２を参照して、本実施例の評価装置１０の動作の全体の流れについて説明する。図２は、本実施例の評価装置１０の動作の全体の流れを示すフローチャートである。

図２に示すように、自律神経機能の評価試験が開始される前であっても、血流計２０が血流情報を出力している場合には、取得部１２は、血流情報を取得してもよい（ステップＳ１１）。評価試験が開始される前に取得された血流情報は、後に詳述する評価部１４が評価パラメータを算出する際に評価部１４によって参照されてもよい。尚、ステップＳ１１の動作は、必ずしも行われなくともよい。

その後、制御部１１は、評価試験の開始を指示する操作が、評価システム１のユーザ（例えば、医師等の医療従事者又は被験者９００等）によって行われたか否かを判定する（ステップＳ１２）。

ステップＳ１２の判定の結果、評価試験の開始を指示する操作が行われていないと判定される場合には（ステップＳ１２：Ｎｏ）、ステップＳ１１の動作が繰り返される。つまり、評価試験の開始の指示が行われるまでは、評価装置１０は実質的に待機することになる。

他方で、ステップＳ１２の判定の結果、評価試験の開始を指示する操作が行われたと判定される場合には（ステップＳ１２：Ｙｅｓ）、続いて、指示部１３は、座位開始時刻Ｔａに座位姿勢をとると共にその後第１維持時間だけ（或いは、第１維持時間以上）座位姿勢を維持し続けることを被験者９００に指示するための姿勢統制指示を出力する（ステップＳ１３）。その結果、被験者９００は、姿勢統制指示に従って、座位開始時刻Ｔａに座位姿勢をとると共に、座位姿勢を第１維持時間以上維持する。尚、「座位姿勢」の一例として、椅子等の器具に被験者９００が座っている姿勢があげられる。「座位姿勢」の他の一例として、被験者９００の体重を折り曲げた下肢で支えるように被験者９００が座っている姿勢（言い換えれば、しゃがみこんでいる姿勢）があげられる。「座位姿勢」の他の一例として、下肢を折り曲げた状態で椅子等の器具に被験者９００が座っている姿勢があげられる。

その後、制御部１１は、座位開始時刻Ｔａを起点として、第１維持時間が経過したか否かを判定する（ステップＳ１４）。つまり、制御部１１は、被験者９００が第１維持時間以上座位姿勢を維持し続けたか否かを判定する。

ここで、後に詳述するように、本実施例では、被験者９００が座位姿勢から立位姿勢に移行した後の血流量の変化の態様に基づいて、被験者９００の自律神経機能が評価される。従って、高精度な（言い換えれば、信頼性の高い）評価を行うという観点からは、立位姿勢に移行する直前の被験者９００の血流量は、十分に安定している（言い換えれば、変動が少ない）ことが好ましい。つまり、立位姿勢に移行する直前の被験者９００の血流量には、座位姿勢を取る前の被験者９００の活動状況（例えば、運動状況）に起因したノイズ成分が重畳されていないことが好ましい。このため、第１維持時間は、被験者９００の血流量が十分に安定するようになるまでに必要な時間に応じて設定されることが好ましい。あるいは、立位姿勢に移行する直前の座位姿勢を、足を台座に乗せる姿勢やしゃがむ姿勢とする場合は、被験者９００は、その姿勢をとる以前に座位姿勢にて所定時間の安静を行った後に足を台座に乗せる姿勢やしゃがむ姿勢をとり、立位姿勢に移行することが好ましい。

ステップＳ１４の判定の結果、第１維持時間が経過していないと判定される場合には（ステップＳ１４：Ｎｏ）、評価装置１０は、第１維持時間が経過するまで待機する。従って、被験者９００は、座位姿勢を維持し続ける。

他方で、ステップＳ１４の判定の結果、第１維持時間が経過したと判定される場合には（ステップＳ１４：Ｙｅｓ）、指示部１３は、立位開始時刻Ｔｓに座位姿勢から立位姿勢に移行すると共にその後第２維持時間だけ（或いは、第２維持時間以上）立位姿勢を維持し続けることを被験者９００に指示するための姿勢統制指示を出力する（ステップＳ１５）。その結果、被験者９００は、姿勢統制指示に従って、立位開始時刻Ｔｓに座位姿勢から立位姿勢に移行すると共に、立位姿勢を第２維持時間以上維持する。尚、「立位姿勢」の一例として、被験者９００が下肢を概ね伸ばした状態で立っている姿勢又は立ち上がった姿勢（つまり、下肢を概ね延ばした状態で被験者９００の体重を足で支えている姿勢）があげられる。

このとき、被験者９００は、被験者９００自身の筋肉を使用することで立ち上がる能動起立を行うことで、座位姿勢から立位姿勢に移行してもよい。或いは、被験者９００は、例えば他者の力を借りて立ち上がる受動起立（例えば、電動ベッドに寝ている被験者９００が、電動ベッドが自動的に起き上がる力を借りて立ち上がる受動起立やヘッドアップチルトを用いた受動起立等）を行うことで、座位姿勢から立位姿勢に移行してもよい。

その後、制御部１１は、立位開始時刻Ｔｓを起点として、第２維持時間が経過したか否かを判定する（ステップＳ１６）。つまり、制御部１１は、被験者９００が第２維持時間以上立位姿勢を維持し続けたか否かを判定する。

ここで、後に詳述するように、本実施例では、被験者９００が立位姿勢に移行することで減少した血流量が回復する態様に基づいて、被験者９００の自律神経機能が評価される。具体的には、立位開始時刻Ｔｓから被験者９００が立位姿勢に移行することで減少した血流量が立位姿勢に移行する前の所定期間中の血流量の平均値Ｆａｖｅ（或いは、任意の所定量）にまで回復した回復時刻Ｔｒ（或いは、所定のカットオフ時刻Ｔｃ）までの間の血流量の変化の態様に基づいて、被験者９００の自律神経機能が評価される。従って、第２維持時間は、立位開始時刻Ｔｓから回復時刻Ｔｒ（或いは、カットオフ時刻Ｔｃ）までの間の時間よりも長いことが好ましい。

ステップＳ１６の判定の結果、第２維持時間が経過していないと判定される場合には（ステップＳ１６：Ｎｏ）、評価装置１０は、第２維持時間が経過するまで待機する。従って、被験者９００は、立位姿勢を維持し続ける。

他方で、ステップＳ１６の判定の結果、第２維持時間が経過したと判定される場合には（ステップＳ１６：Ｙｅｓ）、パラメータ演算部１４１は、ステップＳ１３からステップＳ１６の動作が行われている間に取得部１２によって取得された血流情報に基づいて、自律神経機能を評価するために用いられる評価パラメータを算出する（ステップＳ１７）。尚、ステップＳ１７においてパラメータ演算部１４１によって評価パラメータが算出されることを考慮すれば、図２では明示的に記載していないものの、取得部１２は、ステップＳ１３からステップＳ１６の動作が行われている間も血流情報を取得する。

ここで、図３（ａ）から図３（ｃ）を参照しながら、ステップＳ１３からステップＳ１６の動作が行われている間に取得部１２によって取得される血流情報が示す血流量の一例について説明する。図３（ａ）から図３（ｃ）は、夫々、血流情報が示す血流量の時間波形を示すグラフである。

図３（ａ）は、被験者９００が立位姿勢に移行することで減少した血流量が立位姿勢に移行する前の所定期間中の血流量の平均値Ｆａｖｅにまで回復するために要する時間が相対的に短い場合に取得される血流情報（血流量の時間波形）を示す。このような血流情報は、被験者９００が健常者である（つまり、被験者９００の自律神経機能が悪化していない）場合に取得される。図３（ａ）に示すように、被験者９００が座位姿勢から立位姿勢に移行すると、血流量は減少する。なぜならば、立位姿勢への移行時に重力による血液の下方シフトが生じるからであり、また、被験者９００が座位姿勢をとっている間に生じていた膝関節の曲げ等による血管の圧迫が、被験者９００が立位姿勢に移行することで解消されるからである。従って、被験者９００の姿勢が座位姿勢から立位姿勢に移行するという動作（状態変化）は、被験者９００の血流量を変化させるための状態変化の一例に相当する。その後、被験者９００の血圧調整機能に起因して、血流量は再び増加していく。その結果、血流量は、平均値Ｆａｖｅにまで回復する。図３（ａ）に示す血流量の時間波形は、通常回復型の時間波形に分類される。尚、被験者９００が健常者である場合には、立位開始時刻Ｔｓから回復時刻Ｔｒまでの間の期間の長さは、概ね十数秒程度になる。

尚、図３（ａ）では、説明の便宜上、立位姿勢に移行する前の所定期間中の血流量の平均値Ｆａｖｅが、被験者９００が座位姿勢をとっている所定期間中の血流量の平均値Ｆａｖｅであるものとする。以下の説明においても同様である。

図３（ｂ）は、被験者９００が立位姿勢に移行することで減少した血流量が平均値Ｆａｖｅにまで回復するために要する時間が相対的に長い場合に取得される血流情報（血流量の時間波形）を示す。このような血流情報は、被験者９００の自律神経機能が悪化している場合に取得される。図３（ｂ）に示すように、被験者９００が座位姿勢から立位姿勢に移行すると、血流量は減少する。その後、被験者９００の血圧調整機能に起因して、血流量は再び増加していく。しかしながら、自律神経機能の悪化に起因して血圧調整機能もまた悪化しているがゆえに、血流量が増加していく速度は、被験者９００が健常者である場合に血流量が増加していく速度よりも遅い。従って、被験者９００の自律神経機能が悪化している場合に血流量が平均値Ｆａｖｅにまで回復するために要する時間は、被験者９００が健常者である場合に血流量が平均値Ｆａｖｅにまで回復するために要する時間よりも長くなる。図３（ｂ）に示す血流量の時間波形は、回復遅延型の時間波形に分類される。

図３（ｃ）は、被験者９００が立位姿勢に移行することで減少した血流量が所定時間（例えば、３０秒から６０秒の範囲内の任意の時間）以内に平均値Ｆａｖｅにまで回復しない場合に取得される血流情報（血流量の時間波形）を示す。このような血流情報は、被験者９００の自律神経機能が悪化している場合に取得される。図３（ｃ）に示すように、被験者９００が座位姿勢から立位姿勢に移行すると、血流量は減少する。その後、被験者９００の血圧調整機能に起因して、血流量は再び増加していく。しかしながら、自律神経機能の悪化に起因して血圧調整機能もまた悪化しているがゆえに、血流量が増加していく速度は、被験者９００が健常者である場合に血流量が増加していく速度よりも遅い。その結果、図３（ｃ）に示す例では、血流量は、立位開始時刻Ｔｓから所定時間以内に（図３（ｃ）に示す例では、カットオフ時刻Ｔｃ以前に）平均値Ｆａｖｅにまで回復することはない。図３（ｃ）に示す血流量の時間波形は、不完全回復型の時間波形に分類される。

尚、カットオフ時刻Ｔｃは、自律神経機能が悪化している被験者９００が立位姿勢に移行してから当該被験者９００の血流量が平均値Ｆａｖｅにまで回復するために要した時間に基づいて設定されることが好ましい。例えば、自律神経機能が悪化している複数の被験者９００を対象に立位姿勢に移行してから血流量が平均値Ｆａｖｅにまで回復するために要した時間が計測され、当該計測された時間の平均値が算出され、立位開始時刻Ｔｓから当該算出された平均値に相当する時間だけ進んだ時刻がカットオフ時刻Ｔｃに設定されてもよい。或いは、例えば、立位開始時刻Ｔｓから当該算出された平均値に対して所定のマージン（例えば、当該算出された平均値の標準偏差に基づいて算出されるマージン）を加算若しくは減算することで得られる値に相当する時間だけ進んだ時刻がカットオフ時刻Ｔｃに設定されてもよい。或いは、例えば、立位開始時刻Ｔｓから当該算出された平均値に対して所定の係数を掛け合わせることで得られる値に相当する時間だけ進んだ時刻がカットオフ時刻Ｔｃに設定されてもよい。

再び図２において、評価結果出力部１４２は、ステップＳ１７で算出された評価パラメータに基づいて、被験者９００の自律神経機能を評価する（ステップＳ１８）。例えば、評価結果出力部１４２は、評価パラメータが閾値ＴＨ１未満である場合には、被験者９００の自律神経機能が「優（めまい度：低）」であると判定してもよい。例えば、評価結果出力部１４２は、評価パラメータが閾値ＴＨ１以上であり且つ閾値ＴＨ２（但し、ＴＨ２＞ＴＨ１）未満である場合には、被験者９００の自律神経機能が「良（めまい度：中）」であると判定してもよい。例えば、評価結果出力部１４２は、評価パラメータが閾値ＴＨ２以上である場合には、被験者９００の自律神経機能が「悪（めまい度：高）」であると判定してもよい。

その後、評価結果出力部１４２は、自律神経機能の評価結果を出力する（ステップＳ１８）。例えば、評価結果出力部１４２は、自律神経機能の評価結果を、ディスプレイ３０上に表示される画像（例えば、自律神経機能の評価結果を文字や図形等で示す画像）として出力する。但し、評価結果出力部１４２は、自律神経機能の評価結果を、必ずしも出力しなくともよい。

その後、制御部１１は、再度の評価試験（つまり、再試験）の開始を指示する操作が、評価システム１のユーザによって行われたか否かを判定する（ステップＳ１９）。ステップＳ１９の判定の結果、再度の評価試験の開始を指示する操作が行われたと判定される場合には（ステップＳ１９：Ｙｅｓ）、ステップＳ１３以降の動作が繰り返される。他方で、ステップＳ１９の判定の結果、再度の評価試験の開始を指示する操作が行われていないと判定される場合には（ステップＳ１９：Ｎｏ）、評価システム１の動作が終了する。

（２−２）評価パラメータの算出動作の流れ
続いて、図４及び図５（ａ）から図５（ｃ）を参照しながら、図２のステップＳ１７で行われる評価パラメータの算出動作について説明する。図４は、評価パラメータの算出動作の流れを示すフローチャートである。図５（ａ）から図５（ｃ）は、夫々、評価パラメータの算出動作で用いられる中間パラメータを、血流情報が示す血流量の時間波形上で示すグラフである。

図４に示すように、パラメータ演算部１４１は、図２のステップＳ１３からステップＳ１６の動作が行われている間に取得部１２が取得した血流情報に基づいて、以下に説明するステップＳ１７２からステップＳ１７９に示す動作を実行する。

具体的には、パラメータ演算部１４１は、被験者９００が立位姿勢に移行する前の所定期間中の血流量の平均値Ｆａｖｅ（図５（ａ）から図５（ｃ）参照）を算出する（ステップＳ１７２）。尚、上述したように、説明の便宜上、平均値Ｆａｖｅは、被験者９００が座位姿勢をとっている所定期間中の血流量の平均値Ｆａｖｅであるものとする。

パラメータ演算部１４１は、更に、被験者９００が立位姿勢に移行した後の血流量の最小値Ｆｍｉｎ（図５（ａ）から図５（ｃ）参照）を算出する（ステップＳ１７３）。更に、パラメータ演算部１４１は、血流量が最小値Ｆｍｉｎとなった最小時刻Ｔｍｉｎ（図５（ａ）から図５（ｃ）参照）を算出する（ステップＳ１７３）。例えば、パラメータ演算部１４１は、血流量の微分結果に基づいて、最小値Ｆｍｉｎ及び最小時刻Ｔｍｉｎを算出してもよい。具体的には、パラメータ演算部１４１は、血流量の微分結果がマイナスからプラスに変化する時刻を算出すると共に、当該算出した時刻を最小時刻Ｔｍｉｎに設定してもよい。更に、パラメータ演算部１４１は、算出した最小時刻Ｔｍｉｎにおける血流量を、最小値Ｆｍｉｎに設定してもよい。

その後、パラメータ演算部１４１は、ステップＳ１７２で算出した平均値Ｆａｖｅ及びステップＳ１７３で算出した最小値Ｆｍｉｎに基づいて、被験者９００が立位姿勢に移行したことに起因した血流量の減少量ΔＦ（図５（ａ）から図５（ｃ）参照）を算出する（ステップＳ１７４）。具体的には、パラメータ演算部１４１は、減少量ΔＦ＝平均値Ｆａｖｅ−最小値Ｆｍｉｎという数式を用いて、減少量ΔＦを算出する。

パラメータ演算部１４１は、更に、被験者９００が立位姿勢に移行したことに起因して減少した血流量が平均値Ｆａｖｅにまで回復した回復時刻Ｔｒ（図５（ａ）から図５（ｂ）参照）を算出する（ステップＳ１７５）。例えば、パラメータ演算部１４１は、立位開始時刻Ｔｓよりも後において血流量が平均値Ｆａｖｅと一致した時刻を算出すると共に、当該算出した時刻を回復時刻Ｔｒに設定してもよい。

尚、血流量の時間波形が不完全回復型の時間波形である場合には、被験者９００が立位姿勢に移行することで減少した血流量が平均値Ｆａｖｅにまで回復しない可能性がある（図５（ｃ）参照）。この場合には、パラメータ演算部１４１は、回復時刻Ｔｒを算出することができない。従って、パラメータ演算部１４１は、立位開始時刻Ｔｓよりも後において血流量が平均値Ｆａｖｅと一致しない場合には、回復時刻Ｔｒを無限大（或いは、カットオフ時刻Ｔｃよりも大きい任意の値）に設定してもよい。但し、血流量の時間波形が不完全回復型の時間波形である場合であっても、被験者９００が立位姿勢に移行することで減少した血流量が、カットオフ時刻Ｔｃ以降に平均値Ｆａｖｅにまで回復する場合もある。この場合には、パラメータ演算部１４１は、実際の回復時刻Ｔｒを算出してもよい。但し、この場合の回復時刻Ｔｒは、当然にカットオフ時刻Ｔｃよりも大きくなる。

その後、パラメータ演算部１４１は、回復時刻Ｔｒがカットオフ時刻Ｔｃよりも早いか否か（つまり、ステップＳ１７５で算出した回復時刻Ｔｒがカットオフ時刻Ｔｃ未満であるか否か）を判定する（ステップＳ１７６）。

ステップＳ１７６の判定の結果、回復時刻Ｔｒがカットオフ時刻Ｔｃよりも早いと判定される場合には（ステップＳ１７６：Ｙｅｓ）、血流量の時間波形が通常回復型（図５（ａ）参照）又は回復遅延型（図５（ｂ）参照）の時間波形であると推定される。この場合には、パラメータ演算部１４１は、血流量の時間波形のうち立位開始時刻Ｔｓから回復時刻Ｔｒまでの間の期間中の血流量の変化（推移）を示す波形部分と平均値Ｆａｖｅとが囲む所定領域（図５（ａ）又は図５（ｂ）の網掛け部分）の面積Ｓを算出する（ステップＳ１７７）。言い換えれば、血流量をｙとし、時間をｔとし、血流波形がｙ＝ｆ（ｔ）という関数で表される場合には、パラメータ演算部１４１は、ｙ＝ｆ（ｔ）というグラフと、ｙ＝Ｆａｖｅというグラフと、ｔ＝Ｔｓというグラフと、ｔ＝Ｔｒというグラフとが囲む所定領域の面積を算出する。

パラメータ演算部１４１は、例えば、立位開始時刻Ｔｓから回復時刻Ｔｒまでの間の期間における血流量と平均値Ｆａｖｅとの差分の積分値を算出すると共に、当該積分値を面積Ｓとして取り扱ってもよい。より具体的には、パラメータ演算部１４１は、数式１を用いて面積Ｓを算出してもよい。

他方で、ステップＳ１７６の判定の結果、回復時刻Ｔｒがカットオフ時刻Ｔｃよりも早くないと判定される場合には（ステップＳ１７６：Ｎｏ）、血流量の時間波形が不完全回復型（図５（ｃ）参照）の時間波形であると推定される。この場合には、パラメータ演算部１４１は、血流量の時間波形のうち立位開始時刻Ｔｓからカットオフ時刻Ｔｃまでの間の期間中の血流量の変化（推移）を示す波形部分と平均値Ｆａｖｅとが囲む所定領域（図５（ｃ）の網掛け部分）の面積Ｓを算出する（ステップＳ１７８）。言い換えれば、パラメータ演算部１４１は、ｙ＝ｆ（ｔ）というグラフと、ｙ＝Ｆａｖｅというグラフと、ｔ＝Ｔｓというグラフと、ｔ＝Ｔｃというグラフとが囲む所定領域の面積を算出する。

パラメータ演算部１４１は、例えば、立位開始時刻Ｔｓからカットオフ時刻Ｔｃまでの間の期間における血流量と平均値Ｆａｖｅとの差分の積分値を算出すると共に、当該積分値を面積Ｓとして取り扱ってもよい。より具体的には、パラメータ演算部１４１は、数式２を用いて面積Ｓを算出してもよい。

その後、パラメータ演算部１４１は、ステップＳ１７４で算出した減少量ΔＦを用いてステップＳ１７７又はステップＳ１７８で算出した面積Ｓを正規化することで、正規化面積ＲＡを算出する（ステップＳ１７９）。具体的には、パラメータ演算部１４１は、正規化面積ＲＡ＝面積Ｓ／減少量ΔＦという数式を用いて、面積Ｓを正規化する（つまり、正規化面積ＲＡを算出する）。

但し、パラメータ演算部１４１は、平均値Ｆａｖｅと立位開始時刻Ｔｓから回復時刻Ｔｒ（或いは、カットオフ時刻Ｔｃ）までの間の期間における血流量との差分を波形として示す差分血流波形に基づいて、面積Ｓ（正規化面積ＲＡ）を算出してもよい。この場合、面積Ｓは、差分血流波形と時間軸（図５（ａ）から図５（ｃ）の夫々のグラフにおける横軸）とによって囲まれる領域の面積に相当する。更に、回復時刻Ｔｒは、立位開始時刻Ｔｓ以降に差分血流波形と時間軸とが交わる時刻に相当する。

その後、パラメータ演算部１４１は、ステップＳ１７９で算出した正規化面積ＲＡを、評価パラメータとして評価結果出力部１４２に出力する。その結果、評価結果出力部１４２は、正規化面積ＲＡに基づいて、自律神経機能を評価する。例えば、図５（ａ）及び図５（ｂ）に示すように、血流量の時間波形が通常回復型である場合の正規化面積ＲＡは、血流量の時間波形が回復遅延型である場合の正規化面積ＲＡよりも小さくなる。例えば、図５（ｂ）及び図５（ｃ）に示すように、血流量の時間波形が回復遅延型である場合の正規化面積ＲＡは、血流量の時間波形が不完全回復型である場合の正規化面積ＲＡよりも小さくなる。従って、正規化面積ＲＡと自律神経機能の良否との間には相関がある。従って、評価結果出力部１４２は、正規化面積ＲＡに基づいて、自律神経機能を適切に評価することができる。

（３）技術的効果
続いて、本実施例の評価装置１０によって得られる技術的効果について説明する。

上述したように、評価装置１０は、被験者９００の血流量の時間波形がどのような反応型であっても、「正規化面積ＲＡ」という共通の評価パラメータを用いて、自律神経機能を評価することができる。従って、評価装置１０は、血流量の時間波形の反応型の違いに応じて選択される異なる複数のパラメータを用いて自律神経機能を評価する比較例の評価装置と比較して、反応型毎の評価精度のばらつきを抑制することができる。つまり、血流量の時間波形が通常回復型又は回復遅延型である場合の評価装置１０の評価精度と、血流量の時間波形が不完全回復型である場合の評価装置１０の評価精度とは一致する（或いは、大きく異なることはない）。従って、評価装置１０は、比較例の評価装置と比較して、自律神経機能をより高精度に評価することができる。

加えて、評価装置１０は、面積Ｓを評価パラメータとして用いることに代えて、正規化面積ＲＡを評価パラメータとして用いている。このため、評価装置１０は、被験者９００毎に減少量ΔＦが異なる場合であっても、自律神経機能を適切に評価することができる。以下、図６（ａ）から図６（ｂ）を参照しながら、正規化面積ＲＡを評価パラメータとして用いることで得られる技術的効果について更に説明する。図６（ａ）及び図６（ｂ）は、自律神経機能の悪化度合いと面積Ｓ及び正規化面積ＲＡとの間の相関を示すグラフである。尚、図６（ａ）及び図６（ｂ）に示す相関は、本願発明者等の実験によって得られた相関であり、ここでの「自律神経機能の悪化度合い」とは、より正確には本願発明者らの実験における糖尿病患者の神経障害の悪化度合いの評価結果に基づくものである。

図６（ａ）は、自律神経機能の悪化度合いと面積Ｓとの間の相関を示している。図６（ａ）に示すように、自律神経機能の悪化度合いと面積Ｓとの間には、一定の相関が存在する。しかしながら、図６（ａ）に示すように、自律神経機能が重度に悪化している場合の面積Ｓは、理論的には自律神経機能が中度に悪化している場合の面積Ｓよりも大きくなるはずであるものの、実際には、自律神経機能が中度に悪化している場合の面積Ｓより小さくなってしまっている。この理由は、以下のとおりであると推定される。まず、面積Ｓは、血流量の反応型を決定付ける立位開始時刻Ｔｓから回復時刻Ｔｒまでの期間の長さのみならず、減少量ΔＦにも依存して変動する。従って、被験者９００が健常者である場合（つまり、立位開始時刻Ｔｓから回復時刻Ｔｒまでの期間の長さが十数秒となる場合）であっても、当該被験者９００の減少量ΔＦが相対的に大きい場合には、面積Ｓが相対的に大きくなる。従って、健常者であるはずの被験者９００の血流量の反応型が、回復遅延型又は不完全回復型であると誤って判定される可能性がある。同様に、被験者９００が健常者でない場合であっても、当該被験者９００の減少量ΔＦが相対的に小さい場合には、面積Ｓが相対的に小さくなる。従って、健常者でないはずの被験者９００の血流量の反応型が、通常回復型であると誤って判定される可能性がある。つまり、面積Ｓが減少量ΔＦに依存して変動するがゆえに、自律神経機能の悪化度合いと面積Ｓとの間の相関の精度が悪化する場合がある。

一方で、図６（ｂ）は、自律神経機能の悪化度合いと正規化面積ＲＡとの間の相関を示している。図６（ｂ）に示すように、自律神経機能の悪化度合いと正規化面積ＲＡとの間には、自律神経機能の悪化度合いと面積Ｓとの間の相関よりも強い相関が存在する。これは、面積Ｓを減少量ΔＦで正規化することで正規化面積ＲＡが得られるがゆえに、正規化面積ＲＡが減少量ΔＦに依存して変動しにくくなるからである。従って、評価装置１０は、このような正規化面積ＲＡに基づいて自律神経機能を評価することで、被験者９００の個体差の影響（特に、減少量ΔＦに関連する個体差の影響）を受けることなく自律神経機能を適切に評価することができる。

特に、図６（ａ）及び図６（ｂ）は、夫々、自律神経機能の悪化度合い毎の面積Ｓ及び正規化面積ＲＡの誤差（エラーバー）も合わせて示している。図６（ａ）及び図６（ｂ）に示すように、自律神経機能の悪化度合いに対する正規化面積ＲＡの誤差は、自律神経機能の悪化度合いに対する面積Ｓの誤差よりも小さくなる。誤差が小さくなる理由は、上述した減少量ΔＦでの正規化である。その結果、減少量ΔＦでの正規化により、面積Ｓに基づく自律神経機能の悪化度合いの評価精度と比較して、正規化面積ＲＡに基づく自律神経機能の悪化度合いの評価精度が高くなる。このように正規化面積ＲＡの誤差が相対的に小さくなるという点においても、評価装置１０は、正規化面積ＲＡを評価パラメータとして用いることで、自律神経機能を適切に評価することができる。

但し、図６（ａ）に示すように、自律神経機能の悪化度合いと面積Ｓとの間には、一定の相関が存在することもまた事実である。従って、評価装置１０は、正規化面積ＲＡに加えて又は代えて、面積Ｓに基づいて自律神経機能を評価してもよい。

加えて、面積Ｓを減少量ΔＦで正規化することで正規化面積ＲＡが得られるがゆえに、正規化面積ＲＡの次元は、時間の次元と一致する。つまり、評価装置１０は、被験者９００の血流量の時間波形がどのような反応型であっても、血流量の反応型を決定付ける立位開始時刻Ｔｓから回復時刻Ｔｒまでの期間の長さと同じ次元を有する「正規化面積ＲＡという評価パラメータ」に基づいて、自律神経機能を評価することができる。この意味においても、評価装置１０は、自律神経機能をより高精度に評価することができる。

加えて、評価装置１０は、カットオフ時刻Ｔｃを採用することで、血流量の反応型が不完全回復型となる場合であっても、正規化面積ＲＡを算出することができる。つまり、評価装置１０は、カットオフ時刻Ｔｃを採用することで、被験者９００の血流量の時間波形がどのような反応型であっても、正規化面積ＲＡという共通の評価パラメータに基づいて自律神経機能を評価することができる。

（４）変形例
以下、本実施例の評価システム１の変形例について説明する。

（４−１）第１変形例
初めに、図７（ａ）から図７（ｃ）を参照しながら、第１変形例について説明する。図７（ａ）から図７（ｃ）は、夫々、血流情報が示す血流量の時間波形上で第１変形例の面積Ｓを示すグラフである。

上述した実施例では、評価パラメータ演算部１４１は、血流量の時間波形のうち立位開始時刻Ｔｓから回復時刻Ｔｒ又はカットオフ時刻Ｔｃまでの間の期間中の血流量の変化（推移）を示す波形部分と平均値Ｆａｖｅとが囲む所定領域の面積Ｓを算出する。一方で、第１変形例では、図７（ａ）から図７（ｃ）に示すように、評価パラメータ演算部１４１は、血流量の時間波形のうち最小時刻Ｔｍｉｎから回復時刻Ｔｒ又はカットオフ時刻Ｔｃまでの間の期間中の血流量の変化（推移）を示す波形部分と平均値Ｆａｖｅとが囲む所定領域の面積Ｓを算出する。

このような面積Ｓを用いる第１変形例においても、評価装置１０は、上述した各種効果を享受することができる。

加えて、被験者９００が立位姿勢に移行することで減少した血流量が被験者９００の血圧調整機能に起因して増加することは上述したとおりである。血圧調整機能は、自律神経機能に密接に関連することもまた上述したとおりである。従って、血流量の時間波形のうち被験者９００が立位姿勢に移行することで減少した血流量が増加している期間中の波形部分は、血圧調整機能の悪化度合い（つまり、自律神経機能の悪化度合い）により密接に関連している可能性が相対的に高い。従って、第１変形例では、評価装置１０は、自律神経機能を評価するために、血流量の時間波形のうち自律神経機能の悪化度合いにより密接に関連する波形部分に着目することができる。このため、評価装置１０は、自律神経機能をより高精度に評価することができる。

尚、評価パラメータ演算部１４１は、血流量の時間波形のうち立位開始時刻Ｔｓと最小時刻Ｔｍｉｎとの間の任意の時刻から回復時刻Ｔｒ又はカットオフ時刻Ｔｃまでの間の期間中の血流量の変化（推移）を示す波形部分と平均値Ｆａｖｅとが囲む所定領域の面積Ｓを算出してもよい。評価パラメータ演算部１４１は、血流量の時間波形のうち最小時刻Ｔｍｉｎと回復時刻Ｔｒ又はカットオフ時刻Ｔｃとの間の任意の時刻から回復時刻Ｔｒ又はカットオフ時刻Ｔｃまでの間の期間中の血流量の変化（推移）を示す波形部分と平均値Ｆａｖｅとが囲む所定領域の面積Ｓを算出してもよい。つまり、面積Ｓを規定する時間波形の開始部分を特定する時刻が任意に調整されてもよい。

評価パラメータ演算部１４１は、血流量の時間波形のうち立位開始時刻Ｔｓから立位開始時刻Ｔｓと回復時刻Ｔｒ又はカットオフ時刻Ｔｃとの間の任意の時刻までの間の期間中の血流量の変化（推移）を示す波形部分と平均値Ｆａｖｅとが囲む所定領域の面積Ｓを算出してもよい。つまり、面積Ｓを規定する時間波形の終了部分を特定する時刻が任意に調整されてもよい。

評価パラメータ演算部１４１は、血流量の時間波形のうち立位開始時刻Ｔｓと回復時刻Ｔｒ又はカットオフ時刻Ｔｃとの間の任意の第１時刻から立位開始時刻Ｔｓと回復時刻Ｔｒ又はカットオフ時刻Ｔｃとの間の任意の第２時刻（但し、第２時刻は第１時刻よりも遅い）までの間の期間中の血流量の変化（推移）を示す波形部分と平均値Ｆａｖｅとが囲む所定領域の面積Ｓを算出してもよい。つまり、面積Ｓを規定する時間波形の開始部分及び終了部分の双方を特定する時刻が任意に調整されてもよい。

（４−２）第２変形例
続いて、図８（ａ）から図８（ｃ）を参照しながら、第２変形例について説明する。図８（ａ）から図８（ｃ）は、夫々、血流情報が示す血流量の時間波形上で第２変形例の面積Ｓを示すグラフである。

上述した実施例では、回復時刻Ｔｒは、被験者９００が立位姿勢に移行したことに起因して減少した血流量が平均値Ｆａｖｅにまで回復した時刻である。一方で、第２変形例では、図８（ａ）から図８（ｃ）に示すように、回復時刻Ｔｒは、被験者９００が立位姿勢に移行したことに起因して減少した血流量が平均値Ｆａｖｅよりも小さい中間値Ｆｃにまで回復した時刻である。従って、第２変形例では、評価パラメータ演算部１４１は、血流量の時間波形のうち立位開始時刻Ｔｓから回復時刻Ｔｒ又はカットオフ時刻Ｔｃまでの間の期間中の血流量の変化（推移）を示す波形部分と中間値Ｆｃとが囲む所定領域の面積Ｓを算出する。

中間値Ｆｃは、平均値Ｆａｖｅに対して所定の係数ＲＦ２（但し、０≦ＲＦ２≦１）を掛け合わせることで得られる値であってもよい。係数ＲＦ２は、０．５以上であってもよいし、０．８以上であってもよい。係数ＲＦ２は、０．９以上であることが好ましい。係数ＲＦ２は、０．９５以上であることがより好ましい。

中間値Ｆｃは、被験者９００が立位姿勢に移行することで減少量ΔＦだけ減少した血流量が、減少量ΔＦ×所定の係数ＲＦ１（但し、０≦ＲＦ１≦１）だけ回復したと仮定した場合の血流量であってもよい。つまり、中間値Ｆｃは、平均値Ｆａｖｅ−減少量ΔＦ＋減少量ΔＦ×係数ＲＦ１＝平均値Ｆａｖｅ−減少量ΔＦ（１−係数ＲＦ１）という数式に基づいて定められてもよい。係数ＲＦ１は、０．５以上であってもよいし、０．８以上であってもよい。係数ＲＦ１は、０．９以上であることが好ましい。係数ＲＦ１は、０．９５以上であることがより好ましい。

このような面積Ｓを用いる第２変形例においても、評価装置１０は、上述した各種効果を享受することができる。

（４−３）第３変形例
続いて、図９（ａ）から図９（ｃ）を参照しながら、第３変形例について説明する。図９（ａ）から図９（ｃ）は、夫々、血流情報が示す血流量の時間波形上で第３変形例の面積Ｓを示すグラフである。

上述した実施例では、評価パラメータ演算部１４１は、血流量の時間波形のうち立位開始時刻Ｔｓから回復時刻Ｔｒ又はカットオフ時刻Ｔｃまでの間の期間中の血流量の変化（推移）を示す波形部分と平均値Ｆａｖｅとが囲む所定領域の面積Ｓを算出する。一方で、第３変形例では、図９（ａ）から図９（ｃ）に示すように、評価パラメータ演算部１４１は、血流量の時間波形のうち立位開始時刻Ｔｓから最小時刻Ｔｍｉｎまでの間の期間中の血流量の変化（推移）を示す波形部分と平均値Ｆａｖｅとが囲む所定領域の面積Ｓ１、及び、最小時刻Ｔｍｉｎから回復時刻Ｔｒ又はカットオフ時刻Ｔｃまでの間の期間中の血流量の変化（推移）を示す波形部分と平均値Ｆａｖｅとが囲む所定領域の面積Ｓ２を別個に算出する。その後、パラメータ演算部１４１は、減少量ΔＦを用いて面積Ｓ１を正規化することで正規化面積ＲＡ１を算出し、減少量ΔＦを用いて面積Ｓ２を正規化することで正規化面積ＲＡ２を算出する。

その後、評価結果出力部１４２は、正規化面積ＲＡ１及びＲＡ２に基づいて、自律神経機能を評価する。例えば、評価結果出力部１４２は、正規化面積ＲＡ１に基づいて自律神経機能を評価し、正規化面積ＲＡ２に基づいて自律神経機能を評価し、正規化面積ＲＡ１に基づく評価結果及び正規化面積ＲＡ２に基づく評価結果に基づいて、自律神経機能を評価してもよい。或いは、評価結果出力部１４２は、正規化面積ＲＡ１及びＲＡ２から算出されるパラメータ（例えば、正規化面積ＲＡ１と正規化面積ＲＡ２との和又は比率等）に基づいて、自律神経機能を評価してもよい。

このような面積を用いる第３変形例においても、評価装置１０は、上述した各種効果を享受することができる。

加えて、血流量の時間波形のうち被験者９００が立位姿勢に移行することで減少した血流量が増加している期間中の波形部分が血圧調整機能の悪化度合い（つまり、自律神経機能の悪化度合い）により密接に関連している可能性が相対的に高いことは、第１変形例において上述したとおりである。つまり、正規化面積ＲＡ２は、自律神経機能の悪化度合いに関連している可能性が相対的に高い。一方で、血流量の時間波形のうち被験者９００が立位姿勢に移行することで血流量が減少している期間中の波形部分もまた、血圧調整機能の悪化度合い（つまり、自律神経機能の悪化度合い）にある程度関連している可能性はある。例えば、血流量の時間波形のうち被験者９００が立位姿勢に移行することで血流量が減少している期間中の波形部分は、血圧調整機能が作用し始めるまでに要する時間を示しているとも言える。従って、正規化面積ＲＡ１は、自律神経機能のうち血圧調整機能の立ち上げ（アクティベーション）と言う第１機能部分に関連しており、正規化面積ＲＡ２は、自律神経機能のうち血圧調整機能の立ち上げ後の動作という第２機能部分に関連していると言える。つまり、正規化面積ＲＡ１及びＲＡ２は、夫々、自律神経機能の異なる機能部分の悪化度合いに関連している可能性が相対的に高いと言える。第３変形例では、パラメータ演算部１４１は、自律神経機能のうち異なる機能部分に夫々関連する複数の面積Ｓ（複数の正規化面積ＲＡ）を算出することができる。従って、評価装置１０は、自律神経機能を高精度に評価することができる。

尚、評価パラメータ演算部１４１は、３つ以上の面積Ｓ（言い換えれば、正規化面積ＲＡ）を算出してもよい。つまり、評価パラメータ演算部１４１は、血流量の時間波形のうち立位開始時刻Ｔｓと回復時刻Ｔｒ又はカットオフ時刻Ｔｃとの間における任意の期間＃１中の血流量の変化（推移）を示す波形部分＃１と平均値Ｆａｖｅとが囲む所定領域の面積Ｓ１と、血流量の時間波形のうち立位開始時刻Ｔｓと回復時刻Ｔｒ又はカットオフ時刻Ｔｃとの間における任意の期間＃２中の血流量の変化（推移）を示す波形部分＃２と平均値Ｆａｖｅとが囲む所定領域の面積Ｓ２と、・・・、血流量の時間波形のうち立位開始時刻Ｔｓと回復時刻Ｔｒ又はカットオフ時刻Ｔｃとの間における任意の期間＃ｋ（但し、ｋは、３以上の整数）中の血流量の変化（推移）を示す波形部分＃ｋと平均値Ｆａｖｅとが囲む所定領域の面積Ｓｋとを算出してもよい。

また、第３変形例においても、第２変形例と同様に、回復時刻Ｔｒは、被験者９００が立位姿勢に移行したことに起因して減少した血流量が平均値Ｆａｖｅよりも小さい中間値Ｆｃにまで回復した時刻であってもよい。つまり、評価パラメータ演算部１４１は、第２変形例と同様の算出態様で面積Ｓ１及びＳ２の夫々を算出してもよい。具体的には、図１０（ａ）から図１０（ｃ）に示すように、評価パラメータ演算部１４１は、血流量の時間波形のうち立位開始時刻Ｔｓから最小時刻Ｔｍｉｎまでの間の期間中の血流量の変化（推移）を示す波形部分と中間値Ｆｃとが囲む所定領域の面積Ｓ１、及び、最小時刻Ｔｍｉｎから回復時刻Ｔｒ又はカットオフ時刻Ｔｃまでの間の期間中の血流量の変化（推移）を示す波形部分と中間値Ｆｃとが囲む所定領域の面積Ｓ２を別個に算出してもよい。

また、第３変形例において、血流量が中間値Ｆｃにまで回復した時刻が回復時刻Ｔｒとなる場合には、図１１（ａ）から図１１（ｃ）に示すように、評価パラメータ演算部１４１は、血流量の時間波形のうち立位開始時刻Ｔｓから最小時刻Ｔｍｉｎまでの間の期間中の血流量の変化（推移）を示す波形部分と平均値Ｆａｖｅとが囲む所定領域の面積Ｓ１、及び、最小時刻Ｔｍｉｎから回復時刻Ｔｒ又はカットオフ時刻Ｔｃまでの間の期間中の血流量の変化（推移）を示す波形部分と中間値Ｆｃとが囲む所定領域の面積Ｓ２を別個に算出してもよい。自律神経機能を評価してもよい。

（４−４）第４変形例
続いて、第４変形例について説明する。上述した実施例では、被験者９００の血流量を変化させる（具体的には、一時的に血流量を減少させる）ための状態変化として、被験者９００の姿勢が座位姿勢から立位姿勢に移行するという状態変化が用いられている。しかしながら、被験者９００の血流量を変化させるための状態変化として、被験者９００の状態が任意の第１状態から任意の第２状態に移行するという任意の状態変化が用いられてもよい。以下、被験者９００の血流量を変化させるための状態変化の他の一例について説明する。
被験者９００の血流量を変化させるための状態変化は、被験者９００がベッド上での臥位姿勢から床へ立位姿勢をとるという状態変化であってもよい。この場合、状態変化が生ずる前の被験者９００の第１状態は、被験者９００が臥位姿勢をとっている状態である。一方で、状態変化が生じた後の被験者９００の第２状態は、被験者９００が立位姿勢であるという状態である。

被験者９００の血流量を変化させるための状態変化は、被験者９００の下肢大腿部（或いは、その他任意の部位）を圧迫しているカフが取り外される又は緩められるという状態変化であってもよい。この場合、状態変化が生ずる前の被験者９００の第１状態は、被験者９００の下肢大腿部がカフで圧迫されているという状態である。一方で、状態変化が生じた後の被験者９００の第２状態は、被験者９００の下肢大腿部がカフで圧迫されていないという状態である。

被験者９００の血流量を変化させるための状態変化は、被験者９００の下半身（その他任意の部位）を負圧環境下にさらすという状態変化であってもよい。この場合、状態変化が生ずる前の被験者９００の第１状態は、被験者９００の下半身が負圧環境下にさらされていないという状態である。一方で、状態変化が生じた後の被験者９００の第２状態は、被験者９００の下半身が負圧環境下にさらされているという状態である
被験者９００の血流量を変化させるための状態変化は、被験者９００の交感神経を刺激するという状態変化であってもよい。この場合、状態変化が生ずる前の被験者９００の第１状態は、被験者９００の交感神経が未だ刺激されていないという状態である。一方で、状態変化が生じた後の被験者９００の第２状態は、被験者９００の交感神経が既に刺激されているという状態である。尚、ここで言う刺激とは、被験者９００の体に直接的に触れる物理的な刺激であってもよい。或いは、ここで言う刺激とは、被験者９００の精神に影響を与える精神的な刺激であってもよい。

また、本発明は、請求の範囲及び明細書全体から読み取るこのできる発明の要旨又は思想に反しない範囲で適宜変更可能であり、そのような変更を伴う評価装置、出力方法、及び、コンピュータプログラムもまた本発明の技術思想に含まれる。

１０ 評価装置
１１ 制御部
１２ 取得部
１３ 指示部
１４ 評価部
１４１ パラメータ演算部
１４２ 評価結果出力部
２０ 血流計
３０ ディスプレイ
９００ 被験者
Ｔａ 座位開始時刻
Ｔｓ 立位開始時刻
Ｔｍｉｎ 最小時刻
Ｔｒ 回復時刻
Ｔｃ カットオフ時刻
Ｆａｖｅ 平均値
Ｆｍｉｎ 最小値
Ｆｃ 中間値
ΔＦ 減少量
Ｓ 面積
ＲＡ 正規化面積

Claims (15)

1. 被験者の血流量を取得する取得部と、
   前記血流量を変化させるための状態変化が前記被験者に生じた後の前記血流量の変化を示す時間波形が規定する所定領域の面積に基づいて、前記被験者の生体機能を評価する評価部と
   を備える評価装置。
2. 前記評価部は、前記状態変化に起因した前記血流量の減少量を用いて前記面積を正規化することで得られる正規化結果に基づいて、前記生体機能を評価する
   請求項１に記載の評価装置。
3. 前記所定領域は、前記時間波形のうち前記状態変化が生じた第１時刻から前記第１時刻よりも後の第２時刻までの間の第１期間中の前記血流量を示す少なくとも一部の波形部分が規定する領域を含む
   請求項１又は２に記載の評価装置。
4. 前記所定領域は、前記時間波形のうち前記状態変化が生じた第１時刻よりも後の第３時刻から前記第３時刻よりも後の第２時刻までの間の第２期間中の前記血流量を示す少なくとも一部の波形部分が規定する領域を含む
   請求項１から３のいずれか一項に記載の評価装置。
5. 前記所定領域は、前記時間波形のうち前記状態変化が生じた第１時刻から前記第１時刻よりも後の第３時刻までの間の第３期間中の前記血流量を示す少なくとも一部の波形部分が規定する領域、及び、前記時間波形のうち前記第３時刻から前記第３時刻よりも後の第２時刻までの間の第２期間中の前記血流量を示す少なくとも一部の波形部分が規定する領域を含む
   請求項１から４のいずれか一項に記載の評価装置。
6. 前記第３時刻は、前記状態変化に起因して前記血流量が最も減少した時刻である
   請求項４又は５に記載の評価装置。
7. 前記第２時刻は、前記状態変化に起因して減少した前記血流量が所定量にまで回復した回復時刻である
   請求項３から６のいずれか一項に記載の評価装置。
8. 前記状態変化に起因して減少した前記血流量が所定量にまで回復した回復時刻が、所定のカットオフ時刻よりも前である場合には、前記第２時刻は前記回復時刻であり、
   前記回復時刻が前記カットオフ時刻よりも後である場合には、前記第２時刻は前記カットオフ時刻である
   請求項３から７のいずれか一項に記載の評価装置。
9. 前記カットオフ時刻は、前記生体機能が正常でない前記被験者に前記状態変化が生じた場合における前記回復時刻に基づいて規定される
   請求項８に記載の評価装置。
10. 前記所定量は、前記状態変化が生ずる前の第４期間中の前記血流量の平均値である
    請求項７から９のいずれか一項に記載の評価装置。
11. 前記少なくとも一部の波形が規定する前記領域は、前記平均値と前記少なくとも一部の波形部分とが囲む領域を含む
    請求項１０に記載の評価装置。
12. 前記所定量は、前記状態変化が生ずる前の第４期間中の前記血流量の平均値よりも小さい中間値である
    請求項７から９のいずれか一項に記載の評価装置。
13. 前記少なくとも一部の波形が規定する前記領域は、前記中間値と前記少なくとも一部の波形部分とが囲む領域を含む
    請求項１２に記載の評価装置。
14. 被験者の血流量を取得する取得工程と、
    前記血流量を変化させるための状態変化が前記被験者に生じた後の前記血流量の変化を示す時間波形が規定する所定領域の面積に基づく前記被験者の生体機能の評価結果を出力する出力工程と
    を備える出力方法。
15. コンピュータを請求項１から１３のいずれか一項に記載の評価装置として機能させることを特徴とするコンピュータプログラム。