JP2001276057A - 終末呼気速度定数と動脈血速度定数との関係式中の定数決定方法およびキセノンｃｔ装置 - Google Patents

Abstract

(57)【要約】 【課題】脳血流量を正確に求めることを可能とする。 【解決手段】Ｋ−Ｓ演算手段１２４は、仮定値設定手段
１１４で設定したＫｅ／Ｋａ変換定数γの仮定値γαに
基づいてＫｅ／Ｋａ変換手段１１２で得られた仮算出値
Ｋａｉα、Ｋａｏαから、脳／血液分配係数λの仮算出
値λαを求める。この場合、仮定値設定手段１１４は、
仮定値γαを所望の範囲で変化させている。γ値確定手
段１２６は、仮算出値λαと目標λ値算出手段１２２か
らの目標値λτとを比較することによって、仮算出値λ
αが目標値λτに最も近づくときの仮定値γαを、真の
Ｋｅ／Ｋａ変換定数γとして確定する。

Description

【発明の詳細な説明】

【０００１】

【発明の属する技術分野】この発明は、終末呼気速度定
数と動脈血速度定数との関係式中の定数決定方法および
キセノンＣＴ装置に関し、一層詳細には、キセノンガス
供給装置とＸ線ＣＴ装置等を利用して、脳の血流量を正
確に測定可能とする終末呼気速度定数と動脈血速度定数
との関係式中の定数決定方法およびキセノンＣＴ装置に
関する。

【０００２】

【従来の技術】Ｘ線ＣＴ装置により被検体としての、例
えば、患者の頭部の断層画像を得ながら、ガス吸入装置
から送出されるキセノンガスと酸素ガスとの混合ガスを
呼吸用マスクを通じて前記患者に一定時間吸入させた
後、通常の空気を呼吸させたときの前記断層画像を解析
して患者の頭部の血流を測定する方法が知られている。

【０００３】すなわち、この測定方法によれば、前記混
合ガスが、患者の肺から肺静脈中に吸収され心臓を経由
して動脈血流として頭部の組織に流れ込み、その頭部組
織を経由し静脈血流を通じて心臓にもどされ、該心臓を
経由して肺動脈にもどされる。このときの、前記頭部の
組織中におけるキセノンガス濃度の時間的変化をＸ線Ｃ
Ｔ装置により観察し、組織が正常である頭部のキセノン
ガス濃度の時間的変化と比較することで、前記患者の頭
部診断を行うことができるようになっている。

【０００４】この測定方法を用いて脳の血流量を得るた
めには、脳の組織中におけるキセノンガス濃度ととも
に、動脈中のキセノンガス濃度が必要となる。近時、こ
の動脈中のキセノンガス濃度として、非侵襲的な方法で
検出が可能な終末呼気中のキセノンガス濃度を代用して
いる。

【０００５】

【発明が解決しようとする課題】ところで、本出願の発
明者は、動脈中のキセノンガス濃度と終末呼気中のキセ
ノンガス濃度との相関関係を見出し、これを明らかにし
ている｛「The Effect of Xenon Inhalation Speed on
Cerebral Blood Flow Obtained Using the End-Tidal M
ethod in Xenon-Enhanced CT」、Shigeru Sase、Journa
l of Computer Assisted Tomography、22(5):786-791、
1998参照｝。

【０００６】この相関関係とは、動脈中のキセノンガス
濃度の速度定数が、変換定数を用いて、終末呼気中のキ
セノンガス濃度の速度定数の一次指数関数で表されると
いうものである。

【０００７】この発明は、このような知見を考慮してな
されたものであり、動脈中の血流のキセノンガスの濃度
と、終末呼気中のキセノンガスの濃度との相関関係を利
用することによって、終末呼気中のキセノンガスの濃度
を用いて脳の血流量を正確に求めることを可能とする、
終末呼気速度定数と動脈血速度定数との関係式中の定数
決定方法およびキセノンＣＴ装置を提供することを目的
とする。

【０００８】

【課題を解決するための手段】この発明に係る、キセノ
ンＣＴ検査に係る終末呼気速度定数をＫｅとし、動脈血
速度定数をＫａとしたときの関係式 Ｋａ＝γ×（１−ｅｘｐ（−Ｋｅ／γ）） 中の定数γを決定する方法は、キセノンＣＴ画像上で関
心領域を設定するＡステップと、設定した関心領域でキ
セノン分配係数λが所定の目標値に最も近づく前記定数
γを求めるＢステップとを備えている（請求項１記載の
発明）。このように、キセノン分配係数λを指標として
定数γを求めることによって、該定数γの正確な値を得
ることができる。

【０００９】また、前記Ｂステップでは、前記定数γを
０．２４から７．７の間の所望の範囲で変化させて前記
分配係数λを計算し、前記分配係数λが前記目標値に最
も近づく前記定数γを求めている（請求項２記載の発
明）。このように、考慮すべき定数γの範囲を限定する
ことによって、定数γを求めるための処理時間を短縮す
ることができる。

【００１０】この場合、前記所望の範囲を、０．３から
２．５までの範囲としてもよい（請求項３記載の発
明）。

【００１１】また、前記Ｂステップでは、前記関心領域
に含まれる所定画素毎に、前記分配係数λ値が前記目標
値に最も近づく前記定数γを求め、得られた定数γ値を
平均して目的の定数γ値としている（請求項４記載の発
明）。このため、定数γを一層正確に求めることができ
る。

【００１２】さらに、前記Ａステップでは、前記関心領
域を脳の白質を含む領域に設定し、前記Ｂステップで
は、前記目標値をヘマトクリット値に応じて決定してい
る（請求項５記載の発明）。このため、脳の血流量を求
めるための定数γを正確に求めることができる。

【００１３】この発明に係るキセノンＣＴ装置は、被検
体にキセノンガスを供給するガス供給装置と、前記被検
体の終末呼気中のキセノンガスの濃度（以下、呼気キセ
ノンガス濃度と記す。）を測定するための濃度測定装置
と、前記被検体の被検査部位のキセノンガスの濃度（以
下、被検査部位キセノンガス濃度と記す。）を得るため
に前記被検査部位のＣＴ画像データを取得するＸ線ＣＴ
装置本体と、前記ＣＴ画像データに基づいて前記被検査
部位キセノンガス濃度を求めるとともに、該被検査部位
キセノンガス濃度および前記呼気キセノンガス濃度に基
づいて前記被検査部位の血流量を求めるデータ処理装置
とを備えており、前記データ処理装置は、前記呼気キセ
ノンガス濃度の速度定数（以下、呼気速度定数と記
す。）を動脈中の血流のキセノンガスの濃度（以下、動
脈キセノンガス濃度と記す。）の速度定数（以下、動脈
血速度定数と記す。）に変換するための変換定数γに基
づいて、前記被検査部位と前記被検体の血液との間のキ
セノンガスの分配係数λを求めるとともに、該分配係数
λが所定の目標値に最も近づくときの前記変換定数γを
真値として確定するようにしている（請求項６記載の発
明）。このように、分配係数λを指標として変換定数γ
を求めることによって、該変換定数γの正確な値を得る
ことができる。

【００１４】この場合、前記データ処理装置は、前記変
換定数γを求める変換定数設定手段を備え、前記変換定
数設定手段は、前記変換定数γの仮定値を設定するとと
もに、該仮定値を変化させる仮定値設定手段を有する
（請求項７記載の発明）。

【００１５】また、前記変換定数設定手段は、前記変換
定数γの仮定値に基づいて、前記呼気速度定数から前記
動脈血速度定数の仮算出値を求める仮速度定数設定手段
と、前記動脈血速度定数の仮算出値と前記被検査部位キ
セノンガス濃度とから、前記分配係数λの仮算出値を求
める仮分配係数算出手段と、前記変換定数γの各仮定値
に基づいて得られた前記分配係数λの各仮算出値のう
ち、前記目標値に最も近づくものに対応する前記変換定
数γの仮定値を、前記真値として確定すべき抽出仮定値
として抽出する変換定数抽出手段とを有する（請求項８
記載の発明）。

