(57)【要約】【課題】対象物内の物質を高い精度で特定可能なＸ線ＣＴ装置を提供する。【解決手段】実施形態のＸ線ＣＴ装置は、対象物をスキャンして得られた投影データに基づいて対象物内の画像を表示するものであり、生成部と、変換部と、画像形成部と、特定部とを有する。生成部は、エネルギーが異なるＸ線でそれぞれ対象物をスキャンして複数の投影データを生成する。変換部は、複数の投影データを、複数の基準物質に対応する複数の新たな投影データに変換する。画像形成部は、変換部で変換した複数の新たな投影データをそれぞれ再構成することで、複数の基準物質に対応する複数の基準物質画像を形成する。特定部は、複数の基準物質画像の画素の値の相関に基づいて対象物質を特定する。【選択図】図１
2013081765.tif 000002
【特許請求の範囲】
【請求項１】
対象物をスキャンして得られた投影データに基づいて、前記対象物内の画像を表示するＸ線ＣＴ装置において、
エネルギーが異なるＸ線でそれぞれ対象物をスキャンして複数の投影データを生成する生成部と、
前記複数の投影データを、複数の基準物質に対応する複数の新たな投影データに変換する変換部と、
前記変換部で変換した複数の新たな投影データをそれぞれ再構成することで、複数の基準物質に対応する複数の基準物質画像を形成する画像形成部と、
前記複数の基準物質画像の画素の値の相関に基づいて対象物質を特定する特定部と、
を備えるＸ線ＣＴ装置。
【請求項２】
前記特定部は、前記複数の基準物質画像における対応する画素の値が、予め設定された相関を有するかを判定することにより、前記対象物質を特定するものであることを特徴とする請求項１記載のＸ線ＣＴ装置。
【請求項３】
前記変換部は、前記基準物質に対応するＸ線減弱係数を含む計算式を用いて、前記複数の投影データを、前記複数の基準物質に対応する複数の新たな投影データに変換するものであることを特徴とする請求項１または２に記載のＸ線ＣＴ装置。
【請求項４】
前記特定部により特定された前記対象物質の領域を識別可能にした前記対象物内の画像を生成する表示画像生成部を備えることを特徴とする請求項１乃至３のいずれか一項に記載のＸ線ＣＴ装置。
【請求項５】
前記特定部により特定された前記対象物質の領域の画素を、他の領域の画素に比べて、強調または抑制した表示画像を生成する表示画像生成部を備えることを特徴とする請求項１乃至３のいずれか一項に記載のＸ線ＣＴ装置。
【請求項６】
予め記憶された計算式に基づいて、前記複数の基準物質画像における対応する画素の値から、前記対象物質の濃度情報を求める濃度算出部を備えることを特徴とする請求項１乃至３のいずれか一項に記載のＸ線ＣＴ装置。
【請求項７】
前記濃度情報を表した表示画像を生成する表示画像生成部を備えることを特徴とする請求項６記載のＸ線ＣＴ装置。
【請求項８】
前記画像形成部は、前記複数の投影データを再構成して得られた複数のエネルギー画像を生成するものであり、
前記特定部は、前記複数のエネルギー画像の前記複数の基準物質画像の画素の値の相関に基づいて対象物質を特定するものである、
ことを特徴とする請求項１乃至７のいずれか一項に記載のＸ線ＣＴ装置。
【請求項９】
前記基準物質画像に基づいて、前記対象物における実効原子番号の分布を表す実効原子番号画像、前記対象物質における物質の密度の分布を表す密度画像、及び、単一のエネルギーを有するＸ線で前記対象物質をスキャンした場合を模擬的に再現する単色Ｘ線画像のうちの少なくとも１つを形成する形成部を更に備え、
前記特定部は、前記形成部により形成された画像に基づいて、前記対象物質の特定を行う、
ことを特徴とする請求項１に記載のＸ線ＣＴ装置。
【請求項１０】
エネルギーが異なるＸ線でそれぞれ対象物をスキャンして生成された複数の投影データを、複数の基準物質に基づき変換して得られた複数の新たな投影データを、それぞれ再構成することで形成された、複数の基準物質に対応する複数の基準物質画像を記憶する記憶部と、
前記記憶部に記憶された前記複数の基準物質画像の画素の値の相関に基づいて対象物質を特定する特定部と、
を備える画像処理装置。
【請求項１１】
エネルギーが異なるＸ線でそれぞれ対象物をスキャンして生成された複数の投影データを、複数の基準物質に基づき変換して得られた複数の新たな投影データを、それぞれ再構成することで形成された、複数の基準物質に対応する複数の基準物質画像を記憶する記憶部を有するコンピュータを、
前記記憶部に記憶された前記複数の基準物質画像の画素の値の相関に基づいて対象物質を特定する特定部として機能させるプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明の実施形態は、Ｘ線ＣＴ装置、画像処理装置及びプログラムに関する。
【背景技術】
【０００２】
Ｘ線コンピュータ断層撮影（Ｃｏｍｐｕｔｅｄ  Ｔｏｍｏｇｒａｐｈｙ：ＣＴ）は、対象物をＸ線ビームでスキャンして得られた投影データを再構成することにより、対象物の情報を表す画像を形成する技術である。
【０００３】
Ｘ線ＣＴの応用として、対象物に含まれる物質の種別を特定する技術がある。この技術においては、単一の管電圧によるＸ線ビームで得られた画像を基に物質を判別することは難しいため、異なる２つの管電圧によるエネルギーの異なる２つのＸ線ビームでそれぞれスキャンを行う手法が近年注目を集めている。この手法は「デュアル・エナジー・ＣＴ（Ｄｕａｌ  Ｅｎｅｒｇｙ  ＣＴ）」と呼ばれる。
【０００４】
非特許文献１には、異なる２つの管電圧を適用して２つの画像を形成し、これら画像のＣＴ値の比を用いて物質を特定する技術が記載されている。なお、非特許文献１に記載された全ての事項は、この明細書の一部を構成するものとして援用される。
【０００５】
また、特許文献１には、異なる２つの管電圧を適用して得られた各投影データに示す線減弱係数を、２つの基準物質（たとえば水と骨）の線減弱係数の線形結合として表現することにより、各基準物質の分布を表す画像（基準物質画像）を形成する技術が記載されている。更に、特許文献１には、これら基準物質画像を組み合わせることにより、実効原子番号画像、密度画像及び単色Ｘ線画像を形成する手法も記載されている。なお、特許文献１に記載された全ての事項は、この明細書の一部を構成するものとして援用される。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２００９−２６１９４２号公報
【非特許文献】
【０００７】
【非特許文献１】Ｊｏｈｎｓｏｎ  ＴＲ．  ｅｔ  ａｌ．、“Ｍａｔｅｒｉａｌ  ｄｉｆｆｅｒｅｎｔｉａｔｉｏｎ  ｂｙ  ｄｕａｌ  ｅｎｅｒｇｙ  ＣＴ：ｉｎｉｔｉａｌ  ｅｘｐｅｒｉｅｎｃｅ”、Ｅｕｒ  Ｒａｄｉｏｌ（２００７）、１７、１５１０−１５１７
【発明の概要】
【発明が解決しようとする課題】
【０００８】
しかしながら、ＣＴ値が管電圧に依存し、その結果ＣＴ値の比が管電圧の組み合わせに依存することを考慮すると、従来のデュアル・エナジー・ＣＴにおいては、ＣＴ値比が近い物質を高い精度で判別することは困難であった。
【０００９】
また、従来のデュアル・エナジー・ＣＴにおいては、基準物質画像に基づいてその基準物質の含有量が多いか少ないかを把握することは可能であるが、その物質の種別を特定することは困難であった。
【００１０】
この発明が解決しようとする課題は、対象物に含まれる物質を高い精度で特定することが可能なＸ線ＣＴ装置、画像処理装置及びプログラムを提供することである。
【課題を解決するための手段】
【００１１】
実施形態のＸ線ＣＴ装置は、対象物をスキャンして得られた投影データに基づいて対象物内の画像を表示するものであり、生成部と、変換部と、画像形成部と、特定部とを有する。生成部は、エネルギーが異なるＸ線でそれぞれ対象物をスキャンして複数の投影データを生成する。変換部は、複数の投影データを、複数の基準物質に対応する複数の新たな投影データに変換する。画像形成部は、変換部で変換した複数の新たな投影データをそれぞれ再構成することで、複数の基準物質に対応する複数の基準物質画像を形成する。特定部は、複数の基準物質画像の画素の値の相関に基づいて対象物質を特定する。
【図面の簡単な説明】
【００１２】
【図１】第１の実施形態に係るＸ線ＣＴ装置の概略構成の一例を表すブロック図である。
【図２】第１の実施形態に係るＸ線ＣＴ装置が実行する処理を説明するためのグラフである。
【図３】第１の実施形態に係るＸ線ＣＴ装置が実行する処理を説明するためのグラフである。
【図４】第１の実施形態に係るＸ線ＣＴ装置が実行する処理を説明するためのグラフである。
【図５】第１の実施形態に係るＸ線ＣＴ装置が実行する処理を説明するためのグラフである。
【図６】第１の実施形態に係るＸ線ＣＴ装置の動作の一例を表すフローチャートである。
【図７】第２の実施形態に係るＸ線ＣＴ装置の概略構成の一例を表すブロック図である。
【図８】第２の実施形態に係るＸ線ＣＴ装置が表示する画像の一例を表す概略図である。
【図９Ａ】第２の実施形態に係るＸ線ＣＴ装置が表示する画像の一例を表す概略図である。
【図９Ｂ】第２の実施形態に係るＸ線ＣＴ装置が表示する画像の一例を表す概略図である。
【図９Ｃ】第２の実施形態に係るＸ線ＣＴ装置が表示する画像の一例を表す概略図である。
【図１０】第２の実施形態に係るＸ線ＣＴ装置が表示する画像の一例を表す概略図である。
【図１１】第３の実施形態に係るＸ線ＣＴ装置が表示する画像の一例を表すブロック図である。
【図１２】第３の実施形態に係るＸ線ＣＴ装置が実行する処理を説明するための概略図である。
【図１３】第３の実施形態に係るＸ線ＣＴ装置が実行する処理を説明するための概略図である。
【図１４】第３の実施形態に係るＸ線ＣＴ装置の動作の一例を表すフローチャートである。
【図１５Ａ】Ｘ線ＣＴ装置の変形例を説明するための概略図である。
【図１５Ｂ】Ｘ線ＣＴ装置の変形例を説明するための概略図である。
【図１６】実施形態に係る画像処理装置の概略構成の一例を表すブロック図である。
【発明を実施するための形態】
【００１３】
実施形態に係るＸ線ＣＴ装置、画像処理装置及びプログラムについて図面を参照しつつ説明する。以下の実施形態では被検体（患者）を対象物とするが、対象物はこれに限定されるものではない。
【００１４】
〈第１の実施形態〉
［構成］
この実施形態のＸ線ＣＴ装置の概略構成を図１に示す。このＸ線ＣＴ装置は、ガントリ１と、前処理部２と、再構成処理部３と、データ処理部４とを有する。なお、図示は省略するが、この実施形態のＸ線ＣＴ装置には、一般的なものと同様に、寝台、コンソール、高電圧発生装置などが設けられている。
【００１５】
〔ガントリ１〕
ガントリ１は、被検体をＸ線でスキャンするために用いられる。ガントリ１には、一般的なものと同様に、互いに対向配置されたＸ線管及びＸ線検出器、これらを回転させる回転機構、スリップリング機構、傾斜（チルト）機構、データ収集部（Ｄａｔａ  Ａｃｑｕｓｉｔｉｏｎ  Ｓｙｓｔｅｍ：ＤＡＳ）などが設けられている。ガントリ１は、Ｘ線管及びＸ線検出器を回転させつつ被検体をＸ線でスキャンする。Ｘ線検出器による検出データは、データ収集部によって収集され、前処理部２に送られる。
【００１６】
ガントリ１は、特に、エネルギーの異なる２つのＸ線ビームでそれぞれスキャンを行う手法、つまりデュアル・エナジー・ＣＴを実施することができる。Ｘ線のエネルギーは、高電圧発生装置がＸ線管に印加する電圧（管電圧）に依存する。デュアル・エナジー・ＣＴの手法としては、Ｓｌｏｗ−ｋＶ  Ｓｗｉｔｃｈｉｎｇ方式や、Ｄｕａｌ  Ｓｏｕｒｃｅ方式や、Ｆａｓｔ−ｋＶ  Ｓｗｉｔｃｈｉｎｇ方式などがある。Ｓｌｏｗ−ｋＶ  Ｓｗｉｔｃｈｉｎｇ方式とは、第１の管電圧でスキャンを行った後に、第２の管電圧でスキャンを行う方式（２回転方式）である。Ｄｕａｌ  Ｓｏｕｒｃｅ方式とは、Ｘ線管を２つ備えたガントリを用い、これらＸ線管に異なる管電圧を印加してスキャンを行う方式（２管球方式）である。Ｆａｓｔ−ｋＶ  Ｓｗｉｔｃｈｉｎｇ方式とは、Ｘ線管及びＸ線検出器を回転させながらビューごとに管電圧を切り替える方式（高速スイッチ方式）である。
【００１７】
〔前処理部２〕
前処理部２は、ガントリ１から送られた検出データに対して所定の前処理（画像再構成処理の前に行われる処理）を施す。この前処理としては、データの対数を計算する処理、リファレンス補正、水補正、ビームハードニング補正、体動補正などがある。前処理部２により生成されるデータは投影データと呼ばれる。前処理部２により生成された投影データは、再構成処理部３やデータ処理部４に送られる。ガントリ１及び前処理部２は「生成部」の一例として機能する。
【００１８】
〔再構成処理部３〕
再構成処理部３は、前処理部２により生成された投影データに対して再構成処理を施すことにより、被検体の画像データを生成する。再構成処理は、投影データから被検体のＸ線吸収係数の分布を逆算する演算処理である。この演算処理としては、２次元フーリエ変換法、コンボリューション・バックプロジェクション法、ファンビーム・コンボリューション・バックプロジェクション法などがある。
【００１９】
また、再構成処理部３は、データ処理部４により得られた投影データに再構成処理を施して画像データを生成するように構成されていてもよい。この処理については後述する。
【００２０】
〔データ処理部〕
データ処理部４は、前処理部２により生成された投影データに対して所定のデータ処理を施すことにより、被検体に含まれる物質の種別を特定する。
【００２１】
データ処理部４には、投影データ変換部４１と、画像形成部４２と、物質特定部４３とが設けられている。なお、投影データ変換部４１は「変換部」の一例として、画像形成部４２は「画像形成部」の一例として、そして物質特定部４３は「特定部」の一例として、それぞれ機能する。
【００２２】
（投影データ変換部４１）
投影データ変換部４１は、デュアル・エナジー・ＣＴの手法により得られた第１及び第２の投影データを、あらかじめ決められた２つの基準物質に対応する２つの投影データに変換する。この処理の一例として、投影データ変換部４１は、特許文献１に記載された手法を用いて、第１及び第２の投影データのそれぞれを、２つの基準物質に対応するあらかじめ決められた２つの基準値と、当該２つの投影データとからなる線形結合として表現する。つまり、この線形結合における２つの係数が、目的の２つの投影データとなる。この処理例の詳細については後述する。
【００２３】
（画像形成部４２）
画像形成部４２は、投影データ変換部４１により得られた２つの投影データをそれぞれ再構成することにより、２つの基準物質に対応する２つの基準物質画像を形成する。基準物質画像は、対象物内における線減弱係数の分布を表す。線減弱係数とは、入射Ｘ線が単一厚さの物質を通過するときに減衰するエネルギーの割合を示す。この再構成処理により得られる情報は、２つの基準値の係数がと２つの基準物質画像とからなる線形結合画像である。つまり、この線形結合画像における２つの係数が、目的の２つの基準物質画像（再構成画像）となる。
【００２４】
この再構成処理は、再構成処理部３と同様の要領で行われる。なお、画像形成部４２の代わりに再構成処理部３が当該再構成処理を行うようにしてもよい。その場合、画像形成部４２は不要である。画像形成部４２又は再構成処理部３が実行する処理の例の詳細については後述する。
【００２５】
なお、２つの基準値は、投影データの線減弱係数を線形結合として表すために用いられる。２つの基準値は、２つの基準物質の線減弱係数であってよい。また、２つの基準値は任意の値であってもよい。
【００２６】
（処理例）
上記の処理例として、特許文献１に記載された手法を説明する。前処理部２により生成された第１の投影データ（高エネルギー側）をｇＨで表し、第２の投影データ（低エネルギー側）をｇＬで表す。投影データ変換部４１は、これら投影データｇＨ、ｇＬに対して次式（１）に示す変換を施すことで、２つの投影データＬ１、Ｌ２を生成する。
【００２７】
【数１】
2013081765.tif 000003
【００２８】
ただし：
Ｄは、次式（２）の右辺の２×２行列の行列式；
〈μ〉１，２Ｈ，Ｌは、特許文献１に記載されたエネルギー平均化線減弱係数；
をそれぞれ示す。
【００２９】
【数２】
2013081765.tif 000004
【００３０】
画像形成部４２は、式（２）に示す２つの投影データＬ１、Ｌ２をそれぞれ再構成することにより、２つの基準物質に対応する２つの基準物質画像（次式（３）のｃ１、ｃ２）を形成する。任意の物質の線減弱係数μは、次式（３）に示すように、２つの基準値μ１、μ２と、２つの基準物質画像ｃ１、ｃ２とを用いた線形結合（線形結合画像）として表現される。
【００３１】
【数３】
2013081765.tif 000005
【００３２】
ただし：
Ｅは、Ｘ線のエネルギー；
μ１（Ｅ）は、エネルギーＥにおける第１の基準物質の線減弱係数（基準値）；
μ２（Ｅ）は、エネルギーＥにおける第２の基準物質の線減弱係数（基準値）；
ｃ１（Ｅ，ｘ，ｙ）は、座標（ｘ，ｙ）に位置する画素における第１の基準物質の存在率；
ｃ２（Ｅ，ｘ，ｙ）は、座標（ｘ，ｙ）に位置する画素における第２の基準物質の存在率；
をそれぞれ示す。
【００３３】
基準物質画像（存在率）ｃ１及びｃ２は、任意の物質の線減弱係数μを２つの基準物質の線減弱係数μ１及びμ２の関数として表したときの係数であり、この任意の物質が各基準物質にどれだけ類似しているかを示す指標である。
【００３４】
以下、第１の基準物質が造影剤（ヨウ素濃度５０［ｍｇＩ／ｍｌ］）であり、第２の基準物質が水である場合を例として特に説明する。
【００３５】
（物質特定部４３）
物質特定部４３は、画像形成部４２により形成された２つの基準物質画像の相関に基づいて対象物質の種別を特定する。なお、対象物質とは、この実施形態において種別の特定処理の対象とされる物質を示す。上記処理例を適用する場合、物質特定部４３は、まず、所定の対象物質について、２つの基準値μ１及びμ２に基づきあらかじめ決められた座標系における、２つの基準物質画像ｃ１及びｃ２に対応する座標を決定する。更に、物質特定部４３は、所定の複数の物質についてあらかじめ得られた上記座標系における複数の座標と、前段の処理で決定された２つの基準物質画像ｃ１及びｃ２に対応する座標とに基づいて、対象物質の種別を特定する。
【００３６】
上記座標系については、任意にこれを設定することが可能である。その一例として、図２に示すような、２つの存在率ｃ１及びｃ２を基底とする２次元座標系を適用できる。この座標系は、第１の基準物質としての造影剤の存在率ｃ１を縦軸とし、第２の基準物質としての水の存在率ｃ２を横軸としたものである。
【００３７】
この座標系における座標については、縦軸の座標、横軸の座標の順で記載する。すなわち、この座標系における座標は（ｃ１，ｃ２）と記載される。縦軸における基準の座標Ｐ（１，０）、つまり縦軸方向における基底として、造影剤の存在率ｃ１が１００％に相当するベクトルを適用する。また、横軸における基準の座標Ｑ（０，１）、つまり横軸方向における基底として、水の存在率ｃ２が１００％に相当するベクトルを適用する。
【００３８】
物質特定部４３は、画像形成部４２により得られた線形結合画像における基準物質画像の組（ｃ１，ｃ２）を、この座標系の座標として扱う。
【００３９】
なお、本例では、基準物質画像の組と座標とが同じ表現となるように座標系（基底）を設定しているが、これに限定されるものではない。たとえば、ｃ１やｃ２に係数を掛けることにより、座標軸が伸縮された座標系を適用することができる。また、本例では、縦軸にｃ１を割り当て、横軸にｃ２を割り当てているが、これを逆にしてもよい。更に、本例のような直交座標系に代えて、斜交座標系等の他の座標系を用いることも可能である。すなわち、この実施形態における座標系は、線形結合画像における２つの基準物質画像の相関を表すことができるものであれば十分であり、その具体的な態様は任意である。
【００４０】
物質特定部４３による処理に供される上記「複数の物質」は、任意の物質でよい。また、その個数も任意である。図２に示す例では、造影剤と水が複数の物質に相当する。本例では、「複数の物質」と、基底の生成に用いられる物質とが同じであるが、これには限定されない。その一例として、造影剤と水の２つの物質を「基準物質」とし、造影剤、水、炭酸カルシウム、脂肪及び尿酸の５つの物質を「複数の物質」として適用する場合を後述する。
【００４１】
複数の物質に対応する座標の取得方法についても任意である。その例として、その物質の測定を実際に行って座標を求める方法や、他の物質の線減弱係数をも考慮して座標を算出する方法などがある。
【００４２】
前者の方法は、たとえば、その物質をＸ線でスキャンして投影データを生成し、これに基づく画素毎の線減弱係数を線形結合として表現し、その係数の組に対応する座標を決定することにより行われる。
【００４３】
後者の方法について説明する。各物質の線減弱係数は既知であるから、異なる２つのエネルギーに対応する線減弱係数を上記の式（３）にそれぞれ代入することで、次の連立方程式（４ａ）、（４ｂ）が得られる。
【００４４】
【数４】
2013081765.tif 000006
【００４５】
ただし：
ＥＬｏｗ、ＥＨｉｇｈは、２種類のＸ線エネルギー；
μ（ＥＬｏｗ）は、低い側のエネルギーＥＬｏｗにおける当該物質の線減弱係数；
μ（ＥＨｉｇｈ）は、高い側のエネルギーＥＨｉｇｈにおける当該物質の線減弱係数；
をそれぞれ示す。
【００４６】
この連立方程式（４ａ）、（４ｂ）において未知数はｃ１とｃ２の２つであるから、これを解くことにより、係数の組ｃ１、ｃ２を算出できる。そして、この係数の組に基づいて目的の座標が得られる。
【００４７】
物質特定部４３は、このようにして取得された複数の物質についての複数の座標と、対象物質について物質特定部４３により決定された座標とに基づいて、この対象物質の種別を特定する。この処理の一例として、物質特定部４３は、（１）当該複数の座標に基づく座標系内の領域を求め、（２）この領域に対する対象物質の座標の位置に基づいて、対象物質の種別を特定する。以下、座標系内の領域として、複数の座標を結ぶグラフを用いる場合を説明する。ここで、特に必要な場合を除き、（１）で得られる領域（グラフ等）を実際に表示させる必要はない。以下の説明では、実施形態の理解を支援するために、表示されたグラフを用いているに過ぎない。それ以降に説明される座標系内の領域についても同様である。なお、座標系内の領域は当該座標系における座標の集合に相当し、グラフはこの領域自体やその外周線に位置する座標の集合に相当する。
【００４８】
座標系内の領域としてのグラフを求める処理について説明する。図２に示す例では、上記複数の座標として、造影剤の座標と水の座標が適用されている。物質特定部４３は、目的のグラフとして、これら座標を結ぶ線分Ｌ１を求める。
【００４９】
一般に、線分の式は、２つの座標に基づいてこれらを通る直線の式（つまり傾きと切片）を算出し、この直線のうち当該２つの座標を両端とする部分を抽出することにより得られる。３つ以上の座標を考慮して２つ以上の線分を求める場合には、任意の２つの座標の組み合わせに対して同様の演算を行えばよい。
【００５０】
さて、線分Ｌ１上の座標は、造影剤のヨウ素濃度に対応している。たとえば、線分Ｌ１上の座標Ｐ１、Ｐ２、Ｐ３、Ｐ４、Ｐ５は、ヨウ素濃度２５［ｍｇＩ／ｍｌ］、２０［ｍｇＩ／ｍｌ］、１５［ｍｇＩ／ｍｌ］、１０［ｍｇＩ／ｍｌ］、５［ｍｇＩ／ｍｌ］にそれぞれ対応している。また、前述のように、座標Ｐは５０［ｍｇＩ／ｍｌ］に対応し、座標Ｑは０［ｍｇＩ／ｍｌ］（単なる水）に対応している。つまり、線分Ｌ１上において座標Ｐに近づくとヨウ素濃度が増加し、座標Ｑに近づくとヨウ素濃度が減少する。
【００５１】
このように、複数の物質の１つとして水を採用することで、物質の濃度をグラフとして表すことができる。たとえば造影剤、水、炭酸カルシウム、脂肪及び尿酸の５つの物質が採用される場合におけるグラフの例を図３に示す。
【００５２】
図３には、造影剤の座標Ｐと、水の座標Ｑと、炭酸カルシウムの座標Ｒと、脂肪の座標Ｓと、尿酸の座標Ｔが示されている。また、座標Ｐと座標Ｑとを結ぶ線分Ｌ１は造影剤のヨウ素濃度に相当し、座標Ｒと座標Ｑとを結ぶ線分Ｌ２は炭酸カルシウムの濃度に相当し、座標Ｓと座標Ｑとを結ぶ線分Ｌ３は脂肪の濃度に相当し、座標Ｔと座標Ｑとを結ぶ線分Ｌ４は尿酸の濃度に相当する。各線分Ｌ１〜Ｌ４において、座標Ｑに近づくほど濃度が低くなり、座標Ｑから遠ざかるほど濃度が高くなる。
【００５３】
また、図３から分かるように、４つの線分Ｌ１〜Ｌ４は、互いに座標Ｑでのみ交差している。これは、濃度０（単なる水）に相当する座標Ｑを除いて、線分Ｌ１〜Ｌ４が互いに分離されていること、つまり４つの物質（造影剤、炭酸カルシウム、脂肪及び尿酸）に対応する座標が当該座標系において互いに分離されていることを示している。換言すると、当該座標系における線分の位置が異なれば、それに対応する物質も異なるということになる。
【００５４】
このようなグラフを利用することにより、物質特定部４３は、物質特定部４３により決定された対象物質の座標とグラフとの位置関係に基づいて、この対象物質の種別を特定する。
【００５５】
図３に示す線分Ｌ１〜Ｌ４をグラフとして利用する場合を例として説明すると、物質特定部４３は、まず、対象物質の座標がこれら線分Ｌ１〜Ｌ４のいずれかの上に位置するか判断する。
【００５６】
対象物質の座標が線分Ｌｉ（ｉ＝１〜４のいずれか）上に位置すると判断された場合、物質特定部４３は、この線分Ｌｉにおける当該座標の位置に基づいて、対象物質の種別を特定する。この種別には、対象物質の物質名だけでなく、その濃度（つまり対象物質と水との成分比）も含まれる。
【００５７】
なお、物質名だけを特定するように構成することも可能である。その場合には、対象物質の座標が位置している線分Ｌｉを特定するだけでよい。
【００５８】
一方、対象物質の座標がどの線分Ｌ１〜Ｌ４上にも位置しないと判断された場合、物質特定部４３は、この対象物質はこれら線分Ｌ１〜Ｌ４に相当する物質のいずれにも該当しないとの結果を得る。
【００５９】
３つ以上の基準物質について座標が得られている場合、座標系内の領域として、これら座標を結んでなる多角形を形成することができる。この多角形により規定される領域は、これら基準物質の混合物に相当する。なお、この混合物とは、各基準物質の成分比を０〜１００％とし、全基準物質の成分比の和を１００％とした場合に定義される物質を意味する。よって、この混合物には、３つ以上の基準物質のうちの１つ又は２つのみを成分とするものも含まれる。１つの基準物質のみを成分とするものは多角形の頂点に相当し、２つの基準物質のみを成分とするものは多角形の辺に相当する。
【００６０】
３つ以上の基準物質を用いる場合の例として、造影剤と水と炭酸カルシウムの３つを基準物質とする場合について説明する。この場合、図４に示すように、造影剤、水、炭酸カルシウムに相当する座標Ｐ、Ｑ、Ｒを頂点とする多角形（三角形）Ａが得られる。多角形Ａの辺ＰＱ（線分Ｌ１）及びＱＲ（線分Ｌ２）に相当する領域は、それぞれ前述のように造影剤の濃度及び炭酸カルシウムの濃度を表す。また、辺ＰＲ（線分Ｌ５）に相当する領域は、造影剤と炭酸カルシウムとの混合物における造影剤と炭酸カルシウムとの成分比を表す。辺ＰＲにおいて、頂点Ｐに近づくほど造影剤の成分比が増大し、頂点Ｒに近づくほど炭酸カルシウムの成分比が増大する。
【００６１】
また、多角形Ａの内部領域、つまり多角形Ａ上の領域から線分Ｌ１、Ｌ２、Ｌ５を除いた領域は、これら３つの基準物質全ての混合物に相当する。この内部領域中の座標についても、頂点Ｐに近づくほど造影剤の成分比が増大し、頂点Ｒに近づくほど炭酸カルシウムの成分比が増大する。そして、水に相当する頂点Ｑに近づくほど当該混合物の濃度が減少する。
【００６２】
物質特定部４３は、決定された対象物質の座標が多角形上に位置するか否か判断する。当該座標が多角形の外部に位置すると判断された場合、物質特定部４３は、この対象物質は当該多角形に対応する混合物ではないと判断する。
【００６３】
一方、当該座標が多角形上に位置すると判断された場合、物質特定部４３は、この対象物質は当該混合物であると判断する。更に、物質特定部４３は、この対象物質の座標の位置に基づいて、この対象物質を構成する３つ以上の基準物質の成分比を求める。
【００６４】
なお、以上の例では、Ｘ線ＣＴ装置の測定誤差を考慮していないため、各物質に相当する座標が一意に決定されると仮定している。この測定誤差は、装置の機差（公差等）やノイズの混入などに起因する。以下、測定誤差を考慮する場合における処理について説明する。なお、測定誤差が十分に小さい場合などこれを許容できる場合には、上記の処理を適用すればよい。
【００６５】
測定誤差を考慮する場合の一例について説明する。まず、様々な物質に対して繰り返し測定を行うことにより、各物質についての座標の分布を取得する。この分布情報は、たとえば、このＸ線ＣＴ装置により得られる投影データに混入するノイズの標準偏差情報である。この分布情報は物質特定部４３に記憶される。
【００６６】
造影剤、水及び炭酸カルシウムのそれぞれについての座標の分布の例を図５に示す。造影剤についての座標の分布範囲をＰａ、水についての座標の分布範囲をＱａ、炭酸カルシウムについての座標の分布範囲をＲａでそれぞれ示す。
【００６７】
物質特定部４３は、標準偏差情報に基づいて、当該物質に対応するグラフを含む２次元領域を求める。この処理の一例を説明する。まず、物質特定部４３は、分布範囲Ｐａと分布範囲Ｑａとを結ぶ２本の線分Ｌ１ａ、Ｌ１ｂを求める。線分Ｌ１ａ、Ｌ１ｂとしては、たとえば、分布範囲Ｐａ、Ｑａの双方に接し、かつ互いに交差しないものが用いられる。それにより、分布範囲Ｐａ、Ｑａ及び線分Ｌ１ａ、Ｌ１ｂで囲まれた領域Ｂが得られる。同様にして、分布範囲Ｒａ、Ｑａ及び線分Ｌ２ａ、Ｌ２ｂで囲まれた領域Ｃが得られる。
【００６８】
領域Ｂは、測定誤差を反映した造影剤の濃度の座標の分布範囲として用いられる。領域Ｂにおいて線分Ｌ１上にない座標には、線分Ｌ１に基づいて濃度の値が対応付けられる。この対応付けの例として、線分Ｌ１上の各位置において線分Ｌ１に直交する直線を求め、この直線上に位置する座標の濃度を等しく設定することができる。測定誤差を反映した炭酸カルシウムの濃度の分布範囲を示す領域Ｃについても同様に濃度が設定される。
【００６９】
物質特定部４３は、決定された対象物質の座標と、領域Ｂ、Ｃとの位置関係に基づいて、この対象物質の種別を特定する。たとえば、対象物質の座標が領域Ｂ上に位置する場合、物質特定部４３は、この対象物質は造影剤であると特定し、更に、当該座標と、領域Ｂ上の座標に設定された濃度とに基づいて、この対象物質の濃度を特定する。
【００７０】
［動作］
この実施形態に係るＸ線ＣＴ装置の動作について説明する。このＸ線ＣＴ装置の動作例を図６に示す。なお、装置各部の動作の詳細については上述したので、ここでは簡単な説明にとどめる。
【００７１】
（Ｓ１）
まず、ガントリ１でデュアル・エナジー・ＣＴによる撮影を行う。前処理部２は、ガントリ１により収集されたデータを投影データに変換する。それにより、エネルギーが異なる第１及び第２の投影データが生成される。第１及び第２の投影データは、データ処理部４に送られる。
【００７２】
（Ｓ２）
投影データ変換部４１は、ステップ１で生成された各投影データに基づく画素毎の線減弱係数を２つの基準物質の線減弱係数の線形結合として表現することにより、２つの基準物質による各投影データの分解を行う。このステップ２により、第１及び第２の投影データが、２つの基準物質に対応する２つの投影データに変換される。
【００７３】
（Ｓ３）
更に、画像形成部４２は、ステップ２で得られた線形結合を再構成することにより線形結合画像を形成する。これにより、２つの基準物質画像が得られる。
【００７４】
（Ｓ４）
物質特定部４３は、ステップ３で得られた２つの基準物質画像の相関を求める。上記の例では、２つの基準物質画像の各画素（ｘ、ｙ）について、所定の２次元座標系における対応する座標（ｃ１（ｘ、ｙ）、ｃ２（ｘ、ｙ））が得られる。
【００７５】
（Ｓ５）
物質特定部４３は、ステップ４で得られた相関に基づいて、対象物質を２つの基準物質に分解する。上記の例では、対象物質が造影剤と水などに分解される。
【００７６】
（Ｓ６）
物質特定部４３は、ステップ５での分解結果に基づいて対象物質の種別を特定する。この特定結果は、たとえば図示しないディスプレイに表示される。また、この特定結果は、Ｘ線ＣＴ装置の記憶デバイスやネットワーク上の記憶デバイスに記憶される。
【００７７】
［作用・効果］
この実施形態に係るＸ線ＣＴ装置の作用と効果について説明する。
【００７８】
このＸ線ＣＴ装置は、ガントリ１及び前処理部２により投影データを生成する。特に、デュアル・エナジー・ＣＴ、つまり、エネルギーが異なる第１及び第２のＸ線でそれぞれ対象物をスキャンすることで、第１及び第２の投影データが生成される。
【００７９】
更に、このＸ線ＣＴ装置は、投影データ変換部４１と、画像形成部４２と、物質特定部４３とを有する。投影データ変換部４１は、第１及び第２の投影データを、あらかじめ決められた２つの基準物質に対応する２つの投影データ（新たな投影データ）に変換する。画像形成部４２は、２つの投影データをそれぞれ再構成することで、２つの基準物質に対応する２つの基準物質画像を形成する。物質特定部４３は、２つの基準物質画像の相関に基づいて対象物質の種別を特定する。
【００８０】
投影データ変換部４１が、上記の変換処理として、第１及び第２の投影データのそれぞれを、２つの基準物質に対応するあらかじめ決められた２つの基準値と、２つの投影データとからなる線形結合として表現するように構成されていてもよい。
【００８１】
また、画像形成部４２が、この線形結合を再構成することにより、２つの基準値と２つの基準物質画像とからなる線形結合画像を形成するように構成されていてもよい。
【００８２】
また、物質特定部４３が、対象物質について、２つの基準値に基づきあらかじめ決められた座標系における２つの基準物質画像に対応する座標を決定し、かつ、複数の物質についてあらかじめ得られた座標系における複数の座標と、決定された座標とに基づいて、この対象物質の種別を特定するように構成されていてもよい。
【００８３】
なお、第１及び第２のＸ線のそれぞれのエネルギーについて、２つの基準値は、あらかじめ決められた２つの基準物質の線減弱係数であってよい。更に、座標系は、２つの基準物質の線減弱係数に基づく２つの基底により張られる２次元座標系であってよい。
【００８４】
また、物質特定部４３が、複数の座標に基づく座標系内の領域に対する対象物質の座標の位置に基づいて対象物質の種別を特定するように構成されていてもよい。
【００８５】
また、複数の物質が第１及び第２の物質を含む場合において、物質特定部４３が、座標系における第１の物質の座標と第２の物質の座標とを結ぶ線分を当該座標系内の領域として求め、かつ、対象物質の座標がこの線分上に位置する場合、当該位置に基づいて第１の物質と第２の物質との成分比を求めるように構成されていてもよい。
【００８６】
また、複数の物質が３つ以上の物質を含む場合において、物質特定部４３が、座標系における３つ以上の物質の座標を結ぶ多角形を上記座標系内の領域として求め、対象物質の座標がこの多角形上に位置する場合、当該位置に基づいて３つ以上の物質の成分比を求めるように構成されていてもよい。
【００８７】
ここで、上記の例では、３つ以上の物質の座標を考慮する場合における上記座標系内の領域は、多角形には限定されない。この場合における上記座標系内の領域は、一般に、３つ以上の物質の座標に基づく図形、つまりこれら座標を考慮して形成された任意の図形でよい。たとえば、３つ以上の物質の座標を通過又は内包する図形を上記座標系内の領域として用いることができる。この例では、図形の外周線が直線である必要はなく、また物質の座標が当該図形の外周線上にある必要もない。対象物質の座標が図形上に位置する場合、当該位置に基づいて３つ以上の物質の成分比が算出される。
【００８８】
また、複数の物質のうちの１つが水である場合において、物質特定部４３が、第１及び第２の対象物質について決定された第１及び第２の座標のそれぞれと、水の座標とを結ぶ２つの線分を座標系内の領域として求め、これら線分の位置関係に基づいて第１の対象物質と第２の対象物質とが同種であるか否か特定するように構成されていてもよい。
【００８９】
また、物質特定部４３が、投影データに混入するノイズの標準偏差情報をあらかじめ記憶し、この標準偏差情報に基づいて座標系内の領域を含む２次元領域を求め、この２次元領域に対する対象物質の座標の位置に基づいてこの対象物質の種別を特定するように構成されていてもよい。
【００９０】
このようなＸ線ＣＴ装置によれば、対象物質の特性を所定の座標系における座標として表現し、この座標と複数の物質についての複数の座標とに基づいて対象物質の種別を特定することができる。したがって、ＣＴ値比が近い物質であっても座標を参照することにより高い精度での識別が可能である。
【００９１】
また、上記座標系内の領域を用いる構成を適用することにより、対象物質を構成する物質の成分比を求めることができる。更に、上記座標系内の領域を参照することにより、２つ以上の対象物質がそれぞれ同種であるか異種であるかを判別することができる。
【００９２】
また、ノイズの影響を考慮する構成を適用することにより、より高精度、高確度での物質の特定、成分比の特定、対象物質の分離などが可能となる。
【００９３】
〈第２の実施形態〉
この実施形態のＸ線ＣＴ装置の概略構成を図７に示す。このＸ線ＣＴ装置は、第１の実施形態で説明した構成を含んでいる。このＸ線ＣＴ装置は、ガントリ１と、前処理部２と、再構成処理部３と、データ処理部４と、制御部５と、表示部６と、操作部７とを有する。データ処理部４には、投影データ変換部４１と、画像形成部４２と、物質特定部４３とが設けられている。第１の実施形態と同じ構成部分は、特に言及しない限り第１の実施形態と同様の機能を有する。
【００９４】
基準物質の１つは水であるとする。また、物質特定部４３は、第１の実施形態で説明したように、対象物質の座標と水の座標とを結ぶ線分を求める（図２等を参照）。この線分は、対象物質と水との混合物における成分比、つまり対象物質の濃度を表すものである。
【００９５】
データ処理部４は、対象物（被検体等）の撮影領域の各位置について第１の実施形態の処理を適用することにより、撮影領域における対象物質の分布を求める。この分布情報は再構成処理部３に入力される。なお、分布情報には、撮影領域における対象物質の存在位置を示す情報と、濃度を示す情報とが含まれている。
【００９６】
分布情報と投影データを受けた再構成処理部３は、投影データにおける対象物質の分布領域に対して再構成処理を施すことで、対象物質の分布（たとえば濃度の分布）を表す分布画像を形成する。この処理の一例として、再構成処理部３は、分布情報に示す対象物質の存在位置に基づいて、この存在位置に対応する画素を特定する。そして、再構成処理部３は、投影データに基づいて、特定された画素についてのみ画像を再構成する。他の処理として、再構成処理部３は、投影データに基づいて通常の再構成処理を行い、それにより得られた画像のうち上記存在位置に対応する画素のみを抽出する。これら手法により得られる再構成画像が上記の分布画像である。
【００９７】
制御部５は、対象物質の濃度を表す線分上における位置の指定を受けて、この指定位置に応じた態様で分布画像を表示部６に表示させる。この処理について以下に説明する。
【００９８】
このＸ線ＣＴ装置には、対象物質の濃度を変更するためのユーザインターフェイスが設けられている。このユーザインターフェイスは、ソフトウェアに基づいて表示部６に表示されるものであってもよいし、ハードウェアであってもよい。ハードウェアとしてはダイヤルやスライドレバーなどがある。ソフトウェアを用いる場合、たとえば、ダイヤルやスライドレバーを模した画像を表示部６に表示させ、これを操作部７で操作するように構成できる。
【００９９】
ソフトウェアを用いる場合の具体例として、制御部５は、図２に示すグラフに基づいて、図８に示す画像を表示部６に表示させる。この画像は、図２に示す座標系と、線分画像Ｄと、スライダＥとを含む。座標系の座標軸には、物質名を示す文字列「造影剤」、「水」が付されている。線分画像Ｄは、線分Ｌ１を示している。スライダＥは、線分画像Ｄ上を移動可能とされている（図中の両側矢印を参照）。スライダＥは、たとえば操作部７のマウスによるドラッグ操作によって移動される。
【０１００】
ユーザは、スライダＥを所望の位置に移動させることにより、造影剤の濃度の指定を行う。より具体的に説明すると、線分画像Ｄと線分Ｌ１とはあらかじめ対応付けられており、制御部５は、スライダＥの位置に対応する線分Ｌ１の位置、つまり造影剤の濃度が指定される。制御部５は、指定された濃度に応じて分布画像の表示態様を変更する。表示態様の変更処理としては、画素値（輝度、色等）を変更するものや、模様を変更するものなどがある。
【０１０１】
分布画像が複数存在する場合には、各分布画像について濃度の指定を行うことが可能である。その場合、ユーザは、たとえばポインティングデバイスを用いて、所望の分布画像を指定する。制御部５は、指定された分布画像における所定の濃度に対応する位置にスライダＥを表示させる。ユーザは、スライダＥを所望の位置に移動させる。制御部５は、移動後のスライダＥの位置に応じて当該分布画像の表示態様を変更する。
【０１０２】
濃度の変更処理として、数値を入力する方法を採用することも可能である。その場合、制御部５は、対象物質の濃度情報を入力するための入力領域を表示部６に表示させる。ユーザは、操作部７（たとえばキーボード）を用いて入力領域に所望の濃度の値を入力する。入力領域は、たとえばプルダウンメニューのように、濃度の複数の選択肢から所望のものを選択可能とするオブジェクトであってもよい。
【０１０３】
初期値よりも濃度を高めた場合の分布画像を強調画像と呼ぶことがある。また、初期値よりも濃度を低下させた場合の分布画像を抑制画像と呼ぶことがある。
【０１０４】
強調画像及び抑制画像の例を説明する。図９Ａは、６０ｋｅＶの単色Ｘ線画像としての分布画像（原画像）Ｇ０を示す。なお、図９Ａ（画像を示す以下の図も同様）においては、実際の画像と表示輝度が逆転されている。つまり、実際の画像では対象物質の濃度が高いほど表示輝度が高くなっているが、図９Ａでは濃度が高いほど表示輝度が低くなっている。
【０１０５】
単色Ｘ線画像は次式（５）で定義される。
【０１０６】
【数５】
2013081765.tif 000007
【０１０７】
ただし：
ＣＴ  ｎｕｍｂｅｒは、ＣＴ値；
μは、式（５）に示す対象物質の線減弱係数；
μｗａｔｅｒは、水の線減弱係数；
をそれぞれ示す。
【０１０８】
原画像Ｇ０には、様々な濃度の造影剤の分布画像Ｈ０１、Ｈ０２、Ｈ０３、Ｈ０４、Ｈ０５が描写されている。分布画像Ｈ０１、Ｈ０２、Ｈ０３、Ｈ０４、Ｈ０５は、それぞれ、濃度２５［ｍｇＩ／ｍｌ］、２０［ｍｇＩ／ｍｌ］、１５［ｍｇＩ／ｍｌ］、１０［ｍｇＩ／ｍｌ］、５［ｍｇＩ／ｍｌ］に相当する。
【０１０９】
図９Ｂは、原画像Ｇ０に基づく強調画像Ｇ１を示す。強調画像Ｇ１には、原画像Ｇ０中の分布画像Ｈ０１、Ｈ０２、Ｈ０３、Ｈ０４、Ｈ０５にそれぞれ対応する分布画像Ｈ１１、Ｈ１２、Ｈ１３、Ｈ１４、Ｈ１５が描写されている。各分布画像Ｈ１ｉ（ｉ＝１〜５）は、分布画像Ｈ０ｉに示す濃度を高めた状態を示している。
【０１１０】
図９Ｃは、原画像Ｇ０に基づく抑制画像Ｇ２を示す。抑制画像Ｇ２には、原画像Ｇ０中の分布画像Ｈ０１、Ｈ０２、Ｈ０３、Ｈ０４、Ｈ０５にそれぞれ対応する分布画像Ｈ２１、Ｈ２２、Ｈ２３、Ｈ２４、Ｈ２５が描写されている。各分布画像Ｈ２ｉ（ｉ＝１〜５）は、分布画像Ｈ０ｉに示す濃度を低下させた状態を示している。
【０１１１】
このような強調画像Ｇ１及び抑制画像Ｇ２の生成方法を説明する。まず、制御部５は、造影剤の濃度を示す線分Ｌ１（図２参照）上に座標が位置する画素を特定する。次に、制御部５は、特定された各画素について、式（３）に示す造影剤の存在率ｃ１を１０００に設定し、水の存在率ｃ２を０に設定し、これらを式（２ａ）、（２ｂ）に代入することにより、強調画像Ｇ１を生成する。また、制御部５は、造影剤の存在率ｃ１を０に設定し、水の存在率ｃ２を１０００に設定し、これらを式（２ａ）、（２ｂ）に代入することにより、抑制画像Ｇ２を生成する。制御部５は、生成された強調画像Ｇ１や抑制画像Ｇ２を表示部６に表示させる。
【０１１２】
次に、濃度の算出方法の例を説明する。前述のように、座標系中の線分は、物質と水との混合物を表しており、線分上の位置が濃度に相当している。図２に示す場合において、造影剤の濃度を示す線分Ｌ１は、傾きが−１、切片が１の直線上にあるので、造影剤の濃度は次式で算出することができる。
【０１１３】
【数６】
2013081765.tif 000008
【０１１４】
図１０は、造影剤の濃度を示す線分Ｌ１上に座標が位置する係数ｃ１（ｘ，ｙ）、ｃ２（ｘ，ｙ）を抽出し、式（６）を用いて算出された濃度をカラー表示させた場合の表示態様の一例を模式的に表している。分布画像Ｇには５つの分布画像Ｈ１、Ｈ２、Ｈ３、Ｈ４、Ｈ５が描写されている。各分布画像Ｈｉ（ｉ＝１〜５）は、その濃度に応じた色で表示される。濃度と表示色との関係はカラーバーＪにより示されている。なお、複数の物質の濃度をカラー表示させる場合には、物質毎に異なる色で分布画像を表示させるように構成することも可能である。
【０１１５】
この実施形態のＸ線ＣＴ装置の作用及び効果を説明する。
【０１１６】
この実施形態のＸ線ＣＴ装置は、第１の実施形態の作用に加えて、次に示す作用を有する。なお、物質の１つは溶剤（たとえば水）である。物質特定部４３は、対象物質の座標と溶剤の座標とを結ぶ線分を求める。再構成処理部３は、投影データにおける対象物質の分布領域に対して再構成処理を施すことにより、対象物質の分布を表す分布画像を形成する。制御部５（表示画像生成部）は、線分上における位置の指定を受けて、当該指定位置に応じた態様の分布画像（表示画像）を生成して表示部６に表示させる。
【０１１７】
また、制御部５により線分を表す線分画像が表示部６に表示され、更に、この線分画像上の位置が操作部７を用いて指定されたことを受けて、制御部５が、当該指定位置に対応する線分上の位置を特定し、当該特定位置に応じた態様で分布画像を表示させるように構成することができる。
【０１１８】
また、物質の濃度情報を入力するための入力領域が表示部６に表示され、更に、操作部７を用いて対象物質の濃度情報が入力領域に入力されたことを受けて、制御部５（濃度算出部、表示画像生成部）が、当該濃度情報に対応する線分上の位置を特定し、当該特定位置に応じた態様の分布画像を生成して表示させるように構成することができる。
【０１１９】
また、制御部５が、線分上の指定位置に応じて分布画像の表示色を変更するように構成することが可能である。また、制御部５が、物質特定部４３により特定された対象物質の種別に応じて分布画像の表示色を変更するように構成することが可能である。更に、制御部５が、対象物質の種別と線分上の指定位置の双方に応じて分布画像の表示色を変更するように構成することも可能である。
【０１２０】
このようなＸ線ＣＴ装置によれば、第１の実施形態の効果に加えて、強調画像や抑制画像を形成し表示することができる。また、物質の濃度を視覚的に表現することができるので、濃度の直感的な理解が可能となる。
【０１２１】
なお、上記の例では、対象物質の座標と溶剤の座標とを結ぶ線分を求めて対象物質の分布画像を形成する場合について説明したが、溶剤以外の物質を適用して同様の処理を行うことも可能である。たとえば、前述の「複数の物質」のうちの一の物質の座標と対象物質の座標とを結ぶ線分を求め、更に、投影データにおける対象物質の分布領域に対して再構成処理を施すことにより分布画像を形成する。この分布画像は対象物質及び一の物質の分布を表すものである。つまり、２つの物質の分布を表現する単一の分布画像が得られる。
【０１２２】
〈第３の実施形態〉
第１の実施形態で説明した方法では、第１の物質の座標と水の座標とを結ぶ第１の線分と、第２の物質の座標と水の座標とを結ぶ第２の線分とが近い場合、これら物質を判別することは困難である。一方、非特許文献１に記載された方法では、第１の物質のＣＴ値と第２の物質のＣＴ値とが近い場合、これら物質を判別することは困難である。
【０１２３】
そこで、この実施形態では、これら２つの方法を併用することにより、より確実な物質の判別を可能とする技術を提供する。また、この実施形態では、特許文献１に記載された技術を利用して、実効原子番号画像、密度画像及び単色Ｘ線画像を形成する処理についても説明する。
【０１２４】
この実施形態のＸ線ＣＴ装置の概略構成を図１１に示す。このＸ線ＣＴ装置は、第１の実施形態で説明した構成を含んでいる。また、図示は省略するが、このＸ線ＣＴ装置は、第２の実施形態で説明した構成を更に含んでいてもよい。
【０１２５】
このＸ線ＣＴ装置は、ガントリ１と、前処理部２と、再構成処理部３と、データ処理部４と、画素値比算出部８と、物質特定部９と、特定画像形成部１０と、特定制御部１１とを有する。データ処理部４には、投影データ変換部４１と、画像形成部４２と、物質特定部４３とが設けられている。なお、物質特定部９と特定制御部１１はそれぞれ「特定部」の一例として機能する。
【０１２６】
第１の実施形態と同じ構成部分は、特に言及しない限り第１の実施形態と同様の機能を有する。なお、第１の実施形態で説明した物質特定処理はデータ処理部４が実行し、非特許文献１に記載された物質特定処理は画素値比算出部８及び物質特定部９が実行する。
【０１２７】
再構成処理部３は、デュアル・エナジー・ＣＴにより得られた第１及び第２の投影データのそれぞれに対して再構成処理を施すことにより第１及び第２の画像を形成する。
【０１２８】
画素値比算出部８は、非特許文献１に記載されているように、第１及び第２の画像における画素値（ＣＴ値）の比を算出する。この処理は、第１の画像と第２の画像との間で画素の対応付けを行い、対応付けられた２つの画素のＣＴ値の比を算出するものである。
【０１２９】
物質特定部９は、非特許文献１に記載された方法を用いることで、画素値比算出部８により算出されたＣＴ値比に基づいて、対象物質の種別を特定する。
【０１３０】
以下、具体例として、尿酸、軟骨及び軟組織が混在する状況において尿酸を特定する場合を説明する。これら物質について得られたＣＴ値比の分布の例を図１２に示す。本例では、管電圧８０ｋＶ及び１３５ｋＶでデュアル・エナジー・ＣＴによる撮影を行うものとする。図１２では、管電圧８０ｋＶでのＣＴ値をｘ軸に、管電圧１３５ｋＶでのＣＴ値をｙ軸にとっている。したがって、原点と任意の座標とを通過する直線の傾きがＣＴ値比に相当する。
【０１３１】
分布領域Ｋ１、Ｋ２、Ｋ３は、それぞれ尿酸、軟骨、軟組織についての分布を示している。分布領域Ｋ１、Ｋ２、Ｋ３の位置関係から分かるように、尿酸と軟骨とはＣＴ値比が異なっており互いの判別が可能であるが、尿酸と軟組織ではＣＴ値比が近いため判別が困難である。なお、直線Ｍ１、Ｍ２は、それぞれ尿酸のＣＴ値比の最大値と最小値を表している。
【０１３２】
一方、図１３は、第１の実施形態で説明した方法で得られた、尿酸、軟骨及び軟組織の座標の分布の例を示している。なお、縦軸及び横軸は、それぞれ造影剤及び水に基づき設定されている。分布領域Ｎ１、Ｎ２、Ｎ３は、それぞれ尿酸、軟骨、軟組織についての分布を示している。直線Ｌ６は、尿酸の座標の分布状態を示す直線である。また、直線Ｌ６ａ、Ｌ６ｂは、それぞれ、尿酸の分布領域の広がりに応じた直線の分布範囲の上限及び下限を示している（第１の実施形態を参照）。軟骨及び軟組織については直線の図示を省略しているが、その形態から直線の向きは容易に把握できる。分布領域Ｎ１、Ｎ２、Ｎ３の位置から分かるように、本例では、尿酸と軟骨では近い直線が得られるので物質相互の判別が困難であり、尿酸と軟組織では直線の向きが十分異なるので判別が可能である。
【０１３３】
以上の考察の下、第１の実施形態で説明した方法と非特許文献１に記載された方法とを併用する処理の一例を説明する（図１４を参照）。まず、第１の実施形態と同様に、デュアル・エナジー・ＣＴにより投影データを取得する（Ｓ２１）。
【０１３４】
再構成処理部３，画素値比算出部８及び物質特定部９は、前述の要領で、図１２に示すＣＴ値比の分布情報、つまり分布領域Ｋ１、Ｋ２、Ｋ３を得る（Ｓ２２）。また、データ処理部４は、前述の要領で、図１３に示す座標の分布情報、つまり分布領域Ｎ１、Ｎ２、Ｎ３を得る（Ｓ２３）。これら分布情報は特定制御部１１に送られる。以下、ＣＴ値比の分布情報を第１の分布情報と呼び、座標の分布情報を第２の分布情報と呼ぶことがある。
【０１３５】
特定制御部１１は、第１及び第２の分布情報に基づいて３つの物質をそれぞれ特定する。その処理の例として、特定制御部１１は、まず、図１２に示す直線Ｍ１、Ｍ２と、図１３に示す直線Ｌ６ａ、Ｌ６ｂを求める。なお、直線Ｍ１、Ｍ２や直線Ｌ６ａ、Ｌ６ｂは、今回の測定により得られたデータから作成したものであっても、事前に得られた測定データや規定の値（理論値や標準的な値等）であってもよい。また、特定制御部１１は、各物質について、そのＣＴ値比の分布領域に基づく直線や、その座標の分布領域に基づく直線を求めてもよい。このようにして、直線Ｍ１、Ｍ２により挟まれた特定領域と、直線Ｌ６ａ、Ｌ６ｂにより挟まれた特定領域とが得られる（Ｓ２４、Ｓ２５）。
【０１３６】
続いて、特定制御部１１は、軟骨及び軟組織のそれぞれのＣＴ値比の分布領域Ｋ２、Ｋ３が、直線Ｍ１、Ｍ２により挟まれた特定領域に含まれているか判断する（Ｓ２６）。この処理は、分布領域Ｋ２、Ｋ３自体が当該特定領域に含まれているか判断するものでもよいし、上記で直線を求めた場合にはその直線が当該特定領域に含まれているか判断するものでもよい。
【０１３７】
なお、前者の場合、分布領域Ｋ２、Ｋ３の全体が当該特定領域に含まれているか判断するようにしてもよいし、分布領域Ｋ２、Ｋ３の所定割合が当該特定領域に含まれているか判断するようにしてもよい。
【０１３８】
また、後者の場合、直線の長さによっては少しの傾きの違いでも直線が当該特定領域からはみ出すことが考えられる。その場合、直線の長さの範囲をあらかじめ決めておき、その範囲内において直線が当該特定領域に含まれているときに、この直線は当該特定領域に含まれていると判断するようにしてもよい。
【０１３９】
続いて、特定制御部１１は、軟骨及び軟組織のそれぞれの座標の分布領域Ｎ２、Ｎ３が、直線Ｌ６ａ、Ｌ６ｂにより挟まれた特定領域に含まれているか判断する（Ｓ２７）。この処理は、ＣＴ値比の場合と同様に実行できる。
【０１４０】
次に、特定制御部１１は、ＣＴ値比の分布領域Ｋ２、Ｋ３が特定領域に含まれているかの判断結果と、座標の分布領域Ｎ２、Ｎ３が特定領域に含まれているか判断結果とに基づいて、尿酸、軟骨及び軟組織の判別を行う。
【０１４１】
その具体例として、特定制御部１１は、ＣＴ値比の分布領域が特定領域に含まれており（Ｓ２８：Ｙｅｓ）、かつ、座標の分布領域が特定領域に含まれている（Ｓ２９：Ｙｅｓ）と判断された場合、これら分布領域に対応する物質は尿酸であると判断する（Ｓ３０）。
【０１４２】
また、特定制御部１１は、ＣＴ値比の分布領域が特定領域に含まれており（Ｓ２８：Ｙｅｓ）、かつ、座標の分布領域が特定領域に含まれていない（Ｓ２９：Ｎｏ）と判断された場合、これら分布領域に対応する物質は軟組織であると判断する（Ｓ３１）。
【０１４３】
また、特定制御部１１は、ＣＴ値比の分布領域が特定領域に含まれておらず（Ｓ２８：Ｎｏ）、かつ、座標の分布領域が特定領域に含まれている（Ｓ３２：Ｙｅｓ）と判断された場合、これら分布領域に対応する物質は軟骨であると判断する（Ｓ３３）。
【０１４４】
また、特定制御部１１は、ＣＴ値比の分布領域が特定領域に含まれておらず（Ｓ２８：Ｎｏ）、かつ、座標の分布領域が特定領域に含まれていない（Ｓ３２：Ｎｏ）と判断された場合、これら分布領域に対応する物質は尿酸でも軟骨でも軟組織でもないと判断する。以上で、この処理例の説明は終了となる。
【０１４５】
次に、特定画像形成部１０について説明する。特定画像形成部１０は、実効原子番号画像、密度画像、及び単色Ｘ線画像のうちの少なくとも１つ（特定画像と呼ぶことがある）を形成する。この処理は、第１の実施形態で説明した線形結合（式（３））における係数ｃ１及びｃ２に基づいて形成される、各基準物質の分布を表す画像（基準物質画像）を、特許文献１に記載された方法で組み合わせることにより実行されるものである。
【０１４６】
なお、実効原子番号画像は、対象物における実効原子番号の分布を表す画像である。また、密度画像は、対象物質における物質の密度の分布を表す画像である。また、単色Ｘ線画像は、単一のエネルギーを有するＸ線で対象物質をスキャンした場合を模擬的に再現する画像である。これら画像を形成可能な特定画像形成部１０は、「形成部」の一例として機能する。
【０１４７】
この実施形態の物質特定部４３は、第１の実施形態で説明した処理で得られた複数の対象物質のそれぞれについての特定結果と、これら対象物質のそれぞれについて特定画像形成部１０により形成された実効原子番号画像、密度画像、及び／又は単色Ｘ線画像に基づいて、これら対象物質の種別をそれぞれ特定することが可能である。なお、第１の実施形態で説明した処理は、複数の物質についてあらかじめ得られた座標系における複数の座標と、対象物質について物質特定部４３により決定された座標とに基づいて、対象物質の種別を特定するものである。
【０１４８】
この実施形態のＸ線ＣＴ装置によれば、第１の実施形態の効果を得ることができるとともに、第１の実施形態と同様の座標系における各物質の座標の分布状態と、各物質のＣＴ値比の分布状態の双方に基づいて、物質の種別を特定することができる。したがって、対象物に含まれる物質の特定精度の更なる向上を図ることが可能である。
【０１４９】
また、この実施形態のＸ線ＣＴ装置によれば、実効原子番号画像、密度画像、及び／又は単色Ｘ線画像を形成することができるので、物質の分布状態を視覚的に把握することが可能である。また、座標系を用いた物質特定処理とこれら画像を用いた物質特定処理を併用することにより、対象物に含まれる物質の特定精度の更なる向上を図ることが可能である。
【０１５０】
〈変形例〉
上記実施形態の変形例を説明する。
【０１５１】
第１の変形例を説明する。対象物質に未知の成分が含まれている場合、図５に示す線分Ｌｉのような直線的な相関関係が得られず、領域Ｂ等を適当に決定し難いことがある。たとえば図１５Ａに示すように、３つの対象物質について、３つの座標Ｖ１、Ｖ２、Ｖ３が得られたものとする。これら座標Ｖ１、Ｖ２、Ｖ３は、同一の直線上に位置しない。その場合、この座標系と、座標Ｖ１、Ｖ２、Ｖ３を表示部に表示させる。ユーザは、表示された座標系及び座標Ｖ１、Ｖ２、Ｖ３を参照し、操作部を用いて当該座標系内に領域を任意に設定する。このとき、たとえば図１５Ｂに示すように、座標Ｖ１、Ｖ２、Ｖ３を含むように領域Ｗが設定される。
【０１５２】
領域Ｗは、対象物の未知の成分などに起因する誤差を反映した造影剤の濃度の座標の分布範囲として用いられる。領域Ｗにおいて座標（ｃ１，ｃ２）＝（１，０）と座標（ｃ１，ｃ２）＝（０，１）とを結ぶ線分（図５の線分Ｌ１に相当する）上にない座標に対しては、当該線分に基づいて濃度の値が対応付けられる。この対応付けは、たとえば図５に示す場合と同様にして行われる。物質特定部４３は、たとえば図５に示す場合と同様の手法により、決定された対象物質の座標と、領域Ｗとの位置関係に基づいて、この対象物質の種別を特定する。
【０１５３】
この変形例によれば、対象物質に未知の成分が含まれている等の誤差要因が介在する場合であっても、ユーザの判断によって物質の特定を高い精度で行うことができる。
【０１５４】
第２の変形例を説明する。上記実施形態では、エネルギーが異なる２つのＸ線を用いるデュアル・エナジー・ＣＴの手法が適用されているが、実施形態はこれに限定されるものではない。
【０１５５】
たとえば、エネルギーが異なる３以上のＸ線を用いてスキャンを行うことが可能なＸ線ＣＴ装置においては、生成部は、エネルギーが異なるＮ（３以上）のＸ線でそれぞれ対象物をスキャンしてＮ個の投影データを生成する。変換部は、これらＮ個の投影データを、Ｎ個の基準物質に対応するＮ個の新たな投影データに変換する。画像形成部は、変換部で変換したＮ個の新たな投影データをそれぞれ再構成することで、Ｎ個の基準物質に対応するＮ個の基準物質画像を形成する。特定部は、Ｎ個の基準物質画像の画素の値の相関に基づいて対象物質を特定する。
【０１５６】
この場合に実行される処理においては、上記実施形態では２つである線形結合の基底の個数がＮ個となり、上記実施形態では２次元である座標系の次元がＮ次元となる。つまり、エネルギーが異なるＸ線の数が２である場合と３以上である場合との間に、実質的な相違はない。よって、上記実施形態は、エネルギーが異なる複数（２以上）のＸ線を用いる場合として、次のように一般化できる。
【０１５７】
実施形態に係るＸ線ＣＴ装置は、対象物をスキャンして得られた投影データに基づいて、対象物内の画像を表示する装置であって、生成部と、変換部と、画像形成部と、特定部とを有する。生成部は、エネルギーが異なるＸ線でそれぞれ対象物をスキャンして複数の投影データを生成する。変換部は、複数の投影データを、複数の基準物質に対応する複数の新たな投影データに変換する。画像形成部は、変換部で変換した複数の新たな投影データをそれぞれ再構成することで、複数の基準物質に対応する複数の基準物質画像を形成する。特定部は、複数の基準物質画像の画素の値の相関に基づいて対象物質を特定する。
【０１５８】
このようなＸ線ＣＴ装置によれば、上記実施形態と同様に、対象物に含まれる物質を高い精度で特定することが可能である。更に、エネルギーが異なる３以上のＸ線を用いることで、２つのＸ線を用いる上記実施形態よりも高い精度での物質特定が可能となる。
【０１５９】
上記実施形態で説明した任意の構成や処理についても同様の一般化を行うことができる。その場合、上記実施形態と同様の効果が得られるとともに、物質特定の精度の更なる向上を図ることが可能となる。
【０１６０】
〈物質特定方法〉
実施形態に係る物質特定方法は、たとえば上記実施形態に係るＸ線ＣＴ装置により実行される。実施形態の物質特定方法は、生成ステップと、変換ステップと、画像形成ステップと、特定ステップを含む。生成ステップでは、エネルギーが異なるＸ線でそれぞれ対象物をスキャンして複数の投影データを生成する。変換ステップでは、複数の投影データを、複数の基準物質に対応する複数の新たな投影データに変換する。画像形成ステップでは、変換ステップで変換した複数の新たな投影データをそれぞれ再構成することで、複数の基準物質に対応する複数の基準物質画像を形成する。特定ステップでは、複数の基準物質画像の画素の値の相関に基づいて対象物質を特定する。
【０１６１】
特定ステップは、複数の基準物質画像における対応する画素の値が、予め設定された相関を有するかを判定することにより、対象物質を特定するものであってよい。
【０１６２】
変換ステップは、基準物質に対応するＸ線減弱係数を含む計算式を用いて、複数の投影データを、複数の基準物質に対応する複数の新たな投影データに変換するものであってよい。
【０１６３】
特定ステップにより特定された対象物質の領域を識別可能にした対象物内の画像を生成する表示画像生成ステップを更に含んでいてもよい。
【０１６４】
特定ステップにより特定された対象物質の領域の画素を、他の領域の画素に比べて、強調または抑制した表示画像を生成する表示画像生成ステップを更に含んでいてもよい。
【０１６５】
予め記憶された計算式に基づいて、複数の基準物質画像における対応する画素の値から、対象物質の濃度情報を求める濃度算出ステップを更に含んでいてもよい。
【０１６６】
濃度情報を表した表示画像を生成する表示画像生成ステップを更に含んでいてもよい。
【０１６７】
画像形成ステップにおいて、複数の投影データを再構成して得られた複数のエネルギー画像を生成し、かつ、特定ステップにおいて、複数のエネルギー画像の複数の基準物質画像の画素の値の相関に基づいて対象物質を特定してもよい。
【０１６８】
基準物質画像に基づいて、対象物における実効原子番号の分布を表す実効原子番号画像、対象物質における物質の密度の分布を表す密度画像、及び、単一のエネルギーを有するＸ線で対象物質をスキャンした場合を模擬的に再現する単色Ｘ線画像のうちの少なくとも１つを形成する形成ステップを更に含み、かつ、特定ステップにおいて、形成ステップで形成された画像に基づいて、対象物質の特定を行ってもよい。
【０１６９】
このような物質特定方法によれば、第１の実施形態と同様に、ＣＴ値比が近い物質であっても座標を参照することにより高い精度での識別が可能である。また、対象物質を構成する物質の成分比を求めることや、２つ以上の対象物質がそれぞれ同種であるか異種であるかを判別することができる。また、第１の実施形態のようにノイズの影響を考慮した物質特定方法を適用することにより、より高精度、高確度での物質の特定、成分比の特定、対象物質の分離などが可能となる。
【０１７０】
また、実施形態に係る物質特定方法によれば、第２の実施形態と同様に、上記効果に加えて、強調画像や抑制画像を形成し表示することができる。また、物質の濃度を視覚的に表現することができるので、濃度の直感的な理解が可能となる。
【０１７１】
また、実施形態に係る物質特定方法によれば、第３の実施形態と同様に、上記効果に加えて、各物質の座標の分布状態と各物質のＣＴ値比の分布状態の双方に基づいて、物質の種別を特定することにより、対象物に含まれる物質の特定精度の更なる向上を図ることが可能である。また、実効原子番号画像、密度画像、及び／又は単色Ｘ線画像を形成することができるので、物質の分布状態を視覚的に把握することが可能である。また、座標系を用いた物質特定処理とこれら画像を用いた物質特定処理を併用することにより、対象物に含まれる物質の特定精度の更なる向上を図ることが可能である。
【０１７２】
〈画像処理装置〉
実施形態に係る画像処理装置について説明する。実施形態に係る画像処理装置の構成例を図１６に示す。画像処理装置１００は、院内ＬＡＮ等のネットワークを介して画像保管装置３００に接続されている。画像保管装置３００は、Ｘ線ＣＴ装置２００により形成された画像を記憶している。
【０１７３】
Ｘ線ＣＴ装置２００は、たとえば第１の実施形態のＸ線ＣＴ装置（図１参照）から物質特定部４３を除いた構成を有する。つまり、Ｘ線ＣＴ装置２００は、対象物をスキャンして得られた投影データに基づいて、前記対象物内の画像を表示する装置であって、次の機能を有する：エネルギーが異なるＸ線でそれぞれ対象物をスキャンして複数の投影データを生成する生成機能；複数の投影データを、複数の基準物質に対応する複数の新たな投影データに変換する変換機能；変換機能で変換した複数の新たな投影データをそれぞれ再構成することで、複数の基準物質に対応する複数の基準物質画像を形成する画像形成機能。Ｘ線ＣＴ装置２００は、形成された複数の基準物質画像を、院内ＬＡＮ等のネットワークを介して画像保管装置３００に送る。画像保管装置３００は、たとえば、画像保存通信システム（Ｐｉｃｔｕｒｅ  Ａｒｃｈｉｖｉｎｇ  ａｎｄ  Ｃｏｍｍｕｎｉｃａｔｉｏｎ  Ｓｙｓｔｅｍｓ；ＰＡＣＳ）である。
【０１７４】
画像処理装置１００は、Ｘ線ＣＴ装置２００により形成されて画像保管装置３００に記憶された画像を処理する。なお、他のルートを経由して画像処理装置１００に画像を入力することも可能である。他のルートの例として、画像が記憶されたＤＶＤ等の記録媒体がある。
【０１７５】
画像処理装置１００は、外部から画像を取得する処理を行う画像取得部１０１を有する。画像保管装置３００に記憶された画像を取得する場合、画像取得部１０１は、ネットワークを介して画像保管装置３００と通信する通信部を含む。更に、画像取得部１０１は、たとえば、取得する画像を選択するためのユーザインターフェイス（表示部、操作部等）と、選択された画像を取得するための信号を生成し、通信部を制御してこの信号を画像保管装置３００に送信させる通信制御部を含む。記録媒体を介して画像を取得する場合、画像取得部１０１は、この記録媒体に記録された情報を読み取り可能な読み取り部（ドライブ装置等）を含む。
【０１７６】
画像取得部１０１により取得された画像は、記憶部１０２に記憶される。記憶部１０２は、ハードディスクドライブ等の記憶装置を含んで構成される。この実施形態では、Ｘ線ＣＴ装置２００が上記の要領で形成した画像が記憶部１０２に記憶される。この画像は、エネルギーが異なるＸ線でそれぞれ対象物をスキャンして生成された複数の投影データを、複数の基準物質に基づき変換して得られた複数の新たな投影データを、それぞれ再構成することで形成された、複数の基準物質に対応する複数の基準物質画像である。
【０１７７】
物質特定部１０３は、記憶部１０２に記憶された複数の基準物質画像の画素の値の相関に基づいて対象物質を特定する。この処理は、たとえば第１の実施形態の物質特定部４３と同様にして実行される。物質特定部１０３は「特定部」の一例として機能する。
【０１７８】
画像処理装置１００が表示部を有する場合、物質特定部１０３による処理結果が表示部に表示される。また、画像処理装置１００がたとえばネットワーク上のサーバである場合、画像処理装置１００は、物質特定部１０３による処理結果を、ネットワークを介して所定のユーザ端末に送信する機能（送信部）を有する。また、画像処理装置１００は、物質特定部１０３による処理結果を記録媒体に記録する記録部（ドライブ装置等）を有していてもよい。
【０１７９】
このような画像処理装置１００によれば、対象物に含まれる物質を高い精度で特定することが可能である。なお、画像処理装置１００は、上記実施形態に係るＸ線ＣＴ装置の任意の機能を有していてもよい。以下のいずれかの構成が適用される場合、上記実施形態において説明した、当該構成に対応する作用及び効果が奏される。
【０１８０】
たとえば、物質特定部１０３は、複数の基準物質画像における対応する画素の値が、予め設定された相関を有するかを判定することにより、対象物質を特定するものであってよい。
【０１８１】
また、画像処理装置１００は、物質特定部１０３により特定された対象物質の領域を識別可能にした対象物内の画像を生成する表示画像生成部を備えていてもよい。
【０１８２】
また、画像処理装置１００は、物質特定部１０３により特定された対象物質の領域の画素を、他の領域の画素に比べて、強調または抑制した表示画像を生成する表示画像生成部を備えていてもよい。
【０１８３】
また、画像処理装置１００は、予め記憶された計算式に基づいて、複数の基準物質画像における対応する画素の値から、対象物質の濃度情報を求める濃度算出部を備えていてもよい。これに加え、画像処理装置１００は、濃度情報を表した表示画像を生成する表示画像生成部を備えていてもよい。
【０１８４】
また、物質特定部１０３は、複数のエネルギー画像の複数の基準物質画像の画素の値の相関に基づいて対象物質を特定するものであってよい。
【０１８５】
また、画像処理装置１００は、基準物質画像に基づいて、対象物における実効原子番号の分布を表す実効原子番号画像、対象物質における物質の密度の分布を表す密度画像、及び、単一のエネルギーを有するＸ線で対象物質をスキャンした場合を模擬的に再現する単色Ｘ線画像のうちの少なくとも１つを形成する形成部を更に備え、かつ、物質特定部１０３が、形成部により形成された画像に基づいて対象物質の特定を行うように構成されていてもよい。
【０１８６】
〈プログラム及び記憶媒体〉
上記実施形態で説明した処理を、Ｘ線ＣＴ装置又はこれに含まれるコンピュータに実行させるためのプログラムを構成することができる。また、上記実施形態で説明した処理を、画像処理装置（コンピュータ）に実行させるためのプログラムを、たとえばＤＶＤ等の記憶媒体に記憶させることができる。また、インターネットやＬＡＮ等のネットワークを介して、これらプログラムのいずれかを送信するシステムを構築することができる。
【０１８７】
この発明の実施形態を説明したが、上記の実施形態は例として提示したものであり、発明の範囲を限定することを意図していない。これら新規な実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれるとともに、特許請求の範囲に記載された発明とその均等の範囲に含まれる。
【符号の説明】
【０１８８】
１  ガントリ
２  前処理部
３  再構成処理部
４  データ処理部
４１  投影データ変換部
４２  画像形成部
４３  物質特定部
５  制御部
６  表示部
７  操作部
８  画素値比算出部
９  物質特定部
１０  特定画像形成部
１１  特定制御部
【図１】
2013081765.tif 000009
【図２】
2013081765.tif 000010
【図３】
2013081765.tif 000011
【図４】
2013081765.tif 000012
【図５】
2013081765.tif 000013
【図６】
2013081765.tif 000014
【図７】
2013081765.tif 000015
【図８】
2013081765.tif 000016
【図９Ａ】
2013081765.tif 000017
【図９Ｂ】
2013081765.tif 000018
【図９Ｃ】
2013081765.tif 000019
【図１０】
2013081765.tif 000020
【図１１】
2013081765.tif 000021
【図１２】
2013081765.tif 000022
【図１３】
2013081765.tif 000023
【図１４】
2013081765.tif 000024
【図１５Ａ】
2013081765.tif 000025
【図１５Ｂ】
2013081765.tif 000026
【図１６】
2013081765.tif 000027