特許 6785732 Ｘ線ＣＴ装置及び画像生成方法

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】6785732

(24)【登録日】2020年10月29日

(45)【発行日】2020年11月18日

(54)【発明の名称】Ｘ線ＣＴ装置及び画像生成方法

(51)【国際特許分類】

   A61B 6/03 20060101AFI20201109BHJP

【ＦＩ】

   A61B6/03 373

   A61B6/03 320R

   A61B6/03ZDM

【請求項の数】6

【全頁数】19

(21)【出願番号】特願2017-162278(P2017-162278)

(22)【出願日】2017年8月25日

(65)【公開番号】特開2019-37550(P2019-37550A)

(43)【公開日】2019年3月14日

【審査請求日】2019年12月17日

(73)【特許権者】

【識別番号】000005108

【氏名又は名称】株式会社日立製作所

(74)【代理人】

【識別番号】110000888

【氏名又は名称】特許業務法人 山王坂特許事務所

(72)【発明者】

【氏名】小嶋 進一

(72)【発明者】

【氏名】渡辺 史人

【審査官】 遠藤 直恵

(56)【参考文献】

【文献】 特開２０１４−０１４４４５（ＪＰ，Ａ）

【文献】 国際公開第２０１４／２０３７８３（ＷＯ，Ａ１）

【文献】 特開２００６−１０１９２６（ＪＰ，Ａ）

【文献】 特開２０１１−０８５４７９（ＪＰ，Ａ）

【文献】 特開２０１６−０４２９３６（ＪＰ，Ａ）

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ａ６１Ｂ ６／００−６／１４

(57)【特許請求の範囲】

【請求項1】

寝台上の被写体の周囲を周回するＸ線源から照射され前記被写体を透過したＸ線フォトンを検出する検出器と、
該検出器によって検出されたＸ線フォトンを収集し処理することにより、予め定めた複数のエネルギー範囲について、エネルギー範囲毎に計数値を出力するデータ収集部と、
前記データ収集部により出力された計数値に対して、前記エネルギー範囲毎に定められ、該エネルギー範囲の所定値に比例した重みを乗じて加算処理を行うデータ加算部と、
前記データ加算部により加算されたデータを用いで再構成処理を行うことにより画像を生成する画像生成部と、
前記データ加算部が処理に用いる前記エネルギー範囲の前記所定値を、補正係数によりずらすことにより前記重みを補正する重み補正部と、
サイズの異なる複数種類のファントムに、前記Ｘ線源からＸ線フォトンを照射して前記検出器により検出した計数値、または、シミュレーションにより算出した前記計数値を用いて、前記エネルギー範囲ごとに予め算出した平均エネルギーと、そのエネルギー範囲について予め定められた前記所定値とのずれ量を、前記エネルギー範囲の前記所定値の前記補正係数として、予め記憶する記憶部とを有し、
前記重み補正部は、前記被写体について予め撮像した画像に基づいて、前記被写体のサイズを算出し、前記被写体のサイズに対応する前記ファントムを用いて求めた前記補正係数を前記記憶部から読み込み、読み込んだ前記補正係数を用いて前記エネルギー範囲ごとに前記所定値をずらすことにより前記重みを補正することを特徴とするＸ線ＣＴ装置。

【請求項2】

請求項１に記載のＸ線ＣＴ装置であって、前記所定値は、前記エネルギー範囲の中央値、最大値および最小値のいずれかであることを特徴とするＸ線ＣＴ装置。

【請求項3】

前記記憶部には、入射Ｘ線のエネルギー分布に応じて前記重みを補正するための第２の補正係数が予め記憶され、
前記データ加算部が、加算処理に際して前記補正係数および前記第２の補正係数を用いて重みを補正する請求項１記載のＸ線ＣＴ装置。

【請求項4】

前記記憶部には、前記検出器の特性に応じて前記重みを補正するための第３の補正係数が予め記憶され、
前記データ加算部が、加算処理に際して前記補正係数および前記第３の補正係数を用いて重みを補正する請求項１記載のＸ線ＣＴ装置。

【請求項5】

前記データ収集部が、前記エネルギー範囲を４以上設定し、
各前記エネルギー範囲において、上限値と下限値との差が前記検出器が検出可能なエネルギー範囲の半値幅より大きい請求項１乃至請求項４の何れか１項記載のＸ線ＣＴ装置。

【請求項6】

寝台上の被写体の周囲を周回するＸ線源から照射され前記被写体を透過したＸ線フォトンを検出する工程と、
検出されたＸ線フォトンを収集し処理することにより、予め定めた複数のエネルギー範囲についてエネルギー範囲毎に計数値を出力する工程と、
前記計数値に対して、前記エネルギー範囲毎に定められ、該エネルギー範囲の中央値に比例した重みを乗じて加算処理を行う工程と、
前記加算処理によって得られたデータを用いで再構成処理を行うことにより画像を生成する工程とを有する画像生成方法であって、
前記加算処理を行う工程では、
前記被写体について予め撮像した画像に基づいて、前記被写体のサイズを算出し、
前記被写体のサイズに対応するファントムを用いて求めておいた補正係数を記憶部から読み込み、
読み込んだ前記補正係数を用いて前記エネルギー範囲ごとに前記中央値をずらすことにより前記重みを補正する前記中央値を、補正係数によりずらすことにより前記重みを補正し、
前記記憶部には、サイズの異なる複数種類のファントムに、前記Ｘ線源からＸ線フォトンを照射して前記検出器により検出した計数値、または、シミュレーションにより算出した前記計数値を用いて、前記エネルギー範囲ごとに予め算出した平均エネルギーと、そのエネルギー範囲の中央値とのずれ量を、前記エネルギー範囲の前記中央値の前記補正係数として、予め格納されていることを特徴とする画像生成方法。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、入射したＸ線の光子を計数することにより入射したＸ線のエネルギー情報を取得するフォトンカウンティング（ｐｈｏｔｏｎ ｃｏｕｎｔｉｎｇ）方式の検出器を備えたＸ線ＣＴ装置（Ｐｈｏｔｏｎ Ｃｏｕｎｔｉｎｇ Ｃｏｍｐｕｔｅｄ Ｔｏｍｏｇｒａｐｈｙ・・・以下、ＰＣＣＴ装置と呼ぶ。）及び画像生成方法に関する。

【背景技術】

【0002】

Ｘ線ＣＴ装置は、被写体を挟んで対向配置されたＸ線源とＸ線検出器の対を回転させながら被写体のＸ線透過データ（投影データ）を得、その断層画像（以下、ＣＴ画像とする）を計算により再構成する装置であり、工業用及びセキュリティ用の検査装置や医学用の画像診断装置等として用いられる。

【0003】

このようなＸ線ＣＴ装置では、半導体検出器やシンチレータ検出器などを適用して、検出器に入射したＸ線のエネルギーに比例した電荷量又は光を生じさせ、これを検出して得られた投影データを再構成することにより断層画像を取得している。より具体的な例として、シンチレータ検出器では、シンチレータ検出器に例えば２０ｋｅＶのＸ線が入射した場合、検出器内で２０ｋｅＶに相当する、つまり２０ｋｅＶのエネルギーに比例した光子が発生し、それを光電子増倍管で電気信号に変換して出力する。従って、Ｘ線ＣＴ装置では、入射したＸ線のエネルギーに応じた信号値、延いてはＣＴ値を得ることができる。