【００１６】さらに、前記変換定数設定手段は、前記抽
出仮定値のうち、該抽出仮定値に対応する前記分配係数
λの仮算出値が所定の濾過範囲内に含まれるものを、前
記真値として確定すべき濾過値として求めるフィルタリ
ング手段を有する（請求項９記載の発明）。

【００１７】そして、前記変換定数設定手段は、前記Ｃ
Ｔ画像データに含まれる各画素のデータのうち、所定の
複数の画素に対応するデータに基づいて得られた前記被
検査部位キセノンガス濃度をそれぞれ抽出するＲＯＩデ
ータ抽出手段と、抽出された前記被検査部位キセノンガ
ス濃度に基づいてそれぞれ得られた前記濾過値の平均値
を前記変換定数γの真値として求める平均値算出手段と
を有する（請求項１０記載の発明）。

【００１８】このため、変換定数γを一層正確に求める
ことができる。

【００１９】また、前記仮定値設定手段は、前記変換定
数γの仮定値を０．２４から７．７の間の所望の範囲で
変化させている（請求項１１記載の発明）。このよう
に、考慮すべき仮定値の範囲を限定することによって、
変換定数γを求めるための処理時間を短縮することがで
きる。

【００２０】この場合、前記所望の範囲を、０．３から
２．５までの範囲としてもよい（請求項１２記載の発
明）。

【００２１】また、前記変換定数設定手段は、前記目標
値を決定する目標値算出手段を有し、前記目標値算出手
段は、前記被検査部位を前記被検体の脳とし、かつ、前
記変換定数γを求めるための関心領域を前記脳の白質を
含む領域に設定したとき、前記目標値をヘマトクリット
値に応じて決定する（請求項１３記載の発明）。これに
より、脳の血流量を求めるための変換定数γを正確に求
めることができる。

【００２２】この場合、前記変換定数γに基づいて前記
呼気速度定数から前記動脈血速度定数を求めるための関
係式は、前記呼気速度定数をＫｅとし、前記動脈血速度
定数をＫａとしたとき、 Ｋａ＝γ×（１−ｅｘｐ（−Ｋｅ／γ）） で表される（請求項１４記載の発明）。

【００２３】また、キセノンＣＴ装置は、前記血流量お
よび／または前記分配係数λの分布図を表示する表示装
置を有する（請求項１５記載の発明）。

【００２４】

【発明の実施の形態】以下、この発明の一実施の形態に
ついて図面を参照して説明する。

【００２５】図１は、この一実施の形態の脳血流量測定
装置１０の全体的構成を示している。また、図２は、図
１に示す脳血流量測定装置１０のブロック図を示してい
る。

【００２６】図１および図２において、脳血流量測定装
置１０は、人等の被検体１２に対するキセノンＣＴ検査
を行うためのキセノンＣＴ装置であり、基本的には、被
検体１２の断層映像（キセノンＣＴ画像）を得るＸ線Ｃ
Ｔ装置１４と、被検体１２にキセノン（Ｘｅ）と酸素
（Ｏ２）の混合ガスを供給する混合ガス供給装置１６
とから構成される。

【００２７】Ｘ線ＣＴ装置１４は、Ｘ線ＣＴ装置本体１
８と、該Ｘ線ＣＴ装置本体１８を制御するとともに混合
ガス供給装置１６を制御する制御装置２０とから構成さ
れる。制御装置２０はＸ線ＣＴ装置本体１８により得ら
れた画像データ等を処理するデータ処理装置としても機
能する。なお、制御装置２０を、Ｘ線ＣＴ装置本体１８
を制御する制御装置と、混合ガス供給装置１６を制御す
る制御装置とに物理的に分離して構成することもでき
る。

【００２８】図１に示すように、Ｘ線ＣＴ装置本体１８
には、被検体１２を載せた状態で矢印ＡＢ方向に移動さ
れる移動テーブル２２が上面に配置された被検体載せ台
２４と、円筒状の開口２６が形成されたガントリ３２と
が備えられている。ガントリ３２には、その円筒状の開
口２６の回りを例えば矢印ａ方向に旋回するように構成
されたＸ線管２８（図２参照）と、前記開口２６の回り
の円周上に配置された複数の検出器からなる検出器３０
（図２参照）とが配されている。

【００２９】図１に示すように、混合ガス供給装置１６
は、キセノンガスボンベ３６と酸素ガスボンベ３８と、
これらからのキセノンガスおよび酸素ガスを内部のコン
ピュータ４０の制御のもとに混合する吸入装置本体４２
と、この吸入装置本体４２に一端側が接続され、かつ他
端側が呼吸用マスク４４に接続される導管４６とを有し
ている。

【００３０】この場合、導管４６は、吸気管４６ａと呼
気管４６ｂおよび呼吸用マスク導管４６ｃとから構成さ
れている。そして、呼吸用マスク４４にはキセノンガス
濃度測定センサ（濃度測定装置）４８が取り付けられ、
その濃度測定センサ４８の検出信号がコンピュータ４０
に供給され、コンピュータ４０により呼気中のキセノン
ガスの濃度が計算されるようになっている。

【００３１】混合ガス供給装置１６の全体動作を制御す
るコンピュータ４０は、制御装置２０と電気的に接続さ
れ、相互に通信を行うように構成されている。

【００３２】Ｘ線ＣＴ装置１４の制御装置２０は、図２
に示したように、制御装置及び処理装置として機能する
コンピュータ５０を有し、このコンピュータ５０により
Ｘ線ＣＴ装置本体１８と混合ガス供給装置１６の動作を
制御するとともに、ガントリ３２内の検出器３０により
検出された、被検体１２の被検査部位としての脳５４の
断層画像を構成する画素データを処理し、該断層画像等
を作成する。

【００３３】コンピュータ５０には、さらに、マウス５
１（図１参照）やキーボードを有する操作コンソール５
２と、光磁気ディスク装置や磁気ディスク装置等の外部
記憶装置５５と、カラーＣＲＴ等の表示装置５６とが接
続されている。

【００３４】図１および図２例の脳血流量測定装置１０
において、実際上、操作コンソール５２は、表示装置５
６の画面上に表示され、マウス５１により操作されるマ
ウスポインタにより画面上の該当表示をクリックするこ
とにより該当表示が示す処理の実行を指示する。

【００３５】前記表示装置５６上には、後述するよう
に、コンピュータ５０による処理を通じて、Ｘ線ＣＴ装
置本体１８により得られたＣＴ画素データにより表され
る脳５４の断層画像（いわゆるＣＴ画像）等がカラーあ
るいはモノクロームで表示されるとともに、脳血流量の
画像が表示される。また、表示装置５６の画面上に表示
された画像は、制御装置２０に内蔵されたプリンタによ
りプリントアウトし、カラーまたはモノクロームのハー
ドコピー５７（図１参照）として出力することができる
ようになっている。

【００３６】次にこの実施の形態の動作について、図３
に示すフローチャートに基づき説明する。フローチャー
トの制御主体はコンピュータ５０である。

【００３７】まず、ステップＳ１において、医師等のオ
ペレータが操作コンソール５２を操作し、図４に示すよ
うに、被検体載せ台２４上に被検体１２を載せた状態で
移動テーブル２２を矢印Ｂ方向に移動させ、被検体１２
の脳５４の断層画像が撮影できる位置で停止させる。

【００３８】次いで、ステップＳ２において、図４に示
したように、呼吸用マスク４４が被検体１２の口と鼻部
を覆うように取り付けられる。

