(57)【特許請求の範囲】
【請求項１】
放射線を照射する放射線管と、前記所定の軸まわりと当該軸に沿う方向との２次元的に配列された複数の放射線検出器と、前記放射線管および前記複数の放射線検出器を前記所定の軸まわりに回転、および、当該軸に沿う方向に前記被検部位に対して相対的に移動させる移動手段とを有し、被検部位を通過した放射線を前記被検部位の所定の軸まわりに検出するスキャンを、前記所定の軸方向に沿って異なる位置において複数回行なう放射線検出装置と、
複数回の前記スキャンの各々について当該スキャンによって得られた検出データに基づく計算により前記所定の軸方向における画像生成位置の被検部位の断層像データを再構成し、当該断層像データに対して各々重みを付与し、重みを付与した当該複数の断層像データを合成して前記画像生成位置の合成断層像データを生成する画像再構成手段と
を有し、
前記各々の断層像データに対する重みを付与において、前記断層像データの各画素データのうち、各々の前記スキャンによる前記検出データが存在する領域以外の領域の画素データに付与する前記重みを０とする
放射線計算断層画像装置。
【請求項２】
前記画像生成位置における前記所定の軸から断層像の各画素までの距離および／または前記所定の軸上の所定の位置から前記画像生成位置までの距離に応じて、前記複数の断層像データに各々付与する前記重みを変化させる
請求項１に記載の放射線計算断層画像装置。
【請求項３】
前記断層像データの前記画素データの各々に対して付与する重みを、前記所定の軸を中心とした同心円状に変化させる
請求項１又は２に記載の放射線計算断層画像装置。
【請求項４】
前記画像再構成手段が前記所定の軸方向に沿って生成した複数の前記合成断層像データに基づいて、所望の方向の断層画像を生成する再構成画像生成手段をさらに有する
請求項１〜３のいずれか一項に記載の放射線計算断層画像装置。
【請求項５】
前記放射線はＸ線である、請求項１〜４のいずれか一項に記載の放射線計算断層画像装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、対象物に対する放射線の投影データから対象物の断層像を計算する放射線計算断層画像装置（以下、放射線CT（Computed Topography）装置）、および、断層画像生成方法に関する。
【０００２】
【従来の技術】
放射線CT装置として、たとえば、放射線にX線を用いるX線CT装置が知られている。このX線CT装置として、複数のX線検出器を２次元的に配列した装置が知られている。複数のX線検出器は、被検体に対する所定の軸に沿った方向に幅を持つように配置される。X線検出器の列が軸方向に沿って所定の幅分形成されるため、２次元的に配列したX線検出器は、多列検出器と呼ばれる。
【０００３】
多列検出器を用いた断層像再構成手法として、コーンBP（Back Projection）法が知られている（たとえば、特許文献１参照。）。
多列検出器を有するX線CT装置においてコーンBP法を適用する場合には、たとえば、上述の軸に直交する断面が扇形のコーン状のX線を放射しながら、多列検出器を軸まわりに１回転させる。そして、被検体を通過したX線を多列検出器により検出して被検体の被検部位をスキャンする。１回のスキャン終了後、多列検出器を軸方向に所定ピッチずらして、次のスキャンを行なう。このようなスキャン手法は、アキシャル（axial）スキャンと呼ばれる。X線検出器の列が軸方向に複数存在するため、X線検出器の列の数に応じて、１回のスキャンによって被検体の断層像を複数枚得ることができる。
コーンBP法においては、あるX線検出器の列に対応する断層像を生成する場合に、対象とするX線検出器の列によって得た検出データだけでなく、断層像に対応する被検部位を通過し他のX線検出器の列によって検出されたX線の検出データも利用する。このように画像生成のためのデータが増えることによって、断層像の画質をある程度向上させることができる。
【０００４】
【特許文献１】
特開２００２−３３６２３９号公報
【０００５】
【発明が解決しようとする課題】
しかしながら、従来のコーンBP法では、軸方向の端部の検出器の列に相当する断層像生成において、画像生成のために利用可能な検出データが減少する。このため、特に端部の検出器の列に相当する断層像に対しては、さらなる画質の向上が望まれている。
【０００６】
また、従来のコーンBP法では、被検体の軸まわりの１回のスキャンによって得られた検出データのみを用いて断層像を生成する。このため、アキシャルスキャンを行なう場合に、軸方向のある位置におけるスキャンによって得られた断層像と、他の位置における他のスキャンによって得られた断層像との間ではデータの連続性が悪い。データの連続性が悪いと、コーンBP法によって得られた軸方向に直交する複数の断層像を用いた断面変換により得られる他の断層像の画質に悪影響が及ぶ等の不利益が発生する。不利益の例として、断面変換により得られるMPR（Multi Plain Reformation）画像が十分に滑らかではなく、スキャンとスキャンとの境界に相当する部分において縞状のアーチファクト（artifact）が発生することが挙げられる。
【０００７】
したがって、本発明の目的は、アキシャルスキャンによる断層像再構成において、被検体の軸方向における所定の画像生成位置の断層像の画質を向上可能な放射線計算断層画像装置（放射線CT装置）を提供することにある。また、本発明の目的は、この放射線CT装置において用いられ、画像生成位置の断層像の画質を向上可能な断層画像生成方法を提供することにもある。
【０００８】
【課題を解決するための手段】
本発明に係る放射線計算断層画像装置は、被検部位を通過した放射線を前記被検部位の所定の軸まわりに検出するスキャンを、前記所定の軸方向に沿って異なる位置において複数回行なう放射線検出装置と、複数回の前記スキャンの各々によって得られた検出データに基づいて、前記被検部位の同じ位置に関して複数のデータを求め、当該複数のデータに基づいて前記所定の軸方向における画像生成位置の被検部位の１つの断層像データを再構成する画像再構成手段とを有する。
【０００９】
本発明に係る放射線計算断層画像装置は、被検部位を通過した放射線を前記被検部位の所定の軸まわりに検出するスキャンを、前記所定の軸方向に沿って異なる位置において複数回行なう放射線検出装置と、複数回の前記スキャンの各々について当該スキャンによって得られた検出データに基づく計算により前記所定の軸方向における画像生成位置の被検部位の断層像データを再構成し、再構成した複数の前記断層像データを合成して前記画像生成位置の合成断層像データを生成する画像再構成手段とを有する構成でもよい。
【００１０】
また、本発明に係る断層画像生成方法は、被検部位を通過した放射線を前記被検部位の所定の軸まわりに検出するスキャンを、前記所定の軸方向に沿って異なる位置において複数回行なうステップと、複数回の前記スキャンの各々によって得られた検出データに基づいて、前記被検部位の同じ位置に関して複数のデータを求め、当該複数のデータに基づいて前記所定の軸方向における画像生成位置の被検部位の１つの断層像データを計算し再構成するステップとを有する。
【００１１】
本発明に係る断層画像生成方法は、被検部位を通過した放射線を前記被検部位の所定の軸まわりに検出するスキャンを、前記所定の軸方向に沿って異なる位置において複数回行なうステップと、複数回の前記スキャンの各々について当該スキャンによって得られた検出データに基づいて前記所定の軸方向における画像生成位置の被検部位の断層像データを計算し再構成するステップと、再構成した複数の前記断層像データを合成して前記画像生成位置の合成断層像データを生成するステップとを有する構成でもよい。
【００１２】
本発明においては、被検部位を通過した放射線を放射線検出装置が被検部位の所定の軸まわりに検出して被検部位をスキャンする。放射線検出装置は、被検部位のスキャンを軸方向に沿って異なる位置において複数回行なう。複数回のスキャンの各々に対して、被検部位を通過した放射線の検出データが得られる。画像再構成手段が、複数回のスキャンの各々によって得られた検出データに基づいて、被検部位の同じ位置に関して複数のデータを求める。画像再構成手段は、この複数のデータに基づいて、軸方向における所定の画像生成位置の被検部位の１つの断層像データを計算して再構成する。
その際に、画像再構成手段は、複数回のスキャンの各々について、得られた検出データから画像生成位置の被検部位の断層像データを再構成する。画像再構成手段は、画像生成位置においてスキャンの回数に応じて得られる複数の断層像データを合成して合成断層像データを生成する。この合成断層像データが、画像生成位置における被検部位の最終的な断層像データになる。
【００１３】
【発明の実施の形態】
以下、本発明の実施の形態について、添付図面を参照しながら述べる。本発明における放射線は、X線を含む。以下では、放射線CT装置としてX線CT装置を例に挙げて述べる。
【００１４】
装置構成
図１は、本発明に係るX線CT装置の一実施の形態の装置構成を示す図であり、図２は、図１に示すX線CT装置の要部の構成を示す図である。図１に示すX線CT装置１が、本発明における放射線計算断層画像装置の一実施の形態である。
【００１５】
図１に示すように、本実施の形態に係るX線CT装置１は、走査ガントリ２と、操作コンソール３と、撮影テーブル４とを有する。
また、走査ガントリ２は、X線管２０と、移動部２１と、コリメータ２２と、検出器アレイ２３と、データ収集システム（Data Acquisition System：DAS）２４と、X線コントローラ２５と、コリメータコントローラ２６と、回転部２７と、回転コントローラ２８とを有している。
【００１６】
X線管２０が本発明における放射線管の一実施態様であり、移動部２１と回転部２７と回転コントローラ２８とを含んで本発明における移動手段の一実施態様が構成される。
また、X線管２０と、移動部２１と、検出器アレイ２３と、回転部２７と、回転コントローラ２８とを含んで、本発明における放射線検出装置の一実施態様が構成される。
【００１７】
図２に示すように、X線コントローラ２５は、X線管２０およびX線管移動部２１に接続される。コリメータコントローラ２６は、コリメータ２２に接続される。検出器アレイ２３はDAS２４に接続され、回転コントローラ２８は回転部２７に接続される。
【００１８】
X線管２０は、X線コントローラ２５からの制御信号CTL２５１に基づいて所定強度のX線をコリメータ２２に向かって放射する。
【００１９】
コリメータ２２は、コリメータコントローラ２６からの制御信号CTL２６１に基づいてアパーチャ２２１の開度を調節することによって、X線管２０から放射されたX線の照射範囲を調節する。
本実施の形態においては、後述する検出器アレイの２次元的なX線検出面の全面にX線が入射するようにアパーチャ２２１の開度が調節される。
【００２０】
検出器アレイ２３は、X線検出器を複数用いてi行j列の２次元のマトリクス（アレイ）状に構成したものである。
各X線検出器は、たとえば、シンチレータ（scintillator）とフォトダイオードの組み合わせによって構成する。
行方向には、たとえばi＝１０００個程度のX線検出器が配置されて１つの検出器列を構成する。図２においては、たとえば、j＝８として８列の検出器列を有する検出器アレイ２３を示している。以下では、各検出器列を表わす場合に、検出器アレイ２３の列方向の中心から右側に向かって順に検出器列１a，２a，３a，４a、中心から左側に向かって順に検出器列１b，２b，３b，４bと表わす。
複数の検出器列１a〜４a，１b〜４bは、隣接して互いに平行に配設する。
以上のように、X線検出器をマトリクス状に配置することにより、全体として２次元的なX線検出面２３Sが構成される。
【００２１】
図１および図２に示すように、検出器アレイ２３は、コリメータ２２から所定距離離れて配置する。コリメータ２２と検出器アレイ２３との間が、X線照射空間２９となる。このX線照射空間２９内に被検体６が搬入される。
検出器アレイ２３は、コリメータ２２のアパーチャ２２１の部分から照射され、被検体６が搬入されているX線照射空間２９を通過したX線ビーム５の強度をX線検出面２３Sによって検出する。
【００２２】
DAS２４は、操作コンソール３からの制御信号CTL３０３に基づいて、X線検出面２３Sを構成する個々のX線検出器からのX線強度の検出データを収集し、操作コンソール３に送信する。DAS２４によって収集される検出データは、ローデータ（Raw Data）と呼ばれる。
DAS２４は、ローデータをAD（Analog to Digital）変換して操作コンソール３に送信する。
【００２３】
回転部２７は、回転コントローラ２８からの制御信号に基づいて、X線照射空間２９内の所定の回転軸AXまわりに回転する。
回転部２７にはX線管２０、移動部２１、コリメータ２２、検出器アレイ２３、DAS２４、X線コントローラ２５、およびコリメータコントローラ２６が搭載されており、回転部２７の回転によってこれらの構成要素が相対的な位置関係を保ったまま回転軸AXまわりに回転する。
【００２４】
本実施の形態においては、回転軸AXの方向に沿って、検出器アレイ２３の各検出器列を配列する。そして、この回転軸AXの方向、すなわち検出器アレイ２３の検出器列１a〜４a，１b〜４bの方向を、z方向とする。
なお、X線管２０と検出器アレイ２３とは、X線管２０のX線の放射中心が検出器アレイ２３のz方向の中心に一致するように配置する。
【００２５】
移動部２１は、X線コントローラ２５からの制御信号CTL２５２に基づいて、X線管２０、コリメータ２２、および検出器アレイ２３を、これらの相対的な位置関係を保ったまま、z方向に沿って移動させる。
【００２６】
被検体６は、撮影テーブル４上に載置される。撮影テーブル４は、たとえばモーターによって位置変更が可能になっており、操作コンソール３からの制御信号CTL３０bを受けて、被検体６をX線照射空間２９内に搬入する。被検体６は、たとえば、頭部から脚部への体軸の方向がz方向に一致するように搬入される。
【００２７】
X線コントローラ２５は、操作コンソール３の後述する中央処理装置３０からの制御信号CTL３０１に基づいて、X線管２０が放射するX線の強度を制御するための制御信号CTL２５１、および、移動部２１を制御するための制御信号CTL２５２を出力する。
X線コントローラ２５は、制御信号CTL２５２によって、z方向に沿って移動部２１を所定ピッチ移動させる。
【００２８】
コリメータコントローラ２６は、操作コンソール３からの制御信号CTL３０２に応じて、コリメータ２２に出力する制御信号CTL２６１によってアパーチャ２２１の開度を制御する。
本実施の形態においては、コリメータ２２を通過したX線ビーム５が、検出器アレイ２３のX線検出面２３Sの全面に照射されるようにアパーチャ２２１の開度を調節する。
【００２９】
回転コントローラ２８は、操作コンソール３の中央処理装置３０からの制御信号CTL３０４に基づいて、回転部２７に対し制御信号を出力して回転軸AXまわりに回転させる。
本実施の形態においては、たとえば、回転軸AX上の所定の位置において回転部２７が1回転し、回転軸AXまわりの全ビューにおいて検出器アレイ２３がX線強度を検出することによって１回のスキャンが行なわれる。
【００３０】
図１に示すように、操作コンソール３は、中央処理装置３０と、入力装置３１と、表示装置３２と、記憶装置２３とを有する。
入力装置３１は、X線CT装置１を動かすための操作者（オペレータ）からの指令入力を受け付けて中央処理装置３０に送信する。
表示装置３２は、中央処理装置３０による計算の結果得られたCT画像を表示するためのものである。
記憶装置２３は、中央処理装置３０を介してX線CT装置１を動かすための各種プログラムやパラメータ、CT画像の画像データ等のデータを記憶する。
【００３１】
図２に示すように、中央処理装置３０は、制御部３４と、再構成部３６と、表示部３８と、再構成画像生成部４０とを有する。
再構成部３６が本発明における画像再構成手段の一実施態様であり、再構成画像生成部４０が再構成画像生成手段の一実施態様である。
中央処理装置３０は、１つのハードウェアのみを用いた構成であってもよいし、上述の各構成要素に対して適宜ハードウェアを割り当てたシステムとして構成してもよい。
【００３２】
制御部３４には、DAS２４、入力装置３１、再構成部３６、表示部３８、X線コントローラ２５、コリメータコントローラ２６、回転コントローラ２８、および撮影テーブル４が接続される。
再構成部３６は表示部３８および再構成画像生成部４０にさらに接続されており、また、表示部３８と再構成画像生成部４０とが接続されている。
また、記憶装置２３は、中央処理装置３０の各構成要素から適宜アクセスされる。
【００３３】
制御部３４は、撮影テーブル４に制御信号CTL３０bを出力して被検体６をX線照射空間２９内に搬入させる。
また、制御部３４は、CT画像としての被検体６の断層像の生成に用いる検出データを得るために、入力装置３１を介した操作者からの指令入力に基づいて、操作ガントリ２に制御信号CTL３０aを出力し、被検体６のスキャンを行なわせる。制御部３４からの制御信号CTL３０aには、前述の制御信号CTL３０１，３０２，３０３，３０４が含まれる。
制御部３４は、スキャンによってDAS２４が収集した検出データを受け取る。
【００３４】
再構成部３６は、制御部３４を介して入手した検出データに基づいて、被検体６の被検部位の断層像を再構成する。
制御部３４を介した走査ガントリ２による被検体６のスキャンの例、および、再構成部３６による被検体６の断層像生成については後ほど詳述する。
【００３５】
再構成画像生成部４０は、再構成部３６を介した制御部３４からの指令に応じて、再構成部３６によって再構成された断層像のデータに基づいて断面変換を行ない、再構成部３６によって再構成された断層像とは異なる方向の断層像を生成する。
表示部３８は、制御部３４からの指令信号に応じて、再構成部３６が生成した断層像および／または再構成画像生成部４０が生成した断層像を表示装置３２に表示させる。
【００３６】
断層像生成
以下では、走査ガントリ２による被検体６の被検部位のスキャン、および、このスキャンによって得られた検出データに基づいた断層像生成について詳細に述べる。
【００３７】
図３は本実施の形態における断層像の生成方法について述べるための図であり、（a）はz軸を含む平面においてX線ビームと断層像との関係を示しており、（b）は断層像の斜視図を示している。
また、図４は、図３（a）の意味を説明するための図である。したがって、以下ではまず図４について述べる。
【００３８】
前述のように、X線管２０および検出器アレイ２３は、回転軸AXまわりに回転する。図４においては、この回転軸AXをz軸とする。
図４においては、X線管２０の初期状態の位置を０°方向とし、そこからz軸まわりに１８０°回転したときのX線管２０の位置を１８０°方向としている。また、０°方向および１８０°方向を含む平面がyz平面となるようにy方向を規定している。図４において、紙面に直交する方向がx方向となる。
【００３９】
本実施の形態においては、制御部３４は、z軸方向に沿った所定の位置において１回スキャンを行ない、その後z軸に沿った方向に所定ピッチずれた位置において次のスキャンを行なう、いわゆるアキシャルスキャンを走査ガントリ２に実行させる。図４においては、２回スキャンを行なった例が示されている。
【００４０】
制御部３４は、X線コントローラ２５に制御信号CTL３０１を出力してX線管２０からX線を放射させる。また、制御部３４は、コリメータコントローラ２６に制御信号CTL３０２を出力して、全てのX線検出器列１a〜４a，１b〜４bにX線ビームが照射されるようにコリメータ２２のアパーチャ２２１の開度を調節する。本実施の形態においては、図４に示すように、yz平面におけるX線ビームの断面形状が扇形のいわゆるファンビームが検出器アレイ２３に照射されることになる。
この状態において、制御部３４は回転部２７を回転させてスキャンを行なわせる制御信号CTL３０４を回転コントローラ２８に出力し、検出器アレイ２３によって得られる検出データを収集させる制御信号CTL３０３をDAS２４に出力する。
これにより、z軸まわりの全ビューにおける検出データ、すなわちローデータがDAS２４によって収集される。ローデータは、X線による被検部位の投影データになっている。
【００４１】
再構成部３６は、ローデータに対して逆投影処理を施し、被検部位の断層像を再構成する。
ただし、ローデータをそのまま用いて逆投影したのでは再構成した断層像にぼけやノイズが発生する。このため、再構成部３６はローデータに対してFFT（Fast Fourier Transform：高速フーリエ変換）を施して周波数領域のデータにした後に、ぼけやノイズを低減するためのフィルタリング処理を行ない、IFFT（Inverse Fast Fourier Transform：高速逆フーリエ変換）によって元のローデータの形式に戻すというローデータの加工を行なう。
再構成部３６は、加工したローデータに対して逆投影処理を施し、断層像を再構成する。
【００４２】
再構成部３６は、X線管２０から放射され検出器列１a〜４a，１b〜４bのそれぞれに入射する各X線ビームが回転軸AXに交差する位置における、回転軸AXに直交する断層像を再構成する。
たとえば、図４に示す断層像I１４bは、１回目のスキャンにおいて、０°方向から出射し１８０°方向に位置する検出器列４bに入射するX線ビーム１０＿４b、および、１８０°方向から出射し０°方向に位置する検出器列４bに入射するX線ビーム１８＿４bが回転軸AXに直交する位置における被検部位の断層像を表わしている。同様に、断層像I１４aは、１回目のスキャンにおいて０°方向および１８０°方向から検出器列４aにそれぞれ入射するX線ビーム１０＿４aおよび１８＿４aが回転軸AXに直交する位置の被検部位の断層像を示している。
以上により、本実施の形態においては、検出器列１a〜４a，１b〜４bのそれぞれに対応した８枚の断層像が１回のスキャンによって得られることがわかる。
【００４３】
１回のスキャンが終了後、制御部３４はX線管２０、コリメータ２２、および検出器アレイ２３をz方向にピッチPt分移動させる制御信号CTL３０１をX線コントローラ２４に出力する。移動部２１のz方向への移動終了後に、制御部３４は回転コントローラ２８およびDAS２４にそれぞれ制御信号CTL３０４およびCTL３０３を出力し、２回目のスキャンを行なわせる。
以上のようにして、制御部３４は走査ガントリ２にアキシャルスキャンを実行させる。
また、図４における断層像I２４ｂおよびI２４aは、２回目のスキャンにおける検出器列４bおよび４aに対応する位置の断層像をそれぞれ表わしている。
ピッチPtは、スキャンとスキャンとの境界に位置する断層像同士の距離が所定範囲内に収まるように決める。具体的には、たとえば、断層像I１４aと断層像I２４bとの間の距離が所定範囲内に収まるようにする。
【００４４】
以上のように、回転軸AXを画像の中心として、検出器列１a〜４a，１b〜４bの各々に対応した断層像が再構成されて生成される。図３（a）は、z軸上のある断層像がいずれの検出器列に対応した断層像であるかを直感的に示すために、z軸と一致している回転軸AX上に検出器列１a〜４a，１b〜４bを描いて図４を描き直した図である。X線管２０およびx，y，zの各軸の位置関係は、図４の場合と同じとなる。
なお、本実施の形態においては、図３（a）に示すように、１スキャン目の検出器列４aから連続的に２スキャン目の検出器列４bが位置するようにピッチPtを決めている。これは、各スキャンにおける端部の断層像同士の距離を他の断層像同士の距離とほぼ同じにし、また、各スキャンにおいてz方向に可能な限り広いを撮影し、撮影時間を短縮するためである。
【００４５】
図３（a）において、X線管２０から放射されるX線ビームが、図３（a）に描かれている所定の検出器列を通過することは、回転軸AXを挟んだ位置にあるその検出器列にそのX線ビームが実際に入射することを意味する。たとえば、図３（a）に示す１回目のスキャンにおいて、回転軸AX上の検出器列４aを０°方向から通過するX線ビーム１０＿４aは、１８０°方向にある実際の検出器列４aに実際に入射する。
以上のことから、回転軸AX上に検出器列１a〜４a，１b〜４bを描くことによって、回転軸AXにおける断層像がどの検出器列に対応したものであるかを判断できることがわかる。
【００４６】
たとえば、図３（a）における検出器列４aの位置の断層像を生成することを考える。
簡単に述べると、本実施の形態においてはまず、１スキャン目および２スキャン目の各々の検出データに基づいて、図３（b）に示すように検出器列４aに対応する位置における断層像I１４aおよびI５bの画像データを再構成する。そして、得られた断層像I１４aおよびI５bの画像データを合成することによって、検出器列４aに対応する位置における最終的な断層像Ie４aの画像データを生成する。
なお、図３（b）は、断層像I１４a，Ie４a，I５bがそれぞれxy平面に平行な平面であることを示しており、断層像I１４a，Ie４a，I５b上の各々の画素（ピクセル）Pe４a，Pe，Pe５bは全て同じ位置の被検部位に対応するピクセルである。
以下に断層像Ie４aのデータ生成の詳細を示す。
【００４７】
まず、再構成部３６による１スキャン目の検出データを用いた断層像I１４aのデータの再構成を考える。このとき、たとえばyz平面において考えた場合に、図３（a）から明らかなように、検出器列４aが実際に検出したX線ビーム１０＿４aおよび１８＿４aが通過する被検部位は、回転軸AX上の点Ctのみである。
このため、再構成部３６はコーンBP（Back Projection）法によって検出器列４aに対応する断層像I１４aを再構成する。１スキャン目において、たとえば、図３（a），（b）中のピクセルPeに対応する位置の被検部位を通過するX線ビーム１８＿３aは、検出器列３aによって検出されている。再構成部３６は、断層像I１４aにおいてピクセルPeと同じ被検部位を表わすピクセルPe４aの画像データを、検出器列３aによって得られた検出データを用いて再構成する。
また、断層像I１４aにおいて、点Ctに相当するピクセルからピクセルPe４aまでのピクセルの画像データは、たとえば、１スキャン目の検出器列４aの検出データと検出器列３aの検出データとの補間によって得られたローデータに基づいて再構成する。
以上のようなコーンBP法によって、再構成部３６は断層像I１４aの画像データを再構成する。
【００４８】
再構成部３６は、１スキャン目の場合と同様にして、２スキャン目の検出データを用いて検出器列４aに対応する断層像を再構成する。そのために、２スキャン目の位置における検出器アレイ２３の１スキャン目側の最端部に、図３（a）に示すように仮想的な検出器列５bが存在するものと考える。再構成部３６は、この検出器列５bに対応する断層像I５bを再構成する。
検出器列５bの位置は１スキャン目の検出器列４aの位置と同じであるため、２スキャン目の検出データを用いて、１スキャン目の検出器列４aに対応する断層像を再構成することになる。
【００４９】
２スキャン目における最も外側のX線ビームはX線ビーム２０＿４bおよび２８＿４bであり、検出器列５bに対応する位置において斜線で示す領域Aeを通過するX線ビームは存在せず、検出データを入手することができない。このため、再構成部３６は、たとえば、検出器列４bの検出データを外挿して、断層像I５bのうちの領域Aeに対応するピクセルの画像データを再構成する。その他の領域については、検出器列１b〜４bの検出データおよびこれらの補間データを用いて画像データを再構成する。
以上のようにして、再構成部３６は２スキャン目の検出データを用いたコーンBP法によって断層像I５bのデータを再構成する。
【００５０】
断層像I１４aおよび断層像I５bの画像データを再構成して入手した後に、再構成部３６は各ピクセルについて各々断層像I１４aの画像データと断層像I１５bの画像データとを加算して合成する。たとえば、断層像I１４aの画像データと断層像I１５bの画像データとの平均を計算し、得られた画像データを合成後の画像データとする。再構成部３６は、合成によって生成した合成断層像データに基づいて、最終的に検出器列４aに対応する位置における合成断層像Ie４aを生成する。検出器列４aに対応する位置のCT画像を表示する指令が操作者から入力されたときには、制御部３４によって表示部３８を介して合成断層像Ie４aが表示装置３２に表示される。
【００５１】
以上のように、本実施の形態においては、z方向に沿って異なる位置における複数回のスキャンの検出データを利用して１つの断層像データを生成する。たとえば、本実施の形態において再構成した合成断層像Ie４aのピクセルPeは、１スキャン目の検出器列３aを通過するX線ビーム１８＿３aによるデータと、２スキャン目の検出器列４bを通過するX線ビーム２８＿４b２によるデータとを用いて再構成される。このため、ピクセルPeの画質は、X線ビーム１８＿３aによるデータのみを用いる従来のコーンBP法によって再構成されるピクセルPe４aの画質よりも高くなる。以上から明らかなように、本実施の形態によって得られる合成断層像の画質は、従来のコーンBP法によって得られる断層像の画質よりも高くなる。
【００５２】
また、本実施の形態においては、検出器アレイ２３によって得られるローデータを合成して再構成することにより合成断層像の画像データを生成するのではなく、複数スキャンの各々によって得られる断層像の画像データを加算することにより合成断層像の画像データを生成する。このため、合成断層像の生成においてFFT等の再構成処理が不用であり、計算時間を短縮することができる。また、容易に合成断層像を得ることができる。
なお、検出器アレイ２３にはファンビームが照射されるため、各スキャンにおいては、外側の検出器列に対応する断層像ほど再構成に利用可能な検出データは少なくなる。このため、検出器アレイ２３のz方向に沿った端部側の検出器列に対応する合成断層像ほど、利用可能な検出データが従来よりも多くなり、画質向上の効果は高い。
【００５３】
重み決定方法
ここまでは複数回のスキャンの各々によって得られた断層像データの平均を取って合成断層像データを生成していた。しかし、z方向における合成断層像の位置や合成断層像中の各ピクセルの位置に応じて、合成に用いる複数の断層像データの各々に重みを適切に設定して合成すれば、さらなる画質の向上が期待できる。以下では、各断層像データに付与する重みの決定方法について述べる。
【００５４】
たとえば、図３（a）において検出器列４aに対応する位置の断層像データを再構成することを考える。簡単化のため、図３（a）における点Ctの位置のピクセルおよびピクセルPeについて考える。なお、合成断層像Ie４aのピクセルPeは、１スキャン目の断層像I１４aのピクセルPe４aと、２スキャン目の断層像I５bのピクセルPe５bとを合成して生成される。
１スキャン目において、点Ctに対応する被検部位については、０°方向からのX線ビーム１０＿４aおよび１８０°方向からのX線ビーム１８＿４aの２つのX線ビームによる検出データが存在するため、点Ctに対応するピクセルを高精度に再構成することができる。
１スキャン目において、ピクセルPe４aに対応する被検部位については、０°方向から通過するX線ビームは存在せず、１８０°方向からのX線ビーム１８＿３aによる検出データのみが入手可能であるため、ピクセルPe４aを高精度に再構成することは困難である。
【００５５】
２スキャン目において、点Ctに対応する被検部位については、０°方向からも１８０°方向からも通過するX線ビームが存在しないため、点Ctを通過したX線による実際の検出データを用いて点Ctに対応するピクセルを再構成することはできない。
２スキャン目において、ピクセルPe５bに対応する被検部位については、０°方向から通過するX線ビームは存在せず、１８０°方向からのX線ビーム２８＿４b２による検出データのみが入手可能であるため、ピクセルPe５bを高精度に再構成することは困難である。
以上により、合成断層像の画質を向上させるためには、z軸からの距離に応じて、各スキャンによって得られる各々の断層像データに付与する重みを変化させることが好ましいことが分かる。
【００５６】
特に、仮想的な検出器列に対応した断層像においては、被検部位を実際に通過したX線の検出データを用いることができない領域が存在し、この領域に対しては重みを０とすべきである。たとえば、図３（a）における斜線で示す領域Aeにおいては、断層像I５bの重みは０とする。
以下、図５を参照して、仮想的な検出器列に対応した断層像において重みを０とする領域の大きさを実際に導出する。
【００５７】
図５は、重みを０とする領域を導出するための図であり、基本的には図３（a）と同じである。そのため、同一部分には同一符号を付し、詳細な記述は省略する。
図５に示すように、z軸上において、z軸に投影した１スキャン目のX線管２０の放射中心から生成すべき合成断層像までの距離をd１とする。また、z軸上において、z軸に投影した２スキャン目のX線管２０の放射中心から生成すべき合成断層像までの距離をd２とする。
さらに、X線管２０の放射中心からz軸までの距離をLとする。
距離d１，d２，Lは、走査ガントリ２の装置構成から予め幾何学的に求めることができる。これらの距離の値は、たとえば記憶装置２３に記憶させておく。
【００５８】
ここで、図５における三角形lmnおよび三角形nopを考える。三角形nopにおける辺npの長さはd１である。したがって、三角形lmnにおける辺mnの長さはd２−d１となる。三角形lmnと三角形nopとは相似であるため、各辺の長さの比から、三角形lmnにおける辺lmの長さはL（d２−d１）／d１と求められる。再構成部３６は、生成すべき合成断層像のz方向の位置に応じた距離d１，d２、および距離Lのデータを記憶装置２３から入手し、辺lmの大きさを計算する。辺lmの大きさが、仮想的な検出器列に対応する断層像において、重みを０とするz軸からの距離r０となる。
【００５９】
仮想的な検出器列に対応する断層像において重みを０とする領域においては、合成断層像データ生成に用いる他の検出器列に対応する断層像データに付与する重みは１とする。
合成断層像におけるz軸からの距離rが距離r０よりも大きい領域においては、各スキャンによって再構成する各断層像のデータの各々に付与する重みを、距離rに応じて適宜変化させる。
また、距離r０が距離d１およびd２の関数であることから、各スキャンによって再構成する各断層像のデータの各々に付与する重みは、再構成する各断層像のz方向における位置に応じても変化させることが、画質向上の点から好ましい。
【００６０】
以上のように、本実施の形態においては、複数のスキャンの各々によって再構成する断層像のz方向における位置、および各断層像におけるz軸からの距離に応じて、各断層像の画像データに付与する重みを変化させる。
図６は重みの変化の一例を示すグラフであり、（a）は各検出器列に対応する断層像ごとに、z軸上の点Ctからの距離に応じた重みの変化を示しており、（b）はz軸方向の距離に応じた重みの変化を示している。
【００６１】
図６（a）において、横軸はz軸上の点Ctからの距離rを表わしており、縦軸は重みの大きさを表わしている。
また、検出器列１a〜４aのそれぞれに対応する断層像に付与する重みw１a〜w４a、および、検出器列５bに対応する断層像I５bに付与する重みを一例として挙げてそれぞれプロットしている。
【００６２】
前述のように、距離r０までは重みw５bの大きさは０である。重みw５bの大きさは、距離rが距離r０から増加するにしたがって、たとえば０．５まで増加させる。
重みw４aの大きさは、距離r０までは１であり、距離rが距離r０から増加するにしたがって、たとえば０．５まで減少させる。
このように、z方向の同じ位置における断層像の各々のデータに付与する重みは、合計が１になるようにする。
また、たとえば、検出器アレイ２３のz方向における中心部近傍の検出器列１aに対応する断層像においては、全てのピクセルは、各ピクセルに対応する被検体を通過したX線ビームの実際の検出データを用いて再構成することができる。このため、図６（a）に示すように、検出器列１aに対応する断層像データに対しては、付与する重みw１aの大きさを全て１とする。
【００６３】
重みw２aおよび重みw３aに関しては、検出器アレイ２３における内側の検出器列２aに付与する重みw２aの方が外側の検出器列３aに付与する重みw３aよりも大きさが１となる範囲が大きくなるようにする。また、重みw２aの大きさが重みw３aの大きさ以上となる状態で、距離rが増加するにしたがって、１から減少するようにする。
【００６４】
図６（a）において距離rをr＝reとして一定にした状態における各重みw１a〜w４aの大きさをプロットしたグラフが図６（b）である。図６（b）において、横軸はz軸上におけるX線管２０の放射中心からの距離d１を各検出器列１a〜４aの位置として表わしており、縦軸は重みの大きさを表わしている。
図６（b）に示すように、検出器アレイ２３の外側の検出器列に向かうほど、重みの大きさは小さくなる傾向にある。これは、外側の検出器列に対応する断層像ほど、直接検出したX線強度データを用いて再構成できるピクセルが少なくなるためである。
【００６５】
なお、ある断層像のピクセルの各々に対する重みの大きさは、たとえば、図７に示すように回転軸AXを中心として同心円状に変化させる。これは、X線管２０が回転軸AXのまわりに３６０°回転することに起因して、xy平面におけるX線の照射状態は回転軸AXを中心に同心円状に変化するためである。
図７において、（a）は検出器列４aに対応する断層像I１４aにおける重みの変化を示しており、（b）は検出器列３aに対応する断層像I１３aにおける重みの変化を示している。
【００６６】
図７（a），（b）においては、断層像I１４a，I１３aの各ピクセルに付加する重みw４a，w３aの大きさの分布状態が、たとえば、１〜０．５の範囲で前述のように回転軸AXを中心として同心円状に変化している。
また、各重みの分布状態を示す同心円の大きさは、z方向における距離に応じて変化している。前述のように、検出器アレイ２３の内側の検出器列に対応した断層像ほど実際の検出データを用いて再構成可能な領域が増える。このため、大きい重みが占める領域は、断層像I１４aよりも断層像I１３aの方が大きくなる。たとえば、重みの大きさを１とする領域を規定する半径r0の大きさは、断層像I１３aの方が断層像I１４aよりも大きい。
なお、以上のように検出器列に対応して同心円状に規定した重みのデータのマップは、たとえば記憶装置２３に記憶させておく。再構成部３６は、合成断層像データの計算時に、生成すべき合成断層像の位置における重みのマップを、合成に用いる断層像ごとに記憶装置２３から入手する。再構成部３６は、入手した重みのマップに基づいて各断層像データを重み付けして加算して、合成断層像データを生成する。
【００６７】
図６（a）においては重みの大きさが１または０から非線形に変化する状態を例示したが、たとえば、図６（a）におけるプロットw４a２に示すように、重みの大きさは線形に変化させてもよい。
また、たとえば、距離rが距離r０までは１、距離r０よりも大きくなれば０．５というように、不連続に重みを変化させてもよい。
重みの決め方は、以上のような方法に限らず、たとえば、w５b＋w４a＝１となるように、合成に用いる複数のピクセルの重みの合計が１となるように任意に決めてよい。
【００６８】
以上のように、本実施の形態においては、合成に用いる各断層像データについて、z方向の距離およびz軸からの距離に応じて重みを適宜付与する。これは、各断層像の各ピクセルの再構成の精度に応じて、各ピクセルに付与する重みを変化させることを意味する。このため、たとえば、再構成精度の高いピクセルの画像データを十分に活用して合成断層像データを生成することが可能になり、その結果、高画質の合成断層像を得ることができる。
また、図７に示すように同心円状に重みを変化させれば、検出器列毎の重みのマップ生成が容易である。
【００６９】
断面変換
以上述べてきたように、再構成部３６はz方向の所定の検出器列に対応する位置において、xy平面に平行な合成断層像を生成する。
この合成断層像をz方向に沿って複数枚生成することによって、xy平面とは異なる断面の断層像を生成することができる。以下にこの断面変換の詳細について述べる。
【００７０】
図８（a）は、yz平面においてz方向に複数の合成断層像を生成した状態を示した図であり、図８（b）は複数の合成断層像を利用して生成したyz平面のMPR画像の一例である。
図８（a）においては、図３（a）および図５と同様に、各スキャンにおいて生成される合成断層像の位置を表わすために、検出器列をz軸上に重ねて描いている。
図８（a）に示すように、i番目のスキャンによって、再構成部３６は、各検出器列４b，３b，２b，１b，１a，２a，３a，４aに対応する合成断層像Ie１i〜Ie８iを生成する。これら８枚の合成断層像をまとめてスキャン画像Siとする。
また、再構成部３６は、i＋１番目のスキャンによって各検出器列４b，３b，２b，１b，１a，２a，３a，４aに対応する合成断層像Ie１i+1〜Ie８i+1を生成する。これら８枚の合成断層像をまとめてスキャン画像Si+1とする。
【００７１】
このように、再構成部３６は、１回のスキャン毎に８枚の合成断層像を生成する。
再構成画像生成部４０は、断面変換を行なわせる指令が操作者から制御部３４に入力されると、再構成部３６がスキャン毎に生成した合成断層像の画像データを入手する。そして、再構成画像生成部４０は、入手した複数の合成断層像のピクセルのうち、z軸を含むyz平面におけるピクセルを全てまとめて、表示部３８を介して表示装置３２に表示させる。これにより、図８（b）に示すような、yz平面におけるMPR画像が得られる。
比較のために、従来のコーンBP法を適用して再構成した断層像を用いて断面変換した結果得られるMPR画像を図１０に示す。
図８（b）と図１０とを比較してわかるように、本実施の形態によって得られるMPR画像は、各ピクセル間の変化がなめらかであり、従来のMPR画像において認められるような、スキャン画像間の境界において画像がずれることによる縞状のアーチファクトの発生が抑制されている。また、各々のスキャン画像の画質も向上する。
【００７２】
以上のように、本実施の形態においては、z方向の異なる位置における複数のスキャンの各々の検出データを利用して生成した合成断層像の画像データを用いてMPR画像を生成している。このため、特に、スキャンとスキャンとの境界における合成断層像データの連続性が良くなる。その結果、従来のMPR画像よりも滑らかで高画質のMPR画像を得ることができる。
また、本実施の形態における各合成断層像は、従来のコーンBP法によって得られる断層像よりも高画質であるため、MPR画像を構成する各スキャン画像の画質も向上する。
【００７３】
ここで、本実施の形態における断層像の生成方法について、図９のフローチャートを参照して簡単にまとめる。
まず、z軸方向の所定の位置において、被検部位の投影データを入手するためのスキャンを行なう（ステップST１）。
複数のX線検出器を２次元的に配列した検出器アレイ２３を回転軸AXまわりに回転させることによって、被検体の容積データが入手される。
【００７４】
１回のスキャン終了後には、入手したX線の検出データに基づいて、断層像データを再構成する（ステップST２）。
検出器列毎に得られる検出データにより、z方向に沿って配列される複数の断層像の画像データを１回のスキャンによって再構成する。
また、合成断層像データの生成に用いるために、仮想的な検出器列に対応する断層像のデータも、ステップST１のスキャンによって得られた検出データを用いて再構成しておく。
【００７５】
ステップST２において得られた複数の断層像を各々構成する画像データのそれぞれに対しては、大きさを適宜定めた重みを付与する（ステップST３）。たとえば、z軸に沿った方向距離およびz軸からの距離に応じて予め規定してある重みを、それぞれの画像データに付与する。
【００７６】
スキャンが１回終了する毎に、所望の範囲の画像が得られるだけの回数のスキャンが行なわれたか否かを判断する（ステップST４）。
スキャンの回数が不足していると判断した場合には、z方向にスキャン位置を所定ピッチずらせる（ステップST５）。
そして、ステップST１のスキャンからステップST３の断層像データへの重み付与までのステップを実行する。
ステップST４において、所望の範囲の画像が得られるだけの回数のスキャンが行なわれると判断されるまで、ステップST１からステップST５までが繰り返される。
【００７７】
たとえば、全てのスキャンが終了しz方向に配列した断層像の全ての画像データが得られた後に、z方向における同じ位置に存在する断層像と仮想的な断層像とを合成する（ステップST６）。
z方向における同じ位置に存在する各断層像の画像データを、それぞれの画像データに対して規定された重みを付与して加算して、合成断層像の画像データを計算する。
【００７８】
z方向における全ての合成断層像の画像データの入手後には、必要に応じて断面変換を行ない、MPR画像を生成する（ステップST７）。
合成断層像の各座標を揃えて配列させ、合成断層像の集合における所定の面内のピクセルをまとめて表示することにより、MPR画像が得られる。
以上により、断層像生成の手順は終了する。
【００７９】
なお、本発明は上記の実施の形態に限定されない。
たとえば、スキャンに用いる放射線は、X線に限らず、γ線等の他の放射線であってもよい。
また、実施の形態においては検出器列の数がj＝８の場合を示したが、X線検出器の列の数はこれに限らず、たとえば、j＝１６またはj＝３２でもよい。
上記の実施の形態においては移動部２１がz方向に移動してアキシャルスキャンを行なう形態を示したが、この形態に限らず、回転部２７と被検体６とがz方向において相対的に移動可能であればよい。
また、図９に示すフローチャートは本発明に係る断層画像生成方法の一例を示したものであり、図９のフローチャートにおける各ステップは、特許請求の範囲内において適宜移動可能である。たとえば、全てのスキャンが終了した後に、断層像を一括して再構成してもよい。また、断層像の合成は、全スキャン終了後ではなく合成に用いる複数の断層像を再構成した時点で行なうことが可能である。
さらに、z軸を基準とした距離ではなく、各断層像における被検部位の形状を認識して、被検部位の形状に応じて各ピクセルの重みを変化させてもよい。
MPR画像についても、yz平面に平行な画像だけでなく、xz平面に平行な画像を生成することもできる。
【００８０】
【発明の効果】
以上のように、本発明によれば、アキシャルスキャンによる断層像再構成において、被検体の軸方向における所定の画像生成位置の断層像の画質を向上可能な放射線計算断層画像装置を提供することができる。また、本発明によれば、この放射線計算断層画像装置において用いられ、画像生成位置の断層像の画質を向上可能な断層画像生成方法を提供することもできる。
【図面の簡単な説明】
【図１】本発明に係るX線CT装置の一実施の形態の装置構成を示す図である。
【図２】図１に示すX線CT装置の要部の構成を示す図である。
【図３】本発明の一実施の形態における断層像の生成方法について述べるための図であり、（a）はz軸を含む平面においてX線ビームと断層像との関係を示しており、（b）は断層像の斜視図を示している。
【図４】図３（a）の意味を説明するための図である。
【図５】断層像において重みを０とする領域を導出するための図である。
【図６】重みの変化の一例を示すグラフであり、（a）は各検出器列に対応する断層像ごとに、z軸上の距離に応じた重みの変化を示しており、（b）はz軸方向の距離に応じた重みの変化を示している。
【図７】断層像における重みの変化を示す図であり、（a），（b）はそれぞれ異なるる断層像における重みの分布状態を示している。
【図８】（a）は、yz平面においてz方向に複数の合成断層像を生成した状態を示した図であり、（b）は複数の合成断層像を利用して生成したyz平面のMPR画像の一例である。
【図９】本発明に係る断層画像生成方法の一実施の形態を示すフローチャートである。
【図１０】従来のMPR画像の一例を示す図である。
【符号の説明】
１…X線CT装置
２…走査ガントリ
３…操作コンソール
５，１０＿４a，１８＿３a，１８＿４a，２０＿４b，２８＿４b，２８＿４b２…X線ビーム
６…被検体
２０…X線管
２２…コリメータ
２３…検出器アレイ
２４…データ収集システム
３０…中央処理装置
３４…制御部
３６…再構成部
３８…表示部
４０…再構成画像生成部
１a〜４a，１b〜５b…検出器列
AX…回転軸
【図１】
0004056922.tif 000002
【図２】
0004056922.tif 000003
【図３】
0004056922.tif 000004
【図４】
0004056922.tif 000005
【図５】
0004056922.tif 000006
【図６】
0004056922.tif 000007
【図７】
0004056922.tif 000008
【図８】
0004056922.tif 000009
【図９】
0004056922.tif 000010
【図１０】
0004056922.tif 000011