(57)【要約】【課題】代謝機能指数を非採血的に測定するサーモグラフィシステムと糖尿病のサーモグラフィ臨床診断用ルーラ（Ｒｕｌｅｒ　）の表示装置及び表示方法を提供する。【解決手段】物体表面温度分布のデジタルデータを収録するサーモグラフィ装置を使って、特定環境条件に置いた生体末梢循環部位に冷水負荷を掛ける前後でのサーモグラフィ画像を複数枚取り、幾つかの仮設条件で生体伝熱理論から体内代謝状況と血流状況を推測し、空間領域と時間領域での生体末梢部の伝熱方程式を解析することにより生体内部代謝系の働きを表現する一連のパラメータを算出し、システム同定技術で末梢部の有効熱伝導率と代謝生熱率と言う代謝系における重要な二つのパラメータを持って代謝機能指数を測定し、モニタ画面上に統計的な方法で糖尿病臨床診断用の代謝機能指数を表示する二次元非線形的なルーラを提供する。【選択図】　　　　図１
2004016734.tif 000002
【特許請求の範囲】
【請求項１】
被測定対象物からの赤外線をサーモグラフィ装置を介して検出して、該被測定対象物から温度測定を行なって代謝機能を表示する代謝機能表示方法に於いて、
座標軸上に血液注入率と代謝機能指数をパラメータとし、二次画面上に血糖値を中心とする糖尿病総合判定で得られた統計データに基づいて糖尿病を型分けした非線形的な曲線を表示させて成ることを特徴とするサーモグラフィ装置による代謝機能表示方法。
【請求項２】
被測定対象物から赤外線をサーモグラフィ装置を介して検出して、該被測定対象物から温度測定を行なって代謝機能を表示する代謝機能表示装置に於いて、
座標軸上に血液注入率と代謝機能指数をパラメータとし、二次画面上に血糖値を中心とする糖尿病総合判定で得られた統計データに基づいて糖尿病を型分けした非線形的な曲線を表示させて成ることを特徴とするサーモグラフィによる代謝機能表示装置。
【請求項３】
被測定対象物からの赤外線をサーモグラフィ装置を介して検出して、該被測定対象物の温度測定を行なって代謝機能を表示する代謝機能表示方法に於いて、
上記被測定対象物に温度負荷を掛ける前後の時刻に上記サーモグラフィ装置で撮像した温度画像データを記憶手段に格納させるステップと、
上記記憶手段に格納した温度画像データから上記被測定対象物の内部発熱状態を伝熱方程式を用いて推測し、該伝熱方程式を分解、解析して代謝系の働きを表現する一連のパラメータをコンピュータを介して算出する演算ステップと、
上記演算ステップで得た上記内部発熱分布から代謝機能指数パラメータとなる有効熱伝導率と代謝生熱率とを計測するステップとから成り、
上記代謝機能指数と血液注入率とをパラメータとし、二次画面上に血糖値を中心とする糖尿病総合判定で得られた統計データに基づいて糖尿病を型分けした非線形的な曲線を表示させて成ることを特徴とするサーモグラフィ装置による代謝機能表示方法。
【請求項４】
被測定対象物からの赤外線をサーモグラフィ装置を介して検出して、該被測定対象物の温度測定を行なって代謝機能を表示する代謝機能表示装置に於いて、
上記被測定対象物に温度負荷を掛ける前後の時刻に上記サーモグラフィ装置で撮像した温度画像データを格納させる記憶手段と、
上記記憶手段に格納した温度画像データから上記被測定対象物の内部発熱状態を伝熱方程式を用いて推測し、該伝熱方程式を分解、解析して代謝系の働きを表現する一連のパラメータをコンピュータを介して算出する演算手段と、
上記演算手段で得た上記内部発熱分布から代謝機能指数パラメータとなる有効熱伝導率と代謝生熱率とを計測する計測手段とを具備し、
上記代謝機能指数と血液注入率とをパラメータとし、二次画面上に血糖値を中心とする糖尿病総合判定で得られた統計データに基づいて糖尿病を型分けした非線形的な曲線を表示させて成ることを特徴とするサーモグラフィ装置による代謝機能表示方法。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明はサーモグラフィ装置による代謝機能表示方法及びその表示装置に係わり、特に非採血的な方法で代謝系の働きを表現する一連のパラメータを測定して、糖尿病のサーモグラフィ臨床診断用のルーラを得て、代謝機能を表示する表示装置及び表示方法に関する。
【０００２】
【従来の技術】
糖尿病の臨床診断は従来から、血糖値を計測している。糖尿病は代表的な代謝系の病気である。代謝は生体が活きる上で不可欠的な重要の機能である。特に膵臓はインスリンを分泌し、血液や体液中に放出する。
【０００３】
このインスリンはさらに葡萄糖を分解して、生体の生命活動に要するエネルギーを供給する。分解されない葡萄糖は血液や尿のなかに残り、体外に排出される。血液のなかの葡萄糖の量は血糖値と言う。普通、血糖値で分類する糖尿病には主に１型糖尿病と２型糖尿病と言う２種類があり、１型糖尿病はインスリンが絶対的に不足であり、２型糖尿病はインスリンが相対的に不足している。この分類法は本質的に膵臓のインスリン分泌能力を評価している。体内にインスリンが不足すると、各種の器官的な或いは組織的な病気に罹り易く、各種の糖尿病の併発症が起こる。
【０００４】
実際に臨床上で使用している糖尿病の血糖値の診断では血糖値を通して、間接的にインスリンを測定することに過ぎないため、また血糖値が食事や生理状態や運動など左右されているため、血糖値で血液中のインスリン量を評価する精度が低く、糖尿病の血糖値診断法の正確率は低い。
【０００５】
例えば、１９８５年に世界保健機構（ＷＨＯ）は糖尿病の血糖値診断基準を提案したが、１９９７年にアメリカの糖尿病専門委員会はその修正案を提唱した。しかし、統計によると、二つの提案とも糖尿病診断の正確率は凡そ６０％前後である。
【０００６】
又、サーモグラフィ装置で生体温度を測定する方法は従来から種々提案されている。しかし、それらの測定方法は生体内部に腫瘍など単なる局所発熱機構異常の場合にサーモグラフィ画像のパターン処理によって直接的に判断する場合に利用される。
【０００７】
生体の代謝系には葡萄糖の分解とともに熱が発生する。それによって、糖尿病診断にもサーモグラフィ装置を利用して、生体内部の代謝や末梢血管循環などの障害を診断することができる。
【０００８】
従来技術に於いて、サーモグラフィ装置を用いた糖尿病診断では温度回復率を用いたものがある。温度回復率（ＴＲと記す）は普通、下記（１）式のように計算する。
ＴＲ＝（Ｔｉ　−Ｔ１　）／（Ｔｏ　−Ｔ１　）　　　　（１）
ここにＴｏ　は温度負荷直前の定常状態温度であり、Ｔ１　は温度負荷直後の温度であり、Ｔｉ　は温度負荷後ｉ　時刻の温度である。
【０００９】
生体内部の生熱機構はかなり複雑である。温度回復パターンのみでは簡単に内部状況を判断することができない。この様なＴＲの検査方法では糖尿病の正確な臨床診断ができず、代謝率などの測定もできない。従って、臨床上、糖尿病を診断するには多回採血して、血液中の血糖値の測定やインスリン抵抗試験などの生化学的な方法にたよっていのが現状である。
【００１０】
また、物体表面の温度から物体内部の生熱機能を非接触的に診断するには数学的なモデルが必要であり、従来では逆問題解析用の数学的なモデルを作って、内部熱力学パラメータを計算していた。しかし、その逆問題解析法は色々問題があって、実用に至っていない。
【００１１】
物体内部の熱力学方程式は理論上では種々作られている。生体の熱測定には従来ペネス（Ｐｅｎｎｅｓ）の伝熱方程式が利用されている。Ｐｅｎｎｅｓの伝熱方程式の一般的な表現は次式（２）で示される。
ρｃ∂Ｔ／∂ｔ＝∇・（Ｋ∇Ｔ）＋Ｗｂ　Ｃｂ　ρｂ　（Ｔａ　−Ｔｖ　）Ｑｂ　＋Ｑｍ　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　（２）
ここに、ρは物質比重、ｃ　は物質比熱、∂は微分子、Ｔは温度関数、∇は勾配子、Ｋは物質の伝熱率定数、Ｗｂ　は血液注入率、Ｃｂ　は血液比熱、ρｂ　は血液比重、Ｔａ　は動脈血温度、Ｔｖ　は静脈血温度、Ｑｂ　は血流変数、Ｑｍ　は代謝生熱率である。
【００１２】
この方程式を解くには物体内部各位置の温度を計測することと各物質定数や変数などを算出することが必要であるが、生体内部の温度は非侵入的な手段で測定することが望ましいが、現状では生体内部に侵入的な温度計測方法を用いて、一つ一つの定数を計測しなければならない。又、Ｐｅｎｎｅｓ伝熱方程式を用いてそのままで代謝生熱率などを計測することもできないし、生体代謝機能指数などの代謝系のパラメータを非侵入的な方法で指定することもできない。
【００１３】
生体代謝率などを表示するのは数字的な表示方法とヒストグラフ的な表示方法があるが、本発明で作成した非線形的な表示手段（ルーラ）を用いた特殊な表示に基づいて直観的に表示する方法は提案されていない。
【００１４】
【発明が解決しようとする課題】
従来技術では上述の様に糖尿病を診断する場合に血糖値を使うので精度も低いし、採血を必要とし、診断結果の表示も直観的でない。生体代謝機能を表すこともできない。そのために精度が高く、採血も必要でなく、直観的に診断結果を表示する生体代謝機能指数を計測するシステムの開発が本発明の解決しようとする課題である。
【００１５】
本発明はサーモグラフィ装置を利用し、代謝系にインスリンによる葡萄糖分解時に放出した熱量に注目して、非採血的なサーモグラフィ温度測定方法で、生体内部の代謝系の働きを表現する代謝機能指数と言うパラメータを測定し、糖尿病の臨床診断に新しい参照基準を提供し、表示手段の画面上に二次元の非線形的な糖尿病サーモグラフィ臨床診断用の代謝機能指数ルーラを作り、糖尿病の型分けと重症度を直観的に表示できるシステムを提供するものである。具体的に主に次の要件を解決する。
▲１▼　糖尿病の型分けと重症度を直観的に表示する。
▲２▼　生体内部代謝系の動きを強制的に誘発し、温度分布データとして記録する。
▲３▼　サーモグラフィのデータから生体代謝系における代謝機能指数を計測する。
【００１６】
【課題を解決するための手段】
本発明は被測定対象物１からの赤外線をサーモグラフィ装置２を介して検出して、被測定対象物１から温度測定を行なって代謝機能を表示する代謝機能表示装置及び表示方法に於いて、座標軸上に血液注入率と代謝機能指数をパラメータとし、二次画面上に血糖値を中心とする糖尿病総合判定で得られた統計データに基づいて糖尿病を型分けした非線形的な曲線を表示させて成ることを特徴とするサーモグラフィ装置による代謝機能表示装置及び表示方法としたものである。
【００１７】
本発明のサーモグラフィ装置による代謝機能表示方法及びその表示装置によれば、表示手段の画面上に二次元パラメータを座標軸として作った糖尿病サーモグラフィ臨床診断用の代謝機能指数ルーラ（表示装置）によって臨床診断の主な目的である糖尿病の型分けと重症度を直観的に判断可能なものが得られる。
【００１８】
【発明の実施の形態】
以下、システムブロック及び計測時のフローチャートに基づいて本発明のサーモグラフィ装置による代謝機能表示方法及びその表示装置を図１及び図２に基づいて詳記する。
図１は本発明に用いるシステムブロック図、図２は計測時のフローチャートを示すものである。図１に於いて、１は被測定対象物であり、この被測定対象物１は所定の環境温度設定領域５に配された環境に於いて、測定が成される。
【００１９】
被測定対象物１はサーモグラフィ装置２によって撮像されて、被測定対象物１から発する表面の赤外線温度を測定する。
【００２０】
サーモグラフィ装置２内にはＰＣカード等の記憶手段６を有する。サーモグラフィ装置２で計測された温度データはコンピュータ・システム（以下ＣＰＵと記す）３に供給される。ＣＰＵ３内には通常のＬＣＤ等の表示装置４やＲＯＭ８、ＲＡＭ９等の記憶手段を有し、計測された温度データはＲＡＭ９等に格納され、計測された温度データを基に解析、演算が行なわれる。
【００２１】
環境温度設定領域５に配設された被測定対象物１は恒温冷水槽７中に配された所定温度の冷水に基づいて、冷水負荷を掛けることで被測定対象物１の内部温度を表面に誘出させる様な処理が行なわれる。
【００２２】
次に図２に示すフローチャートに基づいて、本発明の概略的なプロセス動作を説明する。
図２に於いては、被測定対象物１を生体とし、糖尿病の型分けと重症度の判別をサーモグラフィ装置を用いて計測する場合について説明する。
【００２３】
先ず、第１ステップＳ１　に示す様に被測定対象物１となる生体内部代謝系の動きを生体末梢部位に温度負荷を掛ける手段で強制的に体表面へ誘発する。但し、温度負荷用の冷水温度は外部環境温度に対しての差が最適値に成る様に選択して固定し、温度負荷の時間も固定する。体表面に出てくる温度情報は冷水負荷直前と冷水負荷後の多時刻に於いて、サーモグラフィ装置２で記録し、装置２内蔵のＰＣカード６に収容する。
【００２４】
かくすれば第２ステップＳ２　に示す様に生体の例えば、手指中心温度が自動的に検出される。
【００２５】
次の第３ステップＳ３　乃至第５ステップＳ５　では従来のＰｅｎｎｅｓ方程式に基づいて、幾つかの仮定条件をおいてから簡略した伝熱方程式を作る。さらに、伝熱方程式を空間領域と時間領域に分解して解析して各時刻に於ける空間上のパラメータの計算（第３ステップＳ３　）を行い、次に体表面の温度分布から体内の温度分布を推測する。具体的には指の平均温度の時間領域のパラメータの計算（第４ステップＳ４　）が行なわれる。また、システム同定技術で代謝機能指数と言う末梢部の有効熱伝導率と代謝生熱率とを計測（第５ステップＳ５　）する。
【００２６】
次の第６ステップＳ６　では臨床で血糖値を中心とする糖尿病総合診断で得られた統計的なデータに基づいて表示装置４の画面上に二次元の非線形的な曲線で糖尿病のサーモグラフィ臨床診断用の代謝機能指数ルーラを作成する。このルーラを用いて、糖尿病の型分けと重症度を直観的に表示することで、サーモグラフィ装置を使って、代謝機能指数の非採血的な測定システムと糖尿病のサーモグラフィ臨床診断用の代謝機能指数ルーラ（表示装置）を提供可能となる。
【００２７】
上述のフローチャートで説明したサーモグラフィ装置による代謝機能指数の測定方法と糖尿病用の代謝機能指数ルーラに依れば、
【００２８】
▲１▼　冷水負荷を掛けてから生体内部の熱環境を強制的に変更することで代謝機能と共に血流状況も温度刺激によって誘発され、代謝系の活動状況が容易に観測される。生体内部代謝機能が体表面温度状態に相関にがあれば、サーモグラフィ装置で非接触的に計測した体表面の温度情報は体内代謝機能情報が含まれている。本発明ではサーモグラフィ装置を使うので採血する必要がなく、サーモグラフィ装置で撮った温度分布データはデジタルでＰＣカード６に収容してあるから、ＣＰＵ３での処理に便利である。また、冷水温度と環境温度の差及び温度負荷時間を固定することは夏も冬も同じ負荷量と考えることができる。
【００２９】
▲２▼　生体伝熱理論に基づいた熱方程式により、生体内部の熱変化を体表面の温度変化から推測することができる。空間領域における熱方程式の解析により、生体内部の代謝系の働きを表現する一連のパラメータを算出することができる。システム同定は最小二乗誤差法で代謝機能指数としての末梢部の有効熱伝導率と代謝生熱率とのパラメータを測定することができる。
【００３０】
▲３▼　更に、表示装置４の画面上に二次元パラメータを座標軸として作成した糖尿病サーモグラフィ臨床診断用の代謝機能指数ルーラにより、臨床診断の主な目的である糖尿病の型分けと重症度を直観的に判断することができる。
【００３１】
以下、本発明の一形態例を図３乃至図１１によって詳記する。
まず本発明を実施する一般的な外部条件と環境温度について説明する。
サーモグラフィ装置で温度を測定する場合、一般的な外部条件が必要である。
例えば、図３に示す検査室１０内で測定する場合に厚い黒色カーテンで屋外の赤外線を遮断して、ドアを閉め、室内の環境温度が目標値に達したら、クーラ１１や赤外線を発するランプなどを切り、人が移動しないようにしてから測定を開始する。図３は測定室の配置を示している。被測定対象物１の生体を心身ともに安定させる状態で、まず１５分前後の環境温度に順化させる。それから両手１３が心臓と略等しい高さの手置台１２に置き、真上から手背に対してサーモグラフィ装置２で一枚の温度画像を撮影する。続いて両手を恒温冷水槽７の冷水に１分間沈めて温度負荷を掛ける。それから数分ごとに各一枚の温度画像を撮影する。
【００３２】
環境温度については夏と冬に対応し易くするために、通常状態の２２℃年平均室温に対する計測室の環境温度を夏は２５±２℃に冬は２０±２℃とに設定し、計測中の開始時刻と終了時刻の温度変動差は夏と冬で共に±０．　５℃以内に収める。測定対象物の生体の温度順化時間は夏は１０分間であり、冬は２０分間である。
【００３３】
温度負荷用の冷水温度の最適値は環境温度と１５℃の差を採り、夏に１０±１℃であり、冬に５±１℃に設定する。但し、実用上、厳密に上述で決めた環境温度を守ることは難しいので、温度補償を付けてパラメータを計算する。この場合には大きな誤差が出るから、できれば厳密に決めた環境条件で始めからやり直す。
【００３４】
次にサーモグラフィ装置２で代謝機能指数を測定する原理とその熱伝達方程式について説明する。
生体代謝系の膵臓はインスリンを分泌し、血液や体液の中に放出する。このインスリンはさらに葡萄糖を分解し、葡萄糖の分解に伴い熱を生成し、この生成された熱が体表面へ伝達し、体表面の温度として現れる。人体の局所温度は主に血流や局所代謝によるものである。全身代謝による熱エネルギーは血流によって末梢部位の温度に影響し、末梢循環の局所代謝もその部位の温度に影響する。
【００３５】
そして人体内部の、特に微循環部の代謝機能は末梢部位である手の指でその局所温度を判断できる。手の指はサイズが小さく、円柱体モデルとして最適である。
【００３６】
従って測定実施部位には手の指を選び、手の背部位の温度を測定する。手の指内には動静脈が並設している、血流や代謝の動的な状態情報が含まれている。温度測定点は手の背部位の中指の指先中心から手の背部位中央部まで計測する。図４には人の中指１４のサーモグラフィ上の温度計測部位を示している。
【００３７】
熱伝達方程式を作成するために人の手の指の組織を図５の円柱体モデルで考える。この中指モデルで中指１４の中心線１５の複数個所のポイントＰ１　〜Ｐ６　で測定を行なう。手の指の温度分布（Ｔ）　は空間位置の関数であり、経時的に変化する時間の関数である。手の指の温度は心臓から血液注入した熱エネルギーと末梢循環の代謝によって発生した熱エネルギーで生成したものである。また、環境との熱交換で皮膚表面の熱放射により熱エネルギーは発散するものもあり、生体の内部温度と動的な平衡状態になる。
【００３８】
上述の中指モデルの末指部位の中指１４に対して、Ｐｅｎｎｅｓ方程式は血流による熱伝達と組織の熱伝達とを一括的に考えて、有効熱伝導率Ｋｅｆｆ　を定義し、以下の仮設条件が置かれ、θ次元の解析を省略し、円柱座標系で表現すると（３）式の簡略的な方程式が立てられる。
【００３９】
（ａ）指組織は均等媒質と等方向熱特性を持つ物質であり、∂Ｔ／∂θ＝０であり、Ｋｅｆｆ　とＱｍ　は空間的な定数である。
（ｂ）体内の中心温度は定数である。
2004016734.tif 000003
ここにρとＣはそれぞれ指組織の物質比重と比熱であり、Ｋｅｆｆ　は有効熱伝導率であり、Ｑｍ　は末梢部の代謝生熱率である。
【００４０】
（３）　式の空間領域での図５の指の径方向（ｒ）の一般解が（４）　式であり、体表面の境界条件を入れると（５）　式が得られる。
2004016734.tif 000004
【００４１】
（５）　式により、定常状態の指の径ｒ方向の温度分布は距離の自乗に比例することを利用して、（６）　式も成り立つ。
Ｔ（ｒ）＝ａ０　＋ａ１　ｒ＋ａ２　ｒ２　　　　（６）
【００４２】
中心軸１５方向には主に心臓からの血液が末梢へ注入され、皮膚表面に注目し、定常状態の温度分布のみ考えると、中心軸１５方向Ｘの温度分布として（７）　式が導かれる。
Ｔ（ｘ）＝Ｋｘ　＋ｄ　　　　（７）
【００４３】
冷水負荷を掛けると、末梢部の温度が降下し、血管が収縮し、血液注入が減少する。それで、末梢部の温度バランスは崩れ、体温を維持するためには代謝が活発化し、血液の注入も時間的に段々回復し、（３）　式中のＫｅｆｆ　とＱｍ　も変化する。
【００４４】
そして、冷水負荷の終了後に末梢部の温度Ｔ　はＫｅｆｆ　とＱｍ　の変化によって調節される。この場合、解析を簡略するために中指１４の中心軸１５のＸ軸方向の温度分布が線形的な変化である（７）　式から分るように、径方向のＴ（ｒ，ｔ）のみ考えばよい。すると、（３）　式のＫｅｆｆ　とＱｍ　と（６）　式のパラメータａも時間的な関数として考えれば良い。
【００４５】
従って、最終的には体表面の温度分布に対して、（８）式が得られる。
2004016734.tif 000005
さらに、（８）　式は体表面と周囲環境の熱交換を考えると、（９）　式が得られる。
2004016734.tif 000006
（９）　式ではαが体表面と周囲環境の熱交換係数であり、Ｔａ　が環境温度であり、Ａが放射面積であり、εが放射率であり、σがボルツマン定数である。
【００４６】
（８）　式と（９）　式を時間的に離散化すれば（１０）式と（１１）式になる。
2004016734.tif 000007
2004016734.tif 000008
それで、サーモグラフィ装置２で冷水負荷前後における生体末梢部位の各時刻の温度を記録し、環境温度と体温を測定しておけばＫｅｆｆ　とＱｍ　が求められる。
【００４７】
実測の結果を見ると、冷水負荷後に末梢部の代謝生熱率Ｑｍは（１２）式のような関数で表現できる。
ａｅ−ｂｔ　＝Ｑｍ０−Ｑｍ　（ｔ）　　　（１２）
（１２）式のａは定常時の代謝率Ｑｍ０及び冷水負荷直後の代謝率Ｑｍ　で決定され、ｂは冷水負荷撤去後の定常状態に回復する速さを表現している。
【００４８】
Ｋｅｆｆ　はＱｍ　と同じような傾向にある。（３）　式のＫｅｆｆ　の定義からも分るように冷水負荷後の末梢部の熱源は代謝と血液注入の二つであり、Ｋｅｆｆ　は血液注入率を表現し、全身の代謝による生熱効率も表現する。すると、ａとｂが代謝系の働きを表すパラメータである。これで、システム同定技術により、最小二乗法で（１２）式のパラメータａとｂを求めるために（１２）式を（１３）式に変形する。
ｌｎａ−ｂｔ＝ｌｎ〔Ｑｍ０−Ｑｍ　（ｔ）〕　　　（１３）
即ち：Ｙ（ｔ）＝Ａ＋Ｂｔ　　　（１４）
ここに、Ｙ（ｔ）＝ｌｎ〔Ｑｍ０−Ｑｍ　（ｔ）〕、Ａ＝　ｌｎａ、Ｂ　＝ｌｎｂであり、ここでＡとＢを代謝機能指数と定義する。
【００４９】
また、本発明では上式で求めた代謝機能指数ａ及びｂを座標軸とし、表示手段４の画面上に表現し、臨床で血糖値を中心とする糖尿病総合診断で得られた統計的なデータに基づいて非線形的な曲線でルーラ（Ｒｕｌｅｒ　）を作成し、実測のａ　とｂ　の座標値を位置付ければ生体の代謝機能が観測され、糖尿病の臨床診断に糖尿病の型分けと重症度を直観的に判断することができる。
【００５０】
上述の（６）式の指の中心軸の二次元的な分布Ｔ（ｒ）を証明するためにＰｅｎｎｅｓが人の前腕で実測した実測値（℃）と、本発明の（７）式の同定値と計算値とを比較すると下記の表１に示す結果が得られた。図６に表１のデータを縦軸に温度（℃）を横軸に測定点半径ｒをプロットしたものである。
【００５１】
【表１】
2004016734.tif 000009
【００５２】
また、本発明に依り、臨床的に収集した糖尿病患者のデータを用いて、図７に示してある冷水負荷後の温度回復曲線Ｔ（ｒ）も二次元的な分布を証明している。
【００５３】
図７（Ａ）、（Ｂ）は縦軸に温度（℃）を、横軸に冷水負荷時間をとったもので図７（Ａ）は冷水負荷状態を示す曲線を、図７（Ｂ）は冷水負荷後の所定時間経過時の統計的な経時温度分布を示し、曲線１６は二次元的分布を示している。
【００５４】
また、臨床的に実測した一人の糖尿病患者のデータを用いて、本発明の有効熱伝導率Ｋｅｒｒ　と代謝生熱率Ｑｍ　を所定時刻毎に実測した結果を下表２に示す。図８はこれらのデータのグラフを示している。
【００５５】
【表２】
2004016734.tif 000010
【００５６】
表２に於いて、時刻０−は冷水負荷前の時刻であり０＋は冷水負荷後の時刻を表す。また、図８（Ａ）は縦軸に温度（℃）を横軸に温度負荷前後の時刻（分）をとったものであり、図８（Ｂ）は縦軸にＫｅｆｆ　を、図８（Ｃ）は縦軸にＱｍ　をとっている。これらのＴ，Ｋｅｆｆ　及びＱｍ　の曲線１７，１８，１９は常に一定の傾向の曲線を示し、サーモグラフィ装置で糖尿病を診断できるパラメータは代謝機能指数の有効熱伝導率や代謝生熱率を用い得ることを示している。
【００５７】
次に、（１４）式で定義した血液注入率と代謝機能指数のパラメータＡとＢを座標軸として画面上にベゼール関数で臨床的に血糖値を中心とする糖尿病総合診断で得られた統計的なデータに基づいて糖尿病を型分けした非線形的な曲線を図９及び図１０の様に、表示手段の表示パターン曲線２０〜２３となるルーラを作成する。曲線２０から右側の領域が糖尿病の無い健常者の分布領域を示し、曲線２１乃至２３間が２型糖尿病患者の分布領域を示し、曲線２３から左側領域が１型糖尿病患者の分布領域を示している。又、２４は左側に行くに従って糖尿病が重くなる重症度合いを示す矢印である。
【００５８】
従って、本発明のサーモグラフィ装置による代謝機能表示方法及びその装置に依れば被検者の血液注入率と代謝機能指数を測定、演算して、図９に示す表示手段のルーラ上にＣＰＵのキーボード等の入力手段を介して入力させることで、二次元非線形的な糖尿病分類尺度から直観的に糖尿病の類型と健常者を見分けることが可能となる。
【００５９】
図１０は１５人の健常者と９４名の糖尿病患者の測定結果をルーラ上にプロットしたもので、健常者が二次元非線形的な曲線２０〜２３に入ることはないことが証明されている。
【００６０】
図１１の直線２５は縦軸に実測血糖値を横軸に代謝機能指数の自乗和をとった場合の相関関係を示すものであり、これらは血糖値の側面から、本発明の有効性を実証していることが解る。
【００６１】
【発明の効果】
本発明のサーモグラフィ装置の代謝機能表示方法及びその装置によればサーモグラフィを使っているので、採血する必要がなく、サーモグラフィ装置が取得した温度分布データはコンパクトなＰＣカード（ＣＦメモリカード）にデジタル記録してあるので、コンピュータで簡単に処理することができる。また冷水負荷により、生体内部の熱環境を強制的に変更し、代謝機能と血流状況が温度刺激によって誘発され、代謝系の活動状況が容易に観測される。生体伝熱理論に基づいた伝熱方程式により、末梢循環部位で内部の熱状況が体表面の温度状況と繋がり、体表面の温度から体内の生熱状況を推測することができる。空間領域と時間領域における伝熱方程式の解析により、生体内部の代謝系の働きを表現する一連のパラメータを算出することができる。更にシステム同定は代謝機能指数としての末梢部の有効熱伝導率と代謝生熱率とのパラメータを測定することもできる。またモニタ画面上に二次元パラメータを座標軸として作成した糖尿病のサーモグラフィ臨床診断用の代謝機能指数を表示するルーラにより、臨床診断の主な目的である糖尿病の型分けと重症度を直観的に判断することができる。
【図面の簡単な説明】
【図１】本発明のサーモグラフィ装置の代謝機能表示装置のシステムブロック図である。
【図２】本発明のサーモグラフィ装置の代謝機能表示装置の表示の為のフローチャート図である。
【図３】本発明に用いる測定室の略線的配置図である。
【図４】本発明のサーモグラフィ装置の代謝機能表示装置による温度計測部位の平面図である。
【図５】本発明のサーモグラフィ装置の被測定対象物となる人体の手の指の円柱形モデルを示す図である。
【図６】前腕の実測温度と最小自乗法での温度の同定値計算温度との関係を示す曲線図である。
【図７】冷水負荷後の温度回復曲線を示す図である。
【図８】実測した糖尿病患者の熱パラメータの分布曲線図である。
【図９】本発明の糖尿病の診断用スケールとなるルーラの表示パターンである。
【図１０】糖尿病診断結果の分布図である。
【図１１】代謝機能指数自乗和と実測血糖値との相関を示す図である。
【符号の説明】
１‥‥被測定対象物、２‥‥サーモグラフィ装置、３‥‥コンピュータ（ＣＰＵ）、４‥‥表示装置、６‥‥ＰＣカード（記憶手段）、７‥‥恒温冷水槽、８‥‥ＲＯＭ、９‥‥ＲＡＭ
【図１】
2004016734.tif 000011
【図２】
2004016734.tif 000012
【図３】
2004016734.tif 000013
【図４】
2004016734.tif 000014
【図５】
2004016734.tif 000015
【図６】
2004016734.tif 000016
【図７】
2004016734.tif 000017
【図８】
2004016734.tif 000018
【図９】
2004016734.tif 000019
【図１０】
2004016734.tif 000020
【図１１】
2004016734.tif 000021