JP2018139754A - 画像処理装置、ｘ線ｃｔ装置及び画像処理方法 - Google Patents

Abstract

【課題】
演算時間を増大させることなく複数種類の基底物質の物理量を算出する。
【解決手段】
エネルギー分布の異なる２種類以上の投影データに基づいて、第１の基底物質群に含まれる２以上の基底物質の物理量をそれぞれ算出する第１の基底物質変換部と、前記第１の基底物質群に含まれる２以上の基底物質の物理量から、被写体の投影像又は再構成像の少なくとも一方であって、かつ、複数の画像を生成する画像生成部と、複数の前記画像に基づいて、前記第１の基底物質群とは異なる第２の基底物質群に含まれる基底物質の物理量を算出する第２の基底物質変換部と、を備える画像処理装置を提供する。
【選択図】図２

Description

本願発明は、画像処理装置、Ｘ線ＣＴ装置及び画像処理方法に関し、特に、フォトンカウンティング型のＸ線検出器により取得されたデータに対する画像処理を行う画像処理装置、Ｘ線ＣＴ装置及び画像処理方法に関する。

Ｘ線ＣＴ装置は、複数の方向から撮影した被検体のＸ線透過像である投影データからＸ線吸収係数（線減弱計数）を算出し、被検体の断層像である再構成像を得る装置である。従来、投影データを取得するためのＸ線検出器として、積分型のＸ線検出器を用いたＸ線ＣＴ装置が主流であったが、近年、Ｘ線フォトンの個数を計測するフォトンカウンティング型のＸ線検出器を搭載したＸ線ＣＴ装置が開発されつつある（例えば、特許文献１）。この装置では、これまでのＸ線ＣＴ装置では取得できなかったエネルギー毎の仮想単色画像や、基底物質密度画像、実効原子番号画像、電子密度画像、光電効果画像、コンプトン散乱画像など、吸収係数以外の分布を示す再構成像、撮影で用いたスペクトル以外のエネルギー分布のＸ線における吸収係数像など、様々な再構成像（以降、これらの画像を「マルチエネルギー画像」という）を生成することができる。
このようなフォトンカウンティング型Ｘ線ＣＴ装置では、複数のエネルギービン毎に投影データを取得し、変換マップなどの変換パラメータを用いて、投影データの出力値に対応する基底物質の物理量を得ることができる。ここで基底物質の物理量とは、投影データの場合、基底物質の長さや面密度である。基底物質の物理量を用いることで、様々なマルチエネルギー画像を生成することができる。

特開２０１４−２３９８４０号公報

しかしながら、特許文献１のＸ線ＣＴ装置では、変換パラメータが、Ｘ線照射条件毎に必要であり、更に基底物質の組み合わせ（セット）毎に必要となる。従って、撮影前に、想定するＸ線照射条件と基底物質の組み合わせに対する多数の変換パラメータを、実験やシミュレーションを用いて決定しておく必要があり、多くの工程を要する。そして、これら工程で得られた多くの変換パラメータを記憶しておくための大きな記憶容量が必要となる。複数種類の基底物質を設定して夫々画像を生成する場合には、さらに多くの工程と大きな記憶容量を要してしまう。

ところで、基底物質の物理量の変換には、投影データを基に行う方法と、再構成像を基に行う方法とがある。投影データを基に行う方法は、特にＸ線検出素子毎の特性の違いを考慮して変換を行うことができるため、精度良く基底物質の物理量を求めることができるという利点がある。
しかしながら、投影データを基に行う方法によって基底物質の物理量変換を行うと、処理時間を要してしまうため、様々な基底物質を設定して画像を生成するような処理には不向きであった。
一方、再構成像を基に行う方法は、Ｘ線検出素子毎の特性の違いを考慮することができないものの、高速に処理を行うことができるという利点がある。

本発明は上記実情に鑑みてなされたものであり、演算時間を増大させることなく複数種類の基底物質の物理量を算出し、延いては所望の画像を取得することを目的とする。

上記課題を解決するために、本発明は以下の手段を提供する。
本発明の一態様は、エネルギー分布の異なる２種類以上の投影データに基づいて、第１の基底物質群に含まれる２以上の基底物質の物理量をそれぞれ算出する第１の基底物質変換部と、前記第１の基底物質群に含まれる２以上の基底物質の物理量から、被写体の投影像又は再構成像の少なくとも一方であって、かつ、複数の画像を生成する画像生成部と、複数の前記画像に基づいて、前記第１の基底物質群とは異なる第２の基底物質群に含まれる基底物質の物理量を算出する第２の基底物質変換部と、を備える画像処理装置を提供する。

本発明によれば、演算時間を増大させることなく複数種類の基底物質の物理量を算出し、所望の画像を取得することができる。

本発明の実施形態に係る画像処理装置を適用したＸ線ＣＴ装置の概略を示すブロック図である。 本発明の実施形態に係る画像処理装置を適用したＸ線ＣＴ装置における演算部の概略を示すブロック図である。 図１のＸ線ＣＴ装置の検出部の概略を示す参考図である。 図１のＸ線ＣＴ装置の検出部によって検出するフォトン数とエネルギー範囲との関係を示すラフである。 図１のＸ線ＣＴ装置の検出部において生じた電荷に応じた電圧信号を示すグラフである。 本発明の実施形態に係る画像処理装置を適用したＸ線ＣＴ装置における基底物質演算処理の概略を説明するフローチャートである。 本発明の実施形態に係る画像処理装置を適用したＸ線ＣＴ装置における基底物質演算処理を説明する説明図である。 （Ａ）〜（Ｅ）本発明の実施形態に係る画像処理装置を適用したＸ線ＣＴ装置の基底物質演算処理において用いる変換マップの一例である。 （Ａ）〜（Ｃ）本発明の実施形態に係る画像処理装置を適用したＸ線ＣＴ装置の基底物質演算処理において用いる変換マップの一例である。 本発明の実施形態に係る画像処理装置を適用したＸ線ＣＴ装置の表示部に表示されるユーザインターフェース画面の一例を示す参考図である。 本発明の実施形態の変形例１に係る画像処理装置を適用したＸ線ＣＴ装置における基底物質演算処理を説明する説明図である。 本発明の実施形態の変形例２に係る画像処理装置を適用したＸ線ＣＴ装置における基底物質演算処理を説明する説明図である。 本発明に係る画像処理装置のブロック図である。

本発明の一実施形態について、図面を参照して説明する。
本発明に係る画像処理装置は、エネルギー分布の異なる２種類以上の投影データに基づいて、第１の基底物質群に含まれる２以上の基底物質の物理量をそれぞれ算出する第１の基底物質変換部と、第１の基底物質群に含まれる２以上の基底物質の物理量から、被写体の投影像又は再構成像の少なくとも一方であって、かつ、複数の画像を生成する画像生成部と、複数の前記画像に基づいて、第１の基底物質群とは異なる第２の基底物質群に含まれる基底物質の物理量を算出する第２の基底物質変換部と、を備える画像処理装置である。
このような画像処理装置によれば、第１の基底物質群に含まれる２以上の基底物質の物理量に基づいて取得した画像から第２の基底物質群に含まれる基底物質の物理量を算出するので、画像演算時間を増大させることなく複数種類の基底物質の物理量を算出することができる。

以下、より具体的に本発明の実施形態について説明する。
＜実施形態＞
以下、本発明の実施形態に係る画像処理装置を備えたＸ線ＣＴ装置について図面を参照して説明する。

図１に示すように、Ｘ線ＣＴ装置は、撮影系としての、Ｘ線源１００と、Ｘ線検出器１１１と、これらＸ線源１００及び検出器１１１の検出部１０４（後述）を対向配置し所定の回転軸を中心に回転するガントリー回転部１０１と、ガントリー回転部１０１の開口内に配置された寝台天板１０３と、これら撮影系の動作に伴いＸ線検出器１１１が取得した信号を処理する信号処理部１１２とを備えている。