【００３９】そして、ステップＳ３では、図４に示した
測定可能状態において操作コンソール５２が操作される
ことで、制御装置２０のコンピュータ５０からの測定開
始指令が混合ガス供給装置１６のコンピュータ４０およ
びＸ線ＣＴ装置本体１８にそれぞれ送られる。

【００４０】このとき、まず、ステップＳ４において、
Ｘ線ＣＴ装置本体１８による脳５４の断層画像、いわゆ
るベースラインＣＴ画像が撮影されて外部記憶装置５５
に取り込まれる。

【００４１】次に、混合ガス供給装置１６のコンピュー
タ４０による制御の下に、キセノンガスボンベ３６と酸
素ガスボンベ３８から送出されるキセノンガスと酸素ガ
スとが吸入装置本体４２によりキセノンガス３０％酸素
７０％の割合で混合され、吸気管４６ａ、呼吸用マスク
導管４６ｃおよび呼吸用マスク４４を通じて被検体１２
の肺に供給され、かつ被検体１２の肺を通じて排出され
る呼気ガスが呼吸用マスク４４、呼吸用マスク導管４６
ｃおよび呼気管４６ｂを通じて吸入装置本体４２にもど
される。

【００４２】このとき、ステップＳ５において、被検体
１２に対する混合ガスの供給開始時点から、混合ガス中
のキセノンガスの濃度が所定値（この場合、３０％）と
なるようにコンピュータ４０により混合ガス供給装置１
６が制御されて測定が開始されることで、吸入過程、い
わゆるＷａｓｈ−ｉｎが開始される。

【００４３】なお、混合ガス供給装置１６として、例え
ば、この出願人による特公平３−３３３２６号公報に開
示された装置を用いることができ、呼気中のキセノンガ
スの濃度測定が測定開始時点から例えば４０ｍｓ毎に行
われる。

【００４４】そして、ステップＳ６においては、被検体
１２に対する混合ガスの吸気開始時点から、後述する飽
和判断、洗い出し過程、いわゆるＷａｓｈ−ｏｕｔ過程
への切り替え等の処理を行いながら約６０ｓ毎にガント
リ３２中のＸ線管２８からＸ線が被検体１２に対して放
射され、被検体１２を通過したＸ線が検出器３０により
検出されることで脳５４の断層映像が約６０ｓ毎に撮像
され、ＣＴ画素データとしてコンピュータ５０に取り込
まれる。

【００４５】次いで、ステップＳ７において、このＣＴ
画素データからＣＴ値（すなわち、ハウンスフィールド
ユニット［ＨＵ］）が各画素毎に抽出され、さらに、こ
のＣＴ値に基づき脳組織中のキセノンガス濃度が各画素
毎に計算される。なお、この実施の形態において画素の
大きさは約０．５ｍｍ角としているが、適当な大きさに
変更することが可能である。

【００４６】また、各画素のキセノンガス濃度は、移動
平均法を用いて計算される。すなわち、複数個（７×７
個、９×９個、１１×１１個等であり、好ましくは、９
×９個）の画素で構成された測定領域から各画素毎にキ
セノンガス濃度を求め、さらに、測定領域全体における
これらキセノンガス濃度の平均値を、測定領域の例えば
中心に位置する画素のキセノンガス濃度として算出す
る。そして、測定領域を１画素単位または複数の画素単
位（例えば、測定領域単位である９画素単位）の幅で移
動させながら、各画素のキセノンガス濃度を算出する。

【００４７】なお、発明の理解を容易にするために、以
下の説明において、「画素毎に」とあるのは、便宜上、
脳５４を構成する「組織毎にあるいは各組織毎に」と同
じ意味であるものとする。

【００４８】このようにして、図５に示すように、被検
体１２の被検査部位のキセノンガス濃度としての、脳５
４の各組織毎のキセノンガス濃度（脳キセノンガス濃
度）Ｃｂ（Ｔ）［ｍｇ／ｇ］｛Ｔは測定時点［ｓ］で、
Ｔ＝０、６２、１２３、１８３、２４３、３０３、３６
３、４２４、４８４、５４４、…等約６０ｓ（１ｍｉ
ｎ）毎の値をとる。｝を得ることができる。

【００４９】次いで、ステップＳ８において、この脳キ
セノンガス濃度Ｃｂ（Ｔ）の増加の割合が、予め定めて
ある所定値より小さくなった場合には、飽和状態になっ
たものと判断し、続くステップＳ９において、洗い出し
過程が済んでいるかどうかの判断をした後、ステップＳ
１０において、混合ガスの供給を中止し、混合ガスの代
わりに通常の空気を送る、いわゆる洗い出しを行う。

【００５０】さらに、ステップＳ６の処理を約１ｍｉｎ
（６０ｓ）毎、ステップＳ７の処理を４０ｍｓ毎に行い
ながら、洗い出し過程が終了した後には（ステップＳ９
の判断が肯定的になった後には）、ステップＳ１１にお
いて、脳キセノンガス濃度Ｃｂ（Ｔ）が所定値以下にな
ったことが確認されるまで、同様に約１ｍｉｎ毎にステ
ップＳ６、４０ｍｓ毎にステップＳ７の処理を行う。そ
して、脳キセノンガス濃度Ｃｂ（Ｔ）が所定値以下にな
ったとき（ステップＳ１１：ＹＥＳ）、ステップＳ１２
において、キセノンガス濃度測定センサ４８で求めた呼
気の濃度データと脳５４の脳キセノンガス濃度Ｃｂ
（Ｔ）とに基づき、以下に説明するように、脳血流量を
算出する。

【００５１】そして、続くステップＳ１３においては、
算出結果等に基づき表示装置５６上に後述する各種の表
示を行う。

【００５２】次いで、脳５４の組織毎の脳血流量を算出
するためのステップＳ１２の処理について詳しく説明す
る。

【００５３】まず、脳血流量を算出するためのアルゴリ
ズムについて説明する。なお、以下の説明においては、
任意の測定時点ｔ［ｓ］における、終末呼気中のキセノ
ンガスの濃度（呼気キセノンガス濃度）をＣｅ（ｔ）
［ｍｇ／ｍｌ］とし、動脈（この場合は、肺動脈以外の
動脈であり、例えば、頸動脈である。）中の血流のキセ
ノンガスの濃度、すなわち、脳５４に流れ込む動脈血流
のキセノンガスの濃度（動脈キセノンガス濃度）をＣａ
（ｔ）［ｍｇ／ｍｌ］とし、動脈キセノンガス濃度Ｃａ
（ｔ）の速度定数（動脈血速度定数）をＫａ［ｍｉ
ｎ^-1］とし、呼気キセノンガス濃度Ｃｅ（ｔ）の速度定
数（終末呼気速度定数）をＫｅ［ｍｉｎ^-1］とし、脳血
流量をｆ［ｍｌ／ｇ／ｍｉｎ］とし、脳／血液分配係数
（被検体１２の脳５４と血液との間のキセノン分配係
数）をλとする。

【００５４】次の（１）式は、この一実施の形態におい
て脳血流量ｆを求めるために用いられるＫｅｔｙ−Ｓｃ
ｈｍｉｄｔの式である。

【００５５】 Ｃｂ（Ｔ） ＝ｆ×∫Ｃａ（ｔ）×ｅｘｐ（−ｆ×（Ｔ−ｔ）／λ）ｄｔ … （１） ここで、右辺の定積分の変域は［０，Ｔ］である。

【００５６】前記（１）式中の左辺の脳キセノンガス濃
度Ｃｂ（Ｔ）は、測定時点ＴにおいてＸ線ＣＴ装置１４
で取得されたＣＴ画像データに含まれるＣＴ値に基づい
て求められている。

【００５７】一方、前記（１）式中の右辺の定積分の解
は、次の手順で求められる。

