(57)【要約】【課題】  冗長性を生じさせず、演算時間を増大させることなく基底物質の物理量を算出し、所望の画像を取得する。【解決手段】  複数のＸ線検出素子に入射したＸ線フォトンを３以上のエネルギー範囲に分別して計数することにより得られた該エネルギー範囲毎の投影データを取得するデータ取得部と、前記エネルギー範囲毎の投影データに基づいて、前記エネルギー範囲の数よりも少ない数の複数の基底物質の物理量を算出する基底物質演算部と、前記基底物質の物理量を用いて所定の画像を生成する画像生成部と、を備え、前記基底物質演算部が、前記基底物質の物理量の算出に用いる適用エネルギー範囲として、３以上の前記エネルギー範囲のうち、前記基底物質と同数の前記エネルギー範囲を選択し、選択された前記適用エネルギー範囲の投影データに基づいて前記基底物質の物理量を算出する画像処理装置を提供する。【選択図】図１
2018023671.tif 000002
【特許請求の範囲】
【請求項１】
複数のＸ線検出素子に入射したＸ線フォトンを３以上のエネルギー範囲に分別して計数することにより得られた該エネルギー範囲毎の投影データを取得するデータ取得部と、
前記エネルギー範囲毎の投影データに基づいて、前記エネルギー範囲の数よりも少ない数の複数の基底物質の物理量を算出する基底物質演算部と、
前記基底物質の物理量を用いて所定の画像を生成する画像生成部と、を備え、
前記基底物質演算部が、前記基底物質の物理量の算出に用いる適用エネルギー範囲として、３以上の前記エネルギー範囲のうち、前記基底物質と同数の前記エネルギー範囲を選択し、選択された前記適用エネルギー範囲の投影データに基づいて前記基底物質の物理量を算出する画像処理装置。
【請求項２】
前記基底物質演算部において、演算する前記基底物質の数と、選択する前記適用エネルギー範囲の数とが共に２であることを特徴とする請求項１記載の画像処理装置。
【請求項３】
前記基底物質演算部が、前記エネルギー範囲毎の投影データにおけるＸ線フォトンのカウント値に基づいて前記適用エネルギー範囲を選択する請求項１記載の画像処理装置。
【請求項４】
前記基底物質演算部が、前記エネルギー範囲毎の投影データのうち、Ｘ線フォトンのカント値が最大となる前記エネルギー範囲を前記適用エネルギー範囲の一つとして選択する請求項３記載の画像処理装置。
【請求項５】
前記基底物質演算部が、前記エネルギー範囲毎の投影データのうち、Ｘ線フォトンのカント値が大きい順に前記適用エネルギー範囲を選択する請求項３記載の画像処理装置。
【請求項６】
前記基底物質演算部が、前記エネルギー範囲毎の再構成像のうち、雑音レベルが小さい順に前記適用エネルギー範囲を選択する請求項３記載の画像処理装置。
【請求項７】
前記データ取得部が、被検体の情報を含む撮像条件を取得し、
前記基底物質演算部が、前記撮像条件に基づいて前記適用エネルギー範囲を選択する請求項１記載の画像処理装置。
【請求項８】
ユーザにより指定された適用エネルギー範囲の入力を受け付ける入力部をさらに備え、
前記基底物質演算部が、ユーザにより指定され前記入力部に入力されたエネルギー範囲を、前記適用エネルギー範囲として選択する請求項１記載の画像処理装置。
【請求項９】
ユーザによる入力を受け付ける入力部に備え、
前記基底物質演算部が、前記適用エネルギー範囲を、前記エネルギー範囲毎の投影データにおけるＸ線フォトンのカウント値、前記エネルギー範囲毎の再構成における雑音レベル、撮像条件、又はユーザによる前記入力部への入力のうち少なくとも２つに基づいて前記適用エネルギー範囲を選択する請求項１記載の画像処理装置。
【請求項１０】
前記投影データから再構成像を作成する再構成処理部をさらに備え、
前記基底物質演算部は、前記投影データから作成された前記再構成像に基づいて前記基底物質の物理量を算出する請求項１記載の画像処理装置。
【請求項１１】
複数のＸ線検出素子に入射したＸ線フォトンを３以上のエネルギー範囲に分別して計数することにより該エネルギー範囲毎の投影データを生成するＸ線検出部と、
請求項１乃至請求項１０の何れか１項に記載の画像処理装置と、を備えたＸ線撮像装置。
【請求項１２】
複数のＸ線検出素子に入射したＸ線フォトンを３以上のエネルギー範囲に分別して計数することにより得られた該エネルギー範囲毎の投影データを取得するデータ取得工程と、
前記エネルギー範囲毎の投影データに基づいて、前記エネルギー範囲の数よりも少ない数の複数の基底物質の物理量を算出する基底物質演算工程と、
前記基底物質の物理量を用いて所定の画像を生成する画像生成工程と、を備え、
前記基底物質演算工程において、前記基底物質の物理量の算出に用いる適用エネルギー範囲として、３以上の前記エネルギー範囲のうち、前記基底物質と同数の前記エネルギー範囲を選択し、選択された前記適用エネルギー範囲の投影データに基づいて前記基底物質の物理量を算出する画像処理方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本願発明は、画像処理装置、Ｘ線撮像装置及び画像処理方法に関し、特に、フォトンカウンティング型のＸ線検出器により取得されたデータに対する画像処理を行う画像処理装置、Ｘ線撮像装置及び画像処理方法に関する。
【背景技術】
【０００２】
Ｘ線ＣＴ装置は、複数の方向から撮影した被検体のＸ線透過像である投影データからＸ線吸収係数（線減弱計数）を算出し、被検体の断層像である再構成像を得る装置である。従来、投影データを取得するためのＸ線検出器として、積分型のＸ線検出器を用いたＸ線ＣＴ装置が主流であったが、近年、Ｘ線フォトンの個数を計測するフォトンカウンティング型のＸ線検出器を搭載したＸ線ＣＴ装置が開発されつつある（例えば、特許文献１）。この装置では、これまでのＸ線ＣＴ装置では取得できなかったエネルギー毎の疑似単色の再構成像や、原子番号などの分布を示す吸収係数以外の再構成像（以降、これらの画像をマルチエネルギー画像と記す）を生成することができる。このような画像の生成は、例えば、投影データから基底物質の長さや面密度などの物理量を求めることで行われる。
【０００３】
一方、フォトンカウンティング型のＸ線検出器は、複数のエネルギー範囲に分別して、エネルギー範囲毎の投影データを得ることができ、基底物質数よりも多い数のエネルギー範囲に分別されることもある。このとき、すべてのエネルギー範囲を基底物質の物理量の決定に用いると冗長性が生じる。例えば、３つのエネルギー範囲の投影データから２つの基底物質の物理量の決定する場合である。すなわち、２つの基底物質の物理量の決定は２つの投影データで行うことができるが、３つのエネルギー範囲の投影データがある場合、２つの基底物質の物理量は、投影データの選び方で複数の解が存在することになり、冗長となる。これら複数の解が一致する場合も有り得るが、入射するＸ線のフォトン数に揺らぎが存在するため、多くの場合は一致しない。このような場合、例えば、揺らぎをモデル化し共役勾配法などの最適化処理を行って最適な解を求める必要があるが、このような最適化処理は多くの計算量が必要となり、結果の表示に遅延が生じる。特にスキャノ画像や、撮影範囲を確認などのためにリアルタイムでの表示が望ましい投影像、再構成像などの表示の妨げとなる場合もある。
【０００４】
そこで、例えば、特許文献１に記載されたＸ線ＣＴ装置では、エネルギー範囲を加算して、加算後のエネルギー範囲数を基底物質数と同数にし、加算後のエネルギー範囲の投影データを用いて、基底物質の物理量を算出している。この場合、冗長性が生じないため、結果を即時に得ることができる。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２０１４−２３３６３３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、特許文献１のＸ線ＣＴ装置では、複数のエネルギー範囲の投影データを加算するため、加算処理のための演算時間や、加算後のデータの保存のための記憶媒体が必要となってしまう。
また、加算処理後に得た基底物質の物理量に比して、加算処理を行わずに最適化処理を行って求めた基底物質の物理量はその情報量が多いことから、これらから求めた画像（以降、全エネルギー範囲利用画像と記す）が診断に供するために必要になると考えられる。このため、加算して得た画像とは別に、全エネルギー範囲利用画像を生成することが望ましいと考えられる。ところが、全エネルギー範囲利用画像の算出には、加算した画像の生成の際に取得した基底物質の物理量や生成したデータを活用することはできないため、全エネルギー範囲利用画像の生成には、画像生成に必要な演算を改めて行った上で再構成像を生成する必要が生じてしまう。
【０００７】
本発明は上記実情に鑑みてなされたものであり、冗長性を生じさせず、演算時間を増大させることなく基底物質の物理量を算出し、所望の画像を取得することを目的とする。
【課題を解決するための手段】
【０００８】
上記課題を解決するために、本発明は以下の手段を提供する。
本発明の一態様は、複数のＸ線検出素子に入射したＸ線フォトンを３以上のエネルギー範囲に分別して計数することにより得られた該エネルギー範囲毎の投影データを取得するデータ取得部と、エネルギー範囲毎の投影データに基づいて、エネルギー範囲の数よりも少ない数の複数の基底物質の物理量を算出する基底物質演算部と、基底物質の物理量を用いて所定の画像を生成する画像生成部と、を備え、基底物質演算部が、基底物質の物理量の算出に用いる適用エネルギー範囲として、３以上の前記エネルギー範囲のうち、基底物質と同数の前記エネルギー範囲を選択し、選択された前記適用エネルギー範囲の投影データに基づいて基底物質の物理量を算出する画像処理装置を提供する。
このように基底物質の物理量を算出する際に用いる適用エネルギー範囲を、算出する基底物質と同数とするので、基底物質の算出の際に冗長性が生じず所望の画像を少ない演算時間で取得することができる。
【発明の効果】
【０００９】
本発明によれば、冗長性を生じさせず、演算時間を増大させることなく基底物質の物理量を算出し、所望の画像を取得することができる。
【図面の簡単な説明】
【００１０】
【図１】本発明の実施形態に係る画像処理装置を適用したＸ線ＣＴ装置の概略を示すブロック図である。
【図２】本発明の実施形態に係る画像処理装置を適用したＸ線ＣＴ装置における演算部の概略を示すブロック図である。
【図３】図１のＸ線ＣＴ装置の検出部の概略を示す参考図である。
【図４】図１のＸ線ＣＴ装置の検出部において生じた電荷に応じた電圧信号を示すグラフである。
【図５】本発明の実施形態に係る画像処理装置を適用したＸ線ＣＴ装置における演算処理に係るフローチャートである。
【図６】本発明の実施形態に係る画像処理装置を適用したＸ線ＣＴ装置における高速算出モードの基底物質演算処理に係るフローチャートである。
【図７】本発明の実施形態に係る画像処理装置を適用したＸ線ＣＴ装置の高速算出モードの基底物質演算処理において用いる変換マップである。
【図８】本発明の実施形態に係る画像処理装置を適用したＸ線ＣＴ装置における通常算出モードの基底物質演算処理に係るフローチャートである。
【図９】本発明の実施形態に係る画像処理装置を適用したＸ線ＣＴ装置の通常算出モードの基底物質演算処理において用いる変換マップである。
【図１０】本発明の実施形態の変形例２に係る画像処理装置を適用したＸ線ＣＴ装置におけるスキャノ画像生成処理に係るフローチャートである。
【図１１】本発明の実施形態の変形例２に係る画像処理装置を適用したＸ線ＣＴ装置においてスキャノ画像用投影データからスキャノ画像を生成する際の説明図である。
【図１２】本発明の実施形態の変形例３に係る画像処理装置を適用したＸ線ＣＴ装置の表示部に表示されるユーザーインターフェースの一例を示す参考図である。
【図１３】本発明の実施形態の変形例３に係る画像処理装置を適用したＸ線ＣＴ装置における適用エネルギー範囲の選択処理に係るフローチャートである。
【図１４】本発明の実施形態の変形例４に係る画像処理装置を適用したＸ線ＣＴ装置における適用エネルギー範囲の選択処理に係るフローチャートである。
【図１５】本発明の実施形態の変形例５に係る画像処理装置を適用したＸ線ＣＴ装置における適用エネルギー範囲の選択処理に係るフローチャートである。
【図１６】本発明に係る画像処理装置のブロック図である。
【図１７】本発明の実施形態の変形例６に係る画像処理装置を適用したＸ線ＣＴ装置における演算処理に係るフローチャートである。
【発明を実施するための形態】
【００１１】
以下、本発明の一実施形態について、図面を参照して説明する。
本発明に係る画像処理装置は、複数のＸ線検出素子に入射したＸ線フォトンを３以上のエネルギー範囲に分別して計数することにより得られた該エネルギー範囲毎の投影データを取得するデータ取得部と、エネルギー範囲毎の投影データに基づいて、エネルギー範囲の数よりも少ない数の複数の基底物質の物理量を算出する基底物質演算部と、基底物質の物理量を用いて所定の画像を生成する画像生成部と、を備え、基底物質演算部が、基底物質の物理量の算出に用いる適用エネルギー範囲として、３以上の前記エネルギー範囲のうち、基底物質と同数の前記エネルギー範囲を選択し、選択された前記適用エネルギー範囲の投影データに基づいて基底物質の物理量を算出する。
このように構成された画像処理装置や、このような画像処理装置を備えたＸ線ＣＴ装置などのＸ線撮像装置によれば、基底物質の物理量を算出する際に用いる適用エネルギー範囲を、算出する基底物質と同数とするので、基底物質の算出の際に冗長性が生じず所望の画像を少ない演算時間で取得することができる。
【００１２】
以下、より具体的に本発明の実施形態について説明する。
＜実施形態＞
以下、本発明の実施形態に係るＸ線ＣＴ装置について図面を参照して説明する。
【００１３】
図１に示すように、Ｘ線ＣＴ装置は、撮影系としての、Ｘ線源１００と、Ｘ線検出器１１１と、これらＸ線源１００及び検出器１１１の検出部１０４（後述）を対向配置し所定の回転軸を中心に回転するガントリー回転部１０１と、ガントリー回転部１０１の開口内に配置された寝台天板１０３と、これら撮影系の動作に伴いＸ線検出器１１１が取得した信号を処理する信号処理部１１２とを備えている。
【００１４】
Ｘ線源１００は、例えば管電圧で加速した電子ビームをタングステンやモリブデンなどのターゲット金属に衝突させ、その衝突位置（焦点）からＸ線を発生させる。
ガントリー回転部１０１は、Ｘ線源１００及び検出部１０４を互いに対向配置し、所定の回転軸を中心に回転する。ガントリー回転部１０１の中央には、被検体１０２が挿入される開口が設けられ、この開口内に、被検体１０２が寝かせられる寝台天板１０３が配置されている。寝台天板１０３とガントリー回転部１０１とは、所定の方向に相対的に移動可能となっている。
【００１５】
Ｘ線検出器１１１は、入射したＸ線フォトンを検出し、３つのエネルギー範囲に分別して計数を行うフォトンカウンティング方式のＸ線検出素子が複数配置された検出部１０４と、Ｘ線検出素子から出力される投影像を収集する信号収集部１０８とを備えている。検出部１０４の詳細は後述する。
【００１６】
信号処理部１１２は、演算部１０５、表示部１０６、制御部１０７、記憶部１０９及び入力部１１０を備えている。
図２に示すように、演算部１０５は、収集した信号に所定の演算処理を行うため、信号収集部１０８で収集した信号に対して補正処理を行う補正処理部１０５２、マルチエネルギー画像等の再構成像を生成する再構成処理部１０５３、及び補正処理部１０５２と再構成処理部１０５３とを制御する主制御部１０５０を備えている。
【００１７】
補正処理部１０５２は、エア補正部１０５４を含んでいる。
再構成処理部１０５３は、基底物質演算部１０５５、マルチエネルギー画像用投影データ生成部１０５６、及び再構成像生成部１０５７を含んで構成されている。
演算部１０５の計算に用いられるパラメータやデータは、記憶部１０９に保存されており、演算部１０５は必要に応じて記憶部１０９からパラメータ等を読み出し、補正処理、演算処理、画像再構成などの計算を行う。このパラメータやデータには、例えば、エア補正部１０５４が用いるＸ線感度分布やＸ線分布を含む感度・Ｘ線分布データ１４１、基底物質演算部１０５５が面密度値から成る投影データ（以降、面密度投影データと記す）を算出する際や、マルチエネルギー画像用投影データ生成部１０５６がマルチエネルギー画像用の投影データを生成する際に用いるＸ線スペクトル分布と質量吸収係数データなどの計算用データ１４０などが含まれる。
【００１８】
ここで、感度・Ｘ線分布データ１４１は、エネルギー範囲毎に生成して予め記憶部１０９に記憶させておく。感度・Ｘ線分布データ１４１は、例えば、被検体１０２を設けずに、Ｘ線管１００からＸ線を照射してエネルギー毎に投影データを取得し、それらに対してＸ線検出素子毎にビュー方向に加算平均を行い、Ｘ線検出器１０４での出力の平均値によって規格化して生成する。
【００１９】
表示部１０６は、演算部１０５により生成された再構成像などを表示する。
制御部１０７は、Ｘ線源１００の発生駆動源の動作を制御するＸ線制御部、Ｘ線検出器１１１の信号読み出し動作を制御する読み出し制御部、ガントリー回転部１０１の回転と寝台天板１０３の移動を制御する撮影制御部、及びこれら各部全体を制御する全体制御部を備えている。記憶部１０９は、演算部１０５における演算処理に用いられるパラメータやデータ等を記憶している。入力部１１０は、Ｘ線ＣＴ装置における撮影条件等の入力を行う。
【００２０】
演算部１０５及び制御部１０７は、その一部又は全部をＣＰＵ（中央処理装置）、メモリ及び主記憶部１０９を含むシステムとして構築することができ、演算部１０５及び制御部１０７を構成する各部の機能は、予め記憶部に格納されたプログラムをＣＰＵがメモリにロードし、実行することにより実現することができる。また機能の一部は、ＡＳＩＣ（Application Specific Integrated Circuit）やＦＰＧＡ（Field Programmable Gate Array）などのハードウェアで構成することも可能である。
【００２１】
以下の説明においては、特に説明しない限り、上述した撮影系、制御部１０７及び信号処理部１１２を構成する要素は、公知のＸ線ＣＴ装置が備える要素と同様の構成を有し、同様の機能を持つ。
【００２２】
続いて、Ｘ線検出器１１１の検出部１０４について説明する。
検出部１０４は、Ｘ線源１００を略中心とした円弧状に複数配置されており、ガントリー回転部１０１の回転に伴い、Ｘ線源１００との位置関係を保ちながら回転する。なお、図１においては、説明の便宜上、検出部１０４が８個の場合を示したが、実際の装置では、例えば４０個程度である。また検出部１０４の前面にはＸ線グリッド（図示せず）が設置されており、Ｘ線源１００から照射されたＸ線のうち、被検体３００などで散乱されたＸ線が、検出部１０４に入射するのを防ぐ。
【００２３】
検出部１０４は、例えば複数の同一のサイズのＸ線検出素子が、チャネル方向とスライス方向に２次元的に配置され、それぞれの方向で等間隔に配置された構造となっている。Ｘ線検出素子は、図１のＸ線ＣＴ装置において、チャネル方向と回転方向１１３を、スライス方向と回転軸方向１１４とを略一致させて配置されている。
【００２４】
一方、検出部１０４の断面は、例えば、図３に示すように、検出層４０１を挟むように正負の電極４０２，４０３が設けられてＸ線検出素子が構成され、その電極４０２，４０３には、読み出し回路４０５が接続された構造を有する。本実施形態では、正の電極４０２は、各Ｘ線検出素子で共通の構造である。
【００２５】
検出層４０１は、例えばＣｄＴｅ（テルル化カドミウム）、ＣｄＺｎＴｅ（カドミジンクテルル）、Ｓｉ（シリコン）などの半導体材料から成る。Ｘ線は、図２中の矢印４０４で示すように、正の電極４０２側から検出層４０１に入射し、Ｘ線フォトンを検出してそのエネルギーに応じた量の電荷を生じる。そして、検出部１０４は、この電荷に応じたアナログの電気信号を信号収集部１０８に出力する。
【００２６】
検出部１０４からのアナログ信号は、信号収集部１０８の読出し回路４０５に入力される。読出し回路４０５は、入射をトリガーとして読み出し、エネルギー分別とデジタル変換を行う。エネルギー分別は、入射した電荷で生じた電気信号を、所定の閾値により複数のエネルギー範囲に分別する処理である。このとき、発生した電気信号の波高や発生量は、入射したＸ線フォトンのエネルギーに依存するため、Ｘ線フォトンに応じたエネルギー範囲に分別できる。このようにエネルギー範囲に分別してそれぞれでカウントすることで、デジタル信号を得る。信号収集部１０８は、このように得たエネルギー範囲毎のデジタル信号を演算部１０５に出力する。
【００２７】
ここで、読出し回路４０５で行われる分別方法の一例を説明する。ここでは、３つのエネルギー範囲、つまり、所定の低エネルギー閾値未満であるエネルギー範囲（以下、「低エネルギー範囲」という）、所定の低エネルギー閾値以上高エネルギー閾値未満のエネルギー範囲（以下、「中エネルギー範囲」という）、所定の高エネルギー閾値以上のエネルギー範囲（以下、「高エネルギー範囲」という）に判別する例について説明する。
【００２８】
図４は、発生した電荷により生じた電圧信号１２７を示すグラフであり、横軸は時間、縦軸は電圧を示す。また、閾値１３１が低エネルギー閾値を、閾値１３２が高エネルギー閾値を表す。図４の例では、サンプリング時間１２３中にＸ線が入射してパルス出力１２０を生じ、サンプリング時間１２４中にＸ線が入射してパルス出力１２１を生じ、サンプリング時間１２５中にＸ線が入射してパルス出力１２２を生じている。なお、図４では、サンプリングはＸ線が入射するタイミングだけでなく、Ｘ線が入射しない場合（サンプリング時間１２６）でも周期的に行われる場合を示したが、Ｘ線フォトンが入射したタイミングでサンプリングが行われる場合も在り得る。
【００２９】
本実施形態においては、読み出し回路４０５は、入射信号によるＸ線フォトンをエネルギー範囲に分けて分別する前に、Ｘ線フォトンが入射したかを判断する。このため、読み出し回路４０５は、サンプリング毎に、その区間における出力電圧１２７の最大値と、閾値１３０と比較する。ここで閾値１３０は、Ｘ線フォトンの入力の有無を判断する閾値であり、出力電圧１２０が閾値１３０以上の場合には、以下に説明するエネルギー範囲の分別を行い、出力電圧１２０は閾値１３０未満の場合は行わない。このような閾値１３０は、電圧１２７は、Ｘ線が入力しないときにも検出部１０４の回路ノイズによって変動しており、これをＸ線による信号と誤検出しないために必要となる。そのため閾値１３０は、ゼロより大きく低エネルギー閾値１３１より小さな値で設定される。
【００３０】
図４におけるサンプリング時間１２６はＸ線が入射していないが、雑音によって変動が生じた出力電圧１２７の一例であり、出力電圧１２０は低エネルギー閾値１２１未満のため、読み出し回路４０５は、信号のエネルギー範囲の分別を行わない。
【００３１】
Ｘ線が入射して閾値１３０以上の電圧１２７が生じた場合、読み出し回路４０５は、その信号をエネルギー範囲に分別し、それぞれのエネルギー範囲で入射Ｘ線フォトン数をカウントする。この分別は、例えば、出力電圧１２７の最大値を、低エネルギー閾値１３１と高エネルギー閾値１３２と比較することで行う。
【００３２】
例えば、サンプリング時間１２３では、出力電圧１２０は低エネルギー閾値１３１未満のため、低エネルギー範囲と分別する。サンプリング時間１２４では、出力電圧１２１は低エネルギー閾値１３１以上であって、高エネルギー閾値１３２未満のため、中エネルギー範囲と分別する。サンプリング時間１２５では、出力電圧１２０は高エネルギー閾値１３２以上のため、高エネルギー範囲と分別する。
【００３３】
エネルギー範囲が分別されて入射したＸ線フォトン数をエネルギー範囲毎にカウントし、この合計をビュー毎に出力する。サンプリング時間はビューの時間と比較して非常に短く、ビュー間に多数回サンプリングを行うことになる。
このように、入射の有無とエネルギー範囲の分別を行い、信号収集部１０８は、ビュー毎に、エネルギー範囲毎のデジタル信号を生成する。
なおサンプリングでの最大値を用いて分別を行う代わりに、例えば、サンプリング中の出力電圧の積分値を用いてもよく、分別手法は上記手法に限定されない。
【００３４】
このように構成されたＸ線ＣＴ装置では、一般に、以下のように撮影動作が行われる。
まず撮影者が、入力部１１０から撮影条件を入力して実撮影の開始を入力すると、制御部１０７はＸ線源１００からのＸ線の照射と、ガントリー回転部１０１を制御し撮影を開始する。例えば、Ｘ線源１００では、１２０ｋＶの管電圧で電子ビームを加速して、寝台天板１０３に載った被検体３００に向けて照射されＸ線が照射される。被検体３００を透過したＸ線は検出部１０４で検出される。検出部１０４は、Ｘ線検出素子毎に入射Ｘ線のエネルギーに応じた電荷を発生する。信号収集部１０８は、上述の通り、この電荷を高エネルギー範囲、中エネルギー範囲、低エネルギー範囲に分別し、エネルギー範囲毎及びビュー毎にデジタル信号を得て演算部１０５に出力する。
【００３５】
制御部１０７は、ガントリー回転部１０１を回転方向に回転させることで、被検体３００に対するＸ線の照射角度を変化させる。このように回転駆動させながら、ビュー毎に焦点位置を変更させて撮影を繰り返し行い、３６０度分のデジタル信号を取得する。撮影は、例えば０．４度ごとに複数ビューにわたって行われる。このような撮影により、３６０度分のデジタル信号が得られ、３６０度分のデジタル信号を投影データとする。なお、Ｘ線源１００から発生されるＸ線は、各ビューに同期したパルスＸ線でも良いし、連続Ｘ線でも良い。
【００３６】
演算部１０５は、取得した投影データに対し、所定の補正処理や演算処理を行い、マルチエネルギー投影データを生成する。すなわち、補正処理部１０２において、補正処理として、例えば、エア補正を、再構成処理部１０５３では、各演算処理として、例えば、密度画像生成、マルチエネルギー画像用投影データ生成、再構成像生成を行う。
【００３７】
以下、演算部１０５による演算処理について、図５のフローチャートに従って説明する。
図５に示すように、演算部１０５が信号収集部１０８から投影データを受け取ると、ステップＳ１０１において、エア補正部１０５４が信号収集部１０８から受け取った投影データに対してエア補正を行う。エア補正は、例えば、本撮影の事前に計測し生成して記憶部１０９に保存しておいた感度・Ｘ線分布データ１４１を用いて、投影データをエネルギー範囲毎に除することで実現する。エア補正は、各エネルギー範囲で取得した投影データ毎に行う。
【００３８】
ステップＳ１０２では、３つのエネルギー範囲の投影データを用いて、基底物質演算部１０５５において基底物質演算処理を行う。本実施形態に係るＸ線ＣＴ装置では、基底物質の物理量として面密度値を算出することとする。従って、本実施形態における基底物質演算処理は、面密度投影データ算出処理であり、ここでは、撮影範囲の確認等を目的として撮影後に直ちに表示が必要なリアルタイム画像を算出する高速算出モードと、時間をかけて診断に活用することのできる高質な画像を生成するための通常算出モードとの２つの算出モードを有することとして説明する。各モードにおける処理の詳細は後述する。
【００３９】
次のステップＳ１０３では、マルチエネルギー画像用投影データ生成部１０５６において、マルチエネルギー画像用の投影データを生成する。本実施形態では、マルチエネルギー画像の一つである疑似単色画像の投影データを生成することとし、以下、疑似単色画像の投影データ生成方法の一例について説明する。
【００４０】
疑似単色画像は、ある特定の単色のエネルギーの再構成像である。その単色のエネルギーをＥ０、投影データをＰＥ０とし、２つの基底物質を基底物質１と基底物質２としてそれらの投影データの値をδ１とδ２、エネルギーＥ０における質量吸収係数をμｍ１（Ｅ０）とμｍ２（Ｅ０）とすると、基底物質演算部１０５５にて算出した２つの基底物質の面密度投影データの値δ１、δ２を用いて、投影データＰＥ０は、例えば式（１）を用いることで算出することができる。
【００４１】
【数１】
2018023671.tif 000003
【００４２】
ここで、エネルギーＥ０におけるフォトン数Ｓ（Ｅ０）や、各基底物質の質量吸収係数μｍ１（Ｅ０）、μｍ２（Ｅ０）は、シミュレーションやデータベースなどを用いて計算用データ（１４０）として事前に用意しておくことができるため、面密度投影データの値δ１、δ２が求まると、疑似単色画像の投影データが生成できることが分かる。
【００４３】
また面密度投影データの値δ１、δ２は、焦点と検出素子とを結ぶ経路中にある物質の密度を積分したものになるため、以下の説明においては、面密度値という。なお、基底物質とは、例えば、脂肪、水、骨等の被検体を構成する特定の物質を示す。
【００４４】
次のステップＳ１０４では、再構成像生成部１０５７にて、再構成処理を行ってマルチエネルギー画像を取得する。再構成処理は、例えば、従来のＸ線ＣＴ装置で行われている一般的な再構成処理を適用すればよい。このようにして生成したマルチエネルギー画像は、表示部１０６に表示される。
【００４５】
本実施形態では、マルチエネルギー画像として疑似単色画像を生成する例について説明したが、基底物質密度画像、実効原子番号画像、電子密度画像、光電効果画像、コンプトン散乱画像、撮影で用いたスペクトル以外のスペクトルにおける吸収係数像などの画像など、他のさまざまなマルチエネルギー画像を生成することもできる。
【００４６】
続いて、ステップＳ１０２の面密度投影データ算出処理の詳細について説明する。
基底物質演算部１０５５は、高速算出モードの場合には、撮影で得た投影データから、基底物質と同数のエネルギー範囲のデータを抽出し、これらを用いて基底物質の物理量を求める。この基底物質の物理量を算出のために抽出したエネルギー範囲を、以下の説明においては、「適用エネルギー範囲」という。本実施形態においては、基底物質の物理量として、面密度値を求めることとしているため、基底物質演算部１０５５では、面密度値の算出を、全Ｘ線検出素子分及び全ビュー分の投影データに対して行うことで、面密度値データの集合である面密度投影データを得ることができる。
【００４７】
基底物質演算部１０５５における、高速算出モードの演算処理について図６のフローチャートに従って説明する。
ステップＳ２０１で、基底物質演算部１０５５は、投影データを用いて適用エネルギー範囲を決定する。ここで基底物質の物理量の決定は、Ｘ線がより多く入射しているエネルギー範囲を用いることで、ＳＮＲが良く精度が向上するため、Ｘ線フォトンが最も多くカウントされたエネルギー範囲を選択することが望ましい。
【００４８】
従って、ステップＳ２０１では、まず全Ｘ線検出素子及び全ビュー分のカウント値の合計値を、エネルギー範囲毎に算出する。つまり、高エネルギー範囲、中エネルギー範囲、低エネルギー範囲の夫々について、カウント値の合計値を算出する。そして、算出されたエネルギー範囲毎の合計値を用いて、当該合計値の大きい方から順に、基底物質と同数のエネルギー範囲を、適用エネルギー範囲として決定する。ここでは、基底物質が２つであり、中エネルギー範囲と高エネルギー範囲の２つのエネルギー範囲を適用エネルギー範囲として選択したとして説明する。
【００４９】
次のステップＳ２０２において、基底物質演算部１０５５は、投影データ中の１つのＸ線検出素子、ビューに対応するデータ点に対して、基底物質変換処理を行う。
図７に、基底物質変換処理に用いる変換用マップの例を示した。図７（Ａ）に示す中エネルギー範囲の変換用マップ及び図７（Ｂ）に示す高エネルギー範囲の変換用マップは、いずれも基底物質１の面密度値と基底物質２の面密度値との組み合わせから投影データが取り得る値をプロットしたものである。つまり、中エネルギー範囲の変換用マップは中エネルギー範囲の投影データが取り得る値をプロットしたもの、高エネルギー範囲の変換用マップは高エネルギー範囲の投影データが取り得る値をプロットしたものである。
【００５０】
これらの変換用マップは、横軸が基底物質１の面密度値であり、縦軸が基底物質２の面密度値であり、基底物質１と基底物質２との組み合わせから予想される投影データの出力値が記録されたものである。このような変換用マップは、例えばシミュレーションなどにより事前に生成し、計算用データ１４０として記憶部１０９に事前に記憶しておく。曲線１５５は、中エネルギー範囲の変換マップ中で投影データの値が同一である等高線であり、中エネルギー範囲で計測された値と一致するものの集合である。曲線１５６は、高エネルギー範囲の変換マップ中での等高線であり、高エネルギー範囲で計測された値と一致するものの集合である。
【００５１】
これらの変換用マップの等高線から理解されるように、１つのエネルギー範囲の１つの出力値に対して、変換用マップ中に、複数の組み合わせの基底物質１と基底物質２の面密度値が存在する。そこで、曲線１５５と曲線１５６を重ね合わせた変換用マップを生成する（図７（Ｃ））。曲線１５５と曲線１５６との交点１５７は、２つのエネルギー範囲の投影データ値を両立する基底物質１と基底物質２の面密度値であり、交点１５７から基底物質１の面密度値は値１５８、基底物質２の面密度値は値１５９と決定することができる。このようにして、２つのエネルギー範囲の結果から２つの基底物質の面密度値を決定する。
【００５２】
次に、次のステップＳ２０３に進み、基底物質演算部１０５５により基底物質変換処理が完了したか否かを判定する。すなわち、基底物質演算部１０５５は、投影データの全てのデータ点で、言い換えると全てのＸ線検出素子分かつ全てのビューにおける基底物質変換処理が完了しているかを判定し、完了していないと判定された場合には、ステップＳ２０２に戻り、未処理のデータ点について基底物質変換処理を行う。
【００５３】
一方、ステップＳ２０３において、投影データの全てのデータ点に対して基底物質変換処理が完了したと判定された場合には、ステップＳ２０４に進み、得られた基底物質の面密度値の集まりである基底物質の面密度投影データを生成し、処理を終了する。この基底物質の面密度投影データは、基底物質毎に生成する。なお、このようにして生成された基底物質の面密度投影データは、図５のステップＳ１０３におけるマルチエネルギー画像用投影データの生成に用いられる。
【００５４】
続いて、基底物質演算部１０５５における、通常算出モードの演算処理について図８のフローチャートに従って説明する。通常撮影モードは、全てのエネルギー範囲の投影データを用いて基底物質の面密度値を求め、基底物質の面密度投影データを生成するものである。
【００５５】
図８のステップＳ３０１において、基底物質演算部１０５５は、まず投影データ中の１つのＸ線検出素子、ビューに対応するデータ点に対して、基底物質変換処理を行う。
図９に、基底物質変換処理に用いる変換用マップの例を示した。図９（Ａ）は、低エネルギー範囲の変換用マップ、図９（Ｂ）は中エネルギー範囲の変換用マップ、図９（Ｃ）は高エネルギー範囲の変換用マップであり、これらの変換用マップは、いずれも基底物質１の面密度値と基底物質２の面密度値との組み合わせにおいて投影データが取り得る値を示す。各変換用マップは、横軸が基底物質１の面密度値、縦軸が基底物質２の面密度値であり、それらの組み合わせで予想される投影データの出力値が記録されたものである。
【００５６】
これらの変換用マップは、例えばシミュレーションなどにより事前に生成され、計算用データ１４０として記憶部１０９に事前に記憶しておく。曲線１６１は、低エネルギー範囲の変換マップ中で投影データの値が同一である等高線であり、低エネルギー範囲で計測された値と一致するものの集合である。曲線１５５は、中エネルギー範囲の変換マップ中での等高線であり、中エネルギー範囲で計測された値と一致するものの集合である。曲線１５６は、高エネルギー範囲の変換マップ中での等高線であり、高エネルギー範囲で計測された値と一致するものの集合である。
【００５７】
高速算出モードと同様に、１つのエネルギー範囲の１つの出力値に対して、複数の組み合わせの基底物質１と基底物質２の面密度値が存在する。図９（Ｄ）は、曲線１６０と曲線１５５と曲線１５６を重ね合わせた変換用マップである。これらの曲線の交点は２つのエネルギー範囲の投影データ値を両立する基底物質１と基底物質２の面密度値であるが、ここで通常算出モードでは、高速算出モードとは異なり、３つの曲線の交点が位置しない。
【００５８】
図９（Ｄ）の変換用マップでは、２つの曲線の交点１５７−１、交点１５７−２、交点１５７−３の３点が存在し、３つの曲線の共通の交点が存在しない。これは、例えば、照射されたＸ線のフォトン数の揺らぎなどの雑音に起因するものである。すなわち、例えば、１ビュー間に入射したＸ線フォトンが、低エネルギー範囲に属し、かつ、フォトン数が少ない場合、曲線１６１は共通の交点が存在する曲線からはずれるため、共通の交点を結ばなくなることからも理解できる。
【００５９】
このように、３つの曲線に共通する交点が存在しないことから、ステップＳ３０２において、最適化処理を行う。具体的には、雑音を考慮し、最も確からしい基底物質１と基底物質２の面密度値の組み合わせを決定する。この処理では、例えば、それぞれのエネルギー範囲の投影データの出力に対して、各曲線１６１、１５５、１５６はＸ線の量子揺らぎなどによる雑音によってどの程度の誤差を有しているかを求め、その誤差から交点１５７−１、交点１５７−２、交点１５７−３の位置がどの程度の誤差の分布を持つかを求め、最も誤差の少ない位置を最も確からしい基底物質１と基底物質２の面密度値として算出する。このように最適な値を求める方法として、例えば、共役勾配法などを適用することができる。
【００６０】
このように１組の基底物質１と基底物質２の面密度値を決定した後、次のステップＳ３０３に進み、基底物質演算部１０５５により基底物質変換処理が完了したか否かを判定する。すなわち、基底物質演算部１０５５は、投影データの全てのデータ点で、言い換えると全てのＸ線検出素子分かつ全てのビューにおける基底物質変換処理が完了しているかを判定し、完了していないと判定された場合には、ステップＳ３０１に戻り、未処理のデータ点について基底物質変換処理及びステップＳ３０２の最適化処理を行う。
【００６１】
一方、ステップＳ３０３において、投影データの全てのデータ点に対して基底物質変換処理が完了したと判定された場合には、ステップＳ３０４に進み、得られた基底物質の面密度値の集まりである基底物質の面密度投影データを生成し、処理を終了する。この基底物質の面密度投影データは、基底物質毎に生成する。
【００６２】
以上述べたように、高速算出モードは、エネルギー範囲数と基底物質数とが同じであるため、冗長性は無く、一意に基底物質の面密度投影データを求めることができる。このため、撮影直後に撮影範囲を確認するために用いる再構成像のように、リアルタイムでの表示が望ましい再構成像を提供できる。
【００６３】
一方、通常撮影モードは、エネルギー範囲数が基底物質数よりも多いため、冗長性がある。このため、一意に基底物質の面密度投影データを求めることができず、最適化処理を要する。このため、処理に時間を必要として、リアルタイム用の画像生成には適さないが、高速算出モードよりも多くのエネルギー範囲の投影データを用いることができる。
【００６４】
従って、エネルギー情報を用いて、より正確に基底物質の面密度値を決定できる。更に多くのデータを用いているので、高速算出モードよりもＳＮＲも良く、より正確に基底物質の面密度値を決定できる。従って、通常撮影モードにて、診断で使用できる高質な画像を提供できる。
以上のようにＸ線ＣＴ装置を構成することにより、高速に基底物質の物理量を算出し、撮影して直ぐに再構成像を提供することが可能となる。
【００６５】
なお、上記した高速算出モード及び通常撮影モードに関し、各モードによる画像の生成は必要に応じて行うことができる。すなわち、高速算出モード及び通常撮影モードによる画像の生成を所望の順番で双方行うこともでき、何れか一方のみを行うこともできる。
【００６６】
例えば、高速算出モード、に次いで、通常撮影モードの順で画像生成を行い、高速算出モードで得た画像を表示した後、通常撮影モードで得た画像を取得する等、モードの使用についてはユーザが適宜決定することができる。このように、まず高速モードによる画像を表示することで、ユーザは、先に表示された再構成像で撮影直後に撮影範囲を確認することができ、その後に表示された通常モードによる精細な再構成像にて、診断を行うことが可能となる。特に、通常撮影モードで生成した再構成像を表示する際、高速算出モードで生成して先に表示している再構成像と同じ撮影位置のものを、高速算出モードの再構成像を置き換えるように重ねて表示することにより、高速算出モードにより取得した画像がどの位置の画像であるかを認識することが容易となる。
【００６７】
また、本実施形態では、Ｘ線検出器は３つのエネルギー範囲にＸ線フォトンを分別する例について説明したが、エネルギー範囲は３つに限られず、４つ以上であってもよい。同様に、２つの基底物質を同定する例について説明したが、２つ以上で、エネルギー範囲数よりも少ない数の基底物質に分別することもできる。またこのときも、通常算出モードでは、全てのエネルギー範囲を用いて基底物質を算出し、高速算出モードでは、基底物質数と同数の一部のエネルギー範囲を用いて基底物質を算出することができる。
【００６８】
更に、通常算出モードでは、基底物質数よりも多い一部のエネルギー範囲の投影データを選定して用いることができる。通常算出モード及び高速算出モードの両モードとも、一部の複数のエネルギー範囲の投影データを加算して用いることもできる。ただし加算後のエネルギー範囲数は、通常算出モードでは基底物質数よりも多く、高速算出モードでは基底物質数と同数とする。
【００６９】
さらに、本実施形態では、基底物質の物理量として、基底物質演算部１０５５が面密度値を求める例について説明したが、基底物質の物理量は面密度値に限られず、基底物質の長さ（厚さ）であっても良い。基底物質の長さは、基底物質の面密度と密度が分かれば求めることができることから分かるように、基底物質の密度を仮定することや実際に決定することで、求めることができる。
【００７０】
（変形例１）
上述の例では、高速算出モードによる画像生成と通常撮影モードによる画像生成を別個独立の処理として説明したが、これに限られず、例えば、高速算出モードによる画像生成の際に算出したデータを通常算出モードの画像を生成する際に活用することができる。
すなわち、基底物質演算部１０５５は、上述した高速算出モードによる面密度投影データ算出処理において、適用エネルギー範囲を決定して、決定したエネルギー範囲において行った基底物質変換処理の結果を記憶部１０９に記憶しておき、これを通常算出モードによる面密度投影データ算出処理に活用する。
【００７１】
つまり、記憶部１０９には適用エネルギー範囲の基底物質変換処理の結果が記憶されているので、基底物質演算部１０５５はこれをこのまま活用すると共に、適用エネルギー範囲の組み合わせ以外のエネルギー範囲の組み合わせの投影データに対して基底物質変換処理を行う。そして、変換した結果と、記憶部１０９に記憶した変換結果を用いて、最適化処理を行う。最適化処理については上述した処理と同様であるのでここでの説明は省略する。
このように高速算出モードによる画像生成の際に算出したデータを通常算出モードの画像を生成する際に活用することで、通常算出モードにおいて、演算処理を削減することができる。
【００７２】
（変形例２）
上述の実施形態においては、高速算出モードにより再構成画像を生成する例について説明したが、スキャノ画像を生成することもできる。
以下、演算部１０５がスキャノ画像を生成する場合について図１０のフローチャートに従って説明する。
【００７３】
まずステップＳ４０１において、ガントリー回転部１０１を回転せずに、寝台天板１０３を移動して撮影して得た投影データに対して、演算部１０５が信号収集部１０８から投影データを受け取ると、ステップＳ４０１において、エア補正部１０５４が信号収集部１０８から受け取った投影データに対してエア補正を行う。
ステップＳ４０２において、エア補正後の投影データに対して、演算部１０５が、上述した高速算出モードによる面密度投影データ算出処理を行って基底物質の面密度投影データを生成する。
【００７４】
そして、次のステップＳ４０３において、演算部１０５が、スキャノ画像生成処理を行う。スキャノ画像生成処理では、ステップＳ４０２において生成された基底物質の面密度投影データから、まずはスキャノ画像用投影データを生成する。スキャノ画像用投影データＰは、例えば、以下の式（２）を適用することにより算出することができる。
【００７５】
【数２】
2018023671.tif 000004
ただし、δ１、δ２は基底物質の面密度投影データの値、Ｓ（Ｅ）はエネルギーＥにおけるフォトン数、μｍ１（Ｅ）とμｍ２（Ｅ）はエネルギーＥにおける基底物質１と基底物質２の質量吸収係数を表す。
【００７６】
上記式（２）によって得られたスキャノ画像用投影データに基づいて、演算部１０５がスキャノ画像を生成する。具体的には、演算部１０５は、スキャノ画像用投影データとスキャノ画像との画素の中心位置の違いを考慮して、スキャノ画像用投影データの出力値を重み付け加算して、スキャノ画像の画素値を算出する。
【００７７】
スキャノ画像の画素値の算出手法の一例について、図１１に従って説明する。
図１１は、検出部１０４がスライス数を４とし、３ビューのスキャノ画像用投影データから、５スライス分のスキャノ画像１７３を生成する例を説明するための投影データに関する概念図である。このような処理は、スキャノ撮影で得られたスライス方向に位置がずれたデータを合成して、検出部１０４のスライス幅よりも大きな幅の画像を作成するため、更に重複または近いサンプル点でのデータを、スキャノ画像の１点とすることでＳＮＲを向上するために行う。  図１１では１つのチャネルの場合を示しているが、チャネル方向に並ぶ複数のＸ線検出素子の出力に対して行うことは言うまでもない。
【００７８】
図１１において、スキャノ画像用投影データ１７０〜１７２は、夫々１から３ビューで取得したスキャノ画像用投影データであり、位置１７０−ｉ、１７１−ｉ、１７２−ｉ（ｉ＝１，２，３，４）は、各スキャノ画像用投影データ１７０〜１７２のスライスｉのデータが取得された中心位置である。撮影の際に、寝台天板１０３を移動させながらデータを取得するため、各スキャノ画像用投影データ１７０〜１７２のスライスｉのデータが取得された中心位置１７０−ｉ、１７１−ｉ、１７２−ｉ（ｉ＝１，２，３，４）は、スライス方向１１５に対してずれている。
【００７９】
また、スキャノ画像用投影データ１７０〜１７２に基づいて生成されるスキャノ画像１３における位置１７３−ｉ（ｉ＝１，２，３，４，５）は、スキャノ画像１７３のスライスｊのデータの位置を示す。
【００８０】
スキャノ画像１７３は、上述したように、スキャノ画像用投影データとスキャノ画像との画素の中心位置の違いを考慮して、スキャノ画像用投影データの出力値を重み付け加算して算出する。これは、例えば、位置１７３−１のスキャノ画像１７３の画素値は、同一位置である位置１７０−１で取得したスキャノ画像用投影データと、スライス方向の画素幅の半分だけずれた位置で計測した位置１７１−１で取得したスキャノ画像用投影データを、重みをそれぞれ１，０．５として重み付加算する。
【００８１】
同様に、位置１７３−２のスキャノ画像の画素値は、同一位置である位置１７０−２と位置１７２−１で取得したスキャノ画像用投影データの重みを何れも１とし、スライス方向の画素幅の半分だけずれた位置で計測した位置１７１−１と位置１７１−２で取得したスキャノ画像用投影データの重みを何れも０．５として、重み付加算する。他の画素に対しても同様に重み付け加算して、全画素の値を求める。このように重み付加算することで、ＳＮＲの良いスキャノ画像を得ることができる。
【００８２】
次に、重み付けの違いを規格化する。例えば、位置１７３−１のスキャノ画像の画素値は、位置１７０−１と位置１７１−１のスキャノ画像用投影データから生成される。位置１７０−１と位置１７１−１における画素値の重みの合計は１．５であったが、位置１７３−２のスキャノ画像の画素値は、位置１７０−２、位置１７１−１、位置１７１−２、位置１７２−１のスキャノ画像用投影データから生成され、それらの重みの合計は３であった。このため出力差が生じてしまうため、例えばこの重みの合計でそれぞれの合計値を除することで規格化する。
【００８３】
このようにスキャノ画像用投影データ１７０〜１７２に基づいてスキャノ画像１７３を生成することができる。
【００８４】
（変形例３）
上述した第１の実施形態では、演算部１０５が、信号収集部１０８から受け取った投影データを用いて、Ｘ線フォトンのカウント値が大きいエネルギー範囲を順に選択して適用エネルギー範囲を決定する例について説明した。
【００８５】
この他、例えば、予め適用エネルギー範囲の一部を定めておくこともできる。特に、照射Ｘ線のスペクトル内で、実質的に最も高いエネルギー範囲を適用エネルギーの１つとして定めておくことができる。一般に、多くの物質は、エネルギーが高い方が吸収し難くなるため、高エネルギー範囲を適用エネルギー範囲として選択することで、被検体のサイズに依存したＸ線フォトン数の変化が小さく、高速算出モードであっても、安定したＳＮＲで、再構成像などを生成することができるからである。
【００８６】
また、Ｎ個（Ｎは３以上の整数）のエネルギー範囲を適用エネルギー範囲とする際、１個以上（Ｎ−１）個以下のエネルギー範囲を固定しても良いし、照射Ｘ線のスペクトル内で、実質的に最も高いエネルギー範囲から１個以上（Ｎ−１）個以下のエネルギー範囲を適用エネルギー範囲として選択するように予め定めておくことができる。
【００８７】
さらに、投影データに基づいて適用エネルギー範囲を決定するのではなく、他のエネルギー範囲に入射するＸ線フォトン数に関する情報を用いてエネルギー範囲を決定することもできる。例えば、入力部１１０に入力した被検体の情報に基づいて適用エネルギー範囲を定めることができる。
【００８８】
被検体の情報として、具体的には、例えば、大人と子供、体格、体重、撮影部位等が考えられる。大人と子供との比較から検討すると、大人に比して体が小さいためＸ線フォトン数の低減（減弱）が少なく、大人は子供よりも体が大きいためＸ線フォトン数の低減（減弱）が大きいと判断する。減弱が大きいときは、特に低エネルギー範囲でのＸ線フォトン数が大幅に減り、このように大幅に減弱したエネルギー範囲のデータを用いると基底物質の決定精度が低下すると予想される。
【００８９】
従って、子供を撮影する際には、適用エネルギー範囲として低エネルギーと高エネルギー範囲を用い、大人の場合には適用エネルギー範囲として、中エネルギーと高エネルギー範囲を用いる等、体格、体重、撮影部位などの被検体に応じて減弱の大小を判断し、減弱が大きくなると判断できる場合では、そうでない場合よりも高いエネルギー範囲を適用エネルギー範囲として選択することができる。
【００９０】
このように制御することで、被検体による減弱の大きさに応じて適したエネルギー範囲を適用エネルギー範囲として選択することができ、被検体に応じて高速でありながらＳＮＲの向上した画像を生成することができる。
なお、被検体の情報を入力する入力部としては、キーボードや音声入力機器を用いる場合の他にも、体重計や可視光のカメラなどの計測機器等種々のものを適用することができる。
【００９１】
このように、被検体の情報に応じて適用エネルギー範囲を選択する場合は、予め定めた適用エネルギー範囲を記憶部１０９に保存しておくことができる他、適用エネルギー範囲を決定するための条件を記憶部１０９に記憶させておこともできる。又、例えば、各エネルギー範囲のカウント値の合計値を推定して適用エネルギー範囲を決定してもよい。
更に、ユーザが何れのエネルギー範囲を適用エネルギー範囲とするかについて直接選択して決定しても良い。
【００９２】
これらの適用エネルギー範囲の決定するための条件を、ユーザが決定することもできる。
この場合、表示部１０６に表示されるＵＩ（ユーザインターフェース）を介してユーザが適用エネルギー範囲を決定するための条件の入力を行う。図１２に表示部１０６に表示されるＵＩの一例を示した。
図１２に示すＵＩ９００には、適用エネルギー範囲を決定するための条件を入力するための項目として、「撮影条件」、「カウント数」を示す選択バー９０１，９０２の他、直接エネルギー範囲を指定するための「指定」を示す選択バー９０３を備えている。[指定]欄には、領域９０５中にさらに「高エネルギー範囲」、「中エネルギー範囲」、「低エネルギー範囲」を選択することができる選択バー９０４Ａ〜９０４Ｃが設けられ、各エネルギー範囲を選択することができるようになっている。
【００９３】
各選択バー９０１〜９０３は、入力部からの入力で選択が可能である。図１２に示す例では、選択バー９０１が選択された場合に、被検体の情報を用いて適用エネルギー範囲を決定する。選択バー９０２が選択された場合に、投影データに基づいてカウント値から適用エネルギー範囲として決定する。選択バー９０３が選択された場合に、適用エネルギー範囲をユーザが直接指定することができる。すなわち、例えば、選択バー９０３が選択されると、領域９０５の選択がアクティブになり、この中の選択バー９０４−１，９０４−２，９０４−３の内から任意のエネルギー範囲を適用エネルギー範囲として選択することができる。
【００９４】
選択バー９０１〜９０３は、選択されると選択バー中のドットが、白から黒に変わるようになっており、これにより選択されたことをユーザが目視でわかるようになっている。図１２においては、選択バー９０１が選択されている。このような選択は撮影開始前まで可能であり、撮影開始を入力した際に選択されていた選択バーに従った方法で適用エネルギー範囲が決定され、画像が生成される。
【００９５】
具体的に、適用エネルギー範囲の選択処理について図１３のフローチャートに従って説明する。
ステップＳ５０１において、演算部１０５は入力部１１０からの撮影開始の入力と共にユーザにより何れの選択バー９０１〜９０３が選択されたかに係る情報を受けて、ステップＳ５０２に進む。ステップＳ５０２では、ユーザからの入力を受け付けて適用エネルギー範囲の選択方法を決定する。具体的には、ユーザが選択バー９０１を選択している場合には、演算部１０５は「カウント値」が選択されていると判定してステップＳ５０３に進む。また、ユーザが選択バー９０２を選択している場合には、演算部１０５は「撮影条件」が選択されていると判定してステップＳ５０４に進む。さらに、ユーザが選択バー９０３を選択している場合には、演算部１０５は「指定」が選択されていると判定してステップＳ５０５に進む。
【００９６】
ステップＳ５０３では、基底物質演算部１０５５は、信号収集部１０８により受け取った投影データを用いて、エネルギー範囲毎にカウント数の合計値を求め、次のステップＳ５０６において、例えば、合計値が大きい順に適用エネルギー範囲を決定する。
ステップＳ５０４では、基底物質演算部１０５５は、撮影条件に応じた各エネルギー範囲のカウント数の合計値を推定して、ステップＳ５０６において、適用エネルギー範囲を決定する。
【００９７】
ステップＳ５０５では、基底物質演算部１０５５は、選択バー９０３にて「指定」が選択されたことから、適用エネルギー範囲が所定数指定されか否かを判定する。所定数が選択されている場合にはステップＳ５０６に進み、選択されたエネルギー範囲を適用エネルギー範囲として決定する。ステップＳ５０５の判定において、適用エネルギー範囲として所定数が選択されていない場合にはステップＳ５０７に進み、例えば、ＵＩ上に設定すべき項目が不足している旨の表示を行って撮影処理を停止する。ここで、所定数の適用エネルギー範囲とは、本変形例においては基準物質の数と同一数とする。
【００９８】
（変形例４）
なお、上述した変形例３のステップＳ５０７では、選択された適用エネルギー範囲が所定数に満たない場合には、撮影処理を停止する例について説明したが、例えば、適用エネルギー範囲が不足している場合には、その不足分を、撮影条件やカウント数の合計値等に基づいて決定しても良い。すなわち、「指定」を示す選択バー９０３と、「カウント値」を示す選択バー９０１又は「撮影条件」を示す選択バー９０２とを組み合わせて選択することができることとし、不足分については撮影条件又はカウント数の合計値に従って適用エネルギー範囲を決定することができる。
【００９９】
より具体的に、適用エネルギー範囲の選択処理について図１４のフローチャートに従って説明する。
ステップＳ６０１において、演算部１０５は入力部１１０からの撮影開始の入力と共にユーザにより何れの選択バー９０１〜９０３が選択されたかに係る情報を受けて、ステップＳ６０２に進む。ステップＳ６０２では、演算部１０５は、ユーザにより選択バー９０１が選択されている場合には、「カウント値」が選択されていると判定してステップＳ６０３に進み、選択バー９０２が選択されている場合には、「撮影条件」が選択されていると判定してステップＳ６０４に進み、選択バー９０３が選択されている場合には、「指定」が選択されていると判定してステップＳ６０５に進む。
【０１００】
ステップＳ６０３では、基底物質演算部１０５５は、信号収集部１０８により受け取った投影データを用いて、エネルギー範囲毎にカウント数の合計値を求め、次のステップＳ６０６において、例えば、合計値が大きい順に適用エネルギー範囲を決定する。
ステップＳ６０４では、基底物質演算部１０５５は、撮影条件に応じた各エネルギー範囲のカウント数の合計値を推定して、ステップＳ６０６において、適用エネルギー範囲を決定する。
【０１０１】
ステップＳ６０５では、基底物質演算部１０５５は、選択バー９０３にて「指定」が選択されたことから、適用エネルギー範囲が所定数指定されか否かを判定する。所定数が選択されている場合にはステップＳ６０６に進み、選択されたエネルギー範囲を適用エネルギー範囲として決定する。
ステップＳ６０５の判定において、適用エネルギー範囲として所定数が選択されていない場合にはステップＳ６０７に進み、何れの選択バー９０１〜９０３が選択されたかを判定する。つまり、「カウント値」か「撮影条件」のどちらが選択されているか判定する。ステップＳ６０７において、選択バー９０１が選択されている場合には、「カウント値」が選択されていると判定してステップＳ６０８に進み、選択バー９０２が選択されている場合には、「撮影条件」が選択されていると判定してステップＳ６０９に進み、何れも選択されていない場合にはステップＳ６１０に進む。
【０１０２】
ステップＳ６０８では、基底物質演算部１０５５が、不足分のエネルギー範囲について、エネルギー範囲毎にカウント数の合計値を求め、次のステップＳ６０６において、例えば、合計値が大きいものを優先的に適用エネルギー範囲として決定する。
ステップＳ６０９では、基底物質演算部１０５５は、撮影条件に応じた各エネルギー範囲のカウント数の合計値を推定して、ステップＳ６０６において、適用エネルギー範囲を決定する。
【０１０３】
ステップＳ６１０では、設定項目が不足している旨を表示部に表示してユーザに注意喚起し、撮影を停止する。
このように、指定された適用エネルギー範囲と共に、全ての適用エネルギー範囲を決定することができる。
【０１０４】
なお、ステップＳ６１０において「カウント値」と「撮影条件」の両方とも指定されていない場合は、デフォルトとして自動的に「カウント値」と「撮影条件」のどちらか一方が選択されるように設定することもできる。
上述の説明においては、適用エネルギー範囲の選択に関し、「撮影条件」に基づいて選択する例、「カウント数」に基づいて選択する例、ユーザが「指定」する例の３つの例について説明したが、「撮影条件」、「カウント数」及び「指定」のうち、２つから選択することもでき、「撮影条件」、「カウント数」及び「指定」に他の方法を加えたＮ個（Ｎは４以上の整数）から２つ以上（Ｎ−１）の方法から選択するように設定することもできる。
【０１０５】
（変形例５）
本実施形態の高速算出モードでは、投影データを用いてカウント数の合計値の高い順に適用エネルギー範囲を選択し、再構成像を生成する場合について説明した。この他、例えば、適用エネルギー範囲の組み合わせを複数パターン決定し、夫々の組み合わせに対して夫々再構成像を生成し、複数の再構成像の中から最適な再構成像を選択することもできる。
以下、このような処理について、図１５のフローチャートに従って説明する。
【０１０６】
ステップＳ７０１において、基底物質演算部１０５５は、次のステップＳ７０２における基底物質変換処理に用いる適用エネルギー範囲を決定する。本変形例では、例えば、高エネルギー範囲、中エネルギー範囲及び低エネルギー範囲の３つのエネルギー範囲から、適用エネルギー範囲として２つを組み合わせて選択する場合について説明する。従って、適用エネルギー範囲の組み合わせとしては、高エネルギー範囲と中エネルギー範囲との組み合わせ、高エネルギー範囲と低エネルギー範囲の組み合わせ、及び、中エネルギー範囲と低エネルギー範囲の組み合わせの３パターンが考えられ、このパターンを予め定めた順番に従って、順次処理していく。
【０１０７】
次のステップＳ７０２において、基底物質演算部１０５５は、投影データ中の１つのＸ線検出素子、ビューに対応するデータ点に対して、基底物質変換処理を行う。基底物質変換処理は、上述した第１の実施形態における基底物質変換処理と同様に、例えば、図７に示す変換用マップを用いて行われる。
【０１０８】
続いて、ステップＳ７０３において、基底物質演算部１０５５は、投影データの全てのデータ点に対して基底物質変換処理が完了したか否かを判定する。ステップＳ７０３において全てのデータ点に対して基底物質変換処理が完了していないと判定された場合には、ステップＳ７０２に戻り基底物質変換処理を繰り返す。ステップＳ７０３において、全てのデータ点に対して基底物質変換処理が完了したと判定された場合には、ステップＳ７０４に進む。ステップＳ７０４では、基底物質演算部１０５５が、得られた基底物質の面密度値の集まりである基底物質の面密度投影データを生成し、記憶部１０９に記憶する。この基底物質の面密度投影データは、基底物質毎に生成する。
【０１０９】
ステップＳ７０５では、適用エネルギー範囲の全ての組み合わせに対して基底物質変換処理が終了しているか否かを判定する。ステップＳ７０５の判定において、適用エネルギー範囲の全ての組み合わせに対して基底物質変換処理が終了していない場合は、ステップＳ７０１に戻り、未処理の適用エネルギー範囲の組み合わせに対する処理を行う。ステップＳ７０５の判定において、適用エネルギー範囲の全ての組み合わせに対して基底物質変換処理が終了したと判定した場合には、ステップＳ７０６に進み、マルチエネルギー画像用投影データ生成部１０５６により夫々マルチエネルギー画像用投影データを生成し、ステップＳ７０７において、再構成像生成部１０５７が再構成画像を生成する。
【０１１０】
次のステップＳ７０８では、先のステップＳ７０７で生成された適用エネルギー範囲の全ての組み合わせに対する再構成画像のうち、最適な画像を選択する。最適な画像の選択は、公知の最適化手法を用いて行われ、例えば、再構成像中の特定位置のＣＮＲ（コントラスト−ノイズ比）が最も良いものを選択することができる。ここで、ＣＮＲについての特定位置は、予め定めておいてもよく、また、入力部１１０にて指定された位置であっても良い。また、雑音レベルなどの物理量を基に、最も良い再構成像を選択しても良い。最適な画像としては、予め定めた条件に鑑みて適宜選択するように設定することができる。そして、選択された再構成像の投影データのエネルギー範囲が、最終的な適用エネルギー範囲となる。
ステップＳ７０９において、演算部１０５は、選択された最適な画像を表示部１０６に表示させ処理を終了する。
【０１１１】
なお、上述の説明では、全ての通りのエネルギー範囲から複数の適用エネルギー範囲の組み合わせに対して夫々基底物質の投影データを生成したが、必ずしもすべての組み合わせに対して投影データを生成する必要はなく、一部の適用エネルギー範囲の組み合わせについてのみ投影データを生成することとしてもよい。これにより処理を低減し、処理時間を短縮できる。
【０１１２】
また、上述の例では、エネルギー範囲の全ての組み合わせから基底物質の投影データを生成した後に再構成像を取得して最適な画像を選択しているが、適用エネルギー範囲の何れか１つの組み合わせに対して、基底物質変換処理から再構成像の生成までを順次行ってから、最適画像の選択を行ってもよい。つまり、処理の順序は上記例に限られない。
【０１１３】
表示部１０６への画像の表示は、通常撮影モードで生成した再構成像を、高速撮影モードで生成した再構成像に重ねて表示することができる他、２つの再構成像を並列して、また別個独立に表示することもできる。
【０１１４】
（変形例６）
上述した本実施形態における高速算出モードでは、投影データを用いて基準物質の物理量を算出し、再構成像を生成する場合について説明した。この他に、エネルギー範囲毎の投影データから再構成像を作成し、その中から最適な再構成像の組を選定して、基準物質の物理量を求めてマルチエネルギー画像を作成しても良い。本変形例における基準物質の物理量は、例えば密度である。
【０１１５】
このような処理の一例を、図１７を参照して説明する。
演算部１０５が信号収集部１０８から投影データを受け取ると、例えば、図５で説明した場合と同様に、ステップＳ８０１において、エア補正部１０５４が信号収集部１０８から受け取った投影データに対してエア補正を行う。次のステップＳ８０２では、再構成像生成部１０５７にて、再構成処理を行う。再構成処理は、例えば、従来のＸ線ＣＴ装置で行われている一般的な再構成処理を適用すればよい。
【０１１６】
続いて、ステップＳ８０３の密度画像算出処理の詳細について説明する。
基底物質演算部１０５５は、高速算出モードの場合には、ステップＳ８０２において取得した再構成像から、基底物質と同数の適用エネルギー範囲のデータを抽出し、これらを用いて基底物質の物理量を求める。これらの適用エネルギー範囲には、例えば雑音のレベルが低い順に、基準物質数だけのエネルギー範囲を選択する。このような選択方法は、投影データにおいてＳＮＲが良い範囲を選ぶ方法と同じである。
【０１１７】
本実施形態においては、基底物質の物理量として、密度値を求めることとしているため、基底物質演算部１０５５では、密度値の算出を、再構成像の全てのボクセルに対して行うことで、密度値データの集合である再構成像（以降、密度画像と記す）を得ることができる。このような物理量への変換は、例えば図７で説明した面密度値を求める場合と同様に、変換マップを用いる。ただし本変形例での変換マップは、横軸が基底物質１の密度値であり、縦軸が基底物質２の密度値であり、基底物質１と基底物質２との組み合わせから予想される再構成像のボクセル値が記録されたものとなる。
【０１１８】
次に、ステップＳ８０４のマルチエネルギー画像作成処理では、密度画像を用いてマルチエネルギー画像を求める。この一例として、例えば、それぞれの基準物質の密度画像に、所定のエネルギーの質量吸収係数を乗じて和を取ることで、その所定のエネルギーの疑似単色画像を作成する。なお、このマルチエネルギー画像は一例であり、さまざまなマルチエネルギー画像を作成する場合があることは言うまでもない。また、ステップＳ８０４のマルチエネルギー画像作成処理では処理をせず、密度画像をマルチエネルギー画像とすることもできる。更に、密度画像を順投影処理して面密度の物理量から成る投影データを求め、先に記したように、面密度からさまざまなマルチエネルギー画像を求めることもできる。
そして、ステップＳ８０５において、、以上のように作成したマルチエネルギー画像を表示部１０６に表示させる。
【０１１９】
このような処理により、最適な再構成像の組を選定して、基準物質の物理量を求めてマルチエネルギー画像を作成することができる。
本変形例では、物理量が密度の場合について説明したが、物理量は適宜選択することができ、例えば、物理量を存在率とすることもできる。これは、上述した実施形態及びその変形例において、面密度から成る投影データを再構成すると、再構成像において密度が物理量となること、長さから成る投影データを再構成すると存在率となることから、存在率も物理量に成り得ることは容易に理解できる。
【０１２０】
本変形例では、適用エネルギー範囲として、再構成像の雑音のレベルが低い順に、エネルギー範囲を適用エネルギー範囲として選択する場合について説明したが、これは一例であり、本発明を限定するものではない。上述した実施形態及びその変形例のように、例えば、実質的に最も高いエネルギー範囲を適用エネルギー範囲の１つとして事前に定めておき、他の適用エネルギー範囲を再構成像を用いて選択しても良い。更に事前に定めておく適用エネルギー範囲が複数あって、他の１つ以上の適用エネルギー範囲を選択するようにして、全ての適用エネルギー範囲を決定しても良い。
【０１２１】
以上述べた実施形態及びその変形例では、適用エネルギー範囲を決定するに際し、投影データの全データ点を用いてカウント数の合計値を算出する例について説明したが、例えば、一部のＸ線検出素子のデータ点や、一部のビューのデータ点のカウント数の合計値を算出することとしてもよい。
【０１２２】
なお、上述した実施形態及びその変形例では、Ｘ線ＣＴ装置により撮像を行って投影データを取得し、この投影データに対して演算処理を行う例について説明したが、投影データを他の撮像装置から取得して演算処理を行う構成とすることもできる。すなわち、別個のＸ線撮像装置により撮像された投影データのからデータを取得して、これに対して演算処理を行う画像処理装置とすることができる。
【０１２３】
具体的には、図１６に示すように、演算部１０５、表示部１０６、制御部１０７、記憶部１０９及び入力部１１０を備えた画像処理装置とすることができる。図１６の画像処理装置では、別個のＸ線撮像装置で撮影して取得した投影データを、例えば、入力部１１０により入力して記憶部１０９に記憶させ、演算部１０５がこの投影データに基づいて、投影データ、スキャノ画像、再構成像等所望の画像を生成する。
【０１２４】
上述した実施形態及び各変形例においてはＸ線ＣＴ装置について説明したが、例えば、画像再構成処理を行わない装置や、Ｘ線源を有さない装置などにも適用することができる。具体的には、Ｘ線画像診断装置、Ｘ線画像撮影装置、Ｘ線透視装置、マンモグラフィー、デジタルサブトラクション装置、Ｘ線検出器、放射線検出器等にも適用することができる。
【符号の説明】
【０１２５】
１０４・・・検出部、１０５・・・演算部、１０６・・・表示部、１０７・・・制御部、１０８・・・信号収集部、１０９・・・記憶部、１１０・・・入力部、１１１・・・Ｘ線検出器、１１２・・・信号処理部、１４０・・・計算用データ、１４１・・・感度・Ｘ線分布データ、１０５２・・・補正処理部、１０５３・・・再構成処理部、１０５４・・・エア補正部、１０５５・・・基底物質演算部、１０５６・・・マルチエネルギー画像用投影データ生成部、１０５７・・・再構成像生成部
【図１】
2018023671.tif 000005
【図２】
2018023671.tif 000006
【図３】
2018023671.tif 000007
【図４】
2018023671.tif 000008
【図５】
2018023671.tif 000009
【図６】
2018023671.tif 000010
【図７】
2018023671.tif 000011
【図８】
2018023671.tif 000012
【図９】
2018023671.tif 000013
【図１０】
2018023671.tif 000014
【図１１】
2018023671.tif 000015
【図１２】
2018023671.tif 000016
【図１３】
2018023671.tif 000017
【図１４】
2018023671.tif 000018
【図１５】
2018023671.tif 000019
【図１６】
2018023671.tif 000020
【図１７】
2018023671.tif 000021