Ｘ線源１００は、例えば管電圧で加速した電子ビームをタングステンやモリブデンなどのターゲット金属に衝突させ、その衝突位置（焦点）からＸ線を発生させる。
ガントリー回転部１０１は、Ｘ線源１００及び検出部１０４を互いに対向配置し、所定の回転軸を中心に回転する。ガントリー回転部１０１の中央には、被検体１０２が挿入される開口が設けられ、この開口内に、被検体１０２が寝かせられる寝台天板１０３が配置されている。寝台天板１０３とガントリー回転部１０１とは、所定の方向に相対的に移動可能となっている。

Ｘ線検出器１１１は、入射したＸ線フォトンを検出し、複数（例えば、３つ）のエネルギー範囲に分別して計数を行うフォトンカウンティング方式のＸ線検出素子が複数配置された検出部１０４と、Ｘ線検出素子から出力される投影像を収集する信号収集部１０８とを備えている。検出部１０４の詳細は後述する。

信号処理部１１２は、演算部１０５、表示部１０６、制御部１０７、記憶部１０９及び入力部１１０を備えている。
図２に示すように、演算部１０５は、収集した信号に所定の演算処理を行うため、信号収集部１０８で収集した信号に対して、基底物質変換に用いる変換パラメータを決定する変換パラメータ決定部１５０、第１の基底物質群に含まれる複数の基底物質の物理量を求める第１の基底物質変換部１５１、第１の基底物質群に属する基底物質の物理量から投影像を生成する投影像生成部（画像生成部）１５２、再構成処理を行いマルチエネルギー画像等の再構成像を生成する再構成像生成部１５３、投影像生成部１５２において生成された投影像又は再構成像生成部１５３において生成された再構成像から第２の基底物質群に属する基底物質の物理量を算出する第２の基底物質変換部１５４、及びこれら各部を制御する主制御部１５９を備えている。

演算部１０５は、記憶部１０９に記憶された情報や変換パラメータに基づいて種々の演算を行う。演算部１０５は必要に応じて記憶部１０９から各種情報やパラメータ等を読み出し、信号処理、画像処理、及び画像再構成等の演算を行う。
表示部１０６は、ユーザによる各種入力を実現するユーザインターフェース画面（詳細は後述する）を表示したり、演算部１０５により生成されたマルチエネルギー画像等を表示したりする。

制御部１０７は、Ｘ線源１００の発生駆動源の動作を制御するＸ線制御部、Ｘ線検出器１１１の信号読み出し動作を制御する読み出し制御部、ガントリー回転部１０１の回転と寝台天板１０３の移動を制御する撮影制御部、及びこれら各部全体を制御する全体制御部を備えている。
記憶部１０９は、演算部１０５における演算処理に用いられるパラメータやデータ、投影データを取得した際の撮影条件等を記憶している。具体的には、基底物質変換に用いられる撮影条件情報１４０や、基底物質の物理量を求めるための変換パラメータ１４１等が含まれる。基底物質の物理量を算出するための変換パラメータとして、例えば、変換マップ等を記憶している。
入力部１１０は、Ｘ線ＣＴ装置における撮影条件等の入力を行う。

演算部１０５及び制御部１０７は、その一部又は全部をＣＰＵ（中央処理装置）、メモリ及び主記憶部１０９を含むシステムとして構築することができ、演算部１０５及び制御部１０７を構成する各部の機能は、予め記憶部１０９に格納されたプログラムをＣＰＵがメモリにロードし、実行することにより実現することができる。また機能の一部または全部は、ＡＳＩＣ（Application Specific Integrated Circuit）やＦＰＧＡ（Field Programmable Gate Array）などのハードウェアで構成することも可能である。

以下の説明においては、特に説明しない限り、上述した撮影系、制御部１０７及び信号処理部１１２を構成する要素は、公知のＸ線ＣＴ装置が備える要素と同様の構成を有し、同様の機能を持つ。

続いて、Ｘ線検出器１１１の検出部１０４について説明する。
検出部１０４は、Ｘ線源１００を略中心とした円弧状に複数配置されており、ガントリー回転部１０１の回転に伴い、Ｘ線源１００との位置関係を保ちながら回転する。なお、図１においては、説明の便宜上、検出部１０４が８個の場合を示したが、実際の装置では、例えば４０個程度である。また検出部１０４の前面にはＸ線グリッド（図示せず）が設置されており、Ｘ線源１００から照射されたＸ線のうち、被検体３００などで散乱されたＸ線が、検出部１０４に入射するのを防ぐ。

検出部１０４は、例えば複数の同一のサイズのＸ線検出素子が、チャネル方向とスライス方向に２次元的に配置され、それぞれの方向で等間隔に配置された構造となっている。Ｘ線検出素子は、図１のＸ線ＣＴ装置において、チャネル方向と回転方向１１３を、スライス方向と回転軸方向１１４とを略一致させて配置されている。

一方、検出部１０４の断面は、例えば、図３に示すように、検出層４０１を挟むように正負の電極４０２，４０３が設けられてＸ線検出素子が構成され、その電極４０２，４０３には、読み出し回路４０５が接続された構造を有する。本実施形態では、正の電極４０２は、各Ｘ線検出素子で共通の構造である。

検出層４０１は、例えばＣｄＴｅ（テルル化カドミウム）、ＣｄＺｎＴｅ（カドミジンクテルル）、Ｓｉ（シリコン）などの半導体材料から成る。Ｘ線は、図２中の矢印４０４で示すように、正の電極４０２側から検出層４０１に入射し、Ｘ線フォトンを検出してそのエネルギーに応じた量の電荷を生じる。そして、検出部１０４は、この電荷に応じたアナログの電気信号を信号収集部１０８に出力する。

検出部１０４からのアナログ信号は、信号収集部１０８の読出し回路４０５に入力される。読出し回路４０５は、入射をトリガーとして読み出し、エネルギー分別とデジタル変換を行う。エネルギー分別は、入射した電荷で生じた電気信号を、所定の閾値により複数のエネルギー範囲に分別する処理である。このとき、発生した電気信号の波高や発生量は、入射したＸ線フォトンのエネルギーに依存するため、Ｘ線フォトンに応じたエネルギー範囲に分別できる。このようにエネルギー範囲に分別してそれぞれでカウントすることで、デジタル信号を得る。信号収集部１０８は、このように得たエネルギー範囲毎のデジタル信号を演算部１０５に出力する。

ここで、読出し回路４０５で行われる分別方法の一例を説明する。ここでは、３つのエネルギー範囲、つまり、所定の低エネルギー閾値未満であるエネルギー範囲（以下、「低エネルギー範囲」という）、所定の低エネルギー閾値以上高エネルギー閾値未満のエネルギー範囲（以下、「中エネルギー範囲」という）、所定の高エネルギー閾値以上のエネルギー範囲（以下、「高エネルギー範囲」という）に判別する例について説明する。

３つのエネルギー範囲について、図４を用いて説明する。図４に示すグラフにおいて、横軸はエネルギーを、縦軸はフォトン数を示している。図４中の曲線１３３は、検出部１０４のＸ線検出素子に入射するＸ線のスペクトルの一例を示す。閾値１３０〜１３２は、Ｘ線フォトンの入力の有無を判断する閾値であり、詳細は後述する。領域１３６−１が低エネルギー範囲、領域１３７が中エネルギー範囲、領域１３８が高エネルギー範囲である。ただし、詳細は後述するが、閾値１３０未満のときは、読み出し回路４０５は信号読み出しを行わないため、実質的には領域１３６−２が低エネルギー範囲となる。入射したＸ線フォトンは、Ｘ線検出素子毎に、読み出し回路４０５にて、エネルギーに応じたエネルギー範囲に分別されて、エネルギー範囲毎に計数されていく。

