(57)【要約】  工数や保持するデータ量を増大させることなく高品質な画像と低被ばく化とを両立させる。  代表的な撮影条件（以下、代表撮影条件とする）で最適化した逐次修正パラメータのテーブルを保持し、当該テーブルに保持される逐次修正パラメータから、実際の撮影条件において最適化した逐次修正パラメータを決定する。また、パラメータのテーブルとともに、基準とする重みも保持し、画素位置毎の逐次修正パラメータ生成時に当該基準重みも反映する。
2016199716.tif 000002
【特許請求の範囲】
【請求項１】
設定した撮影の条件に応じてＸ線を発生するＸ線発生部と、
被写体を透過後の前記Ｘ線を複数の検出素子で検出し、測定投影データを得るＸ線検出部と、
前記測定投影データから生成したＣＴ画像から順投影計算で求めた計算投影データと前記測定投影データとが等しくなるように、前記ＣＴ画像に対し逐次修正を行う逐次再構成部を有する画像生成部と、を備え、
前記逐次再構成部は、
代表的な撮影条件である代表撮影条件で取得した前記ＣＴ画像について、前記逐次修正によりノイズまたはＸ線量を低減させる割合である低減割合と、当該低減割合を実現する前記逐次修正に用いるパラメータである基準パラメータとの関係をテーブルとして保存するテーブル部と、
実際の撮影時の撮影条件である実撮影条件における所望の前記低減割合を、前記代表撮影条件における低減割合である代表低減割合に変換し、前記テーブルを参照して当該所望の低減割合を実現する前記逐次修正パラメータを決定するパラメータ決定部と、
前記パラメータ決定部で決定した前記逐次修正パラメータを用いて前記逐次修正を行う逐次修正部と、を備えること
を特徴とするＸ線ＣＴ装置。
【請求項２】
請求項１記載のＸ線ＣＴ装置であって、
前記パラメータ決定部は、
前記代表撮影条件と前記実撮影条件との差異による、前記実撮影条件で取得したＣＴ画像のノイズに対する前記代表撮影条件で取得したＣＴ画像のノイズの比であるノイズ増減比を用いて、前記所望の低減割合を前記代表低減割合に変換する代表ノイズ低減割合変換部と、
前記変換後の代表低減割合に対応する基準パラメータを前記テーブルから抽出する計算テーブル参照部と、
前記検出素子毎の重みに基づいて、ＣＴ画像上の画素毎の重みを計算する重み計算部と、
前記画素毎の重みと前記基準パラメータとから、画素毎の逐次修正パラメータを決定するパラメータ変換部と、を備えること
を特徴とするＸ線ＣＴ装置。
【請求項３】
請求項２記載のＸ線ＣＴ装置であって、
前記検出素子毎の重みに基づいて、ＣＴ画像上の予め定めた基準とする領域の画素の重み値を基準重みとして算出する基準重み計算部をさらに備え、
前記パラメータ変換部は、前記基準重みをさらに用い、前記画素毎の逐次修正パラメータを決定すること
を特徴とするＸ線ＣＴ装置。
【請求項４】
請求項２記載のＸ線ＣＴ装置であって、
ユーザが設定する撮影条件と前記代表撮影条件とは、再構成フィルタが異なり、
前記逐次再構成部は、前記再構成フィルタ毎に、同一の対象を前記代表撮影条件で撮影した際の前記ＣＴ画像のノイズを保存するフィルタノイズテーブルをさらに備え、
前記代表ノイズ低減割合変換部は、前記フィルタノイズテーブルを参照し、前記所望の低減割合を前記代表低減割合に変換すること
を特徴とするＸ線ＣＴ装置。
【請求項５】
請求項２記載のＸ線ＣＴ装置であって、
ユーザが設定する撮影条件と前記代表撮影条件とは、測定投影データに対する補正法が異なり、
前記逐次再構成部は、前記補正法毎に、基準の補正法に対するノイズ増減比を保存する補正法ノイズテーブルをさらに備え、
前記代表ノイズ低減割合変換部は、前記補正法ノイズテーブルを参照し、前記所望の低減割合を前記代表低減割合に変換すること
を特徴とするＸ線ＣＴ装置。
【請求項６】
請求項２記載のＸ線ＣＴ装置であって、
ユーザが設定する撮影条件と前記代表撮影条件とは、前記ＣＴ画像生成に用いる前記測定投影データの範囲および逐次修正に用いる前記測定投影データの範囲が異なり、
前記代表ノイズ低減割合変換部は、前記ユーザが設定する撮影条件および前記代表撮影条件それぞれにおける前記ＣＴ画像生成時に用いる測定投影データ範囲と逐次修正に用いる前記測定投影データ範囲とを用い、前記所望の低減割合を前記代表低減割合に変換すること
を特徴とするＸ線ＣＴ装置。
【請求項７】
請求項１記載のＸ線ＣＴ装置であって、
前記逐次再構成部は、前記測定投影データから生成したＣＴ画像の適否を判別する適否判別部をさらに備え、
前記適否判別部は、前記逐次修正において、前記ＣＴ画像が、予め定めた繰り返し回数内で収束しないと判別した場合、不適と判別すること
を特徴とするＸ線ＣＴ装置。
【請求項８】
請求項７記載のＸ線ＣＴ装置であって、
前記逐次再構成部は、前記適否判別部による判別結果が不適な場合、ユーザに通知する判別結果通知部、および、前記適否判別部による判別結果が不適な場合、前記ＣＴ画像を変更する画像変更部の少なくとも一方をさらに備えること
を特徴とするＸ線ＣＴ装置。
【請求項９】
請求項２記載のＸ線ＣＴ装置であって、
前記検出素子毎の重みの値は、一定および各前記検出素子で検出した光子数に比例した値のいずれかであること
を特徴とするＸ線ＣＴ装置。
【請求項１０】
設定した撮影の条件に応じてＸ線を発生するＸ線発生部と、被写体を透過後の前記Ｘ線を複数の検出素子で検出し、測定投影データを得るＸ線検出部と、前記測定投影データから生成したＣＴ画像から順投影計算で求めた計算投影データと前記測定投影データとが等しくなるように、前記ＣＴ画像に対し逐次修正を行う逐次再構成部を有する画像生成部と、を備えるＸ線ＣＴ装置における、前記逐次修正に用いる逐次修正パラメータ決定方法であって、
実際の撮影条件における、前記逐次修正による前記ＣＴ画像のノイズまたはＸ線量を低減させる割合である低減割合を、予め定めた代表的な撮影条件である代表撮影条件における前記低減割合に換算し、
予め保存される、前記代表撮影条件における前記低減割合と、当該低減割合を実現するために前記逐次修正に用いるパラメータである基準パラメータとの関係を参照し、前記逐次修正パラメータを決定すること
を特徴とする逐次修正パラメータ決定方法。
【請求項１１】
設定した撮影の条件に応じてＸ線を発生するＸ線発生部と、
被写体を透過後の前記Ｘ線を複数の検出素子で検出し、測定投影データを得るＸ線検出部と、
前記測定投影データから生成したＣＴ画像から順投影計算で求めた計算投影データと前記測定投影データとが等しくなるように、逐次修正を行う逐次再構成部を有する画像生成部と、を備え、
前記逐次再構成部は、
代表的な撮影条件である代表撮影条件で取得した前記ＣＴ画像について、前記逐次修正によりノイズまたはＸ線量を低減させる割合である低減割合と、当該低減割合を実現する前記逐次修正に用いるパラメータである基準パラメータとの関係をテーブルとして保存するテーブル部と、
実際の撮影時の撮影条件である実撮影条件における所望の前記低減割合を、前記代表撮影条件における低減割合である代表低減割合に変換し、前記テーブルを参照して当該所望の低減割合を実現する前記逐次修正パラメータを決定するパラメータ決定部と、
前記パラメータ決定部で決定した前記逐次修正パラメータを用いて前記逐次修正を行う逐次修正部と、を備えること
を特徴とするＸ線ＣＴ装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、Ｘ線ＣＴ（Ｃｏｍｐｕｔｅｄ  Ｔｏｍｏｇｒａｐｈｙ）装置に係り、特に、ＣＴ画像を逐次的に修正する逐次再構成技術に関する。
【背景技術】
【０００２】
Ｘ線ＣＴ装置は、被写体を多方向から撮影して得た測定投影データから被写体内の各点のＸ線吸収係数（ＣＴ値）を算出し、その分布画像（ＣＴ画像）を得る。被写体の低被曝化と高画質化とを両立させる技術として、逐次再構成法と呼ばれる手法がある（例えば、特許文献１参照）。逐次再構成法では、Ｘ線ＣＴ装置で得た測定投影データと、当該測定データから生成したＣＴ画像から算出した計算投影データとが等しくなるよう、計算投影データまたはＣＴ画像を逐次的に修正することにより、低線量でＣＴ画像のノイズを低減する。【０００３】
逐次再構成法は、従来の解析的にＣＴ値を計算する手法と比較して、反復更新に伴い計算量が増加する。また、膨大なパラメータの最適化等が必要となる。特に、更新中に隣接画素間を平滑化処理してノイズの低減効果を得る。この際、多くのパラメータの設定を必要とする。【０００４】
この計算量を低減する手法として、逐次修正中に出力されるＣＴ画像のノイズ、ここではＣＴ値のばらつきを表す標準偏差（Ｓｔａｎｄａｒｄ  Ｄｅｖｉａｔｉｏｎ、以下、ＳＤとする）の計測値が所望の値に到達した時点で修正を打ち切ったり、パラメータを変更したりする手法がある（例えば、特許文献２参照）。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００６−２５８６８号公報
【特許文献２】特許第４５３５７９５号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
特許文献２に開示の技術によれば、逐次修正中のＣＴ画像においてＳＤを計測する必要がある。ＳＤは、一様なＣＴ値で構成される組織に関心領域（Ｒｅｇｉｏｎ  ｏｆ  Ｉｎｔｅｒｅｓｔ、以下、ＲＯＩとする）を設定して計測する必要があるため、撮影対象が異なるＣＴ値で構成される組織の場合、正確に計測することが困難である。【０００７】
所望のＳＤとするため、撮影前に予め想定されるＸ線管や寝台速度等の撮影条件に対して、最適化した再構成パラメータをテーブルに保有する方法が挙げられる。しかし、撮影条件の数は膨大であるため、テーブルの作成に多大な工数がかかる。また、パラメータを保有するテーブル数も増加する。【０００８】
本発明は、上記事情に鑑みて成されたもので、工数や保持するデータ量を増大させることなく高品質な画像と低被ばく化とを両立させることを目的とする。
【課題を解決するための手段】
【０００９】
本発明は、代表的な撮影条件（以下、代表撮影条件とする）で最適化した逐次修正パラメータのテーブルを保持し、当該テーブルに保持される逐次修正パラメータから、実際の撮影条件において最適化した逐次修正パラメータを決定する。また、パラメータのテーブルとともに、基準とする重みも保持し、画素位置毎の逐次修正パラメータ生成時に当該基準重みも反映する。
【発明の効果】
【００１０】
本発明によれば、工数や保持するデータ量を増大させることなく高品質な画像と低被ばく化とを両立できる。
【図面の簡単な説明】
【００１１】
【図１】第一の実施形態のＸ線ＣＴ装置各部のハードウェアの構成のブロック図である。
【図２】第一の実施形態のＸ線ＣＴ装置の機能ブロック図である。
【図３】第一の実施形態の撮影条件受付画面例を説明するための説明図である。
【図４】第一の実施形態の逐次再構成部の機能ブロック図である。
【図５】第一の実施形態の逐次修正処理のうち、パラメータ決定部およびテーブル部による処理のフローチャートである。
【図６】（ａ）は、第一の実施形態のテーブル計算部による処理を説明するための説明図であり、（ｂ）は、第一の実施形態の基準パラメータ毎のノイズ低減割合の変化の一例を説明するためのグラフである。
【図７】第一の実施形態にパラメータテーブルの一例を説明するための説明図である。
【図８】（ａ）は、第一の実施形態の検出素子の重みを、（ｂ）は、基準重み計算部による計算結果を、それぞれ説明するための説明図である。
【図９】（ａ）は、第一の実施形態の、逐次修正処理における更新回数とノイズとの関係を説明するための説明図である。（ｂ）は、第一の実施形態のノイズテーブルの一例を説明するための説明図である。
【図１０】（ａ）は、第一の実施形態の重み画像例を、（ｂ）は、第一の実施形態の各画素位置の逐次修正パラメータ例を、それぞれ説明するための説明図である。
【図１１】第一の実施形態の逐次修正部による逐次修正処理のフローチャートである。
【図１２】（ａ）〜（ｃ）は、第一の実施形態の実施例の結果を説明するための説明図である。
【図１３】第二の実施形態の逐次修正処理のうち、パラメータ決定部およびテーブル部による処理のフローチャートである。
【図１４】第二の実施形態の補正法ノイズテーブルの一例を説明するための説明図である。
【図１５】（ａ）および（ｂ）は、第二の実施形態のノイズ増減比計算手法を説明するための説明図である。
【図１６】第三の実施形態の逐次修正処理のうち、パラメータ決定部およびテーブル部による処理のフローチャートである。
【図１７】第四の実施形態の逐次再構成部の機能ブロック図である。
【図１８】第四の実施形態の代表ノイズ低減割合変換部の処理を説明するための説明図である。
【図１９】本発明の実施形態の変形例の重み画像を説明するための説明図である。
【発明を実施するための形態】
【００１２】
＜＜第一の実施形態＞＞
本発明の第一の実施形態を、図面を用いて説明する。以下、本発明の実施形態を説明するための全図において、基本的に同一機能を有するものは同一符号を付し、その繰り返しの説明は省略する。
【００１３】
［Ｘ線ＣＴ装置］
まず、本実施形態のＸ線ＣＴ装置について説明する。Ｘ線ＣＴ装置は、被写体を多方向から撮影して得た測定投影データから被写体内の各点のＸ線吸収係数を算出し、Ｘ線吸収係数分布画像を得る装置である。各点のＸ線吸収係数は、空気を−１０００、水を０として規格化される。以後、規格化されたＸ線吸収係数をＣＴ値と呼び、Ｘ線吸収係数分布画像をＣＴ画像と呼ぶ。
【００１４】
図１に、本実施形態のＸ線ＣＴ装置１００のハードウェア構成を示す。また、図２に、本実施形態のＸ線ＣＴ装置１００のソフトウェアが実現する機能のブロック図を示す。
【００１５】
これらの図に示すように、本実施形態のＸ線ＣＴ装置１００は、入力部１１０と、撮影部１２０と、画像生成部１３０とを備える。
【００１６】
［入力部］
入力部１１０は、撮影条件を入力する撮影条件入力部２１１として機能する。本実施形態の撮影条件入力部２１１は、モニタに撮影条件受付画面を表示し、当該画面を介して撮影条件の入力を受け付け、受け付けた撮影条件を設定する。本実施形態の撮影条件受付画面および当該画面を介して受け付ける撮影条件の詳細は後述する。
【００１７】
また、図１に示すように、入力部１１０は、汎用のコンピュータが備えるハードウェア構成を有し、入出力インタフェースであるキーボード１１１やマウス１１２、データを記憶するメモリ１１３やＨＤＤ（Ｈａｒｄ  Ｄｉｓｋ  Ｄｒｉｖｅ）装置１１５、演算処理を行う中央処理装置（ＣＰＵ）１１４等を備える。さらに、入力部１１０は、図示を省略したモニタを備える。入力部１１０内の各構成要素はデータバス１１６によって接続される。
【００１８】
キーボード１１１およびマウス１１２は、撮影条件の入力等を行うために用いられる。なお、撮影条件の入力用として、ペンタブレットやタッチパネル等の他の入力手段を備えていてもよい。キーボード１１１等により入力されたデータは、処理部であるＣＰＵ１１４に受け渡される。
【００１９】
ＣＰＵ１１４は、メモリ１１３、ＨＤＤ装置１１５等に予め格納されている所定のプログラムを展開・起動することにより、撮影条件入力部２１１として機能する。また、ＣＰＵ１１４は、別のプログラムを展開・起動することにより、撮影部１２０に制御信号を送り、撮影部１２０の撮影制御部２２１の一部としても機能する。
【００２０】
なお、ＣＰＵ１１４が実現する全部または一部の機能は、ＡＳＩＣ（Ａｐｐｌｉｃａｔｉｏｎ  Ｓｐｅｃｉｆｉｃ  Ｉｎｔｅｇｒａｔｅｄ  Ｃｉｒｃｕｉｔ）、ＦＰＧＡ（ｆｉｅｌｄ−ｐｒｏｇｒａｍｍａｂｌｅ  ｇａｔｅ  ａｒｒａｙ）などのハードウェアによって実現してもよい。また、各機能の処理に用いる各種のデータ、処理中に生成される各種のデータは、メモリ１１３、ＨＤＤ装置１１５等に格納される。
【００２１】
［撮影部］
撮影部１２０は、図２に示すように、撮影条件入力部２１１を介して入力された撮影条件に基づいて、各部を制御して撮影を実行する撮影制御部２２１と、Ｘ線の照射および検出を行う信号取得部２２２として機能する。撮影部１２０は、これらの機能により、受け付けた撮影条件に応じてＸ線撮影を行い、測定投影データを得る。
【００２２】
信号取得部２２２は、図１に示すように、設定された撮影条件に応じてＸ線を発生するＸ線発生部１２１と、被写体１０１を透過後のＸ線を検出し、測定投影データを得るＸ線検出部１２２と、Ｘ線発生部１２１とＸ線検出部１２２とを搭載して被写体１０１の周囲を回転する機構とを備える。この機構は、ガントリ１２３と、Ｘ線発生部１２１とＸ線検出部１２２とを搭載して被写体１０１の周囲を回転する回転板１２４と、被写体１０１が載置される寝台１２５が挿入される円形の開口部１２６と、を備える。
【００２３】
Ｘ線発生部１２１とＸ線検出部１２２は、一般的なＸ線ＣＴ装置と同様、被写体１０１へのＸ線の照射および被写体１０１を透過したＸ線の検出を実現する。Ｘ線検出部１２２は、シンチレータ及びフォトダイオード等から構成される公知のＸ線検出素子（以下、単に検出素子と呼ぶ）を含み、複数の検出素子が、チャネル方向、すなわち、回転板１２４の主平面に平行な面内でＸ線発生部１２１から等距離の円弧に沿った方向、およびスライス方向、すなわち、被写体１０１の体軸方向にそれぞれ配列されている。
【００２４】
Ｘ線発生部１２１のＸ線発生点とＸ線検出部１２２のＸ線入力面との距離の代表例は、１０００［ｍｍ］である。開口部１２６の直径の代表例は、７００［ｍｍ］である。回転板１２４の１回転の所要時間の代表例は、１．０［ｓ］である。
【００２５】
例えば、チャネル方向の検出素子の数（以下、チャネル数とする）は、１０００個である。各検出素子のチャネル方向のサイズの代表例は、１［ｍｍ］である。回転板１２４が１回転する間の撮影回数は、９００回であり、回転板１２４が０．４度回転する毎に１回の撮影が行われる。撮影を行う際の回転板１２４の角度を、投影角と呼ぶ。なお、各仕様は、上記の値に限定されるものはなく、Ｘ線ＣＴ装置１００の構成に応じて種々変更可能である。
【００２６】
また、撮影制御部２２１は、図１に示すように、回転板１２４の回転動作を制御するガントリ制御器１２７と、Ｘ線発生部１２１およびＸ線検出部１２２の動作を制御するＸ線制御器１２８と、寝台１２５の位置を制御する寝台制御器１２９と、を備える。ガントリ制御器１２７、Ｘ線制御器１２８、および、寝台制御器１２９は、ＣＰＵ１１４の指示により、各部を制御する。
【００２７】
ＣＰＵ１１４は、ユーザから撮影開始の指示を受け付けると、ガントリ制御器１２７に指示をだし、回転板１２４の回転を開始させる。回転板１２４は、駆動モータにより駆動される。回転板１２４の回転が定速状態に入り、かつ、被写体１０１の撮影位置への配置が終了した時点で、ＣＰＵ１１４は、Ｘ線制御器１２８にＸ線の照射タイミングを指示し、Ｘ線検出部１２２にＸ線検出タイミングを指示する。
【００２８】
次に、図２の撮影部１２０は、撮影条件入力部２１１が受け付けた撮影条件に応じたＸ線撮影を行う。ユーザがマウス１１２やキーボード１１１等を用いて撮影開始を指示すると、ＣＰＵ１１４は、撮影制御部２２１の寝台制御器１２９およびガントリ制御器１２７に出力する。寝台制御器１２９は、制御信号を受けて、寝台１２５を回転板１２４の回転軸方向に移動させる制御を行い、被写体１０１の撮影部位が、Ｘ線発生部１２１とＸ線検出部１２２の間のＸ線通過範囲、すなわち、撮影位置に一致した時点で、寝台１２５の移動を停止させる。これにより被写体１０１の撮影位置への配置が完了する。
【００２９】
また、ガントリ制御器１２７は、ＣＰＵ１１４から撮影開始が指示されると同時に駆動モータを介して回転板１２４の回転を開始させる。回転板１２４の回転が定速状態に入り、かつ、被写体１０１の撮影位置への配置が終了した時点で、ＣＰＵ１１４は、Ｘ線制御器１２８にＸ線発生部１２１のＸ線照射タイミング、及び、Ｘ線検出部１２２の撮影タイミングを指示する。Ｘ線制御器１２８は、この指示に従ってＸ線発生部１２１からＸ線を照射させ、Ｘ線検出部１２２は、Ｘ線を検出して撮影を開始する。また、Ｘ線制御器１２８は、例えばユーザが設定したＸ線発生部１２１の管電圧および管電流時間積により、照射するＸ線のエネルギースペクトルと出力量を決定する。
【００３０】
なお、ここでは１種類のエネルギースペクトルを有するＸ線を使用する例について説明したが、本実施形態の構成はマルチエネルギーＣＴにも適用できる。その場合には、例えば、１回転毎または１回転中に管電圧を高速に切り替えて２種類以上のエネルギースペクトルを有するＸ線を照射し、撮影データを取得するように制御する。
【００３１】
また、本実施形態のＸ線ＣＴ装置１００において、被写体とは撮影対象を意味し、被写体１０１と、被写体１０１を支える寝台１２５とを包含する。なお、被写体１０１は、人体に限らず、ファントムや機械等の検査対象の物体であってもよい。
【００３２】
［画像生成部］
画像生成部１３０は、撮影部１２０で検出した信号（取得した測定投影データ）からＣＴ画像を再構成する。ＣＴ画像は、被写体の断層面を体軸方向に重ねて表現される。医療現場において、ＣＴ画像の利用は、正確かつ即時に患者の病状を診断できるため、臨床上有用である。しかし、医師の診断に必要な高い画質を得る条件において、被写体は一定量の被曝を伴う。
【００３３】
低被曝化を実現するためにＸ線量を低くすると、検出した信号に対するノイズの比率が増加する。これにより、誤診断の原因になるライン状のストリークアーチファクトや粒状性のノイズが多く発生する。そのため、低線量撮影時にストリークアーチファクトやノイズを低減することで、画像の高品質化と低被曝化を両立させることが望まれている。
【００３４】
本実施形態の画像生成部１３０は、これを実現するため、逐次再構成法により、最終的な出力画像を生成する。
【００３５】
これを実現するため、本実施形態の画像生成部１３０は、図１に示すように、データ収集システム（Ｄａｔａ  Ａｃｑｕｉｓｉｔｉｏｎ  Ｓｙｓｔｅｍ、以下、ＤＡＳ）１３１と、データを記憶するメモリ１３２と、ＨＤＤ装置１３４と、演算処理を行う中央処理装置（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ、以下、ＣＰＵ）１３３と、処理結果等を表示する表示部（モニタ）１３５とを備える。これらはデータバス１３６によって接続される。
【００３６】
画像生成部１３０は、まず、測定投影データからＣＴ画像（初期画像）を生成する。そして、その初期画像から順投影計算で算出する計算投影データと測定投影データとが等しくなるように、初期画像を逐次再構成（逐次修正）する。
【００３７】
これを実現するため、本実施形態の画像生成部１３０は、信号取得部２２２のＸ線検出部１２２が検出した信号をディジタル信号に変換するＡＤ変換を行う信号収集部２３１と、ディジタル信号に変換された測定投影データを補正する補正処理部２３２と、補正後の測定投影データからＣＴ画像を再構成する逐次再構成部２３３と、再構成したＣＴ画像を出力する画像表示部２３４として機能する。
【００３８】
画像生成部１３０の信号収集部２３１は、Ｘ線検出部１２２の出力信号をディジタル信号に変換し、メモリ１３２に保存する。この信号収集部２３１は、ＤＡＳ１３１により実現される。すなわち、撮影部１２０のＸ線検出部１２２で検出された信号は、信号収集部２３１として機能するＤＡＳ１３１によって収集されて、ディジタル信号に変換され、ＣＰＵ１３３に受け渡される。
【００３９】
補正処理部２３２は、信号収集部２３１による処理後の測定投影データに対し、Ｘ線の検出信号のゼロ値を較正するオフセット補正や、投影角度毎に検出した信号成分のばらつきを補正するリファレンス補正や、検出素子間の感度を補正する公知のエアキャリブレーション処理等の補正を行う。補正後の測定投影データは、逐次再構成部２３３に送られる。
【００４０】
このように，Ｘ線検出部１２２が得た測定投影データは、逐次修正前に補正処理部２３２において補正処理を施される。以後、補正後の測定投影データを、測定投影データと呼ぶ。
【００４１】
逐次再構成部２３３は、測定投影データからＣＴ画像を生成する。このとき、測定投影データから生成したＣＴ画像から順投影計算で求めた計算投影データと測定投影データとの差を低減させるよう、ＣＴ画像に対し逐次修正を行う。本実施形態では、ＣＴ画像の逐次修正は、所望のノイズ低減割合あるいはＸ線量の低減割合を実現するよう行う。ノイズ低減割合あるいはＸ線量の低減割合は、逐次修正時に用いるパラメータ（逐次修正パラメータ）の選択により制御される。
【００４２】
逐次再構成部２３３は、逐次修正において、測定投影データと計算投影データとの差を用いて当該差が小さくなるようにＣＴ画像を修正する計算（第一の計算）と、修正前のＣＴ画像の２以上の画素間のＣＴ値差を用いてＣＴ値差が小さくなるようにＣＴ画像を修正する計算（第二の計算）とを繰り返し行う。
【００４３】
第一の計算として、例えば、Ｌｉｋｅｌｉｈｏｏｄ計算、またはＤａｔａｆｉｔ計算を行う。以下、両者を代表してＬｉｋｅｌｉｈｏｏｄ計算と呼ぶ。また、第二の計算として、Ｐｒｉｏｒ計算、またはＲｅｇｕｌａｒｉｚａｔｉｏｎ計算を行う。以下、両者を代表してＰｒｉｏｒ計算と呼ぶ。逐次修正パラメータは、このＰｒｉｏｒ計算の係数として用いる。なお、逐次修正パラメータは、Ｐｒｉｏｒ計算の係数の代わりにＬｉｋｅｌｉｈｏｏｄ計算の係数に用いてもよい。
【００４４】
本実施形態の逐次再構成部２３３では、代表的な撮影条件（代表撮影条件）において、目標とするノイズ低減割合あるいはＸ線量の低減割合と逐次修正パラメータとの関係を予め算出し、保持する。そして、実際の撮影時は、実際の撮影条件（実撮影条件）で指定された（所望の）ノイズ低減割合あるいはＸ線量の低減割合を、代表撮影条件でのノイズ低減割合あるいはＸ線量の低減割合に換算し、その低減割合に対応する逐次修正パラメータを決定し、逐次修正に用いる。本実施形態の逐次再構成部２３３による逐次修正処理の詳細は、後述する。
【００４５】
画像表示部２３４は、逐次再構成部２３３で逐次修正後のＣＴ画像を、表示する。画像表示部２３４は、モニタ１３５により実現される。
【００４６】
なお、補正処理部２３２および逐次再構成部２３３は、ＣＰＵ１３３が、メモリ１３２、ＨＤＤ装置１３４等に予め格納されている所定のプログラムを展開・起動することにより実現される。
【００４７】
すなわち、ＣＰＵ１３３は、補正処理部２３２として信号に対して補正を行い、逐次再構成部２３３として逐次処理を用いて画像再構成を行う。また、ＨＤＤ装置１３４等にデータは保存され、必要に応じて、データは外部へ入出力される。画像再構成により得たＣＴ画像は、画像表示部２３４として機能する液晶ディスプレイやＣＲＴ等のモニタ１３５に表示される。
【００４８】
なお、ＣＰＵ１３３が実現する全部または一部の機能は、ＡＳＩＣ（Ａｐｐｌｉｃａｔｉｏｎ  Ｓｐｅｃｉｆｉｃ  Ｉｎｔｅｇｒａｔｅｄ  Ｃｉｒｃｕｉｔ）、ＦＰＧＡ（ｆｉｅｌｄ−ｐｒｏｇｒａｍｍａｂｌｅ  ｇａｔｅ  ａｒｒａｙ）などのハードウェアによって実現してもよい。また、各機能の処理に用いる各種のデータ、処理中に生成される各種のデータは、メモリ１３２、ＨＤＤ装置１３４等に格納される。
【００４９】
また、信号収集部２３１は、信号取得部２２２が備えてもよい。この場合、撮影部１２０は、ディジタル信号を出力する。例えば、ネットワークを介して画像生成部１３０を接続する場合、この構成は好適である。
【００５０】
また、入力部１１０と画像生成部１３０とは、独立したハードウェアとしても良いし、これらのハードウェアを共用した構成としても良い。従って、上述のようにＣＰＵ１３３やメモリ１３２やモニタ１３５等は入力部１１０と共用できる。
【００５１】
［撮影条件受付画面］
ここで、撮影条件入力部２１１がモニタに表示し、撮影条件の入力を受け付ける撮影条件受付画面について説明する。本実施形態の撮影条件受付画面４００の例を図３に示す。
【００５２】
本実施形態の撮影条件受付画面４００は、照射するＸ線のエネルギー及び出力量に対応する管電圧、管電流時間積、および、１回転における撮影回数を設定するためのＸ線条件設定領域４１０と、再構成画像の範囲を設定するための再構成範囲設定領域４２０と、逐次再構成に用いる重みを選択するための重み設定領域４３０と、撮影部位を設定するための撮影部位設定領域４４０と、所望の画質を選択するための撮影・画像設定領域４５０と、再構成フィルタを選択するための再構成フィルタ設定領域４６０と、測定投影データの補正方法を選択するための投影データ補正設定領域４７０と、投影データから再構成したＣＴ画像（以下、初期画像と呼ぶ。）または逐次修正に使用する測定投影データの範囲を選択するための投影データ範囲設定領域４８０と、を備える。
【００５３】
ユーザは、モニタに表示される撮影条件受付画面４００を見ながら、マウス１１２やキーボード１１１等を操作して、Ｘ線条件をＸ線条件設定領域４１０に、再構成範囲を再構成範囲設定領域４２０に、重み条件を重み設定領域４３０に、撮影部位を撮影部位設定領域４４０に、所望の画質を撮影・画像設定領域４５０に、再構成フィルタの条件を再構成フィルタ設定領域４６０に、測定投影データの補正方法を投影データ補正設定領域４７０に、初期画像または逐次修正に使用する測定投影データの範囲を投影データ範囲設定領域４８０にそれぞれ設定する。以下、さらに詳しく説明する。
【００５４】
図３では一例として、ユーザによってＸ線条件設定領域４１０に、管電圧値１２０［ｋＶ］、管電流時間積２００［ｍＡｓ］、撮影回数９００［回／回転］が設定されている例を示す。なお、図３では、１種類のエネルギースペクトルを有するＸ線を用いる例について示しているが、２種類以上のＸ線を用いるマルチエネルギーＣＴの場合には、管電圧、管電流時間積、及び撮影回数の項目をＸ線条件設定領域４１０に追加し、Ｘ線の種類ごとに同様に設定する。
【００５５】
再構成範囲設定領域４２０において、ユーザは、画像再構成を行う領域である再構成範囲（Ｆｉｅｌｄ  ｏｆ  Ｖｉｅｗ、以下、ＦＯＶとする）を設定する。ＦＯＶは、ＦＯＶの大きさと中心位置を設定することにより設定される。本実施形態では一例として、ＦＯＶを正方形で定義する。図３の例では、ＦＯＶは一辺７００［ｍｍ］が設定され、ＦＯＶの中心位置は、Ｘ線ＣＴ装置１００の回転中心と等しい、Ｘ＝Ｙ＝Ｚ＝０［ｍｍ］に設定される。ただし、ＦＯＶは、正方形に限ることはなく、円形、長方形、立方体、直方体、球等の任意の形状に設定することも可能である。この場合も本実施形態の構成を適用できる。
【００５６】
重み設定領域４３０は、Ｘ線検出部１２２の各検出素子で検出した測定投影データと順投影計算で算出する計算投影データとの差分に付与する重みの種類を設定する。重みの種類は、全ての検出素子のデータに付与する重みを一定にする「一定値」と、検出素子で検出したＸ線光子数に応じた重みを付与する「検出素子の光子数」の２種類から選択される。なお、Ｘ線光子数とは、例えば公知であるフォトンカウンティングタイプの検出素子で検出したＸ線光子の数や、前述したディジタル信号に変換後の測定投影データから逆対数変換を用いて、検出素子で検出したＸ線光子の数を推定したものを指す。ここでは、「一定値」が選択されている場合を例示する。以下、本明細書では、各検出素子で検出した測定投影データと順投影計算で算出する計算投影データとの差分に付与する重みを、単に、検出素子の重みと呼ぶ。
【００５７】
撮影部位設定領域４４０は、撮影部位の選択を受け付ける。本実施形態では、列挙したＸ線照射対象（頭部、胸部、肺野等の部位や組織）の中から選択を受け付ける、および、Ｘ線照射対象を近似した楕円体の条件として数値の指定を受け付ける、のいずれかにより受け付ける。ここでは、列挙したＸ線照射対象の中から「頭部」が選択されている場合を例示する。
【００５８】
撮影・画像設定領域４５０は、所望の画質の選択を受け付ける。本実施形態では、例えば、列挙したモードの中から選択を受け付けるとともに、低減割合、あるいはノイズ値を数値で受け付ける。
【００５９】
ここでは、モードとして、所望のノイズ低減割合を達成したＣＴ画像を取得するモード（「固定のノイズ低減割合」）、所望のＸ線量の低減割合を達成したＣＴ画像を取得するモード（「固定のＸ線量低減割合」）、および、所望のノイズ値を達成したＣＴ画像を取得するモード（「固定ノイズ値」）が用意され、ユーザは、これらのモードの中から選択するとともに、目標値を指定する。
【００６０】
図３には、所望のノイズ低減割合を達成したＣＴ画像を取得するモードが選択され、その低減割合の目標値は７５％である場合を例示する。これは、後述する逐次修正を用いて、初期画像のノイズを７５％低減させたＣＴ画像を取得することを示す。
【００６１】
なお、「固定のＸ線量低減割合」は、指定したＸ線量低減割合に相当するＸ線量で撮影した場合に、低減前のＸ線量で撮影したデータを解析的に再構成することで得られるＣＴ画像と同等の画像ノイズを取得するモードである。また、「固定のノイズ値」は、逐次修正を用いて所望するノイズ値をもつＣＴ画像を取得するモードである。
【００６２】
本実施形態では、後述するように、逐次修正に入力する初期画像を基準としてノイズの低減割合を計算しているが、初期画像とは別に、逐次修正前にＸ線ＣＴ装置１００から出力するＣＴ画像（以下、装置出力画像）を基準としてもよい。以降、初期画像と装置出力画像は同様に扱うことが可能なため、装置出力画像の記載を省略し、初期画像に統一し、説明する。
【００６３】
再構成フィルタ設定領域４６０は、逐次修正時に用いる初期画像に対し適用する再構成フィルタの種類を受け付ける領域である。再構成フィルタは、ノイズ低減量を決定するための基準となる初期画像の画質を決定する。この再構成フィルタは、部位によって定められていることが多いため、ここでは、部位の指定を受け付けることにより、再構成フィルタの種類を受け付ける。一般に、頭部や肺野等の撮影時に用いる再構成フィルタは、高分解能かつ高ノイズのＣＴ画像を取得する効果がある。一方、腹部の撮影時に用いる再構成フィルタは、低分解能かつ低ノイズのＣＴ画像を取得する効果がある。図３では、頭部用のフィルタが指定された場合を例示する。
【００６４】
投影データ補正設定領域４７０は、初期画像に用いる測定投影データの補正方法の指定を受け付ける領域である。補正効果が大きいほどチャネル方向、投影角度方向、または列方向の検出素子のデータが隣接素子のデータと平滑化処理され、検出素子のデータのノイズを低下できる。ここでは、補正効果をその大きさで段階的に分け、大きさの段階で指定を受け付ける。図３では、補正効果の大きさで補正方法を小、中、大と区分し、大が選択されている場合を例示する。
【００６５】
投影データ範囲設定領域４８０は、初期画像の生成（再構成）および後述する逐次修正それぞれに使用する測定投影データの範囲の指定を受け付ける領域である。ここでは、投影角度等の指定により、受け付ける。例えば、投影角度１８０度分の測定投影データ範囲を使用すれば、高い時間分解能を得ることができる。また、投影角度３６０度分の測定投影データ範囲を使用すると、ノイズを抑えることができる。また、公知であるヘリカルスキャンでは、寝台移動距離［ｍｍ］をＸ線のビーム幅［ｍｍ］にて除算したビームピッチに応じて、ＣＴ画像の各画素位置で収集される測定投影データの範囲が変わる。その他、ＣＴ画像の各画素位置で収集した全てのデータを使用する（最大収集値）といった指定も可能である。図３では、初期画像については、３６０度分の測定投影データ範囲を使用し、逐次修正においては、全てのデータを使用するよう選択がなされている場合を例示する。
【００６６】
なお、撮影条件受付画面４００は、図３の画面構成に限定されない。また、撮影条件受付画面４００で設定を受け付けたＸ線条件、再構成範囲、重みの設定条件、撮影部位の設定条件、撮影・画像条件、再構成フィルタ、投影データの補正方法、および投影データ範囲の組み合わせをＨＤＤ装置１１５に保存してもよい。この場合、次の同条件の撮影時に、ユーザの指示に従って、撮影条件入力部２１１がＨＤＤ装置１１５から読み出し、使用する。この場合、ユーザが、毎撮影時に、Ｘ線条件等を入力する必要はない。また、上記設定条件の組み合わせを予め複数種類、ＨＤＤ装置１１５に保存しておき、ユーザが複数種類の中から選択するよう構成してもよい。
【００６７】
［逐次修正処理］
次に、本実施形態の逐次再構成部２３３による逐次修正処理の詳細を、図４および図５を用いて説明する。
【００６８】
上述のように、本実施形態では、実撮影条件で指定された（所望の）ノイズ低減割合あるいはＸ線量の低減割合を、代表撮影条件でのノイズ低減割合あるいはＸ線量の低減割合に換算し、その低減割合に対応する逐次修正パラメータを、逐次修正に用いる。なお、固定ノイズ値のモードが選択された場合、初期画像における該当範囲のノイズ値を計測し、所望の固定ノイズ値との割合から実撮影条件で指定された（所望の）ノイズ低減割合あるいはＸ線量の低減割合を計算する。そして、低減割合に対応する逐次修正パラメータを、逐次修正に用いる。
【００６９】
撮影条件のうち、逐次修正パラメータに影響を与える撮影条件には、図３の撮影条件受付画面４００を介して設定される撮影条件のうち、再構成フィルタ設定領域４６０を介して受け付ける再構成フィルタ、投影データ補正設定領域４７０を介して受け付ける測定投影データの補正方法、および、投影データ範囲設定領域４８０を介して受け付ける逐次修正に使用する測定投影データの範囲、が含まれる。
【００７０】
以下、本実施形態では、ユーザが設定する撮影条件（実撮影条件）と代表撮影条件とは、再構成フィルタが異なる場合の換算手法を説明する。
【００７１】
逐次再構成部２３３は、前述のように、まず、初期画像を順方向に投影計算（以下、順投影計算とする）することにより計算投影データを求める。求めた計算投影データと測定投影データとが等しくなるように、初期画像に対し逐次修正を行い、最終的にＣＴ画像を得る。逐次修正時に用いる逐次修正パラメータは、画素位置毎に異なる重みに応じた最適な値を用いる。この逐次修正パラメータは、代表的な撮影条件（代表撮影条件）において最適化したパラメータ（基準パラメータ）を、撮影時の撮影条件に最適な値に換算することにより、算出する。
【００７２】
なお、逐次修正パラメータは、所望のノイズ低減割合、または、Ｘ線量低減割合を達成したＣＴ画像を算出するために、逐次修正時に用いるパラメータである。
【００７３】
これを実現するため、本実施形態の逐次再構成部２３３は、図４に示すように、代表撮影条件における、ノイズ低減割合と逐次修正パラメータとの関係をテーブルとして保持するとともに、代表撮影条件の画像における基準となる重みを保持するテーブル部３１０と、実際の撮影時に受け付けた撮影条件に応じて、画素毎の逐次修正に用いる逐次修正パラメータを決定するパラメータ決定部３２０と、決定した逐次修正パラメータを用いて初期画像を逐次修正する逐次修正部３３０と、を備える。
【００７４】
［ノイズの低減割合  Ｘ線量の低減割合］
各部の説明に先立ち、ノイズの低減割合およびＸ線量の低減割合について説明する。
【００７５】
本実施形態において、ノイズの低減割合とは、例えば、以下の式（１）に示すように、公知であるＦｅｌｄｋａｍｐ法等の解析的な再構成法を用いて再構成した初期画像のノイズを基準として、逐次修正後、ＣＴ画像のノイズが低減した割合を百分率で示すものである。
【数１】
2016199716.tif 000003
【００７６】
一方、本実施形態において、Ｘ線量の低減割合とは、前述した解析的な再構成法を用いて再構成したＣＴ画像と同等な画質を取得する条件の下、逐次再構成により低減可能なＸ線量の割合を百分率で示す。本実施形態において、画質は任意領域のノイズを示す標準偏差ＳＤを用いて説明しているが、空間分解能等の別の評価指標を用いても構わない。Ｘ線量の低減割合は、ＳＤからＸ線量を近似的に換算可能であるため、下式（２）で表すことができる。
【数２】
2016199716.tif 000004
【００７７】
以降、ノイズの低減割合とＸ線量の低減割合は、式（１）、および式（２）から換算可能であるため、Ｘ線量の低減割合の記載を省略し、ノイズの低減割合で代表する。
【００７８】
［テーブル部］
テーブル部３１０は、代表的な撮影条件である代表撮影条件で取得したＣＴ画像（初期画像）について、逐次修正によりノイズまたはＸ線量を低減させる割合である低減割合と、当該低減割合を実現する逐次修正に用いるパラメータである基準パラメータとの関係を保存する。低減割合と基準パラメータとの関係は、例えば、ノイズ低減割合を変数とする基準パラメータの関数として予め作成し、メモリ１３２あるいはＨＤＤ装置１３４に保持する。以下、本明細書では、ノイズ低減割合と基準パラメータとの関係を示す関数を、テーブルと呼ぶ。
【００７９】
また、検出素子毎の重みから算出する重み画像は、画素位置毎に重み値が異なる。従って、撮影条件だけでなく、画像上の画素位置に応じても、所望の画質（ノイズ低減割合、空間分解能など）の画像を得るために最適な逐次修正パラメータが異なる。本実施形態のテーブル部３１０は、画素位置毎に最適な逐次パラメータを決定するために、基準となる位置の重み（基準重みと呼ぶ）も算出し、保持する。
【００８０】
これを実現するため、本実施形態のテーブル部３１０は、図４に示すように、テーブル計算部３１１と、基準重み計算部３１２と、計算テーブル保存部３１３と、を備える。各部は、ＣＰＵ１３３がプログラムを実行することにより実現されるが、その一部は、メモリ１３２およびＨＤＤ装置１３４により構成される記憶部上に構築される。
【００８１】
［テーブル計算部］
テーブル計算部３１１は、事前に取得した代表的なファントムを用いて、基準パラメータの大きさとＣＴ画像のノイズ低減割合との関係を特定する関数（テーブル）を作成する（図５のステップＳ１１０１）。テーブルは、実際に当該ファントムのＣＴ画像（初期画像）を逐次修正した結果を用いて算出する。なお、このＣＴ画像は、代表撮影条件で取得する。
【００８２】
テーブル計算部３１１は、得られた初期画像を、複数の異なる逐次修正パラメータである基準パラメータβｂを用いて逐次修正する。そして、初期画像のノイズと、逐次修正後のＣＴ画像のノイズとをそれぞれ計測し、基準パラメータβｂ毎のノイズの低減割合を算出する。ノイズの計測は、予め設定したＲＯＩ内で行う。
【００８３】
ノイズ、および重みを計測するためのＲＯＩは、ＣＴ画像の、Ｘ線ＣＴ装置１００の回転中心の位置に設定する。しかしながら、ＲＯＩの位置は、これに限定されない。周辺の位置に設定しても構わない。また、中心以外の周辺に位置する複数のＲＯＩを設定しても構わない。このとき、ＲＯＩ毎にノイズを計測し、テーブルを作成してもよいし、複数ＲＯＩのノイズの計測値を平均化し、テーブルを作成してもよい。
【００８４】
テーブル計算部３１１の処理を、図６（ａ）を用いて説明する。図６（ａ）は、内部を水で満たした、直径３０ｃｍ、高さ１００ｃｍの円柱ファントムを、代表撮影条件で撮影し、得た断層面のＣＴ画像５１０である。ノイズを計測するＲＯＩ５１１は、上述のように、ＣＴ画像５１０の回転中心に設定される。
【００８５】
テーブル計算部３１１は、まず、初期画像のＲＯＩ５１１内のノイズを計測する。その後、基準パラメータβｂの値を変えて、それぞれ、逐次修正処理を行い、修正後のＣＴ画像のＲＯＩ５１１のノイズを計測する。そして、それぞれ、式（１）に従って、ノイズの低減割合を算出する。
【００８６】
このときの、基準パラメータβｂの値毎のノイズ低減割合の変化の一例を図６（ｂ）に示す。本図において、５１２は、計測値であり、各基準パラメータβｂの値毎の、ノイズ低減割合をプロットしたものである。
【００８７】
テーブル計算部３１１は、計測値５１２のプロット結果から、例えば、最小二乗法等を用いて近似曲線５１３を算出する。そして、この基準パラメータβｂの値と、ノイズ低減割合との関係を示す近似曲線５１３を表す関数を、テーブルとして保持する。このテーブルを用いることにより、任意のノイズ低減割合を実現する基準パラメータβｂを得ることができる。
【００８８】
なお、逐次修正は、上述のように、第一の計算と第二の計算とを繰り返す。その繰返し回数は、予め定められる。そして、同じ逐次修正パラメータであっても、繰返し回数によって、ノイズの低減割合は異なる。従って、テーブル計算部３１１は、例えば、繰返し回数毎に、上記テーブルを作成してもよい。
【００８９】
上述の手順で、テーブル計算部３１１が作成したテーブルの一例を図７に示す。ここでは、繰返し回数（更新回数）毎に作成したものを例にあげて説明する。繰返し回数毎に、テーブルを保持するデータベースを、パラメータテーブルと呼ぶ。本図に示すように、パラメータテーブル７１０は、繰返し回数（更新回数）７１１毎に、テーブル７１３を保持する。このとき、代表撮影条件として用いた撮影条件（本実施形態では、再構成フィルタの種類）７１２も併せて保持してもよい。
【００９０】
なお、ここでは、実際にＸ線ＣＴ装置１００を用いて取得したＣＴ画像を用いて、パラメータテーブル７１０を作成する場合を例にあげて説明したが、パラメータテーブル７１０の作成手法は、これに限定されない。パラメータテーブル７１０は、仮想的なＸ線ＣＴ装置によるシミュレーションデータを用いて作成してもよい。
【００９１】
［基準重み計算部］
基準重み計算部３１２は、検出素子毎の重みに基づいて、ＣＴ画像上の予め定めた基準とする領域（基準領域）の画素の重み値を基準重みとして算出する。具体的には、基準重み計算部３１２は、ユーザが撮影条件として指定可能なＦＯＶ、およびＦＯＶの中心位置の条件に応じて、各検出素子の重みＷ（ｉ）から重み分布（重み画像）ＷＩ（ｊ）を計算し、基準領域の重みＷｂ（以下、基準重みとする）を取得する（図５のステップＳ１１０２）。
【００９２】
なお、検出素子毎の重みは、予め、撮影条件受付画面４００の重み設定領域４３０を介して設定される。また、重み画像ＷＩ（ｊ）は、検出素子の重みＷ（ｉ）に対し、公知である再構成フィルタ処理を含まない単純逆投影処理を施すことにより得る、画素毎の重み値の分布である。
【００９３】
例えば、ユーザが指定できるＦＯＶが、ＦＯＶ＝３００、５００、７００［ｍｍ］の３種類であり、かつ、ＦＯＶの中心位置が、Ｘ＝Ｙ＝Ｚ＝０［ｍｍ］と、Ｘ＝３０［ｍｍ］、Ｙ＝Ｚ＝０［ｍｍ］と、の２種類の場合、合計６通りの基準重みＷｂを取得する。
【００９４】
基準領域は、テーブル計算部３１１でノイズ計測に用いたＲＯＩ５１１と同じ位置が望ましい。中心以外の周辺に位置する複数のＲＯＩを用いても構わない。
【００９５】
図８（ａ）および図８（ｂ）を用いて、基準重み計算部３１２の処理を具体例で説明する。
【００９６】
図８（ａ）に、検出素子の重みＷ（ｉ）分布を示す。なお、ｉは、投影角度毎の検出素子それぞれに振られた検出素子番号である。投影角度が異なる場合は、同一チャネルおよび列の検出素子においても、別の検出素子番号とする。ここでは、上述の撮影条件受付画面４００の重み設定領域４３０で、一定値が選択された場合の重みＷ（ｉ）分布を示す。すなわち、各重みＷ（ｉ）の値は、全て同じ値である。重みＷ（ｉ）の値は、例えば、１とする。もちろん、１以外の値であってもよい。
【００９７】
図８（ｂ）は、図８（ａ）に示す検出素子の重みＷ（ｉ）分布から作成した、重み画像ＷＩ（ｊ）である。ここでは、ＦＯＶが７００［ｍｍ］、ＦＯＶの中心位置がＸ＝Ｙ＝Ｚ＝０［ｍｍ］と指定された際に取得したＣＴ画像の重み画像ＷＩ（ｊ）を示す。重み画像ＷＩ（ｊ）は、各検出素子の重みＷ（ｉ）に対し、ＣＴの画像再構成と同じ処理を施すことにより得る。なお、ｊは、画素番号である。
【００９８】
図８（ａ）に示すように、各検出素子の重みＷ（ｉ）の値が一定値であるため、重み画像ＷＩ（ｊ）は、画像中心付近に、円形状にほぼ一定の画素値を有し、画像周辺が視野外のため、０に近い値を有する画像となる。
【００９９】
基準重み計算部３１２は、基準領域（ＲＯＩ５１４）の画素ｊの重み値ＷＩ（ｊ）を、以下の式（３）を用いて算出する。画素ｊの重みＷＩ（ｊ）は、図８（ｂ）のＲＯＩ５１４内の画素の重み画像を計算した結果である。
【数３】
2016199716.tif 000005
ここで、Ｉは投影角度毎において異なる検出素子番号の総数、Ｌは、総画素数、ｌは、更新対象の画素ｊと検出素子ｉを結ぶライン上にＬ個の画素がある場合、それらの画素の番号を表す。Ｃ（ｉ，ｊ）は、画素ｊが検出素子ｉに寄与する割合を表し、Ｘ線検出部１２２の位置や順投影計算、または、逆投影計算の手法により異なる。
【０１００】
そして、基準重み計算部３１２は、ＲＯＩ５１４内の各画素ｊの重みＷＩ（ｊ）から、基準重みＷｂを算出する。基準重みＷｂは、例えば、各画素ｊの重みＷＩ（ｊ）の平均を算出することにより得る。なお、平均値に限らず、ＷＩ（ｊ）の、最大値、最小値などを基準重みＷｂとしてもよい。
【０１０１】
［計算テーブル保存部］
計算テーブル保存部３１３は、テーブル計算部３１１が計算したテーブル７１３、および、基準重み計算部３１２が算出した基準重みＷｂを、メモリ１３２やＨＤＤ装置１３４等に保存する。保存されたテーブル７１３および基準重みＷｂは、パラメータ決定部３２０が参照する。
【０１０２】
なお、テーブル部３１０の処理は、被写体１０１の撮影前に行う。なお、Ｘ線ＣＴ装置１００の出荷より前であってもよい。
【０１０３】
また、本実施形態では、図７に示すように、実際に逐次修正で用いられる可能性のある各繰返し回数（ここでは、２０回、６０回、・・・）７１１についてテーブル７１３を取得しているが、これに限定されない。２以上の異なる繰返し回数（更新回数）７１１のテーブル７１３を取得しておけばよい。この場合、他の繰返し回数のテーブルは、取得した繰返し回数のテーブル７１３に基づいて、公知の線形補間等で内挿または外挿により得る。これにより、保存するテーブル７１３の数を低減できる。
【０１０４】
本実施形態では、図７に示すように、代表撮影条件の再構成フィルタ７１２を用いて、繰返し回数（更新回数）７１１毎に、基準パラメータとノイズ低減割合とを対応づける関数をテーブル７１３として決定する。しかしながら、代表撮影条件は、これに限定されない。各管電圧、寝台の動作速度等、他の撮影条件または再構成条件に対して、それぞれテーブル（基準パラメータとノイズ低減割合とを対応づける近似曲線）を取得しても良い。なお、更新回数とは、逐次修正の繰返し回数のことである。
【０１０５】
本実施形態では、管電圧等の撮影条件、ＦＯＶ等の再構成条件、および測定投影データの種類に応じて、複数の計算テーブルを保存してもよい。これにより、条件の違いによる逐次修正パラメータの真値との誤差を低減することができる。
【０１０６】
［パラメータ決定部］
パラメータ決定部３２０は、実際の撮影時の撮影条件である実撮影条件における所望の低減割合を、代表撮影条件における低減割合である代表低減割合に変換し、テーブル部３１０が保存する関係（テーブル）を参照して、当該所望の低減割合を実現する逐次修正パラメータを決定する。実撮影条件は、撮影条件受付画面４００で受け付ける。
【０１０７】
すなわち、パラメータ決定部３２０は、実撮影条件として、撮影条件受付画面４００で受け付けた撮影条件で、当該撮影条件受付画面４００の、撮影・画像設定領域４５０を介して受け付けたノイズ低減割合（実低減割合）を実現する逐次修正パラメータを決定する。このとき、パラメータ決定部３２０は、実低減割合を、代表撮影条件のノイズ低減割合（代表低減割合）に換算し、上記テーブル７１３に従って、逐次修正パラメータを得る。
【０１０８】
これを実現するため、本実施形態のパラメータ決定部３２０は、図４に示すように、撮影条件として設定された所望の低減割合（実ノイズ低減割合）を、代表ノイズ低減割合に変換する代表ノイズ低減割合変換部３２１と、テーブル部３１０が保存した関係（テーブル）を参照し、代表ノイズ低減割合に対応する基準パラメータβｂを決定する計算テーブル参照部３２２と、各検出素子の重みに基づいてＣＴ画像の各位置（画素）の重みを計算する重み計算部３２３と、基準パラメータβｂと各位置の重みとを用いて各位置のパラメータに変換するパラメータ変換部３２４と、を備える。
【０１０９】
以下、本実施形態のパラメータ決定部３２０の各部の処理の詳細を、図５を用いて説明する。
【０１１０】
［代表ノイズ低減割合変換部］
代表ノイズ低減割合変換部３２１は、代表撮影条件と実撮影条件との差異による、実撮影条件で取得したＣＴ画像のノイズに対する代表撮影条件で取得したＣＴ画像のノイズの比であるノイズ増減比を用いて、撮像条件として設定された所望の低減割合（実ノイズ低減割合）を、代表撮影条件のテーブル７１３に入力するためのノイズ低減割合（代表低減割合）に変換する（図５のステップＳ１２０１、Ｓ１２０２）。以下、変換処理の詳細を説明する。
【０１１１】
変換処理の説明に先立ち、繰返し回数（更新回数）とノイズ（ＳＤ）との関係を説明する。この関係を図９（ａ）に示す。図９（ａ）は、横軸を繰返し回数（更新回数）、縦軸をノイズ（ＳＤ）とし、繰返し回数（更新回数）の増加に伴うノイズ（ＳＤ）５３１の変化の様子を示すグラフである。
【０１１２】
本図に示すように、ノイズ（ＳＤ）は、繰返し回数（更新回数）が増えるにつれて低減し、最終的には、所定の値５３４に収束する。図９（ａ）には、一例として、繰返し回数（更新回数）が６０回を超えると、ノイズＳＤ５３１が変化しない収束段階となる場合を例示する。なお、繰返し回数（更新回数）０回は初期画像を意味する。
【０１１３】
本実施形態では、同一の逐次修正パラメータを用いる時、初期画像のノイズによらず、逐次修正の収束段階のノイズは等しいものとする。
【０１１４】
まず、代表ノイズ低減割合変換部３２１は、実撮影条件の再構成フィルタで取得したＣＴ画像（初期画像）のノイズ（図９（ａ）の５３２）を用い、代表撮影条件で取得するＣＴ画像（初期画像）のノイズ（図９（ａ）の５３３）を算出する（ステップＳ１２０１）。
【０１１５】
ここで、代表ノイズ低減割合変換部３２１は、まず、実撮影条件の再構成フィルタで取得した初期画像のノイズ（ＳＤＡｉｎｔ）５３２を算出する。算出は、実際に得た初期画像の画素値を用いて行う。そして、撮影・画像設定領域４５０で設定された実低減割合（ＲＲＡ％）を用い、逐次修正後の目標ノイズ（ＳＤＡｔａｒ）を算出する。目標ノイズ（ＳＤＡｔａｒ）は、以下の式（４）によって算出できる。
ＳＤＡｔａｒ＝ＳＤＡｉｎｔ×（１００−ＲＲＡ）/１００  ・・・（４）
【０１１６】
次に、代表ノイズ低減割合変換部３２１は、以下の式（５）に示す、代表撮影条件で取得したＣＴ画像のノイズの、実際の撮影条件で取得したＣＴ画像のノイズに対する比であるノイズ増減比を算出する。
【数５】
2016199716.tif 000006
【０１１７】
本実施形態では、実際の撮影条件と代表撮影条件とは、再構成フィルタのみが異なるものとしている。従って、両画像のノイズの差は、再構成フィルタによる。本実施形態では、図９（ｂ）に示す、再構成フィルタの種類５２１ごとに、ノイズ（ＳＤ）５２２を保持するフィルタノイズテーブル５２０を用意し、フィルタ間ノイズ比を算出し、ノイズ増減比とする。
【０１１８】
フィルタノイズテーブル５２０は、使用を想定される再構成フィルタ毎に、同一の対象を前記代表撮影条件で撮影した際の前記ＣＴ画像のノイズを保存する。ここでは、各再構成フィルタについて、代表的なファントムを用いて、代表撮影条件で撮影し、ＳＤを計測し、その結果を保存する。図９（ｂ）には、使用を想定される再構成フィルタの種類５２１として、腹部用、頭部用、胸部肺野用、Ｒａｍｐ、Ｓｈｅｐｐ−ｌｏｇａｎが保持される場合を例示する。本図に示すように、なお、フィルタノイズテーブル５２０は、予め作成され、ＨＤＤ装置１３４等に保持される。
【０１１９】
代表ノイズ低減割合変換部３２１は、前記フィルタノイズテーブル５２０を参照し、所望の低減割合（実低減割合）を、代表低減割合に換算する。すなわち、ノイズ増減比（ＲＳＤ）として、代表撮影条件の再構成フィルタによるノイズ（ＳＤＲｆｉｌ）の、実撮影条件の再構成フィルタによるノイズ（ＳＤＡｆｉｌ）に対するフィルタ間ノイズ比（ＳＤＲｆｉｌ/ＳＤＡｆｉｌ）を算出する。
【０１２０】
例えば、代表撮影条件の再構成フィルタの種類５２１がＲａｍｐで、実撮影条件の再構成フィルタの種類５２１が頭部用の場合、フィルタノイズテーブル５２０の値を用いて算出されるフィルタ間ノイズ比は、２０／４０＝０．５となる。
【０１２１】
そして、代表ノイズ低減割合変換部３２１は、このノイズ増減比（ＲＳＤ）を用い、式（５）の関係から、代表撮影条件で取得した初期画像のノイズ（ＳＤＲｉｎｔ）を算出する。代表撮影条件で取得した初期画像のノイズ（ＳＤＲｉｎｔ）は、以下の式（６）で算出できる。
ＳＤＲｉｎｔ＝ＲＳＤ×ＳＤＡｉｎｔ  ・・・（６）
【０１２２】
次に、代表ノイズ低減割合変換部３２１は、代表撮影条件で取得した初期画像のノイズ（ＳＤＲｉｎｔ）と、目標ノイズ（ＳＤＡｔａｒ）とを用い、代表撮影条件におけるノイズ低減割合（代表ノイズ低減割合ＲＲＲ％）を算出する（ステップＳ１２０２）。
【０１２３】
代表ノイズ低減割合（ＲＲＲ）は、以下の式（７）で算出できる。
ＲＲＲ＝（１−ＳＤＡｔａｒ/ＳＤＲｉｎｔ）×１００  ・・・（７）
【０１２４】
なお、上記では、代表ノイズ低減割合変換部３２１は、実撮影条件で算出した初期画像のノイズ（ＳＤＡｉｎｔ）から初期画像を代表撮影条件で取得した場合のノイズ（ＳＤＲｉｎｔ）を算出し、目標ノイズ（ＳＤＡｔａｒ）とから、代表低減割合（ＲＲＲ）を算出しているが、この手法に限定されない。
【０１２５】
例えば、代表ノイズ低減割合ＲＲＲは、実撮影条件として撮影・画像設定領域４５０を介して設定されたノイズ低減割合（実ノイズ低減割合）と、ノイズ増減比とを用いて、以下の式（８）により算出できる。
【数８】
2016199716.tif 000007
【０１２６】
上述のように、ノイズ増減比は、フィルタノイズテーブル５２０から算出できる。従って必ずしも、上述のように、実撮影条件の初期画像のノイズ（ＳＤＡｉｎｔ）や、代表撮影条件の初期画像のノイズ（ＳＤＲｉｎｔ）を算出しなくてもよく、式（８）により、直接代表ノイズ低減割合を算出してもよい。
【０１２７】
例えば、実ノイズ低減割合が７５％、代表撮影条件の再構成フィルタの種類５２１がＲａｍｐで、実撮影条件の再構成フィルタの種類５２１が頭部用の場合、代表ノイズ低減割合は、（１−（１−７５／１００）／０．５）×１００＝５０％となる。
【０１２８】
［計算テーブル参照部］
計算テーブル参照部３２２は、実ノイズ低減割合を変換して得た代表ノイズ低減割合に対応する基準パラメータβｂをテーブル部３１０が保存した関係（テーブル）から抽出し、逐次修正パラメータを決定する（図５のステップＳ１２０３）。
【０１２９】
具体的には、計算テーブル参照部３２２は、テーブル部３１０のテーブル計算部３１１が作成したテーブル７１３に従って、代表ノイズ低減割合に対応する基準パラメータβｂを決定する。
【０１３０】
例えば、図７のパラメータテーブル７１０内で、更新回数７１１が６０回に対応づけられるテーブル７１３の近似曲線が、図６（ｂ）の変化態様を示すものとする。計算テーブル参照部３２２は、例えば、代表ノイズ低減割合が５０％と算出された場合、この近似曲線５１３に従って、基準パラメータβｂとして、０．６５を得る。
【０１３１】
［重み計算部］
重み計算部３２３は、検出素子毎の重みＷ（ｉ）に基づいて、ＣＴ画像上の画素毎の重みＷＩ（ｊ）を計算する。（図５のステップＳ１２０４）。重み画像ＷＩ（ｊ）の算出手法は、基準重み計算部３１２による算出手法と同様である。ただし、基準重み計算部３１２は、基準領域部分のみの重み値を算出するが、重み計算部３２３は、画像内の全領域の重み値を算出する。
【０１３２】
従って、例えば、基準重み計算部３１２が、基準重みを算出する際、全領域の重み値を算出し、それをメモリ１３２等に保持しておくよう構成し、重み計算部３２３は、メモリ１３２等から読み出すよう構成してもよい。
【０１３３】
算出した重み画像ＷＩ（ｊ）を図１０（ａ）に示す。重み画像ＷＩ（ｊ）は、ＣＴ画像の各位置（画素）の重み値を画素値とする画像である。ｊは、画素番号を示す。上述のように、本実施形態では、各検出素子の重みＷ（ｉ）の値は一定とする。重み画像ＷＩ（ｊ）は、この重みＷ（ｉ）分布から作成するため、各画素位置（例えば、ｊ１、ｊ２）において、画素値（ＷＩ（ｊ１）、ＷＩ（ｊ２））は、それぞれ異なる。
【０１３４】
［パラメータ変換部］
パラメータ変換部３２４は、画素ｊ毎の重みＷＩ（ｊ）と基準パラメータβｂとから、画素毎の逐次修正パラメータβＩ（ｊ）を決定する。
【０１３５】
パラメータ変換部３２４は、基準重み計算部３１２が算出した基準重みＷｂと、計算テーブル参照部３２２が特定した基準パラメータβｂと、重み計算部３２３が計算した重み画像ＷＩ（ｊ）とを用いて、下式（９）に従って、各画素位置ｊの逐次修正パラメータの分布βＩ（ｊ）を計算する。
【数９】
2016199716.tif 000008
なお、上述のように、基準重みＷｂは、撮影条件で指定されたＦＯＶ、およびＦＯＶの中心位置と同じ条件で算出されたものである。
【０１３６】
パラメータ変換部３２４が計算した、各画素位置ｊの、逐次修正パラメータの分布βＩ（ｊ）を、図１０（ｂ）示す。各画素位置（例えば、ｊ１、ｊ２）において、画素値（ＷＩ（ｊ１）、ＷＩ（ｊ２））が異なるため、パラメータβの値（βＩ（ｊ１）、βＩ（ｊ２））も、それぞれ異なる。
【０１３７】
［逐次修正部］
逐次修正部３３０は、パラメータ決定部３２０で計算した逐次修正パラメータβＩ（ｊ）を用いて逐次修正（逐次再構成）処理を行う（図５のステップＳ１３００）。本実施形態では、逐次修正部３３０は、画素ｊ毎に、画素ｊ毎の逐次修正パラメータβＩ（ｊ）を適用し、逐次修正を行う。
【０１３８】
逐次修正部３３０は、ＣＴ画像を計算により順投影して求めた計算投影データと、測定投影データとが近づくように逐次的にＣＴ画像を修正する。このとき、逐次修正パラメータとして、パラメータ決定部３２０で計算した、各画素位置ｊの逐次修正パラメータβＩ（ｊ）（逐次修正パラメータの分布βＩ（ｊ））を用いてＣＴ画像を逐次修正する。これにより、高精度にノイズを除去したＣＴ画像を生成する。
【０１３９】
画像を修正するアルゴリズムには、公知の逐次再構成手法を用いることができる。ここでは一例として、ＳＰＳ（Ｓｅｐａｒａｂｌｅ−Ｐａｒａｂｏｌｏｉｄａｌ−Ｓｕｒｒｏｇａｔｅ）法を用いる場合について説明する。
【０１４０】
これを実現するため、逐次修正部３３０は、解析的再構成部３３１と、順投影部３３２と、差分部３３３と、逆投影処理部３３４と、Ｐｒｉｏｒ計算部３３５と、画像修正部３３６と、を備える。以下、逐次修正処理ステップＳ１３００における、各部の処理を、図１１に沿って説明する。
【０１４１】
解析的再構成部３３１は、公知のＦｅｌｄｋａｍｐ法等の解析的再構成手法を用いて、補正処理部２３２が補正した測定投影データＲ（ｉ）からＣＴ画像λｋ＝０（ｊ）を得る（ステップＳ１３０１）。なお、ｋは、逐次再構成の繰返し回数（更新回数、修正回数）を示す０以上の整数であり、ｋ＝０は、初期画像を示す。また、ｊは画素番号で、λｋ（ｊ）は、更新回数がｋ回の画像の画素ｊの画素値を示す。
【０１４２】
ＳＰＳ法では、修正回数が（ｋ＋１）回のＣＴ画像λｋ＋１（ｊ）は、ｋ回のＣＴ画像λｋ（ｊ）を用いて、以下の式（１０）で表される。
【数１０】
2016199716.tif 000009
なお、Ｊは、総画素数である。ｉは、検出素子の番号を示し、Ｉは、総検出素子数である。また、Ｗ（ｉ）は、各検出素子の重みである。また、Ｐ１およびＰ２は、それぞれ分子および分母のＰｒｉｏｒの計算式を表す。
【０１４３】
なお、この逐次再構成法は、一般的な２次元（ｘ、ｙ方向）の断層像だけでなく、１次元データ（ｘ方向）、体軸方向ｚに像を重ねあわせた３次元データ（ｘ、ｙ、ｚ方向）、また、さらに時間方向ｔを考慮した４次元データ（ｘ、ｙ、ｚ、ｔ）にも適用可能である。画素番号ｊは、逐次再構成処理の計算対象の画素を表しており、例えば２次元データの場合、画素番号ｊの範囲は、（１×１）、（２×１）、…、（Ｘ×１）、（１×２）、…、（Ｘ×Ｙ）である。Ｘ、Ｙはそれぞれｘ、ｙ方向の画素数を表す。
【０１４４】
以下、式（１０）による逐次再構成処理、すなわち、修正回数がｋ回のＣＴ画像λｋ（ｊ）を修正し、修正回数が（ｋ＋１）回のＣＴ画像λｋ＋１（ｊ）を算出する処理を、ステップ毎に説明する。
【０１４５】
逐次修正部３３０は、繰返し回数（更新回数）が予め定めた回数になるまで、以下の処理を（ステップＳ１３０２、Ｓ１３０７）。
【０１４６】
順投影部３３２は、ＣＴ画像λｋ（ｊ）の画素を順投影処理し、計算投影データＳ（ｉ）を求める（ステップＳ１３０３）。計算投影データＳ（ｉ）は、以下の式（１１）を計算することにより、得られる。
【数１１】
2016199716.tif 000010
【０１４７】
式（１１）において、ｌは、修正対象の画素ｊとｉ番目の検出素子（検出素子ｉ）とを結ぶライン上にあるＬ個の画素の番号を表す。Ｃ（ｉ、ｌ）は、画素ｌが検出素子ｉに寄与する割合を表す。なお、Ｃ（ｉ、ｌ）の値は、検出素子の位置や順投影計算、または逆投影計算の手法によって異なる値が設定される。
【０１４８】
差分部３３３は、以下の式（１２）に従って、測定投影データＲ（ｉ）から、計算投影データＳ（ｉ）を減算し、修正投影データΔＲｋ（ｉ）を求める（ステップＳ１３０４）。
【数１２】
2016199716.tif 000011
【０１４９】
Ｐｒｉｏｒ計算部３３５は、以下の式（１３）および式（１４）に従って、Ｐ１およびＰ２をそれぞれ計算する。なお、ここでは、Ψは、公知であるＧｅｎｅｒａｌｉｚｅｄ−Ｇｅｍａｎ−Ｐｒｉｏｒの１階導関数を表す。Ｐ１、Ｐ２は、それぞれ公知であるＳｕｒｒｏｇａｔｅ  ｆｕｎｃｔｉｏｎを用いた式である。
【数１３】
2016199716.tif 000012
【０１５０】
このときＰｒｉｏｒ計算部３３５は、ＣＴ画像を構成する画素間のＣＴ値の差分値に対して、逐次修正パラメータを乗算後、Ｌｉｋｅｌｉｈｏｏｄ計算後の修正画像に加算する。この処理は、逐次修正中に画素間のＣＴ値差を小さくすることができ、ノイズを低下させる効果をもつ。
【０１５１】
上記式（１３）、式（１４）において、βは、Ｐｒｉｏｒの強度を示す固定の逐次修正パラメータ（逐次修正パラメータ）である。本実施形態では、βの代わりに、パラメータ決定部３２０で計算した画素ｊ毎の逐次修正パラメータβＩ（ｊ）を用いる。また、（λｊｋ−λｍｋ）は、ＣＴ画像λｋ（ｊ）における、２つの画素のＣＴ値の差分値（λｊｋ−λｍｋ）を変数とする関数である。
【０１５２】
逆投影処理部３３４は、以下の式（１５）により、修正投影データΔＲｋ（ｉ）を逆投影処理し、修正画像Δλｋ（ｊ）を生成する（ステップＳ１３０５）。
【数１５】
2016199716.tif 000013
【０１５３】
画像修正部３３６は、以下の式（１６）を計算することにより、修正画像Δλｋ（ｊ）を用いて修正したＣＴ画像λｋ＋１（ｊ）を求め、繰返し回数（更新回数）（ｋ＋１）番目の画像を得る（ステップＳ１３０６）。
【数１６】
2016199716.tif 000014
【０１５４】
逐次修正部３３０は、繰返し回数（更新回数）（ｋ＋１）番目の画像を得ると、修正回数ｋを１インクリメントし（ステップＳ１３０７）、ステップＳ１３０２へ移行し、インクリメント後の修正回数ｋが予め設定された修正回数Ｋと等しくなるまで（ステップＳ１３０２）、解析的再構成部３３１の処理（ステップＳ１３０３）から処理を繰り返す。これにより、インクリメント後の修正回数ｋが予め設定された修正回数Ｋと等しくなるまで、上記各部による処理が繰り返し行われる。
【０１５５】
修正回数ｋがＫに達すると、画像表示部２３４が得られたＣＴ画像ＣＴ画像λｋ（ｊ）をモニタ１３５に表示する。
【０１５６】
以上説明したように、本実施形態のＸ線ＣＴ装置１００は、設定した撮影の条件に応じてＸ線を発生するＸ線発生部１２１と、被写体１０１を透過後の前記Ｘ線を複数の検出素子で検出し、測定投影データを得るＸ線検出部１２２と、前記測定投影データから生成したＣＴ画像から順投影計算で求めた計算投影データと前記測定投影データとが等しくなるように、前記ＣＴ画像に対し逐次修正を行う逐次再構成部２３３を有する画像生成部１３０と、を備え、前記逐次再構成部２３３は、代表的な撮影条件である代表撮影条件で取得した前記ＣＴ画像について、前記逐次修正によりノイズまたはＸ線量を低減させる割合である低減割合と、当該低減割合を実現する前記逐次修正に用いるパラメータである基準パラメータとの関係をテーブルとして保存するテーブル部３１０と、実際の撮影時の撮影条件である実撮影条件における所望の前記低減割合を、前記代表撮影条件における低減割合である代表低減割合に変換し、前記テーブル部３１０が保存する前記関係を参照して当該所望の低減割合を実現する前記逐次修正パラメータを決定するパラメータ決定部３２０と、前記パラメータ決定部３２０で決定した前記逐次修正パラメータを用いて前記逐次修正を行う逐次修正部３３０と、を備える。
【０１５７】
また、前記パラメータ決定部３２０は、前記代表撮影条件と前記実撮影条件との差異による、前記実撮影条件で取得したＣＴ画像のノイズに対する前記代表撮影条件で取得したＣＴ画像のノイズの比であるノイズ増減比を用いて、前記所望の低減割合を前記代表低減割合に変換する代表ノイズ低減割合変換部３２１と、前記変換後の代表低減割合に対応する基準パラメータを前記テーブルから抽出する計算テーブル参照部３２２と、前記検出素子毎の重みに基づいて、ＣＴ画像上の画素毎の重みを計算する重み計算部３２３と、前記画素毎の重みと前記基準パラメータとから、画素毎の逐次修正パラメータを決定するパラメータ変換部３２４と、を備えてもよい。
【０１５８】
さらに、前記テーブル部３１０は、前記検出素子毎の重みに基づいて、ＣＴ画像上の予め定めた基準とする領域の画素の重み値を基準重みとして算出する基準重み計算部３１２をさらに備え、前記パラメータ変換部３２４は、前記基準重みをさらに用い、前記画素毎の逐次修正パラメータを決定してもよい。
【０１５９】
ユーザが設定する撮影条件と前記代表撮影条件とは、再構成フィルタが異なり、前記逐次再構成部２３３は、前記再構成フィルタ毎に、同一の対象を前記代表撮影条件で撮影した際の前記ＣＴ画像のノイズを保存するフィルタノイズテーブル５２０をさらに備え、前記代表ノイズ低減割合変換部３２１は、前記フィルタノイズテーブル５２０を参照し、前記所望の低減割合を前記代表低減割合に変換してもよい。
【０１６０】
このように、本実施形態によれば、実際の撮影条件の低減割合を、代表撮影条件の低減割合に換算し、実撮影条件において最適な逐次修正パラメータを得る。従って、撮影条件毎に、最適な逐次修正パラメータを保持しておく必要がなく、代表撮影条件についてのみ、最適な逐次修正パラメータを保持しておけばよい。従って、本実施形態によれば、テーブルの作成に要する工数、およびテーブル数を増加させることなく、実撮影条件での最適な逐次修正パラメータを得ることができる。
【０１６１】
また、本実施形態では、事前に算出した基準重みＷｂを用い、各位置の重みに応じた逐次修正パラメータβＩ（ｊ）を決定する。これにより、逐次修正部３３０では、各位置におけるＰｒｉｏｒ計算後のパラメータを決定できるため、ＣＴ画像の領域によらず、所望するノイズの低減割合を達成したＣＴ画像を取得することができる。すなわち、逐次修正部３３０では、画素毎に最適な逐次修正パラメータβＩ（ｊ）を用いてＣＴ画像を逐次修正できる。
【０１６２】
この結果、本実施形態の画像生成部１３０では、高精度にノイズを除去したＣＴ画像を生成できる。このＣＴ画像は、逐次再構成により生成されているため、これを投影した計算投影データは、測定投影データとよく一致し、測定投影データを高精度に画像化したＣＴ画像を得ることができる。
【０１６３】
また、本実施形態では、ノイズの低減効果は、測定投影データと逐次修正パラメータの設定値に依存することに着目し、事前に計算したテーブルからＣＴ画像の各位置のパラメータを算出して用いる。これにより、逐次修正処理自体を大きく変更することなく、所望するノイズ低減割合、撮影条件、再構成条件、および測定投影データに応じて決定した最適な逐次修正パラメータβＩ（ｊ）を容易に導入できる。
【０１６４】
以上説明したように、本実施形態のＸ線ＣＴ装置１００は、逐次再構成（逐次修正）にあたり、画素毎に最適な逐次修正パラメータを用いる。このため、画素毎に異なる重みに応じた最適な逐次修正パラメータを用いて再構成を行うことができる。従って、高い精度で画質の制御を行うことができ、所望するノイズの低減割合を達成した高品質なＣＴ画像を得ることができる。
【０１６５】
また、本実施形態では、再構成フィルタに関し、代表的な撮影条件にのみ対応づけて、各画素に最適な逐次修正パラメータを保持する。そして、実際の撮影時は、実際の撮影条件用の値に換算して用いる。従って、メモリ容量も節約できる。
【０１６６】
なお、本実施形態では、一周３６０度分の測定投影データを用いて、ＣＴ画像を再構成する場合を例にあげて説明した。しかしながら、これに限定されない。ＣＴ再構成手法は、例えば、公知であるハーフ再構成または一周３６０度分以上の測定投影データを用いる再構成であってもよい。
【０１６７】
＜実施例＞
本実施形態の有効性を検証するため、内部が水で満たされた直径３８ｃｍ、高さ１００ｃｍの円柱ファントムを撮影した。ファントムは、Ｘ線ＣＴ装置１００の回転中心に設置した。
【０１６８】
図１２（ａ）は、初期画像６１１を示し、図１２（ｂ）は、従来方式で逐次修正して得たＣＴ画像６１２を示し、図１２（ｃ）は、本実施形態の方式で逐次修正して得たＣＴ画像６１３を示す。なお、従来方式は、画素位置ｊによらず同じ逐次修正パラメータβを用いた逐次修正である。また、本実施形態の方式は、画素位置ｊ毎に最適な逐次修正パラメータβＩ（ｊ）を用いた逐次修正である。図１２（ａ）〜図１２（ｃ）の各画像の右下の白枠内は、各再構成画像の中心部を拡大した結果である。また、ウィンドウレベル（以下、ＷＬとする）は０［ＨＵ］、ウィンドウ幅（以下、ＷＷとする）は１００［ＨＵ］とした。
【０１６９】
逐次修正には、公知であるサブセット法を用いたＯＳ−ＳＰＳを用い、繰返し回数（更新回数）を６０回、サブセット数を２４に設定した。図３の重み設定領域４３０で設定する重み条件には、一定値を設定し、所望の撮影・画像条件には、所望するノイズ低減割合として２０％を設定した。
【０１７０】
従来方式により得たＣＴ画像６１２は、初期画像６１１と比較して、ノイズが増加したことがわかる。一方、本実施形態の方式で得たＣＴ画像６１３は、初期画像６１１と比較してノイズが低減していることがわかる。また本実施形態の方式で得たＣＴ画像６１３においては、その中心に設定したＲＯＩ内のノイズを計測した結果、初期画像６１１から、約２０％、ノイズが低減していた。すなわち、２０％のノイズ低減割合を達成した。
【０１７１】
このように、本実施形態の方式により、所望するノイズ低減割合を達成したＣＴ画像を取得できることが示された。
【０１７２】
＜＜第二の実施形態＞＞
次に、本実施形態の第二の実施形態を説明する。第一の実施形態では、ユーザが設定する撮影条件（実撮影条件）と代表撮影条件とは、再構成フィルタが異なる。一方、本実施形態では、ユーザが設定する撮影条件（実撮影条件）と代表撮影条件とは、測定投影データに対する補正法が異なる。
【０１７３】
本実施形態のＸ線ＣＴ装置は、基本的に第一の実施形態のＸ線ＣＴ装置１００と同様の構成を有する。しかしながら、上述のように、実撮影条件と代表撮影条件との間で異なる撮影条件に違いがある。このため、本実施形態では、パラメータ決定部３２０の代表ノイズ低減割合変換部３２１の処理が、第一の実施形態と異なる。以下、本実施形態について、第一の実施形態と異なる構成を中心に説明する。
【０１７４】
図１３に従って、本実施形態のパラメータ決定部３２０の代表ノイズ低減割合変換部３２１を説明する。
【０１７５】
本実施形態のパラメータ決定部３２０の代表ノイズ低減割合変換部３２１は、第一の実施形態同様、実撮影条件で指定されたノイズ低減割合を、代表撮影条件のテーブルに入力するための代表ノイズ低減割合に換算する。
【０１７６】
ただし、本実施形態では、第一の実施形態のフィルタノイズテーブル５２０の代わりに、補正法毎に、基準の補正法に対するノイズ増減比を保存する補正法ノイズテーブルをさらに備える。そして、代表ノイズ低減割合変換部３２１は、補正法ノイズテーブルを参照し、所望の低減割合（実低減割合）を代表低減割合に変換する（ステップＳ３２０１）。
【０１７７】
本実施形態の補正法ノイズテーブル５４０の一例を、図１４に示す。本図に示すように、補正法ノイズテーブル５４０は、測定投影データの補正方法５４１毎に、基準の補正法に対するノイズ増減比５４２を格納する。
【０１７８】
第一の実施形態の図３に示すように、補正方法は、その効果により大、中、小と区分され、その区分により指定を受け付ける。従って、補正法ノイズテーブル５４０も、補正方法５４１として、この区分（補正小、補正中、補正大）が格納される。
【０１７９】
ノイズ増減比５４２は、補正効果小、中、大の各投影データ補正手段を用いて、それぞれ、代表的なファントムのＣＴ画像を取得し、当該ＣＴ画像上で計測したＳＤから算出する。本図に示す補正法ノイズテーブル５４０は、補正効果小の値を１．０とした時のノイズ増減比を格納したものである。
【０１８０】
例えば、撮影条件に用いた投影データの補正手段は補正効果を小、代表撮影条件に用いた投影データの補正手段は補正効果を中とした場合のノイズ増減比として、この補正法ノイズテーブル５４０から、０．７７／１．０＝０．７７が得られる。
【０１８１】
なお、補正法ノイズテーブル５４０は、予め作成され、ＨＤＤ装置１３４等に保持される。
【０１８２】
代表ノイズ低減割合変換部３２１は、このノイズ増減比を用い、第一の実施形態と同様の手法で、式（８）に従って、実ノイズ低減割合を代表ノイズ低減割合に変換する。
【０１８３】
例えば、実撮影条件および代表撮影条件の補正効果がそれぞれ小と中であり、実ノイズ低減割合が７５％の場合、上記式（８）より、本実施形態の代表ノイズ低減割合は（１−（１−０．７５）／０．７７）×１００＝６７．５％となる。
【０１８４】
他の処理は、第一の実施形態と同様であるため、ここでは、説明を省略する。
【０１８５】
以上説明したように、本実施形態のＸ線ＣＴ装置１００は、第一の実施形態同様、Ｘ線発生部１２１と、Ｘ線検出部１２２と、逐次再構成部２３３を有する画像生成部１３０と、を備える。また、前記逐次再構成部２３３は、テーブル部３１０と、パラメータ決定部３２０と、逐次修正部３３０と、を備える。そして、前記パラメータ決定部３２０は、代表ノイズ低減割合変換部３２１と、計算テーブル参照部３２２と、重み計算部３２３と、パラメータ変換部３２４と、を備えてもよい。また、前記テーブル部３１０は、基準重み計算部３１２をさらに備えてもよい。
【０１８６】
また、ユーザが設定する撮影条件と前記代表撮影条件とは、測定投影データに対する補正法が異なり、前記逐次再構成部２３３は、前記補正法毎に、基準の補正法に対するノイズ増減比を保存する補正法ノイズテーブル３０をさらに備え、前記代表ノイズ低減割合変換部３２１は、前記補正法ノイズテーブル５４０を参照し、前記所望の低減割合を前記代表低減割合に変換してもよい。
【０１８７】
このように、本実施形態によれば、実際の撮影条件の測定投影データの補正方法が代表撮影条件のそれと異なる場合でも、所望するノイズ低減割合を達成したＣＴ画像を取得することができる。
【０１８８】
すなわち、本実施形態によれば、実撮影条件が代表撮影条件と異なる場合であっても、第一の実施形態同様、実撮影条件の低減割合を代表撮影条件の低減割合に換算し、最適な逐次修正パラメータを得る。従って、第一の実施形態同様、テーブルの作成に要する工数、およびテーブル数を増加させることなく、実撮影条件での最適な逐次修正パラメータを得ることができる。
【０１８９】
また、第一の実施形態同様、画素毎に最適な逐次修正パラメータβＩ（ｊ）を用いて逐次修正できる。従って、所望のノイズ低減割合を実現した、高品質なＣＴ画像を得ることができる。
【０１９０】
なお、本実施形態によれば、実際の撮影条件のノイズ（ＳＤ）を計算することなく、補正法ノイズテーブル５４０を用いてノイズ増減比を取得する。しかし、これに限定されない。本実施形態においても、第一の実施形態同様、実撮影条件および代表撮影条件それぞれの補正手段を用いてＣＴ画像を取得し、予め定めた領域（ＲＯＩ）のＳＤをそれぞれ計測し、その結果を用いてノイズ増減比を算出してもよい。
【０１９１】
図１５（ａ）および図１５（ｂ）に、それぞれ、実撮影条件および代表撮影条件で指定された補正手段を用いて測定投影データを補正し、その結果から得たＣＴ画像６３０、６４０を示す。ファントム中心のＲＯＩ６３１、６４１は、それぞれＳＤ計測に用いた領域である。この計測結果から、式（５）に従って、ノイズ増減比を計算する。
【０１９２】
なお、ＳＤ計測に用いる領域は、ファントム中心の領域６２１、６３１に限定されない。画像６３０および画像６４０の同位置の領域であればよく、例えば、周辺の領域６３２と６４２、領域６３３と６４３、であってもよい。
【０１９３】
さらに、ノイズ増減比は、各画像で一様に設定する必要はない。位置毎に、独立にノイズ増減比を求めてもよい。すなわち、上記例では、領域６３１、６４１におけるノイズ増減比、領域６３２、６４２におけるノイズ増減比、領域６３３、６４３におけるノイズ増減比を、それぞれ算出してもよい。
【０１９４】
さらに、本実施形態においても、第一の実施形態同様、ＣＴ再構成手法は、例えば、公知であるハーフ再構成または一周３６０度分以上の測定投影データを用いる再構成であってもよい。
【０１９５】
＜＜第三の実施形態＞＞
本発明の第三の実施形態を説明する。第一の実施形態では、ユーザが設定する撮影条件（実撮影条件）と代表撮影条件とは、再構成フィルタが異なる。一方、本実施形態では、ユーザが設定する撮影条件（実撮影条件）と代表撮影条件とは、ＣＴ画像生成に用いる測定投影データの範囲および逐次修正に用いる測定投影データの範囲が異なる。
【０１９６】
本実施形態のＸ線ＣＴ装置は、基本的に第一の実施形態のＸ線ＣＴ装置１００と同様の構成を有する。しかしながら、上述のように、実撮影条件と代表撮影条件との間で異なる撮影条件に違いがある。このため、本実施形態では、パラメータ決定部３２０の代表ノイズ低減割合変換部３２１の処理が、第一の実施形態と異なる。以下、本実施形態について、第一の実施形態と異なる構成を中心に説明する。
【０１９７】
図１６に従って、本実施形態のパラメータ決定部３２０の代表ノイズ低減割合変換部３２１を説明する。
【０１９８】
本実施形態においても、パラメータ決定部３２０の代表ノイズ低減割合変換部３２１は、第一の実施形態同様、実撮影条件で指定されたノイズ低減割合を、代表撮影条件のテーブルに入力するための代表ノイズ低減割合に換算する。ただし、本実施形態では、第一の実施形態のように、フィルタノイズテーブル５２０を備えない。
【０１９９】
すなわち、本実施形態の代表ノイズ低減割合変換部３２１は、実撮影条件および代表撮影条件それぞれにおけるＣＴ画像生成時に用いる測定投影データ範囲と、逐次修正に用いる測定投影データ範囲とを用い、前記所望の低減割合を前記代表低減割合に変換する（ステップＳ４２０１）。
【０２００】
なお、測定投影データ範囲とは、上述のように、初期画像の生成および逐次修正による再構成に用いる測定投影データの範囲であり、投影角度で指定される。
【０２０１】
投影データ範囲が３６０度と指定された場合に得られるＣＴ画像のＳＤを１とすると、投影データがＡＮＧ度と指定された場合に得られるＣＴ画像のＳＤは、以下の式（１７）で表される。
ＳＤ＝√（３６０／ＡＮＧ）  ・・・（１７）
例えば、全ての投影データの範囲に対して、投影角度方向、検出素子のチャンネル方向、または、検出素子のスライス列方向にデータを使用する割合である重みが等しい場合、投影データ範囲が１８０度と指定された場合に得られるＣＴ画像のＳＤは、√２倍に、また、２７０度と指定された場合は、１．１５倍に増加する。一方、ＣＴ画像の投影データの範囲に対して、各方向にデータを使用する割合である重みを変更した場合は、重みを考慮したＳＤの増加量を計算する必要がある。
【０２０２】
代表ノイズ低減割合変換部３２１は、以下の式（１８）を用いて、ノイズ増減比を算出する。
【数１８】
2016199716.tif 000015
なお、初期投影データ範囲は、初期画像生成時の測定投影データ範囲であり、逐次修正投影データ範囲は、逐次修正時に用いる測定投影データ範囲である。また、ヘリカルスキャンでは、ＣＴ画像の各画素位置で収集される測定投影データの範囲が変わるため、ＣＴ画像の中心位置のノイズ増減比を代表値として計算してもよいし、各画素位置のノイズ増減比を計算してもよい。
【０２０３】
例えば、実撮影条件における、初期画像生成に用いる測定投影データ範囲を１８０度、逐次修正に用いる測定投影データ範囲を３６０度とし、代表撮影条件における、初期画像生成に用いる測定投影データ範囲を３６０度、逐次修正に用いる測定投影データ範囲を３６０度とする。この場合、式（１８）により算出されるノイズ増減比は、（１／√２）×（１／１）＝１／√２である。
【０２０４】
代表ノイズ低減割合変換部３２１は、このノイズ増減比を用い、第一の実施形態と同様の手法で、式（８）に従って、実ノイズ低減割合を代表ノイズ低減割合に変換する。
【０２０５】
すなわち、本実施形態の代表ノイズ低減割合変換部３２１は、実撮影条件に用いた初期投影データ範囲に対する、代表撮影条件の初期投影データ範囲の割合、および、代表撮影条件の逐次修正投影データ範囲に対する実撮影条件の逐次修正投影データ範囲の割合を計算し、これらを乗算し、ノイズ低減割合を得る。そして、得られたノイズ低減割合を用いて、代表撮影条件の初期画像のノイズ、および、実撮影条件の逐次修正後のＣＴ画像のノイズ（目標ノイズ）を計算する。
【０２０６】
このとき、代表ノイズ低減割合変換部３２１は、ＣＴ画像の投影データの範囲に対して、投影角度方向、検出素子のチャンネル方向、または、検出素子のスライス列方向にデータを使用する割合である重みを寄与するよう構成してもよい。
【０２０７】
他の処理は、第一の実施形態と同様であるため、ここでは、説明を省略する。
【０２０８】
以上説明したように、本実施形態のＸ線ＣＴ装置１００は、第一の実施形態同様、Ｘ線発生部１２１と、Ｘ線検出部１２２と、逐次再構成部２３３を有する画像生成部１３０と、を備える。また、前記逐次再構成部２３３は、テーブル部３１０と、パラメータ決定部３２０と、逐次修正部３３０と、を備える。そして、前記パラメータ決定部３２０は、代表ノイズ低減割合変換部３２１と、計算テーブル参照部３２２と、重み計算部３２３と、パラメータ変換部３２４と、を備えてもよい。また、前記テーブル部３１０は、基準重み計算部３１２をさらに備えてもよい。
【０２０９】
ユーザが設定する撮影条件と前記代表撮影条件とは、前記ＣＴ画像生成に用いる前記測定投影データの範囲および逐次修正に用いる前記測定投影データの範囲が異なり、前記代表ノイズ低減割合変換部３２１は、前記ユーザが設定する撮影条件および前記代表撮影条件それぞれにおける前記ＣＴ画像生成時に用いる測定投影データ範囲と逐次修正に用いる前記測定投影データ範囲とを用い、前記所望の低減割合を前記代表低減割合に変換してもよい。
【０２１０】
このように、本実施形態によれば、実際の撮影条件で、代表撮影条件と異なる投影データの範囲が使用される場合でも、所望するノイズ低減割合を達成したＣＴ画像を取得することができる。
【０２１１】
すなわち、本実施形態によれば、実撮影条件が代表撮影条件と異なる場合であっても、第一の実施形態同様、実撮影条件の低減割合を代表撮影条件の低減割合に換算し、最適な逐次修正パラメータを得る。従って、第一の実施形態同様、テーブルの作成に要する工数、およびテーブル数を増加させることなく、実撮影条件での最適な逐次修正パラメータを得ることができる。
【０２１２】
また、第一の実施形態同様、画素毎に最適な逐次修正パラメータβＩ（ｊ）を用いて逐次修正できる。従って、所望のノイズ低減割合を実現した、高品質なＣＴ画像を得ることができる。
【０２１３】
なお、本実施形態では、上記式（１８）により、ノイズ増減比を計算したが、この手法に限定されない。例えば、予め、使用する測定投影データ範囲の組み合わせに応じたテーブルを作成し、それを参照するよう構成してもよい。
【０２１４】
なお、公知である重み付き解析的再構成のように、投影角度毎に逆投影に寄与する重みを変える場合、同様にして各投影角度の重みから相対的なＳＤの増減比を計算する必要がある。この場合、投影角度の方向だけでなく、チャネル方向またはスライス方向に重みを変える場合も、同様にして相対的なＳＤの増減比を計算する必要がある。
【０２１５】
＜＜第四の実施形態＞＞
次に、本発明の第四の実施形態を説明する。上記各実施形態では、逐次修正の収束段階のノイズが等しいことを前提とする。しかしながら、実際には、設定される更新回数によっては、収束段階に達しないこともある。そこで、本実施形態では、更新回数を予め定め、実撮影条件で得た初期画像が、その繰返し回数（更新回数）で収束段階のノイズに到達しない場合、逐次修正処理前に、ユーザに通知する。
【０２１６】
本実施形態のＸ線ＣＴ装置は、基本的に第一の実施形態のＸ線ＣＴ装置１００と同一の構成を有する。ただし、上述のようにユーザに通知する必要があるため、本実施形態の逐次再構成部２３３は、図１７に示すように、さらに、適否判別部３４１と判別結果通知部３４２と、を備える。
【０２１７】
[適否判別部]
適否判別部３４１は、実撮影条件に従って得た測定投影データから生成したＣＴ画像（初期画像）の適否を判別する。そして、適否判別部３４１は、逐次修正において、ＣＴ画像（初期画像）が、予め定めた繰り返し回数内で収束しないと判別した場合、不適と判別する。すなわち、与えられた初期画像が、上記各実施形態の手法で逐次修正処理が可能な初期画像であるか否かを判別し、否である場合、その旨、ユーザに通知する。
【０２１８】
判別は、与えられた初期画像のノイズと、収束段階のノイズ(以下、収束ノイズ）とから得られる実ノイズ低減割合に応じて行う。すなわち、この実ノイズ低減割合が、所定回数で収束ノイズとならない場合、適否判別部３４１は、否と判別する。
【０２１９】
適否判別部３４１は、予め定めた前記低減割合の閾値によって適否を判別する。判定の基準とする閾値は、更新回数に応じて予め設定する。設定する閾値は、所定回数で、初期画像のノイズが収束段階のノイズとなる、ノイズ低減割合の上限値（最大閾値ＴＨ％）および下限値（最小閾値ＴＬ％）である。
【０２２０】
適否判別部３４１は、実際に初期画像が与えられた際、当該初期画像のノイズと収束段階のノイズとから、そのノイズ低減割合（Ａ％）を算出し、それが最大閾値ＴＨと最小閾値ＴＬとの間であるか否かを判別し、両閾値外である場合、その旨、ユーザに通知する。すなわち、ＴＬ≦Ａ≦ＴＨの場合、処理可能初期画像と判別し、Ａ＜ＴＬ、または、Ａ＞ＴＨの場合、否と判別する。
【０２２１】
判定に用いる閾値（ＴＨおよびＴＬ）は、例えば、図６（ａ）に示すように、事前に代表的なファントムを用いて、複数の異なる基準パラメータβｂを用いて逐次修正した結果から求める。
【０２２２】
具体例で説明する。図１８は、代表的な逐次修正パラメータを用いた時の、ノイズの変化を示すグラフ５５２である。横軸は各繰返し回数（更新回数）であり、縦軸はノイズを示すである。なお、繰返し回数（更新回数）０回は、初期画像を意味する。また、グラフ５５２では、繰返し回数（更新回数）６０回以上でＳＤが変動しない収束段階となる。収束ノイズは５５４とする。
【０２２３】
このとき、与えられた初期画像のノイズ５５３と収束ノイズ５５４とから算出されるノイズ低減割合Ａ％が、最大閾値ＴＨ％より大きい場合、繰返し回数（更新回数）６０回で収束段階のノイズに到達しない。本実施形態の適否判別部３４１は、Ａ％とＴＨ％とを比較し、判定を行い、ノイズ低減割合が、予め定めた最大閾値ＴＨ％より大きい場合、不適（否）と判別する。
【０２２４】
例えばグラフ５５５は、ノイズ低減割合Ｂ％が最大閾値ＴＨ％以上となる初期画像を用いた為、繰返し回数（更新回数）６０回で収束ノイズ５５４に到達しない。
【０２２５】
反対に、初期画像のノイズが５５６である場合、このノイズ５５６と収束ノイズ５５４とから決定するノイズ低減割合Ｃ％が、最小閾値ＴＬ％以下の場合も、同様に繰返し回数（更新回数）６０回で収束ノイズ５５４に到達しない。この場合も、適否判別部３４１は、不適（否）と判別する。
【０２２６】
［判別結果通知部］
判別結果通知部３４２は、適否判別部３４１よる判別結果が不適な場合、ユーザに通知する。ユーザへの通知は、例えば、エラーメッセージなどの表示とする。しかしながら、これに限定されない。例えば、表示等による出力ではなく、操作を停止するといった処理を行うことにより、ユーザに通知してもよい。
【０２２７】
さらに、初期画像を、閾値内に納まる初期画像に変更するようユーザに通知してもよいし、自動的に変更するよう構成してもよい。
【０２２８】
例えば、自動的に変更するよう構成する場合、逐次再構成部２３３は、判別結果通知部３４２に加え、あるいは、判別結果通知部３４２の代わりに、画像変更部３４３を備える。
【０２２９】
画像変更部３４３は、適否判別部３４１による判別結果が不適な場合、ＣＴ画像（初期画像）を変更する。このとき、画像変更部３４３は、低減割合が閾値内であって、判別対象のＣＴ画像（初期画像）の低減割合に最も近いＣＴ画像に変更する。すなわち、変更する初期画像は、ノイズ低減割合が、閾値内に納まる初期画像で、初期画像のノイズが最も近いものとする。
【０２３０】
また、本実施形態では、１種類の繰返し回数（更新回数；６０回）について、閾値を保持しておく場合を例にあげて説明したが、これに限定されない。例えば、繰返し回数（更新回数）１００回、２００回等、複数種類の繰り返し回数について、それぞれ、閾値を保持しておくよう構成してもよい。
【０２３１】
以上説明したように、本実施形態のＸ線ＣＴ装置１００は、第一の実施形態同様、Ｘ線発生部１２１と、Ｘ線検出部１２２と、逐次再構成部２３３を有する画像生成部１３０と、を備える。また、前記逐次再構成部２３３は、テーブル部３１０と、パラメータ決定部３２０と、逐次修正部３３０と、を備える。そして、前記パラメータ決定部３２０は、代表ノイズ低減割合変換部３２１と、計算テーブル参照部３２２と、重み計算部３２３と、パラメータ変換部３２４と、を備えてもよい。また、前記テーブル部３１０は、基準重み計算部３１２をさらに備えてもよい。
【０２３２】
前記逐次再構成部２３３は、前記測定投影データから生成したＣＴ画像の適否を判別する適否判別部３４１をさらに備え、前記適否判別部３４１は、前記逐次修正において、前記ＣＴ画像が、予め定めた繰り返し回数内で収束しないと判別した場合、不適と判別してもよい。また、前記逐次再構成部２３３は、前記適否判別部３４１による判別結果が不適な場合、ユーザに通知する判別結果通知部３４２、および、前記適否判別部３４１による判別結果が不適な場合、前記ＣＴ画像を変更する画像変更部３４３の少なくとも一方をさらに備えてもよい。
【０２３３】
このように、本実施形態によれば、実撮影条件で得た初期画像が、予め定めた更新回数で、収束段階のノイズに到達しない場合、ユーザに通知する。そして、ユーザは、撮影条件等を変更し、初期画像を、予め定めた更新回数で収束段階のノイズに到達するものに変更することができる。あるいは、初期画像が、自動的に、予め定めた更新回数で収束段階のノイズに到達するものに変更される。
【０２３４】
従って、本実施形態によれば、上記第一から第三の実施形態で説明した手法を用いるにあたり、初期画像が必ず、必要な条件を満たすことになる。よって、本実施形態によれば、確実に所望のノイズ低減割合を実現した、高品質のＣＴ画像を得ることができる。
【０２３５】
特に、画像変更部３４３を備える場合、実撮影条件で得た初期画像に近いノイズをもつ画像を用いて逐次修正を実施することで、実撮影条件の初期画像の画質に近く、かつ所望のノイズ低減割合の画像を取得できる。
【０２３６】
＜重み条件の変形例＞
なお、上記各実施形態では、全ての検出素子の重みを一定にした場合を例にあげて説明した。しかしながら、これに限定されない。例えば、重み設定領域４３０を介して、検出素子で検出した光子数に応じた重みを用いてもよい。
【０２３７】
この場合、重み条件が異なるため、基準重み計算部３１２、重み計算部３２３およびパラメータ変換部３２４の処理が異なる。
【０２３８】
［基準重み計算部］
上記各実施形態では、各検出素子の重みを一定としていたため、重み画像は、ＷＩ（ｊ）と算出される。しかしながら、本変形例では、各検出素子の重みを、光子数に応じたものとする。各検出素子ｉの重みをＷＰ（ｉ）とする。
【０２３９】
本変形例の基準重み計算部３１２は、この各検出素子ｉの重みＷＰ（ｉ）に、上記実施形態同様、ＣＴの画像再構成と同じ処理を施し、基準領域の各画素位置の重みＷＩＰ（ｊ）を得る。そして、上記各実施形態同様、各重み値から得た平均値、最大値、最小値等の値を、基準重みＷｂとする。なお、ｊは、画素位置を表す。
【０２４０】
［重み計算部］
本変形例の重み計算部３２３は、この各検出素子ｉの重みＷＰ（ｉ）に、上記実施形態同様、ＣＴの画像再構成と同じ処理を施し、重み画像ＷＩＰ（ｊ）を得る。
【０２４１】
［パラメータ変換部］
本変形例のパラメータ変換部３２４は、上記各実施形態と同様の手法で算出された基準パラメータβｂと、本変形例の手法で算出した基準重みＷｂと、各画素位置の重み値ＷＩＰ（ｊ）と、を用いて、以下の式（１９）に従って、本変形例の逐次修正パラメータの分布βＩＰ（ｊ）を算出する。
【数１９】
2016199716.tif 000016
【０２４２】
なお、逐次修正部３３０は、この算出された各画素位置ｊの逐次修正パラメータβＩＰ（ｊ）を用い、逐次修正を行う。逐次修正処理の詳細は、上記実施形態と同様である。
【０２４３】
図１９に、水で満たされた直径３０ｃｍ、高さ１００ｃｍの円柱ファントムを撮影した結果から、検出素子の光子数を用いて計算した重み画像ＷＩＰ（ｊ）６２０を示す。ファントムは、回転中心から下に１０ｃｍ移動させた位置に配置した。本図に示すように、各位置ｊ１、およびｊ２の重み画像ＷＩＰ（ｊ１）、およびＷＩＰ（ｊ２）は、異なる値をもつ。
【０２４４】
このように、重み条件として、検出素子の光子数を設定した場合であっても、逐次修正時の重みとして、位置（画素）毎の重み値（各位置のパラメータβ画像βＩＰ（ｊ））を用いることにより、所望のノイズ低減割合を達成したＣＴ画像を取得することができる。
【０２４５】
＜その他の変形例＞
上記実施形態１−４は、組み合わせて用いることができる。組み合わせることにより、実撮影条件と、代表撮影条件と複数の条件が異なる場合であっても、上記各実施形態の手法で、実撮影条件の所望のノイズ低減割合を、代表撮影条件の目標ノイズ低減割合に換算し、簡易な構成で所望のノイズ低減割合を実現することができる。
【０２４６】
また、上記第一から第三の実施形態では、実撮影条件と異なる撮影条件として、再構成フィルタ、測定投影データに対する補正法、測定投影データ範囲を取り上げて説明した。しかし、異なる撮影条件は、これらに限定されない。撮影条件を変化させることによるノイズの変化量を把握可能な撮影条件であれば、どのような撮影条件であっても対応できる。
【０２４７】
また、上記各実施形態および変形例では、計算時間の短縮の為、テーブル部３１０およびパラメータ決定部３２０の処理を、予め実施する場合を例にあげて説明した。しかしながら、テーブル部３１０およびパラメータ決定部３２０による処理は、必ずしも事前に実施する必要はない。逐次修正部３３０による処理の前であれば、ユーザが指定した撮影条件、または再構成条件に応じて、撮影時、または画像生成時に実施してもよい。
【０２４８】
また、例えば、重み計算部３２３およびパラメータ変換部３２４を、Ｐｒｉｏｒ計算部３３５の一部として加えてもよい。このとき、式（９）に示す重み画像ＷＩ（ｊ）は、式（１５）右辺の分母のＰ２以外の部分である、以下の項（２０）と等価である。
【数２０】
2016199716.tif 000017
従って、Ｐｒｉｏｒ計算部３３５では、式（９）のβＩ（ｊ）を、式（１５）（式（１３）および式（１４））のβに代入することにより、重み計算部３２３、およびパラメータ変換部３２４を省略することができる。これにより、重複するＷＩ（ｊ）の重み計算に要する計算量を削減できる。
【０２４９】
また、上記各実施形態および変形例において、入力部１１０および画像生成部１３０は、撮影部１２０を備える本体装置（Ｘ線ＣＴ装置１００）と必ずしも一体に構成する必要はなく、撮影部１２０とは独立した情報処理装置上に構築されてもよい。この場合、当該情報処理装置と、Ｘ線ＣＴ装置１００とは、ネットワークを介して接続してもよい。また、画像生成部１３０のみを、測定投影データを処理する処理装置として、独立した情報処理装置上に構築してもよい。
【０２５０】
この場合、Ｘ線ＣＴ装置１００は、逐次修正を終えたＣＴ画像を、ネットワークアダプタを用いて、ローカルエリアネットワーク、電話回線、インターネット等のネットワークを介して外部の端末に送信可能なよう構成する。
【０２５１】
また、上記各実施形態および変形例では、画質を制御する逐次修正パラメータとして、平滑化項の係数であるβ（βＩ（ｊ））を用いているが、これに限定されない。逐次修正パラメータは、Ｌｉｋｅｌｉｈｏｏｄ計算とＰｒｉｏｒ計算との割合を決定するパラメータであれば、Ｌｉｋｅｌｉｈｏｏｄ計算の係数、または二つの計算の係数を用いてもよい。
【０２５２】
また、上述の式（１０）で示した逐次再構成手法は一例であり、これに限定されない。公知であるＯＳ−ＳＰＳ、ＯＳ−ＳＰＳ−ＴＶ、ＰＷＬＳ、ＯＳ−ＰＷＬＳ、ＡＳＩＲＴ、ＭＳＩＲＴ、ＧＲＡＤＹ、ＣＯＮＧＲ、ＡＲＴ、ＳＡＲＴ、ＳＡＲＴ−ＴＶ、ＯＳ−ＳＡＲＴ、ＯＳ−ＳＡＲＴ−ＴＶ、ＭＬ−ＥＭ、ＯＳ−ＥＭ、ＦＩＲＡ、ＲＡＭＬＡ、ＤＲＡＭＡ等、他の手法を用いてもよい。
【０２５３】
また、上記各実施形態および変形例では、寝台１２５およびガントリ１２３が静止した状態であるコンベンショナルスキャン方式で測定投影データを取得する場合を例にあげて説明したが、これに限定されない。例えば、寝台１２５の動作と停止を一定間隔で順に繰り返しながらコンベンショナルスキャンを行うステップアンドシュート方式や、寝台１２５を動かしながら撮影する螺旋スキャン方式でであってもよい。
【０２５４】
更に、上記各実施形態および変形例では、一例として生体用のＸ線ＣＴ装置１００を例にあげて説明したが、爆発物検査や製品検査等の非破壊検査を目的としたＸ線ＣＴ装置であってもよい。
【０２５５】
また、上記各実施形態および変形例では、一例として公知の第３世代のマルチスライスＸ線ＣＴ装置１００を例にあげて説明したが、Ｘ線ＣＴ装置は、公知の第１、第２、第４世代のＸ線ＣＴ装置であってもよい。また、公知のシングルスライスＸ線ＣＴ装置やエレクトロンビームＣＴであってもよい。
【符号の説明】
【０２５６】
１００…Ｘ線ＣＴ装置、１０１…被写体、１１０…入力部、１１１…キーボード、１１２…マウス、１１３…メモリ、１１４…ＣＰＵ、１１５…ＨＤＤ装置、１１６…データバス、１３２…メモリ、１２０…撮影部、１３３…ＣＰＵ、１２１…Ｘ線発生部、１２２…Ｘ線検出部、１２３…ガントリ、１２４…回転板、１２５…寝台、１２６…開口部、１２７…ガントリ制御器、１２８…Ｘ線制御器、１２９…寝台制御器、１３０…画像生成部、１３１…ＤＡＳ、１３２…メモリ、１３３…ＣＰＵ、１３４…ＨＤＤ装置、１３５…モニタ、１３６…データバス、２１１…撮影条件入力部、２２１…撮影制御部、２２２…信号取得部、２３１…信号収集部、２３２…補正処理部、２３３…逐次再構成部、２３４…画像表示部、３１０…テーブル部、３１１…テーブル計算部、３１２…基準重み計算部、３１３…計算テーブル保存部、３２０…パラメータ決定部、３２１…代表ノイズ低減割合変換部、３２２…計算テーブル参照部、３２３…重み計算部、３２４…パラメータ変換部、３３０…逐次修正部、３３１…解析的再構成部、３３２…順投影部、３３３…差分部、３３４…逆投影処理部、３３５…Ｐｒｉｏｒ計算部、３３６…画像修正部、４００…撮影条件受付画面、４１０…Ｘ線条件設定領域、４２０…再構成範囲設定領域、４３０…重み設定領域、４４０…撮影部位設定領域、４５０…撮影・画像設定領域、４６０…再構成フィルタ設定領域、４７０…投影データ補正設定領域、４８０…投影データ範囲設定領域、５１０…ＣＴ画像、５１１…ＲＯＩ、５１２…計測値、５１３…近似曲線、５１４…ＲＯＩ、５２０…フィルタノイズテーブル、５２１…再構成フィルタの種類、５３１…ノイズ、５３２…実撮影条件による初期画像のノイズ、５３３：代表撮影条件による初期画像のノイズ、５３４：ノイズの収束値、５４０：補正法ノイズテーブル、５４１：測定投影データの補正方法、５４２：ノイズ増減比、５５２：ノイズの変化グラフ、５５３：初期画像のノイズ、５５４：収束ノイズ、５５５：ノイズの変化グラフ、５５６：初期画像のノイズ、６１１…初期画像、６１２…ＣＴ画像、６１３…ＣＴ画像、６２０：重み画像、６３０：ＣＴ画像、６３１：領域、６３２：領域、６３３：領域、６４０：ＣＴ画像、６４１：領域、６４２：領域、６４３：領域、７１０…パラメータテーブル、７１１…更新回数、７１２…再構成フィルタ、７１３…テーブル
【図１】
2016199716.tif 000018
【図２】
2016199716.tif 000019
【図３】
2016199716.tif 000020
【図４】
2016199716.tif 000021
【図５】
2016199716.tif 000022
【図６】
2016199716.tif 000023
【図７】
2016199716.tif 000024
【図８】
2016199716.tif 000025
【図９】
2016199716.tif 000026
【図１０】
2016199716.tif 000027
【図１１】
2016199716.tif 000028
【図１２】
2016199716.tif 000029
【図１３】
2016199716.tif 000030
【図１４】
2016199716.tif 000031
【図１５】
2016199716.tif 000032
【図１６】
2016199716.tif 000033
【図１７】
2016199716.tif 000034
【図１８】
2016199716.tif 000035
【図１９】
2016199716.tif 000036
【国際調査報告】
2016199716.tif 000037
2016199716.tif 000038