【００５８】図６は、実験データ等に基づいて得られた
動脈キセノンガス濃度Ｃａ（ｔ）の特性を示している。
この図６中、時点ｔ＝０〜ｗの期間は、吸入過程（Ｗａ
ｓｈ−ｉｎ）の期間であり、時点ｔ＝ｗ〜９の期間は、
洗い出し過程（Ｗａｓｈ−ｏｕｔ）の期間である。

【００５９】動脈キセノンガス濃度Ｃａ（ｔ）の時間的
変化（測定時点ｔに対する変遷）の特性は、最小自乗法
により、Ｗａｓｈ−ｉｎ期間では例えば次の（２）式の
一次指数関数で近似することができ、Ｗａｓｈ−ｏｕｔ
期間では例えば次の（３）式の一次指数関数で近似する
ことができることが知られている。

【００６０】 Ｃａ（ｔ）＝Ａａ×（１−ｅｘｐ（−Ｋａｉ×ｔ）） … （２） Ｃａ（ｔ）＝Ｃａｗ×ｅｘｐ（−Ｋａｏ×（ｔ−ｗ）） … （３） ここで、Ａａは所定の定数であり、Ｃａｗは時点ｔ＝ｗ
における動脈キセノンガス濃度Ｃａ（ｔ）の値である。
また、Ｋａｉは、Ｗａｓｈ−ｉｎ期間における動脈血速
度定数Ｋａであり、Ｋａｏは、Ｗａｓｈ−ｏｕｔ期間に
おける動脈血速度定数Ｋａである。

【００６１】前記（２）式および（３）式中の動脈血速
度定数Ｋａｉ、Ｋａｏは、呼気キセノンガス濃度Ｃｅ
（ｔ）の呼気速度定数Ｋｅｉ、Ｋｅｏ（Ｋｅｉは、Ｗａ
ｓｈ−ｉｎ期間における呼気速度定数Ｋｅであり、Ｋｅ
ｏは、Ｗａｓｈ−ｏｕｔ期間における呼気速度定数Ｋｅ
である。）から求めることができる。特に、この一実施
の形態においては、本出願の発明者によって明らかにさ
れた、動脈血速度定数Ｋａｉ、Ｋａｏと呼気速度定数Ｋ
ｅｉ、Ｋｅｏの関係式｛次の（４）式および（５）式｝
によって求められる。

【００６２】 Ｋａｉ＝γ×（１−ｅｘｐ（−Ｋｅｉ／γ）） … （４） Ｋａｏ＝γ×（１−ｅｘｐ（−Ｋｅｏ／γ）） … （５） ここで、γは、Ｋｅ／Ｋａ変換定数であり、肺における
ガス交換能を表すパラメータである。

【００６３】そして、前記（４）式および（５）式中の
呼気速度定数Ｋｅｉ、Ｋｅｏは、濃度測定センサ４８で
得られる呼気キセノンガス濃度Ｃｅ（ｔ）のデータ｛図
６に示す動脈キセノンガス濃度Ｃａ（ｔ）と類似の時間
的変化の特性を有する。｝から最小自乗法によって得ら
れる例えば次の（６）式および（７）式からそれぞれ求
められる。

【００６４】 Ｃｅ（ｔ）＝Ａｅ×（１−ｅｘｐ（−Ｋｅｉ×ｔ）） … （６） Ｃｅ（ｔ）＝Ｃｅｗ×ｅｘｐ（−Ｋｅｏ×（ｔ−ｗ）） … （７） ここで、Ａｅは所定の定数であり、Ｃｅｗは時点ｔ＝ｗ
における呼気キセノンガス濃度Ｃｅ（ｔ）の値である。

【００６５】すなわち、動脈血速度定数Ｋａｉ、Ｋａｏ
は、Ｋｅ／Ｋａ変換定数γを特定することによって、実
測値である呼気キセノンガス濃度Ｃｅ（ｔ）に基づいて
求めることができる。

【００６６】ここで、Ｋｅ／Ｋａ変換定数γの値がとり
得る範囲を特定するための手順について説明する。

【００６７】図７は、下田らによって開示された、呼気
速度定数Ｋｅと動脈血速度定数Ｋａとの関係を記録した
実験データ（記号「□」で示す。）を表すグラフである
（「Discrepancy of xenon concentrations between en
d-tidal and blood collection methods in xenon-enha
nced computed tomographic measurement of cerebral
blood flow」、M.Shimoda, et al.、Neuroradiology、3
5:66-68、1992参照）。

【００６８】また、図８は、前記（４）式および（５）
式（この場合、Ｋａ＝Ｋａｉ＝Ｋａｏとする。）の特性
を、γ＝０．５、１、２、４、８の場合毎に表したグラ
フである。

【００６９】そして、図９は、図７のグラフに図８のグ
ラフ｛この場合、γ＝０．３、２．５（実線で示
す。）、０．２４、７．７（破線で示す。）｝が重ね合
わされたグラフである。

【００７０】図９に示すように、記号「□」の位置は、
破線で示すγ＝０．２４および７．７のときの前記
（４）式および（５）式の特性を表す曲線に挟まれた領
域に全て含まれており、特に、実線で示すγ＝０．３お
よび２．５のときの特性を表す曲線に挟まれた領域に多
く含まれている。

【００７１】すなわち、図７の呼気速度定数Ｋｅと動脈
血速度定数Ｋａとの関係に基づいて得られるＫｅ／Ｋａ
変換定数γは、全て、γ＝０．２４〜７．７の範囲に含
まれ、また、その多くは、γ＝０．３〜２．５の範囲に
含まれることとなる。このため、前記（４）式および
（５）式中のＫｅ／Ｋａ変換定数γの値の範囲は、γ＝
０．２４〜７．７の範囲、特に、γ＝０．３〜２．５の
範囲に限定することができる。

【００７２】このように、Ｋｅ／Ｋａ変換定数γの値が
とり得る範囲を特定することによって、後述するように
Ｋｅ／Ｋａ変換定数γの仮定値γαを変化させるべき範
囲を制限することができるため、処理時間の短縮化を図
ることができる。

【００７３】なお、このＫｅ／Ｋａ変換定数γの値の範
囲は、図７に示した実験データを例えば補完する他の実
験データに基づいて変更することが可能である。

【００７４】次に、Ｋｅ／Ｋａ変換定数γの値を特定す
るための手順について説明する。

【００７５】この一実施の形態では、以下に説明するＬ
ＧＣ法（Lambda-Guided Calculation Method）、すなわ
ち、前記（１）式中の脳／血液分配係数λを指標として
Ｋｅ／Ｋａ変換定数γを求める手法を用いて該Ｋｅ／Ｋ
ａ変換定数γを特定することを主な特徴としている。

【００７６】このＬＧＣ法では、まず、所望の範囲（例
えば、γ＝０．３〜２．５の範囲）に含まれるＫｅ／Ｋ
ａ変換定数γの任意の値を仮定値γαとして定めて、こ
の仮定値γαと測定値である呼気キセノンガス濃度Ｃｅ
（ｔ）とから前記（２）式〜（４）式に基づいて動脈キ
セノンガス濃度Ｃａ（ｔ）を求める。そして、この動脈
キセノンガス濃度Ｃａ（ｔ）と各測定時点Ｔにおける測
定値である脳キセノンガス濃度Ｃｂ（Ｔ）とから前記
（１）式の計算処理を行うことによって、前記仮定値γ
αに対応する脳／血液分配係数λの仮算出値λαを求め
る。

【００７７】同様に、仮定値γαを所望の範囲内で変化
させながらこの仮定値γαに対応する仮算出値λαをそ
れぞれ求める。そして、これら仮算出値λαを指標とし
ての後述する目標値λτと比較することによって、この
目標値λτに最も近い仮算出値λαを抽出し、この仮算
出値λαに対応する仮定値γα（すなわち、仮算出値λ
αが目標値λτに最も近づくときの仮定値γα）を真の
Ｋｅ／Ｋａ変換定数γとして確定する。