図５は、発生した電荷により生じた電圧信号１２７を示すグラフであり、横軸は時間、縦軸は電圧を示す。また、閾値１３１が低エネルギー閾値を、閾値１３２が高エネルギー閾値を表す。図５の例では、サンプリング時間１２３中にＸ線が入射してパルス出力１２０を生じ、サンプリング時間１２４中にＸ線が入射してパルス出力１２１を生じ、サンプリング時間１２５中にＸ線が入射してパルス出力１２２を生じている。なお、図５では、サンプリングはＸ線が入射するタイミングだけでなく、Ｘ線が入射しない場合（サンプリング時間１２６）でも周期的に行われる場合を示したが、Ｘ線フォトンが入射したタイミングでサンプリングが行われる場合も在り得る。

本実施形態においては、読み出し回路４０５は、入射信号によるＸ線フォトンをエネルギー範囲に分けて分別する前に、Ｘ線フォトンが入射したかを判断する。このため、読み出し回路４０５は、サンプリング毎に、その区間における出力電圧１２７の最大値と、閾値１３０と比較する。ここで閾値１３０は、Ｘ線フォトンの入力の有無を判断する閾値であり、出力電圧１２０が閾値１３０以上の場合には、以下に説明するエネルギー範囲の分別を行い、出力電圧１２０は閾値１３０未満の場合は行わない。このような閾値１３０は、電圧１２７は、Ｘ線が入力しないときにも検出部１０４の回路ノイズによって変動しており、これをＸ線による信号と誤検出しないために必要となる。そのため閾値１３０は、ゼロより大きく低エネルギー閾値１３１より小さな値で設定される。

図５におけるサンプリング時間１２６はＸ線が入射していないが、雑音によって変動が生じた出力電圧１２７の一例であり、出力電圧１２０は低エネルギー閾値１２１未満のため、読み出し回路４０５は、信号のエネルギー範囲の分別を行わない。

Ｘ線が入射して閾値１３０以上の電圧１２７が生じた場合、読み出し回路４０５は、その信号をエネルギー範囲に分別し、それぞれのエネルギー範囲で入射Ｘ線フォトン数をカウントする。この分別は、例えば、出力電圧１２７の最大値を、低エネルギー閾値１３１と高エネルギー閾値１３２と比較することで行う。

例えば、サンプリング時間１２３では、出力電圧１２０は低エネルギー閾値１３１未満のため、低エネルギー範囲と分別する。サンプリング時間１２４では、出力電圧１２１は低エネルギー閾値１３１以上であって、高エネルギー閾値１３２未満のため、中エネルギー範囲と分別する。サンプリング時間１２５では、出力電圧１２０は高エネルギー閾値１３２以上のため、高エネルギー範囲と分別する。

エネルギー範囲が分別されて入射したＸ線フォトン数をエネルギー範囲毎にカウントし、この合計をビュー毎に出力する。サンプリング時間はビューの時間と比較して非常に短く、ビュー間に多数回サンプリングを行うことになる。
このように、入射の有無とエネルギー範囲の分別を行い、信号収集部１０８は、ビュー毎に、エネルギー範囲毎のデジタル信号を生成する。
なおサンプリングでの最大値を用いて分別を行う代わりに、例えば、サンプリング中の出力電圧の積分値を用いてもよく、分別手法は上記手法に限定されない。

このように構成されたＸ線ＣＴ装置では、一般に、以下のように撮影動作が行われる。
まず撮影者が、入力部１１０から撮影条件を入力して実撮影の開始を入力すると、制御部１０７はＸ線源１００からのＸ線の照射と、ガントリー回転部１０１を制御し撮影を開始する。例えば、Ｘ線源１００では、１２０ｋＶの管電圧で電子ビームを加速して、寝台天板１０３に載った被検体３００に向けて照射されＸ線が照射される。被検体３００を透過したＸ線は検出部１０４で検出される。検出部１０４は、Ｘ線検出素子毎に入射Ｘ線のエネルギーに応じた電荷を発生する。信号収集部１０８は、上述の通り、この電荷を高エネルギー範囲、中エネルギー範囲、低エネルギー範囲に分別し、エネルギー範囲毎及びビュー毎にデジタル信号を得て演算部１０５に出力する。

制御部１０７は、ガントリー回転部１０１を回転方向に回転させることで、被検体３００に対するＸ線の照射角度を変化させる。このように回転駆動させながら、ビュー毎に焦点位置を変更させて撮影を繰り返し行い、３６０度分のデジタル信号を取得する。撮影は、例えば０．４度ごとに複数ビューにわたって行われる。このような撮影により、３６０度分のデジタル信号が得られ、３６０度分のデジタル信号を投影データとする。なお、Ｘ線源１００から発生されるＸ線は、各ビューに同期したパルスＸ線でも良いし、連続Ｘ線でも良い。

演算部１０５は、取得した投影データに対し、所定の補正処理や演算処理を行い、マルチエネルギー画像を生成する。すなわち、画像処理部としての演算部１０５において、第１の基底物質変換、再構成像生成、第２の基底物質変換、マルチエネルギー画像生成等を行う。

以下、演算部１０５による演算処理について、図６のフローチャートに従って説明する。
本実施形態においては、第１の基底物質変換処理において、第１の基底物質群に属する基底物質Ａ及び基底物質Ｂの物理量を算出し、第２の基底物質変換処理においては、第１の基底物質群とは異なる第２の基底物質群に属する基底物質Ｃ及び基底物質Ｄの物理量を算出する例について説明する。また、第１の基底物質群に属する基底物質Ａ及び基底物質Ｂの物理量及び第２の基底物質群に属する基底物質Ｃ及び基底物質Ｄの物理量として面密度を算出する例について説明する。

なお、第１の基底物質群と第２の基底物質群とは互いに異なり、第１の基底物質群と第２の基底物質群の何れか一方にのみ含まれる基底物質が存在する。一方で、第１の基底物質群と第２の基底物質群の夫々に含まれる基底物質数が同数であっても異なる数であってもよく、また、基底物質の一部が共通していてもよい。
また、以下の説明において、理解の容易のため便宜的に、信号収集部１０８から得られる信号を「投影データ」といい、第１の基底物質変換の結果や第２の基底物質変換の結果を用いて生成した投影データを「投影像」という。

図６に示すように、ステップＳ１１において、主制御部１５９がエネルギー分布の異なる２種類以上の投影データを取得する。ステップＳ１２では、第１の基底物質変換部１５１が、ステップＳ１１で取得した投影データに基づいて第１の基底物質変換を行う。ここで、第１の基底物質変換とは、エネルギー分布の異なる２種類以上の投影データから第１の基底物質群に属する基底物質Ａ及び基底物質Ｂの物理量を算出することをいう。

次のステップＳ１３では、投影像生成部１５２又は再構成像生成部１５３が、ステップＳ１２で算出した第１の基底物質群に属する基底物質Ａ及び基底物質Ｂの物理量から投影像又は再構成像の少なくとも一つを生成する。ステップＳ１４では、主制御部１９がステップＳ１３で生成した画像を記憶部１０９に記憶させると共に、必要に応じて表示部１０７に表示させる。

次のステップＳ１５では、主制御部１５９が、ステップＳ１２で算出した基底物質を変更するか否かの判定を行う。つまり、第２の基底物質変換を行うか否かの判定を行う。例えば、ユーザが基底物質の変更を指示した場合、第２の基底物質変換を行うと判定する。基底物質の変更を行うと判定された場合はステップＳ１６に進み、第２の基底物質変換部１５４が、ユーザにより入力された第２の基底物質群に属する基底物質Ｂ及び基底物質Ｃの物理量を、ステップＳ１３で生成された投影像又は再構成像に基づいて算出する。

