(57)【要約】  ノイズ低減効果を維持しつつ処理時間を短縮することが可能な演算装置、X線CT装置、及び画像再構成方法を提供するために、従来の逐次近似画像再構成方法(SPS法)では、繰り返し行われる画像更新の都度、更新画像を順投影し、順投影データを元の被写体投影データと差分しこれを逆投影することで差分画像を得て、順投影処理と逆投影処理とを繰り返し行っていたが、本発明では、演算時間を要する順投影処理と逆投影処理を、更新画像と基準画像との差分という比較的演算量の少ない処理に置き換え、その結果、逐次近似画像再構成処理の演算量が大幅に減少し、処理時間を短縮する事が可能となる。
2016132880.tif 000002
【特許請求の範囲】
【請求項１】
被写体の周囲の各方向から照射され、前記被写体を透過したX線の情報に基づいて被写体投影データを作成する投影データ作成部と、
前記被写体投影データから基準画像を作成する基準画像作成部と、
更新処理におけるノイズ低減強度を調整する規格化係数を画素毎に定義した画像である規格化画像を作成する規格化画像作成部と、
前記基準画像と前記規格化画像を用いて、反復処理を所定回数実行する反復処理部と、を備え、
前記反復処理部は、
前記基準画像と前記更新処理により得た更新画像とを差分し、差分画像を作成する差分画像作成部と、
前記差分画像及び前記規格化画像を用いて前記更新処理を実行し、更新画像を作成する画像更新部と、
を有することを特徴とする演算装置。
【請求項２】
前記基準画像作成部は、フィルタ逆投影法により基準画像を作成することを特徴とする請求項1に記載の演算装置。
【請求項３】
前記基準画像作成部は、前記被写体投影データに対しノイズ低減処理を施して補正投影データを作成し、作成した補正投影データを画像再構成して第1補正断層像を作成し、作成した前記第1補正断層像に鮮鋭化処理を施した第2補正断層像を前記基準画像として作成することを特徴とする請求項1に記載の演算装置。
【請求項４】
前記反復処理部は、
前記差分画像作成部における基準画像及び前記画像更新部における規格化画像として、それぞれサイズを拡大した拡大断層像及び拡大規格化画像を使用して前記差分画像の作成及び前記更新処理を所定回数反復して実行することにより拡大更新画像を作成し、前記拡大更新画像を元のサイズに縮小した縮小画像を作成することを特徴とする請求項1に記載の演算装置。
【請求項５】
前記基準画像と前記更新画像とを所定の比率で合成した画像を、前記差分画像作成部における基準画像として使用することを特徴とする請求項1に記載の演算装置。
【請求項６】
前記規格化画像作成部により作成された規格化画像を前記画像更新部及び前記差分画像作成部における基準画像として前記更新処理を実行することにより更新規格化画像を作成し、作成した前記更新規格化画像を前記画像更新部における規格化画像として使用することを特徴とする請求項1に記載の演算装置。
【請求項７】
前記基準画像上に関心領域を設定する関心領域設定部と、
前記関心領域の内外で異なる規格化係数を設定する設定部と、を備え、
前記規格化画像作成部は、設定された前記規格化係数を画素毎に定義した規格化画像を作成することを特徴とする請求項1に記載の演算装置。
【請求項８】
前記規格化画像作成部は、前記被写体投影データと空気投影データとを加算して逆対数変換係数を乗じて指数をとることにより線量換算データを求め、前記線量換算データを画像再構成することにより、前記規格化画像を作成することを特徴とする請求項1に記載の演算装置。
【請求項９】
前記規格化画像作成部は、前記被写体の断層像のCT値を非負値化した非負値化画像を前記規格化画像とすることを特徴とする請求項1に記載の演算装置。
【請求項１０】
前記規格化画像作成部は、CT値と規格化係数との相関を設定した相関関数に基づいて前記被写体の断層像から規格化画像を算出することを特徴とする請求項1に記載の演算装置。
【請求項１１】
請求項1に記載の演算装置を備えたX線CT装置。
【請求項１２】
演算装置が、
被写体の周囲の各方向から照射され、前記被写体を透過したX線の情報に基づいて被写体投影データを作成する投影データ作成ステップと、
前記被写体投影データから基準画像を作成する基準画像作成ステップと、
更新処理におけるノイズ低減強度を調整する規格化係数を画素毎に定義した画像である規格化画像を作成する規格化画像作成ステップと、
前記基準画像と前記規格化画像を用いて、反復処理を所定回数実行する反復処理ステップと、を含み、
前記反復処理ステップにおいて、
前記基準画像と前記更新処理により得た更新画像とを差分し、差分画像を作成する差分画像作成ステップと、
前記差分画像及び前記規格化画像を用いて前記更新処理を実行し、更新画像を作成する画像更新ステップと、
を含むことを特徴とする画像再構成方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、演算装置、X線CT装置、及び画像再構成方法に関し、詳細には、逐次近似再構成法による画像再構成における高速ノイズ低減処理技術に関する。
【背景技術】
【０００２】
近年、被曝線量の少ないCT検査を実施するため、低線量でもノイズが少ない画像が得られる逐次近似法による画像再構成が行われている。【０００３】
特許文献1には、式(1)に示すように、逐次近似法の一つであるSeparable Paraboloidal Surrogate法(SPS法)で用いられる更新式が開示されている。
2016132880.tif 000003
【０００４】
式(1)において、μは更新画素値、yは元の投影値、βは正規化係数、dは規格化係数、wは画素係数、iは検出素子番号、j及びkは画素番号、nは反復回数である。γは投影値の画素値に変換する行列の要素、aは投影値を画素値に変換する行列の要素であり、Σγijは逆投影処理、Σaijは順投影処理である。ポテンシャル関数(罰則項とも言う)の代替関数Ψの一次成分がφ、2次成分がωである。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】国際公開2013/008702号
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、上述の式(1)に示す更新式には、更新式中に演算負荷の高い行列演算である逆投影処理Σγijと順投影処理Σaijが含まれる。そのため演算量が多く、処理時間が長くなるという問題がある。【０００７】
本発明は、以上の問題点に鑑みてなされたものであり、その目的とするところは、逐次近似画像再構成法において、ノイズ低減効果を維持しつつ処理時間を短縮することが可能な演算装置、X線CT装置、及び画像再構成方法を提供することである。
【課題を解決するための手段】
【０００８】
前述した目的を達成するために本発明は、被写体の周囲の各方向から照射され前記被写体を透過したX線の情報に基づいて被写体投影データを作成する投影データ作成部と、前記被写体投影データから基準画像を作成する基準画像作成部と、更新処理におけるノイズ低減強度を調整する規格化係数を画素毎に定義した画像である規格化画像を作成する規格化画像作成部と、前記基準画像と前記規格化画像を用いて、反復処理を所定回数実行する反復処理部と、を備え、前記反復処理部は、前記基準画像と前記更新処理により得た更新画像とを差分し、差分画像を作成する差分画像作成部と、前記差分画像及び前記規格化画像を用いて前記更新処理を実行し、更新画像を作成する画像更新部と、を有することを特徴とする演算装置である。【０００９】
また、本発明の演算装置を備えたX線CT装置である。【００１０】
また、演算装置が、被写体の周囲の各方向から照射され前記被写体を透過したX線の情報に基づいて被写体投影データを作成する投影データ作成ステップと、前記被写体投影データから基準画像を作成する基準画像作成ステップと、更新処理におけるノイズ低減強度を調整する規格化係数を画素毎に定義した画像である規格化画像を作成する規格化画像作成ステップと、前記基準画像と前記規格化画像を用いて、反復処理を所定回数実行する反復処理ステップと、を含み、前記反復処理ステップにおいて、前記基準画像と前記更新処理により得た更新画像とを差分し、差分画像を作成する差分画像作成ステップと、前記差分画像及び前記規格化画像を用いて前記更新処理を実行し、更新画像を作成する画像更新ステップと、を含むことを特徴とする画像再構成方法である。
【発明の効果】
【００１１】
本発明によれば、演算負荷の高い行列演算である逆投影処理と順投影処理を更新式から省き、演算負荷の低い画像処理に置き換えるため、ノイズ低減効果を維持しつつ処理時間を短縮することが可能な演算装置、X線CT装置、及び画像再構成方法を提供できる。
【図面の簡単な説明】
【００１２】
【図１】X線CT装置1の全体構成を示す外観図
【図２】X線CT装置1の内部構成を示すブロック図
【図３】演算装置202の機能ブロック図
【図４】処理全体の流れを示すフローチャート
【図５】図4のステップS102の逐次近似画像再構成処理の流れを示すフローチャート
【図６】基準画像作成処理(A)の流れを示すフローチャート
【図７】規格化画像作成処理の流れを示すフローチャート
【図８】基準画像作成処理(B)の流れを示すフローチャート
【図９】画像更新処理(A)の流れを示すフローチャート
【図１０】画像更新処理(B)の流れを示すフローチャート
【図１１】ノイズ低減画像表示画面51の例
【図１２】断層像への関心領域60の設定例
【図１３】CT値と規格化係数との相関を示す相関関数80の一例
【発明を実施するための形態】
【００１３】
以下、図面に基づいて本発明の好適な実施形態を詳細に説明する。まず、図1及び図2を参照しながら、X線CT装置1のハードウエア構成を説明する。
【００１４】
図1及び図2に示すように、X線CT装置1は、大別してスキャナ10及び操作ユニット20から構成される。スキャナ10は、ガントリ100及び寝台101を有する。ガントリ100は、X線発生装置102、X線検出装置103、コリメータ104、高電圧発生装置105、データ収集装置106、及びガントリ駆動装置107等を有する。操作ユニット20は、中央制御装置200、入出力装置201、及び演算装置202等を有する。
【００１５】
操作者は、操作ユニット20の入出力装置201を用いて撮影条件や再構成条件等を入力する。撮影条件は、例えばX線ビーム幅、寝台送り速度、管電流、管電圧、撮影範囲(体軸方向範囲)、周回当たりの撮影ビュー数等である。また再構成条件は、例えば関心領域、FOV(Field Of View)、再構成フィルタ関数等である。入出力装置201は、CT画像等を表示する表示装置211、マウス、トラックボール、キーボード、タッチパネル等の入力装置212、データを記憶する記憶装置213等を含む。
【００１６】
中央制御装置200は、CPU(Central Processing Unit)、ROM(Read Only Memory)、RAM(Random Access Memory)等を備えたコンピュータであり、X線CT装置1全体の動作を制御する。中央制御装置200は、撮影条件や再構成条件に基づいて撮影に必要な制御信号をスキャナ10の各装置に送信する。
【００１７】
中央制御装置200からの撮影スタート信号を受けて撮影が開始されると、高電圧発生装置105は中央制御装置200からの制御信号に基づいてX線発生装置102に所定の大きさの管電圧、管電流を印加する。X線発生装置102は、印加された管電圧に応じたエネルギーの電子を陰極から放出し、放出された電子がターゲット(陽極)に衝突することによって電子エネルギーに応じたエネルギーのX線を被写体3に照射する。
【００１８】
X線発生装置102から照射されるX線は、コリメータ104によって照射領域が制限される。コリメータ104の開口幅は中央制御装置200からの制御信号に基づいて制御される。
【００１９】
X線CT装置1は、2次元方向に検出素子が配列されるX線検出装置103を用いるマルチスライスCTと、検出素子が1列すなわち1次元方向(チャネル方向のみ)に配列されるX線検出装置103を用いるシングルスライスCTに大別される。マルチスライスCTでは、X線検出装置103に合わせてX線源であるX線発生装置102から円錐状、もしくは角錐状に広がるX線ビームが照射される。シングルスライスCTでは、X線発生装置102から扇状に広がるX線ビームが照射される。
【００２０】
X線発生装置102から照射され、コリメータ104によって照射領域が制限されたX線は被写体3内の各組織においてX線減弱係数に応じて吸収(減衰)され、被写体3を通過し、X線発生装置102に対向する位置に配置されたX線検出装置103によって検出される。
【００２１】
寝台101は、被写体3を寝載する天板、上下動装置、及び天板駆動装置を備え、中央制御装置200からの制御信号に基づき天板高さを上下に昇降したり、体軸方向へ前後動したり、体軸と垂直方向かつ床面に対し平行な方向(左右方向)へ左右動する。撮影中において寝台101は、中央制御装置200により決定された寝台移動速度及び移動方向で天板を移動させる。
【００２２】
ガントリ駆動装置107は、中央制御装置200からの制御信号に基づいて、ガントリ100の回転盤を周回させる。
【００２３】
X線検出装置103は、例えばシンチレータとフォトダイオードの組み合わせによって構成されるX線検出素子群をチャンネル方向(周回方向)及び列方向(体軸方向)に2次元配列したものである。X線検出装置103は、被写体3を介してX線発生装置102に対向するように配置される。X線検出装置103はX線発生装置102から照射されて被写体を透過したX線量を検出し、データ収集装置106に出力する。
【００２４】
データ収集装置106は、X線検出器103の個々のX線検出素子により検出されるX線量の情報を収集し、デジタル信号に変換し、透過X線情報として操作ユニット20の演算装置202に順次出力する。
【００２５】
演算装置202は、再構成処理装置221、及び画像処理装置222等を備える。
【００２６】
再構成処理装置221は、データ収集装置106により収集された透過X線情報を取得し、画像を再構成するために必要な投影データを作成する。また、再構成処理装置221は投影データを用いて被写体3の断層像(CT画像)を再構成する。
【００２７】
本発明において、再構成処理装置221は、後述する逐次近似画像再構成処理を行う。
これにより、ノイズが低減されたCT画像を高速に再構成する。本発明の演算装置202(再構成処理装置221)が実行する逐次近似画像再構成処理の詳細については後述する。再構成処理装置221は、生成したCT画像を記憶装置213に保存するとともに表示装置211に表示する。
【００２８】
画像処理装置222は、再構成処理装置221により作成され、記憶装置213に記憶されたCT画像に対して画像処理を行う。また、画像処理後の画像を表示装置211に表示するとともに記憶装置213に保存する。
【００２９】
次に、図3を参照しながら、本発明のX線CT装置1の機能構成について説明する。
【００３０】
X線CT装置1は主要な機能構成として、データ記憶部31、透過線量収集部32、投影データ作成部33、基準画像作成部34、規格化画像作成部35、反復処理部30、及び画像表示部38を有する。反復処理部30は画像更新部36及び差分画像作成部37を有する。
【００３１】
データ記憶部31は、予め計測または計算して得られた空気投影データを記憶する領域であり、例えば操作ユニット20の記憶装置213に設けられる。データ記憶部31は、空気投影データを規格化画像作成部35へ送信する。空気投影データとは、被写体3及び寝台101のない状態での撮影により計測された投影データである。
【００３２】
透過線量収集部32は、X線発生装置102から被写体3に照射され、被写体3及び寝台101を透過し、X線検出装置103によって検出されたX線の線量情報である透過X線情報を収集する。透過線量収集部32は、収集した透過X線情報を投影データ作成部33に送信する。
【００３３】
投影データ作成部33は、透過線量収集部32から送信される透過X線情報を取得し、所定の対数変換係数を乗じて対数変換し、画像を再構成するために必要な被写体投影データを作成する。作成された被写体投影データは基準画像作成部34及び規格化画像作成部35に送信される。
【００３４】
基準画像作成部34は、取得した被写体投影データを用いて被写体断層像を再構成し、この被写体断層像を後述する画像更新処理の基準画像とする。再構成された基準画像は、差分画像作成部37へ送信される。本発明の逐次近似画像再構成処理では、この基準画像を繰り返し行われる画像更新処理の初回のみならず、毎回参照する。
【００３５】
規格化画像作成部35は、画像更新処理におけるノイズ低減強度を調整する規格化係数を画素毎に定義した画像である規格化画像を作成する。規格化画像の具体的な作成手順は、例えば規格化画像作成部35は、空気投影データと被写体投影データとを取得して加算し、逆対数変換係数を乗じ、指数をとって線量換算投影データを作成する。この線量換算投影データを画像再構成したものを規格化画像とする。
【００３６】
規格化画像作成部35は、作成した規格化画像を画像更新部36へ送信する。なお、規格化画像は、線量換算投影データから再構成した画像に限らない。例えば、被写体断層像のCT値を非負値化したものを用いてもよいし、被写体断層像の各画素をCT値に応じた規格化係数に変換した画像を用いてもよい。
【００３７】
差分画像作成部37は、画像更新部36から更新画像を取得して基準画像と更新画像の差分をとることにより差分画像を作成する。作成された差分画像は、画像更新部36に送られる。
【００３８】
画像更新部36は、差分画像作成部37から取得した差分画像と規格化画像作成部35から取得した規格化画像とを用いて、逐次近似法に基づく更新式(後述する式(4))を用いて画像更新処理を行う。
【００３９】
反復処理部30は、基準画像作成部34により作成された基準画像及び規格化画像作成部35により作成された規格化画像を用いて反復処理を所定回数実行する。反復処理とは、差分画像作成部37における差分画像の作成と画像更新部36における更新画像の作成とを所定の更新回数だけ繰り返す処理である。所定の更新回数に満たない場合は、画像更新部36は更新画像を差分画像作成部37に送信し、差分画像の作成及び上述の画像更新処理を所定回数反復して実行する。所定の更新回数になった場合は、画像更新を停止する。画像更新部36は、更新処理により得られた画像をノイズ低減画像として画像表示部38に送信する。
【００４０】
画像表示部38は、画像更新部36により更新された画像であるノイズ低減画像を表示装置211に表示する。
【００４１】
本発明の逐次近似画像再構成処理は、従来の逐次近似画像再構成処理(SPS法)と比較して、基準画像と更新画像との差分をとる差分画像作成部37の構成を有する点で異なる。従来は繰り返し行われる画像更新の都度、更新画像を順投影し、順投影データを元の被写体投影データと差分しこれを逆投影することで差分画像を得ていた。つまり、順投影処理と逆投影処理とを繰り返し行っていた。しかし本発明では、演算時間を要する順投影処理と逆投影処理を、更新画像と基準画像との差分という、比較的演算量の少ない処理に置き換えた。このため、差分画像を作成するための処理時間が大幅に削減されることとなる。
【００４２】
次に図4を参照しながら、本発明のX線CT装置1の処理全体の流れを説明する。
【００４３】
まずX線CT装置1は、被写体3に対して位置決め撮影を行う。位置決め撮影では、ガントリ100を回転させずX線の照射方向を固定し、寝台101を所定の速さで移動しながら被写体3及び寝台101を透過したX線量を計測する。X線CT装置1は、位置決め撮影によって得られた透過X線データに基づいて位置決め画像を作成する。
【００４４】
中央制御装置200は、位置決め画像を使用して撮影条件や再構成条件等の各種条件設定を受け付ける。そして、中央制御装置200は設定された各種条件に基づいて本撮影を実行する。本撮影では、ガントリ100を回転させて被写体3の周囲の各方向からX線を照射し、被写体3及び寝台101を透過したX線情報を計測する。この本撮影により、被写体投影データを取得する(ステップS101)。
【００４５】
演算装置202は、取得した被写体投影データを用いて逐次近似画像再構成処理を行う(ステップS102)。逐次近似画像再構成処理の詳細な手順については、後述する。
【００４６】
演算装置202は、ステップS102の逐次近似画像再構成処理によって作成されたノイズ低減画像を表示装置211に表示する(ステップS103)。
【００４７】
次に、上述のステップS102の逐次近似画像再構成処理の手順について図5を参照して説明する。
【００４８】
図5に示すように、演算装置202の再構成処理装置221は、ステップS101の本撮影により得た被写体投影データを用いて、画像更新処理における基準画像として被写体断層像を再構成する(ステップS201)。作成された被写体断層像(基準画像)はRAM等に保持され、ステップS203の画像更新処理において使用される。基準画像の作成方法については、後述する(図6等参照)。
【００４９】
次に再構成処理装置221は、規格化画像を作成する(ステップS202)。具体的には、再構成処理装置221は、記憶装置213に記憶されている空気投影データを用いて、ステップS101で取得した被写体投影データを線量換算投影データに変換する。再構成処理装置221は、線量換算投影データを再構成した画像を規格化画像とする。空気投影データは、予め被写体及び寝台101がガントリ100内にない状態で収集した透過X線情報から作成してもよいし、ガントリ100の幾何学系を考慮して計算により作成してもよい。線量換算投影データは、上述したように、空気投影データと被写体投影データとを加算し、逆対数変換係数を乗じ、指数をとることで得られる。
【００５０】
再構成処理装置221は、ステップS201で作成した基準画像と、ステップS202で作成した規格化画像を用いて、画像更新を反復して行う(ステップS203)。ステップS203では、後述する式(4)に示す演算を所定の回数繰り返し行う(反復処理)。
【００５１】
すなわち、逆投影処理と順投影処理を含まない更新式(式(4))を用いて、所定回数画像を繰り返し更新する。画像更新処理の詳細は、後述する。
【００５２】
ステップS201〜ステップS203の処理により、多くの処理時間を要する逆投影及び順投影処理を省き、画像の差分処理(基準画像と更新画像との差分処理)に置き換えて画像を更新することが可能となる。これにより処理時間を短縮しつつ、ノイズが低減された更新画像(ノイズ低減画像)を得ることが可能となる。
【００５３】
以上の処理手順が本発明の画像再構成処理の手順である。以下の各実施形態では、基準画像や規格化画像についての具体的な態様を説明する。
【００５４】
［第1の実施の形態］
図6及び図7を参照して、本発明の第1の実施の形態を詳細に説明する。
【００５５】
本発明の第1の実施の形態のX線CT装置1は、画像更新処理の基準画像を更新せず、常に同じものを用いる点が特徴である。例えば、フィルタ逆投影法で再構成した画像を基準画像として用いる。
【００５６】
図6は、フィルタ逆投影法により基準画像を作成する場合の基準画像作成処理(A)の手順を説明するフローチャートである。
【００５７】
X線CT装置1は、X線発生装置102から扇形のX線(ファンビーム)を被写体3に対して照射して透過X線情報を収集し、被写体投影データを作成する。このように扇形に照射されたX線により得た投影データを、以下の説明では扇形照射投影データと呼ぶ。再構成処理装置221は、基準画像を作成するにあたり、まず扇形照射投影データを平行照射投影データに変換する投影データ幾何学変換処理を行う(ステップS301)。
【００５８】
投影データ幾何学変換処理では、複数角度(ビュー)の扇形照射投影データを用い、平行照射投影データに変換する。平行照射投影データとは、X線が被写体3に対して平行に照射されたと仮定した場合に得られる投影データである。ステップS301の処理により、後述の逆投影処理(ステップS303)における座標計算の演算量を低減できる。
【００５９】
次に再構成処理装置221は、平行照射投影データをフィルタ処理する(ステップS302)。フィルタ処理では、再構成処理装置221はまず平行照射投影データをフーリエ変換し、得られた周波数空間平行照射投影データに対して周波数フィルタを乗じる。更に逆フーリエ変換してフィルタ平行照射投影データを作成する。
【００６０】
フーリエ変換と逆フーリエ変換には、高速フーリエ変換処理(FFT)と高速逆フーリエ変換処理(IFFT)を用いることが望ましい。これにより投影データのフィルタ処理を高速に実行できる。また、上述の周波数フィルタとしては、例えば正規化ランプフィルタ等を用いることが好適であるが、これに限定されずその他のフィルタを用いてもよい。
【００６１】
再構成処理装置221は、ステップS302で得られたフィルタ平行照射投影データについて逆投影処理を行う(ステップS303)。すなわち、再構成処理装置221は、スキャナ10の回転角度(ビュー)毎に、X線発生装置102の焦点と対象画素の中心を通る直線のX線検出装置103との交点に位置する検出素子の座標を算出して、算出した検出素子位置に対応するフィルタ平行照射投影データの投影値を対象画素に加算する。全画素について、同様の処理を繰り返すと、被写体断層像が作成される。この被写体断層像を、逐次近似画像再構成処理において用いる基準画像とする。
【００６２】
図7は、図5のステップS202の規格化画像作成処理の詳細な手順を説明するフローチャートである。
【００６３】
図7に示すように、まず再構成処理装置221は、投影データ線量換算処理を行う(ステップS401)。投影データ線量換算処理とは、対数変換後のデータである被写体投影データを、対数変換前のデータ(線量換算投影データ)に換算する処理である。
【００６４】
具体的な処理手順としては、再構成処理装置221は、データ収集装置106から得られた扇形照射投影データ(被写体投影データ)と記憶装置213から取得した空気投影データとを加算し、所定の変換係数を乗じて指数をとる。これにより被写体投影データが扇形照射線量換算投影データに変換される。扇形照射線量換算投影データは、ファンビームによって得られ被写体3を透過した透過X線量を示すデータである。
【００６５】
次に、再構成処理装置221は、ステップS401で得られた扇形照射線量換算データに対してステップS301の投影データ幾何学変換処理と同様の処理を行い、平行照射線量換算投影データに変換する(ステップS402)。
【００６６】
再構成処理装置221は、フィルタ平行照射線量投影データの代わりにステップS402で得た平行照射線量換算投影データを用いて、ステップS303と同様の逆投影処理を行う。これにより規格化画像が得られる(ステップS403)。
【００６７】
第1の実施の形態では再構成処理装置221は、図6の処理手順(フィルタ逆投影法)で作成した基準画像と、図7の処理手順で作成した規格化画像とを用いて画像の更新処理(図5のステップS203)を行う。
【００６８】
以下、更新処理の詳細について説明する。
【００６９】
一般的な逐次近似再構成処理(SPS法)では、式(1)を用いて画像更新を行っている。
【００７０】
一方、本発明の第1の実施の形態では、以下の式(2)に示す近似を用いる。
2016132880.tif 000004
【００７１】
つまり、再構成処理装置221は、式(1)の更新式に含まれる行列演算部を以下の式(3)に示す近似を用いて画像の差分演算に置き換える。
2016132880.tif 000005
【００７２】
これにより、更新式は以下の式(4)となる。
2016132880.tif 000006
【００７３】
ここで、μj(0)は基準画像である。
【００７４】
再構成処理装置221は、行列演算部を除去した更新式である式(4)を用いて画像の更新処理を行う。
【００７５】
以上説明したように、従来の逐次近似画像再構成処理(SPS法)では、更新画像を順投影処理して作成する順投影データと、被写体投影データの差分投影データを逆投影処理する演算を含むが、本発明では、被写体投影データの順投影処理結果として基準画像を用いる(式(2)の近似を用いる)。これにより、演算負荷の高い行列演算である順投影処理及び逆投影処理を省き、更新式を演算負荷の低い画像処理(基準画像と更新画像の差分)に置き換えて構成することができる。その結果、逐次近似画像再構成処理の演算量は大幅に減り、処理時間を短縮することが可能となる。
【００７６】
［第2の実施の形態］
次に図8を参照して、本発明の第2の実施の形態を詳細に説明する。
【００７７】
第2の実施の形態では、画像の更新処理の基準画像として、逐次近似投影データ補正処理によって補正した投影データにより作成した画像を用いる点が第1の実施の形態と異なる。
【００７８】
図8に、第2の実施の形態における基準画像作成処理(B)について説明する。
【００７９】
再構成処理装置221は、被写体投影データを取得して、逐次近似投影データ補正処理を施す(ステップS501)。逐次近似投影データ補正処理において再構成処理装置221は、被写体投影データに畳み込みフィルタ(コンボリューションフィルタ)を所定回数繰り返し課す。畳み込みフィルタの積分核(カーネル)には平滑化効果のある加重平均フィルタ等を用いる。
【００８０】
次に再構成処理装置221は、ステップS501の逐次近似投影データ補正処理によって補正された投影データ(被写体補正投影データ)を用いてフィルタ逆投影処理を行うことにより画像を再構成する(ステップS502)。フィルタ逆投影処理は、第1の実施の形態の基準画像作成処理(A)と同様である。すなわち、再構成処理装置221は、被写体補正投影データについて投影データ幾何学変換処理、投影データフィルタ処理、及び逆投影処理を行い、被写体断層像を得る。この被写体断層像を第1補正断層像とする。
【００８１】
次に再構成処理装置221は、ステップS502の処理により得た第1補正断層像に鮮鋭化フィルタを課し、鮮鋭化された被写体断層像を第2補正断層像とする(ステップS503)。鮮鋭化フィルタを課すことで、逐次近似投影データ補正処理で過度に平滑化された画素の画素値を回復させる。鮮鋭化フィルタには、例えばラプラシアンフィルタから構成する加重平均フィルタ等を用いることができる。
【００８２】
再構成処理装置221は、図8の基準画像作成処理(B)により作成した第2補正断層像を図5に示す逐次近似画像再構成処理における基準画像として使用する。また、規格化画像は、第1の実施の形態と同様の処理(図7参照)で作成したものを用いればよい。
【００８３】
再構成処理装置221は、上述の第2補正断層像と規格化画像を用いて、式(4)に示す更新式により画像の更新処理を行う。画像の更新を所定回数繰り返すことにより、ノイズ低減画像を作成する。
【００８４】
以上説明したように、第2の実施の形態では、従来の逐次近似画像再構成処理(SPS法)における被写体投影データの順投影処理結果として、逐次近似投影データ補正処理を施した断層像(第2補正断層像)を基準画像として用いる。これにより、第1の実施の形態の効果に加え、ストリークアーチファクトを低減して画質向上したノイズ低減画像を提供できる。
【００８５】
［第3の実施の形態］
次に図9を参照して、第3の実施の形態を詳細に説明する。
【００８６】
第3の実施の形態では、X線CT装置1は、逐次近似画像再構成処理における更新処理(図5のステップS203)において、画像更新の前に画像拡大処理を行い、画像更新の後に画像縮小処理を行う。なお、基準画像及び規格化画像は第1または第2の実施の形態と同様の手法で作成されたものを使用すればよい。
【００８７】
図9に、第3の実施の形態における画像更新処理(A)の手順を示す。
【００８８】
まず再構成処理装置221は、画像の拡大処理を行う(ステップS601)。
【００８９】
この拡大処理において再構成処理装置221は、寝台101対し垂直、かつガントリ100に対し平行な断面(アキシャル面)において基準画像の画像サイズを拡大した拡大基準画像を作成する。例えば、画像サイズを2倍に拡大すると、画素サイズは元の画像の2分の1になる。再構成処理装置221は、規格化画像についても同様にサイズを拡大して拡大規格化画像を作成する。
【００９０】
次に、再構成処理装置221は、拡大基準画像と拡大規格化画像とを用いて、第1の実施の形態と同様に式(4)に示した更新処理を繰り返し行う。これにより拡大ノイズ低減画像が作成される(ステップS602)。
【００９１】
再構成処理装置221は、アキシャル面において拡大ノイズ低減画像の画像サイズを縮小し、元の基準画像と同じ画像サイズにする(ステップS603)。
【００９２】
以上説明したように、第3の実施の形態では、拡大した画像を用いて更新処理を行うため、第1の実施の形態の効果に加え、ノイズ低減画像中の小さな物体の形状を維持することができる。これにより、画質向上したノイズ低減画像を作成することが可能となる。
【００９３】
［第4の実施の形態］
次に、本発明の第4の実施の形態を説明する。
【００９４】
第4の実施の形態の逐次近似画像再構成処理では、更新式(4)の基準画像μj(0)として合成画像を用いる。合成画像は、第1の実施の形態の手法(フィルタ逆投影処理)で作成した基準画像(以下、第1基準画像と呼ぶ)と更新画像とを所定の比率で合成した画像とする。
【００９５】
再構成処理装置221は、まず基準画像及び規格化画像を第1の実施の形態と同様に作成する。
【００９６】
図5のステップS203の更新処理において、再構成処理装置221は、式(4)に示す更新式の基準画像μj(0)を以下の式(5)のように合成画像に置き換える。
2016132880.tif 000007
【００９７】
ここで、νj(0)は、n回目の更新時の画素番号jに対する合成係数(重み)である。この合成係数は、操作者が予め指定した任意の値としてもよいし、撮影条件等に基づいて決定してもよい。基準画像と合成する更新画像は1回前の更新画像を用いる。
【００９８】
以上のように、第4の実施の形態の逐次近似画像再構成処理では、基準画像と更新画像とを所定の合成係数をかけて合成した画像である合成画像を、更新式(4)における基準画像であるμj(0)に置換して更新処理を行う。これにより、更新式(4)における差分画像「μj(0)−μj(n)」が大きくならないようにできる。
【００９９】
第4の実施の形態の処理により、第1の実施の形態の効果に加え、ノイズ低減効果を高めることができ、画質向上したノイズ低減画像を提供できる。
【０１００】
［第5の実施の形態］
次に図10を参照して、第5の実施の形態を詳細に説明する。
【０１０１】
第5の実施の形態の逐次近似画像再構成処理では、規格化画像としてノイズ低減した規格化画像を用いる。
【０１０２】
再構成処理装置221は、まず基準画像及び規格化画像を第1または第2の実施の形態と同様に作成する。
【０１０３】
画像更新処理では、再構成処理装置221は、図10に示す手順の画像更新処理(B)を行う。
【０１０４】
すなわち、再構成処理装置221は、例えば図7の手順で作成した規格化画像を、以下の式(6)に示すように所定の回数だけ更新する(ステップS701)。この処理の結果、ノイズが低減された規格化画像を得る。以下、ステップS701の処理により得た画像をノイズ低減規格化画像と呼ぶ。
2016132880.tif 000008
【０１０５】
式(6)は、上述の更新式(4)において、以下の式(7)による近似を適用したものである。
2016132880.tif 000009
【０１０６】
再構成処理装置221は、基準画像とノイズ低減規格化画像を用いて、第1の実施の形態の更新処理(式(4))と同様に、画像更新を行う。これにより、ノイズ低減画像を得る。
【０１０７】
以上により、第5の実施の形態では、規格化画像としてノイズ低減した規格化画像を使用し、規格化画像を更新処理してノイズを低減する。これにより、第1の実施の形態の効果に加え、ノイズ低減効果を高めることができ、画質向上したノイズ低減画像を提供できる。
【０１０８】
［第6の実施の形態］
次に、図11〜図12を参照して、第6の実施の形態を詳細に説明する。
【０１０９】
第6の実施の形態では、逐次近似画像再構成処理において、X線CT装置1は、被写体断層像上で操作者が指定した関心領域の内外で異なる規格化係数を設定可能とする。規格化画像は、画像更新時に画素毎に画像のノイズ低減強度を調整するために用いる画像であるため、操作者が指定した関心領域のノイズ低減強度を操作者の操作によって自在に設定可能とする。
【０１１０】
再構成処理装置221は、まず基準画像を第1または第2の実施の形態と同様に作成する。
【０１１１】
第6の実施の形態では、規格化画像作成処理において、操作者による関心領域の指定を受け付ける。そのため、再構成処理装置221は表示装置211に基準画像を表示する。
【０１１２】
図11は、関心領域を指定するための表示画面51の一例である。表示画面51には、被写体断層像(基準画像)61が表示される。
【０１１３】
図11に示すように表示された被写体断層像61に対して、操作者がマウス等の入力装置212を使用して、図12に示すように関心領域60を指定する。再構成処理装置221は、操作者が指定した関心領域60の内側と外側との境界の座標情報をもとに規格化画像の各画素の規格化係数を設定する。
【０１１４】
このとき、再構成処理装置221は、関心領域60内の規格化係数を大きく設定し、関心領域60外の規格化係数が小さくなるように設定することが望ましい。規格化係数を大きくすると、その部位のノイズ低減強度を大きくすることができる。また、関心領域60の内と外の境界付近については、滑らかに連続するように規格化係数を設定することが望ましい。また、関心領域60からの距離に応じて規格化係数を滑らかに変化させてもよい。
【０１１５】
なお、関心領域60は一箇所に限定されず、複数設定してもよい。
【０１１６】
以上のように第6の実施の形態によれば、関心領域の内外で異なる規格化係数を設定することが可能となる。これにより、例えば、関心領域のノイズ低減の度合いを大きくするなど、ノイズ低減の強度を自在に設定可能となる。
【０１１７】
［第7の実施の形態］
次に、第7の実施の形態を詳細に説明する。
【０１１８】
第7の実施の形態の逐次近似画像再構成処理では、組織毎にCT値が異なることを考慮して、組織に応じた規格化画像を用いる。また、規格化画像として、被写体断層像を以下の手順で補正処理した画像を用いるものとする。
【０１１９】
逐次近似画像再構成処理において、まず、再構成処理装置221は、第1または第2の実施の形態と同様に基準画像(被写体断層像)を作成する。
【０１２０】
次に、再構成処理装置221は規格化画像を作成する。第7の実施の形態では、上述したように規格化画像として、被写体断層像を補正処理したものを用いる。以下、被写体断層像の補正について説明する。
【０１２１】
被写体断層像の画素値はCT値であるため負数を含むことがある。負数が含まれていると規格化画像として適さないので、再構成処理装置221は、被写体断層像の画素値の非負値化処理を行い、負数を含まない非負値化画像を作成する。画素値の最小値が負数の場合、各画素値から最小負数値を差し引くことで非負値化処理を行ってもよい。第7の実施の形態では、この非負値化画像を規格化画像として用いる。被写体断層像が負数を含んでいない場合は、被写体断層像をそのまま規格化画像として使用することもできる。
【０１２２】
また、被写体断層像の補正処理の別の例としては、CT値と規格化係数との相関を設定した相関関数80を用いた補正等が挙げられる。
【０１２３】
すなわち、再構成処理装置221は、図13に示すようなCT値と規格化係数との相関を示す関数(相関関数80)を使用して、被写体断層像の各画素を画素値に応じた規格化係数に変換する。このようにCT値が規格化係数に変換された画像を、更新処理における規格化画像とする。
【０１２４】
なお、この場合においても、規格化係数が非負値になるような相関関数80を使用する。
【０１２５】
このとき、ノイズ低減効果を相対的に高めたい組織に該当するCT値の規格化係数が大きくなるような相関関数80を用いることが望ましい。図13の例は、軟部組織に該当するCT値の規格化係数が大きくなるような相関関数80を示している。
【０１２６】
再構成処理装置221は、上述したように被写体断層像を非負値化した画像、または相関関数80により変換した画像を、画像更新処理における規格化画像として使用する。再構成処理装置221は、作成した規格化画像を表示装置211に表示し、操作者に示すようにしてもよい。
【０１２７】
更新処理は、基準画像と、上述のように作成した基準画像とを用いて、第1の実施の形態と同様に行う。
【０１２８】
以上説明したように、第7の実施の形態では、被写体断層像のCT値情報を基に規格化係数を設定した画像を規格化画像として利用する。このため、第1の実施の形態の効果に加え、更に、組織に応じたノイズ低減効果を得ることが可能となる。これにより操作者が注目したい組織の画質が向上したノイズ低減画像を提供できる。
【０１２９】
以上、本発明の好適な実施形態について説明したが、本発明は、上述の実施形態に限定されるものではない。当業者であれば、本願で開示した技術的思想の範疇内において、各種の変更例または修正例に想到し得ることは明らかであり、それらについても当然に本発明の技術的範囲に属するものと了解される。
【産業上の利用可能性】
【０１３０】
本発明によれば、ノイズ低減効果を維持しつつ処理時間を短縮することが可能な演算装置、X線CT装置、及び画像再構成方法を提供できる。
【符号の説明】
【０１３１】
1  X線CT装置、3  被写体、10  スキャナ、20  操作ユニット、100  ガントリ、101  寝台、102  X線発生装置、103  X線検出装置、104  コリメータ、105  高電圧発生装置、106  データ収集装置、107  ガントリ駆動装置、200  中央制御装置、201  入出力装置、202  演算装置、211  表示装置、212  入力装置、213  記憶装置、221  再構成処理装置、222  画像処理装置、30  反復処理部、31  データ記憶部、32  透過線量収集部、33  投影データ作成部、34  基準画像作成部、35  規格化画像作成部、36  画像更新部、37  差分画像作成部、38  画像表示部、51  表示画面、60  関心領域、61  基準画像、80  相関関数
【図１】
2016132880.tif 000010
【図２】
2016132880.tif 000011
【図３】
2016132880.tif 000012
【図４】
2016132880.tif 000013
【図５】
2016132880.tif 000014
【図６】
2016132880.tif 000015
【図７】
2016132880.tif 000016
【図８】
2016132880.tif 000017
【図９】
2016132880.tif 000018
【図１０】
2016132880.tif 000019
【図１１】
2016132880.tif 000020
【図１２】
2016132880.tif 000021
【図１３】
2016132880.tif 000022
【手続補正書】
【提出日】平成29年7月7日(2017.7.7)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
被写体の周囲の各方向から照射され、前記被写体を透過したX線の情報に基づいて被写体投影データを作成する投影データ作成部と、
前記被写体投影データから基準画像を作成する基準画像作成部と、
更新処理におけるノイズ低減強度を調整する規格化係数を画素毎に定義した画像である規格化画像を作成する規格化画像作成部と、
前記基準画像と前記規格化画像を用いて、反復処理を所定回数実行する反復処理部と、を備え、
前記反復処理部は、
前記基準画像と前記更新処理により得た更新画像とを差分し、差分画像を作成する差分画像作成部と、
前記差分画像及び前記規格化画像を用いて前記更新処理を実行し、更新画像を作成する画像更新部と、
を有することを特徴とする演算装置。
【請求項２】
前記基準画像作成部は、フィルタ逆投影法により基準画像を作成することを特徴とする請求項1に記載の演算装置。
【請求項３】
前記基準画像作成部は、前記被写体投影データに対しノイズ低減処理を施して補正投影データを作成し、作成した補正投影データを画像再構成して第1補正断層像を作成し、作成した前記第1補正断層像に鮮鋭化処理を施した第2補正断層像を前記基準画像として作成することを特徴とする請求項1に記載の演算装置。
【請求項４】
前記反復処理部は、
前記差分画像作成部における基準画像及び前記画像更新部における規格化画像として、それぞれサイズを拡大した拡大断層像及び拡大規格化画像を使用して前記差分画像の作成及び前記更新処理を所定回数反復して実行することにより拡大更新画像を作成し、前記拡大更新画像を元のサイズに縮小した縮小画像を作成することを特徴とする請求項1に記載の演算装置。
【請求項５】
前記基準画像と前記更新画像とを所定の比率で合成した画像を、前記差分画像作成部における基準画像として使用することを特徴とする請求項1に記載の演算装置。
【請求項６】
前記規格化画像作成部により作成された規格化画像を前記画像更新部及び前記差分画像作成部における基準画像として前記更新処理を実行することにより更新規格化画像を作成し、作成した前記更新規格化画像を前記画像更新部における規格化画像として使用することを特徴とする請求項1に記載の演算装置。
【請求項７】
前記基準画像上に関心領域を設定する関心領域設定部と、
前記関心領域の内外で異なる規格化係数を設定する設定部と、を備え、
前記規格化画像作成部は、設定された前記規格化係数を画素毎に定義した規格化画像を作成することを特徴とする請求項1に記載の演算装置。
【請求項８】
前記規格化画像作成部は、前記被写体投影データと空気投影データとを加算して逆対数変換係数を乗じて指数をとることにより線量換算データを求め、前記線量換算データを画像再構成することにより、前記規格化画像を作成することを特徴とする請求項1に記載の演算装置。
【請求項９】
前記規格化画像作成部は、前記被写体の断層像のCT値を非負値化した非負値化画像を前記規格化画像とすることを特徴とする請求項1に記載の演算装置。
【請求項１０】
前記規格化画像作成部は、CT値と規格化係数との相関を設定した相関関数に基づいて前記被写体の断層像から規格化画像を算出することを特徴とする請求項1に記載の演算装置。
【請求項１１】
請求項1に記載の演算装置を備えたX線CT装置。
【請求項１２】
被写体の周囲の各方向から照射され、前記被写体を透過したX線の情報に基づいて被写体投影データを作成する投影データ作成ステップと、
前記被写体投影データから基準画像を作成する基準画像作成ステップと、
更新処理におけるノイズ低減強度を調整する規格化係数を画素毎に定義した画像である規格化画像を作成する規格化画像作成ステップと、
前記基準画像と前記規格化画像を用いて、反復処理を所定回数実行する反復処理ステップと、を含み、
前記反復処理ステップにおいて、
前記基準画像と前記更新処理により得た更新画像とを差分し、差分画像を作成する差分画像作成ステップと、
前記差分画像及び前記規格化画像を用いて前記更新処理を実行し、更新画像を作成する画像更新ステップと、
を含むことを特徴とする画像再構成方法。
【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】０００５
【補正方法】変更
【補正の内容】
【０００５】
【特許文献１】国際公開第2013/008702号
【手続補正３】
【補正対象書類名】明細書
【補正対象項目名】００２７
【補正方法】変更
【補正の内容】
【００２７】
本発明において、再構成処理装置221は、後述する逐次近似画像再構成処理を行う。
これにより、ノイズが低減されたCT画像を高速に再構成する。本発明の演算装置202(再構成処理装置221)が実行する逐次近似画像再構成処理の詳細については後述する。再構成処理装置221は、生成したCT画像を記憶装置213に保存するとともに表示装置211に表示する。
【手続補正４】
【補正対象書類名】図面
【補正対象項目名】図７
【補正方法】変更
【補正の内容】
【図７】
2016132880.tif 000023
【手続補正書】
【提出日】平成29年8月28日(2017.8.28)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
被写体の周囲の各方向から照射され、前記被写体を透過したX線の情報に基づいて被写体投影データを作成する投影データ作成部と、
前記被写体投影データから基準画像を作成する基準画像作成部と、
更新処理におけるノイズ低減強度を調整する規格化係数を画素毎に定義した画像である規格化画像を作成する規格化画像作成部と、
前記基準画像と前記規格化画像を用いて、反復処理を所定回数実行する反復処理部と、を備え、
前記反復処理部は、
前記基準画像と前記更新処理により得た更新画像とを差分し、差分画像を作成する差分画像作成部と、
前記差分画像及び前記規格化画像を用いて前記更新処理を実行し、更新画像を作成する画像更新部と、
を有することを特徴とする演算装置。
【請求項２】
前記基準画像作成部は、フィルタ逆投影法により基準画像を作成することを特徴とする請求項1に記載の演算装置。
【請求項３】
前記基準画像作成部は、前記被写体投影データに対しノイズ低減処理を施して補正投影データを作成し、作成した前記補正投影データを画像再構成して第1補正断層像を作成し、作成した前記第1補正断層像に鮮鋭化処理を施した第2補正断層像を前記基準画像として作成することを特徴とする請求項1に記載の演算装置。
【請求項４】
前記反復処理部は、
前記差分画像作成部における基準画像及び前記画像更新部における規格化画像として、それぞれサイズを拡大した拡大断層像及び拡大規格化画像を使用して前記差分画像の作成及び前記更新処理を所定回数反復して実行することにより拡大更新画像を作成し、前記拡大更新画像を元のサイズに縮小した縮小画像を作成することを特徴とする請求項1に記載の演算装置。
【請求項５】
前記基準画像と前記更新画像とを所定の比率で合成した画像を、前記差分画像作成部における基準画像として使用することを特徴とする請求項1に記載の演算装置。
【請求項６】
前記規格化画像作成部により作成された規格化画像を前記画像更新部及び前記差分画像作成部における基準画像として前記更新処理を実行することにより更新規格化画像を作成し、作成した前記更新規格化画像を前記画像更新部における規格化画像として使用することを特徴とする請求項1に記載の演算装置。
【請求項７】
前記基準画像上に関心領域を設定する関心領域設定部と、
前記関心領域の内外で異なる規格化係数を設定する設定部と、をさらに備え、
前記規格化画像作成部は、前記規格化係数を画素毎に定義した規格化画像を作成することを特徴とする請求項1に記載の演算装置。
【請求項８】
前記規格化画像作成部は、前記被写体投影データと空気投影データとを加算して逆対数変換係数を乗じて指数をとることにより線量換算データを求め、前記線量換算データを画像再構成することにより、前記規格化画像を作成することを特徴とする請求項1に記載の演算装置。
【請求項９】
前記規格化画像作成部は、前記被写体の断層像のCT値を非負値化した非負値化画像を前記規格化画像とすることを特徴とする請求項1に記載の演算装置。
【請求項１０】
前記規格化画像作成部は、CT値と規格化係数との相関を設定した相関関数に基づいて前記被写体の断層像から規格化画像を算出することを特徴とする請求項1に記載の演算装置。
【請求項１１】
請求項1に記載の演算装置を備えたX線CT装置。
【請求項１２】
被写体の周囲の各方向から照射され、前記被写体を透過したX線の情報に基づいて被写体投影データを作成する投影データ作成ステップと、
前記被写体投影データから基準画像を作成する基準画像作成ステップと、
更新処理におけるノイズ低減強度を調整する規格化係数を画素毎に定義した画像である規格化画像を作成する規格化画像作成ステップと、
前記基準画像と前記規格化画像を用いて、反復処理を所定回数実行する反復処理ステップと、を含み、
前記反復処理ステップにおいて、
前記基準画像と前記更新処理により得た更新画像とを差分し、差分画像を作成する差分画像作成ステップと、
前記差分画像及び前記規格化画像を用いて前記更新処理を実行し、更新画像を作成する画像更新ステップと、
を含むことを特徴とする画像再構成方法。
【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】００２４
【補正方法】変更
【補正の内容】
【００２４】
データ収集装置106は、X線検出装置103の個々のX線検出素子により検出されるX線量の情報を収集し、デジタル信号に変換し、透過X線情報として操作ユニット20の演算装置202に順次出力する。
【手続補正３】
【補正対象書類名】図面
【補正対象項目名】図７
【補正方法】変更
【補正の内容】
【図７】
2016132880.tif 000024
【国際調査報告】
2016132880.tif 000025
2016132880.tif 000026