【0004】

ところで、近年、フォトンカウンティング方式の検出器を搭載したＸ線ＣＴ装置であるＰＣＣＴ装置が知られている（例えば、特許文献１及び特許文献２）。ＰＣＣＴ装置では、フォトンカウンティング方式の検出器により、Ｘ線源から照射されて被写体を透過したＸ線の光子（Ｘ線フォトン）を検出素子毎に計数（カウント）し、これによって得られるカウント値に基づいてＸ線ＣＴ画像を再構成する。ＰＣＣＴ装置では、例えば、Ｘ線が透過した被写体の内部組織を構成する元素を推定可能なスペクトラムを取得し、元素レベルの違いが詳細に描出されたＸ線ＣＴ画像を得ることができる。また、ＰＣＣＴ装置では、物質毎に存在するｋ−エッジ付近のエネルギー範囲のＸ線のみを抽出して画像化し、ｋ−エッジが存在する元素のみの画像を作成して診断に用いることができる。

【先行技術文献】

【特許文献】

【0005】

【特許文献1】特開２０１５−１９８８３３号公報

【特許文献2】国際公開２０１２／０２９４９６号

【発明の概要】

【発明が解決しようとする課題】

【0006】

上述のようにシンチレータ検出器や半導体検出器等を適用したＸ線ＣＴ装置では、入射したＸ線のエネルギーに応じた信号値、延いてはＣＴ値を得ることができる。
しかしながら、ＰＣＣＴ装置では、予め設定したエネルギー幅におけるカウント数を得ることができるに過ぎず、各カウント値について詳細なエネルギー情報を得ることができない。例えば、エネルギー範囲が４０〜６０ｋｅＶに設定されており、その範囲に所定の時間内に５カウント得られた場合、それらの各々のエネルギーがどの程度であるかの情報は得られない。従って、ＰＣＣＴ装置では、エネルギーに応じた情報が得られないことから、シンチレータ検出器等を適用した一般的なＸ線ＣＴ装置で得られるＣＴ値を得ることができない。
このため、これまでのＸ線ＣＴ装置による検査での知見を生かすために、ＰＣＣＴ装置においても、従来と同等のＣＴ値を得ることが求められている。

【0007】

ここで、仮想単色画像等からＣＴ値を求めることも考えられるが、基底物質に分解した時点で選択した基底物質の組み合わせで該当する物質の減衰曲線が表される必要がある。その結果、複数の基底物質を選んでも、該当物質がそれらの組み合わせで表現できなければ誤差が発生し、その物質のＣＴ値を精度よく求めることができない。よって、ＰＣＣＴ装置においても、正確なＣＴ値を得るために、投影データを求めることが求められている。

【0008】

本発明は上記実情に鑑みてなされたものであり、フォトンカウンティング方式の検出器を適用したＸ線ＣＴ装置において、精度よく投影データ延いてはＣＴ値を取得することを目的とする。

【課題を解決するための手段】

【0009】

上記課題を解決するために、本発明は以下の手段を提供する。
本発明は、寝台上の被写体の周囲を周回するＸ線源から照射され前記被写体を透過したＸ線フォトンを検出する検出器と、該検出器によって検出されたＸ線フォトンを収集し処理することにより、予め定めた複数のエネルギー範囲について、エネルギー範囲毎に計数値を出力するデータ収集部と、前記データ収集部により出力された計数値に対して、前記エネルギー範囲毎に定められ、該エネルギー範囲の中央値に比例した重みを乗じて加算処理を行うデータ加算部と、前記データ加算部により加算されたデータを用いで再構成処理を行うことにより画像を生成する画像生成部と、を備えたＸ線ＣＴ装置を提供する。

【発明の効果】

【0010】

本発明によれば、フォトンカウンティング方式の検出器を適用したＸ線ＣＴ装置において、精度よく投影データ、延いてはＣＴ値を取得することができる。

【図面の簡単な説明】

【0011】

【図1】本発明の第１の実施形態に係るＸ線ＣＴ装置の概略構成を示すブロック図である。

【図2】本発明の第１の実施形態に係るＸ線ＣＴ装置の検出器の一部を示すブロック図である。

【図3】本発明の第１の実施形態に係るＸ線ＣＴ装置の演算部の概略構成を示すブロック図である。

【図4】本発明の第１の実施形態に係るＸ線ＣＴ装置における撮像処理の流れを示すフローチャートである。

【図5】（ａ）は水の減弱係数を示すグラフであり、（ｂ）は１００ｍｍの水を透過した場合において、全てのエネルギー範囲での入力Ｘ線の量を１とした場合に検出される量を示すグラフである。

【図6】本発明の第２の実施形態に係るＸ線ＣＴ装置における撮像処理の流れを示すフローチャートである。

【図7】Ｘ線管のエネルギー分布を示し、特に、設定管電圧に依存するピークと、Ｘ線管のターゲット材質に起因する複数のピークを示すグラフである。

【発明を実施するための形態】

【0012】

本発明の実施形態に係るＸ線ＣＴ装置は、寝台上の被写体の周囲を周回するＸ線源から照射され被写体を透過したＸ線フォトンを検出する検出器と、検出器によって検出されたＸ線フォトンを収集し処理することにより、予め定めた複数のエネルギー範囲について、エネルギー範囲毎に計数値を出力するデータ収集部と、データ収集部により出力された計数値に対して、エネルギー範囲毎に定められ、エネルギー範囲の中央値に比例した重みを乗じて加算処理を行うデータ加算部と、データ加算部により加算されたデータを用いで再構成処理を行うことにより画像を生成する画像生成部と、を備えている。

【0013】

以下、より詳細に本発明の実施形態について、図面を参照してより詳細に説明する。
（第１の実施形態）
本実施形態に係るＸ線ＣＴ装置は、フォトンカウンティング方式の検出器を備えたフォトンカウンティングＣＴ装置（ＰＣＣＴ装置）であり、被写体を透過したＸ線に由来する光子（Ｘ線フォトン）を検出器において計数する。個々のＸ線フォトンは、異なるエネルギーを有しており、検出器では、Ｘ線フォトンを予め定めたエネルギー帯毎に弁別して計数する。これにより、当該エネルギー帯毎のＸ線フォトンの数、すなわちＸ線強度を得る。

【0014】

図１に示すように、本実施形態に係るＰＣＣＴ装置１００は、ＵＩ部２００と、計測部３００と、演算部４００と、を備えている。
ＵＩ部２００は、キーボードやマウス等からなる入力部２１０と、表示部（モニタ）やプリンタ等の出力部２２０とを備え、ユーザからの入力を受け付け、演算部４００による処理結果をユーザに提示する。表示部として、液晶ディスプレイやＣＲＴ（Ｃａｔｈｏｄｅ Ｒａｙ Ｔｕｂｅ）等を適用することができる他、表示部がタッチパネル機能を有し、入力装置２１０として機能するようにすることもできる。