ステップＳ１７では、第２の基底物質群に属する基底物質Ｂ及び基底物質Ｃの物理量を用いて、再構成像生成部１５３が、マルチエネルギー画像を生成し、ステップＳ１８で生成したマルチエネルギー画像を表示部１０７に表示させて上記処理を終了する。また、ステップＳ１５において基底物質の変更を行わないと判定された場合は、上記処理を終了する。

上述した演算部１０５による基底物質変換に係る演算処理の詳細及び各処理によって生成されるデータ又は画像について、図７に示す説明図を用いてより詳細に説明する。図７中、実線は処理の内容を示し、破線は実線で示す処理を実行することによって生成されるデータ又は画像を示す。

図７に示すように、ステップＳ１００からステップＳ１０６までの処理は、第１の基底物質の物理量を算出してマルチエネルギー画像を取得する処理であり、ステップＳ１０７からＳ１１０までの処理は、第２の基底物質の物理量を算出して、新たにマルチエネルギー画像を取得する処理である。

ステップＳ１００では、ユーザが、入力部１１０によりＸ線ＣＴ装置に対して、撮影準備として種々の設定の入力と共に撮影開始指示の入力を行い、撮影を開始し、演算部１０５が信号収集部１０８から投影データ２００を取得する。撮影準備としては、例えばＸ線源１００のＸ線管電圧やＸ線管電流、Ｘ線照射範囲、焦点サイズ、Ｘ線フィルタの種類や有無、ボータイフィルタの種類や有無、撮影時間、スキャン時間、撮影範囲などの撮影条件の設定や、基底物質の種類、基底物質の物理量、生成する再構成像のサイズや種類（例えば、マルチエネルギー画像として、仮装単色画像を選択する等）などの設定が含まれる。
設定した撮影条件は、記憶部１０９に記憶される。以下の説明においては、撮影条件情報として、特に、Ｘ線管電圧を記憶部１０９に記憶させた場合について説明する。

ステップＳ１０１では、変換パラメータ決定部１５０が、信号収集部１０８から取得した投影データ２００に対して第１の基底物質変換を行う際に適用する変換パラメータを決定する。ここで、投影データは、低エネルギー範囲、中エネルギー範囲、及び高エネルギー範囲の３つのエネルギー範囲のそれぞれで取得する。基底物質変換処理に適用する変換パラメータは、予め本撮影の前に計測して生成して記憶部１０９に保存しておき、ステップＳ１０１では、記憶部１０９に保存した撮影条件から、適用する変換パラメータを決定する。

基底物質変換処理に適用する変換パラメータとしては、例えば、図８（Ａ）〜（Ｃ）に示す変換用マップを用いることができる。図８（Ａ）に示す変換マップは、基底物質Ａの面密度値（横軸）と基底物質Ｂの面密度値（縦軸）が決まった時に得られる低エネルギー範囲の出力値を表すものであり、曲線１６１はある出力値の取りうる基底物質Ａの面密度値と基底物質Ｂの面密度値との組み合わせの集合を示す曲線（以下、「等高線」という）の一例を表す。これは同じ低エネルギー範囲の出力値に対して、複数の基底物質Ａと基底物質Ｂの面密度値の組み合わせがあることを意味する。

同様に、図８（Ｂ）に示す変換マップは、基底物質Ａの面密度値（横軸）と基底物質Ｂの面密度値（縦軸）が決まった時に得られる中エネルギー範囲の出力値を表すものであり、曲線１５５はある出力値の等高線の一例を表す。同様に、図８（Ｃ）に示す変換マップは、基底物質Ａの面密度値（横軸）と基底物質Ｂの面密度値（縦軸）が決まった時に得られる高エネルギー範囲の出力値を表すものであり、曲線１５６はある出力値の等高線の一例を表す。

このような変換用マップのセットは、Ｘ線照射状況によってそれぞれ異なる。すなわち、撮影条件として、Ｘ線管電圧８０ｋＶp及び１２０ｋＶｐが設定可能である場合、それぞれのＸ線管電圧で異なる変換用マップのセットが必要となる。従って、撮影条件毎、すなわちＸ線管電圧毎の変換マップのセットを変換パラメータとして記憶部１０９に保存しておく。Ｓ１０１にて、変換パラメータ決定部１５０が、撮影条件、ここでは、Ｘ線管電圧に基づいて変換用マップの１つのセットを、各Ｘ線検出素子、各ビューにおいてそれぞれ決定する。

このように第１の基底物質変換処理に適用する変換パラメータを決定した後、次のステップＳ１０２では、第１の基底物質変換部１５１により、低エネルギー範囲、中エネルギー範囲及び高エネルギー範囲の３つのエネルギー範囲における投影データと、ステップＳ１０１で決定した第１の基底物質変換処理に適用する変換パラメータとしての変換マップのセットを用いて、第１の基底物質群に含まれる各基底物質の物理量を算出し、第１の基底物質群の基底物質画像を生成する。本実施形態では、物理量として面密度を算出する。面密度の算出は各Ｘ線検出素子、各ビューでそれぞれ行い、第１の基底物質群に含まれる各基底物質の面密度から成る基底物質画像（投影像２０１）をそれぞれ取得する。すなわち、基底物質画像として、基底物質毎の投影像、本実施形態では基底物質Ａと基底物質Ｂの投影像２０１をそれぞれ生成する。

この面密度値の決定方法について、図８（Ｄ）を用いて説明する。
図８（Ｄ）は、図８（Ａ）〜（Ｃ）を重ねて表示したものであり、１つのエネルギー範囲のデータでは複数の第１の基底物質の面密度値の組み合わせがあるが、複数のエネルギー範囲を用いることで、それらの交点が全てのエネルギー範囲の結果を満たす値となり、最適な第１の基底物質群の面密度値となる。図８（Ｄ）中、曲線１６１、１５５、１５６がそれぞれのエネルギー範囲で計測された結果である場合、図８（Ｄ）の交点１５７が最適な値となり、交点１５７の横軸の値が最適な基底物質Ａの面密度値、縦軸の値が最適な基底物質Ｂの面密度値となる。

ただし、Ｘ線の線量や量子数の揺らぎなどにより、投影データの出力値は誤差を有するため、変換マップから一意に解を得られない場合がある。図８（Ｅ）に、このような一意に解を得られない場合の一例を示している。図８（Ｅ）に示すように、曲線１５５と曲線１５６の交点が点１５７−１、曲線１６１と曲線１５５の交点が点１５７−２、曲線１６１と曲線１５６の交点が点１５７−３となり、解は一意とならない。この場合、最も可能性が高い解を推定する必要がある。

解を推定するために、例えば、シミュレーションにより交点１５７−１、１５７−２、１５７−３の付近で基底物質Ａと基底物質Ｂのそれぞれの面密度の値を変化させて、再構成像を生成した際に最もアーチファクトが少なくなる点を選択することができる。このような解の推定を行うための推定部（図示せず）を、第１の基底物質変換部１０５３が具備するように構成することもできる。ただし、このような処理は演算量が多く、投影データに対して行う場合、全Ｘ線検出素子、全ビューで行うことになるため、多くの時間が必要となる場合がある。

次にステップＳ１０３では、投影像生成部１５２により、基底物質画像２０１に基づいてマルチエネルギー画像の投影像２０２を生成する。この投影像２０２は、ステップＳ１００で設定したマルチエネルギー画像（以下、「指定マルチエネルギー画像」という）の投影像と、第２の基底物質変換で使用するマルチエネルギー画像（以下、「第２の基底物質変換用マルチエネルギー画像」という）の投影像とを含んでいる。

ここで、第２の基底物質変換用マルチエネルギー画像は、撮影に先立って予め決定されているものであり、第２の基底物質変換に必要な複数の画像から成る。本実施形態では、指定マルチエネルギー画像として、撮影で使用したものと異なるＸ線管電圧８０ｋＶpでの画像であり、第２の基底物質変換用マルチエネルギー画像として５０ｋｅＶと１００ｋｅＶの仮想単色画像とした場合の例について説明する。

