JP4731571B2

JP4731571B2 - 反復再構成用の打切り補償アルゴリズム

Info

Publication number: JP4731571B2
Application number: JP2007546239A
Authority: JP
Inventors: ジンハンイェ; マリーケイダービン
Original assignee: Koninklijke Philips NV; Koninklijke Philips Electronics NV
Current assignee: Koninklijke Philips NV
Priority date: 2004-12-17
Filing date: 2005-12-05
Publication date: 2011-07-27
Anticipated expiration: 2025-12-05
Also published as: CN101080651B; WO2006064404A2; DE602005020844D1; US8013307B2; US20090310746A1; EP1828808B1; WO2006064404A3; ATE465425T1; JP2008524575A; CN101080651A; EP1828808A2

Description

本発明は、診断画像化システム及び方法に関する。本発明は、減衰補償をもつ単光子放出コンピュータ断層撮影（ＳＰＥＣＴ）システムと関連した特別な用途を見出し、特にこのようなシステムに関連して記載されるだろう。本発明は、ポジトロン放出断層法システム（ＰＥＴ）及びコンピュータ断層撮影システム（ＣＴ）等の他の画像化システムにも適用可能であると理解されるだろう。

核医療画像化は、患者を画像化するために放射能源を採用する。典型的に放射性薬品が、患者に注入される。放射性薬品化合物は、予測可能な割合及び特性エネルギーでガンマ線崩壊を経る放射性同位元素を含む。１以上の放射検出器が、放出された放射線を監視及び記録するために患者に隣接して位置される。時々、該検出器は、複数の方向からの放出された放射線を監視するために、患者の周りを回転又は角度を割り出される。検出された位置及びエネルギーのような情報に基づいて、体内の放射性薬品分布は決定され、該分布の画像は、循環系、選択された器官又は組織への放射性薬品の取り込み等を研究するために再構成される。

典型的に、反復再構成技術では、再構成されたデータの新たなボリュームの見積もりが生成される毎に、画像データの前に再構成されたボリュームが、検出器の面に順投影（forward project）される。順投影されたデータは、実際の投影データと比較される。もし再構成された画像が完全だったならば、これら２つの投影データは、一致し、違いはないだろう。しかしながら、画像が作られるとき、典型的に違い又は誤りが存在する。誤り又はその逆数が、ボリュームの画像を訂正するために画像ボリュームに逆投影（back project）される。

これらの技術はうまくいくが、打ち切り誤り（truncation errors）を生じる傾向にある。すなわち、対象が全ての検出器の位置で検出器により完全に見られていない場合、ボクセルのいくつかを再構成するデータは、視野のいくつかのみに現れる。したがって、興味のある主要な部位のデータは、典型的に完全にサンプリングされるが、取り囲む組織は、低密度でサンプリングされるか、又は、視野のわずかな部分のみでサンプリングされる。各反復サイクルでは、逆投影された誤り値は、完全にサンプリングされた又はサンプリング不足の領域を含め、全体の画像ボリュームに及ぶ。検出器に比較のためのデータが存在しなかった比較からの誤り値の逆投影を許可することは、対象が低密度でサンプリングされた領域でクリッピングするような、誤り及びアーチファクトとなり得る。

本発明は、上で参照された問題及び他のものを克服する、新規で改善された画像化装置及び方法を提供する。

本発明の一態様にしたがって、画像化システムが開示される。少なくとも１つの放射線検出器が、被験体からの放出及び透過放射線の少なくとも１つを検出及び測定するために、被験体受入開口部に隣接して配置され、該検出器は、放射線を受け取るために被験体の周りで可動であり、複数の投影角で測定された投影データを生成する。視野（ＦＯＶ）手段は、各投影角での放射線検出器の視野に属する複数のピクセルを決定する。画像プロセッサは、前記視野内のみで検出された放射線を、画像表示に反復して再構成する。画像表示は、画像メモリに反復的に再構成される。

本発明の他の態様によると、画像化の方法が開示される。複数の投影角での被験体からの放出及び透過放射線の少なくとも１つが、検出及び測定される。測定された投影データは、複数の投影角で生成される。視野にある複数のピクセルは、放射線が各投影角で放射線検出器により検出される範囲内で決定される。該視野内で検出された放射線は、画像メモリ内の画像表示に反復的に再構成される。

本発明の１つの利点は、削減されたアーチファクトにある。

他の利点は、再構成された画像のクリッピングの削減にある。