【0015】

計測部３００は、後述する演算部４００のＣＰＵ４０１による制御に従って、被写体１０１にＸ線を照射し、被写体１０１を透過したＸ線フォトンを計測する。計測部３００は、Ｘ線照射部３１０、Ｘ線検出部３２０、ガントリ３３０、制御部３４０、データ収集部４０４、及び被写体１０１を載置する寝台１０２を備えている。

【0016】

ガントリ３３０の中央には、被写体１０１と、被写体１０１を載置する寝台１０２とを配置するための円形の開口部３３１が設けられている。ガントリ３３０の内部には、後述するＸ線管３１１及びＸ線検出器３２１を搭載する回転板３３２と、回転板３３２を回転させるための駆動機構が配置される。以下の説明において、開口部３３１の周方向をｘ方向、径方向をｙ方向、それらに直交する方向をｚ方向とする。一般にｚ方向は、被写体１０１の体軸方向となる。

【0017】

通常の撮像時は、回転板３３２が回転し、かつ寝台１０２がｚ方向に動くことで被写体１０１の様々な方向からＸ線を照射し、被写体１０１を透過したＸ線をＸ線検出器３２１で検出してデータを取得する。一方、スキャノグラム撮像時は、回転板３３２は回転せず、寝台１０２のみをｚ方向に１回移動させて終了する、もしくはその撮像の後Ｘ線の照射方向を変化（例えば９０度回転板３３２を回転）させた後撮像する、を規定回数くりかえし、ある一方向からの照射のデータ取得を角度を変えて複数回行ってデータを得る。

【0018】

回転板３３２の回転の所要時間は、ユーザがＵＩ部２００を介して入力したパラメータに依存する。本実施形態では、例えば、回転の所要時間を１．０ｓ／回とする。計測部３００による１回転における撮像回数は、例えば、９００回であり、回転板３３２が０．４度回転する毎に１回の撮像が行われる。なお、各仕様はこれらの値に限定されるものはなく、ＰＣＣＴ装置１００の構成に応じて種々変更可能である。

【0019】

Ｘ線照射部３１０は、Ｘ線を発生し、発生したＸ線を被写体１０１に照射する。Ｘ線照射部３１０は、Ｘ線管３１１と、Ｘ線フィルタ３１２と、ボウタイ（ｂｏｗｔｉｅ）フィルタ３１３と、を備えている。

【0020】

Ｘ線管３１１は、後述する照射制御器３４１の制御に従って供給される高電圧により、被写体１０１にＸ線ビームを照射する。照射されるＸ線ビームは、ファン角及びコーン角を持って広がる。Ｘ線ビームは、後述するガントリ３３０の回転板３３２の回転に伴って、被写体１０１に照射される。

【0021】

Ｘ線フィルタ３１２は、Ｘ線管３１１から照射されたＸ線のＸ線量を調節する。すなわち、Ｘ線のスペクトルを変化させる。本実施形態のＸ線フィルタ３１２は、Ｘ線管３１１から被写体１０１へ照射されるＸ線が、予め定めたエネルギー分布となるよう、Ｘ線管３１１から照射されたＸ線を減衰させる。Ｘ線フィルタ３１２は、被写体１０１である患者の被ばく量を最適化するために用いられる。このため、必要なエネルギー帯の線量が強くなるよう設計される。

【0022】

ボウタイフィルタ３１３は、周辺部の被ばく量を抑える。被写体１０１である人体の断面が楕円形であることを用い、中心付近の線量を強くし、周囲の線量を低くして被ばく量を最適化するために用いられる。

【0023】

Ｘ線検出部３２０は、Ｘ線フォトンが入射する毎に、当該Ｘ線フォトンのエネルギー値を計測可能な信号を出力する。Ｘ線検出部３２０は、Ｘ線検出器３２１を備えている。
図２に、Ｘ線検出器３２１の一部を例示する。本実施形態のＸ線検出器３２１は、複数の検出素子３２２、カウンティング回路３２４、及びＸ線検出器３２１への入射方向を制限するコリメータ３２３と、を備えている。
なお、製作を容易にするために平面状の検出器（検出器モジュール）を複数作成し、平面の中心部分が円弧になるように配置して疑似的に円弧状に配列し、Ｘ線検出器３２１としてもよい。

【0024】

各検出素子３２２に入射したＸ線は、カウンティング回路３２４により、１つのＸ線フォトンが入射する毎に、１パルスの電気信号（アナログ信号）に変換される。変換された電気信号は、データ収集部４０４に入力される。

【0025】

検出素子３２２には、入射したＸ線フォトンを、直接電気信号に変換する、例えば、ＣｄＴｅテルル化カドミウム（cadmium telluride）系の半導体素子を用いる。なお、検出素子３２２は、Ｘ線を受けて蛍光を発するシンチレータ（Scintillator）及び蛍光を電気に変換するフォトダイオードを用いてもよい。

【0026】

Ｘ線検出器３２１の検出素子３２２の数（チャンネル数）は、例えば、１０００個である。各検出素子の、ｘ方向のサイズは、例えば、１ｍｍである。
また、例えば、Ｘ線管３１１のＸ線発生点と、Ｘ線検出器３２１のＸ線入射面との距離は、１０００ｍｍである。ガントリ３３０の開口部３３１の直径は、７００ｍｍである。
なお、ガントリと同様に各仕様はこれらの値に限定されるものはなく、ＰＣＣＴ装置１００の構成に応じて種々変更可能である。

【0027】

制御部３４０は、Ｘ線管３１１からのＸ線の照射を制御する照射制御器３４１、回転板３３２の回転駆動を制御するガントリ制御器３４２、寝台１０２の駆動を制御する寝台制御器３４３、Ｘ線検出器３２１におけるＸ線検出を制御する検出制御器３４４を備えている。これら各部は、後述する演算部４００の計測制御部４２０による制御に従って動作する。

【0028】

演算部４００は、ＰＣＣＴ装置１００の動作全体を制御し、計測部３００で取得したデータを収集し、処理することにより撮像を行う。図１及び図３に示すように、演算部４００は、計測部３００に撮像を実行させて得られたデータの各種処理を行う中央処理装置（以下、「ＣＰＵ」という）４０１と、処理の際の作業領域となるメモリ４０２と、予め必要なプログラムや処理に用いるデータ、処理中に生成されるデータ、処理の結果得られるデータ等を記憶したＨＤＤ（Hard disk drive）装置４０３を備えている。なお、演算部４００による処理結果は、ＵＩ部２００の出力装置２２０にも出力される。

【0029】

ＣＰＵ４０１は、図３に示すように、撮像条件設定部４１０と、計測制御部４２０と、データ加算部４３０と、補正部４４０と、画像生成部４５０を含み、ＣＰＵ４０１に含まれる各部によって、被写体の位置決め画像を取得し、位置決め画像に基づいてＸ線照射量を決定する自動露光制御を行う。