ここで、第２の基底物質変換用マルチエネルギー画像の仮想単色画像の投影像２０２の生成方法の一例について説明する。
単色のエネルギーをＥ_０、投影像をＰ_Ｅ０とし、２つの基底物質を基底物質Ａと基底物質Ｂとしてそれらの投影像の値をδ_１とδ_２、エネルギーＥ_０における質量吸収係数をμｍ_１（Ｅ_０）とμｍ_２（Ｅ_０）、エネルギーＥ_０におけるフォトン数Ｓ（Ｅ_０）とすると、基底物質演算部１０５５にて算出した２つの基底物質の面密度投影像の値δ_１、δ_２を用いて、投影像Ｐ_Ｅ０は、例えば、式（１）を用いることで算出することができる。

エネルギーＥ_０に５０ｋｅＶの場合と１００ｋｅＶの場合を求めることで、それぞれの仮想単色画像の投影像２０２が生成できる。

続いて、指定マルチエネルギー画像として、Ｘ線管電圧８０ｋＶｐのときの画像の生成方法の一例について説明する。
Ｘ線管電圧８０ｋＶｐである場合は、例えば、フォトン数Ｓ（Ｅ_０）として、Ｘ線管電圧８０ｋＶｐのときに照射されるＸ線スペクトルのフォトン数を用い、式（１）のエネルギーＥ_０を０ｋｅＶから８０ｋｅＶまでの仮想単色画像を求め、それらを積分することで仮想単色画像の投影像を得ることができる。ここでは、例えば、１ｋｅＶ毎に仮想単色画像を求めて、式（２）のように加算することで、Ｘ線管電圧８０ｋＶｐの際の投影像Ｐ_{８０ｋＶｐ}を生成できる。

指定マルチエネルギー画像として、は、仮想単色画像、基底物質密度画像、実効原子番号画像、電子密度画像、光電効果画像、及びコンプトン散乱画像など、吸収係数以外の分布を示す再構成像などの様々な画像を生成することができる。

また、本実施形態では、第２の基底物質変換用マルチエネルギー画像として、５０ｋｅＶと１００ｋｅＶの仮想単色画像について説明したが、Ｘ線管電圧の組み合わせはこれに限定されるものではなく、様々なエネルギーの仮想単色画像の組み合わせを適用することができる。２種類のエネルギーの組み合わせだけでなく、３種以上の組み合わせとすることもできる。更に、異なる複数のＸ線管電圧の画像の場合でも良い。また、エネルギー範囲毎の再構成像、すなわち、低エネルギー範囲、中エネルギー範囲、高エネルギー範囲の一部または全てのエネルギー範囲に対応する再構成像のセットでも良い。更に実効原子番号画像、電子密度画像、光電効果画像、コンプトン散乱画像のような様々なマルチエネルギー画像のセットであっても良い。

次にステップＳ１０４では、再構成像生成部１５３が、ステップＳ１０３で求めた投影像２０２に対して再構成処理を行い、再構成像を生成する。再構成像生成部１５３は、まず指定マルチエネルギー画像の投影像を再構成して指定マルチエネルギー画像を生成し、ステップＳ１０５にて表示する。次に、第２の基底物質変換用マルチエネルギー画像の投影像を再構成して、第２の基底物質変換用マルチエネルギー画像を生成し、ステップＳ１０６にて記憶部１０９に保存する。

再構成処理には、例えば、フェルドカンプ画像再構成法などの様々な技術を適用することができる。

なお、ステップＳ１０５において、演算部１０５は、表示部１０７に指定マルチエネルギー画像のみを表示することもできる他、指定マルチエネルギー画像と共に第２の基底物質変換用マルチエネルギー画像を表示することもできる。また、ステップＳ１０６では、第２の基底物質変換用マルチエネルギー画像の他に、指定マルチエネルギー画像を保存することもできる。

次にステップＳ１０７において、第２の基底物質変換部１５４は、入力部１１０を介して、ユーザによる基底物質の変更の有無の入力を受け付ける。第２の基底物質変換部１５４は、ユーザによる基底物質の変更を行わない旨の入力を受け付けた場合に、処理を終了する。一方、第２の基底物質変換部１５４は、基底物質の変更を行う旨の入力を受け付けた場合に、次のステップＳ１０８に進む。

ステップＳ１０８では、第２の基底物質変換部１５４が、算出する基底物質や生成条件の入力を、入力部１１０を介して受け付ける。ここで、ユーザにより入力部１１０を介して入力された基底物質は第２の基底物質群に含まれる基底物質（仮に、基底物質Ｃ及び基底物質Ｄとする）であり、第１の基底物質群に含まれる基底物質Ａ及び基底物質Ｂとは異なるものである。また、生成条件には、例えば生成するマルチエネルギー画像のサイズや種類などが含まれる。

ステップＳ１０８における基底物質や生成条件の入力が終了すると、ステップＳ１０９に進み、第２の基底物質変換部１５４が、記憶部１０９に記憶された第２の基底物質変換用マルチエネルギー画像を読み出し、第２の基底物質の物理量を算出する。
以下、第２の基底物質群に含まれる基底物質Ｃ及び基底物質Ｄの物理量の算出手法の一例について、図９を用いて説明する。本実施形態では、第２の基底物質群に含まれる基底物質Ｃ及び基底物質Ｄの密度を算出し、第２の基底物質変換用マルチエネルギー画像として、エネルギー５０ｋｅＶとエネルギー１００ｋｅＶの仮想単色画像を生成する例について説明する。

基底物質Ｃ及び基底物質Ｄの密度の算出には、例えば、第１の基底物質の面密度を用いた場合と同様に変換マップを用いる。この変換マップは、それぞれの第２の基底物質変換用マルチエネルギー画像毎に、予め記憶部１０９に記憶される。本実施例では、図９（Ａ）の変換マップがエネルギー５０ｋｅＶの仮想単色画像用であり、図９（Ｂ）の変換マップがエネルギー１００ｋｅＶの仮想単色画像用の場合である。図９（Ａ）の横軸は基底物質Ｃの密度値、縦軸は基底物質Ｄの密度値を示す。本第２の基底物質変換用マルチエネルギー画像は、５０ｋｅＶと想定する入射Ｘ線のスペクトルが決まっているので、基底物質Ｃと基底物質Ｄの質量吸収係数が求まり、基底物質Ｃの密度値と基底物質Ｄの密度値との組み合わせが決まる毎に再構成像の画素値（以降、再構成像の画素をボクセル、その画素値をボクセル値と記す）が得られる。図９（Ａ）及び（Ｃ）中、曲線１５８は、あるボクセル値の取りうる基底物質Ｃの密度値と基底物質Ｄの密度値の組み合わせの集合を示す曲線(以下、「等高線」という)である。

図９（Ｂ）に示すエネルギー１００ｋｅＶの仮想単色画像の変換マップにおいて、曲線１５９は、あるボクセル値の取りうる基底物質Ｃの密度値と基底物質Ｄの密度値の組み合わせを示す曲線（等高線）である。図９（Ｃ）は、図９（Ａ）と図９（Ｂ）を重ねたものであり、曲線１５８と曲線１５９の交点１６０が、最適な基底物質Ｃと基底物質Ｄの密度値となる。

上述の例では、エネルギー５０ｋｅＶ及びエネルギー１００ｋｅＶの２つの再構成像を用いて第２の基底物質変換を行う例について説明したが、３つ以上の再構成像を用いてもよい。このとき、図８で示すように、一意に解が得られない場合には、例えば、シミュレーションにより解を推定することが好ましい。再構成像のボクセル数は、投影データに比べて少ないため、投影データから推定する場合に比して短時間で処理することができる。
ステップＳ１０９では、このような処理を再構成像の全ボクセルで行い、密度値の集合である密度画像２０４を生成する。この密度画像２０４は、基底物質Ｃと基底物質Ｄの双方について生成する。