【００７８】このように、指標としての脳／血液分配係
数λの目標値λτに基づいてＫｅ／Ｋａ変換定数γを求
めることができるのは、脳５４の特に白質部分におい
て、脳／血液分配係数λと動脈血速度定数Ｋａｉ、Ｋａ
ｏとの間に次のような相関関係があるためである。

【００７９】図１０は、脳／血液分配係数λと動脈血速
度定数Ｋａｉ、Ｋａｏとの関係を表すグラフである。動
脈血速度定数Ｋａｉ、Ｋａｏが正しい値（ＲＥＡＬ値）
であるとき、白質部分における脳／血液分配係数λの値
は約１．５となることが知られている。一方、動脈血速
度定数Ｋａｉ、Ｋａｏが過大評価または過小評価された
とき、脳／血液分配係数λの値は１．５から離れる。

【００８０】このため、脳／血液分配係数λの仮算出値
λαが約１．５となるときの動脈血速度定数Ｋａｉ、Ｋ
ａｏの値を正しい値とみなし、さらに、このような正し
い値の動脈血速度定数Ｋａｉ、Ｋａｏを前記（４）式ま
たは（５）式から得ることが可能な仮定値γαを真のＫ
ｅ／Ｋａ変換定数γと判断することができる。

【００８１】なお、白質部分における脳／血液分配係数
λの正確な値は、後述するように、計算によって求める
ことができる。

【００８２】次に、コンピュータ４０における、図３の
前記ステップＳ１２の処理（脳血流量算出処理）を行う
ための手段（データ処理手段）について説明する。

【００８３】図１１は、コンピュータ４０のデータ処理
手段１００の概略的な構成を示す機能ブロック図であ
る。

【００８４】データ処理手段１００は、真のＫｅ／Ｋａ
変換定数γを求めるとともに、このＫｅ／Ｋａ変換定数
γに基づいて動脈血速度定数Ｋａｉ、Ｋａｏを求めるγ
設定手段（変換定数設定手段）１０２と、このγ設定手
段１０２で求められた動脈血速度定数Ｋａｉ、Ｋａｏに
基づいて脳血流量ｆおよび脳／血液分配係数λを求める
ｆ・λ算出手段１０４と、これらＫｅ／Ｋａ変換定数
γ、脳血流量ｆおよび脳／血液分配係数λ等に基づいて
出力用の画像データ等を形成する出力処理手段１０６と
を備えている。

【００８５】図１２は、γ設定手段１０２の概略的な構
成を示す機能ブロック図である。γ設定手段１０２は、
Ｋｅ算出手段１１０、Ｋｅ／Ｋａ変換手段（仮速度定数
設定手段）１１２、仮定値設定手段１１４およびＫａ記
憶手段１１６を備えている。

【００８６】Ｋｅ算出手段１１０は、まず、混合ガス供
給装置１６の濃度測定センサ４８（図１参照）で所定時
間（４０ｍｓ）毎に検出され、外部記憶装置５５に記憶
されている呼気キセノンガス濃度Ｃｅ（ｔ）を該外部記
憶装置５５から読み出す。そして、この呼気キセノンガ
ス濃度Ｃｅ（ｔ）から最小自乗法によって、前記（６）
式および（７）式として示した近似式（特に、これら式
中の呼気速度定数Ｋｅｉ、Ｋｅｏ）を求める。

【００８７】Ｋｅ／Ｋａ変換手段１１２は、Ｋｅ算出手
段１１０からの呼気速度定数Ｋｅｉ、Ｋｅｏと、仮定値
設定手段１１４からのＫｅ／Ｋａ変換定数γの仮定値γ
αとに基づいて前記（４）式および（５）式の計算処理
を行うことによって、動脈血速度定数Ｋａｉ、Ｋａｏの
仮算出値Ｋａｉα、Ｋａｏαを求める。

【００８８】この場合、仮定値設定手段１１４は、後述
するＫ−Ｓ演算手段１２４からの制御信号に従って、仮
定値γαを、例えば０．３から２．５まで、例えば０．
０１ずつ増加させる。

【００８９】Ｋａ記憶手段１１６には、仮定値設定手段
１１４からの仮定値γαと、これに対応するＫｅ／Ｋａ
変換手段１１２からの仮算出値Ｋａｉα、Ｋａｏαがそ
れぞれ蓄積される。

【００９０】また、γ設定手段１０２は、ＲＯＩデータ
抽出手段１２０、目標λ値（目標値）算出手段１２２、
Ｋ−Ｓ演算手段（仮分配係数算出手段）１２４およびγ
値確定手段１２６を備えている。

【００９１】図１３は、脳５４の断層画像β上に設定さ
れた関心領域ＲＯＩを示している。

【００９２】図１２のＲＯＩデータ抽出手段１２０は、
図１３に示すように、オペレータによって操作コンソー
ル５２を通じて設定された関心領域ＲＯＩ（図１３中、
○印で示した領域参照）に含まれる画素を特定する。

【００９３】この場合、オペレータは、表示装置５６の
画面１２８上に表示された脳５４の断層画像β（特に、
大脳基底核が出現する断層面の画像）から、特定の領域
（好ましくは、白質が多く存在する前頭葉白質部βｗが
含まれるような領域）を円で囲む操作によって、関心領
域ＲＯＩの設定作業を行う。

【００９４】このように、前頭葉白質部βｗを囲むよう
に関心領域ＲＯＩを設定することによって、関心領域Ｒ
ＯＩの設定作業にオペレータの主観が反映されてしまう
ことを防止できる。なお、画像処理によって、前頭葉白
質部βｗの位置を検出することにより、関心領域ＲＯＩ
を自動的に設定するようにしてもよい。

【００９５】また、ＲＯＩデータ抽出手段１２０は、図
１２に示すように、Ｘ線ＣＴ装置本体１８の検出器３０
（図２参照）で所定の時間間隔（約６０ｓ）の測定時点
Ｔ毎に検出され、外部記憶装置５５に記憶されている脳
キセノンガス濃度Ｃｂ（Ｔ）から、関心領域ＲＯＩに含
まれる画素に対応するものを抽出する。

【００９６】この場合、脳キセノンガス濃度Ｃｂ（Ｔ）
は、図２に示すように、Ｘ線ＣＴ装置本体１８からのＣ
Ｔ画像データに基づいて制御装置２０で予め計算され
て、外部記憶装置５５に記憶されているものとする。ま
た、脳キセノンガス濃度Ｃｂ（Ｔ）を求める際には、上
述したように、ノイズを除去するのに適切な、例えば、
９×９個の画素範囲における移動平均処理が施されてい
る。

【００９７】図１２に示すように、Ｋ−Ｓ演算手段１２
４は、ＲＯＩデータ抽出手段１２０からの脳キセノンガ
ス濃度Ｃｂ（Ｔ）（測定時点Ｔの情報も含んでいる。）
と、Ｋｅ／Ｋａ変換手段１１２からの動脈血速度定数Ｋ
ａｉ、Ｋａｏの仮算出値Ｋａｉα、Ｋａｏαとに基づい
て、前記（１）式〜（３）式を用いた計算処理を行う。
なお、前記（２）式および（３）式中のＫａｉ、Ｋａｏ
は、Ｋａｉα、Ｋａｏαに置き換えて計算を行う。

【００９８】この計算処理によって得られた脳／血液分
配係数λの仮算出値λαは、これに対応する仮定値設定
手段１１４からの仮定値γαとともに、γ値確定手段１
２６に供給される。なお、この計算処理では、実際に
は、脳血流量ｆも求められるが、この脳血流量ｆはγ値
確定手段１２６では考慮されない。

【００９９】また、Ｋ−Ｓ演算手段１２４は、１つの仮
定値γαに基づく仮算出値λαの算出処理が終わると、
次の仮定値γα（０．０１だけ値が増加した仮定値γ
α）を出力させるための制御信号を仮定値設定手段１１
４に対して出力する。