【0030】

なお、ＣＰＵ４０１に含まれる各部の機能は、ＣＰＵ４０１が予めＨＤＤ装置４０３等の記憶装置に格納されたプログラムを読み込んで実行することによりソフトウエアとして実現することができる。また、ＣＰＵ４０１に含まれる各部が実行する動作の一部又は全部を、ＡＳＩＣ（application specific integrated circuit）やＦＰＧＡ（field-programmable gate array）等のハードウェアにより実現することもできる。

【0031】

撮像条件設定部４１０は、ユーザからＵＩ２００を介して入力された撮像条件を設定する。例えば、撮像条件設定部４１０は、撮像条件を受け付ける受付画面を表示装置２２０に表示させ、受付画面を介してユーザによる撮像条件の入力を受け付け、設定する。
撮像条件には、例えば、Ｘ線管３１１の管電流、管電圧、被写体１０１の撮像範囲、Ｘ線フィルタ３１２の形状、ボウタイフィルタ３１３の形状、分解能等が含まれる。なお、撮像条件は、必ずしも毎回ユーザが入力する必要はない。例えば、予め典型的な撮像条件をＨＤＤ装置４０３等に保存しておき、ＣＰＵ４０１が読み出して用いてもよい。

【0032】

計測制御部４２０は、撮像条件設定部４１０により設定された撮像条件に従って、制御部３４０を制御し計測を実行する。
以下、計測制御部４２０の制御に従って行われるＣＴ像撮像方法について具体的に説明する。
計測制御部４２０は、寝台制御器３４３に対し、寝台１０２を回転板３３２に対して垂直な方向に移動させ、回転板３３２の撮像位置が指定された撮像位置と一致した時点で移動を停止するように指示する。これにより、被写体１０１の配置が完了する。

【0033】

計測制御部４２０は、寝台制御器３４３への指示と同じタイミングで、ガントリ制御器３４２に対して駆動モーターを動作させ、回転板３３２に対して回転を開始するよう指示を行う。

【0034】

この指示により回転板３３２が回転を開始し回転板３３２の回転が定速状態になり、かつ被写体１０１の配置が終了すると、計測制御部４２０は、照射制御器３４１に対しＸ線管３１１のＸ線照射タイミングを指示すると共に、検出制御器３４４に対してＸ線検出器３２１の撮像タイミングを指示する。これにより、計測制御部４２０は、Ｘ線の照射及び検出器によるＸ線フォトンの検出を行い、計測を開始する。計測制御部４２０が、これらの指示を所定回数繰り返すことで計測部３００による撮像範囲全体の計測を行う。なお、公知のヘリカルスキャン（Helical Scan）のように、寝台１０２を移動させながら撮像を行うよう制御してもよい。

【0035】

次に、スキャノグラムの撮像方法を示す。計測制御部４２０は、寝台制御器３４３に対し寝台１０２を回転板３３２に対して垂直な方向に移動させ、回転板３３２の撮像位置が指定された撮像位置と一致した時点で移動を停止するように指示する。これにより、被写体１０１の配置が完了する。

【0036】

次に、計測制御部４２０は、寝台制御器３４３へ移動の指示のみを出し、ガントリ制御器３４２に対しては何も命令しない。よって、回転板３３２は回転せずに寝台１０２のみが移動する。この状態で計測制御部４２０は、照射制御器３４１に対しＸ線管３１１のＸ線照射タイミングを指示し、検出制御器３４４に対しＸ線検出器３２１の撮像タイミングを指示する。これにより、計測制御部４２０は、Ｘ線の照射及び検出器によるＸ線フォトンの検出を行い、計測を開始する。そして、寝台１０２の移動距離がユーザの指定した距離、もしくは寝台１０２の移動限界距離に達すると撮像を終了する。

【0037】

この作業を撮像と撮像の間に回転板３３２の角度を変化させてＸ線管３１１、及びＸ線検出器３２１の照射−計測方向を変化させて繰り返すことで複数の方向における透過データを取得することができる。このようにして検出器３２１が検出したＸ線フォトンに係る透過データは、データ収集部４０４により収集される。

【0038】

データ収集部４０４によるデータの収集は、例えば、スキャノグラムの場合以下のように行われる。上述のように、回転板３３２は静止した状態でＸ線を照射しながら寝台１０２を移動させる。その間、例えば０．１ｍｓ毎にＸ線検出器３２１で検出され、予め定めた複数のエネルギー範囲に弁別されたそれぞれのエネルギー幅のカウント数を得る。

【0039】

ここで、Ｘ線照射範囲が広い場合は、Ｘ線照射範囲を通過する間に複数データ収集が実施される。その場合は同じ位置のデータに関して平均値を使用することでデータの精度を高める。このような方法で収集したデータを保存、もしくは補正部４４０や画像生成部４５０の演算実施後に保存する。

【0040】

なお、必要に応じてＸ線のｚ方向の照射範囲を限定してもよい。ｚ方向の照射範囲を限定する方法としては例えばＸ線管への電圧印加をＯＮ／ＯＦＦする方法や、Ｘ線管の前にＸ線を減衰、吸収するシャッターを設けシャッターのＯＮ／ＯＦＦでＸ線を制御してもよい。また、回転板３３２を例えば９０度回転させるなどをし、別の角度からＸ線を照射したデータも作成してもよい。

【0041】

このように、データ収集部４０４は、エネルギー範囲毎に、Ｘ線フォトンの数を計数する。計数して得られる結果は、Ｘ線フォトンのエネルギー値（単位ｋｅＶ）の分布を示す。従って、データ収集部４０４は、これにより、Ｘ線検出器３２１で検出したＸ線のエネルギー分布（スペクトル）を得る。データ収集部４０４は、得られた結果を、計数情報としてＣＰＵ４０１に出力する。
データ収集部４０４のエネルギー範囲の詳細については後述する。

【0042】

データ加算部４３０は、データ収集部４０４で収集したエネルギー範囲毎のデータについて重みを付けて加算する。データ加算部４３０の詳細については後述する。
補正部４４０は、データ収集部４０４により収集されたデータに対し、所定の補正データに従って補正処理を実施する。ここで行う補正処理は、例えば、リファレンス補正回路のリニアリティ補正、対数変換処理、オフセット処理、感度補正、ビームハードニング補正、水ファントムキャリブレーション、ＣＴ値補正などである。

【0043】

補正データとしては、例えば各エネルギー範囲のデータから基準となる物質を弁別する際に必要となる、基準物質の各エネルギー範囲におけるＸ線減弱係数である。他にも、ファントムのサイズに応じて補正を行うファントムキャリブレーションを実施する場合の各エネルギー範囲における補正データやＣＴ値補正データなども設定したエネルギー範囲ごとに必要となる場合がある。補正データの作成は、たとえばＵＩ部２００から補正データ取得命令を実施して、水の入ったアクリル容器等の補正用ファントムを撮像したデータを、補正部４４０で補正データ作成用のプログラムにより計算することで求める。

【0044】