次のステップＳ１１０では、再構成像生成部１５３により、密度画像２０４からマルチエネルギー画像の生成を行う。以下、マルチエネルギー画像の一例として、仮想単色画像を生成する例について説明する。
エネルギーＥにおける吸収係数（線減弱係数）μ（Ｅ）を、再構成像の全ボクセルについて算出する。吸収係数μ（Ｅ）は、以下の式（３）に従って算出する。

ただし、μｍ_ｎ（Ｅ）は、基底物質の質量吸収係数、ｃ_ｎは基底物質の密度である。基底物質の密度ｃ_ｎは、密度画像２０４の値となり、質量吸収係数μｍ_ｎ（Ｅ）は、例えば文献値を用いる。また、ｎは基底物質の種類を表し、基底物質Ｃのときｎ＝３、基底物質Ｄのときはｎ＝４である。

以上のように吸収係数を、再構成像の全ボクセルで求め、更にＣＴ値への変換することで、マルチエネルギー画像（再構成像）２０５として仮想単色画像を生成することができる。次のステップＳ１１１では、生成したマルチエネルギー画像２０５を表示する。

この後、ステップＳ１０７に戻り、第２の基底物質変換部１５４は、入力部１１０を介して、ユーザによる基底物質の変更の有無の入力を受け付ける。第２の基底物質変換部１５４は、ユーザによる基底物質の変更を行わない旨の入力を受け付けた場合に、処理を終了する。一方、第２の基底物質変換部１５４は、基底物質の変更を行う旨の入力を受け付けた場合に、ステップＳ１０８以降の処理を繰り返し、再度マルチエネルギー画像２０５を生成して表示する。このような処理により、第２の基底物質を変更して、それぞれのマルチエネルギー画像を繰り返し取得することができる。

本実施形態によれば、第１の基底物質変換に際して、第１の基底物質群に含まれる基底物質の組み合わせが取り得るＸ線照射条件に応じた変換用パラメータを予め記憶部１０９に記憶しておくことで基底物質の物理量を算出することができる。
そして、第２の基底物質変換では、既に得られた投影データもしくは再構成像を用いて第２の基底物質群に含まれる基底物質について変換、すなわち基底物質の物理量を算出するので、投影データ及び再構成像はＸ線照射条件に依存せず、１つのＸ線照射条件に応じた変換用パラメータを記憶部１０９に記憶しておくのみで基底物質変換処理を行うことができる。

すなわち、１回の基底物質変換で予め定めた物理量を求める場合、基底物質群として取り得る基底物質の組み合わせに応じた変換用パラメータが、取り得るＸ線照射条件に応じた分だけ必要となり、膨大な変換用パラメータが必要となる。これに比して、上述した本実施形態によれば、基底物質を変換するための変換パラメータを予め用意する工数を低減することができるのは勿論、膨大な変換パラメータを記憶させる必要がない。

また、基底物質の変換処理を行う都度被写体の撮影を行う必要がなく、演算部１０５における計算コストを抑制することができる。
さらに、投影データから基底物質の物理量を算出するには処理時間を要してしまうが、再構成像から基底物質の物理量を算出する場合には、再構成像のボクセル数が投影データに比して少ないため高速に処理することができる。

このように本実施形態によれば、当初算出した第１の基底物質群から変更して、第２の基底物質群に含まれる基底物質の物理量を算出する場合であっても、撮影のし直しや投影データからの多くの演算処理を行う必要がなく、短時間で所望の基底物質の物理量を示すマルチエネルギー画像を取得することができる。

なお、上述した実施形態においては、基底物質を変更するか（第２の基底物質変換を行うか）否かに拘わらず、第２の基底物質変換に用いるマルチエネルギー画像を生成する例について説明したが、第２の基底物質変換に用いるマルチエネルギー画像の生成のタイミングは適宜決定することができる。例えば、基底物質の変更指示の入力後に、第２の基底物質変換用マルチエネルギー画像を生成し、これを用いて第２の基底物質変換処理及びマルチエネルギー画像の生成を行い、変更された基底物質のマルチエネルギー画像を生成することもできる。

（表示部１０６におけるユーザインターフェースの例）
以下、上述した実施形態に係るＸ線ＣＴ装置の表示部１０６に表示されるユーザインターフェース画面の一例について説明する。
図１０に、表示部１０６に表示するユーザインターフェース画面の一例を示す。図１０のユーザインターフェース画面９００は、撮影開始入力ボタン９１１、撮影中止ボタン９１２、撮影条件設定領域９１３、基底物質設定領域９０１、画像種別入力ボタン９０７、基底物質変更ボタン９０３、終了ボタン９０２及び画像表示領域９０４，９０５，９０６を備えている。

撮影条件設定領域９１３は、撮影の前に、撮影条件を設定し入力するための領域である。図１０の例では、Ｘ線照射条件として、Ｘ線管電圧、Ｘ線管電流、Ｘ線照射範囲（照射スライス幅）、焦点サイズ、Ｘ線フィルタ（種類と有無）、ボータイフィルタ（種類と有無）、撮影時間が表示されている。
撮影開始入力ボタン９１１は、ユーザにより押下されることにより撮影開始指示を受け付ける。撮影開始以降、撮影を中止する場合は、撮影中止入力部９１２を押下することにより撮影の停止指示を受け付ける。

画像表示領域９０４，９０５，９０６には、例えば、第１の基底物質から作成したマルチエネルギー画像が画像表示部９０４に表示され、第１の基底物質から作成したマルチエネルギー画像が画像表示部９０５に表示される。

基底物質変更ボタン９０３は、基底物質を変更するか否かの入力を受け付ける。基底物質変更ボタン９０３がユーザにより押下されると、ユーザは続いて基底物質設定領域９０１にて第２の基底物質を設定し、画像種別入力ボタン９０７を介してマルチエネルギー画像の種類を設定する。図１０の例では、マルチエネルギー画像の種類として４０ｋｅＶの仮想単色画像を選択した例を示す。基底物質設定領域９０１では、２セット以上の組み合わせの基底物質を設定可能であり、その中の１つのセットを設定する。図１０の例では、基底物質として水と骨を選択した例を示している。
基底物質を変更しない場合は、終了ボタン９０２を押下する。

次に、基底物質変換の処理開始指示部９０３にて処理開始を指示すると、演算部１０５により、再設定された第２の基底物質を算出して得られたマルチエネルギー画像を生成し、画像表示領域９０６に表示される。

繰り返し第２の基底物質を変更して、マルチエネルギー画像を生成する場合は、上記の設定及び入力を繰り返し行うことで、何度で画像の生成及び表示を行うことができる。
生成された画像は、画像表示領域９０４，９０５，９０６に表示する他、別の画像表示領域を設けて表示することもできる。撮影や基底物質の変換に必要な情報を確認しながら入力できると共に、基底物質の種類や数などが互いに異なる画像を、一画面中に同時に表示することで複数の画像を比較して両者の違いを容易に把握することができるというメリットがある。

上述の例では、第１の基底物質群が２つの基底物質（水とカルシウム）からなる１セットを含み、第２の基底物質群が２つの基底物質（水と骨）からなる１セットを含む場合について説明したが、各基底物質群に含まれる基底物質の組み合わせは適宜変更することができる。ただし、第１の基底物質群に含まれる基底物質数が多いと、それに応じて必要となる変換用パラメータ数が増大するため、第１の基底物質群に含まれる基底物質数は、第２の基底物質群に含まれる基底物質数よりも少ないことが好ましい。