他の利点は、再構成された画像のぼやけの削減にある。

本発明の他の利点及び利益は、以下の好ましい実施例の詳細な説明を読み且つ理解することで、当業者には明らかになるだろう。

本発明は、様々なコンポーネント及びコンポーネントの配置並びに、様々なステップ及びステップの構成を取り得る。図面は、好ましい実施例を描写する目的のみのためにあり、本発明を制限するものとして解釈されるべきではない。

図１を参照すると、核画像化システム１０は、典型的に、回転ガントリ１４をサポートする固定ガントリ１２を含む。１以上の検出ヘッド１６は、興味ある領域又は検査の領域１８から放出される放射線を検出するために回転ガントリ１４により担持される。各検出ヘッドは、例えば光電子増倍管、フォトダイオード等のような、シンチレータ及び光感受性要素のような検出要素又は検出器２０の二次元アレイを含む。ＣＺＴ素子のようなＸ線／電気直接変換器も考えられる。各ヘッド１６は、各放射線応答を、検出器の面上の位置（ｘ，ｙ）及びエネルギー（ｚ）のデジタル信号表示に変換する回路２２を含む。検出器２０上の事象の位置は、名目上ｘ及びｙ座標と呼ぶ二次元（２Ｄ）直交座標系で分解及び／又は決定される。しかしながら、他の座標系も考慮される。コリメータ２４は、検出器２０の各要素が放射線を受けることができる方向及び角度の広がりを制御する、すなわち検出器２０は既知の線に沿ってのみ放射線を受けることができる。したがって、放射線が検出される検出器２０の決定された位置及びカメラ１６の角度の位置は、各放射線の事象が起こる名目的な線を規定する。

典型的に、画像化される対象は、１以上の放射性薬品又は放射性同位元素を注入され、コーチ２６によりサポートされた検査領域１８に位置される。このような同位元素の数少ない例は、Ｔｃ−９９ｍ，Ｇａ−６７，及びＩｎ−１１１である。対象内部の放射性薬品の存在は、対象から放出放射線を生成する。放射線は、１以上の選択された投影方向で投影放出データを収集するために、検査領域１８の周りで角度を割り出されるか、又は回転することができる検出ヘッド１６により検出される。投影放出データ、例えば検査領域１８（例えば角度位置レゾルバ（resolver）２８から得られる）の周りの各検出ヘッド１６の位置（ｘ，ｙ），エネルギー（ｚ），角度位置（θ）が、放出データメモリ３０に格納される。

放出データは、患者の人体の吸収特性の変化により生じる誤りを含むので、しばしば透過放射線源が、放出データを訂正する付加的な減衰情報を提供するために用いられる。一実施例では、１以上の放射線源５０が、検出ヘッド１６から検査領域１８を横切って据付けられる。オプションとして、放射線源５０は、放射線源５０からの透過放射線が、放出データを補完するために、検査領域１８の反対側の対応する検出ヘッド１６に向かって該ヘッドにより受けられるように、検出ヘッド１６間に又は検出ヘッド１６に据付けられる。好ましくは、コリメータ２４は、透過放射線を視準合わせする。例えば該コリメータ２４は、検出器２０に（平行なビーム構造に対して）垂直な線、又は、放射線源と検出器との間の他の直接の経路に沿って進行しない、透過放射線の部分を受けることからシンチレータ２０を制限する。代替として、他の視準形態が採用される及び／又は該視準は、前記放射線源で起こりうる。

好ましくは、放射線源５０は、対応する各検出ヘッド１６の軸の長さをそれぞれ伸ばす線源である。好ましくは、該線源は、放射性核種で満たされ、端で密閉される細い鉄の管の形態を取る。代替として、放射線源５０は、棒状源、点源、平面長方形状源、円盤源、フラッド（flood）源、放射性核種で満たされた管若しくは容器、又はＸ線管のようなアクティブ放射線生成器でもあり得る。

引き続き図１を参照し、更に図３，４を参照すると、透過スキャンは、透過放射線源５０からの透過放射線も検出ヘッド１６により受けられ、透過投影データが生成されるように実施される。典型的に放出投影データは、異なるエネルギーのとき以外、同時に収集される。透過及び放出データは、対象５２の大きさが典型的に検出器２０の視野よりも大きいため、しばしば打ち切られ、結果として打ち切られた投影Ｐ１，Ｐ２となる。検出器２０の視野に含まれない対象５２の打ち切られた部分５４は、切り取られる。以下でより詳細に議論されるように、打切り部分５４は、打切り部分５４が打ち切られない他の角度で得られた投影から補償される。