それらの補正データを予めエネルギー範囲の設定をスキャノグラム等から撮像前に予め設定することで撮像前での演算が可能となり、画像作成時間の短縮が可能となる。なお、これらの補正データは、撮像毎に生成することができる他、被写体１０１が変わっても変化しないデータである場合には、予め補正データを収集することもできる。

【0045】

画像生成部４５０は、データ加算部４３０で加算したデータから、Ｘ線ＣＴ画像を再構成する。画像は、例えば、Ｘ線フォトン数に対し、Ｌｏｇ変換を行い、再構成する。再構成には、ＦｅｌｄＫａｍｐ法、逐次近似法など、各種の公知の手法を用いることができる。

【0046】

データ収集部４０４は、Ｘ線検出器３２１が検出したＸ線に由来するフォトン（Ｘ線フォトン）を、予め定めた複数のエネルギー範囲についてエネルギー範囲毎に計数し、当該エネルギー範囲毎の計数情報を得る。

【0047】

データ収集部４０４は、Ｘ線検出器３２１が検出したＸ線フォトン１つ１つのエネルギー値を取得し、そのエネルギー値に応じてエネルギー範囲毎に設けられたエネルギービン（Ｂｉｎ）の計数結果に加算する。エネルギービンは、エネルギー範囲毎に設定される記憶領域である。なお、エネルギー範囲は、ある最小エネルギーからＸ線管３１１の最大エネルギーまでのエネルギー範囲を、所定のエネルギー幅ΔＢで区切ったものである。

【0048】

ここで、ＣＴ値を精度良く求めるにはエネルギー幅の設定が重要となる。エネルギー範囲が大きいほど誤差が発生する可能性が高くなる。レイトレースシミュレーションにより１６５ｃｍの５ｍｍ厚さのポリエチレン容器に入った水におけるＣＴ値の誤差を求めたところ、エネルギー範囲を１つに設定した場合は最大で約８ほどのＣＴ値のずれが発生したが、エネルギー範囲を５に設定したところ最大でも１以下となり、４に設定すると２以下となった。

【0049】

ＣＴ値は２以上ずれなければ臨床上問題ないと考えられているため、エネルギー範囲を４以上で設定すれば問題ないと考えられる。一方、エネルギー範囲を細かく設定しても検出器のエネルギー分解能以上に設定した場合は、エネルギーの弁別の正しさが問題となるため細かくしすぎる必要はない。

【0050】

ＣｄＴｅ半導体検出器を常温で用いた場合は数ｋｅＶのエネルギー分解能となる。エネルギー分解能は温度に依存することから一義的ではないものの、２ｋｅＶのエネルギー分解能で、管電圧が１２０ｋＶである場合、１２０/２より、エネルギー範囲は６０となる。しかし、６０以上のエネルギー範囲を使用しても画質の向上は期待できない上に、データの転送量は多くなることからシステムが複雑化する。そのため、本実施形態ではエネルギー範囲を４以上６０以下とすることが望ましいと考えられる。また、各エネルギー範囲において、上限値と下限値との差が検出器が検出可能なエネルギー範囲の半値幅より大きいことが望ましい。

【0051】

従って、本実施形態において、データ収集部４０４はエネルギー幅ΔＢを、例えば２０ｋｅＶとし、最小エネルギーを４０ｋｅＶ、最大エネルギーを１４０ｋｅＶとする。つまり、エネルギー範囲４０ｋｅＶ〜１４０ｋｅＶを、Ｂ１（４０〜６０ｋｅＶ）、Ｂ２（６０〜８０ｋｅＶ）、Ｂ３（８０〜１００ｋｅＶ）、Ｂ４（１００〜１２０ｋｅＶ）、Ｂ５（１２０〜１４０ｋｅＶ）の５つのエネルギー範囲に区分する。

【0052】

なお、エネルギー幅ΔＢはＢ１〜Ｂ５で同じ値でなくてもよいし、一部のエネルギーのデータを重複させるように設定したり、特定のエネルギーのデータを検出しないように設定したりすることも可能である。ＤＡＳは、検出したＸ線フォトンのエネルギー値に応じて、該当するエネルギー範囲に対応づけて設けられたエネルギービンの計数結果に加算する。

【0053】

（データ加算部４３０の詳細）
上述のように、ＣＴ値を求める場合において、ＰＣＣＴ装置では、特定のエネルギー幅のカウント数が得られるに留まり、詳細なエネルギー情報を得ることができないことから精度の高いＣＴ値を求めることが困難である。一方、フォトンカウンティング方式ではない検出器を適用したＸ線ＣＴ装置では、エネルギーの大きさに比例したＣＴ値を得ることができる。これは、実質的にＸ線のエネルギーに応じた重みがかかって加算された値ということができる。そこで、ＰＣＣＴ装置においても、データ加算部４３０が、データ収集部４０４から得られた各エネルギー範囲の計数結果に対して、Ｘ線の減衰率やエネルギー分布を考慮して適切な重みを乗じることでエネルギーの大きさに応じたＣＴ値を得ることができる。

【0054】

Ｘ線の減衰率（単位距離進む間に散乱、吸収されるＸ線の割合）はＫエッジなどの特殊な状況を除いてエネルギーが小さいほど高くなる。また、Ｋエッジの効果も限定的で総じてみればエネルギーが小さいほど高くなる。そのため、エネルギー方向に均一分布しているＸ線を照射した場合、被写体透過後は高エネルギー側のＸ線のほうが多くなる。これは、あるエネルギー範囲を切り取った時に平均エネルギーが被写体透過前に比べて被写体透過後の方が高くなることを意味する。この効果に鑑みて、重みをかけることでより正確なＣＴ値を得ることができる。

【0055】

さらに、Ｘ線管から照射されるＸ線も上記のようなエネルギー方向に均一分布しているわけではない。たとえば１２０ｋＶの管電圧の場合、Ｘ線のエネルギーは６０ｋｅＶ前後にピークを持つ分布となる。その分布により、高エネルギー側では被写体透過前の時点でエネルギー範囲の中央値よりも低エネルギー側にシフトし、低エネルギー側では中央値よりも高エネルギー側にシフトしている。これらの効果を含めて重みをかけることでさらに正確になる。

【0056】

フォトンカウンティング方式の検出器では１つの検出器の中で信号がすべて減衰すれば全Ｘ線エネルギーが該当検出器の電気信号に変換されるが、一部の信号が別の検出器で減衰したり、あるＸ線の信号のエネルギーを求めているときにさらにＸ線が入るなどして別のエネルギー範囲の信号として計測されることが発生する。これらの効果を含めて重みを計算すると上記よりもさらに正確なＣＴ値を得ることができる。

【0057】

データ加算部４３０により、上記した各特性を踏まえて決定した重みを乗じた加算処理を実施し、これを画像生成部４５０において再構成等の処理を行うことで、従来型のＸ線ＣＴ装置とほぼ同等の投影データを得ることが可能となる。そして、ここで得た投影データを用いて、従来と同等のデータ処理を実施することで、従来型のＸ線ＣＴ装置と同等のＣＴ値を得ることが可能となる。