【０１００】γ値確定手段１２６には、Ｋ−Ｓ演算手段
１２４からの仮算出値λαとともに、目標λ値算出手段
１２２からの脳／血液分配係数λの目標値λτも供給さ
れている。

【０１０１】この目標値λτは、正常（健常）な脳５４
の白質部分の脳／血液分配係数λであり、操作コンソー
ル５２を通じてオペレータによって入力された、例えば
被検体１２のヘマトクリット値（赤血球容積率）Ｈ
［％］に基づいて求められる。この場合、ヘマトクリッ
ト値Ｈは、被検体１２の血液から得ることができる。

【０１０２】目標値λτは、実際には、図１４のλ−Ｈ
グラフ１３０をヘマトクリット値Ｈで検索することによ
って求められる。このλ−Ｈグラフ１３０は、ベール
（N.Veall）らによって明らかにされた、脳５４の白質
部分における脳／血液分配係数λとヘマトクリット値Ｈ
との関係を記録したものである（「The Partition of T
race Amounts of Xenon Between Human Blood and Brai
n Tissues at 37℃」、N.Veall, et al.、Phys. Med. B
iol.、1965、Vol.10、No.3、375-380参照）。

【０１０３】この場合、図１４のλ−Ｈグラフ１３０の
特性を表す計算式や、前記特性が記録されたテーブルを
用いて、ヘマトクリット値Ｈから目標値λτを得るよう
にしてもよい。

【０１０４】なお、脳５４の白質部分の脳／血液分配係
数λは、約１．５であることが知られている。従って、
この値（約１．５）を目標値λτとして設定するように
してもよい。

【０１０５】図１２に示すように、γ値確定手段１２６
は、Ｋ−Ｓ演算手段１２４からの各仮定値γαに対応す
る仮算出値λαを、目標λ値算出手段１２２からの目標
値λτと比較することによって、真のＫｅ／Ｋａ変換定
数γを特定する処理を行う。なお、この処理は、仮定値
設定手段１１４から最後の仮定値γα（例えば、２．
５）が出力されたことを伝える制御信号に基づいて開始
するようにしてもよい。

【０１０６】図１５は、γ値確定手段１２６における、
目標値λτを指標として真のＫｅ／Ｋａ変換定数γを特
定するための手段｛γ値（変換定数）抽出手段１４４、
フィルタリング手段１４６、平均値算出手段１４８｝を
示す機能ブロック図である。

【０１０７】図１６に示すように、γ値確定手段１２６
においては、Ｋ−Ｓ演算手段１２４からの各仮定値γα
に対応する仮算出値λαが、画素番号ｐ（ｐ＝１、２、
…）毎にテーブル１４０ｐとして分類されている。さら
に、これらテーブル１４０ｐ中においては、仮算出値λ
αが各仮定値γα毎に分類されている。

【０１０８】図１５および図１６に示すように、γ値抽
出手段１４４は、目標値λτを用いて、各テーブル１４
０ｐ中から目標値λτと最も値が近い仮算出値λαを検
索する。そして、検索された仮算出値λαとこれに対応
する仮定値γα（すなわち、仮算出値λαが目標値λτ
に最も近づくときの仮定値γα）を、それぞれ抽出仮算
出値λα′、抽出仮定値γα′として抽出する。

【０１０９】図１７に示すように、この抽出処理を経
て、抽出仮算出値λα′および抽出仮定値γα′が画素
番号ｐ毎に分類されたテーブル１４２が作られる。

【０１１０】図１５および図１７に示すように、フィル
タリング手段１４６は、テーブル１４２中の画素番号
ｐ、抽出仮算出値λα′および抽出仮定値γα′の組か
ら、抽出仮算出値λα′が所定の濾過範囲から外れてい
る組を除去する。この場合、この濾過範囲は、目標値λ
τに基づいて設定された範囲であり、例えば、目標値λ
τに対して該目標値λτの１０％の値を加減算して得ら
れた値（λτ±λτ×１０％）を上下限値とした範囲と
することが好ましい。

【０１１１】このようなフィルタリング処理を行うこと
によって、脳５４の白質部分に対応する抽出仮算出値λ
α′および抽出仮定値γα′の組が濾過値として抽出さ
れる。

【０１１２】続いて、平均値算出手段１４８は、フィル
タリング手段１４６におけるフィルタリング処理を経た
テーブル１４２′中の抽出仮定値γα′の平均値（すな
わち、関心領域ＲＯＩ中の白質部分に対応する抽出仮定
値γα′の平均値）を求める。そして、この平均値を、
真のＫｅ／Ｋａ変換定数γ（真値）として確定する。

【０１１３】なお、図１２のγ値確定手段１２６におい
ては、Ｋｅ／Ｋａ変換定数γの仮定値γαを振動させ
て、仮算出値λαが目標値λτに収束したとき（すなわ
ち、仮算出値λαが目標値λτに基づく収束条件を満足
したとき）に得られた仮定値γαを抽出仮定値γα′と
するようにしてもよい。

【０１１４】図１２に示すように、γ値確定手段１２６
で得られたＫｅ／Ｋａ変換定数γは、Ｋａ記憶手段１１
６および出力処理手段１０６にそれぞれ供給される。

【０１１５】Ｋａ記憶手段１１６は、仮定値設定手段１
１４からの各仮定値γαと、γ値確定手段１２６からの
Ｋｅ／Ｋａ変換定数γとを比較して、Ｋｅ／Ｋａ変換定
数γと一致する仮定値γαを抽出する。そして、この仮
定値γαに対応する（すなわち、この仮定値γαを用い
て求められた）Ｋｅ／Ｋａ変換手段１１２からの仮算出
値Ｋａｉα、Ｋａｏαを、真の動脈血速度定数Ｋａｉ、
Ｋａｏとしてｆ・λ算出手段１０４に供給する。

【０１１６】図１１に示すように、ｆ・λ算出手段１０
４は、γ設定手段１０２からの動脈血速度定数Ｋａｉ、
Ｋａｏと、外部記憶装置５５からの脳キセノンガス濃度
Ｃｂ（Ｔ）｛各測定時点ＴにＸ線ＣＴ装置１４で得られ
た脳キセノンガス濃度Ｃｂ（Ｔ）｝とに基づいて、各画
素毎に前記（１）式〜（３）式を用いた計算処理を行
う。そして、この計算処理によって得られた各画素の脳
血流量ｆおよび脳／血液分配係数λを出力処理手段１０
６に供給する。

【０１１７】出力処理手段１０６は、ｆ・λ算出手段１
０４からの脳血流量ｆおよび脳／血液分配係数λ、並び
にγ設定手段１０２からのＫｅ／Ｋａ変換定数γに基づ
いて、表示装置５６に表示させるための、または、プリ
ンタ５８からハードコピー５７として出力させるための
表示画像データ（後述するｆマップＭｆおよびλマップ
Ｍλ等の表示画像データ）を形成する。

【０１１８】図１８〜図２１は、表示装置５６の画面１
２８を示している。

【０１１９】図１８および図１９は、表示装置５６の画
面１２８上にカラーで表示された、脳５４のｆマップ
（脳血流量ｆの分布図）Ｍｆａ、Ｍｆｂを示している。
また、図２０および図２１は、表示装置５６の画面１２
８上にカラーで表示された、脳５４のλマップ（脳／血
液分配係数λの分布図）Ｍλａ、Ｍλｂを示している。