分類器６０は、放出投影データ及び透過投影データを、その相対的なエネルギーに基づいて分類する。放出データは放出データメモリ３０に格納され、透過データは、透過メモリ６２に格納される。ＦＯＶ決定手段６４は、以下で詳細に議論されるであろう従来技術のよく知られた方法の１つにより、どのデータが対応する検出器２０の視野に収集されるか決定する。

引き続き図１を参照し、更に図２を参照すると、データプロセッサ７０は、ＦＯＶ透過データから３Ｄ透過放射線画像又は減衰マップ７４を反復的に再構成し、一方、画像プロセッサ７２は、ＦＯＶ放出データから３Ｄ放出画像を反復的に再構成する。

好ましくは、データプロセッサ７０は、反復的にベイジアン反復透過勾配アルゴリズム（ＢＩＴＧＡ）を実行し、一方、画像プロセッサ７２は、最尤推定期待値最大化アルゴリズム（ＭＬＥＭ）を実行する。再構成プロセスの最初の反復の準備では、減衰マップメモリ７４及び画像メモリ７６は、メモリ７４，７６を仮定してロードすることにより初期化されるか、又は減衰マップ及び放出マップの第１見積もりをする。減衰マップの第１見積もりは、オプションとして被験体を含む既定の輪郭の内側では一様な減衰値であり、外側ではゼロであることを特徴とする。同様に第１放出マップ見積もりは、オプションとして、輪郭の内側では一様な値であり、外側ではゼロであることを特徴とする。代替として、付加的な先験的な情報の利用可能性は、第１見積もりをより正確にすることを可能にする。

第１順投影機８０は、透過減衰マップ７４から投影透過データを生成する。第１比較手段８２は、訂正要素を決定するために、測定された透過データを順投影されたデータと比較する。ＦＯＶ決定手段６４は、被験体のどの部分かが従来技術の方法の１つで知られた視野よりも広いかどうかを決定する。例えばＦＯＶ決定手段６４は、いくつかのピクセルがゼロにされるかどうかを決定する。好ましくは、一連の基準スキャンが生成され、そこから視野が決定される。代替としてＦＯＶ決定手段６４は、隣接したピクセル間の値の激しい変化を決定するため、正規化された透過投影データの各列を端から中心まで検索する。もし現在のピクセルと前のピクセルとの間の違いが、事前に特定された閾値よりも大きいならば、ＦＯＶ決定手段６４は、打切りが存在すると結論付ける。フィルタ手段８４は、視野データの外をフィルタリングすることにより、打切りの決定点を越える全てのピクセルを訂正要素マトリクスから除外する。例えば、再構成データの次の反復は、ＦＯＶの外からの値を誤って訂正されない。視野に対応する訂正要素は、第１訂正メモリ８６に格納される。

引き続き図１及び２を参照し、更に図５を参照すると、第１角度位置α１で投影する間、第１打切り部分９０は、投影Ｐ４により検出されないが、該第１部分９０は、第２角度位置α２で取られた投影Ｐ３により検出される。同様に投影が第２角度位置α２で取られる間、第２打切り部分９２は、投影Ｐ３により検出されず、第２打切り部分９２は、第１角度位置α１で取られる投影Ｐ４により検出される。部分９４は、説明された投影Ｐ３及びＰ４の両方で打ち切られる。打切りデータ補償手段１００は、視野ＦＯＶの外の打ち切られた部分９０，９２を補償するために、異なる角度位置で取られた打ち切られていない投影のデータを供給する。第１逆投影機１０２は、訂正要素を減衰マップメモリ７４に逆投影する。（非ゼロ値の）逆投影は、実際のデータＰ１，Ｐ２が収集された領域に限定される。減衰マップ７４の要素は、減衰要素メモリ１０４に格納される各ボクセルの減衰要素を表す。