上述の例では、第２の基底物質変換用マルチエネルギー画像を保存する例について説明したが、これに限られず、例えば、基底物質変換用マルチエネルギー画像の投影像を保存することもできる。具体的には、第２の基底物質変換用マルチエネルギー画像の投影像２０２を保存しておき、再構成像作成Ｓ１０４を行わない状態で基底物質の変更の有無Ｓ１０７の入力を待つようにし、基底物質の変更の有無Ｓ１０７で「有」が選択された時のみ、その投影像に再構成演算２０３を行って第２の基底物質変換用マルチエネルギー画像２０３を作成し、これを用いて基底物質の変更の有無Ｓ１０７以降の処理を行っても良い。

ただし、指定マルチエネルギー画像の作成及び表示Ｓ１０５は、図７に示した手順と同様に、基底物質の変更の有無Ｓ１０７の前に完了させておく。このようにすることで、基底物質の変更の有無Ｓ１０７に移るまでの時間を短縮でき、特に基底物質の変更の有無Ｓ１０７で「無」を選択する場合の処理時間を短縮できる。

同様に、例えば基底物質画像２０１を保存しておき、投影像作成Ｓ１０３及び再構成像作成Ｓ１０４は行わない状態で基底物質の変更の有無Ｓ１０７の入力を待つようにし、基底物質の変更の有無Ｓ１０７で「有」が選択された時のみ、その基底物質画像２０１に投影像作成Ｓ１０３及び再構成演算２０３を行って第２の基底物質変換用マルチエネルギー画像２０３を作成し、これを用いて基底物質の変更の有無Ｓ１０７以降の処理を行っても良い。

ただし、指定マルチエネルギー画像の作成及び表示Ｓ１０５は、図７に示した手順と同様に、基底物質の変更の有無Ｓ１０７の前に完了させておく。このようにすることで、同様に、基底物質の変更の有無Ｓ１０７に移るまでの時間を短縮できる。

＜変形例１＞
図７において説明した演算部１０５による演算処理では、第２の基底物質変換部１０５８において、再構成像から第２の基底物質変換処理を行う例について説明したが、第２の基底物質変換処理は、例えば、第１の基底物質群に含まれる基底物質の投影像である基底物質画像２０１や、基底物質画像から作成したマルチエネルギー画像の投影像２０２を用いて、第２の基底物質変換を行うこともできる。

この場合、第２の基底物質変換を行った後に、第２の基底物質群に含まれる各基底物質の物理量の投影像が得られ、基底物質の投影像からマルチエネルギー画像の投影像を生成してもよく、さらにマルチエネルギー画像の投影像に対して再構成処理を行ってマルチエネルギー画像を取得することもできる。
この場合の演算部１０５による基底物質変換に係る演算処理の詳細及び生成されるデータ又は画像の例を、図１１に従って説明する。図１１中、実線は処理の内容を示し、破線は実線で示す処理を実行することによって生成されるデータ又は画像を示す。

図１１に示すように、ＳステップＳ２００において、ユーザが、入力部１１０によりＸ線ＣＴ装置に対して、撮影準備として種々の設定の入力と共に撮影開始指示の入力を行い、撮影を開始する。次のステップＳ２０１において変換パラメータを決定し、ステップＳ２０２において第１の基底物質変換を行って基底物質画像（投影像）２０１を生成する。続いて、ステップＳ２０３において、投影像生成部１０５５が、基底物質毎の投影像に基づいて、マルチエネルギー画像の投影像２０２を生成する。このとき投影像として、指定マルチエネルギー画像用と、第２の基底物質変換用マルチエネルギー画像用とを生成する。

ステップＳ２０４では、ステップＳ２０３で生成された指定マルチエネルギー画像用の投影像に対して再構成像生成部１５３により再構成処理を行ってマルチエネルギー画像を生成し、次のステップＳ２０５でこれを表示部１０６に表示させる。一方、投影像生成部１５１は、生成した第２の基底物質変換用マルチエネルギー画像の投影像を記憶部１０９に保存する（ステップＳ２０６）。

次にステップＳ２０７において、第２の基底物質変換部１５４が、入力部１１０を介して、ユーザによる基底物質の変更の有無の入力を受け付ける。第２の基底物質変換部１５４は、ユーザによる基底物質の変更を行わない旨の入力を受け付けた場合には、処理を終了する。一方、第２の基底物質変換部１５４が、基底物質の変更を行う旨の入力を受け付けた場合には、次のステップＳ２０８に進む。

ステップＳ２０８では、第２の基底物質変換部１５４が、算出する基底物質や基底物質画像の生成条件の入力を、入力部１１０を介して受け付ける。入力された情報に基づいて、次のステップＳ２０９において、第２の基底物質変換部１５４により第２の基底物質変換を行って基底物質画像（投影像）２０６を生成する。ステップＳ２０９における第２の基底物質変換は、記憶部１０９に保存されている第２の基底物質変換用マルチエネルギー画像用の投影像を用いて行う。これは、例えば上述した変換マップ等を適用して行うことができる。

次にステップＳ２１０において、再構成像生成部１５３により、基底物質画像（投影像）２０６を用いて、マルチエネルギー画像の投影像２０７を生成し、次のステップＳ２１１において再構成処理を行って、マルチエネルギー画像２０８を生成し、ステップＳ２１２において主制御部１５９が表示部１０６にマルチエネルギー画像２０８を表示させる。その後、ステップＳ１０７に戻り、再度の基底物質の変更の有無について入力を受け付ける。

第２の基底物質変換部１５４は、ユーザによる基底物質の変更を行わない旨の入力を受け付けた場合に、処理を終了する。一方、第２の基底物質変換部１５４は、基底物質の変更を行う旨の入力を受け付けた場合に、ステップＳ２０８以降の処理を繰り返し、再度マルチエネルギー画像２０８を生成して表示する。このような処理により、第２の基底物質を変更して、それぞれのマルチエネルギー画像を繰り返し取得することができる。

＜変形例２＞
上述した実施形態では、再構成像から第２の基底物質変換処理を行う例について、変形例１では、マルチエネルギー画像の投影像２０２から第２の基底物質変換を行う例について説明した。以下、図１２に従って、第１の基底物質群に含まれる基底物質の投影像である基底物質画像２０１から第２の基底物質変換を行う例について説明する。

図１２示すように、ＳステップＳ３００において、ユーザが、入力部１１０によりＸ線ＣＴ装置に対して、撮影準備として種々の設定の入力と共に撮影開始指示の入力を行い、撮影を開始する。次のステップＳ３０１において変換パラメータを決定し、ステップＳ３０２において第１の基底物質変換を行って基底物質画像（投影像）２０１を生成する。

続いて、ステップＳ３０３において、投影像生成部１５２が、基底物質毎の投影像に基づいて、マルチエネルギー画像の投影像２０２を生成する。このとき投影像として、指定マルチエネルギー画像用の投影像を生成する。

ステップＳ３０４では、ステップＳ２０３で生成された指定マルチエネルギー画像用の投影像に対して再構成像生成部１５３により再構成処理を行ってマルチエネルギー画像を生成し、次のステップＳ３０５でこれを表示部１０６に表示させる。

ステップＳ３０７では基底物質の変更の有無を入力し、基底物質の変更を行わない場合は終了となる。一方、基底物質の変更を行う場合は、次にステップＳ３０８にて基底物質や作成条件入力を行う。この情報を基に、ステップＳ３０９において、第２の基底物質変換用マルチエネルギー画像の投影像を生成する。

この第２の基底物質変換用マルチエネルギー画像の投影像２０９を用いて、ステップＳ３１０からステップＳ３１３の処理を行うことでマルチエネルギー画像２０８を生成し表示させることができる。その後、ステップＳ１０７に戻って、再度の基底物質の変更の有無を入力する。

このように、第２の基底物質変換用マルチエネルギー画像やその投影像を作成する前に、第１の基底物質変換後のマルチエネルギー画像を表示させることができるため、参照したい画像をより早期に表示させることができる。