【０１２０】この場合、図１８および図２０は、それぞ
れ、呼気速度定数Ｋｅｉ、Ｋｅｏから動脈血速度定数Ｋ
ａｉ、Ｋａｏを得る際にＫｅ／Ｋａ変換定数γによる補
正が施された場合｛すなわち、前記（４）式および
（５）式に基づいて動脈血速度定数Ｋａｉ、Ｋａｏを得
た場合｝のｆマップＭｆａおよびλマップＭλａを示し
ており、一方、図１９および図２１は、それぞれ、Ｋｅ
／Ｋａ変換定数γによる補正が施されなかった場合｛す
なわち、前記（４）式および（５）式において、Ｋａｉ
＝Ｋｅｉ、かつ、Ｋａｏ＝Ｋｅｏとした場合｝のｆマッ
プＭｆｂおよびλマップＭλｂを示している。

【０１２１】なお、図１８〜図２１例のこれらｆマップ
Ｍｆａ、ＭｆｂおよびλマップＭλａ、Ｍλｂは、被検
体１２としての３２歳の健常な男性から得られたもので
ある。

【０１２２】図１８に示す、補正が施されたｆマップＭ
ｆａにおいては、図１９に示す、補正が施されなかった
ｆマップＭｆｂに比べて、脳血流量ｆが高い値となって
いる。また、図２０に示す、補正が施されたλマップＭ
λａにおいては、図２１に示す、補正が施されなかった
λマップＭλｂに比べて、脳／血液分配係数λが低い値
となっている。

【０１２３】すなわち、Ｋｅ／Ｋａ変換定数γによる補
正が施されなかった場合には、図１０に示すように、動
脈血速度定数Ｋａｉ、Ｋａｏが過大評価されるととも
に、脳／血液分配係数λが過大評価され、さらに、脳血
流量ｆが過小評価されている。

【０１２４】これに対して、Ｋｅ／Ｋａ変換定数γによ
る補正が施された場合には、動脈血速度定数Ｋａｉ、Ｋ
ａｏおよび脳／血液分配係数λの過大評価が解消され、
適正な値の脳血流量ｆが得られている。

【０１２５】図１８に示すように、補正が施されたｆマ
ップＭｆａから得られる、脳１００ｇ当たりの血流量
は、脳５４の左半球において約４９．７ｍｌ／１００ｇ
／ｍｉｎであり、右半球において約５３．２ｍｌ／１０
０ｇ／ｍｉｎである。これらは、成人の脳１００ｇ当た
りの平均的な血流量（約５４ｍｌ／１００ｇ／ｍｉｎ）
に非常に近い値であり、このことからも、図１８のｆマ
ップＭｆａに表されている脳血流量ｆの値が適正なもの
であることがわかる。

【０１２６】なお、図１８〜図２１に示したｆマップＭ
ｆａ、ＭｆｂおよびλマップＭλａ、Ｍλｂは、それぞ
れ、Ｒ（赤）、Ｇ（緑）、Ｂ（青）の３色でカラー表示
されているが、さらに多くの色（例えば、１３色）でカ
ラー表示を行うようにしてもよい。

【０１２７】このように、この発明の一実施の形態にお
いては、脳／血液分配係数λを指標としてＫｅ／Ｋａ変
換定数γを求め、さらに、呼気速度定数Ｋｅｉ、Ｋｅｏ
から動脈血速度定数Ｋａｉ、Ｋａｏを求める際に、この
Ｋｅ／Ｋａ変換定数γによる補正を施すようにしてい
る。このため、動脈キセノンガス濃度Ｃａ（ｔ）の代わ
りとして用いた呼気キセノンガス濃度Ｃｅ（ｔ）から、
脳血流量ｆを正確に求めることができる。

【０１２８】この場合、Ｋｅ／Ｋａ変換定数γは、目標
値λτに基づく抽出処理、およびフィルタリング処理を
経て得られている。このため、Ｋｅ／Ｋａ変換定数γの
正確な値を得ることができる。

【０１２９】また、Ｋｅ／Ｋａ変換定数γは、関心領域
ＲＯＩを対象とした平均化処理を経て得られている。こ
のため、Ｋｅ／Ｋａ変換定数γの一層正確な値を得るこ
とができる。

【０１３０】さらに、指標としての目標値λτは、被検
体１２の血液から得られるヘマトクリット値Ｈに基づい
て求められるため、得られたＫｅ／Ｋａ変換定数γの値
の正確性を一層確実なものとすることができる。

【０１３１】なお、上述の実施の形態においては、キセ
ノンガスの吸入過程の終了を呼気中のキセノンガス濃度
の飽和により判定し（図３のフローチャート中、ステッ
プＳ８参照）、空気の吸入に切り替えて行われるキセノ
ンガスの洗い出し過程の終了をキセノンガス濃度が所定
値以下になることにより判定しているが（図３のフロー
チャート中、ステップＳ１１参照）、吸入過程を所定時
間（例えば、約４ｍｉｎ）で終了し、洗い出し過程を所
定時間（例えば、約５ｍｉｎ）で終了するという時間管
理で行うようにしてもよい。

【０１３２】また、上述の実施の形態は、頭部の診断ば
かりではなく、頭部のように動脈血流が流入して静脈血
流として流出する通常の器官（臓器）、例えば、胃、
腸、膵臓、肝臓などの被検査部位にも適用することがで
きる。

【０１３３】

【発明の効果】以上説明したように、この発明によれ
ば、脳／血液分配係数λを指標として変換定数γを求め
るようにしているため、脳血流量を正確に求めることが
できる。

【図面の簡単な説明】

【図１】この発明の一実施の形態の全体構成を示す斜視
的模式図である。

【図２】この発明の一実施の形態の構成を示すブロック
図である。

【図３】この発明の一実施の形態の動作説明に供される
フローチャートである。

【図４】キセノンガスが供給されている被検体である患
者の脳をＸ線ＣＴ装置により撮影している状態を示す側
面視的模式図である。

【図５】脳キセノンガス濃度の変化を示す特性図であ
る。

【図６】動脈キセノンガス濃度の変化を示す特性図であ
る。

【図７】呼気速度定数と動脈血速度定数の関係を表すグ
ラフである。

【図８】呼気速度定数と動脈血速度定数との関係式の特
性を表すグラフである。

【図９】図７のグラフに図８のグラフが重ね合わされた
グラフである。

【図１０】脳／血液分配係数と動脈血速度定数の関係を
表すグラフである。

【図１１】コンピュータのデータ処理手段の概略的な構
成を示す機能ブロック図である。

【図１２】図１１のデータ処理手段のγ設定手段の概略
的な構成を示す機能ブロック図である。

【図１３】脳の断層画像上に設定された関心領域を示す
図である。

【図１４】λ−Ｈテーブルを示す図である。

【図１５】図１２のデータ処理手段を構成するγ値確定
手段における、目標値を指標として真のＫｅ／Ｋａ変換
定数を特定するための手段を示す機能ブロック図であ
る。

【図１６】各仮定値に対応する仮算出値が画素番号毎に
分類されたテーブルを示す図である。

【図１７】抽出仮算出値および抽出仮定値が画素番号毎
に分類されたテーブルを示す図である。

【図１８】表示装置の画面上に表示された、脳のｆマッ
プを示す図である。

【図１９】表示装置の画面上に表示された、脳のｆマッ
プを示す図である。

【図２０】表示装置の画面上に表示された、脳のλマッ
プを示す図である。

【図２１】表示装置の画面上に表示された、脳のλマッ
プを示す図である。

【符号の説明】