【0058】

したがって、データ加算部４３０は、先ず、重みを計算する。例えば、重みとして、本実施形態では各エネルギー範囲における中央値を使う。具体的には、Ｂ１（４０〜６０ｋｅＶ）の場合には５０ｋｅＶであり、Ｂ２（６０〜８０ｋｅＶ）の場合には７０ｋｅＶとなる。重みは、このように予め設定した重みを適用することができる他、被写体に併せて重みを都度算出してもよい。

【0059】

続いて、このように構成されたＸ線ＣＴ装置における撮像処理について、図４のフローチャートに従って説明する。
まず、撮像条件設定部４１０が、ＵＩ部２００を介して、ユーザから撮像条件を受け付ける（ステップＳ１０１）。ここで入力を受け付ける撮像条件には、管電圧、管電流、Ｘ線フィルタ３１２の厚み、形状、ボウタイフィルタ３１３の形状などがある。
次にスキャノグラムを撮像する（ステップＳ１０２）。撮像は、データ収集部４０４の説明で示したように回転板を回転させずに収集する。

【0060】

次に、データ収集部４０４は、エネルギー範囲を設定又は変更する（ステップＳ１０３）。本実施形態においては上述したＢ１（４０〜６０ｋｅＶ）、Ｂ２（６０〜８０ｋｅＶ）、Ｂ３（８０〜１００ｋｅＶ）、Ｂ４（１００〜１２０ｋｅＶ）、Ｂ５（１２０〜１４０ｋｅＶ）の範囲を設定する。
次にこの設定に従って、補正部４４０及び画像生成部４５０において使用する補正データの作成を行う（ステップＳ１０４）。本実施形態ではエネルギー範囲を５つ設定しているため、これら５つの設定に必要な補正データを作成する。

【0061】

次に、計測制御部４２０は、ステップＳ１０２で設定された撮像条件に従って、撮像を実行し（ステップＳ１０５）、データ収集部４０４は、回転板を回転させたＣＴ撮像の方法に従ってデータを収集する。

【0062】

次に収集したデータを加算する（ステップＳ１０６）。それぞれのエネルギー範囲において設定された重みを乗算して加算する。本実施形態ではたとえば各エネルギー範囲の中間値を重みに用いる。なお、加算後のデータのみを保存してもよいし、必要に応じて加算前のデータも残しておいてもよい。

【0063】

その後、補正部４４０は、データ収集部４０４が収集したデータを補正する（ステップＳ１０７）。ここで、ステップＳ１０４において予め作成した補正データを用いて補正することで画像生成時間の短縮が可能となる。補正部４４０では、そのまま収集したカウントの情報で処理してもよいし、一旦基底物質の距離の情報に変換してから処理を実施してもよい。

【0064】

画像生成部４５０は、ステップＳ１０７において補正されたデータを用いて画像を生成する（ステップＳ１０８）ここでも、必要に応じてステップＳ１０４で作成した補正データを用いて画像生成を実施する。ステップＳ１０８で生成された画像をＨＤＤ装置４０３に保存し、例えば入力装置２１０のモニタに画像を表示することで撮像のシーケンスが完了する（ステップＳ１０９）。その後はユーザが画像を用いて診断し、必要に応じて画像解析を実施する。

【0065】

このように本実施形態によれば、データ加算部４３０において、各エネルギー範囲のカウント値に対する加算処理の際に重みを乗じ、その後再構成等の処理を行うことで、フォトンカウンティング方式の検出器を用いていない従来型のＸ線ＣＴ装置において得られる投影データと同等の投影データを得ることができる。従って、この投影データに従って所定の処理を行うことにより精度の高い所望のＣＴ値を得ることができる。

【0066】

なお、上述の例では、データ加算部４３０は、エネルギー範囲の中央値を重みとして使用したが、エネルギー範囲の間であれば別の値を使用することも可能である。例えば、そのエネルギー範囲の最大値を使用したり、最小値を使用したりすることも可能である。また、エネルギー範囲毎に使用する重みの条件を変化させてもよい。

【0067】

例えば、Ｂ４（１００〜１２０ｋｅＶ）、Ｂ５（１２０〜１４０ｋｅＶ）のエネルギー範囲では最小値の１００ｋｅＶ、１２０ｋｅＶをそれぞれ使用し、Ｂ３（８０−１００ｋｅＶ）では中央値の９０ｋｅＶを使用、そしてＢ１（４０〜６０ｋｅＶ）、Ｂ２（６０〜８０ｋｅＶ）ではそれぞれ最大値の６０ｋｅＶ、８０ｋｅＶを使ってもよい。

【0068】

（第２の実施形態）
上述した第１の実施形態に係るＸ線ＣＴ装置においては、各エネルギー範囲に対して、中央値を重みとして適用していたが、例えば、エネルギー範囲毎に、被写体の透過距離に合わせて重みを最適化させてもよい。

【0069】

本実施形態に係るＸ線ＣＴ装置は、上述した第１の実施形態におけるＸ線ＣＴ装置と同一の構成を有し、演算部４００に含まれる各部が行う処理が異なるのみである。従って、本実施形態に係るＸ線ＣＴ装置を構成する各部について、第１の実施形態におけるＸ線ＣＴ装置の各構成と同一の符号を付し、夫々詳細な説明は省略する。

【0070】

以下に説明する本実施形態では、被写体の透過距離に応じてエネルギー範囲毎に設定された重みを変化させる。つまり、本実施形態では、データ加算部４３０において、予め定めた重み（例えば、各エネルギー範囲における中央値）に対して、被写体の透過距離による変化を考慮した値を加味することで重みを補正し、補正された重みを適用して加算処理を行う。

【0071】

このため、本実施形態に係るＸ線ＣＴ装置は、予めＨＤＤ４０３に重みを補正するための補正係数を保持しており、補正係数に基づいて補正した重みを用いて加算処理を行う。予め定めた重みが、各エネルギー範囲の中央値である場合には、重みを補正するための補正係数として、例えば、エネルギー範囲の中央値に対するずれをＨＤＤ４０３に記憶しておく。

【0072】

ＨＤＤ４０３には、例えば、以下の手順で算出された、エネルギー範囲の中央値に対する「ずれ」が補正係数として記憶されている。
被写体１０１の撮像に先立って、いくつかの水ファントムを準備し撮像する。水ファントムはたとえばアクリルの容器の中に水を満たしたファントムで、ほとんどが水や水と近い密度の物質で構成される人体を模擬する。水ファントムは例えば直径１００ｍｍ、２００ｍｍ３００ｍｍ、４００ｍｍの４つの水ファントムを準備する。

【0073】

これらのファントムを撮像する際に、エネルギー範囲は各エネルギー範囲をより細かく設定して撮像する。例えば、Ｂ１（４０〜６０ｋｅＶ）において、４０−４２ｋｅＶ、４２−４４ｋｅＶ、…と２ｋｅＶ刻みで設定する。これを本実施形態では４つのファントム、５つのエネルギー範囲において２ｋｅＶ刻みで撮像する。

【0074】