引き続き図１及び２を参照すると、画像プロセッサ７２は、画像メモリ７６内に３Ｄ画像表示を反復的に再構成し、各反復は、順投影操作及び逆投影操作を含む。第２順投影機１１０は、画像メモリ７６に格納された再構成された放出画像から投影放出データを生成する。より詳細には、放出データが沿って受信される各線に対して、第２順投影機１１０は、減衰要素メモリ１０４に格納された透過減衰要素を通じて対応する線の投影を計算する。減衰訂正手段１１２は、同じ線に沿う減衰要素投影を補償するために、受信された放出データを調整する。第２比較手段１１４は、線に沿う減衰訂正測定放出データと、同じ線に沿う投影されたデータとを訂正要素を決定するために比較する。フィルタ手段８４'は、打ち切り決定点を越える全てのピクセルを、視野データの外をフィルタリングすることにより再構成から除外する。視野に対応する訂正要素は、第２訂正メモリ１１６に格納される。第２バックプロジェクタ１１８は、訂正要素を再構成された画像に逆投影する、すなわち再構成された画像が訂正要素で更新される。打切りデータ補償手段１００は、視野ＦＯＶの外側の打ち切られた部分９０，９２を補償するために、異なる角度位置で得られた打ち切られていない投影のデータを供給する。視野の外側の領域のいかなる決定された訂正も捨てられる。

一実施例では、打切りデータ補償手段１００は、再構成された画像内の打ち切られた側の境界線を滑らかにするために、先験的な知識をアルゴリズムに適用する。再構成された画像内のいくつかのピクセルは、もし投影が非常に打ち切られるならば、わずかな投影角度のみで検出器ＦＯＶを「視る」。この場合、再構成された画像の打ち切られた側に、いくらか不明瞭さがつくられる。先験的知識をアルゴリズムに適用することは、視野の外側になることを確かめられる再構成された対象のピクセルをゼロにすることにより、打ち切られた側のあいまいさを軽減することを手助けする。

もちろん、データプロセッサ７０及び画像プロセッサ７２が、様々な他の勾配及び非勾配反復アルゴリズム並びに透過及び放出反復アルゴリズムの両方を含む他の反復アルゴリズムを実行することも考慮される。各連続した反復は、最も新たに更新された放出画像で実行される。再構成された放出データ及び測定された放出データが収束する、すなわち訂正が事前に選択されたレベルより下に落ち込むとき、反復再構成プロセスは終了する。

一実施例では、母集団の先験的な画像（すなわち「ナレッジセット」）が収集される。該先験的な画像は、現在検査中の患者からではなく、ＣＴスキャン、ＰＥＴ透過スキャン、他のＳＰＥＣＴ透過スキャン等から得られる他の被験体の画像のような、類似の構造の横断面であり得る。他のオプションとして初期の先験的な画像は、同じ被験体を異なるモダリティで得る。

ビデオプロセッサ１３０は、スライス、投影、３Ｄレンダリング及び他の画像情報を３Ｄ画像メモリ７６から取ってきて、ビデオモニタ１３２、プリンタ、記憶メディア等のような１以上の人が視られるディスプレイにディスプレイするための画像表示を適切にフォーマットする。もしビデオプロセッサが選択された画像情報を再構成時に繰り返し取ってくるならば、ディスプレイは、再構成された画像が最終的な画像に収束するように各反復でよりはっきりとなるだろう。

引き続き図１及び２を参照すると、反復透過勾配アルゴリズムは、以下のように表すことができる。

ただし、
μ^ｎ _ｊは、ｊ番目のピクセルのｎ番目の反復での減衰係数であり、
ｆ_ｉは、ｉ番目の検出器の基準スキャン値であり
ｌ_ｉ，ｊは、ｊ番目のピクセル内でｉ番目の検出器から延びる線のセグメントの長さであり、
Ｋ_ｉ，ｊは、ｊ番目の検出器からｉ番目の検出器までの一群のピクセルであり、
Ｙ_ｉは、ｉ番目の検出器での透過カウントである。

事前空間変化ガンマを含むベイジアン反復透過勾配アルゴリズム（ＢＩＴＧＡ）は、次のように表すことができる。

ただし、ｒ_ｊは、ピクセル値に「好ましい」値であり、ａ_ｊは、その重みである。

好ましくは、ベイジアン反復透過勾配アルゴリズムは、各ピクセルの値を水例えば組織又は空気例えば肺の中の空気の減衰係数に収束するように高めるために、先験的な知識を使用する修正された前のブロックを含む。

ただし、σは、事前の範囲を規定する。