１０…脳血流量測定装置 １２…被検体 １４…Ｘ線ＣＴ装置 １６…混合ガス供
給装置 １８…Ｘ線ＣＴ装置本体 ２０…制御装置 ２８…Ｘ線管 ３０…検出器 ４８…濃度測定センサ ５０…コンピュー
タ ５２…操作コンソール ５４…脳 ５６…表示装置 １００…データ処
理手段 １０２…γ設定手段 １０４…ｆ・λ算
出手段 １１０…Ｋｅ算出手段 １１２…Ｋｅ／Ｋ
ａ変換手段 １１４…仮定値設定手段 １２０…ＲＯＩデ
ータ抽出手段 １２２…目標λ値算出手段 １２４…Ｋ−Ｓ演
算手段 １２６…γ値確定手段 １３０…λ−Ｈグ
ラフ １４４…γ値抽出手段 １４６…フィルタ
リング手段 １４８…平均値算出手段 Ｃａ（ｔ）…動脈
キセノンガス濃度 Ｃｂ（Ｔ）…脳キセノンガス濃度 Ｃｅ（ｔ）…呼気
キセノンガス濃度 Ｋａ…動脈血速度定数 Ｋｅ…終末呼気速
度定数 ｆ…脳血流量 λ…脳／血液分配
係数 λτ…目標値 γ…Ｋｅ／Ｋａ変
換定数 γα…仮定値 ＲＯＩ…関心領域 Ｈ…ヘマトクリット値 Ｍｆａ、Ｍｆｂ…
ｆマップ Ｍλａ、Ｍλｂ…λマップ

Claims (15)

【特許請求の範囲】
1. 【請求項１】キセノンＣＴ検査に係る終末呼気速度定数
   をＫｅとし、動脈血速度定数をＫａとしたときの関係式 Ｋａ＝γ×（１−ｅｘｐ（−Ｋｅ／γ）） 中の定数γを決定する方法において、 キセノンＣＴ画像上で関心領域を設定するＡステップ
   と、 設定した関心領域でキセノン分配係数λが所定の目標値
   に最も近づく前記定数γを求めるＢステップと、 を備えることを特徴とする定数決定方法。
2. 【請求項２】請求項１記載の定数決定方法において、 前記Ｂステップでは、前記定数γを０．２４から７．７
   の間の所望の範囲で変化させて前記分配係数λを計算
   し、前記分配係数λが前記目標値に最も近づく前記定数
   γを求めることを特徴とする定数決定方法。
3. 【請求項３】請求項２記載の定数決定方法において、 前記所望の範囲は、０．３から２．５までの範囲とする
   ことを特徴とする定数決定方法。
4. 【請求項４】請求項１〜３のいずれか１項に記載の定数
   決定方法において、 前記Ｂステップでは、前記関心領域に含まれる所定画素
   毎に、前記分配係数λ値が前記目標値に最も近づく前記
   定数γを求め、得られた定数γ値を平均して目的の定数
   γ値とすることを特徴とする定数決定方法。
5. 【請求項５】請求項１〜４のいずれか１項に記載の定数
   決定方法において、 前記Ａステップでは、前記関心領域を脳の白質を含む領
   域に設定し、 前記Ｂステップでは、前記目標値をヘマトクリット値に
   応じて決定することを特徴とする定数決定方法。
6. 【請求項６】被検体にキセノンガスを供給するガス供給
   装置と、 前記被検体の終末呼気中のキセノンガスの濃度（以下、
   呼気キセノンガス濃度と記す。）を測定するための濃度
   測定装置と、 前記被検体の被検査部位のキセノンガスの濃度（以下、
   被検査部位キセノンガス濃度と記す。）を得るために前
   記被検査部位のＣＴ画像データを取得するＸ線ＣＴ装置
   本体と、 前記ＣＴ画像データに基づいて前記被検査部位キセノン
   ガス濃度を求めるとともに、該被検査部位キセノンガス
   濃度および前記呼気キセノンガス濃度に基づいて前記被
   検査部位の血流量を求めるデータ処理装置と、 を備えたキセノンＣＴ装置において、 前記データ処理装置は、前記呼気キセノンガス濃度の速
   度定数（以下、呼気速度定数と記す。）を動脈中の血流
   のキセノンガスの濃度（以下、動脈キセノンガス濃度と
   記す。）の速度定数（以下、動脈血速度定数と記す。）
   に変換するための変換定数γに基づいて、前記被検査部
   位と前記被検体の血液との間のキセノンガスの分配係数
   λを求めるとともに、該分配係数λが所定の目標値に最
   も近づくときの前記変換定数γを真値として確定するこ
   とを特徴とするキセノンＣＴ装置。
7. 【請求項７】請求項６記載のキセノンＣＴ装置におい
   て、 前記データ処理装置は、前記変換定数γを求める変換定
   数設定手段を備え、 前記変換定数設定手段は、前記変換定数γの仮定値を設
   定するとともに、該仮定値を変化させる仮定値設定手段
   を有することを特徴とするキセノンＣＴ装置。
8. 【請求項８】請求項７記載のキセノンＣＴ装置におい
   て、 前記変換定数設定手段は、 前記変換定数γの仮定値に基づいて、前記呼気速度定数
   から前記動脈血速度定数の仮算出値を求める仮速度定数
   設定手段と、 前記動脈血速度定数の仮算出値と前記被検査部位キセノ
   ンガス濃度とから、前記分配係数λの仮算出値を求める
   仮分配係数算出手段と、 前記変換定数γの各仮定値に基づいて得られた前記分配
   係数λの各仮算出値のうち、前記目標値に最も近づくも
   のに対応する前記変換定数γの仮定値を、前記真値とし
   て確定すべき抽出仮定値として抽出する変換定数抽出手
   段と、 を有することを特徴とするキセノンＣＴ装置。
9. 【請求項９】請求項８記載のキセノンＣＴ装置におい
   て、 前記変換定数設定手段は、前記抽出仮定値のうち、該抽
   出仮定値に対応する前記分配係数λの仮算出値が所定の
   濾過範囲内に含まれるものを、前記真値として確定すべ
   き濾過値として求めるフィルタリング手段を有すること
   を特徴とするキセノンＣＴ装置。
10. 【請求項１０】請求項９記載のキセノンＣＴ装置におい
    て、 前記変換定数設定手段は、 前記ＣＴ画像データに含まれる各画素のデータのうち、
    所定の複数の画素に対応するデータに基づいて得られた
    前記被検査部位キセノンガス濃度をそれぞれ抽出するＲ
    ＯＩデータ抽出手段と、 抽出された前記被検査部位キセノンガス濃度に基づいて
    それぞれ得られた前記濾過値の平均値を前記変換定数γ
    の真値として求める平均値算出手段と、 を有することを特徴とするキセノンＣＴ装置。
11. 【請求項１１】請求項７〜１０のいずれか１項に記載の
    キセノンＣＴ装置において、 前記仮定値設定手段は、前記変換定数γの仮定値を０．
    ２４から７．７の間の所望の範囲で変化させることを特
    徴とするキセノンＣＴ装置。
12. 【請求項１２】請求項１１記載のキセノンＣＴ装置にお
    いて、 前記所望の範囲は、０．３から２．５までの範囲とする
    ことを特徴とするキセノンＣＴ装置。
13. 【請求項１３】請求項７〜１２のいずれか１項に記載の
    キセノンＣＴ装置において、 前記変換定数設定手段は、前記目標値を決定する目標値
    算出手段を有し、 前記目標値算出手段は、前記被検査部位を前記被検体の
    脳とし、かつ、前記変換定数γを求めるための関心領域
    を前記脳の白質を含む領域に設定したとき、前記目標値
    をヘマトクリット値に応じて決定することを特徴とする
    キセノンＣＴ装置。
14. 【請求項１４】請求項６〜１３のいずれか１項に記載の
    キセノンＣＴ装置において、 前記変換定数γに基づいて前記呼気速度定数から前記動
    脈血速度定数を求めるための関係式は、前記呼気速度定
    数をＫｅとし、前記動脈血速度定数をＫａとしたとき、 Ｋａ＝γ×（１−ｅｘｐ（−Ｋｅ／γ）） で表されることを特徴とするキセノンＣＴ装置。
15. 【請求項１５】請求項６〜１４のいずれか１項に記載の
    キセノンＣＴ装置において、 さらに、前記血流量および／または前記分配係数λの分
    布図を表示する表示装置を有することを特徴とするキセ
    ノンＣＴ装置。