次に、このようにして得られた撮像データに基づいて、各直径の水ファントムにおける各エネルギー範囲の中央値に対するずれを計算する。より具体的には、２ｋｅＶ刻みで設定したエネルギー範囲の夫々の計数値に基づいて、中央値に対するエネルギー範囲のずれを計算する。

【0075】

例えば、４０−６０ｋｅＶのエネルギー範囲のファントム透過後のエネルギーに対するずれは、次のように計算する。本実施形態では４０−４２ｋｅＶ、４２−４４ｋｅＶ、…と２ｋｅＶ刻みのエネルギー範囲に対して夫々、（（４０−４２ｋｅＶのカウント数）×４１ｋｅＶ＋（４２−４４ｋｅＶのカウント数）×４３ｋｅＶ＋...＋（５８−６０ｋｅＶのカウント数）×５９ｋｅＶ）/（４０−６０ｋｅＶのカウント数）として平均エネルギーを算出し、その値と中央値とのずれを計算する。
ファントムを透過すると高エネルギー側のＸ線ほど透過するため、エネルギー分布は通常高エネルギー側にシフトする。

【0076】

ここで、図５（ａ）に水の減弱計数を示す。１００ｍｍの水を透過した場合、全てのエネルギー範囲での入力Ｘ線の量を１とした場合に検出される量は図５（ｂ）のようになる。このように入力のエネルギースペクトルが、エネルギーによらず一定の場合は高エネルギー側にシフトする。その平均エネルギーを各ファントムについて計算する。

【0077】

例えば図５（ｂ）に示すようなカウントが得られた場合の平均エネルギーは５０．８ｋｅＶとなり、０．８ｋｅＶがずれ量となる。なお、本実施形態では実測して補正係数を求める方法を示したが、たとえばモンテカルロシミュレーションを用いたり、理論的な減衰量を用いたりして計算することもできる。

【0078】

以下、本実施形態に係るＸ線ＣＴ装置における撮像処理について、図６のフローチャートに従って説明する。

【0079】

まず、ステップＳ２０１において、撮像条件設定部４１０が、ＵＩ部２００を介して、ユーザから撮像条件を受け付け、ステップＳ２０２においてスキャノグラムを撮像する。ステップＳ２０３で、エネルギー範囲を設定又は変更し、設定されたエネルギー範囲に応じてステップＳ２０４で補正データの作成を行う。

【0080】

次のステップＳ２０５において、データ加算部４３０において加算処理を行う際に用いる重みを補正する。すなわち、データ加算部４３０は、ステップＳ２０２において撮像したスキャノグラムから被写体のサイズを算出し、その被写体のサイズに近い撮像データに対する補正係数をＨＤＤ４３０から読み込み、これを用いて重みを補正する。ここで、データ加算部４３０は、被写体のサイズに近い撮像データに対する補正係数を使用せずに、２つのファントムのデータの値を線形補間してもよいし、数点の補正係数を関数フィッティングすることにより重みを補正してもよい。

【0081】

ここで、データ加算部４３０において重みを補正する際に、ＨＤＤ４０３に記憶された補正係数を算出するときに適用されたエネルギー範囲がスキャノグラムを撮像したエネルギー範囲と変わっている場合には、以下のように補正係数を決定する。
例えば４０−５０ｋｅＶの範囲が設定された場合は、２ｋｅＶ刻みで撮像した４０−４２、４２−４４、...４８−５０の５つの撮像データの値からもう一度、対応範囲におけるファントム透過後のエネルギー範囲のずれを算出して補正係数とする。

【0082】

さらに、より被写体のサイズを正確に求めるためには、スキャノグラムのデータを用いずに、重みの補正に先立って本撮像を行い、本撮像によって得られたデータから被写体１０１のサイズを求めてもよい。

【0083】

なお、本実施形態における説明では、予め定めた重みをエネルギー範囲の中央値としているため、この中央値に対するずれを補正係数としたが、ずれ量に限られず、計算した平均エネルギー補正係数としてもよい。

【0084】

ステップＳ２０６では、ステップＳ２０２で設定された撮像条件に従って、計測制御部４２０により撮像を実行し、データ収集部４０４によってデータを収集する。収集されたデータはデータ加算部４３０に出力され、ステップＳ２０７において、ステップＳ２０５において、それぞれのエネルギー範囲において設定され補正された重みを乗算して加算される。

【0085】

その後、補正部４４０により、データ収集部４０４により収集されたデータが補正データに従って補正される（ステップＳ２０８）。なお、補正部４４０は、データ収集部４０４が収集したカウントのそのままの情報で処理してもよいし、一旦基底物質の距離の情報に変換してから処理を実施してもよい。

【0086】

画像生成部４５０は、ステップＳ２０８において補正されたデータを用いて画像を生成し（ステップＳ２０９）、生成された画像をＨＤＤ装置４０３に保存し、例えば入力装置２１０のモニタに画像を表示することで撮像のシーケンスを完了させる（ステップＳ２１０）。

【0087】

このように本実施形態によれば、データ加算部４３０において、各エネルギー範囲のカウント値に対する加算処理の際に重みを乗じ、その後再構成等の処理を行うことで、フォトンカウンティング方式の検出器を用いていない従来型のＸ線ＣＴ装置において得られる投影データと同等の投影データを得ることができる。従って、この投影データに従って所定の処理を行うことにより精度の高い所望のＣＴ値を得ることができる。そして、データ加算部４３０における加算処理の際に用いる重みを、被写体の大きさに伴うエネルギー範囲のずれに応じて補正することにより、より精度の高いＣＴ値を算出可能である。

【0088】

（第３の実施形態）
上述した第１の実施形態に係るＸ線ＣＴ装置においては、各エネルギー範囲の中央値を重みとして適用していたが、例えば、入射Ｘ線のエネルギー分布に合わせて重みを最適化させることもできる。

【0089】

本実施形態に係るＸ線ＣＴ装置は、上述した第１の実施形態及び第２の実施形態におけるＸ線ＣＴ装置と同一の構成を有し、演算部４００に含まれる各部が行う処理が異なるのみである。従って、本実施形態に係るＸ線ＣＴ装置を構成する各部について、第１の実施形態におけるＸ線ＣＴ装置の各構成と同一の符号を付し、夫々詳細な説明は省略する。

【0090】

以下に説明する本実施形態では、入射Ｘ線のエネルギー分布に合わせてエネルギー範囲毎に設定された重みを変化させる。つまり、本実施形態では、データ加算部４３０において、予め定めた重み（例えば、各エネルギー範囲における中央値）に対して、入射Ｘ線のエネルギー分布による変化を考慮した値を加味することで重みを補正し、補正された重みを適用して加算処理を行う。

【0091】

このため、本実施形態に係るＸ線ＣＴ装置は、予めＨＤＤ４０３に重みを補正するための補正係数を保持しており、補正係数に基づいて補正した重みを用いて加算処理を行う。予め定めた重みが、各エネルギー範囲の中央値である場合には、重みを補正するための補正係数として、例えば、エネルギー範囲の中央値に対するずれをＨＤＤ４０３に記憶しておく。

