"生体の生体機能の異なる複数の部位にそれぞれ光を照射し、出射する光を受光する複数の生体用プローブと、前記複数の生体用プローブによって検出された光情報を入力し、演算、制御又は記憶を行う装置本体とを有し、近赤外分光法を利用して生体機能を診断する生体機能診断装置であって、前記装置本体は、前記複数の各生体用プローブからの光情報に基づいて、酸化型ヘモグロビンの濃度変化量と脱酸化型ヘモグロビンの濃度変化量又はこれらの関係から導きだされる各種パラメータを含む生理的指標を算出する算出部と、前記算出部によって算出された前記生体の異なる部位に対応する生理的指標同士を比較可能に調整する調整部と、前記調整部によって調整された情報に基づいて画像表示する表示部とを有する、ことを特徴とする生体機能診断装置。前記調整部は、生体機能の異なる複数の部位における同一の生理的指標同士を比較可能に調整することを特徴とする請求項1に記載の生体機能診断装置。前記調整部は、生体機能の異なる複数の部位における異なる生理的指標同士を比較可能に調整することを特徴とする請求項1に記載の生体機能診断装置。前記算出部は、生体機能の異なる複数の部位における前記生理的指標同士の相関係数を算出し、前記表示部は、前記相関係数に関する情報を表示する、ことを特徴とする請求項1乃至3のいずれか1つの項に記載の生体機能診断装置。前記算出部は、生体機能の異なる複数の部位における同一の生理的指標同士を乗算し、前記表示部は、前記乗算値に関する情報を表示する、ことを特徴とする請求項1乃至4のいずれか1つの項に記載の生体機能診断装置。前記算出部は、生体機能の異なる複数の部位における異なる生理的指標同士を乗算し、前記表示部は、前記乗算値に関する情報を表示する、ことを特徴とする請求項1乃至4のいずれか1つの項に記載の生体機能診断装置。前記生体の生体機能の異なる複数の部位は、少なくとも脳の部位と筋肉の部位とを含むことを特徴とする請求項1乃至6のいずれか1つの項に記載の生体機能診断装置。前記生体の生体機能の異なる複数の部位は、少なくとも脳の異なる2つの部位を含むことを特徴とする請求項1乃至6のいずれか1つの項に記載の生体機能診断装置。前記算出部は、関心領域ROIにおける血液量BVの濃度変化量ΔBVを式1で算出する、ことを特徴とする請求項1乃至8のいずれか1つの項に記載の生体機能診断装置。ΔBVΔHbΔHbO2式1ここで、ΔHbは、脱酸化型ヘモグロビンの濃度変化量、ΔHbO2は酸化型ヘモグロビンの濃度変化量である。前記算出部は、脳と筋肉における血液量BVの濃度変化量ΔBV同士の関係を示す二次元ベクトルのダイアグラムから、任意の関心領域ROI中の筋肉血液量最大時間、脳血液量最大時間を算出する、ことを特徴とする請求項9に記載の生体機能診断装置。前記算出部は、関心領域ROIにおける酸素交換量OEの濃度変化量ΔOEを式2で算出する、ことを特徴とする請求項1乃至10のいずれか1つの項に記載の生体機能診断装置。ΔOEΔHb―ΔHbO2式2ここで、ΔHbは、脱酸化型ヘモグロビンの濃度変化量、ΔHbO2は酸化型ヘモグロビンの濃度変化量である。前記算出部は、脳と筋肉における酸素交換量OEの濃度変化量ΔOE同士の関係を示す二次元ベクトルのダイアグラムから、任意の関心領域ROI中の筋肉酸素消費最大時間、脳酸素消費最大時間を算出する、ことを特徴とする請求項11に記載の生体機能診断装置。前記算出部は、関心領域ROIにおける酸素交換量OEの濃度変化量ΔOEと血液量BVの濃度変化量ΔBVの比である比eを式3で算出する、ことを特徴とする請求項1乃至12のいずれか1つの項に記載の生体機能診断装置。e=ΔOEΔBVΔHb―ΔHbO2ΔHbΔHbO2式3ここで、ΔHbは、脱酸化型ヘモグロビンの濃度変化量、ΔHbO2は酸化型ヘモグロビンの濃度変化量である。前記算出部は、E角を式3aで算出することを特徴とする請求項13に記載の生体機能診断装置。E角=Arctan式3a前記算出部は、生体の異なる2つの部位のうち一方の部位の第1のE角を縦軸、他方の部位の第2のE角を横軸として経時的にプロットして得られる二次元ダイアグラムにおいて、第1のE角と第2のE角のなす角度であるEハイブリッド角EH角を、式3bで算出する、ことを特徴とする請求項14に記載の生体機能診断装置。EH角=Arctan式3b前記算出部は、2つの異なる部位における同一の生理的指標同士を乗算し、一方の部位の第1の乗算値を縦軸、他方の部位の第2の乗算値を横軸として経時的にプロットして得られる二次元ダイアグラムにおいて、第1の乗算値と第2の乗算値のなす角度である第1のハイブリッド角H1角を、式3cで算出する、ことを特徴とする請求項1乃至15のいずれか1つの項に記載の生体機能診断装置。H1角=Arctan式3c前記算出部は、2つの異なる部位における異なる生理的指標同士を乗算し、一方の部位の第3の乗算値を縦軸、他方の部位の第4の乗算値を横軸として経時的にプロットして得られる二次元ダイアグラムにおいて、第3の乗算値と第4の乗算値のなす角度である第2のハイブリッド角H2角を、式3dで算出する、ことを特徴とする請求項1乃至16のいずれか1つの項に記載の生体機能診断装置。H2角=Arctan式3d前記算出部は、脳と筋肉の複数の関心領域ROIsにおける運動時の脳筋肉酸素負担率BM比1とスカラーPL値1パワー値を二次元ダイアグラム上からと式4で算出する、ことを特徴とする請求項1乃至17のいずれか1つの項に記載の生体機能診断装置。BM比1=[筋肉の酸素交換量MOE変化量]/[脳の酸素交換量COE変化量]式4前記算出部は、脳と筋肉の複数の関心領域ROIsにおける運動時の脳筋肉血液量負担率BM比2とスカラーPL値2パワー値を二次元ダイアグラム上からと式5で算出する、ことを特徴とする請求項1乃至18のいずれか1つの項に記載の生体機能診断装置。BM比2=[筋肉の血液量MBV変化量]/[脳の血液量CBV変化量]式5前記算出部は、横軸をROI中の酸化型Hb量O、縦軸をROI中の脱酸化型Hb量Dとする二次元ダイアグラムから、酸素飽和度YをOD平面上の傾きY角として、式6で算出する、ことを特徴とする請求項1乃至19のいずれか1つの項に記載の生体機能診断装置。酸素飽和度Y=1-ArctanY角式6前記算出部は、横軸をROI中の酸化型Hb量O、縦軸をROI中の脱酸化型Hb量Dとする二次元ダイアグラムから、計測開始点E0の総ヘモグロビン量DOに対する酸素交換量D-Oの割合として定義されるE比血液量対酸素交換量の比を式7で算出する、ことを特徴とする請求項1乃至20のいずれか1つの項に記載の生体機能診断装置。E=D-O/DO式7前記算出部は、横軸をROI中の酸化型Hb量O、縦軸をROI中の脱酸化型Hb量Dとし、酸素飽和度Yの変化と酸化型Hbの変化量ΔO、脱酸化型Hbの変化量ΔDの関係を示すグラフから、酸素飽和度Yの変化ΔYを式8で算出する、ことを特徴とする請求項1乃至21のいずれか1つの項に記載の生体機能診断装置。酸素飽和度の変化ΔY=ArctanΔY角式8前記算出部は、関心領域ROIにおける酸素飽和度の推定変化量ΔYを実測値である血液量の最大変化量MaxΔBVに対応する任意の推定値の最大ヘマトクリットの変化MaxΔHtから式9でROI中の計測開始点の血液量BV0を算出する、ことを特徴とする請求項1乃至22のいずれか1つの項に記載の生体機能診断装置。MaxΔHt=MaxMaxΔBVBV0式9前記MaxΔBVは、複数のROI中で最大の変化を選択することを特徴とする請求項23に記載の生体機能診断装置。任意のROIに対して推定ヘマトクリット値を設定することを特徴とする請求項23に記載の生体機能診断装置。前記算出部は、関心領域ROIにおける酸素飽和度の推定変化量ΔYをROI中の計測開始点の血液量BV0を算出し、式10または式11に代入して求め、前記調整部は、異なる部位におけるΔYの時系列の変化を示すグラフを前記表示部で表示できるように調整する、ことを特徴とする請求項1乃至25のいずれか1つの項に記載の生体機能診断装置。ΔY=ΔOE2=ΔHb2ΔHbO22式12前記調整部は、異なる部位同士の前記生理的指標同士の積の変化を示すグラフを前記表示部で表示できるように調整する、ことを特徴とする請求項1乃至28のいずれか1つの項に記載の生体機能診断装置。前記表示部は、縦軸を前記生理的指標、横軸を時間として時系列変化を示すグラフを表示することを特徴とする請求項1乃至29のいずれか1つの項に記載の生体機能診断装置。前記表示部は、異なる2つの前記生理的指標の一方を縦軸、他方を横軸として、経時的にプロットして得られる二次元ダイアグラムを表示することを特徴とする請求項1乃至30のいずれか1つの項に記載の生体機能診断装置。前記表示部は、前記生理的指標の値の大きさに応じて色分けして表示された画像を表示することを特徴とする請求項1乃至31のいずれか1つの項に記載の生体機能診断装置。前記表示部は、安静時、前記生体への課題の負荷時、回復時の前記生理的指標の変化を同一のグラフ上で表示することを特徴とする請求項1乃至32のいずれか1つの項に記載の生体機能診断装置。前記表示部は、安静時、前記生体への課題の負荷時、回復時の前記生理的指標の変化をベクトル空間上のヘモグロビン関連のパラメータの波動、及び回転運動としての動態として、その方向と力の変化をそれぞれ異なるグラフ、または画像表示で表示することを特徴とする請求項1乃至33のいずれか1つの項に記載の生体機能診断装置。前記請求項1乃至34のいずれか1つの項に記載の生体機能診断装置の装置本体の処理を実行させることを特徴とするプログラム。従来例1においては次のような課題があった。1脳磁図やPET、fMRIは、仰向けで安静の状態で計測する必要があり、運動トレーニングのための脳計測には適さない。2脳波は、大脳の局在診断が大雑把で、脳の第一次運動野や補足運動野の働きを区別するためには、何百回もの計測が必要であり、1回の運動施行で脳波の変化を論じる感度と精度がない。3NIRSは脳と筋肉のトレーニング後の回復状態もリアルタイムに計測することができるが、NIRSが光の計測であるので、それぞれの部位を計測できても、光路長の違いや計測領域の体積評価が困難であった。そのため、脳と筋肉のNIRSを計測しても、酸化型ヘモグロビンと脱酸化型ヘモグロビンの変化の増減を定性的に比較できても定量性に乏しかった。41つのモダリティー医療機器:modalityを使って、同時に、CMRO2、CBF、CBV、OEFを計測することができないので、CMRO2酸素消費率を算出するために必要なCBF、CBVの値は、別々の実験で計測して、かつ複数のパラメータに対して推測値を用いなければならなかった。5酸素代謝に関係するCMRO2、OEFの単位は毎分あたりの定量であり、ミリ秒オーダーでリアルタイムにモニターすることはできなかった。6NIRSによる脳機能計測では、従来、頭皮からの計測によって脳実質までの距離が一定でないために発光と受光のプローブ間距離を等間隔に配置してもROI関心領域:regionofinterestにおけるヘモグロビン量光量や課題を開始する時点での各部位の組織酸素濃度が異なっていたとしても、一定として算出されてきた。この点は、NIRSによる脳機能計測だけではなく、脳機能計測全般に言える問題点であった。なお、頭皮からの計測によって脳実質までの距離が一定でないために、発光と受光のプローブ間距離を等間隔に配置してもROIにおけるヘモグロビン量光量が異なる。ヒトの頭部計測では、大人の場合、頭皮から脳実質までの距離が、ほとんどの部位で頭皮上から脳実質まで1cm前後であるので図41参照、受光と発光プローブの距離は、プローブ配列のように最低でも1cm以上、2cm程度確保されていることが望ましい。複数部位を計測するプローブは、頭部に用いる場合には、半径は10-20mm以上確保されることで、脳実質からの信号を含むことができる。受光と発光のプローブ間距離を1cm~2cm以上確保したままで、高密度にプローブを複数配列することにより、複数の計測領域では、各ROIをより均一にすることができるが、ROIのサイズを定量することができなかった。従って、複数の部位からNIRS計測する上での不均一なROIの問題を軽減する測定技術が望まれていた。従来例2においては次のような課題があった。7筋電図は、筋収縮の強さ、弱さを計測するためには有効であるが、筋肉が動かない状態を区別して診断評価することができなかった。8筋電図では、筋肉を動かしていない時の脳の状態を同時に計測できないので、運動が終わって休息している最中の筋肉と脳の状態を同時に電気的に計測して診断判断できなかった。従来例3においては次のような課題があった。9筋肉の加圧トレーニングは、筋肉にゴムバンドを巻くことで、筋肉の低酸素を引き起こして、効果的に筋トレーニングする方法である。しかし本発明者は、加圧トレーニングが広範囲に脳の血圧を上げて、脳血管に過度の圧を負荷するので、盲目的な加圧は危険性があることが考えられる。なお、脳圧をモニターする方法は、特許文献1に開示されている。10従来、筋肉は脳からの指令によって動いていると理解されてきた。しかし、筋肉疲労、あるいは、脳疲労のいずれでも筋肉は動かなくなるはずである。ところが、筋肉と脳の関係を計測して、筋肉と脳の疲労、回復の度合いをみて、トレーニングを開始したり、止めたりすることは困難であった。11筋肉にどれだけ負荷を加えることで、脳の特定の番地部位をトレーニングすることができるかは脳を見ながら筋肉トレーニングする必要があったが、従来の技術では困難であった。12脳の疲労によって筋肉が動かない時に筋トレをしても効果的に作用していない可能性があったり、脳への負荷がない状態での筋肉トレーニングのレベルを設定することが困難であった。13筋肉と大脳皮質を結ぶ皮質脊髄路は、M1以外に、運動前野、一次感覚野からも投射される。筋肉トレーニングで、脳区分用番地のどこを鍛えているのか、筋肉を脳とのかかわりで、どの状態のトレーニング負荷なのか分からないために、トレーニングメニューとトレーニングの開始、停止、休息のタイミングも適切に行われていなかった。14脳と筋肉のNIRSの同時測定装置が無かった。同時に測定する時には、頭皮上のプローブと筋肉でのプローブの形状が異なり、計測されたデータと解析、診断法が見つかっていなかった。15脳と筋肉の機能の関係をリアルタイムに画像化して表示する方法も見つかっていなかった。本発明は、上記課題を解決するためになされたものであり、NIRSを用いて脳と筋肉、脳の異なる箇所、筋肉の異なる箇所等生体機能の異なる複数の部位の生体機能を同時に計測し、診断し、画像化して表示することが可能な生体機能診断装置及びプログラムを提供することを目的とする。本発明の生体機能診断装置は、生体の生体機能の異なる複数の部位にそれぞれ光を照射し、出射する光を受光する複数の生体用プローブと、前記複数の生体用プローブによって検出された光情報を入力し、演算、制御又は記憶を行う装置本体とを有し、近赤外分光法を利用して生体機能を診断する生体機能診断装置であって、前記装置本体は、前記複数の各生体用プローブからの光情報に基づいて、酸化型ヘモグロビンの濃度変化量と脱酸化型ヘモグロビンの濃度変化量又はこれらの関係から導きだされる各種パラメータを含む生理的指標を算出する算出部と、前記算出部によって算出された前記生体の異なる部位に対応する生理的指標同士を比較可能に調整する調整部と、前記調整部によって調整された情報に基づいて画像表示する表示部とを有することを特徴とするものである。前記調整部は、生体機能の異なる複数の部位における同一の生理的指標同士を比較可能に調整してもよい。前記調整部は、生体機能の異なる複数の部位における異なる生理的指標同士を比較可能に調整してもよい。前記算出部は、生体機能の異なる複数の部位における前記生理的指標同士の相関係数を算出し、前記表示部は、前記相関係数に関する情報を表示する、ものでもよい。前記算出部は、生体機能の異なる複数の部位における同一の生理的指標同士を乗算し、前記表示部は、前記乗算値に関する情報を表示する、ものでもよい前記算出部は、生体機能の異なる複数の部位における異なる生理的指標同士を乗算し、前記表示部は、前記乗算値に関する情報を表示する、ものでもよい。前記生体の生体機能の異なる複数の部位は、少なくとも脳の部位と筋肉の部位とを含んでもよい。前記生体の生体機能の異なる複数の部位は、少なくとも脳の異なる2つの部位を含んでもよい。前記算出部は、関心領域ROIにおける血液量BVの濃度変化量ΔBVを式1で算出する、ものでもよい。ΔBVΔHbΔHbO2式1ここで、ΔHbは、脱酸化型ヘモグロビンの濃度変化量、ΔHbO2は酸化型ヘモグロビンの濃度変化量である。前記算出部は、脳と筋肉における血液量BVの濃度変化量ΔBV同士の関係を示す二次元ベクトルのダイアグラムから、任意の関心領域ROI中の筋肉血液量最大時間、脳血液量最大時間を算出する、ものでもよい。前記算出部は、関心領域ROIにおける酸素交換量OEの濃度変化量ΔOEを式2で算出する、ものでもよい。ΔOEΔHb―ΔHbO2式2ここで、ΔHbは、脱酸化型ヘモグロビンの濃度変化量、ΔHbO2は酸化型ヘモグロビンの濃度変化量である。前記算出部は、脳と筋肉における酸素交換量OEの濃度変化量ΔOE同士の関係を示す二次元ベクトルのダイアグラムから、任意の関心領域ROI中の筋肉酸素消費最大時間、脳酸素消費最大時間を算出する、ものでもよい。前記算出部は、関心領域ROIにおける酸素交換量OEの濃度変化量ΔOEと血液量BVの濃度変化量ΔBVの比である比eを式3で算出する、ものでもよい。e=ΔOEΔBVΔHb―ΔHbO2ΔHbΔHbO2式3ここで、ΔHbは、脱酸化型ヘモグロビンの濃度変化量、ΔHbO2は酸化型ヘモグロビンの濃度変化量である。前記算出部は、E角を式3aで算出するものでもよい。E角=Arctan式3a前記算出部は、生体の異なる2つの部位のうち一方の部位の第1のE角を縦軸、他方の部位の第2のE角を横軸として経時的にプロットして得られる二次元ダイアグラムにおいて、第1のE角と第2のE角のなす角度であるEハイブリッド角EH角を、式3bで算出する、ものでもよい。EH角=Arctan式3b前記算出部は、2つの異なる部位における同一の生理的指標同士を乗算し、一方の部位の第1の乗算値を縦軸、他方の部位の第2の乗算値を横軸として経時的にプロットして得られる二次元ダイアグラムにおいて、第1の乗算値と第2の乗算値のなす角度である第1のハイブリッド角H1角を、式3cで算出する、ものでもよい。H1角=Arctan式3c前記算出部は、2つの異なる部位における異なる生理的指標同士を乗算し、一方の部位の第3の乗算値を縦軸、他方の部位の第4の乗算値を横軸として経時的にプロットして得られる二次元ダイアグラムにおいて、第3の乗算値と第4の乗算値のなす角度である第2のハイブリッド角H2角を、式3dで算出する、ものでもよい。H2角=Arctan式3d前記算出部は、脳と筋肉の複数の関心領域ROIsにおける運動時の脳筋肉酸素負担率BM比1とスカラーPL値1パワー値を二次元ダイアグラム上からと式4で算出する、ものでもよい。BM比1=[筋肉の酸素交換量MOE変化量]/[脳の酸素交換量COE変化量]式4前記算出部は、脳と筋肉の複数の関心領域ROIsにおける運動時の脳筋肉血液量負担率BM比2とスカラーPL値2パワー値を二次元ダイアグラム上からと式5で算出する、ものでもよい。BM比2=[筋肉の血液量MBV変化量]/[脳の血液量CBV変化量]式5前記算出部は、横軸をROI中の酸化型Hb量O、縦軸をROI中の脱酸化型Hb量Dとする二次元ダイアグラムから、酸素飽和度YをOD平面上の傾きY角として、式6で算出する、ものでもよい。酸素飽和度Y=1-ArctanY角式6前記算出部は、横軸をROI中の酸化型Hb量O、縦軸をROI中の脱酸化型Hb量Dとする二次元ダイアグラムから、計測開始点E0の総ヘモグロビン量DOに対する酸素交換量D-Oの割合として定義されるE比血液量対酸素交換量の比を式7で算出する、ものでもよい。E=D-O/DO式7前記算出部は、横軸をROI中の酸化型Hb量O、縦軸をROI中の脱酸化型Hb量Dとし、酸素飽和度Yの変化と酸化型Hbの変化量ΔO、脱酸化型Hbの変化量ΔDの関係を示すグラフから、酸素飽和度Yの変化ΔYを式8で算出する、ものでもよい。酸素飽和度の変化ΔY=ArctanΔY角式8前記算出部は、関心領域ROIにおける酸素飽和度の推定変化量ΔYを実測値である血液量の最大変化量MaxΔBVに対応する任意の推定値の最大ヘマトクリットの変化MaxΔHtから式9でROI中の計測開始点の血液量BV0を算出する、ものでもよい。MaxΔHt=MaxMaxΔBVBV0式9前記MaxΔBVは、複数のROI中で最大の変化を選択するものでもよい。任意のROIに対して推定ヘマトクリット値を設定してもよい。前記算出部は、関心領域ROIにおける酸素飽和度の推定変化量ΔYをROI中の計測開始点の血液量BV0を算出し、式10または式11に代入して求め、前記調整部は、異なる部位におけるΔYの時系列の変化を示すグラフを前記表示部で表示できるように調整する、ものでもよい。ΔY=ΔOE2=ΔHb2ΔHbO22式12前記調整部は、異なる部位同士の前記生理的指標同士の積の変化を示すグラフを前記表示部で表示できるように調整する、ものでもよい。前記表示部は、縦軸を前記生理的指標、横軸を時間として時系列変化を示すグラフを表示してもよい。前記表示部は、異なる2つの前記生理的指標の一方を縦軸、他方を横軸として、経時的にプロットして得られる二次元ダイアグラムを表示してもよい。前記表示部は、前記生理的指標の値の大きさに応じて色分けして表示された画像を表示してもよい。前記表示部は、安静時、前記生体への課題の負荷時、回復時の前記生理的指標の変化を同一のグラフ上で表示してもよい。前記表示部は、安静時、前記生体への課題の負荷時、回復時の前記生理的指標の変化をベクトル空間上のヘモグロビン関連のパラメータの波動、及び回転運動としての動態として、その方向と力の変化をそれぞれ異なるグラフ、または画像表示で表示してもよい。本発明のプログラムは、前記生体機能診断装置の装置本体の処理を実行させることを特徴とするものである。本発明によれば、次のような優れた効果を奏する。1脳と筋肉、脳の異なる箇所、筋肉の異なる箇所等生体機能の異なる複数の部位でNIRSにより同時計測する際に、ROIのサイズが不均一であっても、異なる部位同士の生体機能を比較して診断することができる。2該当する運動によって、もっとも強く酸素負荷あるいは血圧負荷が起こる脳の部位を経時的に、診断、表示して、運動トレーニングの継続、中止、ギアチェンジを行うことができる。3運動負荷に比べて、脳の部位への負荷が少ない状態でトレーニングができる。4筋肉への負荷をぎりぎりの状態まで強くした状態で脳へのトレーニングができる。5脳と筋肉の回復期を、以下の4分類に診断して、トレーニングの開始、休息継続、ギアチェンジを行うことができる。脳筋肉回復済み、回復済み未回復未回復回復済み未回復未回復回復済み6筋肉疲労のステージと脳疲労のステージが異なっているので、ダンベル等の負荷を変える事で、脳トレと筋トレの比重を調整できる。7筋肉と脳の疲労を見ながら、トレーニングの負荷量を決めることができることで、効率良く筋トレと脳トレを同時に行うことができる。8運動のインターバルトレーニングを効果的に行うことができる。9高血圧や脳血管に問題のあるひとのリハビリを有効に行うことができる。10手で介助したときに脳と筋肉における酸素消費と血液量の負担率BM比が極端に変化するので、それをみて、効率よくリハビリを行うことができる。11反磁性体である酸化型ヘモグロビンの変化を検出できないために、fMRIでは常反磁性体である脱酸化型ヘモグロビンの変化だけを用いて酸素飽和度の変化を推定して脳機能マッピングを行ってきた。しかし脱酸化型ヘモグロビンの変化だけを用いた方法では、正確な酸素飽和度の変化が得られていないことが明らかとなったので、酸化型ヘモグロビンの変化と脱酸化型ヘモグロビンの変化の両方を同時に用いたより精度の高い酸素飽和度の時系列変化を用いた脳機能計測ができる。12さらに、酸化型ヘモグロビンと脱酸化型ヘモグロビンをベクトル空間上で波や運動として変化するダイナミクスを表す指標として捉えることで、酸化型ヘモグロビンの変化と脱酸化型ヘモグロビンの変化の両方の変化の関係を表す複数の指標を創出することができる。13脳と筋肉の機能の関係をリアルタイムに画像化して表示することができる。以下、本発明の実施の形態について図面を参照して説明する。生体機能診断装置の概略図1は本発明の実施形態例に係る生体機能診断装置の構成を示すブロック図、図2は本発明の実施形態例に係る生体機能診断装置の使用例を概略的に示す説明図である。本発明の実施形態例に係る生体機能診断装置Kは、図1及び図2に示すように、生体の生体機能の異なる複数の部位にそれぞれ光を照射し、出射する光を受光する複数の生体用プローブ1と、複数の生体用プローブ1によって検出された光情報を入力し、演算、制御又は記憶を行う装置本体2とを有し、近赤外分光法を利用して生体機能を診断するために用いられる。各生体用プローブ1は、生体の任意の計測部位組織に光を照射する少なくとも2以上の発光素子発光ダイオード1a...と、計測部位からの透過光、反射光あるいは散乱光等、生体と相互作用した後の光を受光する少なくとも2以上の受光素子フォトダイオード1b...とで構成されている。なお、生体用プローブ1は、脳、腕の筋肉、顎の筋肉、指先の筋肉、口腔の歯肉等に設置して計測されるが、設置箇所に応じて生体用プローブ1の形状、計測面積、設置方法等が設定される。発光素子1aと受光素子1bとの間隔は、脳の計測の場合には1.5cmから3cm程度、腕の筋肉や顎の筋肉の計測の場合には2-3cm程度、指先の筋肉の計測の場合には5-10mm程度、口腔の歯肉の計測の場合には2-3.5mm程度である。また、生体用プローブ1を腕の筋肉に配置する際、腕に対して平行に配置する場合には、特定の筋肉の縦断面の伸びと収縮に伴う筋肉の酸素代謝を計測することができる。筋肉の縦断面の方向では、筋肉の縦方向の伸び縮みが反映されるからである。生体用プローブ1を腕に対して垂直に配置する場合には、特定の筋肉の横断面の伸びと収縮に伴う筋肉の酸素代謝を計測することができる。筋肉の横断面では、筋肉の屈曲伸展に伴って、横断面の面積変化が反映されるからである。また、生体用プローブ1は、その発光素子1aと受光素子1bとの距離を等間隔に複数マトリックス状に配列するのが好ましい。ただ、選択した筋肉運動と最も関連性のある脳部位を、複数の脳の計測部位から、もっとも酸素代謝の変化増加、減少を引き起こす場所と脳血液量の変化増加、減少を引き起こす場所、もっとも相関性の高い部位を同定する目的のため、発光素子1aと受光素子1bとの距離は、必ずしも等間隔である必要なく、ランダムの配置することも可能である。また、一端、関連性の強い部位が判明すれば、さらに、高密度の発光素子1aと受光素子1bの距離も短い生体用プローブ1の配列をして、該当部位を精密に同定する2段階計測法も可能である。すなわち、従来問題であった、計測部位ごとの光のS/Nの違い、光路長の違い、生体用プローブ1にはさまれたサンプリングする領域の大きさの違いの影響を排除して、筋肉と脳の関連性を調べることができる。装置本体2は、発光素子1a...の発光光量を調節する光量調節部3と、任意の受光素子1b...を選択的に有効化又は無効化し、全部の計測感度を調節する選択調節部4と、受光素子1b...からの信号を増幅するゲイン制御可能な信号増幅部5と、信号増幅部5の出力を数値化するA/D変換部6と、各部の制御処理やA/D変換部6の出力に基づいて所定の演算処理を実行する制御部7と、A/D変換部6の出力、各部の制御用データあるいは演算結果等の記憶に使用される記憶部8と、A/D変換部6の出力結果や演算結果等に基づく表示を行う表示部9とを有する。なお、装置本体2には、各種データを印刷したり、通信ネットワークを介してデータを送受信する機能を有してもよい。制御部7は、複数の各生体用プローブ1からの光情報に基づいて、酸化型ヘモグロビンの濃度変化量と脱酸化型ヘモグロビンの濃度変化量又はこれらの関係から導きだされるパラメータを算出する算出部10と、算出部10によって算出された生体の異なる部位に対応する変化量又はパラメータ同士を比較可能に調整して、表示部9に表示できるようにする調整部11とを有する。生体用プローブ1の発光素子1a...は、波長730nmの光を照射するものと、波長850nmの光を照射するものの二種類が用意されるなお、この光の波長の数値は一例であり、これに限定されるものではなく、3波長以上組み合わせてもよい。これらは、例えば列方向に交互に配置されるが、その他のパターンを検討するに当たっては、組織中での波長に依存する減衰を考慮して、受光光量をバランスよく計測できるような配置にすることが重要である。全ての発光素子1a...は、装置本体2の光量調節部3に接続されており、全体的にあるいはそれぞれ独立的に発光光量の調節が可能である。一方、全ての受光素子1b...は、装置本体2の選択調節部4を介して信号増幅部5に接続されており、それぞれの受光素子1bから出力される受光信号は、全部あるいは一部が選択調節部4で選択調節された状態で信号増幅部5に出力され、ここで増幅される。そして、増幅された受光信号は、A/D変換部6で数値化されて制御部7に出力される。制御部7は、A/D変換部6から入力されたデジタルデータをローパスフィルタにかけてノイズ除去処理を実行した後、この処理データ以下、受光光量というをタイムテーブル的に記憶部8に記憶させる。また、制御部7は、得られた受光光量に基づき、以下に説明する演算処理を実行する。まず、波長730nmの吸光度OD730を式13により、そして、波長850nmの吸光度OD850を式14により算出すると共に、該算出結果をタイムテーブル的に記憶部8に記憶させる。OD730=log10I0730I730式13OD850=log10I0850I850式14I0730:波長730nmの発光光量I730:波長730nmの受光光量I0850:波長850nmの発光光量I850:波長850nmの受光光量ここで、酸化型ヘモグロビンの濃度変化量と脱酸化型ヘモグロビンの濃度変化量と吸光度変化量との間には、式15,式16の関係があることが公知理論によりわかっている。ΔOD730=a1Δ[HbO2]a1’Δ[Hb]式15ΔOD850=a2Δ[HbO2]a2’Δ[Hb]式16ΔOD730:波長730nmの吸光度変化量ΔOD850:波長850nmの吸光度変化量Δ[HbO2]:酸化型ヘモグロビンの濃度変化量Δ[Hb]:脱酸化型ヘモグロビンの濃度変化量a1,a1’,a2,a2’:吸光度係数従って、この公知の連立方程式から、式17,式18が求められる。Δ[HbO2]=a{ΔOD730-a1’a2’ΔOD850}式17Δ[Hb]=aa2a2’{a1a2ΔOD850-ΔOD730}式18a=a2’a1a2’-a1’a2≒11あるいはその近傍値そこで、波長730nmの吸光度変化量ΔOD730及び波長850nmの吸光度変化量ΔOD850を求めた上で、酸化型ヘモグロビンの濃度変化量Δ[HbO2]を式17により、そして、脱酸化型ヘモグロビンの濃度変化量Δ[Hb]を式18により算出すると共に、該算出結果をタイムテーブル的に記憶部8に記憶させる。なお、総ヘモグロビンの濃度変化量Δ[totalHb]は、式19で表される。Δ[totalHb]=Δ[HbO2]Δ[Hb]式19ところで、組織に対する刺激によって誘発される毛細血管中の酸化型ヘモグロビン及び脱酸化型ヘモグロビンの各濃度変化量の変化態様は、その増減の組み合わせによって以下の9パターンを示す。1Δ[HbO2]増加Δ[Hb]増加2Δ[HbO2]増加Δ[Hb]減少3Δ[HbO2]増加Δ[Hb]ゼロ4Δ[HbO2]減少Δ[Hb]増加5Δ[HbO2]減少Δ[Hb]減少6Δ[HbO2]減少Δ[Hb]ゼロ7Δ[HbO2]ゼロΔ[Hb]増加8Δ[HbO2]ゼロΔ[Hb]減少9Δ[HbO2]ゼロΔ[Hb]ゼロ実際のところ、組織の代謝活動は、刺激の印加条件や安静状態の生理的状態の違いにより経時的に上記のパターンが変化している。毛細血管中の酸化型ヘモグロビンから組織中に酸素を取り込むための血流代謝活動として、毛細血管のΔ[Hb]とΔ[HbO2]が変動する。そこで、本発明においては、酸化型ヘモグロビンの濃度変化量ΔHbO2]と脱酸化型ヘモグロビンの濃度変化量Δ[Hb]に基づいて導き出される各種パラメータを制御部7の算出部10により算出する。また、脳と筋肉、脳の異なる箇所、筋肉の異なる箇所等生体機能の異なる複数の部位の生体機能を同時に計測し、診断し、画像化して表示するために、調整部11は、算出部10によって算出された生体の異なる部位に対応する濃度変化量又はパラメータ同士を比較可能に調整して、表示部9に表示できるようにする。例えば、調整部11はグラフの縦軸や横軸の目盛りの単位、目盛りの間隔、画像のサイズ、形状、色等を調整する。なお、例えば筋肉と脳では変化量やパラメータの大きさが異なるので、目盛りの単位や間隔等を変えると直線傾き等も変わる。しかし、符号や二次元ダイアグラム上の軌跡の向き右回り、左回り等は変わることはなく、二次元ダイアグラムの象限も変わらないので、複合指標で可視化したり、動態を観察することは可能である。図3は、本発明の実施形態例に係る生体機能診断装置の動作を説明するためのフローチャートである。図3に示すように、まず、生体の所定の部位、例えば左脳の第一次運動野M1の部位を生体用プローブ1を用いて計測するステップS1と同時に、生体の他の部位、例えば、右腕の筋肉二頭筋等を生体用プローブ1を用いて計測するステップS2。次いで、制御部7の算出部10によって、複数の各生体用プローブ1からの光情報に基づいて、酸化型ヘモグロビンの濃度変化量と脱酸化型ヘモグロビンの濃度変化量又はこれらの関係から導きだされるパラメータを算出するステップS3。次いで、制御部7の調整部11によって、算出部10により算出された生体の異なる部位に対応する濃度変化量又はパラメータ同士を比較可能に調整しステップS4、各種グラフや画像等を表示部9に表示するステップS5。時系列変化を示すグラフによる表示図4は9kgのダンベルを持ち上げて負荷をかける課題を与えた際の脳左脳の第一次運動野と筋肉右腕の二頭筋での脱酸化型ヘモグロビンと酸化型ヘモグロビンの時系列変化を示し、Aは横軸が時間、右縦軸が脳左脳の第一次運動野の脱酸化型ヘモグロビンの濃度変化量、左縦軸が筋肉右腕での脱酸化型ヘモグロビンの濃度変化量のグラフ、Bは横軸が時間、右縦軸が脳左脳の第一次運動野の酸化型ヘモグロビンの濃度変化量、左縦軸が筋肉右腕での酸化型ヘモグロビンの濃度変化量のグラフである。図5は所定期間ごとに課題を与えた際の脳左脳の第一次運動野と筋肉右腕の二頭筋でのΔHb-ΔHbO2とΔHbΔHbO2の時系列変化を示し、Aは横軸が時間、右縦軸が脳左脳の第一次運動野M1のΔHb-ΔHbO2、左縦軸が筋肉右腕の二頭筋でのΔHb-ΔHbO2のグラフ、Bは横軸が時間、右縦軸が脳左脳の第一次運動野のΔHbΔHbO2、左縦軸が筋肉右腕の二頭筋でのΔHbΔHbO2のグラフである。図4及び図5からわかるように、脳と筋肉での各濃度変化量等の時系列データを並べて表示することにより、課題負荷によって濃度変化量等が増減することがわかる。しかし、脳と筋肉の相互関係を診断することは困難である。そこで、本発明の実施形態例に係る生体機能診断装置Kでは、二次元ダイアグラムを用いて生体の異なる部位に対応する濃度変化量又はパラメータ同士を比較することにより、異なる部位の相互関係を診断できるようにした。二次元ダイアグラムによる表示図6は、横軸を脳脳の第1次運動野M1の血液量のCBV変化量ΔHbΔHbO2、縦軸を筋肉腕の二頭筋の血液量のMBV変化量ΔHbΔHbO2として、同時計測して経時的にプロットして得られる二次元ダイアグラムである。ここで、CBVはCerebral脳のBloodVolumeの略、MBVはMuscular筋肉のBloodVolumeの略である。ここで、ΔHbは脱酸化型ヘモグロビンの濃度変化量、ΔHbO2は酸化型ヘモグロビンの濃度変化量である。筋肉運動に対する刺激によって誘発される脳の毛細血管中の酸化型ヘモグロビン及び脱酸化型ヘモグロビンの各濃度変化量の変化態様は、その増減の組み合わせによって以下の9パターンを示す。筋肉の総Hb量脳の総Hb量Δ[HbO2]ΔHbΔ[HbO2]ΔHb1増加増加第1象限2増加ゼロ3増加減少第4象限4減少増加第2象限5減少ゼロ6減少減少第3象限7ゼロ増加8ゼロゼロ9ゼロ減少図6で、第1象限は脳と筋肉の血圧増大を示し、筋肉がパンピング、脳の血圧が上がっている象限である。第2象限は脳血圧増大と筋肉の虚血、収縮増大を示し、脳と筋肉への負荷が十分確保されている象限である。第3象限は脳と筋肉の虚血を示し、脳にはまだ余力があり、筋肉への負荷が十分確保されている象限である。第4象限は筋肉の血圧増大が脳血圧低下より大きいことを示し、脳に予備力があり、筋肉がパンピングしている象限である。図7A-Eは、横軸を脳左脳の第1次運動野M1の周囲の血液量のCBV変化量ΔHbΔHbO2、縦軸を筋肉右腕の二頭筋の血液量のMBV変化量ΔHbΔHbO2として、同時計測して経時的にプロットして得られる二次元ダイアグラムであり、Aはダンベルを持ち上げてない場合、Bは2.5kgのダンベルを持ち上げた場合、Cは4.5kgのダンベルを持ち上げた場合、Dは7.0kgのダンベルを持ち上げた場合、Eは9.5kgのダンベルを持ち上げた場合を示し、Fは説明図である。なお、縦軸と横軸の単位はmmollである。図8は、横軸を脳脳の第1次運動野M1の酸素交換量のCOE変化量ΔHb-ΔHbO2、縦軸を筋肉腕の二頭筋の酸素交換量のMOE変化量ΔHb-ΔHbO2として、同時計測して経時的にプロットして得られる二次元ダイアグラムである。ここで、COEはCerebral脳のOxygenExchangeの略、MOEはMuscular筋肉のOxygenExchangeの略である。筋肉運動に対する刺激によって誘発される脳の毛細血管中の酸化型ヘモグロビン及び脱酸化型ヘモグロビンの各濃度変化量の変化態様は、その増減の組み合わせによって以下の9パターンを示す。筋肉の酸素消費脳の酸素消費ΔHb-Δ[HbO2]ΔHb-Δ[HbO2]1増加増加第1象限2増加ゼロ3増加減少第4象限4減少増加第2象限5減少ゼロ6減少減少第3象限7ゼロ増加8ゼロゼロ9ゼロ減少図8で、第1象限は脳と筋肉の活動増大を示し、脳と筋肉が効率よく働いている象限である。第2象限は脳が筋肉よりも活動していることを示し、もっと筋肉に負荷を与えられる象限である。第3象限は脳と筋肉の活動低下を示し、十分に脳と筋肉が回復している象限である。第4象限は筋肉が脳よりも活動していることを示し、もっと脳に負荷を与えられる象限である。図9A-Eは、横軸を脳左脳の第1次運動野M1の酸素交換量のCOE変化量ΔHb-ΔHbO2、縦軸を筋肉右腕の二頭筋の酸素交換量のMOE変化量ΔHb-ΔHbO2として、同時計測して経時的にプロットして得られる二次元ダイアグラムであり、Aはダンベルを持ち上げてない場合、Bは2.5kgのダンベルを持ち上げた場合、Cは4.5kgのダンベルを持ち上げた場合、Dは7.0kgのダンベルを持ち上げた場合、Eは9.5kgのダンベルを持ち上げた場合を示し、Fは説明図である。なお、縦軸と横軸の単位はmmollである。図9Eに示すように、脳酸素消費が増加し負荷により酸素交換のBM比1運動時の脳筋肉酸素負担率は0に近くなり時計回転になる部位は、神経活動が活発化して鍛える対象となる。ここで、BM比1は、脳筋肉酸素負担率BM比1=[筋肉の酸素交換量MOE変化量]/[脳の酸素交換量COE変化量]式20で定義される。図10A-Eは、横軸を脳左脳の第1次運動野M1の周囲の酸素交換量のCOE変化量ΔHb-ΔHbO2、縦軸を筋肉右腕の二頭筋の酸素交換量のMOE変化量ΔHb-ΔHbO2として、同時計測して経時的にプロットして得られる二次元ダイアグラムであり、Aはダンベルを持ち上げてない場合、Bは2.5kgのダンベルを持ち上げた場合、Cは4.5kgのダンベルを持ち上げた場合、Dは7.0kgのダンベルを持ち上げた場合、Eは9.5kgのダンベルを持ち上げた場合を示し、Fは説明図である。なお、縦軸と横軸の単位はmmollである。図10Eに示すように、半時計周りへシフトする部位は高酸素化で、血流増加を引き起こすトレーニング対象となる。図11は、横軸を脳脳の第1次運動野M1の脱酸化型ヘモグロビンの濃度変化量ΔHb、縦軸を筋肉腕の二頭筋の脱酸化型ヘモグロビンの濃度変化量ΔHbとして、同時計測して経時的にプロットして得られる二次元ダイアグラムである。筋肉運動に対する刺激によって誘発される脳の毛細血管中の脱酸化型ヘモグロビンの濃度変化量の変化態様は、その増減の組み合わせによって以下の9パターンを示す。筋肉のHb脳のHbΔHb]ΔHb1増加増加第1象限2増加ゼロ3増加減少第4象限4減少増加第2象限5減少ゼロ6減少減少第3象限7ゼロ増加8ゼロゼロ9ゼロ減少図12は、横軸を脳脳の第1次運動野M1の脱酸化型ヘモグロビンの濃度変化量ΔHbO2、縦軸を筋肉腕の二頭筋の脱酸化型ヘモグロビンの濃度変化量ΔHbO2として、同時計測して経時的にプロットして得られる二次元ダイアグラムである。筋肉運動に対する刺激によって誘発される脳の毛細血管中の脱酸化型ヘモグロビンの濃度変化量の変化態様は、その増減の組み合わせによって以下の9パターンを示す。筋肉のHb脳のHbΔHbO2ΔHbO21増加増加第1象限2増加ゼロ3増加減少第4象限4減少増加第2象限5減少ゼロ6減少減少第3象限7ゼロ増加8ゼロゼロ9ゼロ減少図13A及びBは9kgのダンベルを持ち上げて負荷をかける課題を与えた際の脳左脳の第一次運動野と筋肉右腕の二頭筋での血液量と酸素交換量の時系列変化安静時、ダンベル負荷時、回復時を示す二次元ダイアグラムである。より詳細には、図13Aは横軸を脳左脳の第1次運動野M1の血液量のCBV変化量ΔHbΔHbO2、縦軸を筋肉右腕の二頭筋の血液量のMBV変化量ΔHbΔHbO2として、同時計測して経時的にプロットして得られる二次元ダイアグラムである。図13Bは、横軸を脳左脳の第1次運動野M1の酸素交換量のCOE変化量ΔHb-ΔHbO2、縦軸を筋肉右腕の二頭筋の酸素交換量のMOE変化量ΔHb-ΔHbO2として、同時計測して経時的にプロットして得られる二次元ダイアグラムである。ダンベルの負荷をかけた後、血液量は回復しているものの図13A参照、酸素交換量は未回復であることがわかる図13B参照。図14A-Fは、横軸を脳左脳の第1次運動野M1の酸素交換量のCOE変化量ΔHb-ΔHbO2、縦軸を筋肉右腕の二頭筋の酸素交換量のMOE変化量ΔHb-ΔHbO2として、同時計測して経時的にプロットして得られる二次元ダイアグラムである。なお、図中、D-OはΔHb-ΔHbO2を示す。図14Aは最初のトレーニング時に計測した安静中、Bは運動課題中、Cは課題後の回復期の二次元ダイアグラムを示す。トレーニングの内容は、ダンベル14.5kgを持ち上げる課題を行うことである。図14Aの安静中では、PL値計測開始からのスカラー値、距離が小さく、その最大値は、図14Bに示す運動課題中のPL値より小さいことがわかる。図14Cに示す課題後の回復期では、図14Bに示す運動課題中と同様に第1象限に座標が存在するが、BM比傾きの平均が低下し、筋肉より脳での酸素負荷が課題終了後も継続していることがわかる。図14Dはトレーニング開始後3日目に計測した安静中、Eは運動課題中、Fは課題後の回復期の二次元ダイアグラムを示す。図14D-Fに示すように、トレーニング後3日目の計測結果では、運動課題中のBM比1の平均傾きが増加し、PL値1が小さくなった。回復期では、安静中まで戻るようになり、より脳での酸素消費が少ない状態でダンベル14.5kgを持ち上げられるようになったトレーニング効果が診断できる。図15A-Fは、脳と筋肉での血液量DOと酸素交換量D-Oの時系列変化を示す二次元ダイアグラムである。図15A-Cは、横軸を脳左脳の第1次運動野M1の酸素交換量のCOE変化量ΔHb-ΔHbO2、縦軸を筋肉右腕の二頭筋の酸素交換量のMOE変化量ΔHb-ΔHbO2として、同時計測して経時的にプロットして得られる二次元ダイアグラムである。なお、図中、D-OはΔHb-ΔHbO2を示す。図15Aは安静時、Bは9kgのダンベルを持ち上げる負荷をかけている課題中、Cは回復期の二次元ダイアグラムを示す。図15D-Fは、横軸を脳脳の第1次運動野M1の血液量のCBV変化量ΔHbΔHbO2、縦軸を筋肉腕の二頭筋の血液量のMBV変化量ΔHbΔHbO2として、同時計測して経時的にプロットして得られる二次元ダイアグラムである。なお、図中、DOはΔHbΔHbO2を示す。図15Dは安静時、Eは9kgのダンベルを持ち上げる負荷をかけている課題中、Fは回復期の二次元ダイアグラムであり、それぞれ図15A-Cに時系列的に対応している。図15A及びDでは、計測開始点原点ゼロとして、脳と筋肉の安静状態をプロットしている。安静時間の軌跡のばらつきを原点からの最大距離を半径rの円として描く。すると、円の中と外で、安静時とみなす軌跡と活動中あるいは回復中とみなす時間帯が分かる。円の外にでた時間帯が活動増加時であり、ダンベル運動後でも円の中に入ってなければ、回復していないと判定できる。実際に、図15CとFとを比較すると血液量の変化DOは回復しているのに酸素交換量D-Oは未回復であることが分かる。これは、脳と筋肉の血圧および血流状態は安静時まで回復しているが、脳と筋肉のうち、脳の軸で未回復であるので、ダンベル運動を止めても脳の細胞がまだ興奮状態にあることを示している。また、図15Bでは、ダンベル運動では、脳と筋肉は右肩あがりに線形の関係で酸素交換が増強していくことを示している。図15E及びFでは、脳の血液量は上昇するが筋肉は低下したまま変化しないで、回復後筋肉の血流は即座に回復することを示している。このように、ダンベル運動で同時に筋肉と脳から血液量DOと酸素交換量D-Oの時系列変化を計測することで、運動によってその運動を止めても、脳がまだ回復過程にあるのか、を2つの指標を比べながら診断できる。脳と筋肉ではどちらが先に回復しやすいのかが分かり、次の運動負荷を与える際に、軌跡を見ながら、未回復時点で開始するのか、十分な回復をみてからか、血液量DOと酸素交換量D-Oのどちらの指標が回復した時に運動を開始するのかを判定して運動プログラムを組むことができる。酸素交換量D-Oが未回復な状態で開始すれば、同じダンベル運動でも、さらに回復が遅れたり、より強い負荷を脳あるいは筋肉にあたえるトレーニングを実行できる。一方、脳と筋肉の疲労回復を血液量DOと酸素交換量D-Oの軌跡で確認して、開始することで1回ごとの回復時間を診断できる。徐々に脳も筋肉も疲れてくるので、1回ごとの回復時間の伸びをみながら、運動継続を判断できる。例えば、5回の運動の回復時間を合計したり、1回ごとの回復時間を計算し表示して、次回の運動と比較して、運動練習の効果を診断できる。図16はダンベル運動による脳と筋肉の同時回復時間の計測を示すグラフである。安静時に脳と筋肉の軌跡をプロットして、最大半径を決めることで基準設定すると、図16に示すように、例えば、5回の運動の回復時間を合計したり、1回ごとの回復時間を計算し表示して、次回の運動と比較して、運動練習の効果を診断できる。上記の回復時間から、初回よりも3日後に回復期に効果が出ていることがわかる。図17Aは安静時の酸素交換量D-Oの変化を示す二次元ダイアグラム、Bは安静時の血液量DOの変化を示す二次元ダイアグラムである。課題後に未回復のこともあるので、図17A及びBに示すように、回復時点を最大ΔL値の二分の1ΔL/2として設定することで、回復時間を定義することもできる。図18A-Dは複数の脳部位Bと上腕二頭筋からの筋肉Mの酸素交換量D-Oの変化を示す二次元ダイアグラムである。図18の二次元ダイアグラム上には、複数の脳部位のチャンネルch1-Ch5について、計測開始点から経時的に軌道をベクトル表示されている。図18から、いつ、どの脳部位が特定の筋肉運動でもっとも酸素を使うか、安静時の状態まで回復したか、脳と筋肉との間でどちらがより酸素消費の負担が強くなったか等が二次元ダイアグラム上のBM比1ベクトルの方向と計測開始点からのスカラーパワー値:PL値1で診断できる。すなわち、BM比1が一定であれば、パワー値PL値1が大きければ大きいほど脳と筋肉に加わる酸素負荷影響が大である。象限によるがBM比1、あるいはBM角がゼロに近くなるほど脳の負担が増大する。図19A-Dは複数の脳部位Bと上腕二頭筋からの筋肉Mの血液量DOの変化を示す二次元ダイアグラムである。図19の二次元ダイアグラム上には、複数の脳部位のチャンネルch1-Ch5について、計測開始点から経時的に軌道をベクトル表示されている。図19から、いつ、どの脳部位が特定の筋肉運動でもっとも血液量が変化することで局所の血圧が変化するか、安静時の状態まで回復したか、脳と筋肉との間でどちらがより血液量変化の負担が強くなったか等が、ダイアグラム上のBM比2ベクトルの方向この場合のBM比2は、運動時の脳筋肉血液量負担率を表すと計測開始点からのスカラーパワー値:PL値2で診断できる。ここで、運動時の脳筋肉血液量負担率BM比2は、脳筋肉血液量負担率のBM比2=[筋肉の血液量MBV変化量]/[脳の血液量CBV変化量]式21で定義される。すなわち、BM比2が一定であれば、パワー値PL値が大きければ大きいほど脳と筋肉に加わる血液量変化、すなわち血圧の影響が大である。象限によるがBM比2、あるいはBM角がゼロに近くなるほど脳の負担が増大する。図20A-DはチャンネルaChaの脳部位Bと上腕二頭筋からの筋肉Mの酸素交換量D-Oの変化を示す二次元ダイアグラムである。図20の二次元ダイアグラム上には、15秒ごとにベクトル表示されている。図20Aでは、安静時の軌跡の範囲各象限の最大値がわかる。図20Bでは、負荷された課題によってD-O指標の第1象限にベクトル表示されており、脳と筋肉に酸素消費が起こっていることがわかる。図20Cでは、課題中には第4象限にベクトル表示が移行しており、脳での酸素消費が低下し、筋肉の酸素消費が増加していることがわかる。この15秒間では、脳への負荷が軽減されたことがわかる。図20Dでは、課題終了により急激に第3象限にベクトル表示が移行しており、脳と筋肉の酸素消費がともに低下したことが診断できる。図21A-DはチャンネルaChaとチャンネルbChbでの脳部位Bと上腕二頭筋からの筋肉Mの酸素交換量D-Oの変化を示す二次元ダイアグラムである。図21の二次元ダイアグラム上には、15秒ごとにベクトル表示されている。図21Aでは、安静時の軌跡の範囲各象限の最大値がわかる。図21Bでは、負荷された課題によってChaではDO指標の第1象限にベクトル表示されており、脳と筋肉に酸素消費が起こっていることがわかる。Chbは第4象限にベクトル表示されており、Chaの方が、より脳への負担が強く、Chbでは脳より筋肉への負担が強いことがわかる。図21Cに示す課題中において、Chaでは第4象限にベクトル表示が移行しており、脳での酸素消費が低下し、筋肉の酸素消費が増加していることがわかる。この15秒間では、脳への負荷が軽減された。一方、Chbでは、第1象限にベクトル表示が移行し、筋肉の酸素消費は横ばいで脳の酸素負荷が増大したことがわかる。図21Dでは、課題終了で急激に第3象限にベクトル表示が移行しており、ChaとChbについて脳と筋肉の酸素消費がともに低下したことが診断できる。図22A-CはチャンネルaChaとチャンネルbChbでの脳部位Bと上腕二頭筋からの筋肉Mの血液量DOの変化を示す二次元ダイアグラムである。図22の二次元ダイアグラム上には、15秒ごとにベクトル表示されている。図22Aでは、安静時の軌跡の範囲各象限の最大値がわかる。図22Bでは、負荷された課題によってChaではDO指標の第1象限にベクトル表示を示し、脳と筋肉に血液量の増加が起こっていることがわかる。Chaに比べてChbの方が、より脳への血液量の増加が強いことがわかる。図22Cでは、課題終了で、急激に第3象限にベクトル表示が移行した。脳の血液量は低下し、脳血圧は低下したが、筋肉は運動課題開始時よりもパンピングを起こして膨らんでいることが判断できる。次に、複数の部位からNIRS計測する上での不均一なROIの問題を軽減する手法について説明する。複数の部位からNIRS計測する上での不均一なROI関心領域:regionofinterestの問題を軽減する手法―1手法―1では、OD平面上に軌道を描くΔD脱酸化型Hbの変化量,ΔO酸化型Hbの変化量をプロットする二次元ダイアグラムを使って、ED-O、最大血液量の変化、ヘマトクリットの変化ΔHt=ΔBVBV0から毛細血管内のΔYを見積もることで、脳と筋肉などNIRSの信号変化量の違いがあっても同じ定量指標で比べることができる。E比血液量対酸素交換量比は総ヘモグロビン量DOに対する酸素交換量D-Oの割合として定義される。ROIのS/Nに依存しにくい、毛細血管内酸素飽和度の変化ΔYを計測することで部位間の比較の精度を高めることができる。図23は、e比とΔYの時系列データを利用した脳と筋肉の同時測定を説明するためのグラフである。図23において、計測開始時からある時点までのROIにおける酸素飽和度の変化ΔYは、酸化型ヘモグロビン量の変化と脱酸化型ヘモグロビン量の変化を伴う。したがって、計測開始時EoO0,D0から総ヘモグロビン量血液量の変化ΔDΔOと酸素交換量の変化ΔDΔOが起こる。そこで、総ヘモグロビン量の変化ΔDΔOと脱酸化型量の変化ΔDΔOの関係を調べるためにe比を定義して調べた。e比血液量変化対酸素交換量変化の比は総ヘモグロビンの変化量ΔDΔOに対する酸素交換量ΔDΔOの割合として定義される。e比=ΔDΔO式22ここで、k=ΔDΔOなお、E比は、段落0188に記載したように、絶対値の比として使っており、e比は、そこからの変化の比として使っており、両者を区別している。従来、頭皮からの計測によって脳実質までの距離が一定でないために入射と受光のプローブ間距離を等間隔に配置してもROI関心領域:regionofinterestにおけるヘモグロビン量光量が異なる。酸化型ヘモグロビン量Oと脱酸化型ヘモグロビン量DはROIから得られる値である。KDO式23Kは、図33Aに示すOD平面上の傾きである。K比の単位はROI中の各濃度の割合であり、ROIが複数あり、ROIの大きさが異なってもお互いに比較したり、同時に計算することが可能なパラメータである。ROI中の酸化型ヘモグロビン量Oと脱酸化型ヘモグロビン量Dは、NIRSを使ったヘモグロビン計測法の中でも、TRStimeresolvedspectroscopy時間分解法で計測される。入射と受光のプローブの1つのペアによって計測されるROIの血液量BVは、血漿成分を除けば、総ヘモグロビン量に近似される。したがって、酸化型ヘモグロビン量Oと脱酸化型ヘモグロビン量Dの和からBVDO式24と記述できる。ROIの酸化型ヘモグロビンと脱酸化型ヘモグロビンの濃度差は、OED-O式25と記述できる。そこで、ROI中の酸化型ヘモグロビンと脱酸化型ヘモグロビン濃度差と濃度和の比Eを新しく定義するとEOEBV式26K比と同様に、e比の単位はROI中の各濃度の割合であるので、ROIが複数あり、ROIの大きさが異なってもお互いに比較したり、同時に計算することが可能なパラメータである。6と23の式からEはKの関数としてE式27と記述できることが分かる。ROIの酸素飽和度をY0≦Y≦1とするとYO式281YD式29が成り立つことから式式301YK式31式6と式28、式29からEはYの関数としてE12Y式32すなわち、Y2式33が成り立つ。これより、NIRSの計測で得られた酸化型ヘモグロビン量Oと脱酸化型ヘモグロビン量Dから得られるe比によって酸素飽和度Yが求まる。このYはROIの血管内血液中の酸素濃度とみなすことができる。このことから、ROIの大きさに依存することなく、複数の計測から複数のYの時系列データが求められる。計測開始点の酸素飽和度Y0から計測開始後点の酸素飽和度Y1の変化ΔYを求める。EはYの1次関数であるので、下記のように求められる。ΔY=22=2=ΔE2式34=は以下のA,B、Cの3つの場合が起こり得る。ここで、AΔBV0の場合OD軸の二次元平面上ではD-O軸に平行な変化ΔY=式10aが成り立つΔYはΔOEの線形である。NIRS計測から得られる酸化型ヘモグロビンと脱酸化型ヘモグロビンの変化から算出される。また、安静時では、ΔBVよりもΔOEが大きく揺らぐことがこの式から理解できる。安静時ROIBV0とするとプローブの設置の仕方で、一見ΔYが変化しているように見える。図24A及びBは、安静時のΔBVとΔOEの時系列の変化を示すグラフである。図24Aに示すデータでは、安静時のΔBVとΔOEの標準偏差の比は8.1であり、Bに示すデータでは、安静時のΔBVとΔOEの標準偏差の比は1.2である。安静時のΔBVとΔOEの標準偏差の比が2.0以上のとき、ROIのHbが少なく、十分な光量が確保されていないので、S/Nが悪いと判断できる。すなわち、ΔBVとΔOEの標準偏差の比を2以下に保つようにプローブセットすることで、複数の計測部位の均一なROIを保つことができる。ここで、ROIの酸素飽和度の変化をΔY1<ΔY<1とするとΔY=OD軸の二次元平面上ではOD軸に平行な変化ΔY=式36=E0ΔD式37すなわち、OD軸に平行に変化することでΔYが変化する。ΔYは、ΔDまたはΔOに逆比例した関数である。BV0は計測対象のROIと考えられるので、図33Bに示すように、計測開始点からのΔYの時系列変化は、ΔO,またはΔDに依存して変化する。すなわち、BV0に任意の数値を代入することで、時系列変化がわかる。CΔOEとΔBVのいずれも0でない場合計測時のOとDを仮定することで、ΔYが求まる。すなわち、BV0に任意の数値を代入することで、時系列変化がわかる。ΔY=ΔY=0の場合E0E1ΔOEΔBVとなる。ここで、式4のΔYの中のBV0とΔBVの関係は、38式で表される。ΔHt=ΔBVBV0式38HtはROI中の血液量、すなわちヘマトクリットで、ΔHtはヘマトクリットの変化を表す。動脈、静脈のヘマトクリットは通常04から045である。毛細血管では、ROIが小さいほど0.2から0.8までに変化する。つまり、ΔHtは最大0.25から4までの値をとりうる。CW法で、光路長PL値が分からなくとも、ΔHtは38の式よりヘマトクリットの変化を最大ΔBVの値MaxΔBVから仮設定する方法を用いることで、ΔYの時系列テータを推定算出することができる方法を考案した。38の式から、ROI中のヘマトクリットが低下し、ΔBVが低下する場合の脳反応や筋肉反応が存在することが解釈できる。推定計算方法として、最大ΔHtMaxΔHt=ΔBVの最大変化の信号強度ROI中のBVの信号強度式39を用いて計算する。このとき、最大ΔHtを-1から1の範囲で任意選択肢し、代入することで、ROI中のBVの信号強度を算出できる。図25は上腕二頭筋の収縮時のΔOEとΔBVの時系列変化を示すグラフである。仮定1として計測開始時のY00.7、E0=-0.4仮定2として、ヘマトクリット変化=0.25とすると39式からROI中のBVの信号強度=ΔBVの最大変化の信号強度として得られる。上記より35式は、ΔY=-3e14067eΔOEΔBVとなる。Y0=0.5酸素飽和度50のときと仮定する。仮に毛細血管の酸素飽和度約50-60と見積もるDOEo=0ΔY=2=-ΔE12ΔOEΔBV1ΔHt1=-e1ΔHt1式40PL値光路長が分からなくとも、ΔHtはヘマトクリットの変化を最大ΔBVの値から推定するモデルを用いることで、ΔYの時系列テータを算出することができる方法を考案した。ΔHt=実測最大ΔBVの絶対値BV式41BV=実測最大ΔBV絶対値ΔHt毛細血管では、ROIが小さいほど0.2から0.8までに変化すなわち、ΔHtはおおよそマイナスー4から4までの値をとりうる。したがって式26の計測開始時をY005とするとΔHtへのから、eの関数によって、ΔYは変化することがわかる。図26はE値とY値を一定にして算出した筋肉運動における推定ΔYの時系列変化を示すグラフであり、横軸は時間秒、縦軸はΔYである。図26より、ヘマトクリット比がΔYに大きく影響することがわかった。この結果から、毛細血管が血球の数を変えることで、急激な酸素濃度変化に対応していることを示唆できる。一方で、ΔHt=2の仮説では、ΔY=-1.2となり、起こりえないことも推定できる。すなわち、この推定モデルでは、ΔHtの変動幅の上限も推定できる。図27はヘマトクリットを一定にして算出した筋肉運動における推定ΔYの時系列変化を示すグラフであり、横軸は時間秒、縦軸はΔYである。図27より、ヘマトクリット比が一定の時、計測開始時Y値には、ΔYに大きく影響しないことがわかった。この結果から、毛細血管が血球の数を変えることで、急激な酸素濃度変化に対応していることを示唆できる。このことは、E0、Y05として、簡略計算できることを示している。すなわち、式10はΔY=とし、酸素飽和度Yとの関係を示すグラフ、Bは、ROI中の酸化型Hbの変化量ΔO、ROI中の脱酸化型Hbの変化量ΔDとY角の変化量ΔY角との関係を示すグラフである。図33Aにおいて、上の傾きY角とすると、酸素飽和度Y=1-ArctanY角式42となる。また、図33Bにおいて、酸素飽和度Yの変化は、ΔY角に依存しているので、酸素飽和度の変化ΔY=ArctanΔY角式43として得られる。これより、酸素飽和度の変化は、酸化型Hbの変化量ΔOと脱酸化型Hbの変化量ΔDのいずれかが変化すれば起こりうる。ところで従来、fMRIでは、常磁性体であるΔDの変化を検出して、反磁性体であるΔOの変化を検出していないので、ΔDが変化しない酸素飽和度の変化は検出していなかった。また、従来のNIRSの手法では、ΔDとΔOを別々に扱ってきただけでなく、ΔOの変化を主要な指標にしてきたので、fMRIやNIRSでは酸素飽和度の変化は精度を欠いていた。図33Cに示すようにOD平面上の座標変化ΔOとΔDによって酸素飽和度の変化が決まるので、ΔDとΔOから独立して酸素飽和度を算出すると正確なΔY角が得られず計測開始点からの変化が低下と上昇で全く逆の結果となるので、酸素飽和度変化を誤診する原因となる。fMRIで用いられるΔDの変化だけでは正確な酸素飽和度の変化が得られないことを示している。複数の部位からNIRS計測する上での不均一なROIの問題を軽減する手法-2手法-2では、L値を使って2つのリアルタイムノーマリゼーションによって単位円を二次元ダイアグラム上に描くものである。ここで、D脱酸化型Hb,O酸化型Hb,ΔD脱酸化型Hbの変化量,ΔO酸化型Hbの変化量として、ΔL値とは、ΔL2脱酸化型Hbの変化量2酸化型Hbの変化量2=12ΔD-ΔO2ΔDΔO2式44で定義される。D脱酸化型Hb,O酸化型Hbが定量できる場合は、ΔL値ではなく、L値として、L脱酸化型Hb2酸化型Hb2=12D-O2DO2式45L値で補正してもよい。まず、計測部位間のL値をリアルタイムに算出して、複数点の中の最大L値を使って、任意秒ごとにノーマリゼーションを行う。または、各部位における最大L値を使って、任意秒ごとにノーマリゼーションを行う。次いで、最大L値による二次元調整補正を例えば5秒ごとに補正することで、全体を等倍率で評価できるので、一見して、どの部位が最大L値を示すチャンネルに対してXの変化リアルタイムの変化をしているかわかる。例えば、筋肉と脳のΔLmaxの比を算出して、一方をN倍して二次元ダイアグラムの傾きを求める。筋肉ΔLmax脳ΔLmax=N倍式46この場合、筋、脳を複数点で計測してもそれそれで最大L値を算出する。図34A-Cは筋肉のデータΔD、ΔOから最大L値を選んで単位円補正する手順を説明するためのグラフである。まず、データΔD、ΔOの時系列のグラフら最大L値を選択する図34A参照。次いで、ΔDΔOΔD-ΔOを算出して時系列のグラフを作成する図34B参照。次いで、横軸ΔDΔO、縦軸ΔD-ΔOとして、最大L値が半径となる単位円を作成する図34C参照。従来ROIのサイズ測定対象の大きさか異なるために、部位間に変化の違いが起こっても小さな変化がROIのサイズの違いによるものか、反応の小ささによるものか、どちらか判定することができなかった。そこで、最大L値の補正をROIごとに設定する場合と、複数のROIの中で、最大L値を選んで設定する場合とが考えられる。最大L値の補正をROIごとに設定する場合は、筋肉と脳、口腔と脳など部位の違いが大きい場合に有効である。複数のROIの中で、最大L値を選んで設定する場合は、脳内のみ、皮膚のみなど、プローブ間距離が均一でかつROIの定量が可能な場合に有効である。図35A-Cは脳のデータΔD、ΔOから最大L値を選んで単位円補正する手順を説明するためのグラフである。まず、データΔD、ΔOの時系列のグラフら最大L値を選択する図35A参照。次いで、ΔDΔOとΔD-ΔOを算出して時系列のグラフを作成する図35B参照。次いで、横軸ΔDΔO、縦軸ΔD-ΔOとして、最大L値が半径となる単位円を作成する図35C参照。単位円補正から4つのベクトル方向のうちどの方向に依存した揺らぎの違いと変化の大きさがわかる。その揺らぎの大きさから、安静時最大のΔL値と課題終了時のΔL値を重ねて描くことで、課題によって賦活化された時点が全体のか100倍することで分かる。図34のデータと図35のデータから、筋肉と脳のΔLmaxの比を算出して、脳の値をN倍して二次元ダイアグラムの傾きを求める。N倍は式46を用いる。図36Aはある部位のデータAΔD、ΔO、ΔDΔO、ΔD-ΔOの時系列のグラフ、BはデータAについての横軸ΔBVΔDΔO、縦軸ΔOEΔD-ΔOにして時系列変化課題前、課題中、課題後を示す二次元ダイアグラム、Cはある部位のデータBΔD、ΔO、ΔDΔO、ΔD-ΔOの時系列のグラフ、DはデータBについての横軸ΔBVΔDΔO、縦軸ΔOEΔD-ΔOにして時系列変化課題前、課題中、課題後を示す二次元ダイアグラムである。上記データA,BについてL値の最大値Lmaxを算出し、L値max1となるようなデータ補正振幅値の標準化を行うために、L値の逆数を算出し、この逆数を標準化の補正値とする。次いで、各指標O,D,DO,D-Oに、補正値を乗算し、標準化したデータを2D二次元平面単位円グラフに描く。図37AはデータAのLmax標準化後の2Dグラフ、BはデータBのLmax標準化後の2Dグラフである。また、BVやOEで標準化する場合にも、上記と同様の手順で行われる。ただし、この場合は単位円グラフにはならない。図38AはデータAのBVmax標準化後の2Dグラフ、BはデータAのOEmax標準化後の2Dグラフ、CはデータBのBVmax標準化後の2Dグラフ、DはデータBのOEmax標準化後の2Dグラフである。また、リアルタイムでLmaxを標準化することも可能である。まず、キャリブレーションタイム例えば3sを設定する。次いで、測定開始後、3sごとに計測開始時からのLmaxを上記と同様の手順で算出し、単位円にグラフを描く。Lmaxが更新されない場合には、同様の補正値変換でグラフがリアルタイムで表示され続ける。Lmaxが更新された場合には、それまでに表示されていたデータごと縮尺を更新し、新たな補正値変換後のグラフを表示する。図39AはデータAについてリアルタイムでLmaxを標準化して単位円を作成する手順を示すグラフ、BはデータBについてリアルタイムでLmaxを標準化して単位円を作成する手順を示すグラフである。複数の部位からNIRS計測する上での不均一なROIの問題を軽減する手法-3手法-3では、複数の部位例えば、脳と筋肉のデータの積をとる方法である。脳の測定値と腕の測定値とを二次元ダイアグラムで表示する場合、変化量や計測条件が異なるために、調整する必要がある。その1つが、脳の変化×筋肉の変化の時系列表示である。筋肉運動に対して、異なった反応を示す脳の部位を選び出すこともできる。筋肉が疲れると脳の働きが強まるが、積の最大を求めることで筋肉と脳の働きが相互に効率的に引き出せる運動やその運動の時間帯を知ることができる。トレーニングの効果を診断するのに有効である。図40Aは脳と筋肉の血液量DOの積の時系列変化を示すグラフ、Bは脳と筋肉の酸素交換量D-Oの積の時系列変化を示すグラフである。図40から、2回目の運動時では、血液量を変えないで、酸素交は換が増加していることから、脳と体の運動効率よくなっていることが分かる。各種生理的指標を複合的に可視化して画像化する手法-1ここで、本発明者は、各種生理的指標変化量やパラメータを複合的に可視化して画像化する手法をハイブリッドイメージングと呼んでいる。このハイブリッドイメージングには、脳と筋肉等複数の部位から共通の生理的指標を複合的に可視化して画像化する第1の手法と、複数の生理的指標を複合的に可視化して画像化する第2の手法とがある。図42は、ハイブリッドイメージングの第1の手法の手順を説明するためのフローチャートである。まず、本発明の実施形態例に係る生体機能診断装置を用いて、例えばO酸化型Hb変化量,D脱酸化型Hb変化量,OE酸素交換量変化量,BV血液量変化量,e比血液量変化対酸素交換量変化の比、L値事象内距離等の生理的指標の値をそれぞれ脳と筋肉から求めるステップS11。次いで、二次元上に、各生理的指標のダイアグラムを描くステップS12。次いで、上記のダイアグラムの傾き、すなわち相関係数を部位ごとに算出部10により算出するステップS13。次いで、相関係数の値を-10から10で色分けして表示部9に表示するステップ14。次いで、課題前、課題中、課題後に分けて画像化するステップS15。これによって、回復期の脳代謝と筋代謝の関係もわかる。これが、各生理的指標の相関係数を利用した脳と筋肉活動のハイブリッドイメージングとなる。図43は、各種生理的指標を説明するための二次元ダイアグラムである。複数のヘモグロビンをベクトル空間上の軌道として考え、ベクトルの向き方向とその力の大きさスカラーを考慮することで、筋肉と脳の力関係をリアルタイムに計測、画像化することができる。酸化型ヘモグロビンの変化量△Oと脱酸化型ヘモグロビンの変化量△Dをベクトル空間上で計測開始点からの波運動や回転運動のダイナミクスとして捉えることで、酸化型ヘモグロビンの変化と脱酸化型ヘモグロビンの変化の両方の変化の関係を表す複数の指標を創出することができる。図43中、△OE,△BV,△L、e比、E角は、生理的指標であり、それぞれ下の式により算出される。△OE△D△O△BV△D△O△L2△D2△O2e比=△OE△BVE角=Arctan、ch21は右脳の第1次運動野である。ch6,7,8,10,11,12はch9を囲み、ch18,19,20,22,23,24はch21を囲み、第1次運動野M1の周囲の運動野である。第1次運動野では、筋肉運動でも荷重、重さに依存して働くので酸素消費が起こるが、その周囲では血液の供給が起こっても酸素を消費しないので異なった反応になる。図45Aは図44におけるチャンネルch9のダイアグラムを示す詳細図、BはOEの経時的な移動回転をベクトル表示した図である。ここで、縦軸は筋肉上腕二頭筋の酸素交換量の濃度変化量MOE、横軸は脳第1次運動野の酸素交換量の濃度変化量COEである。各軸の単位はmmollである。図45Aからわかるように、指標OEについては、MOEとCOEの最大値の時点が異なり、MOE20秒の時点で最大がCOE36秒の時点で最大より16秒ほど早く最大値をとる。この2つの時間差は、筋肉と脳の活動性に時間的の違いがあることを示す。また、図45Bに示すベクトル表示により、運動課題中でも脳と筋肉の酸素代謝の関係が変化していることを診断することができる。図46Aは図44におけるチャンネルch10のダイアグラムを示す詳細図、BはOEの移動回転を経時的にベクトル表示した図である。ここで、縦軸は筋肉上腕二頭筋の酸素交換量の濃度変化量MOE、横軸は脳第1次運動野周囲の酸素交換量の濃度変化量COEである。各軸の単位はmmollである。図46Aからわかるように、指標OEについては、MOEとCOEの最大値の時点が異なり、MOE25秒の時点で最大がCOE36秒の時点で最大より11秒ほど早く最大値をとる。この2つの時間差は、筋肉と脳の活動性に時間的の違いがあることを示す。また、図46Bに示すベクトル表示により、運動課題中でも脳と筋肉の酸素代謝の関係が変化していることを診断することができる。図45A及び図46A中の式及びR2については、相関係数ピアソンの積率相関係数の二乗であり、指定された区間で最小二乗法によって、近似直線を引き、回帰直線を求める。相関係数rが1,-1に近いほど相関は強くなり、0に近いほど相関は弱くなる。図47は、筋肉と脳からBV血液量の生理的指標を取り出して、二次元上に複合化して描いた二次元ダイアグラムの例である。ここで、縦軸は筋肉上腕二頭筋の血液量の濃度変化量MBV、横軸は脳第1次運動野や第1次運動野周囲の血液量の濃度変化量CBVである。各軸の単位はmmollである。また、系列1は課題前の5秒間の軌跡、系列2は課題中の36秒間の軌跡、系列3は課題後の55秒間の軌跡をそれぞれ示す。ここで、課題は14.5kgのダンベルを上げることである。図48Aは図47におけるチャンネルch9のダイアグラムを示す詳細図、BはBVの経時的な移動一次線形をベクトル表示した図である。ここで、縦軸は筋肉上腕二頭筋の血液量の変化量MBV、横軸は脳第1次運動野の血液量の濃度変化量CBVである。各軸の単位はmmollである。図48Aからわかるように、指標BVについては、MBVとCBVの最大値の時点が一致している。また、図48Bに示すベクトル表示により、課題中は同じ向きを示し、課題終了とともに正反対の向きを示すことがわかる。図49は、14.5kgのダンベルを両手の持ち上げる課題の最中36秒間の右上腕二頭筋の1カ所と脳の14箇所の部位左脳のch6-ch12、右脳のch18-ch24で△OEと△BVを指標として相関係数をグラフ化したものであり、Aは筋肉△OEと脳△OEの相関係数を示し、Bは筋肉△BVと脳△BVの相関係数を示し、Cは筋肉△OEと脳△BVの相関係数を示し、Dは筋肉△BVと脳△OEの相関係数を示す。図49Aでは、△OEは脳の2箇所で相関係数0.6以上の正の相関を示しているが、ほかの部位は負の相関を示した。これは、筋肉の酸素消費の上昇と一致して脳でも酸素消費が上昇する部位はCh9とch21の第1次運動野M1で限定していると診断できる。図49Bでは、△BVは筋肉と脳の14箇所で相関係数0.8以上を示している。筋肉の血液量の増加と脳全体の血液量の増加が連動して上昇していることを示す。筋肉圧も脳圧も同時に上昇することを意味する。図49Cでは、筋肉の△OEと脳の△BVは相関係数0.6前後の値を示している。図49Dでは、筋肉の△BVと脳の△OEはCh9とch21で正の相関を示し、ほかの脳部位で負の相関を示した。これは、筋肉の充血、パンピングと連動してch9、Ch21の第1次運動野M1で酸素消費が増加することを示している。さらに、脳の部位では、酸素消費が起こらない部位でも、酸素消費を低下させて連動して活動していることを示している。図50は筋肉と脳のE角,L値を指標として相関係数をグラフ化したものであり、Eは筋肉Eと脳Eの相関係数を示し、Fは筋肉Lと脳Lの相関係数を示し、Gは筋肉Eと脳Lの相関係数を示し、Hは筋肉Lと脳Eの相関係数を示す。図50Eでは、右脳のch21と左脳のch9は、いずれも第1次運動野M1であるが、筋肉のE角と正相関したのは右脳の第1次運動野のCh21であると診断できる。筋肉の酸素消費効率は脳の酸素消費効率と強く関係していると診断される。図50Fでは、筋肉と脳のL値同士は、06以上の相関を全部で示したので、この指標では部位の特徴がなく同じように反応する。図50Gでは、筋肉のE角と脳のL値は相関性が全くなく、独立した現象であることが診断される。すなわち、脳の働きの大きさと筋肉の酸素消費効率とは関係なかったことがわかる。図50Hでは、筋肉のL値と脳のE角はCh9とch21で正の相関を示し、ほかの脳部位で負の相関を示し、筋肉の動きが、脳の酸素消費効率と関係していることが診断できる。図51は筋肉と脳のOΔO,DΔDを指標として相関係数をグラフ化したものであり、Iは筋肉△Dと脳△Dの相関係数を示し、Jは筋肉△Oと脳△Oの相関係数を示し、Kは筋肉△Dと脳△Oの相関係数を示し、Lは筋肉△Oと脳△Dの相関係数を示す。図51Iでは、筋肉と脳のDの変化量は、ch10を除き06以上の相関性を示すことがわかる。図51Jでは、筋肉と脳ではOの変化は0.6の相関性で一致性が乏しいことがわかる。図51Kでは、筋肉とDと脳のOの変化量が強く相関していると診断できる。特に、第1次運動野M1とその周囲では相関性の符号正負が異なった。図51Lでは、筋肉のOと脳のDでは、ch10のみが高い相関を示した。Ch10は、第1次運動野M1とは異なった動態を示しながら筋肉運動と関係していることがわかる。図52は、筋肉と脳における各生理的指標O,D,OE,BV,L値,E角同士の相関係数の時系列変化を-10から10で色分けして表示した図である。色分けについては、例えば赤色相関係数10、オレンジ色相関係数075、黄色相関係数05、緑色相関係数00、水色相関係数-075、青色相関係数-10のように相関係数に応じて色が変化するようにしている。図52において、各マスの中はchが7個配置されたものを可視化しており、2つ並んでいる左側のマスは左脳の7つのチャンネルch6-ch12、右側のマスは右脳の7つのチャンネルch18-ch24をそれぞれ示す。図52の図は、生体の異なった部位間の機能を同時に可視化するハイブリッドイメージングであり、脳は筋肉に指令しているだけでなく筋肉と脳のフィードバックシステムが診断できる。すなわち、筋肉と脳の相関係数の関係から、脳の複数部位と筋肉運動の指標がそれぞれ影響しあって動いていることが診断できる。このことから、筋肉の指標を見て運動量、運動パターン、時間を調節しながら行うことで、脳の部位に選択的に酸素消費を起こさせ、脳活動、脳トレーニングができる。従来の筋電図等では、筋肉が動かない時の代謝状態が評価できないが、本発明では、脳と筋肉の同時計測によって回復過程の脳と筋肉の状態評価が診断できる。なお、図52において、酸素交換量OEの図が一番顕著に変化がでている。OEの図では、脳の第1次運動野M1と筋肉の両方でOEの値が高い。図49Aのグラフを参照すると、筋肉が動くために脳には酸素を使って指令を出す部位があり、それが、第1次運動野M1であると診断できる。すなわち、オレンジ色の部位が最も酸素消費をしている脳部位であり、かつ筋肉運動と密接に関係していることがわかる。このとき第1次運動野の周囲が正反対の青色を示しているので、M1以外の脳部位は筋肉の動きと連動して、酸素を使わない正反対の動きを示したと診断できる。このように、必ずしも1つの脳部位で筋肉に動きを支えているのではなく、直接関わらない部位は間接的に脳活動をうまくいくように酸素を使わない動きをしていることがわかる。この傾向は、回復期でもある。一方、血液量BVを参照すると、ほとんど同様に高血圧を示している。E角を使った脳と筋肉の二次元ダイアグラム筋肉と脳のE角同士がなす角は、以下の式47でEハイブリッド角EH角、ダブルE角として定義される。Eハイブリッド角=Arctan=ArctanMOExCBV式47図53は、横軸を脳脳の第1次運動野M1のE角度、縦軸を筋肉腕の二頭筋のE角度として、同時計測して経時的にプロットして得られる二次元ダイアグラムである。図53に示す二次元ダイアグラムから得られる筋肉と脳のE角の比からなるEH角によって、脳と筋肉の相互関係が8つの象限に区分できて、定量的に、かつ経時的に診断、画像化できる。図53で、第1象限では、脳も筋肉もヘモグロビン変化に対する酸素交換効率が高い。第1象限Bでは筋肉が脳より、第1象限Aでは脳が筋肉より効率よく活動している。第2象限では、筋肉の活動は低下し脳が効率よく活動している第3象限では、脳も筋肉も活動が低下している。第4象限では、脳の活動は低下し、筋肉が効率よく活動している図54は筋肉と脳のそれぞれのE角を用いたch6-ch12の二次元ダイアグラムである。ここで、縦軸は筋肉のE角度、横軸は脳の筋肉度である。また、系列1は課題前の5秒間の軌跡、系列2は課題中ダンベルも持ち上げる等の36秒間の軌跡、系列3は課題後の55秒間の軌跡をそれぞれ示す。ここで、課題は14.5kgのダンベルを上げることである。図55はch6-ch12のEH角の時系列変化を示すグラフである。ここで、横軸は時間秒、縦軸はEH角度である。また、横軸に沿った太線は課題中の期間を示す。図55から、Ch9は課題中45-90度のEH角を示し、第1象限Bであり、他の部位は課題中90-135度のEH角を示し、第4象限Bであることがわかる。図56はch18-ch24のEH角の時系列変化を示すグラフである。ここで、横軸は時間秒、縦軸はEH角度である。また、横軸に沿った太線は課題中の期間を示す。図56から、Ch21は課題中45-90度のEH角を示し、第1象限Bであり、他の部位は課題中90-135度のEH角を示し、第4象限Bであることがわかる。このようにEH角は、それぞれの臓器、あるいは部位の酸素交換効率を二次元上に示した状態で、その2つの角度の相互関係を読み取り診断することができる。この筋肉と脳のE角を用いた二次元ベクトルのダイアグラム上の筋肉と脳のE角の比から、脳と筋肉の相互関係が8つの象限に区分でき、定量的に、かつ経時的に診断、画像化できる。また、象限ごとに、脳と筋肉の働きの関係性が明らかになる。また、課題の負荷が酸素効率からみて脳と筋肉に十分作用しているか否かをみて運動を変えることができる。図57Aは、横軸を脳脳の第1次運動野周囲のE角度、縦軸を筋肉腕の二頭筋のE角度として、9.5kgのダンベルをあげる前と後で同時計測して経時的にプロットして得られる二次元ダイアグラムであり、Bは、横軸を脳脳の第1次運動野のE角度、縦軸を筋肉腕の二頭筋のE角度として、9.5kgのダンベルをあげる前と後で同時計測して経時的にプロットして得られる二次元ダイアグラムである。図中、○は課題中のデータ、●は安静時と課題後のデータである。図57A及びBからわかるように、ダンベルの重量が上がると運動野とその周囲の部位では90度からマイナス90度で囲まれた範囲にプロットが集まる。特に、9.5kgのダンベル運動で、脳の運動野ではゼロと90度で囲まれた領域に集まることがわかる。課題中に90度内に入る確率が高くなったことで、90度内は多血であるので、脳圧も筋肉圧も上昇して、血液量が上がった状態で運動をしていることを示す。逆に言えば、この枠外では、まだ、それほど血液量を増やす調節反応をしなくとも対応できる負荷だったが、ダンベル9kgを持ち上げる運動では、血液量を増やす調節反応まで課題がきつくなってきて、しっかりと筋肉も脳も使えている状態にあることがわかる。各種生理的指標を複合的に可視化して画像化する手法-2図58は、ハイブリッドイメージングの第2の手法の手順を説明するためのフローチャートである。まず、本発明の実施形態例に係る生体機能診断装置を用いて、例えばO酸化型Hb変化量,D脱酸化型Hb変化量,OE酸素交換量変化量,BV血液量変化量,e比血液量変化対酸素交換量変化の比、L値事象内距離等の生理的指標の値をそれぞれ脳と筋肉から求めるステップS21。次いで、算出部110により各生理的指標の乗算を部位ごとに行うステップS22。次いで、乗算値を時系列に表示部に表示するステップS23。次いで、部位ごとの乗算値を値の大きさで色分けして表示部9に表示するステップS24。次いで、例えばOE同士の積算値とBV同士の積算値をOE-BVの二次元ダイアグラム上に描くステップS25。これが、各生理的指標の乗算値を利用した脳と筋肉活動のハイブリッドイメージングとなる。筋肉と脳のOE,BVの積を使った2次元ダイアグラム診断脳と筋肉を同時に計測したことで、OEとBVの両者の積によって示される最大値を示すことができる。OEの積では、右上腕二頭筋を動かすことに寄与している運動領域が明確になる乗算したOEとBVのなす角は、以下の式48でハイブリッド角と定義でされる。ハイブリッド角式48ハイブリッド角H角は、脳と筋肉が同時に働いている最中の酸素交換効率を表し、ゼロから90度までより効率よく働いかを表す。図59は、14.5kgのダンベルを持ち上げる課題を行った際における筋肉腕の二頭筋と左脳の第1次運動野及び周囲ch6-ch12の血液量の濃度変化量BV同士の積MBVCBV及び酸素交換量の濃度変化量OE同士の積MOECOEの時系列変化を示すグラフである。横軸は時間秒、縦軸はBV、OEの積である。横軸に沿った太線は課題中の期間を示す。図60は、14.5kgのダンベルを持ち上げる課題を行った際における筋肉腕の二頭筋と右脳の第1次運動野及び周囲ch18-ch24の血液量の濃度変化量BV同士の積MBVCBV及び酸素交換量の濃度変化量OE同士の積MOECOEの時系列変化を示すグラフである。横軸は時間秒、縦軸はBV、OEの積である。横軸に沿った太線は課題中の期間を示す。図61は筋肉のMBVとMOEの時系列変化を示すグラフである。横軸は時間秒、縦軸はMBV、MOEの値である。横軸に沿った太線は課題中の期間を示す。図62は、14.5kgのダンベルを持ち上げる課題を行った際における筋肉腕の二頭筋の酸素交換量の濃度変化量OEと右脳の第1次運動野及び周囲ch18-ch24の血液量の濃度変化量BVの積MOECBVの時系列変化を示すグラフである。横軸は時間秒、縦軸はMOECBVの積である。横軸に沿った太線は課題中の期間を示す。図62では、すべてのchから、筋肉のOEと脳のBVが運動中上昇して線形性を示して正の相関を示しており、筋肉運動が脳全体の脳血液量を増やすことがわかる。図63は、14.5kgのダンベルを持ち上げる課題を行った際における筋肉腕の二頭筋の血液量の濃度変化量BVと右脳の第1次運動野及び周囲ch18-ch24の酸素交換量の濃度変化量OEの積MBVCOEの時系列変化を示すグラフである。横軸は時間秒、縦軸はMBVCOEの積である。横軸に沿った太線は課題中の期間を示す。図63では、第1次運度野M1のCh21について筋肉の血液量MBVと脳の酸素交換量COEが運動中に線形性を示し正の相関を示しており、第1次運度野M1の周囲は負の相関を示した。図64は、14.5kgのダンベルを持ち上げる課題を行った際における筋肉腕の二頭筋と左脳の第1次運動野及び周囲ch6-ch12の血液量の濃度変化量BV同士の積MBVCBVを横軸、酸素交換量の濃度変化量OE同士の積MOECOEを縦軸にして、経時的にプロットして得られる二次元ダイアグラムである。なお、図中、pretaskは課題前の5秒間の軌跡、taskは課題中の36秒間の軌跡、posttaskは課題後の55秒間の軌跡を示す。図64では、左脳の第1次運動野のch9におけるハイブリッド角H角が課題中に最大酸素交換効率を示し。ch9の周囲はむしろ負の酸素交換効率を示し、酸素交換が低い関係を示している。また、2つの臓器が連動して働いているとき最も効率の良い酸素消費が行われた時間20秒がわかる。筋肉と脳のOE,BVの積を使った二次元ダイアグラムによる診断脳と筋が消費した負荷量生理的な酸素消費負荷量MOECOEあるいは、血液負荷量MBVCBVが、生体の中枢と末梢を総合したエネルギー負荷量を表すので、これらをリアルタイムで表示する可視化には十分な意義がある。1乗算されたOEと乗算されたBVの増加ピークが同時でない、2筋肉のOEと、脳のOEの増加ピークが同じでない、という結果も明らかである。この診断には次のような効果を奏する。1時系列的には、課題中のBVとOEの増加曲線が第1次運動野M1左脳のch9または右脳のch21で一致せず、2つが独立的に調節されているということBVよりOEが早く増加することがわかる。2筋肉のOE増加が、脳よりも早く、筋肉のOE増加が、脳のOE増加のトリガーになっているかもしれないということ等が、同時に計測してわかる。乗算によるイメージング画像化をすることで、脳と筋肉が消費した負荷量生理的な酸素消費負荷量が診断できる。3脳と筋肉の二次元平面上のBVxOEは、ベクトルの外積面積で、脳と筋肉が出力したパワーが面積が得られる。D-Oなら、脳と筋の両方のOE負荷量を反映した中枢と末梢を総合したOE負荷量が乗算したことで新たに診断できる。図65Aは、14.5kgのダンベルを持ち上げる課題を行った際における筋肉腕の二頭筋と左脳の第1次運動野及び周囲ch6-ch12の血液量の濃度変化量BV同士の積MBVCBVの時系列変化を色分けして表示した図、Bは酸素交換量の濃度変化量OE同士の積MOECOEの時系列変化を色分けして表示した図である。図65において、各マスの中はchが7個配置されたものを可視化している。色分けについては、例えば図65Aでは、赤色積の値1000、オレンジ色積の値750、黄色積の値500、緑色積の値00、水色積の値-500、青色積の値-1000のように積の値の大きさに応じて色が変化するようにしている。図65Bでは、例えば赤色積の値2000、オレンジ色積の値1500、黄色積の値1000、緑色積の値00、水色積の値-1000、青色積の値-2000のように積の値の大きさに応じて色が変化するようにしている。図65の図は、生体の異なった部位間の機能を同時に可視化するハイブリッドイメージングである。図66Aは、14.5kgのダンベルを持ち上げる課題を行った際における筋肉腕の二頭筋と左脳の第1次運動野及び周囲ch6-ch12のL値同士の積の時系列変化を色分けして表示した図、Bはハイブリッド角H角の時系列変化を色分けして表示した図である。図66において、各マスの中はchが7個配置されたものを可視化している。色分けについては、例えば図66Aでは、赤色積の値2500、オレンジ色積の値1750、黄色積の値1250、緑色積の値00、水色積の値-1250、青色積の値-2500のように積の値の大きさに応じて色が変化するようにしている。図66Bでは、例えば赤色角度90、オレンジ色角度65、黄色角度45、緑色角度00、水色角度-45、青色角度-90のようにH角の角度の大きさに応じて色が変化するようにしている。図66の図は、生体の異なった部位間の機能を同時に可視化するハイブリッドイメージングである。図66Aに示すように、L値同士の値が高い部位が、第1次運動野M1に一致していて、かつ運動を継続すると徐々に上昇し。終了すると一気に平坦化するので、脳と筋か同時に活動していることとどれだけ強く活動していることがわかる。30秒と示された画像かピークをしめしている。図66Bに示すように、ハイブリッド角H角は、脳と筋が酸素を消費してもっとも効率よく働いている時間を選び出すことで診断することができる。この場合は運動開始後、20秒である事がわかる。すなわち、H角は異なった部位、臓器同士が連動して働く際の酸素交換効率を定量的に診断、測定する指標として優れている。プログラム図1に示す本発明の実施形態例に係るプログラム12は、上記生体機能診断装置Kの制御部7の上記の処理を実行させることを特徴とするものである。このプログラム12は、磁気ディスク、CD-ROM、半導体メモリ等の記録媒体に記録されていてもよく、通信ネットワークを介してダウンロードされるものでもよい。本発明は、上記実施の形態に限定されることはなく、特許請求の範囲に記載された技術的事項の範囲内において、種々の変更が可能である。例えば、明細書及び図面に記載された各種生理的指標変化量、パラメータの種類、調整部11による生理的指標の調整の程度、比較調整すべき生理的指標同士の組み合わせ、生理的指標同士の相関係数、生理的指標同士の乗算の組み合わせ等は例示であり、これに限定されるものではない。また、生体の生体機能の異なる複数の部位とは、調整部11によって生理的指標同士を調整しなければならないほど生体機能が異なっている部位であり、筋肉と脳だけでなく、脳の異なる部位、筋肉の異なる部位例えば腕と足等も含まれる。さらに、3以上の部位の生理的指標を算出、調整し、三次元の立体的な画像で表示部9に表示してもよい。"