(57)【要約】【課題】非採血により、被験者の所定部位の血流量を推定する。【解決手段】被験者にＲＩ（Ｒａｄｉｏ  Ｉｓｏｔｏｐｅ）を含む薬剤を投与してから所定の計測時間の間計測した、肺ダイナミック画像１１１を記憶する記憶部１１と、肺ダイナミック画像１１１に基づいて、計測時間経過後の測定基準時刻における薬剤血中濃度Ｃａを推定するＣａ推定部２０と、測定時刻よりも後の撮影時間帯に、被験者の所定部位を撮影したＳＰＥＣＴ画像１１７を記憶する記憶部１１と、撮影時間帯に撮影されたＳＰＥＣＴ画像１１７をに変換するＳＰＥＣＴ変換部３３と、Ｃａ推定部２０により推定された薬剤血中濃度Ｃａと、測定基準時刻相当のＳＰＥＣＴ画像とに基づいて、所定部位の血流量ＣＢＦを推定するＣＢＦ算出部３５と、を備える。【選択図】  図１
2008185499.tif 000002
【特許請求の範囲】
【請求項１】
被験者にＲＩ（Ｒａｄｉｏ  Ｉｓｏｔｏｐｅ）を含む薬剤を投与してから所定の計測時間の間計測した、前記被験者からのＲＩ情報を記憶する記憶手段と、
前記ＲＩ情報に基づいて、前記薬剤投与から前記計測時間経過後の測定基準時刻までの薬剤血中濃度を推定する薬剤血中濃度推定手段と、
前記測定時刻よりも後の撮影時間帯に、前記被験者の所定部位を撮影した３次元画像を記憶する記憶手段と、
前記撮影時間帯に撮影された前記３次元画像を、前記測定基準時刻相当の３次元画像に変換する変換手段と、
前記薬剤血中濃度推定手段により推定された前記薬剤血中濃度と、前記変換手段により変換された前記測定基準時刻相当の３次元画像とに基づいて、前記所定部位の血流量を推定する血流量推定手段と、を備えた血流量の推定システム。
【請求項２】
前記薬剤血中濃度推定手段は、前記ＲＩ情報に基づいて前記被験者の薬剤体内循環量を推定し、予め定められている関数を用いて、前記推定された薬剤体内循環量を前記薬剤血中濃度へ変換することを特徴とする請求項１記載の血流量の推定システム。
【請求項３】
前記薬剤血中濃度推定手段は、前記薬剤を投与してから前記計測基準時刻までの間に、前記被験者から動脈持続採血法により採血した所定量の血液から観測される放射線量を推定することを特徴とする請求項１記載の血流量の推定システム。
【請求項４】
前記ＲＩ情報は、前記被験者の肺から放出される放射線量の時間変化を示す情報であり、
前記薬剤血中濃度推定手段は、前記肺から放出される放射線量の時間変化情報に基づいて、前記薬剤を投与してから前記測定基準時刻までの間に肺から流出した薬剤量を推定し、前記肺から流出した薬剤量を予め定められている関数を用いて前記薬剤血中濃度へ変換することを特徴とする請求項１記載の血流量の推定システム。
【請求項５】
前記ＲＩ情報は、異なる時刻に撮影された前記被験者の肺の複数のプラナ画像であり、
前記薬剤血中濃度推定手段は、各プラナ画像における肺領域の画素の合計画素値の時間変化を算出し、当該合計画素値の時間変化に基づいて肺から流出した薬剤量を推定し、前記肺から流出した薬剤量を、実測値に基づいて予め定められている関数を用いて前記薬剤血中濃度へ変換することを特徴とする請求項１記載の血流量の推定システム。
【請求項６】
前記所定部位を撮影した３次元画像は、脳のＳＰＥＣＴ（Ｓｉｎｇｌｅ  Ｐｈｏｔｏｎ  Ｅｍｉｓｓｉｏｎ  Ｃｏｍｐｕｔｅｄ  Ｔｏｍｏｇｒａｐｈｙ）画像であり、
前記血流量推定手段は、前記被験者の脳血流量を推定することを特徴とする請求項１記載の血流量の推定システム。
【請求項７】
被験者にＲＩ（Ｒａｄｉｏ  Ｉｓｏｔｏｐｅ）を含む薬剤を投与してから所定の計測時間の間、前記被験者からのＲＩ情報を計測するステップと、
前記計測されたＲＩ情報に基づいて、前記薬剤投与から前記計測時間経過後の測定基準時刻までの薬剤血中濃度を推定するステップと、
前記測定時刻よりも後の撮影時間帯に、前記被験者の所定部位の３次元画像を撮影するステップと、
前記撮影時間帯に撮影された前記３次元画像を、前記測定基準時刻相当の３次元画像に変換するステップと、
前記推定された薬剤血中濃度と、前記変換された前記測定基準時刻相当の３次元画像とに基づいて、前記所定部位の血流量を推定するステップと、を備えた血流量の推定方法。
【請求項８】
所定部位の血流量を推定するためのコンピュータプログラムであって、
被験者にＲＩ（Ｒａｄｉｏ  Ｉｓｏｔｏｐｅ）を含む薬剤を投与してから所定の計測時間の間計測して得たＲＩ情報に基づいて、前記薬剤投与から前記計測時間経過後の測定基準時刻までの薬剤血中濃度を推定するステップと、
前記測定時刻よりも後の撮影時間帯に撮影した、前記被験者の所定部位の３次元画像を、前記測定基準時刻相当の３次元画像に変換するステップと、
前記推定された薬剤血中濃度と、前記変換された前記測定基準時刻相当の３次元画像とに基づいて、前記所定部位の血流量を推定するステップと、をコンピュータに実行させる、薬剤血中濃度推定のためのコンピュータプログラム。
【請求項９】
被験者にＲＩ（Ｒａｄｉｏ  Ｉｓｏｔｏｐｅ）を含む薬剤を投与してから所定の計測時間の間計測して得たＲＩ情報を記憶した記憶手段と、
前記記憶手段に記憶されたＲＩ情報に基づいて、前記薬剤投与から前記計測時間経過後の測定基準時刻までの薬剤血中濃度を推定する薬剤血中濃度推定手段と、を備えた薬剤血流濃度の推定装置。
【請求項１０】
薬剤血中濃度推定のためのコンピュータプログラムであって、
被験者にＲＩ（Ｒａｄｉｏ  Ｉｓｏｔｏｐｅ）を含む薬剤を投与してから所定の計測時間の間計測して得たＲＩ情報を記憶するステップと、
前記記憶されたＲＩ情報に基づいて、前記薬剤投与から前記計測時間経過後の測定基準時刻までの薬剤血中濃度を推定するステップと、をコンピュータに実行させる、薬剤血中濃度推定のためのコンピュータプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ＲＩ（Ｒａｄｉｏ  Ｉｓｏｔｏｐｅ）を含む薬剤を投与した被験者からの放射能分布を計測して、被験者からの採血を行わずに所定部位の血流量を推定する技術に関する。
【背景技術】
【０００２】
ヒトの組織血流を推定するための手法として、マイクロスフェア法が知られている。マイクロスフェア法では、ＲＩ（Ｒａｄｉｏ  Ｉｓｏｔｏｐｅ）を含む薬剤（例えば、臨床用脳血流製剤）を被験者に投与し、この被験者からの放射能の分布及び測定基準時刻までの薬物動脈血中量を測定することによって、所定部位（例えば、脳）の血流量を推定する（例えば、非特許文献１）。
【非特許文献１】「ＳＰＥＣＴ機能画像  −定量化の基礎と臨床−」８６−９３ページ、１９９８年４月６日第１刷発行、編集者：西村恒彦、発行所：株式会社メジカルビュー社、発売所：株式会社グロービュー社
【発明の開示】
【発明が解決しようとする課題】
【０００３】
マイクロスフェア法によって血流量を推定するには、測定基準時刻までの薬剤血中濃度Ｃａを知る必要があり、従来は動脈持続採血法によって薬剤投与後の被験者から動脈血を採取していた。
【０００４】
しかしながら、この動脈採血は、被験者にとって非常に大きな負担となっているので、非採血による血流量測定手法の開発が望まれていた。
【０００５】
そこで、本発明の目的は、非採血により、被験者の所定部位の血流量を推定するための技術を提供することである。
【０００６】
本発明の別の目的は、マイクロスフェア法によって血流量を推定するために、採血を行わずに薬剤血中濃度Ｃａを推定することである。
【課題を解決するための手段】
【０００７】
本発明の一実施態様に従う血流量の推定システムは、被験者にＲＩ（Ｒａｄｉｏ  Ｉｓｏｔｏｐｅ）を含む薬剤を投与してから所定の計測時間の間計測した、前記被験者からのＲＩ情報を記憶する記憶手段と、前記ＲＩ情報に基づいて、前記薬剤投与から前記計測時間経過後の測定基準時刻までの薬剤血中濃度を推定する薬剤血中濃度推定手段と、前記測定時刻よりも後の撮影時間帯に、前記被験者の所定部位を撮影した３次元画像を記憶する記憶手段と、前記撮影時間帯に撮影された前記３次元画像を、前記測定基準時刻相当の３次元画像に変換する変換手段と、前記薬剤血中濃度推定手段により推定された前記薬剤血中濃度と、前記変換手段により変換された前記測定基準時刻相当の３次元画像とに基づいて、前記所定部位の血流量を推定する血流量推定手段と、を備える。
【０００８】
好適な実施形態では、前記薬剤血中濃度推定手段は、前記ＲＩ情報に基づいて前記被験者の薬剤体内循環量を推定し、予め定められている関数を用いて、前記推定された薬剤体内循環量を前記薬剤血中濃度へ変換するようにしてもよい。
【０００９】
好適な実施形態では、前記薬剤血中濃度推定手段は、前記薬剤を投与してから前記計測基準時刻までの間に、前記被験者から持続採血法により採血した所定量の血液から観測される放射線量を推定するようにしてもよい。
【００１０】
好適な実施形態では、前記ＲＩ情報は、前記被験者の肺から放出される放射線量の時間変化を示す情報であり、前記薬剤血中濃度推定手段は、前記肺から放出される放射線量の時間変化情報に基づいて、前記薬剤を投与してから前記測定基準時刻までの間に肺から流出した薬剤量を推定し、前記肺から流出した薬剤量を予め定められている関数を用いて前記薬剤血中濃度へ変換することを特徴とするようにしてもよい。
【００１１】
好適な実施形態では、前記ＲＩ情報は、異なる時刻に撮影された前記被験者の肺の複数のプラナ画像であり、前記薬剤血中濃度推定手段は、各プラナ画像における肺領域の画素の合計画素値の時間変化を算出し、当該合計画素値の時間変化に基づいて肺から流出した薬剤量を推定し、前記肺から流出した薬剤量を、実測値に基づいて予め定められている関数を用いて前記薬剤血中濃度へ変換するようにしてもよい。
【００１２】
好適な実施形態では、前記所定部位を撮影した３次元画像は、脳のＳＰＥＣＴ（Ｓｉｎｇｌｅ  Ｐｈｏｔｏｎ  Ｅｍｉｓｓｉｏｎ  Ｃｏｍｐｕｔｅｄ  Ｔｏｍｏｇｒａｐｈｙ）画像であり、前記血流量推定手段は、前記被験者の脳血流量を推定するようにしてもよい。
【００１３】
本発明の一実施態様に従う薬剤血流濃度の推定装置は、被験者にＲＩ（Ｒａｄｉｏ  Ｉｓｏｔｏｐｅ）を含む薬剤を投与してから所定の計測時間の間計測して得たＲＩ情報を記憶した記憶手段と、前記記憶手段に記憶されたＲＩ情報に基づいて、前記薬剤投与から前記計測時間経過後の測定基準時刻までの薬剤血中濃度を推定する薬剤血中濃度推定手段と、を備える。
【発明を実施するための最良の形態】
【００１４】
以下、本発明の一実施形態に係る血流量推定システムについて、図面を参照して説明する。
【００１５】
図１は、本実施形態に係る血流量推定システムの構成図である。
【００１６】
血流量推定システムは、血流量推定装置１と、撮影装置２とを含む。そして、血流量推定装置１が、撮影装置２で撮影された画像データを処理して、非採血で脳などの所定部位の血流量を推定する。
【００１７】
撮影装置２は、例えば、ガンマカメラ型のＳＰＥＣＴ画像撮影装置である。撮影装置２では、ＲＩを含む薬剤が投与された被験者の２次元画像（プラナ画像）、及び３次元画像（ＳＰＥＣＴ画像）を撮影することができる。
【００１８】
血流量推定装置１は、例えば汎用的なコンピュータシステムにより構成され、以下に説明する血流量推定装置１内の個々の構成要素または機能は、例えば、コンピュータプログラムを実行することにより実現される。
【００１９】
血流量推定装置１は、撮影装置２で撮影した画像データを記憶する画像データ記憶部１１と、所定の関数を記憶する関数記憶部１３と、薬剤血中濃度Ｃａを推定するＣａ推定部２０と、マイクロスフェア法に従う処理を行うマイクロスフェア処理部３０とを備える。
【００２０】
画像データ記憶部１１には、被験者にＲＩを投与した後、以下に説明する撮影プロトコールに従って撮影することによって得た画像データが記憶される。本実施形態では、例えば、肺ダイナミック画像１１１と、脳スタティック画像１１３、１１５と、ＳＰＥＣＴ画像１１７とが画像データ記憶部１１に記憶される。
【００２１】
図２は、撮影プロトコール及び対応する画像データを示す。
【００２２】
撮影プロトコールとは、薬剤を投与した時刻を基準とし、この時刻からの経過時間ごとの撮影箇所及び撮影方法を規定するものである。
【００２３】
本実施形態の撮影プロトコールは、まず、ＲＩを含む薬剤投与から所定の計測時間（例えば３分間）の間、被験者の体内における放射能の分布状態を把握するためのＲＩ情報を取得する。すなわち、このプロセスが、例えば図２の（１）及び（２）に相当する。
（１）薬剤投与からの１分間（０分００秒から１分００秒）は、撮影装置２が肺のプラナ画像をダイナミック収集する。このときの各画像の収集時間は１秒である。これにより、６０枚の肺ダイナミック画像１１１が得られる。
（２）薬剤投与から１分００秒から３分００秒までの間も、肺のプラナ画像をダイナミック収集する。このときの各画像の収集時間は１０秒である。これにより、１２枚の肺ダイナミック画像１１１が得られる。
【００２４】
次に、所定の測定基準時刻（例えば、薬剤投与から５分後）における、測定対象の所定部位の２次元画像を取得する。すなわち、このプロセスが、例えば図２の（３）に相当する。
（３）薬剤投与から４分３０秒から５分３０秒までの間は、脳のプラナ画像をスタティック収集する。これにより、脳スタティック画像１１３が得られる。この脳スタティック画像１１３は、薬剤投与から５分後のプラナ画像であるから、以降、ＰＬ５と呼ぶことがある。
【００２５】
次に、測定対象部位の３次元画像を撮影する。ここで、ＳＰＥＣＴ収集をして３次元画像を得るためには２０分程度の時間を要すること、及び、体内の薬剤の分布が安定している時間帯に撮影することが好ましいので、この撮影は測定基準時刻よりも後に撮影（図２の（５）に相当）が行われる。なお、このＳＰＥＣＴ画像の撮影の前後には、このＳＰＥＣＴ画像を修正するための２次元画像の撮影（図２の（４）（６）に相当）が行われる。
（４）薬剤投与から１２分３０秒から１３分３０秒までの間は、脳のプラナ画像をスタティック収集する。これにより、脳スタティック画像１１５ａが得られる。
（５）薬剤投与から１５分００秒から３５分００秒までの間は、脳を３次元的にＳＰＥＣＴ収集する。撮影装置２は、このＳＰＥＣＴ収集したデータに基づいてＳＰＥＣＴ画像１１７を合成する。ＳＰＥＣＴ画像１１７は３次元画像であるから、ＳＰＥＣＴ画像のデータはボクセル値により構成されている。
なお、このＳＰＥＣＴ画像１１７は、薬剤投与後２０分を中心とする時間に撮影されたものであるから、ＳＰＥＣＴ２０と呼ぶことがある。
（６）薬剤投与から３６分３０秒から３７分３０秒までの間は、脳のプラナ画像をスタティック収集する。これにより、脳スタティック画像（ＰＬ２０ｂ）１１５ｂが得られる。
【００２６】
ここで、脳スタティック画像１１５ａ及び１１５ｂは、いずれも、薬剤投与後２０分に相当するプラナ画像を取得するために用いられる。そこで、以降、脳スタティック画像１１５ａをＰＬ２０ａ、脳スタティック画像１１５ｂをＰＬ２０ｂと呼ぶことがある。
【００２７】
再び図１を参照すると、Ｃａ推定部２０は、肺ダイナミック画像１１１に基づいて、薬剤投与から所定時間が経過した測定基準時刻における被験者の薬剤血中濃度Ｃａを推定する。すなわち、本実施形態では、従来は被験者から採血した血液を用いて実測することによって得られた薬剤血中濃度Ｃａを、採血することなく画像データから推定する。なお、動脈持続採血法により実測される薬剤血中濃度Ｃａは、薬剤投与からの経過時間ｔの関数であるＣａ（ｔ）を用いると以下の式で表される。
【００２８】
【数１】
2008185499.tif 000003
【００２９】
Ｃａ推定部２０は、図１に示すように、肺ダイナミック画像１１１から肺領域のＲＯＩを定めるＲＯＩ決定部２１と、肺から流出し、体内を循環している薬剤量（薬剤体内循環量）ΔＬを算出するΔＬ算出部２３と、薬剤体内循環量ΔＬを薬剤血中濃度Ｃａに変換する変換処理部２５とを備える。
【００３０】
ＲＯＩ決定部２１は、例えば以下に示すような手法で、肺ダイナミック画像１１１における肺のＲＯＩを決定する。
【００３１】
図３は、表示装置４に表示される肺のＲＯＩ入力画面の例を示す。
【００３２】
例えば、ＲＯＩ決定部２１は、同図に示すように、肺ダイナミック画像１１１のうちの一枚を表示装置４に表示させる。そして、ユーザが入力装置３を用いて肺領域を囲むと、ＲＯＩ決定部２１は、この囲まれた領域を肺のＲＯＩ５１として決定する。
【００３３】
ΔＬ算出部２３は、全肺ダイナミック画像１１１における肺のＲＯＩ５１内の全画素の画素値（カウント値）に基づいて、薬剤体内循環量ΔＬを算出する。例えば、ΔＬ算出部２３は、各肺ダイナミック画像１１１のそれぞれについて、肺ＲＯＩ５１内の全画素の画素値の総和を算出する。これを、横軸を薬剤投与からの経過時間、縦軸を各ダイナミック画像１１１における肺ＲＯＩ５１内画素の画素値の総和としてグラフ化したタイムアクティビティカーブを生成する。
【００３４】
図４にタイムアクティビティカーブの例を示す。本実施形態では、薬剤投与から３分間の肺ダイナミック画像１１１を取得したので、薬剤投与から３分間のタイムアクティビティカーブが得られる。
【００３５】
ΔＬ算出部２３は、図４に示すような、薬剤投与から３分間のタイムアクティビティカーブ（実線）に基づいて、薬剤投与から５分、つまり測定基準時刻までの間に肺から流出した薬剤量を推定する。例えば、ΔＬ算出部２３は、肺ＲＯＩ内画素値の合計が最大（Ｌｐｅａｋ）となる時刻（Ｔｐｅａｋ）以降のタイムアクティビティカーブを指数関数で近似し、薬剤投与から３分以降の画素値の総和を予測する。そして、Ｌｐｅａｋと、この近似した指数関数に基づいて算出した５分後の画素値の総和（Ｌ５）との差を求めることにより、薬剤体内循環量ΔＬを算出する。
【００３６】
次に、改めて図１を参照すると、変換処理部２５は、所定の関数を用いて、上のようにして求めた薬剤体内循環量ΔＬを薬剤血中濃度Ｃａに変換する。例えば、本実施形態では、変換処理２５は、関数記憶部１３に記憶されている関数を用いて、薬剤体内循環量ΔＬを薬剤血中濃度Ｃａに変換する。ここで、まず関数記憶部１３に記憶されている関数について説明する。
【００３７】
関数記憶部１３には、事前に予め計測により求められた薬剤血中濃度Ｃａと、薬剤体内循環量ΔＬとの関係を示す関数が記憶されている。例えば、本実施形態では、この関数は上記のような薬剤体内循環量ΔＬと、動脈持続採血法により被験者から採血した血液の放射線量を実測して得た薬剤血中濃度Ｃａとの関係を示す関数である。
【００３８】
図５は、この関数の説明図である。すなわち、ある被験者について、肺ダイナミック画像１１１の撮影と同時に、動脈持続採血法による採血を５分間行う。そして、ダイナミック画像１１１からは上述の要領で薬剤体内循環量ΔＬを算出する。一方、採血した血液については、放射線量をカウントして薬剤血中濃度Ｃａを実測する。そして、最小自乗法などの所定の手法を用いて、薬剤体内循環量ΔＬ及び薬剤血中濃度Ｃａの関係を示す関数を導く。本実施形態では、例えば、ΔＬ及びＣａのそれぞれの対数が１次関数で近似でき、
ｌｎ（Ｃａ）＝Ａ×ｌｎ（ΔＬ）＋Ｂ
となる。
【００３９】
なお、薬剤血中濃度Ｃａは、動脈持続採血法により採血した全血を用いて測定した値であってもよいし、採血した血液をオクタノール抽出してから測定した値であってもよい。
【００４０】
再び図１を参照すると、変換処理部２５は、ΔＬ算出部２３が算出した薬剤体内循環量ΔＬを関数記憶部１３に記憶されている関数を適用して、薬剤血中濃度Ｃａに変換する。
【００４１】
これにより、事前に関数を決定しておけば、被験者からの採血を行うことなく画像データに基づいて薬剤血中濃度Ｃａを推定することができる。
【００４２】
マイクロスフェア処理部３０は、ＳＰＥＣＴ画像１１７（ＳＰＥＣＴ２０）を用いて、測定基準時刻である薬剤投与から５分後のＳＰＥＣＴ画像に変換を行って、薬剤投与から５分後の脳内の血流量を推定する。ここで、薬剤投与から５分後のＳＰＥＣＴ画像をＳＰＥＣＴ５と呼ぶ。
【００４３】
マイクロスフェア処理部３０は、脳スタティック画像１１３，１１５に基づいて脳のＲＯＩを定めるＲＯＩ決定部３１と、ＳＰＥＣＴ２０をＳＰＥＣＴ５に変換するＳＰＥＣＴ変換部３３と、測定基準時刻の脳の血流量ＣＢＦ（Ｃｅｒｅｂｒａｌ  Ｂｌｏｏｄ  Ｆｌｏｗ）を算出するＣＢＦ算出部３５とを備える。
【００４４】
ＲＯＩ決定部３１は、脳スタティック画像１１３，１１５における脳のＲＯＩを決定する。
【００４５】
図６は、表示装置４に表示される脳のＲＯＩ入力画面の例を示す。
【００４６】
例えば、ＲＯＩ決定部３１は、同図に示すように、脳スタティック画像１１３，１１５のいずれかを表示装置４に表示させる。そして、ユーザが入力装置３を用いて脳領域を囲むと、ＲＯＩ決定部３１は、この囲まれた領域を脳のＲＯＩ５３として決定する。各スタティック画像ごとにＲＯＩを指定してもよいし、同一被験者の画像については、一つのスタティック画像で囲ったＲＯＩを他のスタティック画像に適用してもよい。
【００４７】
ＳＰＥＣＴ変換部３３は、式（１）により、薬剤投与後２０分を中心とする時間帯に撮影されたＳＰＥＣＴ２０を、薬剤投与後５分の状態に相当するＳＰＥＣＴ５に変換する。例えば、ＳＰＥＣＴ変換部３３は、ＳＰＥＣＴ２０の各ボクセル値に対して、薬剤投与後２０分のプラナ画像ＰＬ２０における脳ＲＯＩ５３内の全画素の平均値と、薬剤投与後５分のプラナ画像ＰＬ５における脳ＲＯＩ５３内の全画素の平均値との比を用いて、ＳＰＥＣＴ２０の各ボクセル値を、ＳＰＥＣＴ５の各ボクセル値に変換する。
【００４８】
【数２】
2008185499.tif 000004
【００４９】
ここで、薬剤投与後２０分のプラナ画像を直接撮影することができないので、ＰＬ２０ａとＰＬ２０ｂとの平均を、薬剤投与後２０分のプラナ画像として扱っている。
【００５０】
ＣＢＦ算出部３５は、Ｃａ推定部２０で推定された薬剤投与から５分間の薬剤血中濃度Ｃａ及びＳＰＥＣＴ変換部３３で算出されたＳＰＥＣＴ５のボクセル値を用いて、式（２）により、ボクセルごとにＣＢＦを算出する。
【００５１】
【数３】
2008185499.tif 000005
ＣＣＦ（Ｃｒｏｓｓ  Ｃａｌｉｂｒａｔｉｏｎ  Ｆａｃｔｏｒ）：ボクセル値を容積に変換する換算係数
【００５２】
なお、ＣＣＦは、被験者ごとの変換処理を行うたびに入力を受け付けるようにしてもよい。
【００５３】
これにより、測定基準時刻である薬剤投与から５分後の脳の血流量を、被験者から採血することなく推定することができる。
【００５４】
次に、上記構成を備えるシステムを用いて、被験者の脳血流を推定するための処理手順について、図７に示すフローチャートに従って説明する。
【００５５】
まず、事前処理として、関数記憶部１３に記憶されるべき関数を生成する（Ｓ１１）。
【００５６】
次に、被験者に薬剤を投与し、図２に示す撮影プロトコールに従って、撮影装置２を用いて撮影を行う（Ｓ１２）。撮影後、肺ダイナミック画像１１１，脳スタティック画像１１３，１１５、及び所定の合成処理を行って合成されたＳＥＰＣＴ画像１１７が画像データ記憶部１１に格納される（Ｓ１３）。
【００５７】
これ以降、血流量推定装置１が、画像データ記憶部１１に格納された画像データに基づいて脳血流量を推定する。すなわち、ＲＯＩ決定部２１が、肺ダイナミック画像１１１における肺のＲＯＩを決定する（Ｓ１４）。そして、ΔＬ算出部２３が、肺ダイナミック画像１１１に基づいて、図４に示すタイムアクティビティカーブの近似式を求め、これから薬剤体内循環量ΔＬを算出する（Ｓ１５）。変換処理部２５が、関数記憶部１３の関数を用いて、薬剤体内循環量ΔＬを薬剤血中濃度Ｃａに変換する（Ｓ１６）。
【００５８】
次に、ＲＯＩ決定部３１が、脳スタティック画像１１３，１１５おける脳のＲＯＩを決定する（Ｓ１７）。ここで決定された脳のＲＯＩを用いて、ＳＰＥＣＴ変換部３３は、ＳＰＥＣＴ２０をＳＰＥＣＴ５に変換する（Ｓ１８）。そして、ＣＢＦ算出部３５は、ＳＰＥＣＴ５を用いて、ボクセル単位に脳の血流量ＣＢＦを算出する（Ｓ１９）。
【００５９】
以上の処理により、被験者ごとの採血を行わずに、撮影装置２で撮影した画像データから脳の血流量を推定することができる。
【００６０】
上述した本発明の実施形態は、本発明の説明のための例示であり、本発明の範囲をそれらの実施形態にのみ限定する趣旨ではない。当業者は、本発明の要旨を逸脱することなしに、他の様々な態様で本発明を実施することができる。
【００６１】
例えば、上記実施形態ではＳＰＥＣＴによる画像を利用しているが、ＰＥＴ（Ｐｏｓｉｔｏｒｏｎ  Ｅｍｉｓｓｉｏｎ  Ｔｏｍｏｇｒａｐｈｙ）による画像を利用することもできる。
【図面の簡単な説明】
【００６２】
【図１】本発明の一実施形態に係る血流量推定システムの構成図である。
【図２】撮影プロトコール及び対応する画像データを示す。
【図３】表示装置４に表示される肺のＲＯＩ入力画面の例を示す。
【図４】タイムアクティビティカーブの例を示す。
【図５】ΔＬ→Ｃａ変換用関数の説明図である。
【図６】表示装置４に表示される脳のＲＯＩ入力画面の例を示す。
【図７】被験者の脳血流を推定するための処理手順を示すフローチャートである。
【符号の説明】
【００６３】
１           血流量推定装置
２           撮影装置
３           入力装置
４           表示装置
１１    画像データ記憶部
１３       関数記憶部
２０       Ｃａ推定部
２１       ＲＯＩ決定部
２３       ΔＬ算出部
２５       変換処理部
３０       マイクロスフェア処理部
３１       ＲＯＩ決定部
３３       ＳＰＥＣＴ変換部
３５       ＣＢＦ算出部
【図１】
2008185499.tif 000006
【図２】
2008185499.tif 000007
【図３】
2008185499.tif 000008
【図４】
2008185499.tif 000009
【図５】
2008185499.tif 000010
【図６】
2008185499.tif 000011
【図７】
2008185499.tif 000012