典型的に検出器マトリクスは、実際の検出器の視野よりも広い領域をカバーする。検出マトリクス内の視野の外の打ち切られた部分５４に位置されたピクセルは、ゼロにされ、上で議論されたようにＦＯＶ決定手段６４により切り取られる。打ち切りの決定点を越える全てのピクセルは、再構成に含まれず、例えば投影データは、再構成から誤ったゼロの値を排除するために打ち切られる。

視野からのみピクセルを再構成する修正されたベイジアン反復透過勾配アルゴリズムは、次のように表すことができる。

検出器ＦＯＶの外側でどの事象も欠落しないほど十分に対象が小さいならば、もちろん投影データは、打ち切られない。

本発明は、好適な実施例に関して記載される。明らかに、修正及び変形は、上記詳細な説明を読み、理解することで他人にもなされるだろう。本発明が、請求項又はそれと等価のものの範囲にある限り、このような修正及び変更を全て含むとして解釈されるべきであるということが、意図される。

図１は、画像化システムの図式的な説明である。図２は、該画像化システムの一部の詳細な図式的な説明である。図３は、平行ボア視準での対象の１つの側の打ち切りを図式的に示す。図４は、拡大又はファンビーム視準での対象の一つの側の打切りを図式的に示す。図５は、検出器の第１角度位置の視野の外側であるが、検出器の第２角度位置の視野内にある打切り部分を図式的に示す。

Claims

被験体からの放出及び透過放射線の少なくとも１つを検出及び測定するために、被験体受入開口に隣接して位置された少なくとも１つの放射線検出器であって、前記放射線を受けるために、前記被験体の周りを動くことができ、複数の投影角度で測定された投影データを生成する、放射線検出器と、
各投影角度での前記放射線検出器の視野に属する複数のピクセルを決定するＦＯＶ手段と、
前記視野内のみで検出された前記放射線を画像表示に反復的に再構成する画像プロセッサと、
前記画像表示が反復して再構成される画像メモリと、
を有する、画像化システム。
前記画像プロセッサが、
対応する前記画像メモリから前記画像表示を投影する順投影機と、
前記投影された画像表示と前記測定された投影データとを比較し、該比較に基づいて一群の訂正要素を決定する手段と、
前記決定された視野内の前記ピクセルのみに対応する前記訂正要素を前記画像表示に逆投影する逆投影機と、
を含む、請求項１に記載のシステム。
前記画像プロセッサが、
訂正要素を逆投影する前に、視野データの外を除去するフィルタ手段
を更に含む、請求項２に記載のシステム。
前記逆投影機が、前記視野ピクセルの外側を逆投影しない、請求項３に記載のシステム。
前記打ち切られたデータが打ち切られていないような異なる角度で得られる前記投影から、打ち切られていないデータを供給することと、
先験的な知識を使用することと、
の少なくとも１つにより、打ち切られたデータを補償する手段
を更に含む、請求項１に記載のシステム。
前記画像プロセッサが、
反復透過再構成勾配アルゴリズムと、
反復透過再構成非勾配アルゴリズムと、
反復放出再構成勾配アルゴリズムと、
反復放出再構成非勾配アルゴリズムと、
の少なくとも１つを実行する、請求項１に記載のシステム。
前記画像プロセッサが、反復透過勾配アルゴリズムを実行し、前記視野の外のピクセルが、