上述した演算部１０５による演算処理において、例えば、Ｘ線検出素子の感度特性の補正や、欠陥素子の出力値を推定する補正処理などを適宜追加することができる。これらの処理により、アーチファクトの低減やＣＴ値の定量性向上などが実現でき、より高画質な画像を得ることが可能となる。
また、例えば、撮影できるＸ線照射条件が１通りの場合などは、変換パラメータが一義的に定まるので、変換パラメータを決定する処理が不要となる。

また、上記実施形態及びその変形例では、撮影条件として、Ｘ線管電圧を記憶して用いる例について説明したが、例えば、Ｘ線管電流、Ｘ線照射範囲（チャネル方向やスライス方向など）、焦点サイズ、Ｘ線フィルタやボータイフィルタなどの種類や有無、撮影時間、等を１つまたは複数用いることがきる。

上述の実施形態及び変形例では、演算部１０５による処理の流れとして、基底物質のマルチエネルギー画像を表示後、第２の基底物質の情報を入力する例について説明したが、例えば、撮影前に、第２の基底物質の情報や、第２の基底物質のマルチエネルギー画像の生成の有無を撮影前に入力することもできる。また、第１の基底物質から作成したマルチエネルギー画像を表示後に、自動的に第２の基底物質のマルチエネルギー画像を生成して表示することもできる。

以上の説明において、第１の基底物質の物理量については面密度、第２の基底物質の物理量については密度としたが、他のエネルギーに依存しない物理量とすることもできる。例えば、第１の基底物質の物理量を基底物質の長さとすることもできる他、第２の基底物質変換を再構成像にて行う場合は、密度以外にも、基底物質の存在率とすることもできる。さらに、第２の基底物質変換を投影像にて行う場合は、面密度や基底物質の長さとすることもできる。

エネルギー範囲についても、上述の例に限らず、２つのエネルギー範囲とすることもできる他、４つ以上のエネルギー範囲とすることもできる。

上記実施形態及び変形例において、Ｘ線を２つ以上のエネルギー範囲に分別して計数を行うフォトンカウンティング方式の検出器を搭載したＸ線ＣＴ装置について説明したが、Ｘ線ＣＴ装置が、複数のＸ線フォトンのエネルギーを積分して出力する積分型のＸ線検出器であって、２つ以上の異なるＸ線エネルギー分布で被写体を撮影して得た複数の投影データから再構成像を得るデュアルエネルギーＸ線ＣＴ装置や、マルチエネルギーＸ線ＣＴ装置であってもよい。

このとき２つ以上の異なるＸ線エネルギー分布で投影データを得る方法として、Ｘ線検出器内で２つ以上の異なるエネルギー分布にＸ線を分離してそれぞれで検出する方式や、Ｘ線源から２つ以上のスペクトルのＸ線を照射して、それぞれで投影データを得る方式であってもよい。Ｘ線の照射方法として、半回転、１回転、複数回転毎に照射Ｘ線のスペクトルを切り替えても良いし、１ビューや複数ビュー毎に切り替えても良い。またスライス方向に異なるエネルギー分布を照射し、ヘリカル撮影を行うことで、同じパスで異なるエネルギー分布の投影データを得ても良い。

なお、上述した実施形態及びその変形例では、Ｘ線ＣＴ装置により撮像を行って投影データを取得し、この投影データに対して演算処理を行う例について説明したが、投影データを他の撮像装置から取得して演算処理を行う構成とすることもできる。すなわち、別個のＸ線撮像装置により撮像された投影データのからデータを取得して、これに対して演算処理を行う画像処理装置とすることができる。

具体的には、図１３に示すように、演算部１０５、表示部１０６、制御部１０７、記憶部１０９及び入力部１１０を備えた画像処理装置とすることができる。図１３の画像処理装置では、別個のＸ線撮像装置で撮影して取得した投影データを、例えば、入力部１１０により入力して記憶部１０９に記憶させ、演算部１０５がこの投影データに基づいて、投影データ、投影像、マルチエネルギー画像等所望の画像を生成する。

上述した実施形態及び各変形例においてはＸ線ＣＴ装置について説明したが、例えば、画像再構成処理を行わない装置や、Ｘ線源を有さない装置などにも適用することができる。具体的には、Ｘ線画像診断装置、Ｘ線画像撮影装置、Ｘ線透視装置、マンモグラフィー、デジタルサブトラクション装置、Ｘ線検出器、放射線検出器等にも適用することができる。

１０４・・・検出部、１０５・・・演算部、１０６・・・表示部、１０７・・・制御部、１０８・・・信号収集部、１０９・・・記憶部、１１０・・・入力部、１１１・・・Ｘ線検出器、１１２・・・信号処理部、１４０・・・撮影条件情報、１４１・・・変換パラメータ、１５０・・・変換パラメータ決定部、１５１・・・第１の基底物質変換部、１５２・・・投影像生成部、１５３・・・再構像生成部、１５４・・・第２の基底物質変換部、１５９・・・主制御部

Claims (8)

1. エネルギー分布の異なる２種類以上の投影データに基づいて、第１の基底物質群に含まれる２以上の基底物質の物理量をそれぞれ算出する第１の基底物質変換部と、
   前記第１の基底物質群に含まれる２以上の基底物質の物理量から、被写体の投影像又は再構成像の少なくとも一方であって、かつ、複数の画像を生成する画像生成部と、
   複数の前記画像に基づいて、前記第１の基底物質群とは異なる第２の基底物質群に含まれる基底物質の物理量を算出する第２の基底物質変換部と、を備える画像処理装置。
2. 前記画像生成部が、前記第２の基底物質変換部により算出された第２の基底物質群に含まれる基底物質の物理量に基づいて再構成像を生成する請求項１記載の画像処理装置。
3. 前記投影データを撮影した際の撮影条件と、前記投影データと前記基底物質の物理量との関係を示す変換パラメータとを記憶する記憶部を備え、
   前記第１の基底物質変換部が、前記記憶部に記憶された撮影条件に基づいて基底物質の物理量を算出する際に用いる変換パラメータを決定し、該変換パラメータを用いて基底物質の物理量を算出する請求項１記載の画像処理装置。
4. 前記撮影条件が、前記投影データを撮影した際のＸ線管電圧、Ｘ線管電流、Ｘ線照射範囲、焦点サイズ、Ｘ線フィルタの種類、ボータイフィルタの種類、及び撮影時間の少なくとも１つに係る情報を含む請求項３記載の画像処理装置。
5. 前記画像生成部が、生成された画像を表示部に表示させる請求項１記載の画像処理装置。
6. 前記画像生成部が、前記２以上の基底物質の物理量から、所定のエネルギーの単色Ｘ線を照射した場合に得られる仮想単色Ｘ線画像の投影像を、異なる２以上のエネルギーの前記単色Ｘ線についてそれぞれ演算により求め、前記２以上のエネルギーのそれぞれの仮想単色Ｘ線画像の投影像からそれぞれ画像を生成し、
   前記第２の基底物質変換部が、前記仮想単色Ｘ線画像の投影像からそれぞれ生成した前記画像に基づいて、前記前記第２の基底物質群に含まれる前記第２の基底物質の物理量を算出する請求項１記載の画像処理装置。
7. 請求項１乃至請求項６の何れか１項記載の画像処理装置を備えたＸ線ＣＴ装置。
8. エネルギー分布の異なる２種類以上の投影データに基づいて、第１の基底物質群に含まれる２以上の基底物質の物理量をそれぞれ算出する第１の基底物質変換工程と、
   前記第１の基底物質群に含まれる２以上の基底物質の物理量から、被写体の投影像又は再構成像の少なくとも一方であって、かつ、複数の画像を生成する画像生成工程と、
   複数の前記画像に基づいて、前記第１の基底物質群とは異なる第２の基底物質群に含まれる基底物質の物理量を算出する第２の基底物質変換工程と、を備える画像処理方法。

Images