(57)【要約】【課題】Ｘ線ＣＴシステムを利用した血流速度の変化を測定する技術を提供すること。【解決手段】マルチスライス型のＸ線ＣＴシステムにおいて、測定対象の血管に造影剤を投与した後にシネスキャンを実施する。シネスキャンの実施中、造影剤は生理的に血管内を流れて、各検出器列を通過するので、造影剤が第１の検出器列に対応するスライス位置に到達した時刻と、その造影剤が後に別のスライス位置に到達した時刻を測定して、両者の時間差を造影剤が距離Ｌを流れるのにかかる時間として求め、これに基づき血流速度を求める。これを、それぞれの隣接スライス位置間に対して行い、求めた血流速度どうしの差分をとることで血流速度の変化を計算する。【選択図】  図３
2004208714.tif 000002
【特許請求の範囲】
【請求項１】
被検体載置用のテーブルの搬送方向に少なくとも３列のＸ線検出器列を備えたＸ線ＣＴシステムであって、
所定位置で前記テーブルの搬送を止めたままスキャンを繰り返し行うシネスキャン手段と、
前記シネスキャンによって収集された投影データに基づき、少なくとも２箇所の異なる検出器列間で、各検出器列に対応するスライス位置の断層像を再構成する再構成手段と、
前記再構成手段によって得られた各断層像に基づき、被検体の血管に投与された造影剤が各断層像のスライス位置に到達した時刻を推定する推定手段と、
前記推定手段によって推定された各時刻の時間差に基づいて、各検出器列間における血流速度をそれぞれ算出する算出手段と、
を有することを特徴とするＸ線ＣＴシステム。
【請求項２】
前記推定手段は、
前記再構成手段で得られた各断層像について測定対象の血管の領域を特定する手段と、
各断層像について前記測定対象の血管の領域における平均ＣＴ値を計算する手段と、
前記平均ＣＴ値がピークとなるスキャンの実施時刻を、前記造影剤が当該スライス位置に到達した時刻として推定する手段と、
を含むことを特徴とする請求項１に記載のＸ線ＣＴシステム。
【請求項３】
前記算出手段は、
２つのスライス位置間における血管経路の距離を求める距離計算手段を含み、前記血管経路の距離を前記時間差で割り算することで前記血流速度を算出することを特徴とする請求項１に記載のＸ線ＣＴシステム。
【請求項４】
前記距離計算手段は、前記２つのスライス位置それぞれの断層像における前記血管の位置の距離を検出し、その距離と前記２つのスライス位置の間隔とに基づいて前記血管経路の距離を求めることを特徴とする請求項３に記載のＸ線ＣＴシステム。
【請求項５】
被検体載置用のテーブルの搬送方向に少なくとも３列のＸ線検出器列を備えたＸ線ＣＴシステムの制御方法であって、
所定位置で前記テーブルの搬送を止めたままスキャンを繰り返し行うシネスキャンステップと、
前記シネスキャンによって収集された投影データに基づき、少なくとも２箇所の異なる検出器列間で、各検出器列に対応するスライス位置の断層像を再構成する再構成ステップと、
前記再構成ステップによって得られた各断層像に基づき、被検体の血管に投与された造影剤が各断層像のスライス位置に到達した時刻を推定する推定ステップと、
前記推定ステップによって推定された各時刻の時間差に基づいて、各検出器列間における血流速度をそれぞれ算出する算出ステップと、
を有することを特徴とするＸ線ＣＴシステムの制御方法。
【請求項６】
前記推定ステップは、
前記再構成ステップで得られた各断層像について測定対象の血管の領域を特定するステップと、
各断層像について前記測定対象の血管の領域における平均ＣＴ値を計算するステップと、
前記平均ＣＴ値がピークとなるスキャンの実施時刻を、前記造影剤が当該スライス位置に到達した時刻として推定するステップと、
を含むことを特徴とする請求項５に記載のＸ線ＣＴシステムの制御方法。
【請求項７】
前記算出ステップは、
２つのスライス位置間における血管経路の距離を求める距離計算ステップを含み、
前記血管経路の距離を前記時間差で割り算することで前記血流速度を算出することを特徴とする請求項５に記載のＸ線ＣＴシステムの制御方法。
【請求項８】
前記距離計算ステップは、前記２つのスライス位置それぞれの断層像における前記血管の位置の距離を検出し、その距離と前記２つのスライス位置の間隔とに基づいて前記血管経路の距離を求めることを特徴とする請求項７に記載のＸ線ＣＴシステムの制御方法。
【請求項９】
被検体載置用のテーブルの搬送方向に少なくとも３列のＸ線検出器列を備えたＸ線ＣＴシステムにおける操作コンソールに、
所定位置で前記テーブルの搬送を止めたままスキャンを繰り返し行うようＸ線ＣＴシステムにおけるガントリに指示する指示ステップ、
前記指示ステップに応じて前記ガントリより転送されてきた投影データに基づき、少なくとも２箇所の異なる検出器列間での、各検出器列に対応するスライス位置の断層像を再構成する再構成ステップ、
前記再構成ステップによって得られた各断層像に基づき、被検体の血管に投与された造影剤が各断層像のスライス位置に到達した時刻を推定する推定ステップ、
前記推定ステップによって推定された各時刻の時間差に基づいて、各検出器列間における血流速度をそれぞれ算出する算出ステップ、
を実行させるためのプログラム。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、Ｘ線ＣＴシステムおよびその制御方法に関する。
【０００２】
【従来の技術】
Ｘ線ＣＴ（Computerized Tomography）システムは、被検体に複数方向からＸ線を投影するスキャンを行い、被検体を透過した各方向からのＸ線より得られる投影データに基づいて画像再構成処理を行うことによって、診断部位の断層像を提供する。
【０００３】
従来より、被検体の血管に造影剤を投与した後にＸ線ＣＴシステムでスキャンを行うと血管が明瞭に描出されることが知られており、これを利用して血流速度の測定を行うことができる。Ｘ線ＣＴシステムを利用した血流速度の測定方法としては、例えば、Ｘ線検出器を複数列備えたいわゆるマルチスライス型のＸ線ＣＴシステムにおいてシネスキャンを行い、２つの検出器列のそれぞれに対応するスライス位置の断層像に基づきそれぞれのスライス位置に造影剤が到達した時刻を推定し、それら到達時刻の時間差から血流速度を求めるものが提案されている（例えば特許文献１を参照。）。
【０００４】
【特許文献１】
特開平１０−２８６２５２号公報
【０００５】
【発明が解決しようとする課題】
しかしながら、上記のようにＸ線ＣＴシステムを利用して血流速度を測定することはできるものの、各検出器間の血流速度の変化を示すことはできなかった。
【０００６】
また、血流速度の測定精度、および血流速度変化の測定精度をさらに高めたいという要請もある。
【０００７】
そこで、本発明は、Ｘ線ＣＴシステムを利用した血流速度の変化を示す技術を提供することを目的とする。
【０００８】
また、本発明は、Ｘ線ＣＴシステムを利用した血流速度の測定精度を高めることを目的とする。
【０００９】
【課題を解決するための手段】
本発明の一側面によれば、被検体載置用のテーブルの搬送方向に少なくとも３列のＸ線検出器列を備えたＸ線ＣＴシステムであって、所定位置で前記テーブルの搬送を止めたままスキャンを繰り返し行うシネスキャン手段と、前記シネスキャンによって収集された投影データに基づき、少なくとも２箇所の異なる検出器列間で、各検出器列に対応するスライス位置の断層像を再構成する再構成手段と、前記再構成手段によって得られた各断層像に基づき、被検体の血管に投与された造影剤が各断層像のスライス位置に到達した時刻を推定する推定手段と、前記推定手段によって推定された各時刻の時間差に基づいて、各検出器列間における血流速度をそれぞれ算出する算出手段と、を有することを特徴とするＸ線ＣＴシステムが提供される。
【００１０】
本発明の好適な実施形態によれば、前記推定手段は、前記再構成手段で得られた各断層像について測定対象の血管の領域を特定する手段と、各断層像について前記測定対象の血管の領域における平均ＣＴ値を計算する手段と、前記平均ＣＴ値がピークとなるスキャンの実施時刻を、前記造影剤が当該スライス位置に到達した時刻として推定する手段と、を含むことが好ましい。
【００１１】
また、本発明の好適な実施形態によれば、前記算出手段は、２つのスライス位置間における血管経路の距離を求める距離計算手段を含み、前記血管経路の距離を前記時間差で割り算することで前記血流速度を算出することが好ましい。
【００１２】
さらに、本発明の好適な実施形態によれば、前記距離計算手段は、前記２つのスライス位置それぞれの断層像における前記血管の位置の距離を検出し、その距離と前記２つのスライス位置の間隔とに基づいて前記血管経路の距離を求めることが好ましい。
【００１３】
本発明の他の側面によれば、被検体載置用のテーブルの搬送方向に少なくとも３列のＸ線検出器列を備えたＸ線ＣＴシステムの制御方法であって、所定位置で前記テーブルの搬送を止めたままスキャンを繰り返し行うシネスキャンステップと、前記シネスキャンによって収集された投影データに基づき、少なくとも２箇所の異なる検出器列間で、各検出器列に対応するスライス位置の断層像を再構成する再構成ステップと、前記再構成ステップによって得られた各断層像に基づき、被検体の血管に投与された造影剤が各断層像のスライス位置に到達した時刻を推定する推定ステップと、前記推定ステップによって推定された各時刻の時間差に基づいて、各検出器列間における血流速度をそれぞれ算出する算出ステップと、を有することを特徴とするＸ線ＣＴシステムの制御方法も提供される。
【００１４】
【発明の実施の形態】
以下、図面を参照して実施形態について詳細に説明する。
【００１５】
図１は、実施形態におけるＸ線ＣＴシステムの構成を示す図である。
【００１６】
図示の如く、本システムは、被検体へのＸ線照射と被検体（患者）を透過したＸ線を検出するガントリ１００と、ガントリ１００に対して各種動作設定を行うとともに、ガントリ１００から出力されてきたデータに基づいて断層像を再構成し出力（表示）する操作コンソール２００を含む構成である。
【００１７】
ガントリ１００は、その全体の制御をつかさどるメインコントローラ１をはじめ以下の構成を備える。
【００１８】
２ａおよび２ｂは操作コンソール２００との通信を行うためのインタフェース、３はテーブル１１上に横たえた被検体を図面に垂直な方向（一般に患者の体軸の方向に一致する方向であり、以下、ｚ軸方向という。）に搬送するための空洞部を有する回転部であり、その内部には、Ｘ線発生源であるＸ線管１４（Ｘ線管コントローラ５により駆動が制御される）、Ｘ線の照射範囲を制限するための開口を有するコリメータ６、コリメータ６のｚ軸方向の開口幅を調整するための開口制御モータ７（開口制御モータドライバ８により駆動が制御される）が設けられている。
【００１９】
また、回転部３には、コリメータ６および空洞部を経由してきたＸ線管１４からのＸ線を検出するための複数（例えば1,000個）の検出チャネルを有するＸ線検出部４、および、Ｘ線検出部４の各検出チャネルの出力に基づき投影データとして収集するデータ収集部１５も備える。Ｘ線管１４およびコリメータ６とＸ線検出部４とは、互いに空洞部を挟んで、すなわち、被検体を挟んで、対向する位置に設けられる。回転部３は、その位置関係を維持した状態で空洞部の周りを回転するように構成されている。この回転は、回転モータドライバ１０からの駆動信号により駆動される回転モータ９によって行われる。また、被検体を載置するテーブル１１は、ｚ軸方向への搬送がなされるが、その駆動は、テーブルモータドライバ１３からの駆動信号により駆動されるテーブルモータ１２によって行われる。
【００２０】
メインコントローラ１は、インタフェース２ａを介して受信した各種コマンドの解析を行い、それに基づいて上記のＸ線管コントローラ５、開口制御モータドライバ８、回転モータドライバ１０、テーブルモータドライバ１３、およびデータ収集部１５に対し、各種制御信号を出力することになる。
【００２１】
また、データ収集部１５で収集されたデータは、インタフェース２ｂを介して操作コンソール２００に送出される。
【００２２】
一方、操作コンソール２００は、いわゆるワークステーションであり、図示するように、装置全体の制御を司るＣＰＵ５１、ブートプログラム等を記憶しているＲＯＭ５２、主記憶装置として機能するＲＡＭ５３をはじめ、以下の構成を備える。
【００２３】
ＨＤＤ５４は、ハードディスク装置であって、ここにＯＳのほか、ガントリ１００に各種指示を与えたり、ガントリ１００より受信したデータに基づく断層像の再構成および表示、ならびに血流速度の測定を実行させるための画像処理プログラムが格納されている。また、ＶＲＡＭ５５は表示しようとするイメージデータを展開するメモリであり、ここにイメージデータ等を展開することでＣＲＴ５６に表示させることができる。５７および５８はそれぞれ、各種設定を行うためのキーボードおよびマウスである。また、５９および６０はガントリ１００と通信を行うためのインタフェースであり、それぞれガントリ１００のインタフェース２ａおよび２ｂに接続される。
【００２４】
実施形態におけるＸ線ＣＴシステムの構成は概ね上記のとおりであるが、次にＸ線管１４、コリメータ６、Ｘ線検出部４の構造を図２を用いてより詳しく説明する。
【００２５】
図２はＸ線管１４、コリメータ６、Ｘ線検出部４の要部構成図である。これらの構成要素は回転部３の所定の基部に支持されて図示のような位置関係を維持している。
【００２６】
同図において、Ｘ線管１４は、集束電極およびフィラメントを内蔵する陰極スリーブ（sleeve）４２と、回転するターゲット４３とをハウジング４１に内蔵した構造であり、焦点ｆからＸ線を放射する。そして、Ｘ線管１４から放射されたＸ線がコリメータ６が形成するスリット６ａを通過させることによって、所定のＸ線照射角（ファン角）のファンビームが形成される。上記したＸ線検出部４は、ファン角に依存した長さにわたる個数（例えば、1,000個）の検出チャネルで一列の検出器を形成し、この検出器がｚ軸方向（テーブル１１の搬送方向に一致する。）に例えば８列配列された構成である。ここでは、各検出器列には中央から外側に向けて１，２，３，４の番号を付すとともに、右４列の検出器列にはＡ、左４列の検出器列にはＢの符号を併記する。つまり、Ｘ線検出器４は、中央から１Ａ，１Ｂ，２Ａ，２Ｂ，３Ａ，３Ｂ，４Ａ，４Ｂの８列の検出器列で構成される。これにより、いわゆるマルチスライスＸ線ＣＴを実現する。
【００２７】
かかる構成のＸ線ＣＴシステムにおいて、投影データの収集は例えば次のように行われる。
【００２８】
まず、被検体を回転部３の空洞部に位置させた状態でｚ軸方向の位置を固定し、Ｘ線管１４からのＸ線ビームを被検体に照射し（Ｘ線の投影）、その透過Ｘ線をＸ線検出部４で検出する。そして、この透過Ｘ線の検出を、Ｘ線管１４とＸ線検出部４を被検体の周囲を回転させながら（すなわち、投影角度（ビュー角度）を変化させながら）複数N（例えば、N=1,000）のビュー方向で、360度分行う。
【００２９】
検出された各透過Ｘ線は、データ収集部１５でディジタル値に変換されて投影データとしてインタフェース２ｂを介して操作コンソール２００に転送される。これら一連の工程を１つの単位として１スキャンとよぶ。そして、順次ｚ軸方向にスキャン位置を所定量移動して、次のスキャンを行っていく。このようなスキャン方式はアキシャルスキャン方式とよばれるが、実施形態におけるガントリ１００は、テーブル１１を所望の位置に固定したまま所定時間スキャンを繰り返す、いわゆるシネスキャンを行うこともできる。このシネスキャンは造影剤を投与した血管の血流解析に好適である。
【００３０】
操作コンソール２００は、ガントリ１００から転送されてくる投影データをＨＤＤ５４に格納するとともに、例えば、所定の再構成関数とたたみ込み演算を行い、バックプロジェクション処理により断層像を再構成する。ここで、操作コンソール２００は、スキャン処理中にガントリ１００から順次転送されてくる投影データからリアルタイムに断層像を再構成し、常に最新の断層像をＣＲＴ５６に表示させることが可能である。さらに、ＨＤＤ５４に格納されている投影データを呼び出して改めて画像再構成を行わせることも可能である。
【００３１】
続いて、実施形態における血流速度の変化の測定処理の概要を、図３の模式図を用いて説明する。
【００３２】
図３は、実施形態におけるＸ線ＣＴシステムによるスキャンを行うときの被検体とＸ線検出器４とを模式的にｚ軸に沿って表示したものである。ここで、それぞれ隣接する検出器列間の距離をｄ、検出器列４Ｂ−４Ａ間の距離をＬとする。スキャンを行うと、スライス位置ｚ１〜ｚ８の、各検出器列に対応した断層像を再構成することができる。
【００３３】
いま、概ねｚ軸に沿って伸びている測定対象の血管に、図面左側から造影剤を投与した後にシネスキャンを実施する。シネスキャンの実施中、造影剤は生理的に血管内を流れて、検出器列４Ｂに対応するスライス位置ｚ１に到達し、順次隣のスライス位置を通過してｚ８に到達することになる。そうすると、造影剤がスライス位置ｚ１に到達した時刻と、スライス位置ｚ８に到達した時刻が分かれば、両者の時間差が造影剤が距離Ｌを流れるのにかかる時間として求められ、これによりマルチスライスにおける両端のスライス間の血流速度を求めることができる。各スライス位置に造影剤が到達したことは、そのスライス位置の断層像の各スキャンにおけるＣＴ値の変化により判定することが可能である。このような血流速度の測定によれば、造影剤が初めに所定のスライス位置から別のスライス位置に到達するまでの時間だけシネスキャンを実施すればよく、投与したすべての造影剤がスキャン位置を通過するまで実施する必要はない。
【００３４】
同様に、隣接する２つのスライス位置それぞれに造影剤が到達した時刻が分かれば、両者の時間差が造影剤が距離ｄを流れるのにかかる時間として求められ、これにより隣接スライス間の血流速度を求めることができる。そうすると、例えば、ｚ１−ｚ２間の血流速度、ｚ２−ｚ３間の血流速度、ｚ３−ｚ４間の血流速度、…、ｚ７−ｚ８間の血流速度をそれぞれ求めることができ、例えば、縦軸を血流速度、横軸をスライス位置とするグラフ上に、各々の血流速度をプロットし、更に隣り合うプロットを結んで折れ線表示すると、血流速度の変化を知ることができる。
【００３５】
上記したような血流速度変化の測定を実現する、実施形態における処理をフローチャートを用いて詳しく説明する。
【００３６】
図４は、実施形態におけるＸ線ＣＴシステムによる血流速度変化の測定処理の概略を示すフローチャートである。このフローチャートに対応するプログラムは、ＨＤＤ５４にインストールされている画像処理プログラムに含まれ、ＲＡＭ５３にロードされてＣＰＵ５１によって実行されるものである。
【００３７】
まず、血流速度の測定を行うスライス位置またはその付近で少なくとも１回のアキシャルスキャンを行うようガントリ１００に指示を出す（ステップＳ１）。このステップには実際には、スライス厚、スキャン開始／終了位置、Ｘ線管４に与える電流値、ガントリ３の回転速度等のスキャン条件を設定するいわゆるスキャン計画工程が含まれるが、このスキャン計画そのものは公知のものであるので、その詳細は省略する。ガントリ１００はこれに応じて、アキシャルスキャンを実行する。
【００３８】
そして、ガントリ１００より転送されてくる投影データを受信し（ステップＳ２）、その投影データに基づき所定の画像再構成アルゴリズム（例えば、フィルタード・バックプロジェクション）に従い断層像の再構成を行い、その断層像をＣＲＴ５６に表示する（ステップＳ３）。ここでオペレータは、キーボード５７またはマウス５８を用いて、表示された断層像の着目したい血管周辺領域をＲＯＩ（Region Of Interest）として指定することができる。ステップＳ４では、キーボード５７またはマウス５８からの入力に基づいて、ＲＯＩを設定する。
【００３９】
次に、オペレータは被検体に造影剤を投与する。その後、造影剤が着目部位に到達すると思われる時刻の少し前の時刻で、オペレータによるキーボード５７またはマウス５８の入力に基づき、シネスキャンの開始をガントリ１００に指示する（ステップＳ５）。ガントリ１００はこの指示に応じてシネスキャンを実行する。このステップＳ５は、シネスキャンを行う位置やシネスキャンの継続時間等の設定工程を含む。ここではシネスキャンの継続時間は、上述のとおり、造影剤が初めに検出器列４Ｂに対応するスライス位置ｚ１から検出器列４Ａに対応するスライス位置ｚ８に到達するのに必要十分と思われる時間を設定すればよい。
【００４０】
そして、ガントリ１００より転送されてくる投影データを受信し（ステップＳ６）、その投影データに基づいて血流速度の計算を行う（ステップＳ７）。このステップＳ７の血流速度の計算は、それぞれの隣接スライス位置間に対して行い、例えば求めた血流速度どうしの差分をとることで血流速度の変化を計算する（ステップＳ８）。
【００４１】
図５は、上記ステップＳ７の血流速度の計算ステップの詳細を示すフローチャートである。
【００４２】
まず、断層像１枚分の投影データを取得し（ステップＳ７１）、次に、各検出器列に対応するスライス位置の断層像を再構成する（ステップＳ７２）。その後、各断層像について、ステップＳ４で設定したＲＯＩから、測定対象の血管を特定する（ステップＳ７３）。
【００４３】
このステップ７３は具体的には、例えば図６のフローチャートに示すような方法が考えられる。一般に、造影されていない血管のＣＴ値はおよそ４０〜５０くらいであるが、造影された血管は１００〜２００程度のＣＴ値を有する。そこで、ステップＳ７３１で、ＲＯＩ内のＣＴ値について、１００以上２００未満のＣＴ値を１、それ以外を０とする２値化処理を行い、ステップＳ７３２で、連続領域抽出のためにラベリングを行う。そして、ステップＳ７３３で、各ラベルの面積Ｓを測定し、大動脈を想定して、例えば、Ｓが１．８平方センチメートル以上２．２平方センチメートル未満のラベル領域を、測定対象の血管とする。
【００４４】
そして、各断層像について、測定対象の血管領域の平均ＣＴ値を計算する（ステップＳ７４）。続くステップＳ７５では、再構成すべき投影データがなくなったか否かを判定し、まだある場合にはステップＳ７６で次の断層像に係る投影データを取得した後ステップＳ７２に戻り、処理を繰り返す。そして、すべての投影データについてステップＳ７２〜Ｓ７４の処理を終えデータがなくなった時点で、ステップＳ７７に進む。
【００４５】
上述したとおり、一般には、造影されていない血管のＣＴ値はおよそ４０〜５０くらいであるが、造影された血管は１００〜２００程度のＣＴ値を有する。したがって、一連のシネスキャンを行う間、そのスライス位置にまだ造影剤が届いていないとき、および、造影剤が通過した後は、そのスライス位置での測定対象の血管領域の平均ＣＴ値は４０〜５０程度の低い値となり、造影剤が到達したときは、その平均ＣＴ値が１００〜２００程度の値のピークを有するようになる。後続のスライス位置でもある程度の時間遅れで同様の平均ＣＴ値の変動を観測することができる。
【００４６】
これを図に示したものが、図７である。同図において、横軸がシネスキャンの継続時間で、縦軸が対象血管領域の平均ＣＴ値を示している。同図に示すように、シネスキャンを開始すると、まず、検出器列４Ｂに対応するスライスにおける平均ＣＴ値にピークが現れ、続いて検出器列３Ｂに対応するスライス、検出器列２Ｂに対応するスライス、・・・、以下、検出器列４Ａに対応するスライスの順で、平均ＣＴ値にピークが現れる。
【００４７】
上記のとおり、このピークをもってそのスライス位置に造影剤が到達した時刻を検出することができる。これにより、例えば、検出器列４Ｂに対応するスライス位置に造影剤が到達したことを示すピーク時刻と、隣接する検出器列３Ｂに対応するスライス位置に造影剤が到達したことを示すピーク時刻の時間遅れΔｔを求めれば、血流速度を求めることができる。
【００４８】
そこで、ステップＳ７７で、いま対象とする隣接する２つの検出器列（例えば検出器列４Ｂおよび３Ｂ）それぞれの断層像について、対象血管領域の平均ＣＴ値がピークとなるときのシネスキャンの実施時刻を測定し、両者の時間差Δｔを求める。そして、次式により血流速度Ｆを計算する（ステップＳ７８）。
【００４９】
Ｆ＝ｄ／Δｔ
ただし、ｄは上記のとおり隣接検出器列間の距離である。
【００５０】
このように、隣接する２つの検出器列それぞれに対応するスライス位置の断層像に基づいて、それぞれのスライス位置に造影剤が到達した時刻を推定し、それら到達時刻の時間差から血流速度が求められる。
【００５１】
このステップＳ７では、上記の処理を例えばすべての検出器列について行い、隣接する検出器列間４Ｂ−３Ｂ、３Ｂ−２Ｂ、２Ｂ−１Ｂ、１Ｂ−１Ａ、１Ａ−２Ａ、２Ａ−３Ａ、３Ａ−４Ａのそれぞれに係る血流速度を測定する。
【００５２】
以上説明した処理によって、血流速度の変化を求めることができる。
【００５３】
ところで、上記の例では、血流速度Ｆを、Ｆ＝ｄ／Δｔにより求めた。隣接検出器列間の距離ｄをΔｔで割るという演算は、測定対象の血管がｚ軸に沿って伸びているという前提の下でのことである。
【００５４】
しかし厳密には、測定対象の血管がｚ軸と平行に直線的に伸びているとは限らない。測定対象の血管がｚ軸に対して斜めに伸びている場合には、上記の測定方法で求めた血流速度の誤差が大きくなってしまう。
【００５５】
そこで、本実施形態では、２つのスライス位置それぞれの断層像に基づき測定対象の血管の両スライス間における長さを精度よく求め、その長さに基づいて血流速度を計算することができる。
【００５６】
以下、図８の模式図を参照して詳しく説明する。
【００５７】
図８は、検出器列４Ｂに対応するスライス１と、その隣の検出器列３Ｂに対応するスライス２とを模式的に示している。ｚ軸方向に形成される各スライス面において、横方向をｘ軸、上下方向をｙ軸とする。また、隣接するスライス位置の間隔（スライス厚）は、隣接検出器列間の距離ｄで近似する。
【００５８】
いま、スライス１においては、測定対象の血管の中心Ｐ１が（ｘ１，ｙ１）の位置で検出され、スライス２においては、測定対象の血管の中心Ｐ２が（ｘ２，ｙ２）の位置で検出された場合を考える。ここで、測定対象の血管経路を直線Ｇで近似することにすれば、血管経路Ｇの距離ΔＬは次のように求めることができる。
【００５９】
まず、スライス１における血管の中心Ｐ１をスライス２に射影して、スライス２におけるＰ１とＰ２の距離ΔＰを求める。ΔＰは次式により求めることができる。
【００６０】
ΔＰ＝sqrt（（ｘ２−ｘ１）２＋（ｙ２−ｙ１）２）
ただし、sqrt（ｎ）は、ｎの平方根である。
【００６１】
そして、次式により、血管経路Ｇの距離ΔＬを求める。
【００６２】
ΔＬ＝sqrt（ΔＰ２＋ｄ２）
【００６３】
このようにして、隣接するスライス間における距離ΔＬが計算された。このΔＬを、血流速度Ｆの計算に用いる。すなわち、Ｆ＝ΔＬ／Δｔとして計算すれば、さらに高い精度で血流速度を求めることができる。
【００６４】
なお、以上の実施形態におけるＸ線ＣＴシステムは８列の検出器列を備えたものであったが、本発明は少なくとも３列の検出器列を備えていれば実現することが可能であり、８列に限定されるものではない。
【００６５】
また、以上説明した実施形態によれば、本発明はＸ線ＣＴシステムにおける操作コンソール２００の制御処理によって実現されたが、操作コンソール２００の構成自体は汎用の画像処理装置（ワークステーションやパーソナルコンピュータ等）で実現できるものであるので、ソフトウェアを同装置にインストールし、それでもって実現することが可能である。
【００６６】
したがって、本発明の機能処理をコンピュータで実現するための、コンピュータにインストールされるプログラム自体および、そのプログラムを格納した記録媒体そのものも本発明を実現するものである。つまり、本発明の特許請求の範囲には、本発明の機能処理を実現するためのプログラム自体および、そのプログラムを格納したコンピュータ読み取り可能な記録媒体も含まれうる。
【００６７】
プログラムを供給するための記録媒体としては、例えば、フレキシブルディスク、光ディスク（CD-ROM、CD-R、CD-RW、DVD等）、光磁気ディスク、磁気テープ、メモリカード等がある。
【００６８】
その他、プログラムの供給方法としては、インターネットを介して本発明のプログラムをファイル転送によって取得する態様も含まれる。
【００６９】
【発明の効果】
以上説明したように、本発明によれば、Ｘ線ＣＴシステムを利用した血流速度の変化を測定する技術を提供することができる。
【００７０】
また、本発明によれば、Ｘ線ＣＴシステムを利用した血流速度の測定精度を高めることができる。
【図面の簡単な説明】
【図１】実施形態におけるＸ線ＣＴシステムの構成を示す図である。
【図２】実施形態におけるＸ線管、コリメータ、Ｘ線検出部の要部構成図である。
【図３】実施形態における血流速度の測定処理の概要を説明するための図である。
【図４】実施形態におけるＸ線ＣＴシステムによる血流速度の測定処理の概略を示すフローチャートである。
【図５】血流速度の計算ステップの詳細を示すフローチャートである。
【図６】実施形態における測定対象の血管を特定する処理を示すフローチャートである。
【図７】各検出器列に対応するスライスにおける対象血管領域の平均ＣＴ値の変動を表す図である。
【図８】隣接するスライス間における測定対象の血管の距離を計算する過程を説明するための図である。
【図１】
2004208714.tif 000003
【図２】
2004208714.tif 000004
【図３】
2004208714.tif 000005
【図４】
2004208714.tif 000006
【図５】
2004208714.tif 000007
【図６】
2004208714.tif 000008
【図７】
2004208714.tif 000009
【図８】
2004208714.tif 000010