ただし、
μ^ｎ _ｊは、ｊ番目のピクセルのｎ番目の反復での減衰係数であり、
ｆ_ｉは、ｉ番目の検出器の基準スキャン値であり、
ｌ_ｉ，ｊは、前記ｊ番目のピクセル内で前記ｉ番目の検出器から延びる線のセグメントの長さであり、
Ｋ_ｉ，ｊは、ｊ番目の検出器から前記ｉ番目の検出器までの一群のピクセルであり、
Ｙ_ｉは、ｉ番目の検出器での透過カウントである
前記再構成から除外される、請求項６に記載のシステム。
前記視野の外のピクセルを再確認し、該再確認された視野の外のピクセルを前記再構成画像においてゼロにするために、補償手段が先験的な知識を用いる、打ち切られたデータを補償する手段
を更に含む、請求項７に記載のシステム。
前記被験体からの放出投影を検出するＰＥＴスキャナ及びＳＰＥＣＴスキャナの少なくとも１つ
を更に含む、請求項１に記載のシステム。
複数の投影角度における被験体からの放出及び透過放射線の少なくとも１つを検出及び測定するステップと、
前記複数の投影角度で測定された投影データを生成するステップと、
前記放射線が、各投影角度で放射線検出器により検出される視野に属する複数のピクセルを決定するステップと、
前記視野で検出された前記放射線を、画像メモリの画像表示に反復して再構成するステップと、
を有する、画像化する方法。
前記再構成するステップが、
前記画像表示を投影するステップと、
前記投影された画像表示と前記測定された投影データとを比較するステップと、
該比較に基づいて、一群の訂正要素を決定するステップと、
前記決定された視野内のピクセルのみに対応する前記訂正要素を、前記画像表示に逆投影するステップと、
を含む、請求項１０に記載の方法。
前記再構成するステップが、
前記視野の外のデータを逆投影の前に除去するステップ
を更に含む、請求項１１に記載の方法。
前記複数のピクセルを決定するステップが、
基準スキャンを生成するステップ
を含む請求項１０に記載の方法。
他の投影からの打ち切られていないデータを供給することと、
先験的な知識を用いることと、
の少なくとも１つにより、打ち切られたデータを補償するステップ
を更に含む、請求項１０に記載の方法。
前記再構成するステップが、
反復透過再構成勾配アルゴリズムと、
反復透過再構成非勾配アルゴリズムと、
反復放出再構成勾配アルゴリズムと、
反復放出再構成非勾配アルゴリズムと、
の少なくとも１つを実行するステップを含む、請求項１０に記載の方法。
前記再構成するステップが、反復透過勾配アルゴリズムを実行するステップを含み、前記視野の外のピクセルが、

ただし、
μ^ｎ _ｊは、ｊ番目のピクセルのｎ番目の反復での減衰係数であり、
ｆ_ｉは、ｉ番目の検出器の基準スキャン値であり、
ｌ_ｉ，ｊは、ｊ番目のピクセル内でｉ番目の検出器から延びる線のセグメントの長さであり、
Ｋ_ｉ，ｊは、ｊ番目の検出器からｉ番目の検出器までの一群のピクセルであり、
Ｙ_ｉは、ｉ番目の検出器での透過カウントである、
再構成から排除される、請求項１５に記載の方法。
前記視野の外のピクセルを先験的な知識を用いることにより再確認するステップと、
前記再確認された視野の外のピクセルを、前記再構成された画像内でゼロにするステップと、
を更に含む、請求項１６に記載の方法。
請求項１１の方法を実行する、少なくとも１つのＰＥＴスキャナ及びＳＰＥＣＴスキャナ。
投影が放射線検出器により検出される被験体からの測定された放出投影を格納する放出メモリと、
データが前記放射線検出器により検出される前記被験体からの測定された透過データを格納する透過メモリと、
前記放射線検出器の視野を決定するプロセッサと、
透過減衰マップを格納する減衰マップメモリと、
放出マップを格納する画像メモリと、
（ａ）前記減衰マップメモリから前記減衰マップを投影することを介して、前記減衰マップを生成するステップと、
（ｂ）前記投影された減衰マップと前記透過データとを比較するステップと、
（ｃ）該比較に基づいて、一群の訂正要素を決定するステップと、
（ｄ）前記決定された視野内の前記データのみに対応する前記訂正要素を、前記減衰マップメモリの前記減衰マップに逆投影するステップと、
（ｅ）前記投影された減衰マップが前記透過データに収束するまで、ステップ（ａ）から（ｄ）を反復して繰り返すステップと、
を実行するようにプログラムされるプロセッサと、
を含む、診断画像化システム。
（ｆ）前記画像メモリから前記放出マップを投影することを介して、前記放出マップを生成するステップと、
（ｇ）前記投影された放出データを、前記減衰マップから計算された前記減衰訂正要素で訂正するステップと、
（ｈ）前記減衰マップから計算された前記減衰訂正要素で、前記測定された放出投影を訂正するステップと、
（ｉ）前記投影された放出データを前記測定された放出データと比較するステップと、
（ｊ）該比較に基づいて、一群の訂正要素を決定するステップと、
（ｋ）前記決定された視野内の前記データのみに対応する前記訂正要素を、前記画像メモリの前記放出マップに逆投影するステップと、
（ｌ）前記減衰訂正投影された放出データが、前記測定された放出データに収束するまで、ステップ（ｆ）から（ｋ）を反復して繰り返すステップと、
を実行するようにプログラムされたプロセッサと、
前記放出マップから再構成された画像表示をレンダリングする、ヒト可読ディスプレイと、
を更に含む、請求項１９に記載のシステム。