【0092】

ＨＤＤ４０３には、例えば、以下の手順で算出された、エネルギー範囲の中央値に対する「ずれ」が補正係数として記憶されている。
被写体１０１の撮像に先立って、ファントムのない状態で撮像を行う。撮像の際に、エネルギー範囲は各エネルギー範囲をより細かく設定して撮像する。例えば、Ｂ１（４０〜６０ｋｅＶ）において、４０−４２ｋｅＶ、４２−４４ｋｅＶ、…と２ｋｅＶ刻みで設定し、撮像を行う。そして、得られた撮像データを用いて各エネルギー範囲の中央値のずれを検出器毎に計算する。

【0093】

検出器毎に計算する理由は、主として二つある。１つ目は、Ｘ線の物理現象としてよく知られたヒール効果により被写体の体軸方向にエネルギースペクトルの変化があることである。また、２つ目は、Ｘ線管回転方向に対して被写体１０１への被ばくを低減するために、撮像中心は薄く、撮像視野の外側に行くほど厚い金属（たとえばアルミ製）を用い、あまり体を通らない部分に対してＸ線照射を抑制するフィルタを用いるため、体軸方向、Ｘ線管回転方向の双方にＸ線の分布が変化するからである。

【0094】

補正係数を算出するためにファントムがない状態で撮像を行っていることから、撮像によって得られたデータにより被写体１０１に入射する前のＸ線の各エネルギー範囲における中央値とのずれを求めることができる。

【0095】

一般的にＸ線管のエネルギー分布は、設定管電圧に依存するピークと、Ｘ線管のターゲット材質に起因する複数のピーク（材質によってはない場合も存在する）がある（図７参照）。設定管電圧に依存するピークより高い範囲は、中心から見て低エネルギー側にずれており、逆にピークより低い部分は、中心から見て高エネルギー側にずれている。
また、Ｘ線照射を抑制するフィルタによってもエネルギー分布が変化する。例えば、上述した水ファントムをおいた撮像のように、高エネルギー側が透過し、低エネルギー側がより散乱、吸収されるので、エネルギー分布が変化する。

【0096】

従って、上記した設定管電圧に依存するずれ及びフィルタ等によるずれ等の複合的な要因に加え、他のＸ線管やＸ線照射部分の構造に伴うＸ線のエネルギースペクトルの変化に応じた補正係数を求めることができる。なお、本実施形態では、ファントムがない状態で撮像を行うことにより、すなわち、実測により補正係数を求める方法を示したが、たとえばモンテカルロシミュレーションを用いたり、理論的な減衰量を用いたりして、計算により補正係数を算出することもできる。

【0097】

なお、被写体毎に重みを算出しなおす必要がない場合には、エネルギー分布を考慮して得られた補正係数を用いて補正した重みを予め記憶させておき、これをデータ加算部４３０において、各エネルギー範囲のカウント値に対する加算処理の際に適用することもできる。

【0098】

このように本実施形態によれば、データ加算部４３０において、各エネルギー範囲のカウント値に対する加算処理の際に重みを乗じ、その後再構成等の処理を行うことで、フォトンカウンティング方式の検出器を用いていない従来型のＸ線ＣＴ装置において得られる投影データと同等の投影データを得ることができる。従って、この投影データに従って所定の処理を行うことにより精度の高い所望のＣＴ値を得ることができる。そして、データ加算部４３０における加算処理の際に用いる重みを、エネルギー分布に起因するずれに応じて補正することにより、より精度の高いＣＴ値を算出することができる。

【0099】

なお、エネルギー分布に起因するずれを考慮した補正計数に加えて、上述した第２の実施形態における被写体の透過距離に応じた補正係数をＨＤＤ４０３に記憶させておくこともできる。この場合、エネルギー分布を考慮しながら、且つ、被写体の透過距離に合わせて重みを補正することができるので、重みをより最適化し、精度の高いＣＴ値を算出することができる。

【0100】

（第４の実施形態）
上述したように、加算処理時に適用する重みについて、エネルギー分布に起因するずれを考慮したり、被写体の透過距離を考慮する等によって最適化させることができる他、検出器における検出の揺らぎを考慮して最適化させることもできる。

【0101】

検出器での検出の揺らぎとは、例えば、フォトンカウンティング方式の検出器におけるパイルアップなどがある。すなわち、フォトンカウンティング方式の検出器では、Ｘ線の入射量が多くなると、パイルアップにより数え落としが発生するため、これを補正することが好ましい。

【0102】

このため、ＨＤＤ４０３に予めパイルアップを考慮した補正係数を記憶させておく。補正係数は以下のように算出する。
被写体の撮像前に、ファントムの無い状態で撮像を行う。この撮像時のエネルギー範囲は、上述した第２の実施形態と同様に、各エネルギー範囲をより細かく設定して撮像する。例えば、Ｂ１（４０〜６０ｋｅＶ）において、４０−４２ｋｅＶ、４２−４４ｋｅＶ、…と２ｋｅＶ刻みで設定し、管電流を変化させて撮像する。

【0103】

次に、この撮像により得られた撮像データを用いて各エネルギー範囲の中央値のずれを検出器ごとに算出し、算出されたずれを補正係数としてＨＤＤ４０３に記憶しておく。検出器ごとにずれを算出する理由は、検出器特性にばらつきがあることが懸念されるためである。管電流が増大してＸ線入射量が増えると、カウント数が正比例で増加する。しかし、パイルアップが発生すると、正比例の関係から逸脱する。この正比例からの逸脱に伴い、各エネルギー範囲の信号量が減少する。この分布の変化を計測し、その減少した信号量を補正することでより正確な補正が可能となる。

【0104】

パイルアップによる揺らぎ、すなわち、数え落としを補正する場合は、撮像データから単位時間当たりのカウント数を計算し、その結果を用いて補正を実施する必要がある。このため、重みの補正は、撮像後に行われる。

【符号の説明】

【0105】

１００・・・ＰＣＣＴ装置、１０１・・・被写体、１０２・・・寝台、２００・・・ＵＩ部、２１０・・・入力装置、２２０・・・出力装置、３００・・・計測部、３１０・・・Ｘ線照射部、３１１・・・Ｘ線管、３１２・・・Ｘ線フィルタ、３１３・・・ボウタイフィルタ、３２０・・・Ｘ線検出部、３２１・・・Ｘ線検出器、３２２・・・検出素子、３２３・・・コリメータ、３２４・・・カウンティング回路、３２５・・・積分回路、３３０・・・ガントリ、３３１・・・開口部、３３２・・・回転板、３４０・・・制御部、３４１・・・照射制御器、３４２・・・ガントリ制御器、３４３・・・テーブル制御器、３４４・・・検出制御器、４００・・・演算部、４０１・・・中央処理装置、４０２・・・メモリ、４０３・・・ＨＤＤ装置、４０４・・・データ収集部、４１０・・・撮像条件設定部、４２０・・・計測制御部、４３０・・・データ加算部、４４０・・・補正部、４５０・・・画像生成部

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】