Ｘ線ＣＴシステム
【課題】高周波域成分を強調する再構成関数を用いて画像再構成された断層像に対して局所領域毎に画像特徴量依存型の画像フィルタが適用され、画質が読影用に最適化された断層像を表示可能としつつ、高周波域成分が欠落していない、画像フィルタ適用前の断層像の表示や調整にも対応可能とする。
【解決手段】高周波域成分が強調された断層像に対して画像フィルタを適用する画像処理を、断層像の画像再構成時に行うのではなく、断層像の表示時に行うようにし、断層像保存時は、画像フィルタ適用前の断層像を保存する。
【解決手段】高周波域成分が強調された断層像に対して画像フィルタを適用する画像処理を、断層像の画像再構成時に行うのではなく、断層像の表示時に行うようにし、断層像保存時は、画像フィルタ適用前の断層像を保存する。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、Ｘ線ＣＴ（ｃｏｍｐｕｔｅｄ ｔｏｍｏｇｒａｐｈｙ）システム（ｓｙｓｔｅｍ）に関し、さらに詳しくは、Ｘ線ＣＴ撮影によって得られた被検体の断層像を、画質改善して表示するＸ線ＣＴシステムに関するものである。
【背景技術】
【０００２】
従来、Ｘ線ＣＴ撮影によって得られた被検体の断層像を表示するＸ線ＣＴシステムとして、次のような処理を行うシステムが知られている。
【０００３】
まず、走査ガントリ（ｇａｎｔｒｙ）により被検体を走査して収集されたＸ線投影データ（ｄａｔａ）に基づいて、被検体の断層像の高周波域成分を強調する再構成関数を用いて被検体の断層像を画像再構成する。このとき、画像再構成される断層像が表す被検体の部位に応じて、画像再構成に用いる再構成関数を切り換える。例えば、断層像が表す部位によっては読影で注目すべき組織の種類が異なる場合があるので、その部位が頭部、胸部、腹部等のいずれであるか等によって、軟部用標準再構成関数、骨部用再構成関数、高周波域成分強調用再構成関数等を使い分ける。断層像が表す部位の判定には、例えば、断層像全体の画像特徴量や操作者によって入力された情報に基づいて行う。なお、このような高周波域成分を強調する処理は断層像の画質向上のために行われるものであるが、画像再構成という画像生成上最も早い段階で断層像の高周波域成分を強調するのは、一旦、高周波域成分の情報が失われると、その成分の情報を復元することが困難になるためである。
【０００４】
次に、その断層像に対して画像フィルタ（ｆｉｌｔｅｒ）を適用し、その断層像を読影に適した画質に調整する。具体的には、断層像の各局所領域に対して、当該局所領域の画像特徴量に応じて画像の高周波域成分の強調度が変化する画像フィルタであって画像特徴量と強調度との対応関係がそれぞれ異なる複数の画像フィルタを、局所領域が表す解剖学的な部位の種類に応じて切り換えて適用する。例えば、断層像の局所領域を構成する画素のＣＴ値の平均値とその局所領域の高周波域成分の強調度との関係を規定するゲイン（ｇａｉｎ）曲線を予め複数種類用意しておき、断層像上で注目する局所領域を走査しながら順次設定し、その注目する局所領域が、予め決められた複数の解剖学的な部位の種類、例えば、肺野部、軟部、骨部のうちいずれに該当するかを判定し、上記複数のゲイン曲線の中からその判定された部位の種類に対応したゲイン曲線を選択し、注目する局所領域に対してその選択したゲイン曲線で定義された画像フィルタを適用するという画像処理を実行する。例えば、上記判定において、局所領域が軟部であると判定された場合には、高周波域成分の強調度が軟部組織に相当するＣＴ値の範囲において比較的小さくなるようなゲイン曲線が選択され、また、局所領域が肺野部であると判定された場合には、高周波域成分の強調度が肺野に相当するＣＴ値の範囲において比較的大きくなるようなゲイン曲線が選択される。
【０００５】
このようにして、局所領域毎に当該局所領域の解剖学的な部位の種類に対応したゲイン曲線の画像フィルタが適用され、画質が読影用に最適化された処理済断層像が得られる。
【０００６】
その後、その処理済断層像は表示処理部に送られ、モニタ（ｍｏｎｉｔｏｒ）画面に表示される。また、処理済断層像は、操作者の要求によりもしくは自動で保存処理部に送られ、その処理済断層像がハードディスク（ｈａｒｄ ｄｉｓｋ）等の記憶装置やＤＶＤ−ＲＯＭ等の記憶媒体に画像データとして保存される。
【０００７】
なお、上記の如く画像再構成して得られた直後の断層像は、上述の通り、高周波域成分が強調された断層像であるため、そのままでは読影に不適な画像である。したがって、上記Ｘ線ＣＴシステムにおいては、高周波域成分を強調する画像再構成と、その画像再構成によって得られた断層像の画質を最適化するための、断層像に対する画像フィルタの適用とは、常に同時に行われる。
【０００８】
このようなＸ線ＣＴシステムによれば、断層像が表す被検体の部位に応じて、また、その個々の局所領域の解剖学的な部位の種類に応じて、適当な再構成関数や画像フィルタが適用され、画質が読影用に最適化された断層像を表示したり、保存したりすることができるため、操作者が、断層像に適用する再構成関数および画像フィルタのゲイン曲線を種々切り換えながら、断層像を読影に適した画質に調整するという煩雑な作業が不要であり、診断効率を向上させることができるという利点がある。
【０００９】
なお、断層像に対して画像特徴量に応じて画像フィルタを切り換えながら適用することにより、断層像の画質を最適化する手法に関する技術は、例えば、特許文献１，２等に開示されている。
【００１０】
【特許文献１】特開平９−１６４１３４号公報
【特許文献２】特開平１０−２３４７２５号公報
【発明の開示】
【発明が解決しようとする課題】
【００１１】
上記Ｘ線ＣＴシステムは、上述の通り、被検体の断層像を、高周波域成分を強調して再構成する処理と、その断層像に画像フィルタを適用する処理とを行って、読影用に画質が最適化された断層像を得、その後、当該断層像を表示したり保存したりしている。
【００１２】
しかしながら、このようなＸ線ＣＴシステムでは、断層像に画像フィルタを適用する際に、場合によっては少なくともその一部が平滑化されるため、高周波域成分の一部の情報が失われ復元が困難になった断層像を保存することになり、後から画像フィルタ適用前の断層像を表示させたい場合や画像フィルタ適用前の断層像に別の画像処理を施して加工したい場合などに対応できないという問題がある。
【００１３】
本発明は、上記事情に鑑み、高周波域成分を強調する画像再構成によって得られた断層像に対して、局所領域毎に、当該領域の解剖学的な部位の種類に応じて、画像特徴量に依存した画像フィルタを切り換えながら適用することにより画質が読影用に最適化された断層像を表示可能としつつ、高周波域成分が欠落していない画像フィルタ適用前の断層像を表示させたい場合や当該断層像に別の画像処理を施して加工したい場合などに対応できるＸ線ＣＴシステムを提供することを目的とするものである。
【課題を解決するための手段】
【００１４】
第１の観点では、本発明は、被検体を走査して収集されたＸ線投影データに基づいて、前記被検体の断層像の高周波域成分を強調する再構成関数を用いて、前記被検体の断層像を画像再構成する画像再構成手段と、前記画像再構成された断層像を記憶装置または記憶媒体に保存する保存処理手段と、前記保存手処理段により保存された断層像を読み出す画像読出手段と、前記読み出された断層像の少なくとも一部における各局所領域に対して、該局所領域の画像特徴量に応じて画像の高周波域成分の強調度が変化する画像フィルタを適用する画像処理であって、前記画像特徴量と前記強調度との対応関係を表すゲイン曲線を、該局所領域が表す解剖学的な部位の種類に応じて切り換える画像処理を施すことにより処理済断層像を得る画像処理手段、および、該処理済断層像を表示する表示処理手段を有する画像表示手段とを備えたことを特徴とするＸ線ＣＴシステムを提供する。
【００１５】
第２の観点では、本発明は、被検体を走査して収集されたＸ線投影データに基づいて、前記被検体の断層像の高周波域成分を強調する再構成関数を用いて、前記被検体の断層像を画像再構成し、前記画像再構成された断層像を記憶装置または記憶媒体に保存し、前記画像保存手段により保存された断層像を読み出し、前記読み出された断層像のうち少なくとも一部における各局所領域に対して、該局所領域の画像特徴量に応じて画像の高周波域成分の強調度が変化する画像フィルタを適用する画像処理であって、前記画像特徴量と前記強調度との対応関係を表すゲイン曲線を、該局所領域が表す解剖学的な部位の種類に応じて切り換える画像処理を施すことにより処理済断層像を得、該処理済断層像を表示することを特徴とするＸ線ＣＴ画像表示方法を提供する。
【００１６】
第３の観点では、本発明は、前記保存処理手段が、前記画像再構成された断層像とともに、コンピュータを前記画像表示手段として機能させるためのプログラム（ｐｒｏｇｒａｍ）を該断層像と関連付けて保存するものであり、前記画像表示手段が、前記断層像とともに前記プログラムを読み出して該プログラムを実行するコンピュータであることを特徴とする上記第１の観点のＸ線ＣＴシステムを提供する。
【００１７】
ここで、高周波域成分とは、断層像を表す画像における高周波域に属する成分を意味するものであり、例えば、空間周波数が、１０〜１５（ｌｉｎｅ ｐａｉｒ／ｃｍ）程度の成分を考えることができる。
【００１８】
また、断層像を保存するとは、その断層像のデータをデータ処理のために暫定的に記憶するような形態ではなく、その断層像のデータを利用可能に残すことを目的として保存することを意味するものである。例えば、操作者の操作に応じて、断層像のデータのタイトル（ｔｉｔｌｅ）や保存場所を特定する情報を表示させたり、他の記憶装置や記憶媒体へコピー（ｃｏｐｙ）したりすることができるように保存することが考えられる。
【００１９】
また、局所領域とは、１つの断層像または複数の連続する断層像で構成される３次元像の中の部分的な領域を意味するものであり、例えば、１または互いに隣接する２次元または３次元の複数画素で構成される領域を考えることができる。
【００２０】
また、画像特徴量とは、画像上の特徴を表す物理量を意味するものであり、例えば、画像を構成する各画素のＣＴ値の平均値やＣＴ値のヒストグラム（ｈｉｓｔｏｇｒａｍ）等を考えることができる。
【００２１】
また、解剖学的な部位の種類としては、例えば、心臓、肝臓、脳等の軟部、肋骨、背骨、頭骨等の骨部、肺に対応する肺野、血管や器官に注入される造影剤等を考えることができる。
【００２２】
なお、画像再構成としては、従来の２次元画像再構成の他、３次元画像再構成と称される、Ｘ線ビームの厚み成分すなわちコーンビーム（ｃｏｒｎ ｂｅａｍ）Ｘ線のコーン角を考慮した画像再構成であり、フェルドカンプ・アルゴリズム（ｆｅｌｄｋａｍｐ ａｌｇｏｒｉｔｈｍ）で代表されるコーンビーム画像再構成アルゴリズムを用いた画像再構成を考えることができる。
【００２３】
また、画像フィルタとしては、例えば、局所領域の略中央に位置する注目画素の値を、その局所領域を構成する各画素のＣＴ値にそれぞれ所定の係数を掛けて積和演算した値に置換する空間フィルタ、ｓｏｂｅｌフィルタ等を考えることができる。
【発明の効果】
【００２４】
本発明のＸ線ＣＴシステムでは、高周波域成分を強調して画像再構成された断層像に対して、局所領域毎に、当該領域の解剖学的な部位の種類に応じて、当該領域の画像特徴量と高周波域成分の強調度との対応関係を表すゲイン曲線が異なる画像フィルタを切り換えて適用する画像処理手段を、画像再構成手段ではなく画像表示手段が具備し、断層像を保存する手段は、画像フィルタ適用前の断層像を保存するようにしている。すなわち、本発明のＸ線ＣＴシステムでは、高周波域成分が強調された断層像に対して画像特徴量依存型の画像フィルタを適用する画像処理を、断層像の画像再構成と切り離して、断層像を表示する際に行い、断層像の保存は、断層像の画像再構成と切り離さずに行うようにしている。
【００２５】
これにより、断層像を保存する際に画像フィルタ適用前の断層像を保存し、断層像を表示する際に断層像に画像フィルタを適用して画像フィルタ適用後の断層像を表示することができ、高周波域成分強調する画像再構成によって得られた断層像に対して、局所領域毎に、当該領域の解剖学的な部位の種類に応じて、画像特徴量に依存した画像フィルタを切り換えながら適用することにより画質が読影用に最適化された断層像を表示可能としつつ、高周波域成分が欠落していない画像フィルタ適用前の断層像を表示させたい場合や当該断層像に別の画像処理を施して加工したい場合などに対応できる。
【発明を実施するための最良の形態】
【００２６】
以下、図面を参照して本発明を実施するための最良の形態を説明する。
【００２７】
（第１の実施形態）
本発明の第１の実施形態について説明する。
【００２８】
図１は、本発明の第１の実施形態によるＸ線ＣＴシステム１の構成ブロック（ｂｌｏｃｋ）図である。図１に示すように、Ｘ線ＣＴシステム１は、Ｘ線ＣＴ装置２と画像解析用ワークステーション（ｗｏｒｋ ｓｔａｔｉｏｎ）６とを備え、Ｘ線ＣＴ装置２は、走査ガントリ３、撮影テーブル（ｔａｂｌｅ）４および操作コンソール（ｃｏｎｓｏｌｅ）５により構成される。
【００２９】
走査ガントリ３は、撮影対象である被検体を走査（以下、スキャン（ｓｃａｎ）ともいう）して被検体の各ビュー（ｖｉｅｗ）のＸ線検出器データを収集するものである。この走査ガントリ３は、後述の制御インターフェース（ｉｎｔｅｒｆａｃｅ）５１を通じて制御され、被検体のスキャンが行われる。走査ガントリ３は、被検体にＸ線を照射するＸ線照射装置３１、被検体を透過したＸ線の強度分布を検出するＸ線検出装置３２、および、Ｘ線検出装置３２で検出されたＸ線の強度分布に基づいてＸ線検出器データを収集するデータ収集部３３を有する。なお、走査ガントリ３と制御インターフェース５１との接続については図示を省略する。
【００３０】
撮影テーブル４は、被検体を載置し、被検体を走査ガントリ３のＸ線照射装置３１とＸ線検出装置３２との間の空間に搬入するものである。この撮影テーブル４は、後述の制御インターフェース５１を通じて制御され、被検体の搬入が行われる。なお、撮影テーブル４と制御インターフェース５１との接続については図示を省略する。
【００３１】
操作コンソール５は、制御インターフェース５１、データ収集バッファ（ｂｕｆｆｅｒ）５２、記憶装置５３、データ処理装置（画像再構成手段）５４、データ保存装置（保存処理手段）５５、表示装置５６および操作装置５７を有する。
【００３２】
制御インターフェース５１は、被検体の断層像が画像再構成されるように、操作者による操作装置５７の操作に応じて、走査ガントリ３、撮影テーブル４、データ処理装置５４等を制御するものである。
【００３３】
データ収集バッファ５２は、制御インターフェース５１からの制御により、走査ガントリ３のデータ収集部３３で収集されたＸ線検出器データを受け取って記憶装置５３に入力するものである。
【００３４】
記憶装置５３は、データ収集バッファ５２からＸ線検出器データを入力され記憶するものである。
【００３５】
データ処理装置５４は、走査ガントリ３により被検体を走査して得られた複数ビューのＸ線検出器データ（Ｘ線投影データ）に基づいて、画像再構成したときに得られる断層像の高周波域成分が強調されるような再構成関数を用いて、被検体の断層像を画像再構成するものである。具体的には、データ処理装置５４は、データ収集バッファ５２を通じて記憶装置５３に記憶されたＸ線検出器データを基に、前処理、ビームハードニング（ｂｅａｍｈａｒｄｅｎｉｎｇ）補正処理、ｚフィルタ重畳処理、再構成関数重畳処理、３次元逆投影処理および後処理を行って、被検体の断層像を画像再構成する。データ処理装置５４は、例えば、コンピュータ（ｃｏｍｐｕｔｅｒ）等によって構成される。
【００３６】
前処理は、Ｘ線検出器データを投影データに変換する処理であり、具体的には、Ｘ線検出器データに対して、オフセット（ｏｆｆｓｅｔ）補正、対数変換、Ｘ線線量補正、感度補正を行う処理である。
【００３７】
ビームハードニング補正処理は、ビームハードニング現象の影響で、収集したデータに現れた非線形性を補正する処理であり、本実施形態では、前処理で得られた投影データに対して、その非線形性を補正する処理となる。ビームハードニング補正処理は、例えば、数式（１）にしたがって行われる。
【００３８】
【数１】
【００３９】
ここで、Ｄ11(view,j,i)はビームハードニング補正された投影データ、Ｄ1(view,j,i)は前処理で得られた投影データ、Ｂ0(j,i)，Ｂ1(j,i)，・・・は係数、“●（ドット（ｄｏｔ））”は乗算演算子を表す。また、View＝１〜Ｎview（ビューの数），ｊ＝１〜Ｎrow，（検出器の列の数），ｉ＝１〜Ｎch（検出器のチャネルの数）である。
【００４０】
ｚフィルタ重畳処理は、断層像において再構成の中心部と周辺部とでスライス厚を略一様にするための補正処理である。一般に、断層像では、再構成中心部に比べ周辺部の方が、スライス厚が厚くなる。一方、Ｘ線検出装置３２は、被検体の搬送方向と垂直な方向に並んだＸ線検出器アレイが被検体の搬送方向（以下、ｚ方向ともいう）に列配された多列Ｘ線検出器である。そこで、各ビュー角度の投影データについて、再構成中心部付近に対応したチャネルのデータに対しては、ｚ方向のデータのうち断層像の再構成に寄与するデータの範囲が狭くなるような係数のフィルタを重畳し、再構成周辺部付近に対応したチャネルのデータに対しては、ｚ方向のデータのうち断層像の再構成に寄与するデータの範囲が広くなるような係数のフィルタを重畳する。ｚフィルタ重畳処理は、例えば、数式（２）で表すような列方向フィルタサイズが５列のフィルタを用いて、数式（３）にしたがって行われる。
【００４１】
【数２】
【００４２】
ここで、Ｄ12(view,j,i)は補正されたＸ線検出器データを表す。
【００４３】
再構成関数重畳処理は、ビームハードニング補正処理およびｚフィルタ重畳処理が施された投影データをフーリエ（ｆｏｕｒｉｅｒ）変換により周波数領域に変換し、その周波数領域に対して再構成関数を掛けた後、さらに逆フーリエ変換を行う処理である。使用する再構成関数は、断層像が表す被検体の部位に応じて切り換えられ。例えば、断層像が表す部位によっては読影で注目すべき組織の種類が異なる場合があるので、その部位が頭部、胸部、腹部等のいずれであるか等によって、軟部用標準再構成関数、骨部用再構成関数、高周波域成分強調用再構成関数等を使い分ける。なお、いずれの再構成関数も断層像の高周波域成分が強調されるような関数である。
【００４４】
図２は、再構成関数重畳処理で使用する再構成関数の一例である、軟部用標準再構成関数ｋ１、骨部用再構成関数ｋ２および高周波域成分強調用再構成関数ｋ３の各々のついて、空間周波数と高周波域成分強調度との関係を示した図である。
【００４５】
再構成関数重畳処理は、例えば、数式（４）にしたがって行われる。
【００４６】
【数３】
【００４７】
ここで、Ｄ12はzフィルタ重畳処理後の投影データ、Ｄ13は再構成関数重畳処理後の投影データ、Kernel（j）は再構成関数、“＊”は重畳（コンボリューション（ｃｏｎｖｏｌｕｔｉｏｎ））演算子を表す。
【００４８】
３次元逆投影処理は、再構成関数重畳処理が施された投影データに対して３次元逆投影を行い、逆投影データを求める処理である。
【００４９】
後処理は、逆投影データに対してＣＴ値変換等を行い、規格化された、被検体の断層像に変換する処理である。
【００５０】
表示装置５６は、データ処理装置５４から出力される被検体の断層像やその他の情報を表示するものである。表示装置５６は、例えば、液晶モニタ等で構成される。
【００５１】
データ保存装置５５は、データ処理装置５４から出力される被検体の断層像を、操作者の要求により、または、自動で、記憶装置５３あるいはＤＶＤ−ＲＯＭ等の記憶媒体５３ｍに画像データとして、操作者による処理が可能な形態で保存するものである。
【００５２】
操作装置５７は、ポインティングデバイス（ｐｏｉｎｔｉｎｇ ｄｅｖｉｃｅ）、キーボード（ｋｅｙｂｏａｒｄ）等で構成される。操作装置５７は、操作者によって操作され、各種の指示や情報等をデータ処理装置５４に入力する。操作者は表示装置５６、データ保存装置５５および操作装置５７を使用してインタラクティブ（ｉｎｔｅｒａｃｔｉｖｅ）にＸ線ＣＴ装置２を操作する。
【００５３】
画像解析用ワークステーション６は、操作コンソール５により保存された断層像を読み出して断層像の画質を最適化して表示するものである。画像解析用ワークステーション６は、データ読出部（画像読出手段）６１、画像処理部（画像処理手段）６２、表示処理部（表示処理手段）６３および記憶部６４を有する。なお、画像処理部６２および表示処理部６３は、本発明における画像表示手段として機能するものである。
【００５４】
データ読出部６１は、操作コンソール５のデータ保存装置５５により保存された断層像を、記憶装置５３あるいは記憶媒体５３ｍから読み出し、記憶部６４に記憶させるものである。
【００５５】
画像処理部６２は、データ読出部６１により読み出され記憶部６４に記憶された断層像（以下、読出断層像という）のうち少なくとも一部における各局所領域に対して、その局所領域の画像特徴量に応じて画像の高周波域成分の強調度が変化する画像フィルタを適用する画像処理であって、その画像特徴量と強調度との対応関係を表すゲイン曲線を、その局所領域が表す解剖学的な部位の種類に応じて切り換えながら画像フィルタを適用することにより処理済断層像を得るものである。なお、局所領域毎にその領域が表す部位に適した画像フィルタを適用する手法については、例えば、特願２００６−０２３５０４に詳しく記載されている。
【００５６】
画像処理部６２は、画像フィルタ制御部６２１、部位判定部６２２、ゲイン曲線選択部６２３および画像フィルタ適用部６２４を有する。
【００５７】
画像フィルタ制御部６２１は、読出断層像のうち少なくとも一部である、画像処理の対象領域上で所定の仮想枠を走査しながら、画像フィルタを適用する局所領域を注目領域として順次設定するものである。
【００５８】
部位判定部６２２は、画像フィルタ制御部６２１により設定された注目領域の画像特徴量、すなわち、注目領域の画像を構成する画素の画素値であるＣＴ値に基づいて、当該注目領域の画像が表す解剖学的な部位の種類が、予め設定された複数の種類のうちいずれに該当するかを判定するものである。本実施形態では、注目領域を、肺野部、軟部および骨部のうちいずれに該当するかを判定する。例えば、局所領域の画像を構成する画素のＣＴ値の平均値（以下、平均ＣＴ値という）が−２００以下であれば肺野部、−２００〜＋２００の範囲内であれば軟部（縦隔）、＋２００以上であれば骨部と判定する。
【００５９】
ゲイン曲線選択部６２３は、予め設定された複数のゲイン曲線の中から、部位判定部６２２により判定された部位の種類に対応したゲイン曲線を選択するものである。具体的には、注目領域の画像が表す解剖学的な部位の種類が軟部であると判定された場合には、高周波域成分の強調度が、軟組織に相当するＣＴ値の範囲において比較的小さくなるようなゲイン曲線を選択し、また、注目領域の画像が表す解剖学的な部位の種類が、肺野部あるいは骨部であると判定された場合には、高周波域成分の強調度が、その判定された肺野あるいは骨のいずれかに相当するＣＴ値の範囲において比較的大きくなるようなゲイン曲線を選択する。なお、これら複数のゲイン曲線を特定する情報は、上記の記憶部６３に記憶されている。
【００６０】
図３は、ゲイン曲線の一例を示した図である。図示のゲイン曲線では、肺野および骨に相当するＣＴ値の範囲では高周波域成分の強調度を比較的大きくし、軟組織に相当するＣＴ値の範囲では高周波域成分の強調度を比較的小さくするように規定されている。このようなゲイン曲線で定義された画像フィルタを画像に適用すると、肺野および骨の領域では高周波域成分が保持されるか、より強調され、縦隔または軟組織の領域では高周波域成分が抑制され平滑化された画像が得られることになる。
【００６１】
図４は、画像フィルタの一例を示した図である。図示の画像フィルタは、３×３画素の領域を局所領域として、この局所領域の画像を構成する各画素に対応したフィルタ係数を、εを用いた所定の係数とする画像フィルタを示している。この画像フィルタは、ε＜０とすれば、高周波域成分強調フィルタとなり、ε＞０とすれば、平滑化フィルタとなる。例えば、ε＝１／９とすれば、この画像フィルタは平滑化フィルタとなり、当該画像フィルタを適用した画像は、数式（５）にしたがって導出される。
【００６２】
【数４】
【００６３】
ここで、ｇ1(x,y)は画像フィルタ適用後の画像における座標(x,y)に位置する画素の値（ＣＴ値）、ｇ(x,y)は画像フィルタ適用前の画像における座標(x,y)に位置する画素の値（ＣＴ値）を表す。なお、画像フィルタは図４に示すものの他、例えば、図５に示すようなものであってもよい。図５に示す画像フィルタは、５×５画素の領域を局所領域として、この局所領域の画像を構成する各画素に対応したフィルタ係数を、ε1，ε2およびε3を用いた所定の係数とするものである。
【００６４】
画像フィルタ適用部６２４は、ゲイン曲線選択部６２２により選択されたゲイン曲線で定義された画像フィルタを注目領域の画像に適用するものである。
【００６５】
画像フィルタ制御部６２１は、読取断層像のうち画像処理の対象領域すべてについて画像フィルタが適用されて、画像処理が施された処理済断層像が得られるように、組織判定部６２２、ゲイン曲線選択部６２３および画像フィルタ適用部６２４を制御するものである。
【００６６】
表示処理部６３は、液晶モニタ等で構成される不図示の表示部を有し、上述のようにして得られた処理済断層像を表示部の画面に表示させるものである。
【００６７】
次に、第１の実施形態に係るＸ線ＣＴシステム１の動作について説明する。
【００６８】
図６は、第１の実施形態に係るＸ線ＣＴシステム１における処理の流れを説明するフローチャート（ｆｌｏｗ ｃｈａｒｔ）である。
【００６９】
操作者による操作コンソール５の操作装置５７の操作に応答して、撮影対象である被検体が載置された撮影テーブル４が移動し、被検体が走査ガントリ３のＸ線照射空間に搬入される。そして、いわゆるスキャンが行われる。例えば、撮影テーブル４を停止させた状態でＸ線照射装置およびＸ線検出装置を回転させることにより、コンベンショナルスキャン（ｃｏｎｖｅｎｔｉｏｎａｌ ｓｃａｎ）、別名、アキシャルスキャン（ａｘｉａｌ ｓｃａｎ）が行われる。あるいは、Ｘ線照射装置およびＸ線検出装置の回転に並行して、撮影テーブル４を被検体の体軸方向に連続的に移動させることにより、Ｘ線照射装置３１およびＸ線検出装置３２は、被検体に関して相対的に、被検体を包囲する螺旋状の軌道に沿って旋回することになり、いわゆるヘリカルスキャン（ｈｅｌｉｃａｌ ｓｃａｎ）が行われる。
【００７０】
ここで、スキャン１回転当たりに複数（例えば１０００程度）のビューのＸ線検出器データが、Ｘ線検出装置３２−データ収集部３３−データ収集バッファ５２の系列によって収集され、記憶装置５３に記憶される（ステップ（ｓｔｅｐ）Ｓ１）。
【００７１】
次に、データ処理装置５４は、記憶装置５３に記憶されたＸ線検出器データを基に、前処理、ビームハードニング補正処理、ｚフィルタ重畳処理、再構成関数重畳処理、３次元逆投影処理および後処理を行って、高周波域成分が強調された被検体の断層像を生成し、その画像データをデータ保存装置５５および表示装置５６に出力する（ステップＳ２）。
【００７２】
被検体の断層像が生成されると、表示装置５６はデータ処理装置５４から出力された画像データに基づいて被検体の断層像を表示し、データ保存装置５５は同画像データを記憶装置５３に保存する（ステップＳ３）。このとき、画像データは、被検体を特定する情報が付帯され、操作者が処理可能な態様で保存される。
【００７３】
画像解析用ワークステーション６は、データ読出部６１により操作コンソール５の記憶装置５３に保存された被検体の断層像の画像データを読み出し、記憶部６４に記憶させる（ステップＳ４）。
【００７４】
画像処理部６２の画像フィルタ制御部６２１は、データ読出部６１により読み出された読出断層像のうち画像処理の対象領域において、局所領域である注目領域を設定する（ステップＳ５）。次に、部位判定部６２２が、その注目領域における平均ＣＴ値に基づいて、その注目領域の画像が表す解剖学的な部位の種類が、「軟部」であるか否かを判定する（ステップＳ６）。ここで、注目領域の画像が表す部位が「軟部」であると判定された場合には、ゲイン曲線選択部６２３は、記憶部６４に記憶されている、複数のゲイン曲線を特定する情報の中から、「軟部用ゲイン曲線」の情報を読み込んで取得し、画像フィルタ適用部６２４が、この「軟部用ゲイン曲線」で定義された画像フィルタを注目領域の画像に適用する（ステップＳ７）。
【００７５】
ステップＳ６において、注目領域の画像が表す部位が「軟部」でないと判定された場合には、部位判定部６２２は注目領域の画像が表す部位が「肺野部」であるか否かを判定する（ステップＳ８）。ここで、注目領域の画像が表す部位が「肺野部」であると判定された場合には、ゲイン曲線選択部６２３が記憶部６４に記憶されている、複数のゲイン曲線を特定する情報の中から、「肺野部用ゲイン曲線」を読み込んで取得し、画像フィルタ適用部６２４が、この「肺野部用ゲイン曲線」で定義された画像フィルタを注目領域の画像に適用する（ステップＳ９）。
【００７６】
ステップＳ８において、注目領域の画像が表す部位が「肺野部」でないと判定された場合には、注目領域の画像が表す部位は「骨部」であるとみなし、ゲイン曲線選択部６２３が記憶部６４に記憶されている複数のゲイン曲線を特定する情報の中から、「骨部用ゲイン曲線」を読み込んで取得し、画像フィルタ適用部６２４が、この「骨部用ゲイン曲線」で定義された画像フィルタを注目領域の画像に適用する（ステップＳ１０）。
【００７７】
注目領域の画像に対して画像フィルタが適用されると、画像フィルタ制御部６２１が、読取断層像における画像処理の対象領域上で、次の注目領域として設定すべき局所領域が残っているか否かを判定する（ステップＳ１１）。次の注目領域として設定すべき局所領域が残っている場合にはステップ５に戻り、新たな注目領域を設定する。一方、次の注目領域として設定すべき局所領域が残っていない場合には、画像処理を終了する。
【００７８】
このようにして処理済断層像が得られると、表示処理部６３は、その処理済断層像を表示部の画面に表示する（ステップＳ１２）。そして、読影者は、画像解析用ワークステーション６の表示部に表示された処理済断層像を見て読影することになる。
【００７９】
（第２の実施形態）
次に、本発明の第２の実施形態について説明する。
【００８０】
本発明の第２の実施形態によるＸ線ＣＴシステムは、図１に示すような、第１の実施形態によるＸ線ＣＴシステムと略同様の構成であるが、部位判定部６２２の動作が、第１の実施形態とは異なるものである。部位判定部６２２は、第１の実施形態においては、注目領域の画像が表す部位を、肺野部、軟部および骨部の３種類のうちいずれに該当するか判定していたが、この第２の実施形態では、注目領域の画像が表す部位を、肺野部、軟部、骨部および造影剤の４種類のうちいずれに該当するか判定するものである。
【００８１】
Ｘ線ＣＴ撮影では、被検体の血管に造影剤を注入して撮影する場合があるが、この造影剤を表す画像に対しては、造影剤に適した画像フィルタが適用されたり、あるいは、画像フィルタが適用されないようにしたりすることが望ましい場合がある。ところが、造影剤に対応するＣＴ値は、骨に対応するＣＴ値と近接しており、ＣＴ値が約＋２００〜＋５００の範囲において区別が難しく、単純な平均ＣＴ値の閾値判定では両者を正確に判別することが難しい。
【００８２】
その一方で、骨と造影剤とでは、内部構造や物質の均一性に違いがあり、その違いがＣＴ値の空間的な分布の違いとして現れる。
【００８３】
図７（ａ）は、血管（大動脈）に注入された造影剤を表す画像におけるＣＴ値の空間的な変化を示した図であり、図７（ｂ）は、骨を表す画像におけるＣＴ値の空間的な変化を示した図である。例えば、大動脈が造影剤で造影された場合には、図７（ａ）に示すように、平均ＣＴ値が＋５００程度であるとすると、その大動脈の断面におけるＣＴ値のばらつきは、ＣＴ値の変化幅として２０〜３０程度である。このＣＴ値のばらつきは、造影剤の均一性やビームハードニング等による影響が考えられる。また、骨の場合には、図７（ｂ）に示すように、一般的に、表面が骨膜で覆われ、内部には緻密質、更に中心部には海綿質が存在しており、ＣＴ値は、緻密質では比較的高くなり、海綿質では比較的低くなる。
【００８４】
そこで、本実施形態では、部位判定部６２２は、注目領域における平均ＣＴ値の閾値判定だけでなく、上記のような骨と造影剤のＣＴ値の空間的な分布に関する特徴の違いを利用して、注目領域の画像が表す組織を肺野部、軟部、骨部および造影剤の４種類のうちいずれに該当するか判定する。以下、その判定手法の一例を説明する。
【００８５】
まず、読出断層像における各画素に対してＣＴ値の閾値判定を行い、所定値以上のＣＴ値を有する画素領域、すなわち、骨または造影剤に相当するＣＴ値を有する画素領域が抽出された画像Ｐ０を得る。そして、この画像Ｐ０において、２次元的に連続した領域（以下、２次元連続領域という）毎に番号ｚ（ｚ＝１，２，・・・）を付与するラベリング処理を行う。次に、画像Ｐ０に対して、各２次元連続領域Ｒｚを所定画素幅分だけ収縮させる収縮用論理フィルタを、すべての２次元連続領域がなくなるまで段階的に適用し、各段階における画像を、Ｐ１，Ｐ２，・・・として得る。そして、隣接する２画像の組合せ毎に２画像間で減算処理を行い、差分画像Ｐ01，Ｐ12，・・・を得る。すなわち、画像Ｐ０から画像Ｐ１を減算して差分画像Ｐ01を得、画像Ｐ１から画像Ｐ２を減算して差分画像Ｐ12を得るといった減算処理を繰り返す。このようにして得られた複数の差分画像Ｐ01，Ｐ12，・・・に基づいて、図８に示すように、各２次元連続領域Ｒｚ毎に、所定画素幅で層状になった殻状領域Ｒz1，Ｒz2，・・・を認識する。そして、各２次元連続領域Ｒｚ毎に、この層状になった各殻状領域Ｒzi（ｉ＝１，２，・・・）に対して、殻状領域を構成する画素のＣＴ値の平均値Ｍziおよび標準偏差Ｓziを求め、さらに、これらＣＴ値の平均値Ｍziおよび標準偏差Ｓziの標準偏差ＳＭz，ＳＳzを算出し、これらの値が所定の閾値Ｔｈ以上である場合、例えば、ＳＭz＋ＳＳz≧Ｔｈを満たす場合には、その２次元連続領域Ｒｚを骨領域と判別し、上記式を満たさない場合には、その２次元連続領域Ｒｚを造影剤領域と判別する。
【００８６】
このようにして、読出断層像における骨領域と造影剤領域とを区別して抽出し、その位置を把握しておく。後は、設定された注目領域における平均ＣＴ値に基づいて、骨部または造影剤のいずれかであると判定された場合に、その注目領域の位置に基づいて、骨部であるか造影剤であるかを判定するようにすればよい。
【００８７】
次に、第２の実施形態によるＸ線ＣＴシステム１の動作について説明する。
【００８８】
図９は、第２の実施形態によるＸ線ＣＴシステム１における処理の流れを示すフローチャートである。なお、ステップＳ２１〜Ｓ２４は、第１の実施形態におけるステップＳ１〜Ｓ４と同様の処理になるので、ここでは説明を省略する。
【００８９】
データ読出部６１により断層像が読み出されると、画像処理部６２の部位判定部６２２が、上述の判定手法により、読出断層像における骨領域と造影剤領域を区別して抽出し、それぞれの位置を予め把握する（ステップＳ２５）。次に、画像フィルタ制御部６２１は、その読出断層像のうち画像処理の対象領域において、局所領域である注目領域を設定する（ステップＳ２６）。部位判定部６２２は、その注目領域における平均ＣＴ値に基づいて、その注目領域の画像が表す解剖学的な部位の種類が、「軟部」であるか否かを判定する（ステップＳＳ２７）。ここで、注目領域の画像が表す部位が「軟部」であると判定された場合には、ゲイン曲線選択部６２３は、記憶部６４に記憶されている、複数のゲイン曲線を特定する情報の中から、「軟部用ゲイン曲線」の情報を読み込んで取得し、画像フィルタ適用部６２４が、この「軟部用ゲイン曲線」で定義された画像フィルタを注目領域の画像に適用する（ステップＳ２８）。
【００９０】
ステップＳ２７において、注目領域の画像が表す部位が「軟部」でないと判定された場合には、部位判定部６２２は注目領域の画像が表す部位が「肺野部」であるか否かを判定する（ステップＳ２９）。ここで、注目領域の画像が表す部位が「肺野部」であると判定された場合には、ゲイン曲線選択部６２３が記憶部６４に記憶されている、複数のゲイン曲線を特定する情報の中から、「肺野部用ゲイン曲線」を読み込んで取得し、画像フィルタ適用部６２４が、この「肺野部用ゲイン曲線」で定義された画像フィルタを注目領域の画像に適用する（ステップＳ３０）。
【００９１】
ステップＳ２９において、注目領域の画像が表す部位が「肺野部」でないと判定された場合には、注目領域の画像が表す部位は「骨部または造影剤」であるとみなし、さらに、その注目領域の位置に基づいて注目領域が先に求めた骨領域上にあるか否かを判定し（ステップＳ３１）、肯定される場合には、注目領域を「骨部」とみなして、ゲイン曲線選択部６２３が記憶部６４に記憶されている複数のゲイン曲線を特定する情報の中から、「骨部用ゲイン曲線」を読み込んで取得し、画像フィルタ適用部６２４が、この「骨部用ゲイン曲線」で定義された画像フィルタを注目領域の画像に適用する（ステップＳ３２）。一方、ステップＳ３１において否定される場合には、注目領域を「造影剤」とみなして、ゲイン曲線選択部６２３が記憶部６４に記憶されている複数のゲイン曲線を特定する情報の中から、「造影剤用ゲイン曲線」を読み込んで取得し、画像フィルタ適用部６２４が、この「造影剤用ゲイン曲線」で定義された画像フィルタを注目領域の画像に適用する（ステップＳ３３）。
【００９２】
注目領域の画像に対して画像フィルタが適用されると、画像フィルタ制御部６２１が、読取断層像における画像処理の対象領域上で、次の注目領域として設定すべき局所領域が残っているか否かを判定する（ステップＳ３４）。次の注目領域として設定すべき局所領域が残っている場合にはステップ２６に戻り、新たな注目領域を設定する。一方、次の注目領域として設定すべき局所領域が残っていない場合には、画像処理を終了する。
【００９３】
このようにして処理済断層像が得られると、表示処理部６３は、その処理済断層像を表示部の画面に表示する（ステップＳ３５）。
【００９４】
このように、上記第１および第２の実施形態によるＸ線ＣＴシステムでは、高周波域成分を強調して画像再構成された断層像に対して、局所領域毎に、当該領域の解剖学的な部位の種類に応じて、当該領域の画像特徴量と高周波域成分の強調度との対応関係を表すゲイン曲線が異なる画像フィルタを切り換えて適用する画像処理手段を、画像再構成手段（例えば、操作コンソール）ではなく画像表示手段（例えば、画像解析用ワークステーション）が具備し、断層像を保存する手段は、画像フィルタ適用前の断層像を保存するようにしている。すなわち、本発明のＸ線ＣＴシステムでは、高周波域成分が強調された断層像に対して画像特徴量依存型の画像フィルタを適用する画像処理を、断層像の画像再構成と切り離して、断層像を表示する際に行い、断層像の保存は、断層像の画像再構成と切り離さずに行うようにしている。
【００９５】
これにより、断層像を保存する際に画像フィルタ適用前の断層像を保存し、断層像を表示する際に断層像に画像フィルタを適用して画像フィルタ適用後の断層像を表示することができ、高周波域成分を強調する画像再構成によって得られた断層像に対して、局所領域毎に、当該領域の解剖学的な部位の種類に応じて、画像特徴量に依存した画像フィルタを切り換えながら適用することにより画質が読影用に最適化された断層像を表示可能としつつ、高周波域成分が欠落していない画像フィルタ適用前の断層像を表示させたい場合や当該断層像に別の画像処理を施して加工したい場合などに対応できる。
【００９６】
なお、上記の実施形態では、被検体の断層像を画像再構成する手段と、画像再構成によって得られた断層像に画像フィルタを適用して処理済断層像を得る手段とを、Ｘ線ＣＴ装置２の操作コンソール５と画像解析用ワークステーション６とに分離して設けているが、例えば、これら両手段を分離せず、操作コンソール５に１つにまとめて設けるようにしてもよい。
【００９７】
また、上記の実施形態では、画像解析用ワークステーション６が、断層像に画像フィルタを適用して処理済断層像を得る画像処理部６２、および、処理済断層像を表示する表示処理部６３を予め有しているが、例えば、図１０に示すように、操作コンソール５の画像保存部Ｘ線ＣＴ装置２が、上記のような画像処理部６２および表示処理部６３を有しない、画像解析用ワークステーション等のコンピュータに対して、断層像の画像データと、コンピュータをこのような画像処理部および表示処理部として機能させるための画像表示ツール（ｔｏｏｌ）であるプログラム（ｐｒｏｇｒａｍ）とを共に転送し、コンピュータ側でそのプログラムを実行するようにしてもよい。このようにすれば、断層像に画像フィルタを適用して処理済断層像を得る画像処理部およびその処理済断層像を表示する表示処理部とを有しない、画像解析用ワークステーション等のコンピュータであっても、同様の画像処理を行うことができる。
【００９８】
また、上記の実施形態では、断層像に適用する画像フィルタとして２次元の画像空間フィルタを用いているが、例えば、被検体の搬送方向すなわち被検体の体軸方向に連続して並ぶ複数の断層像における当該体軸方向の画素を考慮した、３次元の画像空間フィルタを用いるようにしてもよい。近年、Ｘ線検出器として、多列Ｘ線検出器やフラットパネルＸ線検出器に代表されるマトリクス構造の２次元Ｘ線検出器が用いられるようになり、被検体の体軸方向に連続した複数の断層像からなる３次元画像において、被検体の搬送方向のスライス厚を充分薄くすることができ、被検体の断層面方向であるｘ，ｙ軸および被検体の体軸方向であるｚ軸の各方向に対して画素の分解能が均一な３次元画像を生成することが可能になっている。換言すると、３次元画像の画素分解能の等方性、すなわち、アイソトロピック（ｉｓｏｔｒｏｐｉｃ）が実現されている。したがって、スライス（ｓｌｉｃｅ）厚が充分に薄い３次元画像については、３次元の画像空間フィルタを用いた方が画質の最適化には高い効果を奏する場合がある。特に、画像ノイズ（ｎｏｉｓｅ）改善という観点では、より多くの画素を用いた方が、よい結果が得られるので、２次元の画像空間フィルタよりも、より多くの画素を用いる３次元の画像空間フィルタを用いる方が、効果が大きい。
【００９９】
また、上記の実施形態では、断層像に適用する画像フィルタとして、ＣＴ値に依存した画像フィルタを用いているが、さらに、断層像のスライス厚に依存した画像フィルタを用いる場合も考えられる。
【０１００】
上記の実施形態は、本発明を実施するための最良の形態の一例であり、本発明はこれらの実施形態に限定されない。すなわち、本発明は、本発明の趣旨を逸脱しない限り、あらゆる変更・追加が可能である。
【０１０１】
また、コンピュータを本発明のＸ線ＣＴシステムあるいはその各手段として機能させるためのプログラムおよびそのプログラムが記録された記録媒体も、本発明の実施形態の一例である。
【図面の簡単な説明】
【０１０２】
【図１】発明を実施するための最良の形態例である医用画像表示システムのブロック図
【図２】各種再構成関数による空間周波数と高域強調度との関係を示す図
【図３】ＣＴ値と高周波域成分の強調または抑制の程度を規定する画像フィルタ曲線の一例を示す図
【図４】断層像に適用する画像フィルタの一例を示す図
【図５】断層像に適用する他の画像フィルタの一例を示す図
【図６】第１の実施形態によるＸ線ＣＴシステムにおける処理の流れを示すフローチャート
【図７】骨の画像および造影剤の画像におけるＣＴ値の空間的な変化を示す図
【図８】骨の画像および造影剤の画像において分割された層状の殻状領域を示す図
【図９】第２の実施形態によるＸ線ＣＴシステムにおける処理の流れを示すフローチャート
【図１０】Ｘ線ＣＴ装置から断層像と画像表示ツールとを共にコンピュータに転送する実施形態の概念を示す図
【符号の説明】
【０１０３】
１ Ｘ線ＣＴシステム
２ Ｘ線ＣＴ装置
３ 走査ガントリ
４ 撮影テーブル
５ 操作コンソール
６ 画像解析用ワークステーション
３１ Ｘ線照射装置
３２ Ｘ線検出装置
３３ データ収集部
５１ 制御インターフェース
５２ データ収集バッファ
５３ 記憶装置
５４ データ処理装置（画像再構成手段）
５５ データ保存装置（保存処理手段）
５６ 表示装置
５７ 操作装置
６１ データ読出部（画像読出手段）
６２ 画像処理部（画像処理手段）
６３ 表示処理部（表示処理手段）
６４ 記憶部
６２１ 画像フィルタ制御部
６２２ 部位判定部
６２３ ゲイン曲線選択部
６２４ 画像フィルタ適用部
【技術分野】
【０００１】
本発明は、Ｘ線ＣＴ（ｃｏｍｐｕｔｅｄ ｔｏｍｏｇｒａｐｈｙ）システム（ｓｙｓｔｅｍ）に関し、さらに詳しくは、Ｘ線ＣＴ撮影によって得られた被検体の断層像を、画質改善して表示するＸ線ＣＴシステムに関するものである。
【背景技術】
【０００２】
従来、Ｘ線ＣＴ撮影によって得られた被検体の断層像を表示するＸ線ＣＴシステムとして、次のような処理を行うシステムが知られている。
【０００３】
まず、走査ガントリ（ｇａｎｔｒｙ）により被検体を走査して収集されたＸ線投影データ（ｄａｔａ）に基づいて、被検体の断層像の高周波域成分を強調する再構成関数を用いて被検体の断層像を画像再構成する。このとき、画像再構成される断層像が表す被検体の部位に応じて、画像再構成に用いる再構成関数を切り換える。例えば、断層像が表す部位によっては読影で注目すべき組織の種類が異なる場合があるので、その部位が頭部、胸部、腹部等のいずれであるか等によって、軟部用標準再構成関数、骨部用再構成関数、高周波域成分強調用再構成関数等を使い分ける。断層像が表す部位の判定には、例えば、断層像全体の画像特徴量や操作者によって入力された情報に基づいて行う。なお、このような高周波域成分を強調する処理は断層像の画質向上のために行われるものであるが、画像再構成という画像生成上最も早い段階で断層像の高周波域成分を強調するのは、一旦、高周波域成分の情報が失われると、その成分の情報を復元することが困難になるためである。
【０００４】
次に、その断層像に対して画像フィルタ（ｆｉｌｔｅｒ）を適用し、その断層像を読影に適した画質に調整する。具体的には、断層像の各局所領域に対して、当該局所領域の画像特徴量に応じて画像の高周波域成分の強調度が変化する画像フィルタであって画像特徴量と強調度との対応関係がそれぞれ異なる複数の画像フィルタを、局所領域が表す解剖学的な部位の種類に応じて切り換えて適用する。例えば、断層像の局所領域を構成する画素のＣＴ値の平均値とその局所領域の高周波域成分の強調度との関係を規定するゲイン（ｇａｉｎ）曲線を予め複数種類用意しておき、断層像上で注目する局所領域を走査しながら順次設定し、その注目する局所領域が、予め決められた複数の解剖学的な部位の種類、例えば、肺野部、軟部、骨部のうちいずれに該当するかを判定し、上記複数のゲイン曲線の中からその判定された部位の種類に対応したゲイン曲線を選択し、注目する局所領域に対してその選択したゲイン曲線で定義された画像フィルタを適用するという画像処理を実行する。例えば、上記判定において、局所領域が軟部であると判定された場合には、高周波域成分の強調度が軟部組織に相当するＣＴ値の範囲において比較的小さくなるようなゲイン曲線が選択され、また、局所領域が肺野部であると判定された場合には、高周波域成分の強調度が肺野に相当するＣＴ値の範囲において比較的大きくなるようなゲイン曲線が選択される。
【０００５】
このようにして、局所領域毎に当該局所領域の解剖学的な部位の種類に対応したゲイン曲線の画像フィルタが適用され、画質が読影用に最適化された処理済断層像が得られる。
【０００６】
その後、その処理済断層像は表示処理部に送られ、モニタ（ｍｏｎｉｔｏｒ）画面に表示される。また、処理済断層像は、操作者の要求によりもしくは自動で保存処理部に送られ、その処理済断層像がハードディスク（ｈａｒｄ ｄｉｓｋ）等の記憶装置やＤＶＤ−ＲＯＭ等の記憶媒体に画像データとして保存される。
【０００７】
なお、上記の如く画像再構成して得られた直後の断層像は、上述の通り、高周波域成分が強調された断層像であるため、そのままでは読影に不適な画像である。したがって、上記Ｘ線ＣＴシステムにおいては、高周波域成分を強調する画像再構成と、その画像再構成によって得られた断層像の画質を最適化するための、断層像に対する画像フィルタの適用とは、常に同時に行われる。
【０００８】
このようなＸ線ＣＴシステムによれば、断層像が表す被検体の部位に応じて、また、その個々の局所領域の解剖学的な部位の種類に応じて、適当な再構成関数や画像フィルタが適用され、画質が読影用に最適化された断層像を表示したり、保存したりすることができるため、操作者が、断層像に適用する再構成関数および画像フィルタのゲイン曲線を種々切り換えながら、断層像を読影に適した画質に調整するという煩雑な作業が不要であり、診断効率を向上させることができるという利点がある。
【０００９】
なお、断層像に対して画像特徴量に応じて画像フィルタを切り換えながら適用することにより、断層像の画質を最適化する手法に関する技術は、例えば、特許文献１，２等に開示されている。
【００１０】
【特許文献１】特開平９−１６４１３４号公報
【特許文献２】特開平１０−２３４７２５号公報
【発明の開示】
【発明が解決しようとする課題】
【００１１】
上記Ｘ線ＣＴシステムは、上述の通り、被検体の断層像を、高周波域成分を強調して再構成する処理と、その断層像に画像フィルタを適用する処理とを行って、読影用に画質が最適化された断層像を得、その後、当該断層像を表示したり保存したりしている。
【００１２】
しかしながら、このようなＸ線ＣＴシステムでは、断層像に画像フィルタを適用する際に、場合によっては少なくともその一部が平滑化されるため、高周波域成分の一部の情報が失われ復元が困難になった断層像を保存することになり、後から画像フィルタ適用前の断層像を表示させたい場合や画像フィルタ適用前の断層像に別の画像処理を施して加工したい場合などに対応できないという問題がある。
【００１３】
本発明は、上記事情に鑑み、高周波域成分を強調する画像再構成によって得られた断層像に対して、局所領域毎に、当該領域の解剖学的な部位の種類に応じて、画像特徴量に依存した画像フィルタを切り換えながら適用することにより画質が読影用に最適化された断層像を表示可能としつつ、高周波域成分が欠落していない画像フィルタ適用前の断層像を表示させたい場合や当該断層像に別の画像処理を施して加工したい場合などに対応できるＸ線ＣＴシステムを提供することを目的とするものである。
【課題を解決するための手段】
【００１４】
第１の観点では、本発明は、被検体を走査して収集されたＸ線投影データに基づいて、前記被検体の断層像の高周波域成分を強調する再構成関数を用いて、前記被検体の断層像を画像再構成する画像再構成手段と、前記画像再構成された断層像を記憶装置または記憶媒体に保存する保存処理手段と、前記保存手処理段により保存された断層像を読み出す画像読出手段と、前記読み出された断層像の少なくとも一部における各局所領域に対して、該局所領域の画像特徴量に応じて画像の高周波域成分の強調度が変化する画像フィルタを適用する画像処理であって、前記画像特徴量と前記強調度との対応関係を表すゲイン曲線を、該局所領域が表す解剖学的な部位の種類に応じて切り換える画像処理を施すことにより処理済断層像を得る画像処理手段、および、該処理済断層像を表示する表示処理手段を有する画像表示手段とを備えたことを特徴とするＸ線ＣＴシステムを提供する。
【００１５】
第２の観点では、本発明は、被検体を走査して収集されたＸ線投影データに基づいて、前記被検体の断層像の高周波域成分を強調する再構成関数を用いて、前記被検体の断層像を画像再構成し、前記画像再構成された断層像を記憶装置または記憶媒体に保存し、前記画像保存手段により保存された断層像を読み出し、前記読み出された断層像のうち少なくとも一部における各局所領域に対して、該局所領域の画像特徴量に応じて画像の高周波域成分の強調度が変化する画像フィルタを適用する画像処理であって、前記画像特徴量と前記強調度との対応関係を表すゲイン曲線を、該局所領域が表す解剖学的な部位の種類に応じて切り換える画像処理を施すことにより処理済断層像を得、該処理済断層像を表示することを特徴とするＸ線ＣＴ画像表示方法を提供する。
【００１６】
第３の観点では、本発明は、前記保存処理手段が、前記画像再構成された断層像とともに、コンピュータを前記画像表示手段として機能させるためのプログラム（ｐｒｏｇｒａｍ）を該断層像と関連付けて保存するものであり、前記画像表示手段が、前記断層像とともに前記プログラムを読み出して該プログラムを実行するコンピュータであることを特徴とする上記第１の観点のＸ線ＣＴシステムを提供する。
【００１７】
ここで、高周波域成分とは、断層像を表す画像における高周波域に属する成分を意味するものであり、例えば、空間周波数が、１０〜１５（ｌｉｎｅ ｐａｉｒ／ｃｍ）程度の成分を考えることができる。
【００１８】
また、断層像を保存するとは、その断層像のデータをデータ処理のために暫定的に記憶するような形態ではなく、その断層像のデータを利用可能に残すことを目的として保存することを意味するものである。例えば、操作者の操作に応じて、断層像のデータのタイトル（ｔｉｔｌｅ）や保存場所を特定する情報を表示させたり、他の記憶装置や記憶媒体へコピー（ｃｏｐｙ）したりすることができるように保存することが考えられる。
【００１９】
また、局所領域とは、１つの断層像または複数の連続する断層像で構成される３次元像の中の部分的な領域を意味するものであり、例えば、１または互いに隣接する２次元または３次元の複数画素で構成される領域を考えることができる。
【００２０】
また、画像特徴量とは、画像上の特徴を表す物理量を意味するものであり、例えば、画像を構成する各画素のＣＴ値の平均値やＣＴ値のヒストグラム（ｈｉｓｔｏｇｒａｍ）等を考えることができる。
【００２１】
また、解剖学的な部位の種類としては、例えば、心臓、肝臓、脳等の軟部、肋骨、背骨、頭骨等の骨部、肺に対応する肺野、血管や器官に注入される造影剤等を考えることができる。
【００２２】
なお、画像再構成としては、従来の２次元画像再構成の他、３次元画像再構成と称される、Ｘ線ビームの厚み成分すなわちコーンビーム（ｃｏｒｎ ｂｅａｍ）Ｘ線のコーン角を考慮した画像再構成であり、フェルドカンプ・アルゴリズム（ｆｅｌｄｋａｍｐ ａｌｇｏｒｉｔｈｍ）で代表されるコーンビーム画像再構成アルゴリズムを用いた画像再構成を考えることができる。
【００２３】
また、画像フィルタとしては、例えば、局所領域の略中央に位置する注目画素の値を、その局所領域を構成する各画素のＣＴ値にそれぞれ所定の係数を掛けて積和演算した値に置換する空間フィルタ、ｓｏｂｅｌフィルタ等を考えることができる。
【発明の効果】
【００２４】
本発明のＸ線ＣＴシステムでは、高周波域成分を強調して画像再構成された断層像に対して、局所領域毎に、当該領域の解剖学的な部位の種類に応じて、当該領域の画像特徴量と高周波域成分の強調度との対応関係を表すゲイン曲線が異なる画像フィルタを切り換えて適用する画像処理手段を、画像再構成手段ではなく画像表示手段が具備し、断層像を保存する手段は、画像フィルタ適用前の断層像を保存するようにしている。すなわち、本発明のＸ線ＣＴシステムでは、高周波域成分が強調された断層像に対して画像特徴量依存型の画像フィルタを適用する画像処理を、断層像の画像再構成と切り離して、断層像を表示する際に行い、断層像の保存は、断層像の画像再構成と切り離さずに行うようにしている。
【００２５】
これにより、断層像を保存する際に画像フィルタ適用前の断層像を保存し、断層像を表示する際に断層像に画像フィルタを適用して画像フィルタ適用後の断層像を表示することができ、高周波域成分強調する画像再構成によって得られた断層像に対して、局所領域毎に、当該領域の解剖学的な部位の種類に応じて、画像特徴量に依存した画像フィルタを切り換えながら適用することにより画質が読影用に最適化された断層像を表示可能としつつ、高周波域成分が欠落していない画像フィルタ適用前の断層像を表示させたい場合や当該断層像に別の画像処理を施して加工したい場合などに対応できる。
【発明を実施するための最良の形態】
【００２６】
以下、図面を参照して本発明を実施するための最良の形態を説明する。
【００２７】
（第１の実施形態）
本発明の第１の実施形態について説明する。
【００２８】
図１は、本発明の第１の実施形態によるＸ線ＣＴシステム１の構成ブロック（ｂｌｏｃｋ）図である。図１に示すように、Ｘ線ＣＴシステム１は、Ｘ線ＣＴ装置２と画像解析用ワークステーション（ｗｏｒｋ ｓｔａｔｉｏｎ）６とを備え、Ｘ線ＣＴ装置２は、走査ガントリ３、撮影テーブル（ｔａｂｌｅ）４および操作コンソール（ｃｏｎｓｏｌｅ）５により構成される。
【００２９】
走査ガントリ３は、撮影対象である被検体を走査（以下、スキャン（ｓｃａｎ）ともいう）して被検体の各ビュー（ｖｉｅｗ）のＸ線検出器データを収集するものである。この走査ガントリ３は、後述の制御インターフェース（ｉｎｔｅｒｆａｃｅ）５１を通じて制御され、被検体のスキャンが行われる。走査ガントリ３は、被検体にＸ線を照射するＸ線照射装置３１、被検体を透過したＸ線の強度分布を検出するＸ線検出装置３２、および、Ｘ線検出装置３２で検出されたＸ線の強度分布に基づいてＸ線検出器データを収集するデータ収集部３３を有する。なお、走査ガントリ３と制御インターフェース５１との接続については図示を省略する。
【００３０】
撮影テーブル４は、被検体を載置し、被検体を走査ガントリ３のＸ線照射装置３１とＸ線検出装置３２との間の空間に搬入するものである。この撮影テーブル４は、後述の制御インターフェース５１を通じて制御され、被検体の搬入が行われる。なお、撮影テーブル４と制御インターフェース５１との接続については図示を省略する。
【００３１】
操作コンソール５は、制御インターフェース５１、データ収集バッファ（ｂｕｆｆｅｒ）５２、記憶装置５３、データ処理装置（画像再構成手段）５４、データ保存装置（保存処理手段）５５、表示装置５６および操作装置５７を有する。
【００３２】
制御インターフェース５１は、被検体の断層像が画像再構成されるように、操作者による操作装置５７の操作に応じて、走査ガントリ３、撮影テーブル４、データ処理装置５４等を制御するものである。
【００３３】
データ収集バッファ５２は、制御インターフェース５１からの制御により、走査ガントリ３のデータ収集部３３で収集されたＸ線検出器データを受け取って記憶装置５３に入力するものである。
【００３４】
記憶装置５３は、データ収集バッファ５２からＸ線検出器データを入力され記憶するものである。
【００３５】
データ処理装置５４は、走査ガントリ３により被検体を走査して得られた複数ビューのＸ線検出器データ（Ｘ線投影データ）に基づいて、画像再構成したときに得られる断層像の高周波域成分が強調されるような再構成関数を用いて、被検体の断層像を画像再構成するものである。具体的には、データ処理装置５４は、データ収集バッファ５２を通じて記憶装置５３に記憶されたＸ線検出器データを基に、前処理、ビームハードニング（ｂｅａｍｈａｒｄｅｎｉｎｇ）補正処理、ｚフィルタ重畳処理、再構成関数重畳処理、３次元逆投影処理および後処理を行って、被検体の断層像を画像再構成する。データ処理装置５４は、例えば、コンピュータ（ｃｏｍｐｕｔｅｒ）等によって構成される。
【００３６】
前処理は、Ｘ線検出器データを投影データに変換する処理であり、具体的には、Ｘ線検出器データに対して、オフセット（ｏｆｆｓｅｔ）補正、対数変換、Ｘ線線量補正、感度補正を行う処理である。
【００３７】
ビームハードニング補正処理は、ビームハードニング現象の影響で、収集したデータに現れた非線形性を補正する処理であり、本実施形態では、前処理で得られた投影データに対して、その非線形性を補正する処理となる。ビームハードニング補正処理は、例えば、数式（１）にしたがって行われる。
【００３８】
【数１】
【００３９】
ここで、Ｄ11(view,j,i)はビームハードニング補正された投影データ、Ｄ1(view,j,i)は前処理で得られた投影データ、Ｂ0(j,i)，Ｂ1(j,i)，・・・は係数、“●（ドット（ｄｏｔ））”は乗算演算子を表す。また、View＝１〜Ｎview（ビューの数），ｊ＝１〜Ｎrow，（検出器の列の数），ｉ＝１〜Ｎch（検出器のチャネルの数）である。
【００４０】
ｚフィルタ重畳処理は、断層像において再構成の中心部と周辺部とでスライス厚を略一様にするための補正処理である。一般に、断層像では、再構成中心部に比べ周辺部の方が、スライス厚が厚くなる。一方、Ｘ線検出装置３２は、被検体の搬送方向と垂直な方向に並んだＸ線検出器アレイが被検体の搬送方向（以下、ｚ方向ともいう）に列配された多列Ｘ線検出器である。そこで、各ビュー角度の投影データについて、再構成中心部付近に対応したチャネルのデータに対しては、ｚ方向のデータのうち断層像の再構成に寄与するデータの範囲が狭くなるような係数のフィルタを重畳し、再構成周辺部付近に対応したチャネルのデータに対しては、ｚ方向のデータのうち断層像の再構成に寄与するデータの範囲が広くなるような係数のフィルタを重畳する。ｚフィルタ重畳処理は、例えば、数式（２）で表すような列方向フィルタサイズが５列のフィルタを用いて、数式（３）にしたがって行われる。
【００４１】
【数２】
【００４２】
ここで、Ｄ12(view,j,i)は補正されたＸ線検出器データを表す。
【００４３】
再構成関数重畳処理は、ビームハードニング補正処理およびｚフィルタ重畳処理が施された投影データをフーリエ（ｆｏｕｒｉｅｒ）変換により周波数領域に変換し、その周波数領域に対して再構成関数を掛けた後、さらに逆フーリエ変換を行う処理である。使用する再構成関数は、断層像が表す被検体の部位に応じて切り換えられ。例えば、断層像が表す部位によっては読影で注目すべき組織の種類が異なる場合があるので、その部位が頭部、胸部、腹部等のいずれであるか等によって、軟部用標準再構成関数、骨部用再構成関数、高周波域成分強調用再構成関数等を使い分ける。なお、いずれの再構成関数も断層像の高周波域成分が強調されるような関数である。
【００４４】
図２は、再構成関数重畳処理で使用する再構成関数の一例である、軟部用標準再構成関数ｋ１、骨部用再構成関数ｋ２および高周波域成分強調用再構成関数ｋ３の各々のついて、空間周波数と高周波域成分強調度との関係を示した図である。
【００４５】
再構成関数重畳処理は、例えば、数式（４）にしたがって行われる。
【００４６】
【数３】
【００４７】
ここで、Ｄ12はzフィルタ重畳処理後の投影データ、Ｄ13は再構成関数重畳処理後の投影データ、Kernel（j）は再構成関数、“＊”は重畳（コンボリューション（ｃｏｎｖｏｌｕｔｉｏｎ））演算子を表す。
【００４８】
３次元逆投影処理は、再構成関数重畳処理が施された投影データに対して３次元逆投影を行い、逆投影データを求める処理である。
【００４９】
後処理は、逆投影データに対してＣＴ値変換等を行い、規格化された、被検体の断層像に変換する処理である。
【００５０】
表示装置５６は、データ処理装置５４から出力される被検体の断層像やその他の情報を表示するものである。表示装置５６は、例えば、液晶モニタ等で構成される。
【００５１】
データ保存装置５５は、データ処理装置５４から出力される被検体の断層像を、操作者の要求により、または、自動で、記憶装置５３あるいはＤＶＤ−ＲＯＭ等の記憶媒体５３ｍに画像データとして、操作者による処理が可能な形態で保存するものである。
【００５２】
操作装置５７は、ポインティングデバイス（ｐｏｉｎｔｉｎｇ ｄｅｖｉｃｅ）、キーボード（ｋｅｙｂｏａｒｄ）等で構成される。操作装置５７は、操作者によって操作され、各種の指示や情報等をデータ処理装置５４に入力する。操作者は表示装置５６、データ保存装置５５および操作装置５７を使用してインタラクティブ（ｉｎｔｅｒａｃｔｉｖｅ）にＸ線ＣＴ装置２を操作する。
【００５３】
画像解析用ワークステーション６は、操作コンソール５により保存された断層像を読み出して断層像の画質を最適化して表示するものである。画像解析用ワークステーション６は、データ読出部（画像読出手段）６１、画像処理部（画像処理手段）６２、表示処理部（表示処理手段）６３および記憶部６４を有する。なお、画像処理部６２および表示処理部６３は、本発明における画像表示手段として機能するものである。
【００５４】
データ読出部６１は、操作コンソール５のデータ保存装置５５により保存された断層像を、記憶装置５３あるいは記憶媒体５３ｍから読み出し、記憶部６４に記憶させるものである。
【００５５】
画像処理部６２は、データ読出部６１により読み出され記憶部６４に記憶された断層像（以下、読出断層像という）のうち少なくとも一部における各局所領域に対して、その局所領域の画像特徴量に応じて画像の高周波域成分の強調度が変化する画像フィルタを適用する画像処理であって、その画像特徴量と強調度との対応関係を表すゲイン曲線を、その局所領域が表す解剖学的な部位の種類に応じて切り換えながら画像フィルタを適用することにより処理済断層像を得るものである。なお、局所領域毎にその領域が表す部位に適した画像フィルタを適用する手法については、例えば、特願２００６−０２３５０４に詳しく記載されている。
【００５６】
画像処理部６２は、画像フィルタ制御部６２１、部位判定部６２２、ゲイン曲線選択部６２３および画像フィルタ適用部６２４を有する。
【００５７】
画像フィルタ制御部６２１は、読出断層像のうち少なくとも一部である、画像処理の対象領域上で所定の仮想枠を走査しながら、画像フィルタを適用する局所領域を注目領域として順次設定するものである。
【００５８】
部位判定部６２２は、画像フィルタ制御部６２１により設定された注目領域の画像特徴量、すなわち、注目領域の画像を構成する画素の画素値であるＣＴ値に基づいて、当該注目領域の画像が表す解剖学的な部位の種類が、予め設定された複数の種類のうちいずれに該当するかを判定するものである。本実施形態では、注目領域を、肺野部、軟部および骨部のうちいずれに該当するかを判定する。例えば、局所領域の画像を構成する画素のＣＴ値の平均値（以下、平均ＣＴ値という）が−２００以下であれば肺野部、−２００〜＋２００の範囲内であれば軟部（縦隔）、＋２００以上であれば骨部と判定する。
【００５９】
ゲイン曲線選択部６２３は、予め設定された複数のゲイン曲線の中から、部位判定部６２２により判定された部位の種類に対応したゲイン曲線を選択するものである。具体的には、注目領域の画像が表す解剖学的な部位の種類が軟部であると判定された場合には、高周波域成分の強調度が、軟組織に相当するＣＴ値の範囲において比較的小さくなるようなゲイン曲線を選択し、また、注目領域の画像が表す解剖学的な部位の種類が、肺野部あるいは骨部であると判定された場合には、高周波域成分の強調度が、その判定された肺野あるいは骨のいずれかに相当するＣＴ値の範囲において比較的大きくなるようなゲイン曲線を選択する。なお、これら複数のゲイン曲線を特定する情報は、上記の記憶部６３に記憶されている。
【００６０】
図３は、ゲイン曲線の一例を示した図である。図示のゲイン曲線では、肺野および骨に相当するＣＴ値の範囲では高周波域成分の強調度を比較的大きくし、軟組織に相当するＣＴ値の範囲では高周波域成分の強調度を比較的小さくするように規定されている。このようなゲイン曲線で定義された画像フィルタを画像に適用すると、肺野および骨の領域では高周波域成分が保持されるか、より強調され、縦隔または軟組織の領域では高周波域成分が抑制され平滑化された画像が得られることになる。
【００６１】
図４は、画像フィルタの一例を示した図である。図示の画像フィルタは、３×３画素の領域を局所領域として、この局所領域の画像を構成する各画素に対応したフィルタ係数を、εを用いた所定の係数とする画像フィルタを示している。この画像フィルタは、ε＜０とすれば、高周波域成分強調フィルタとなり、ε＞０とすれば、平滑化フィルタとなる。例えば、ε＝１／９とすれば、この画像フィルタは平滑化フィルタとなり、当該画像フィルタを適用した画像は、数式（５）にしたがって導出される。
【００６２】
【数４】
【００６３】
ここで、ｇ1(x,y)は画像フィルタ適用後の画像における座標(x,y)に位置する画素の値（ＣＴ値）、ｇ(x,y)は画像フィルタ適用前の画像における座標(x,y)に位置する画素の値（ＣＴ値）を表す。なお、画像フィルタは図４に示すものの他、例えば、図５に示すようなものであってもよい。図５に示す画像フィルタは、５×５画素の領域を局所領域として、この局所領域の画像を構成する各画素に対応したフィルタ係数を、ε1，ε2およびε3を用いた所定の係数とするものである。
【００６４】
画像フィルタ適用部６２４は、ゲイン曲線選択部６２２により選択されたゲイン曲線で定義された画像フィルタを注目領域の画像に適用するものである。
【００６５】
画像フィルタ制御部６２１は、読取断層像のうち画像処理の対象領域すべてについて画像フィルタが適用されて、画像処理が施された処理済断層像が得られるように、組織判定部６２２、ゲイン曲線選択部６２３および画像フィルタ適用部６２４を制御するものである。
【００６６】
表示処理部６３は、液晶モニタ等で構成される不図示の表示部を有し、上述のようにして得られた処理済断層像を表示部の画面に表示させるものである。
【００６７】
次に、第１の実施形態に係るＸ線ＣＴシステム１の動作について説明する。
【００６８】
図６は、第１の実施形態に係るＸ線ＣＴシステム１における処理の流れを説明するフローチャート（ｆｌｏｗ ｃｈａｒｔ）である。
【００６９】
操作者による操作コンソール５の操作装置５７の操作に応答して、撮影対象である被検体が載置された撮影テーブル４が移動し、被検体が走査ガントリ３のＸ線照射空間に搬入される。そして、いわゆるスキャンが行われる。例えば、撮影テーブル４を停止させた状態でＸ線照射装置およびＸ線検出装置を回転させることにより、コンベンショナルスキャン（ｃｏｎｖｅｎｔｉｏｎａｌ ｓｃａｎ）、別名、アキシャルスキャン（ａｘｉａｌ ｓｃａｎ）が行われる。あるいは、Ｘ線照射装置およびＸ線検出装置の回転に並行して、撮影テーブル４を被検体の体軸方向に連続的に移動させることにより、Ｘ線照射装置３１およびＸ線検出装置３２は、被検体に関して相対的に、被検体を包囲する螺旋状の軌道に沿って旋回することになり、いわゆるヘリカルスキャン（ｈｅｌｉｃａｌ ｓｃａｎ）が行われる。
【００７０】
ここで、スキャン１回転当たりに複数（例えば１０００程度）のビューのＸ線検出器データが、Ｘ線検出装置３２−データ収集部３３−データ収集バッファ５２の系列によって収集され、記憶装置５３に記憶される（ステップ（ｓｔｅｐ）Ｓ１）。
【００７１】
次に、データ処理装置５４は、記憶装置５３に記憶されたＸ線検出器データを基に、前処理、ビームハードニング補正処理、ｚフィルタ重畳処理、再構成関数重畳処理、３次元逆投影処理および後処理を行って、高周波域成分が強調された被検体の断層像を生成し、その画像データをデータ保存装置５５および表示装置５６に出力する（ステップＳ２）。
【００７２】
被検体の断層像が生成されると、表示装置５６はデータ処理装置５４から出力された画像データに基づいて被検体の断層像を表示し、データ保存装置５５は同画像データを記憶装置５３に保存する（ステップＳ３）。このとき、画像データは、被検体を特定する情報が付帯され、操作者が処理可能な態様で保存される。
【００７３】
画像解析用ワークステーション６は、データ読出部６１により操作コンソール５の記憶装置５３に保存された被検体の断層像の画像データを読み出し、記憶部６４に記憶させる（ステップＳ４）。
【００７４】
画像処理部６２の画像フィルタ制御部６２１は、データ読出部６１により読み出された読出断層像のうち画像処理の対象領域において、局所領域である注目領域を設定する（ステップＳ５）。次に、部位判定部６２２が、その注目領域における平均ＣＴ値に基づいて、その注目領域の画像が表す解剖学的な部位の種類が、「軟部」であるか否かを判定する（ステップＳ６）。ここで、注目領域の画像が表す部位が「軟部」であると判定された場合には、ゲイン曲線選択部６２３は、記憶部６４に記憶されている、複数のゲイン曲線を特定する情報の中から、「軟部用ゲイン曲線」の情報を読み込んで取得し、画像フィルタ適用部６２４が、この「軟部用ゲイン曲線」で定義された画像フィルタを注目領域の画像に適用する（ステップＳ７）。
【００７５】
ステップＳ６において、注目領域の画像が表す部位が「軟部」でないと判定された場合には、部位判定部６２２は注目領域の画像が表す部位が「肺野部」であるか否かを判定する（ステップＳ８）。ここで、注目領域の画像が表す部位が「肺野部」であると判定された場合には、ゲイン曲線選択部６２３が記憶部６４に記憶されている、複数のゲイン曲線を特定する情報の中から、「肺野部用ゲイン曲線」を読み込んで取得し、画像フィルタ適用部６２４が、この「肺野部用ゲイン曲線」で定義された画像フィルタを注目領域の画像に適用する（ステップＳ９）。
【００７６】
ステップＳ８において、注目領域の画像が表す部位が「肺野部」でないと判定された場合には、注目領域の画像が表す部位は「骨部」であるとみなし、ゲイン曲線選択部６２３が記憶部６４に記憶されている複数のゲイン曲線を特定する情報の中から、「骨部用ゲイン曲線」を読み込んで取得し、画像フィルタ適用部６２４が、この「骨部用ゲイン曲線」で定義された画像フィルタを注目領域の画像に適用する（ステップＳ１０）。
【００７７】
注目領域の画像に対して画像フィルタが適用されると、画像フィルタ制御部６２１が、読取断層像における画像処理の対象領域上で、次の注目領域として設定すべき局所領域が残っているか否かを判定する（ステップＳ１１）。次の注目領域として設定すべき局所領域が残っている場合にはステップ５に戻り、新たな注目領域を設定する。一方、次の注目領域として設定すべき局所領域が残っていない場合には、画像処理を終了する。
【００７８】
このようにして処理済断層像が得られると、表示処理部６３は、その処理済断層像を表示部の画面に表示する（ステップＳ１２）。そして、読影者は、画像解析用ワークステーション６の表示部に表示された処理済断層像を見て読影することになる。
【００７９】
（第２の実施形態）
次に、本発明の第２の実施形態について説明する。
【００８０】
本発明の第２の実施形態によるＸ線ＣＴシステムは、図１に示すような、第１の実施形態によるＸ線ＣＴシステムと略同様の構成であるが、部位判定部６２２の動作が、第１の実施形態とは異なるものである。部位判定部６２２は、第１の実施形態においては、注目領域の画像が表す部位を、肺野部、軟部および骨部の３種類のうちいずれに該当するか判定していたが、この第２の実施形態では、注目領域の画像が表す部位を、肺野部、軟部、骨部および造影剤の４種類のうちいずれに該当するか判定するものである。
【００８１】
Ｘ線ＣＴ撮影では、被検体の血管に造影剤を注入して撮影する場合があるが、この造影剤を表す画像に対しては、造影剤に適した画像フィルタが適用されたり、あるいは、画像フィルタが適用されないようにしたりすることが望ましい場合がある。ところが、造影剤に対応するＣＴ値は、骨に対応するＣＴ値と近接しており、ＣＴ値が約＋２００〜＋５００の範囲において区別が難しく、単純な平均ＣＴ値の閾値判定では両者を正確に判別することが難しい。
【００８２】
その一方で、骨と造影剤とでは、内部構造や物質の均一性に違いがあり、その違いがＣＴ値の空間的な分布の違いとして現れる。
【００８３】
図７（ａ）は、血管（大動脈）に注入された造影剤を表す画像におけるＣＴ値の空間的な変化を示した図であり、図７（ｂ）は、骨を表す画像におけるＣＴ値の空間的な変化を示した図である。例えば、大動脈が造影剤で造影された場合には、図７（ａ）に示すように、平均ＣＴ値が＋５００程度であるとすると、その大動脈の断面におけるＣＴ値のばらつきは、ＣＴ値の変化幅として２０〜３０程度である。このＣＴ値のばらつきは、造影剤の均一性やビームハードニング等による影響が考えられる。また、骨の場合には、図７（ｂ）に示すように、一般的に、表面が骨膜で覆われ、内部には緻密質、更に中心部には海綿質が存在しており、ＣＴ値は、緻密質では比較的高くなり、海綿質では比較的低くなる。
【００８４】
そこで、本実施形態では、部位判定部６２２は、注目領域における平均ＣＴ値の閾値判定だけでなく、上記のような骨と造影剤のＣＴ値の空間的な分布に関する特徴の違いを利用して、注目領域の画像が表す組織を肺野部、軟部、骨部および造影剤の４種類のうちいずれに該当するか判定する。以下、その判定手法の一例を説明する。
【００８５】
まず、読出断層像における各画素に対してＣＴ値の閾値判定を行い、所定値以上のＣＴ値を有する画素領域、すなわち、骨または造影剤に相当するＣＴ値を有する画素領域が抽出された画像Ｐ０を得る。そして、この画像Ｐ０において、２次元的に連続した領域（以下、２次元連続領域という）毎に番号ｚ（ｚ＝１，２，・・・）を付与するラベリング処理を行う。次に、画像Ｐ０に対して、各２次元連続領域Ｒｚを所定画素幅分だけ収縮させる収縮用論理フィルタを、すべての２次元連続領域がなくなるまで段階的に適用し、各段階における画像を、Ｐ１，Ｐ２，・・・として得る。そして、隣接する２画像の組合せ毎に２画像間で減算処理を行い、差分画像Ｐ01，Ｐ12，・・・を得る。すなわち、画像Ｐ０から画像Ｐ１を減算して差分画像Ｐ01を得、画像Ｐ１から画像Ｐ２を減算して差分画像Ｐ12を得るといった減算処理を繰り返す。このようにして得られた複数の差分画像Ｐ01，Ｐ12，・・・に基づいて、図８に示すように、各２次元連続領域Ｒｚ毎に、所定画素幅で層状になった殻状領域Ｒz1，Ｒz2，・・・を認識する。そして、各２次元連続領域Ｒｚ毎に、この層状になった各殻状領域Ｒzi（ｉ＝１，２，・・・）に対して、殻状領域を構成する画素のＣＴ値の平均値Ｍziおよび標準偏差Ｓziを求め、さらに、これらＣＴ値の平均値Ｍziおよび標準偏差Ｓziの標準偏差ＳＭz，ＳＳzを算出し、これらの値が所定の閾値Ｔｈ以上である場合、例えば、ＳＭz＋ＳＳz≧Ｔｈを満たす場合には、その２次元連続領域Ｒｚを骨領域と判別し、上記式を満たさない場合には、その２次元連続領域Ｒｚを造影剤領域と判別する。
【００８６】
このようにして、読出断層像における骨領域と造影剤領域とを区別して抽出し、その位置を把握しておく。後は、設定された注目領域における平均ＣＴ値に基づいて、骨部または造影剤のいずれかであると判定された場合に、その注目領域の位置に基づいて、骨部であるか造影剤であるかを判定するようにすればよい。
【００８７】
次に、第２の実施形態によるＸ線ＣＴシステム１の動作について説明する。
【００８８】
図９は、第２の実施形態によるＸ線ＣＴシステム１における処理の流れを示すフローチャートである。なお、ステップＳ２１〜Ｓ２４は、第１の実施形態におけるステップＳ１〜Ｓ４と同様の処理になるので、ここでは説明を省略する。
【００８９】
データ読出部６１により断層像が読み出されると、画像処理部６２の部位判定部６２２が、上述の判定手法により、読出断層像における骨領域と造影剤領域を区別して抽出し、それぞれの位置を予め把握する（ステップＳ２５）。次に、画像フィルタ制御部６２１は、その読出断層像のうち画像処理の対象領域において、局所領域である注目領域を設定する（ステップＳ２６）。部位判定部６２２は、その注目領域における平均ＣＴ値に基づいて、その注目領域の画像が表す解剖学的な部位の種類が、「軟部」であるか否かを判定する（ステップＳＳ２７）。ここで、注目領域の画像が表す部位が「軟部」であると判定された場合には、ゲイン曲線選択部６２３は、記憶部６４に記憶されている、複数のゲイン曲線を特定する情報の中から、「軟部用ゲイン曲線」の情報を読み込んで取得し、画像フィルタ適用部６２４が、この「軟部用ゲイン曲線」で定義された画像フィルタを注目領域の画像に適用する（ステップＳ２８）。
【００９０】
ステップＳ２７において、注目領域の画像が表す部位が「軟部」でないと判定された場合には、部位判定部６２２は注目領域の画像が表す部位が「肺野部」であるか否かを判定する（ステップＳ２９）。ここで、注目領域の画像が表す部位が「肺野部」であると判定された場合には、ゲイン曲線選択部６２３が記憶部６４に記憶されている、複数のゲイン曲線を特定する情報の中から、「肺野部用ゲイン曲線」を読み込んで取得し、画像フィルタ適用部６２４が、この「肺野部用ゲイン曲線」で定義された画像フィルタを注目領域の画像に適用する（ステップＳ３０）。
【００９１】
ステップＳ２９において、注目領域の画像が表す部位が「肺野部」でないと判定された場合には、注目領域の画像が表す部位は「骨部または造影剤」であるとみなし、さらに、その注目領域の位置に基づいて注目領域が先に求めた骨領域上にあるか否かを判定し（ステップＳ３１）、肯定される場合には、注目領域を「骨部」とみなして、ゲイン曲線選択部６２３が記憶部６４に記憶されている複数のゲイン曲線を特定する情報の中から、「骨部用ゲイン曲線」を読み込んで取得し、画像フィルタ適用部６２４が、この「骨部用ゲイン曲線」で定義された画像フィルタを注目領域の画像に適用する（ステップＳ３２）。一方、ステップＳ３１において否定される場合には、注目領域を「造影剤」とみなして、ゲイン曲線選択部６２３が記憶部６４に記憶されている複数のゲイン曲線を特定する情報の中から、「造影剤用ゲイン曲線」を読み込んで取得し、画像フィルタ適用部６２４が、この「造影剤用ゲイン曲線」で定義された画像フィルタを注目領域の画像に適用する（ステップＳ３３）。
【００９２】
注目領域の画像に対して画像フィルタが適用されると、画像フィルタ制御部６２１が、読取断層像における画像処理の対象領域上で、次の注目領域として設定すべき局所領域が残っているか否かを判定する（ステップＳ３４）。次の注目領域として設定すべき局所領域が残っている場合にはステップ２６に戻り、新たな注目領域を設定する。一方、次の注目領域として設定すべき局所領域が残っていない場合には、画像処理を終了する。
【００９３】
このようにして処理済断層像が得られると、表示処理部６３は、その処理済断層像を表示部の画面に表示する（ステップＳ３５）。
【００９４】
このように、上記第１および第２の実施形態によるＸ線ＣＴシステムでは、高周波域成分を強調して画像再構成された断層像に対して、局所領域毎に、当該領域の解剖学的な部位の種類に応じて、当該領域の画像特徴量と高周波域成分の強調度との対応関係を表すゲイン曲線が異なる画像フィルタを切り換えて適用する画像処理手段を、画像再構成手段（例えば、操作コンソール）ではなく画像表示手段（例えば、画像解析用ワークステーション）が具備し、断層像を保存する手段は、画像フィルタ適用前の断層像を保存するようにしている。すなわち、本発明のＸ線ＣＴシステムでは、高周波域成分が強調された断層像に対して画像特徴量依存型の画像フィルタを適用する画像処理を、断層像の画像再構成と切り離して、断層像を表示する際に行い、断層像の保存は、断層像の画像再構成と切り離さずに行うようにしている。
【００９５】
これにより、断層像を保存する際に画像フィルタ適用前の断層像を保存し、断層像を表示する際に断層像に画像フィルタを適用して画像フィルタ適用後の断層像を表示することができ、高周波域成分を強調する画像再構成によって得られた断層像に対して、局所領域毎に、当該領域の解剖学的な部位の種類に応じて、画像特徴量に依存した画像フィルタを切り換えながら適用することにより画質が読影用に最適化された断層像を表示可能としつつ、高周波域成分が欠落していない画像フィルタ適用前の断層像を表示させたい場合や当該断層像に別の画像処理を施して加工したい場合などに対応できる。
【００９６】
なお、上記の実施形態では、被検体の断層像を画像再構成する手段と、画像再構成によって得られた断層像に画像フィルタを適用して処理済断層像を得る手段とを、Ｘ線ＣＴ装置２の操作コンソール５と画像解析用ワークステーション６とに分離して設けているが、例えば、これら両手段を分離せず、操作コンソール５に１つにまとめて設けるようにしてもよい。
【００９７】
また、上記の実施形態では、画像解析用ワークステーション６が、断層像に画像フィルタを適用して処理済断層像を得る画像処理部６２、および、処理済断層像を表示する表示処理部６３を予め有しているが、例えば、図１０に示すように、操作コンソール５の画像保存部Ｘ線ＣＴ装置２が、上記のような画像処理部６２および表示処理部６３を有しない、画像解析用ワークステーション等のコンピュータに対して、断層像の画像データと、コンピュータをこのような画像処理部および表示処理部として機能させるための画像表示ツール（ｔｏｏｌ）であるプログラム（ｐｒｏｇｒａｍ）とを共に転送し、コンピュータ側でそのプログラムを実行するようにしてもよい。このようにすれば、断層像に画像フィルタを適用して処理済断層像を得る画像処理部およびその処理済断層像を表示する表示処理部とを有しない、画像解析用ワークステーション等のコンピュータであっても、同様の画像処理を行うことができる。
【００９８】
また、上記の実施形態では、断層像に適用する画像フィルタとして２次元の画像空間フィルタを用いているが、例えば、被検体の搬送方向すなわち被検体の体軸方向に連続して並ぶ複数の断層像における当該体軸方向の画素を考慮した、３次元の画像空間フィルタを用いるようにしてもよい。近年、Ｘ線検出器として、多列Ｘ線検出器やフラットパネルＸ線検出器に代表されるマトリクス構造の２次元Ｘ線検出器が用いられるようになり、被検体の体軸方向に連続した複数の断層像からなる３次元画像において、被検体の搬送方向のスライス厚を充分薄くすることができ、被検体の断層面方向であるｘ，ｙ軸および被検体の体軸方向であるｚ軸の各方向に対して画素の分解能が均一な３次元画像を生成することが可能になっている。換言すると、３次元画像の画素分解能の等方性、すなわち、アイソトロピック（ｉｓｏｔｒｏｐｉｃ）が実現されている。したがって、スライス（ｓｌｉｃｅ）厚が充分に薄い３次元画像については、３次元の画像空間フィルタを用いた方が画質の最適化には高い効果を奏する場合がある。特に、画像ノイズ（ｎｏｉｓｅ）改善という観点では、より多くの画素を用いた方が、よい結果が得られるので、２次元の画像空間フィルタよりも、より多くの画素を用いる３次元の画像空間フィルタを用いる方が、効果が大きい。
【００９９】
また、上記の実施形態では、断層像に適用する画像フィルタとして、ＣＴ値に依存した画像フィルタを用いているが、さらに、断層像のスライス厚に依存した画像フィルタを用いる場合も考えられる。
【０１００】
上記の実施形態は、本発明を実施するための最良の形態の一例であり、本発明はこれらの実施形態に限定されない。すなわち、本発明は、本発明の趣旨を逸脱しない限り、あらゆる変更・追加が可能である。
【０１０１】
また、コンピュータを本発明のＸ線ＣＴシステムあるいはその各手段として機能させるためのプログラムおよびそのプログラムが記録された記録媒体も、本発明の実施形態の一例である。
【図面の簡単な説明】
【０１０２】
【図１】発明を実施するための最良の形態例である医用画像表示システムのブロック図
【図２】各種再構成関数による空間周波数と高域強調度との関係を示す図
【図３】ＣＴ値と高周波域成分の強調または抑制の程度を規定する画像フィルタ曲線の一例を示す図
【図４】断層像に適用する画像フィルタの一例を示す図
【図５】断層像に適用する他の画像フィルタの一例を示す図
【図６】第１の実施形態によるＸ線ＣＴシステムにおける処理の流れを示すフローチャート
【図７】骨の画像および造影剤の画像におけるＣＴ値の空間的な変化を示す図
【図８】骨の画像および造影剤の画像において分割された層状の殻状領域を示す図
【図９】第２の実施形態によるＸ線ＣＴシステムにおける処理の流れを示すフローチャート
【図１０】Ｘ線ＣＴ装置から断層像と画像表示ツールとを共にコンピュータに転送する実施形態の概念を示す図
【符号の説明】
【０１０３】
１ Ｘ線ＣＴシステム
２ Ｘ線ＣＴ装置
３ 走査ガントリ
４ 撮影テーブル
５ 操作コンソール
６ 画像解析用ワークステーション
３１ Ｘ線照射装置
３２ Ｘ線検出装置
３３ データ収集部
５１ 制御インターフェース
５２ データ収集バッファ
５３ 記憶装置
５４ データ処理装置（画像再構成手段）
５５ データ保存装置（保存処理手段）
５６ 表示装置
５７ 操作装置
６１ データ読出部（画像読出手段）
６２ 画像処理部（画像処理手段）
６３ 表示処理部（表示処理手段）
６４ 記憶部
６２１ 画像フィルタ制御部
６２２ 部位判定部
６２３ ゲイン曲線選択部
６２４ 画像フィルタ適用部
【特許請求の範囲】
【請求項１】
被検体を走査して収集されたＸ線投影データに基づいて、前記被検体の断層像の高周波域成分を強調する再構成関数を用いて、前記被検体の断層像を画像再構成する画像再構成手段と、
前記画像再構成された断層像を記憶装置または記憶媒体に保存する保存処理手段と、
前記保存処理手段により保存された断層像を読み出す画像読出手段と、
前記読み出された断層像の少なくとも一部における各局所領域に対して、該局所領域の画像特徴量に応じて画像の高周波域成分の強調度が変化する画像フィルタを適用する画像処理であって、前記画像特徴量と前記強調度との対応関係を表すゲイン曲線を、該局所領域が表す解剖学的な部位の種類に応じて切り換える画像処理を施すことにより処理済断層像を得る画像処理手段、および、該処理済断層像を表示する表示処理手段を有する画像表示手段とを備えたことを特徴とするＸ線ＣＴシステム。
【請求項１】
被検体を走査して収集されたＸ線投影データに基づいて、前記被検体の断層像の高周波域成分を強調する再構成関数を用いて、前記被検体の断層像を画像再構成する画像再構成手段と、
前記画像再構成された断層像を記憶装置または記憶媒体に保存する保存処理手段と、
前記保存処理手段により保存された断層像を読み出す画像読出手段と、
前記読み出された断層像の少なくとも一部における各局所領域に対して、該局所領域の画像特徴量に応じて画像の高周波域成分の強調度が変化する画像フィルタを適用する画像処理であって、前記画像特徴量と前記強調度との対応関係を表すゲイン曲線を、該局所領域が表す解剖学的な部位の種類に応じて切り換える画像処理を施すことにより処理済断層像を得る画像処理手段、および、該処理済断層像を表示する表示処理手段を有する画像表示手段とを備えたことを特徴とするＸ線ＣＴシステム。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【公開番号】特開２００８−１４８９７０（Ｐ２００８−１４８９７０Ａ）
【公開日】平成２０年７月３日（２００８．７．３）
【国際特許分類】
【出願番号】特願２００６−３４０７９５（Ｐ２００６−３４０７９５）
【出願日】平成１８年１２月１９日（２００６．１２．１９）
【出願人】（３０００１９２３８）ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー (1,125)
【Ｆターム（参考）】
【公開日】平成２０年７月３日（２００８．７．３）
【国際特許分類】
【出願日】平成１８年１２月１９日（２００６．１２．１９）
【出願人】（３０００１９２３８）ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー (1,125)
【Ｆターム（参考）】
[ Back to top ]