JPH075263A - コンプトン・フリー・ガンマ線映像 - Google Patents

コンプトン・フリー・ガンマ線映像

Info

Publication number
JPH075263A
JPH075263A JP3139186A JP13918691A JPH075263A JP H075263 A JPH075263 A JP H075263A JP 3139186 A JP3139186 A JP 3139186A JP 13918691 A JP13918691 A JP 13918691A JP H075263 A JPH075263 A JP H075263A
Authority
JP
Japan
Prior art keywords
energy
photons
compton
scattered
counts
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP3139186A
Other languages
English (en)
Other versions
JP3123024B2 (ja
Inventor
Gideon Berlad
ベルラド ジデオ
Dov Maor
マオル ドブ
Yigal Shrem
シュレム イガル
Adrian Soil
ソイル アドリアン
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Elscint Ltd
Original Assignee
Elscint Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Elscint Ltd filed Critical Elscint Ltd
Publication of JPH075263A publication Critical patent/JPH075263A/ja
Application granted granted Critical
Publication of JP3123024B2 publication Critical patent/JP3123024B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/16—Measuring radiation intensity
    • G01T1/161—Applications in the field of nuclear medicine, e.g. in vivo counting
    • G01T1/164—Scintigraphy
    • G01T1/1641—Static instruments for imaging the distribution of radioactivity in one or two dimensions using one or several scintillating elements; Radio-isotope cameras
    • G01T1/1647—Processing of scintigraphic data
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/16—Measuring radiation intensity
    • G01T1/161—Applications in the field of nuclear medicine, e.g. in vivo counting
    • G01T1/164—Scintigraphy
    • G01T1/1641—Static instruments for imaging the distribution of radioactivity in one or two dimensions using one or several scintillating elements; Radio-isotope cameras
    • G01T1/1642—Static instruments for imaging the distribution of radioactivity in one or two dimensions using one or several scintillating elements; Radio-isotope cameras using a scintillation crystal and position sensing photodetector arrays, e.g. ANGER cameras
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/16—Measuring radiation intensity
    • G01T1/17—Circuit arrangements not adapted to a particular type of detector
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/42—Arrangements for detecting radiation specially adapted for radiation diagnosis
    • A61B6/4208—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector
    • A61B6/4258—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector for detecting non x-ray radiation, e.g. gamma radiation

Landscapes

  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Molecular Biology (AREA)
  • High Energy & Nuclear Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Biomedical Technology (AREA)
  • Optics & Photonics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Medical Informatics (AREA)
  • General Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Nuclear Medicine (AREA)
  • Measurement Of Radiation (AREA)

Abstract

(57)【要約】 【目的】散乱された光子が最終映像に及ぼす影響を大幅
に削減し、取得後数秒以内に事実上コンプトン・フリー
である映像を得ることによって、記録される映像を定性
的にも定量的にも改善すること。 【構成】位置当たり、また測定されたエネルギー当たり
について記憶された検出カウントから位置についてのエ
ネルギー・スペクトルを作る。次に、システム依存性コ
ンプトン・エネルギー分布と位置を提供するためシステ
ムのエネルギー拡散関数によるコンプトン散乱光子のエ
ネルギー分布の繰り込みによって、非散乱光子が作る位
置当たりのカウントを決定する。

Description

【発明の詳細な説明】
【０００１】
【産業上の利用分野】この発明はガンマ線カメラ映像シ
ステムに関し、さらに特定すれば、コンプトン散乱イベ
ントの分量が非常に少ない映像、すなわち、事実上コン
プトン・フリーである映像を得るための方法とシステム
に関する。
【０００２】
【従来の技術】人体を貫通する場合、ガンマ線光子は、
コンプトン効果によって散乱する特定の確率をもってい
る。そのような散乱は光子の運動の方向と光子のエネル
ギーを変更させる。散乱された１個の光子がガンマ線カ
メラ映像中に記録されると、散乱された光子から偽りの
位置情報が誘導される。
【０００３】原則として、散乱された光子は無視しなけ
ればならない。ただし、散乱された光子を無視するため
の効率的で有効な基準に到達することは容易でない。例
えば、散乱の過程において光子は自己が有するエネルギ
ーの一部を喪失するが、典型的なガンマ線カメラのエネ
ルギー解像度は、散乱された光子と散乱されない光子と
の間の記録されたエネルギーに大量の重複を許容する解
像度であるため、エネルギー・レベル基準は有効でな
い。
【０００４】
【発明が解決しょうとする課題】この発明の目的は、散
乱された光子が最終映像に及ぼす影響を大幅に削減し、
取得後数秒以内に事実上コンプトン・フリーである映像
を得ることによって、記録された映像を定性的にも定量
的にも改善することである。やはり、散乱された光子が
最終映像に及ぼす影響を減少させ、または、減少させる
ことを意図した先行技術の方法が多数ある。ただし、既
知のすべての先行技術の方法とシステムのうちには、散
乱された光子が最終映像に及ぼす影響が大幅に削減され
ており、最終映像に、定性的にも定量的にも、事実上、
コンプトン散乱の汚染がないことが示されている、高品
質であって実際的な方法は知られていない。
【０００５】先行技術の一部を示す：
【０００６】（ａ）デュアル・ウインドウ取得 ピーク・エネルギーの周辺の標準全エネルギー・ウイン
ドウを使用するデータ取得に加えて、第２映像が『低エ
ネルギー』で取得される。その後、第２映像には実験定
数が乗算され、散乱のイベントによって最終映像上にも
たらされた効果を控除するため、第１映像から減算され
る。
【０００７】（ｂ）アクセルソンらの方法（論文『コン
プトン散乱光子輻射コンピュータ断層撮影法の減算』
Ｊ．Ｎｕｃ．Ｍｅｄ．２５，ｐｐ．４９０−４９４，１
９８４参照）
【０００８】この方法は空間deconvolution 方法であ
り、コンプトン散乱によって生じた対数テールを点拡散
関数から抽出する。対数テールの大きさが映像全体を通
しての定数であると考えられる。
【０００９】この方法は米国特許第４−７８０，８２３
号に記載されている。基本的に、この特許において説明
されているのは、映像の全区域上のフライ上に適用され
る定空間フィルタである。この方法はコンプトン散乱イ
ベントをまったく削減しないが、散乱によってぼんやり
になった縁を鋭利にする傾向がある。
【００１０】コーラルらの方法 この方法は米国特許第４，８３９，８０８号に示されて
いる。この特許において、エネルギー・スペクトルは、
コンプトン散乱スペクトルを代表すると仮定される第３
オーダーの多項式と組み合わされた光電ピーク先鋭カー
ブを伴った各与えられた場所に適合される。その適合の
結果は、各点におけるイベントを散乱部分と非散乱部分
に分離すると仮定される。相当な統計を得るため、ま
た、計算時間を節約するため、実際には、エネルギー・
スペクトルが粗い空間光子上で解析される。その後、さ
らに精細な映像光子に関する値を得るため、コンプトン
散乱分数の値は内挿される。
【００１１】前記の方法全部が有する基本的欠点は、一
般に、コンプトン散乱イベントを説明するためにエネル
ギー・スペクトルの本質的な物理的性質を利用していな
いことである。その結果、前記のコンプトン散乱補正法
は、すべて、本質的に、コンプトン散乱イベントを正確
に計算し、減算することができない。
【００１２】前記の方法（ａ）と（ｂ）は、共にまった
く経験的であって非局所的である。その結果、どちらの
方法もコンプトン散乱に起因する人工物の克服に成功し
ておらず、その結果、実際には使用されていない。
【００１３】加重取得モジュール方法は、コンプトン散
乱イベントを非散乱イベントから分離しようと試みさえ
しない。むしろ、その方法は、散乱によってぼんやりに
された縁をいくらか鋭利にする。この方法は、定性的映
像においては限定的に使用されるが、散乱補正映像の主
な用途であると期待されている定量的ＳＰＥＣＴ映像に
おいては事実上まったく無用である。
【００１４】コーラルらの方法は、いくつかの短所を抱
えている。主な短所は、コンプトン散乱分布を適合させ
るために使用された第３オーダー多項式がコンプトン散
乱イベントの物理的分布を真に定義していないという事
実である。第３オーダー多項式は、特に、光電ピークを
重複させるスペクトルの高いエネルギー・コンプトン散
乱しきい値において限定されていない。重複の高エネル
ギー部分は最も関連がある部分であり、その部分におい
てエネルギー・ウインドウの従来の方法でコンプトン散
乱をなくすことは困難である。
【００１５】映像あたり数千回もの適合を行う場合に、
過誤の多い映像量と人工物をもたらす多数の不正確なコ
ンポジット・カーブの使用を防止することは、事実上不
可能である。
【００１６】
【課題を解決するための手段】この発明は、米国特許第
４，８３９，８０８号に記載されたコーラルらのシステ
ムにいくらか似ている。この発明は、エネルギー・スペ
クトルを局所的に解析し、エネルギー形状が既知であっ
て大きさが未知である光電ピーク成分と、形状も大きさ
も未知であるコンプトン散乱成分とで構成された『試
行』関数と適合させようと試みる。ただし、コーラルら
が使用したようなコンプトン成分エネルギー・スペクト
ルに適合させるために任意の多項式を使用する代りに、
その後においてコンプトン成分エネルギー・スペクトル
を作るために使用されるコンプトン多重散乱関数を誘導
するためにコンプトン過程の真の物理的特性が使用され
る。従って、この適合は、光電ピーク成分とコンプトン
成分の両方についての定量情報を提供する。
【００１７】非常に詳細に見ると、未知のもの、すなわ
ち、光電ピーク成分の大きさ、コンプトン多重散乱成分
の形状と大きさを決定するため、下記の『入力』が使用
される。 （１）ピクセル当たりの測定されたエネルギー・スペク
トル。この中には散乱された光子と非散乱光子に起因す
るカウントが含まれる。 （２）測定されたシステム・エネルギーの拡散関数。同
位元素中心線のためには、この拡散関数は光電ピーク・
エネルギー形状である。
【００１８】試行関数のコンプトン成分の形状は、下記
の段階を使用して誘導される：
【００１９】１）仁科−ケルビン方程式。この方程式
は、光子と電子の物理的量子散乱が、与えられたエネル
ギーからの単独の相互作用で、さらに低いエネルギーに
散乱させるための光子に関する確率分布に変換されるこ
とを説明する。
【００２０】２）繰り込みの反復によって、高いオーダ
ーの相互作用に関する確率分布が作られる。
【００２１】３）まだ未知である値の局所的係数によっ
て、様々なオーダーの重ね合わせが作られる。
【００２２】４）物理的なコンプトン多重散乱成分を作
るため、この重ね合わせが測定されたシステム・エネル
ギー散乱関数に繰り込まれる。
【００２３】この試行関数を、取得されデータの測定済
みエネルギー・スペクトルに局所的に適合させることに
より、多重散乱コンプトン係数の値と光電ピークのマグ
ニチュードさが得られる。このことによって、取得され
たデータからコンプトン汚染を除去することができるよ
うになる。
【００２４】多重散乱確率のセットの標準直交化によっ
て、この局所的適合の実行を促進させることが望まし
い。この標準直交化は、例えば、グラハム−シュミット
手順によって行うことができる。実際には、コンプトン
多重散乱確率の係数を解く代りに、標準直交化された関
数の係数が決定されるべき未知数である。
【００２５】この発明の一代替システムは、やはり、コ
ーラルらが使用した方法よりさらに正確な方法で補償す
るため、計算に関する統計を改善する問題を扱う。この
発明のその代替システムにおいては、『準局所的』区域
上における、すなわち、四角形（２ｎ＋１）×（２ｍ＋
１）［ｎとｍは１，２，３．．．］のピクセル内におけ
るコンプトン散乱データを集計するに際して、コンプト
ン分布の円滑さが利用される。
【００２６】四角形のピクセルから得られた値は、前記
四角形内のピクセルの数によって除算され、前記四角形
の中心ピクセルに割り当てられる。この計算が、各ピク
セルについて反復される。この『準局所的』方法はコー
ラルらが使用した内挿法よりも正確である。さらに、こ
の方法の計算時間はコーラルらが内挿法で要した計算時
間とほぼ同じである。望ましい実際のアルゴリズムは、
処理時間がその四角形の中のピクセルの数とは事実上無
関係であるようにする、連続的に最新化される四角形の
アルゴリズムである。
【００２７】局所的または準局所的に測定されたエネル
ギー・スペクトルへの『試行』関数の実際の適合は、多
数の方法で行うことができる。望ましい実施例において
は、最小自乗適合または最大公算方法のいずれかが使用
される。最大公算方法は最上自乗適合より多くの時間を
消費するが、最大公算方法は、ガンマ線処理におけるポ
アソンの統計を明確に考慮する。
【００２８】この発明の望ましい局面にしたがって、ガ
ンマ線映像システムによって作られた映像に及ぼされる
コンプトン散乱光子の影響を取り除く方法が提供され
る。前記方法は下記の段階で構成される：
【００２９】源から来るガンマ線輻射に応答してガンマ
線検出器に衝突する光子をカウントとして検出する段
階；
【００３０】検出器上における光子の衝突の位置にした
がって、検出された各カウントについてのＸ，Ｙ座標位
置を決定する段階；
【００３１】衝突する光子のエネルギーを測定する段
階；
【００３２】共にエネルギー範囲上にまたがっており、
各ウインドウがエネルギー範囲の一部上で延伸している
複数のエネルギー・ウインドウを使用する段階；
【００３３】検出された各カウントを、衝突する光子の
測定されたエネルギーにしたがって前記各エネルギー・
ウインドウに割り当てられたメモリ・マトリックス中の
位置に置く段階；前記メモリ・マトリックスは、それぞ
れ、Ｘ，Ｙ座標で構成される単独のエネルギー・ウイン
ドウに対応するデジタル映像である。
【００３４】各Ｘ，Ｙ座標位置について各メモリ・マト
リックス中に集積されたカウントを使用することによっ
て、Ｘ，Ｙ位置当たりの測定されたエネルギー・スペク
トルを決定する段階；前記測定されたエネルギー・スペ
クトルには、散乱された光子と散乱されない光子の両方
が作ったカウントが含まれる。
【００３５】コンプトン多重散乱成分を使用して、コン
プトン多重散乱光子のエネルギー分布を計算する段階；
この場合、『多重散乱』には少なくとも第１のオーダー
の散乱が含まれる。
【００３６】コンプトン多重散乱光子のシステム依存性
エネルギー分布を得るため、多重散乱光子の計算された
エネルギー分布をガンマ線映像システムのエネルギー拡
散関数に繰り込ませる段階；
【００３７】下記によって試行関数を作る段階；コンプ
トン多重散乱光子のシステム依存性エネルギー分布を各
オーダーの散乱について１個の未知の係数と合計するこ
と、および、非散乱光子の影響を係数として光電ピーク
・エネルギー形状の時間に加算すること
【００３８】未知のＸ，Ｙ係数を解決し、したがって、
非散乱光子のＸ，Ｙ位置についてのカウントを決定する
ため、作られた試行関数を、Ｘ，Ｙ座標位置について測
定されたエネルギー・スペクトルに適合させること；お
よび
【００３９】コンプトン・フリーである映像を提供する
ため、ピクセルについての非散乱光子のカウントを使用
すること。
【００４０】この発明の特長に従うと、散乱光子のエネ
ルギー分布の計算には、第１オーダーの散乱光子のエネ
ルギー分布を決定するための仁科−ケルビン方程式の第
１オーダー確率分布への変換が含まれ、もっと高いオー
ダーの散乱光子分布が考慮される場合は、もっと高いオ
ーダーの散乱光子分布を得るための、第１オーダー確率
分布の繰り込みが含まれる。
【００４１】この発明の付加された特長にしたがうと、
各エネルギー・ウインドウについてのエネルギー分布を
平均するか統合することによって、エネルギー分布から
個別の関数のセットが得られる。
【００４２】この発明のまた別の付加された特長は、計
算数を減少させるための、個別関数の直交関数セットへ
の変換である。
【００４３】この発明のまた別の特長に従うと、エネル
ギー・スペクトルについての個別関数の取得は、『大き
な』ピクセルを使用して準局所的に達成される。比較的
円滑なコンプトン分布について、この『大きな』ピクセ
ルを使用することは、正規のピクセルそれぞれにおける
統計カウントの低さを補償する。
【００４４】この発明のまた別の特長に従うと、適合は
最小自乗適合を使用して行われる。代替方法として、最
大公算解決法が使用される。
【００４５】この発明のまた別の特長は、試行関数に使
用されたすべてのエネルギー形状の起源がデータから独
立していることである。したがって、コンプトン・フリ
ーである映像を作るデータ処理に先立って、そのような
形状を記憶または計算することができる。
【００４６】
【一般的説明】図１の１１は、コンプトン・フリーであ
る映像を作るためのこの発明のガンマ線映像システム１
１をブロック図の形態で示す。図１は、エネルギー・ス
ペクトル測定段階１２、試行関数準備段階１４と、コン
プトン・フリーである映像１６を提供するカーブ適合ま
たは計算段階１５から成っている。
【００４７】エネルギー・スペクトル測定段階１２は、
ガンマ線輻射検出器１７から成っている。ガンマ線輻射
検出器１７は、１８において示されているような、ガン
マ線輻射検出器１７の面に衝突する光子に応答して電気
信号を提供する。検出器の面に光子が衝突すると、導体
１９，２１，２２上に電気信号が提供される。導体１
９，２１，２２は直接、座標コンピュータ２３に向けら
れており、座標コンピュータ２３はガンマ線輻射検出器
１７上における光子１８の衝突のＸ，Ｙ位置を決定す
る。
【００４８】導体２２，２４は、光子のエネルギーを示
す電気信号を運ぶ。エネルギーを示す電気信号はエネル
ギー（Ｚ）補正回路２５に与えられる。エネルギー処理
回路２６は、検出されたエネルギーの範囲をシステム・
オペレータがあらかじめ定めたエネルギー・ウインドウ
の数に分割する。
【００４９】エネルギーが特定の限度以内である場合、
エネルギー補正回路は、座標コンピュータ２３がＸ，Ｙ
位置を決定することができるようにする許可信号を導体
２７上に送る。この情報は、映像補正器とデジタイザ回
路３１に向けられる。映像補正器とデジタイザ回路３１
は、各イベントのＸ，Ｙ座標を補正し、デジタル化す
る。多数のイベントに関する情報は、光子のエネルギー
に応じて複数のマトリックス３２中に置かれる。各マト
リックスはメモりであって、例えば、ウインドウ＃１の
場合は２２ｋｅｖから２５ｋｅｖまでのウインドウ、ウ
インドウ＃２の場合は２５ｋｅｖから２８ｋｅｖまでの
ウインドウというように、特定のエネルギー・ウインド
ウについて、Ｘ，Ｙ位置についてのイベントのカウント
を保持する。ウインドウは、Ｗ１，Ｗ２，Ｗ３からＷｎ
までとして示されている。この場合、ｎは所定のエネル
ギー・ウインドウの数である。
【００５０】このように、マトリックスは検出器上のイ
ベントの位置に対応するＸ，Ｙ位置に分割される。Ｘ，
Ｙ位置は、また、最終映像中のピクセルにも対応する。
映像予備処理装置３３は、各ウインドウからのピクセル
毎のデータを受信し、ブロック３４において示されてい
る通り、ピクセルについて測定されたか取得されたエネ
ルギー・スペクトルＮE を計算する。この取得されたエ
ネルギー・スペクトルには、非散乱光子と散乱光子の両
方のカウントが含まれる。
【００５１】図１の試行関数段階１４は、散乱イベント
と非散乱イベントを含む理論的すなわち試行エネルギー
分布ｎ（Ｘ，Ｙ，Ｅ）を作る。この場合：
【数１】
【００５２】上記の式において、 ε＝Ｅ／ｍ−ｃ2 ：電子静止エネルギーｍ−ｃ2 を単位
とする光子エネルギー Ｐ（ε’，ε）＝エネルギーε’におけるシステム・エ
ネルギー拡散関数。その形状は測定されており既知であ
る。 п（εo ，ε）＝原エネルギーεo からエネルギーεま
でｍ回散乱された光子によって作られたイベントのエネ
ルギー分布（すなわち、ｍ回散乱された光子のエネルギ
ー確率分布の形状） ＮP （ｘ，ｙ）＝散乱しなかった光子によって作られた
イベントの空間分布（カウント／ピクセル） Ｑm （ｘ，ｙ）＝ｍ回散乱した光子によって作られたイ
ベントの空間分布（カウント／ピクセル） εo ＝光子の原エネルギー ε＝測定された光子のエネルギー ε’＝光子の中間エネルギー
【００５３】この発明の重要な目的は、ＮP （ｘ，
ｙ）、すなわち、散乱しなかったイベントの空間分布を
決定することである。
【００５４】ピクセルについてのイベントのカウントを
決定するため、ブロック１５は、試行分布ｎ（ｘ，ｙ，
ε）の計算された値を測定された値ＮＥ（ｘ，ｙ）に適
合させる。非散乱空間分布の知識があれば、１６に示さ
れる通り、コンプトン・フリーである映像が作られる。
【００５５】図１の試行関数準備段階１４において発生
する計算の詳細が図２に示されている。さらに特定する
と、図２に示されている通り、図１と図２のブロック４
１中に示されたシステム・エネルギー拡散関数上の値
は、図１と図２のブロック３６中に入力される。システ
ムのエネルギー拡散関数は既知であると仮定される。シ
ステムのエネルギー拡散関数は、一度測定された後、シ
ステムのメモリ中に保持される。測定は、既知のエネル
ギーのガンマ線輻射の源の提供と、例えば、図１の機器
１１によるガンマ線輻射の検出によって容易に行われ
る。検出は、エネルギー源と検出器の間のコンプトン散
乱媒体なしで行われる。この検出によって、ブロック４
１に示されたように、検出器のエネルギー解像度に起因
する単一エネルギー源に関するエネルギー拡散関数が提
供される。
【００５６】準備ブロック３６は、πm 、すなわち、散
乱された光子のエネルギー分布を各散乱オーダーについ
て計算する。
【００５７】この計算は、拡散された非偏光光子の様々
なオーダーを誘導するための仁科−ケルビン方程式を使
用して行われる。すなわち、
【００５８】
【数２】
【００５９】さらに高いオーダーの散乱は、前記方程式
を使用し、反復繰り込みによって誘導される：
【数３】
【００６０】前記方程式が回帰的には解決されること、
すなわち、高いオーダーの方程式は、それぞれ、それよ
り低いオーダーの方程式の知識を必要とすることに注目
しなければならない。
【００６１】コンプトン散乱光子のエネルギー分布は、
各オーダーの散乱について、システムから独立したカー
ブを提供する。ただし、システムから独立したカーブ
は、Ｃ m （εo ，ε）によって示されたシステム依存性
コンプトン多重散乱エネルギー分布を提供するため、シ
ステム・エネルギー拡散関数の作用を受ける。ｃm （ε
o ，ε）分布の形状は、システム・エネルギー拡散関数
Ｐ（ε’，ε）でπm を繰り込むことによって得られ
る。すなわち、
【数４】
【００６２】この方程式のセットは、システム・エネル
ギー拡散関数による作用を受けた後の、各オーダーの拡
散についてのコンプトン・エネルギー分布の形状を提供
する。図２は、π1 およびそのπ2 ．．．πm につい
て、ボックス４２，４３，４４中で仁科−ケルビン方程
式を使用し、πm 値においてもたらされる計算を示す。
【００６３】ブロック４２，４３，４４中のπ1 ，π2
とπm の形状が、それぞれ、ブロック４６，４７，４８
の中で、ブロック４１のシステム・エネルギー拡散関数
に繰込まれ、それによって、形状Ｃ１，Ｃ２，Ｃ３を提
供していることが示されている。π1 ，π2 等とＣ１，
Ｃ２等を決定するための計算が回帰的であることが、π
1 ，π2 等からとＣ１，Ｃ２等から進む矢印によって示
されている。
【００６４】この発明によって、計算数を抜本的に減少
させる方法が提供される。この計算の減少は、セットＣ
m （εo ，ε）の標準直交化によって達成される。標準
直交化は、グラハム−シュミット手順を使用して直交関
数（ベクトル）を作ることによって提供される：
【００６５】
【数５】
【００６６】上記において、エネルギー上の合計（積
分）は下記の式で定義される：
【数６】
【００６７】アレー・セット｛Θm ｝は次式に従うこと
に注目しなければならない：
【数７】
【００６８】標準直交化はコンピュータ４９内で行わ
れ、結果、すなわち、Θ1 ，Θ2 ，Θ m は、例えば、ブ
ロック５１，５２，５３中に示される。
【００６９】Θk を使用して、コンプトン総和（ＥＱ
（４））を書き直すことができる：
【数８】
【００７０】上記の式において、
【数９】
【００７１】［標準直交ベース｛Θm ｝を使用すると、
どのようなベクトルＶでもΘm ’のアレーの重ね合わせ
によって示される： ｖ＝Σ＜ｖ．Θm ＞Θm （１０）］
【００７２】試行分布は、今や、下記になる： ｎ（Ｘ，Ｙ；Ｅ）＝Ｎｐ（Ｘ，Ｙ）Ｐ（Ｅo ，Ｅ）＋Ｃ（Ｘ，Ｙ；Ｅo ，Ｅ ） （１１）
【００７３】上記の式において、
【数１２】
【００７４】以降、既知のエネルギー拡散関数Ｐは、＜
Ｐ＞＝１のように標準直交化される。
【００７５】望ましい一つの実施例においては、最小自
乗適応が使用される。さらに特定すると、方程式（１）
の試行関数がｎ（Ｘ，Ｙ；Ｅ）であり、ブロック３４か
らの多重ウインドウ取得の結果がＮE である場合、非散
乱光子を原因とするカウント数ＮE （Ｘ，Ｙ）のための
解決であって、各ピクセルについての相違の自乗の総和
を最低にする解決δ（Ｘ，Ｙ）が求められる。 δ（Ｘ，Ｙ）＝＜［ｎ（Ｘ，Ｙ；Ｅ）−ＮE （Ｘ，Ｙ）］2 ＞ （１３）
【００７６】さらに特定すると、ブロック１５におい
て、下記の『適合』作業が行われる。すなわち、
【数１４】
【００７７】前記方程式の解決が下記であることを示す
ことができる： Ｎｐ（Ｘ，Ｙ）＝＜ＮE （Ｘ，Ｙ）．Ｊ（Ｅ）＞ （１５） ｑk （Ｘ，Ｙ）＝＜ＮE （Ｘ，Ｙ）．Ｇk （Ｅ）＞ （１６）
【００７８】上記の式において、
【００７９】
【数１７】 Ｇk （Ｅ）＝Θk （Ｅo ，Ｅ）−＜ＰΘk ＞・Ｊ（Ｅ） （１８）
【００８０】Ｊ（Ｅ）とＧk （Ｅ）はデータから独立し
ているので、Θ1 ，Θ2 ．．．ΘkとＰを使用する計算
Ｊ（Ｅ）を示すブロック５４によって、図２に示されて
いる通り、演繹的に誘導することができる。『ピクセル
当たり』作業が必要とするのは下記のスカラー積の評価
だけである： Ｎｐ＝＜ＮE ・Ｊ（Ｅ）＞ （１９）
【００８１】したがって、この適合は、ピクセル当たり
の、コンプトン拡散のないカウントＮp を決定する。
【００８２】エネルギー・ウインドウ当たりピクセル当
たりの統計の低さを補償するため、大きな空間域上の散
乱分布の相対安定度が『準局所的』解決の使用による利
点に与えられる。さらに特定すると、拡張された、すな
わち、『大きな』ピクセルが使用される。したがって、
最上部のピクセルが（Ｘ，Ｙ）である場合、空間ウイン
ドウＷは、面積ｓ＝（２Ｗ＋１）2 の式（２０）として
定義される： Ｘ0 −Ｗ＜Ｘ＜Ｘ0 ＋Ｗ3 Ｙ0 −Ｗ＜Ｙ＜Ｙ0 ＋Ｗ （２０）
【００８３】ウインドウｓ全体のコンプトン成分ＣS を
計算する場合『ピクセル当たり』コンプトン活動度Ｃ1
は、その平均によって近似させることができる： Ｃ1 ＝ＣS ／ｓ （２１）
【００８４】空間ウインドウｓ（座標（Ｘ，Ｙ）の周囲
で対称）において測定された活動度はＮE S （Ｘ，Ｙ）
で示され、（Ｘ，Ｙ）−ピクセル活性度はＮE 2 （Ｘ，
Ｙ）で示される。局所的な光電ピークを発見するため、
１個のパラメータの適合が行われる。局所的な光電ピー
クが式（２２）によって与えられることを示すことがで
きる：
【００８５】
【数２２】
【００８６】式（２２）において、
【００８７】
【数２３】
【００８８】
【数２４】
【００８９】Ｎ1 p （Ｘ，Ｙ）を解くと、コンプトン効
果がない映像のカウント／ピクセルが与えられる。
【００９０】別の適合法として、最大公算方法がある。
測定された活動度（ＮE ）が与えられた場合、試行関数
ｎ（Ｅ）のパラメータに関する同時ポアゾン確率が最大
化される。すなわち、式（２５）であるｎ（Ｅ）が発見
される。
【００９１】
【数２５】 または、Ｐは正なのであるから、 Ｉｎｐ＝＜ＮE Ｉｎ ｎ（Ｅ）−ｎ（Ｅ）−ＩｎＮE ！＞＝最大 （２６） である。
【００９２】最大公算解決のためには、 ＜ｎ（Ｅ）＞＝＜ＮE ＞ であることを示すことができる。この式は、ｎ（Ｅ）関
数からＮp を取り除くことができる。
【００９３】以下、（Ｘ，Ｙ）は明白には示されず、例
えば、ｎ（Ｘ，Ｙ；Ｅ）＝ｎ（Ｅ）となる。したがっ
て、方程式（１１）と（１２）から下記が導かれる：
【００９４】
【数２７】
【００９５】ｑk に関するＩｎＰの導関数を計算し、得
られた方程式を最大条件が要求する通り０に設定する
と、下記の方程式が得られる：
【００９６】
【数２８】
【００９７】これは非直線結合方程式であり、閉じた形
態で解決することはできない。多重グラジェント方法を
使用すると、ｑk ｓに関する対話解決を得ることができ
る。δｑk が対話の前のｑk と対話の後のｑ’k との間
の差であるとすると、 δｑk ＝ｑ’k −ｑk
【００９８】方程式の結合セットは直線化され、解決可
能になる：
【数２９】
【００９９】アレーｑk に関する解決の適切な収束が行
われた後、コンプトン効果がない活動度ｎp を式（３
０）から得ることができる：
【０１００】
【数３０】
【０１０１】さらに別の適合は、部分的最大公算解決で
ある。最小自乗解決はコンプトン成分の近似関数構造を
提供する。ただし、結合分布を最適化するためには、コ
ンプトン・イベント分数に対する光電ピークの比率を変
更することによって、ポアゾン統計を導入することが望
まれる。したがって、試行関数ｎ（Ｅ）は、式（３１）
である： ｎ（Ｅ）＝＜ＮE ＞［ｆp ．Ｐ＋（１−ｆp ）Ｃ］ （３１）
【０１０２】式（３１）において、
【数３２】 であり、Ｃは１に正規化された最小自乗コンプトン解決
である。すなわち、 Ｃ＝Ｃ／＜Ｃ＞ （３３）
【０１０３】ここで、最大公算方程式が、１個のパラメ
ータ、すなわち、光電ピーク関数ｆ p に関して最大化さ
れる。ひとたびｆp が計算されると、その方程式を使用
して、散乱のないイベント分布Ｎp を発見することがで
きる： Ｎp ＝ｆp ・＜ＮE ＞ （３４）
【０１０４】ｆp に関する微分から得られる最適化方程
式は下記となる：
【０１０５】
【数３５】
【０１０６】この式は、対話型ニュートン−ラフソン法
によって解くことができる： ｆ’p ＝ｆp ＋＜ＮE ．Ｔ＞／＜ＮE ．Ｔ2 ＞ （３６） 式（３６）において、
【数３７】
【０１０７】Ｎp の値を得ることに関係するまた別の方
法には、半局所的最大公算適合解決法が関与する。準局
所的解決法の場合は、下記の通り実行される：先ず、ピ
クセルＸo ，Ｙo を囲むＳと呼ばれる正方形に関する解
が得られる。すなわち、
【数３８】
【０１０８】完全最大公算法または部分的な最大公算法
のいずれかにより、ウインドウ全体のコンプトン効果が
ない成分がひとたび得られると、すなわち、Ｎ
p S （Ｘ，Ｙ）が既知となると、 Ｃ1 ＝ＣS ／ｓ （３９） であると仮定することによって、（Ｘ，Ｙ）−ピクセル
・コンプトン効果がない活動性を得るため、遅いコンプ
トン空間変動が使用される：
【数４０】
【０１０９】作動時、この発明のシステムはエネルギー
・スペクトルを局所的に解析し、非散乱光電ピーク関数
と、コンプトン散乱スペクトルを代表する関数の組み合
せから成る試行関数に適合させる。コンプトン散乱スペ
クトルを代表する関数は、仁科−ケルビン公式を使用し
て誘導される。したがって、コンプトン散乱スペクトル
の形状は、以前に使用された任意の多項式の場合とは異
なり、固有の性質として、コンプトン散乱の真の相対論
的分布を反映する。
【０１１０】仁科−ケルビン公式は、多重散乱コンプト
ン分布πm を作るために回帰的に使用される。その後、
Ｃm を得るため、各πm はシステム・エネルギー拡散関
数に繰り込まれる。計算のために要求される個別アレー
を得るため、繰り込まれた関数は個別のウインドウにつ
いて平均化または積分される。必要な計算の数を減少さ
せるため、および、この発明の方法が、取得の後数秒以
内に事実上コンプトン・フリーである映像を提供するこ
とができることを保証するため、個別関数Ｃ1のセット
が標準直交化されることが望ましい。標準直交化された
関数の係数は、ピクセル当たりの散乱されたカウントを
提供するパラメータである。このパラメータは、光電ピ
ーク成分と散乱成分から成る最終関数を、散乱された光
子と非散乱光子の両方を含む測定されたエネルギー分布
に適合させることによって決定される。局所的（および
準局所的）適合は適合関数の係数を最適化する。
【０１１１】この発明の独自のアプローチは、決定しな
ければならないパラメータが物理的コンプトン散乱関数
の係数であることである。前記関数はどの点においても
適合が正しいことを保証する正確な高エネルギーしきい
行動を有している。
【０１１２】また、この発明の方法が、コンプトン分布
（準局所的方法）の円滑さという長所を利用する方法に
よって、コンプトン適合のための計算のための統計を改
善することが望ましい。すなわち、望ましい実施例にお
いて、データは、適合のため、（２ｎ＋１）平方のピク
セル上で総和される（ｎは整数である）。その後、値
は、その正方形の中央のピクセルのみに与えられる。
【０１１３】未知の係数、すなわち非散乱イベントの場
合のピクセル当たりの量を解決するため（および、散乱
イベントのピクセル当たりの量が望まれる場合）、最小
自乗適合を使用することが望ましい。ただし、最大公算
適合を採用するいくつかの変わった方法も説明されてい
る。
【０１１４】具体的な実施例に関してこの発明を説明し
て来たが、前記の説明は、添付の特許請求の範囲によっ
て定義されたこの発明の範囲を限定するものとしてでは
なく、例としてなされたものに過ぎないことを理解しな
ければならない。
【０１１５】前記のこの発明の目的と特長は、その追加
目的と特長と共に、添付図面を参照してなされる説明に
照らして考慮すれば、最もよく理解されるであろう。
【図面の簡単な説明】
【図１】コンプトン・フリー映像を提供するガンマ放射
線映像システムを示すブロック図
【図２】図１に示されたシステムにおいて使用される準
備、計算と作動の詳細を示す図
【符号の説明】
１１ ブロック図 １４ 試行関数準備
段階 １５ カーブ適合または計算段階 １６ コンプトン効
果がない映像 １７ ガンマ線輻射検出器 １８ ガンマ線輻射検出器の面に衝突する光子 １９ 導体 ２１ 導体 ２２ 導体 ２３ 座標コンピュ
ータ ２４ 導体 ２５ エネルギー
（Ｚ）補正回路 ２６ エネルギー処理回路 ２７ 導体 ２８ ２９ ３１ 映像補正器とデジタイザ回路 ３２ マトリックス ３３ 映像予備処理
装置 ３４ ブロック ３６ 準備ブロック ４１ ブロック ４２ ブロック ４３ ブロック ４４ ブロック ４６ ブロック ４７ ブロック ４８ ブロック ４９ コンピュータ ５１ ブロック ５２ ブロック ５３ ブロック ５４ ブロック
─────────────────────────────────────────────────────
【手続補正書】
【提出日】平成３年９月５日
【手続補正１】
【補正対象書類名】明細書
【補正対象項目名】００５１
【補正方法】変更
【補正内容】
【００５１】図１の試行関数段階１４は、散乱イベント
と非散乱イベントを含む理論的すなわち試行エネルギー
分布ｎ（Ｘ，Ｙ，Ｅ）を作る。この場合：
【数１】
【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】００５２
【補正方法】変更
【補正内容】
【００５２】上記の式において、 ε＝Ｅ／ｍe ｃ2 ：電子静止エネルギーｍe ｃ2 を単位
とする光子エネルギー Ｐ（ε’，ε）＝エネルギーε’におけるシステム・エ
ネルギー拡散関数。その形状は測定されており既知であ
る。 п（εo ，ε）＝原エネルギーεo からエネルギーεま
でｍ回散乱された光子によって作られたイベントのエネ
ルギー分布（すなわち、ｍ回散乱された光子のエネルギ
ー確率分布の形状） ＮP （ｘ，ｙ）＝散乱しなかった光子によって作られた
イベントの空間分布（カウント／ピクセル） Ｑm （ｘ，ｙ）＝ｍ回散乱した光子によって作られたイ
ベントの空間分布（カウント／ピクセル） εo ＝光子の原エネルギー ε＝測定された光子のエネルギー ε’＝光子の中間エネルギー
【手続補正３】
【補正対象書類名】明細書
【補正対象項目名】００６５
【補正方法】変更
【補正内容】
【００６５】
【数５】
【手続補正４】
【補正対象書類名】明細書
【補正対象項目名】００６８
【補正方法】変更
【補正内容】
【００６８】標準直交化はコンピュータ４９内で行わ
れ、結果、すなわち、Θ1 ，Θ2 ……Θm は、例えば、
ブロック５１，５２，５３中に示される。
【手続補正５】
【補正対象書類名】明細書
【補正対象項目名】００７１
【補正方法】変更
【補正内容】
【００７１】［標準直交ベース｛Θm ｝を使用すると、
どのようなベクトルＶでもΘm ’のアレーの重ね合わせ
によって示される：
【手続補正６】
【補正対象書類名】明細書
【補正対象項目名】００７６
【補正方法】変更
【補正内容】
【００７６】さらに特定すると、ブロック１５におい
て、下記の『適合』作業が行われる。すなわち、
【数１４】
【手続補正７】
【補正対象書類名】明細書
【補正対象項目名】００７９
【補正方法】変更
【補正内容】
【００７９】
【数１７】 Ｇk （Ｅ）＝Θk （Ｅo ，Ｅ）−＜ＰΘk ＞Ｊ（Ｅ） （１８）
【手続補正８】
【補正対象書類名】明細書
【補正対象項目名】００８２
【補正方法】変更
【補正内容】
【００８２】エネルギー・ウインドウ当たりピクセル当
たりの統計の低さを補償するため、大きな空間域上の散
乱分布の相対安定度が『準局所的』解決の使用による利
点に与えられる。さらに特定すると、拡張された、すな
わち、『大きな』ピクセルが使用される。したがって、
最上部のピクセルが（Ｘ，Ｙ）である場合、空間ウイン
ドウＷは、面積ｓ＝（２Ｗ＋１）2 の式（２０）として
定義される： Ｘ0 −Ｗ＜Ｘ＜Ｘ0 ＋Ｗ ， Ｙ0 −Ｗ＜Ｙ＜Ｙ0 ＋Ｗ （２０）
【手続補正９】
【補正対象書類名】明細書
【補正対象項目名】００８４
【補正方法】変更
【補正内容】
【００８４】空間ウインドウｓ（座標（Ｘ，Ｙ）の周囲
で対称）において測定された活動度はＮE S （Ｘ，Ｙ）
で示され、（Ｘ，Ｙ）−ピクセル活性度はＮE 1 （Ｘ，
Ｙ）で示される。局所的な光電ピークを発見するため、
１個のパラメータの適合が行われる。局所的な光電ピー
クが式（２２）によって与えられることを示すことがで
きる：
【手続補正１０】
【補正対象書類名】明細書
【補正対象項目名】００８５
【補正方法】変更
【補正内容】
【００８５】
【数２２】
【手続補正１１】
【補正対象書類名】明細書
【補正対象項目名】００８８
【補正方法】変更
【補正内容】
【００８８】
【数２４】
【手続補正１２】
【補正対象書類名】明細書
【補正対象項目名】００９１
【補正方法】変更
【補正内容】
【００９１】
【数２５】 または、Ｐは正なのであるから、 ｌｎＰ＝＜ＮE ｌｎ ｎ（Ｅ）−ｎ（Ｅ）−ｌｎＮE ！＞＝最大 （２６） である。
【手続補正１３】
【補正対象書類名】明細書
【補正対象項目名】００９５
【補正方法】変更
【補正内容】
【００９５】ｑk に関するｌｎＰの導関数を計算し、得
られた方程式を最大条件が要求する通り０に設定する
と、下記の方程式が得られる：
【手続補正１４】
【補正対象書類名】明細書
【補正対象項目名】００９６
【補正方法】変更
【補正内容】
【００９６】
【数２８】
【手続補正１５】
【補正対象書類名】明細書
【補正対象項目名】００９７
【補正方法】変更
【補正内容】
【００９７】これは非直線結合方程式であり、閉じた形
態で解決することはできない。多重グラジェント方法を
使用すると、ｑk ｓに関する対話解決を得ることができ
る。δｑk が対話の前のｑk と対話の後のｑ’k との間
の差であるとすると、 δｑk ＝ｑk −ｑ’k
【手続補正１６】
【補正対象書類名】明細書
【補正対象項目名】０１０５
【補正方法】変更
【補正内容】
【０１０５】
【数３５】
【手続補正１７】
【補正対象書類名】明細書
【補正対象項目名】０１０６
【補正方法】変更
【補正内容】
【０１０６】この式は、対話型ニュートン−ラフソン法
によって解くことができる： ｆ’p ＝ｆp ＋＜ＮE ．τ＞／＜ＮE ．τ2 ＞ （３６） 式（３６）において、
【数３７】
【手続補正１８】
【補正対象書類名】明細書
【補正対象項目名】０１０７
【補正方法】変更
【補正内容】
【０１０７】Ｎp の値を得ることに関係するまた別の方
法には、半局所的最大公算適合解決法が関与する。準局
所的解決法の場合は、下記の通り実行される：先ず、ピ
クセルＸo ，Ｙo を囲むＳと呼ばれる正方形に関する解
が得られる。すなわち、
【数３８】
───────────────────────────────────────────────────── フロントページの続き (72)発明者 イガル シュレム イスラエル国 ハイファ， シムソン 17 エー (72)発明者 アドリアン ソイル イスラエル国 ハイファ， ネベシャナ ン， アバ ヒレル シルバー 113／42

Claims (45)

    【特許請求の範囲】
  1. 【請求項１】ガンマ線映像システムが作る映像に及ぼす
    コンプトン散乱光子の効果の影響を減少させる方法であ
    って、 源から来るガンマ線輻射に応答してガンマ線検出器に衝
    突する光子をカウントとして検出する段階；検出器上の
    衝突位置に従い、各カウントについてＸ，Ｙ座標位置を
    決定する段階；衝突する光子のエネルギーを測定する段
    階；共にエネルギー範囲をまたいでおり、それぞれが前
    記エネルギー範囲の一部分上を延伸する複数のエネルギ
    ー・ウインドウを使用する段階；測定された衝突光子の
    測定されたエネルギーに従って、検出された各カウント
    をメモリ・マトリックス上に置く段階；前記マトリック
    スは空間Ｘ，Ｙ座標から成っており、それぞれが１個の
    エネルギー・ウインドウに対応するデジタル映像であ
    り、前記検出されたカウントは、検出器上の衝突の位置
    に従って前記Ｘ，Ｙ座標に置かれる；メモリ・マトリッ
    クス内の各Ｘ，Ｙ座標位置におけるカウントを集積する
    段階；マトリックス当たりの位置当たりで集積されたカ
    ウントを使用し、各Ｘ，Ｙ位置の測定されたエネルギー
    ・スペクトルを作る段階；前記測定されたエネルギー・
    スペクトルには散乱された光子と散乱されない光子とが
    作るカウントが含まれる；散乱されない光子のみが作る
    Ｘ，Ｙ位置当たりのカウントを得る段階；前記散乱され
    ない光子のみが作るＸ，Ｙ位置当たりのカウントを得る
    段階は、１個または複数のオーダーのコンプトン散乱光
    子のエネルギー分布を計算する段階を含んでいる；位置
    当たりのシステム依存性コンプトン・エネルギー分布を
    得るため、システムのエネルギー拡散関数でコンプトン
    散乱光子のエネルギー分布を繰り込む段階；システム・
    エネルギー拡散関数と非散乱光子の未知の係数との積、
    および、散乱光子の未知の係数とシステム依存性コンプ
    トン・エネルギー分布との積とから成る試行関数を作る
    段階；非散乱光子の未知の係数を解くことによって非散
    乱光子のカウントを得るため、測定されたエネルギー・
    スペクトルに前記試行関数を局所的に適合させる段階；
    および事実上コンプトン・フリーである映像を作るため
    に、非散乱光子のカウントを使用する段階；から成るこ
    とを特徴とする方法。
  2. 【請求項２】ガンマ線映像システムが作る映像に及ぼす
    コンプトン散乱光子の効果の影響を減少させる方法であ
    って、 源から来るガンマ線輻射に応答してガンマ線検出器に衝
    突する光子をカウントとして検出する段階；検出器上の
    光子の衝突位置に従い、各カウントについてＸ，Ｙ座標
    位置を決定する段階；衝突する光子のエネルギーを測定
    する段階；共にエネルギー範囲上にまたがる複数のエネ
    ルギー・ウインドウを使用する段階；測定された衝突光
    子のエネルギーに従って、検出された各カウントを、複
    数のメモリ・マトリックスのうち前記複数のエネルギー
    ・ウインドウの１個に対応する１個中に置く段階；前記
    マトリックスはそれぞれ空間Ｘ，Ｙ座標から成ってお
    り、前記マトリックスはそれぞれがデジタル映像であっ
    て１個のエネルギー・ウインドウに対応しており、前記
    Ｘ，Ｙ座標は検出器上の光子の衝突の位置と映像中のピ
    クセルの位置に対応し、；各マトリックスの各Ｘ，Ｙ座
    標位置において集積されたカウントを使用して、Ｘ，Ｙ
    位置当たりの測定されたエネルギー・スペクトルを提供
    する段階；前記測定されたエネルギー・スペクトルには
    散乱された光子と散乱されない光子とが作るカウントが
    含まれる；物理相対論的コンプトン特性を使用して、１
    個または複数のオーダーのコンプトン散乱光子のエネル
    ギー分布を計算する段階；１個または複数のオーダーの
    コンプトン散乱光子のシステム依存性コンプトン・エネ
    ルギー分布を得るため、１個または複数のオーダーのコ
    ンプトン散乱光子の計算されたエネルギー分布をシステ
    ムのエネルギー拡散関数で繰り込む段階；未知のＸ，Ｙ
    係数が付いた、１個または複数のオーダーのコンプトン
    散乱光子の、繰り込まれたシステム依存性エネルギー分
    布と、未知のＸ，Ｙ係数が付いた光電ピークエネルギー
    形状の影響とを合計することによって、試行関数を作る
    段階；非散乱光子の未知の係数を解くことによって非散
    乱光子のカウントを得るため、測定されたエネルギー・
    スペクトルに前記試行関数を局所的に適合させる段階；
    およびコンプトン・フリーである映像を提供するため
    に、ピクセル当たりの非散乱光子のカウントを使用する
    段階；から成ることを特徴とする方法。
  3. 【請求項３】散乱光子のエネルギー分布を計算する前記
    段階が、 第１オーダーの拡散光子のエネルギー分布を決定するた
    め、仁科−ケルビン方程式を第１オーダー光子拡散確率
    分布に変換する段階；およびさらにオーダーが高い分散
    光子の前記分布が使用される場合は、さらにオーダーが
    高い散乱光子の確率分布を得るため、第１オーダーの分
    散光子のエネルギー分布を繰り込む段階；から成ること
    を特徴とする、請求項２に記載の方法。
  4. 【請求項４】各エネルギー・ウインドウについて計算さ
    れたコンプトン・エネルギー分布を平均することによっ
    て、コンプトン多重散乱光子のエネルギー分布から個別
    の関数のセットを取得する段階を含むことを特徴とす
    る、請求項３に記載の方法。
  5. 【請求項５】計算の数を減少させる段階を含むことを特
    徴とする、請求項４に記載の方法。
  6. 【請求項６】計算の数を減少させる段階が、個別の関数
    の正規直交関数セットへの変換から成ることを特徴とす
    る、請求項５に記載の方法。
  7. 【請求項７】適合段階が、 散乱光子の局所的カウントを得るために準局所ピクセル
    を使用すること；および局所的光電ピーク・カウントを
    決定するため１回のパラメータ適用を行うこと；を含む
    ことを特徴とする、請求項６に記載の方法。
  8. 【請求項８】準局所ピクセルを使用する段階が、 （２ｎ＋１）×（２ｍ＋１）の四角形ピクセルから成る
    大きなピクセルを使用すること［ｎとｍは共に正の整数
    であり、大きなピクセルの値は大きなピクセル中のピク
    セル数によって分割され、中心ピクセルに与えられ
    る］；およびこの方法ですべてのピクセルを評価するこ
    と；から成ることを特徴とする、請求項７に記載の方
    法。
  9. 【請求項９】試行関数を測定されたエネルギー・スペク
    トルに適合させる段が最小自乗適合を使用して行われる
    ことを特徴とする、請求項２に記載の方法
  10. 【請求項１０】試行関数を測定されたエネルギー・スペ
    クトルに適合させる段階が最大公算適合を使用して行わ
    れることを特徴とする、請求項２に記載の方法
  11. 【請求項１１】試行関数を測定されたエネルギー・スペ
    クトルに適合させる段階が部分的最大公算適合を使用し
    て行われることを特徴とする、請求項２に記載の方法
  12. 【請求項１２】試行関数を測定されたエネルギー・スペ
    クトルに適合させる段階が最小自乗と最大公算適合の組
    み合せを使用して行われることを特徴とする、請求項２
    に記載の方法
  13. 【請求項１３】ガンマ線映像システムで事実上コンプト
    ン・フリーである患者の映像を作る方法であって、 既知のシステム・エネルギー拡散関数に、非散乱光子で
    あるため未知であるカウントを乗じること；システム依
    存性のコンプトン適合関数のセットを提供するため、光
    子の物理的相互作用の確率から誘導される散乱光子の少
    なくとも１個のエネルギー分布を既知のシステム・エネ
    ルギー拡散関数を繰り込むこと：およびシステム依存性
    コンプトン適合関数に、散乱光子のカウントに相当する
    未知数を乗じることの各段階から成る試行関数を決定す
    る段階；散乱した光子と散乱しない光子に起因する総カ
    ウントを含む患者のエネルギー・スペクトルを局所的に
    測定する段階；および事実上コンプトン・フリーである
    映像を作るために非散乱光子に起因するカウント数を使
    用する段階；から成ることを特徴とする方法。
  14. 【請求項１４】エネルギー・スペクトルを測定する段階
    を終えた後数秒以内に、非散乱光子に起因する未知のカ
    ウント数を決定することが含まれることを特徴とする、
    請求項１３に記載の方法。
  15. 【請求項１５】試行エネルギー関数を決定する段階に、
    仁科−ケルビン方程式を第１オーダー光子散乱確率エネ
    ルギー分布への変換、および、第１オーダー散乱光子の
    分布を得るための変換された仁科−ケルビン方程式の使
    用が含まれることを特徴とする、請求項１３に記載の方
    法。
  16. 【請求項１６】高いオーダーが考慮される場合、高いオ
    ーダーの散乱光子のエネルギー分布を得るため、第１オ
    ーダー散乱光子のエネルギー分布を繰り込む段階を含む
    ことを特徴とする、請求項１５に記載の方法。
  17. 【請求項１７】各マトリックスについてコンプトン・エ
    ネルギー分布を統合することによって、計算されたコン
    プトン・エネルギー分布から個別の関数のセットを得る
    段階を含むことを特徴とする、請求項１６に記載の方
    法。
  18. 【請求項１８】エネルギー関数を測定する段階に、 システムの光子検出器に衝突する光子のＸ，Ｙ位置とエ
    ネルギーを決定する段階；衝突光子のカウントを集積す
    るため、衝突位置に対応するＸ，Ｙピクセルを有する複
    数のマトリックスを使用する段階；前記各マトリックス
    は複数のエネルギー・ウインドウの１個に対応し、前記
    複数のエネルギー・ウインドウはエネルギー範囲をまた
    いでいる；およびマトリックス中の集積カウントから測
    定されたエネルギー・スペクトルを作る段階；が含まれ
    ることを特徴とする、請求項１７に記載の方法。
  19. 【請求項１９】個別の関数を標準直交化する段階が含ま
    れることを特徴とする、請求項１８に記載の方法。
  20. 【請求項２０】標準直交化のためのグラハム−シュミッ
    ト手順の使用が含まれることを特徴とする、請求項１９
    に記載の方法。
  21. 【請求項２１】測定されたエネルギー・スペクトルに試
    行関数を適合する段階に最小自乗適合の使用が含まれる
    ことを特徴とする、請求項１３に記載の方法。
  22. 【請求項２２】測定されたエネルギー・スペクトルに試
    行関数を適合する段階に最大公算解決の使用が含まれる
    ことを特徴とする、請求項１３に記載の方法。
  23. 【請求項２３】適合の段階が最小自乗適合と最大公算適
    合の組み合わせを含むことを特徴とする、請求項１３に
    記載の方法。
  24. 【請求項２４】ガンマ線映像システムが作る映像に及ぼ
    すコンプトン散乱光子の効果の影響を減少させるシステ
    ムであって、 源から来るガンマ線輻射に応答してガンマ線検出器に衝
    突する光子をカウントとして検出する手段；検出器上の
    光子の衝突位置に従い、各カウントについてＸ，Ｙ座標
    位置を決定する手段；衝突する光子のエネルギーを測定
    する手段；共にエネルギー範囲上をまたがり、それぞれ
    が前記エネルギー範囲の一部分上を延伸する複数のエネ
    ルギー・ウインドウを使用する手段；測定された衝突光
    子のエネルギーに従って、検出された各カウントを前記
    複数のメモリ・マトリックスの１個のメモリ・マトリッ
    クス上に挿入する手段；前記マトリックスは空間Ｘ，Ｙ
    座標から成っており、前記検出されたカウントは検出器
    上の衝突の位置に従ってＸ，Ｙ座標位置において挿入さ
    れる；前記メモリ・マトリックスは各Ｘ，Ｙ座標位置に
    おいてカウントを集積する；マトリックス当たり位置当
    たりで集積されたカウントを使用し、各Ｘ，Ｙ位置の測
    定されたエネルギー・スペクトルを作る手段；前記測定
    されたエネルギー・スペクトルには散乱された光子と散
    乱されない光子が作るカウントが含まれる；散乱されな
    い光子のみが作るＸ，Ｙ位置当たりのカウントを得る手
    段；散乱されない光子のみが作るＸ，Ｙ位置当たりのカ
    ウントを得る手段には、コンプトン散乱光子の１個また
    は複数のオーダーのエネルギー分布を計算する手段が含
    まれる；システム依存性コンプトン・エネルギー分布を
    得るため、コンプトン散乱光子の１個または複数のオー
    ダーのエネルギー分布をシステムのエネルギー拡散関数
    を繰り込むための手段；システム・エネルギー拡散関数
    と非散乱光子の未知の係数との積、および、散乱光子の
    未知の係数とシステム依存性コンプトン・エネルギー分
    布との積との和から成る試行関数を作る手段；非散乱光
    子の未知の係数を解くことによって非散乱光子のカウン
    トを得るため、測定されたエネルギー・スペクトルに前
    記試行関数を部分的に適合させる手段；および事実上コ
    ンプトン・フリーである映像を作るために、非散乱光子
    のカウントを使用する手段；から成ることを特徴とする
    システム。
  25. 【請求項２５】作られた映像に及ぼすコンプトン散乱光
    子の影響を減少させるガンマ線映像システムであって、 源から来るガンマ線輻射に応答してガンマ線検出器に衝
    突する光子をカウントとして検出する手段；検出器上の
    光子の衝突位置に従い、検出された各カウントについて
    Ｘ，Ｙ座標位置を決定する手段；衝突する光子のエネル
    ギーを測定する手段；エネルギー範囲に広がり、それぞ
    れが１個のエネルギー・ウインドウに対応する複数のメ
    モリ・マトリックス；測定された衝突光子のエネルギー
    に従って、検出された各カウントを前記複数のメモリ・
    マトリックスの１個のメモリ・マトリックス上に置く手
    段；前記各メモリ・マトリックスはＸ，Ｙ座標から成っ
    ており、前記各マトリックスはＸ，Ｙ座標から成ってお
    り、前記各マトリックスはデジタル映像であり、前記
    Ｘ，Ｙ座標は検出器上の光子の衝突の位置と映像上の１
    個のピクセルの位置に対応する；各マトリックスの各
    Ｘ，Ｙ座標位置において集積されたカウントを使用し、
    Ｘ，Ｙ位置あたりの測定されたエネルギー・スペクトル
    を提供する手段；前記測定されたエネルギー・スペクト
    ルには散乱された光子と散乱されない光子が作るカウン
    トが含まれる；物理相対論コンプトン特性を使用して、
    コンプトン多重散乱成分のエネルギー分布を計算する手
    段；前記多重散乱には、少なくとも第１オーダーの散乱
    が含まれる；未知の局所的係数が付いた多重散乱光子の
    システム依存性コンプトン・エネルギー分布を得るた
    め、システム・エネルギー拡散関数が付いたコンプトン
    多重散乱光子の計算されたエネルギー分布を繰り込むた
    めの手段；少なくとも第１オーダーの拡散についての、
    未知の局所係数が付いた、繰り込まれた多重散乱光子エ
    ネルギー分布と未知の局所係数が付いた光電ピーク・エ
    ネルギー形状の影響とを合計する手段を含む試行関数を
    作るための手段；非散乱光子の未知の係数を解くことに
    よって非散乱光子の位置当たりのカウントを決定するた
    め、測定されたエネルギー・スペクトルに前記試行関数
    を部分的に適合させる手段；およびコンプトン・フリー
    である映像を作るために、ピクセル当たりの非散乱光子
    のカウントを使用する手段；から成ることを特徴とする
    ガンマ線映像システム。
  26. 【請求項２６】前記散乱された光子のエネルギー分布を
    計算するための手段が、 第１オーダーの散乱光子のエネルギー分布を決定するた
    め、第１オーダーの光子散乱確率分布に変換された仁科
    −ケルビン方程式を使用する手段；および高いオーダー
    の散乱光子の分布を使用する場合は、高いオーダーの散
    乱光子の確率分布を得るため、第１オーダーの散乱光子
    のエネルギー分布を繰り込むための手段；から成ること
    を特徴とする請求項２５に記載のシステム。
  27. 【請求項２７】計算された各エネルギー・ウインドウの
    コンプトン・エネルギー分布を平均することによって、
    コンプトン多重散乱光子のエネルギー分布から個別の関
    数のセットを得るための手段を含むことを特徴とする、
    請求項２６に記載のシステム。
  28. 【請求項２８】計算数を減少させるための手段が含まれ
    ることを特徴とする、請求項２７に記載のシステム。
  29. 【請求項２９】計算数を減少させるための手段が、個別
    の関数を関数の標準直交セットに変換するための手段か
    ら成ることを特徴とする、請求項２８のシステム。
  30. 【請求項３０】適合のための手段に、散乱した光子の局
    所的カウントを得るための準局所ピクセルと局所的光電
    ピーク・カウントを決定するために単独パラメータ適合
    を行うための手段が含まれることを特徴とする、請求項
    ２８のシステム。
  31. 【請求項３１】準局所ピクセルが、 四角形（２ｎ＋１）×（２ｍ＋１）［ｎとｍは正の整
    数］で構成された大きなピクセル；前記大きなピクセル
    の値を前記大きなピクセル中のピクセル数で除算するた
    めの手段；前記除算の商をピクセル値として中心ピクセ
    ルに与えるための手段；およびこの方法ですべてのピク
    セルを評価するための手段；とから成ることを特徴とす
    る、請求項３０のシステム。
  32. 【請求項３２】測定されたエネルギー・スペクトルに試
    行関数を適合させる手段が最小自乗適合から成ることを
    特徴とする、請求項２５のシステム。
  33. 【請求項３３】測定されたエネルギー・スペクトルに試
    行関数を適合させる手段が最大公算適合から成ることを
    特徴とする、請求項２５のシステム。
  34. 【請求項３４】測定されたエネルギー・スペクトルに試
    行関数を適合させる手段が最小自乗適合と最大公算適合
    を組み合わせた手段から成ることを特徴とする、請求項
    ２５のシステム。
  35. 【請求項３５】事実上コンプトン効果のない患者の映像
    をガンマ線映像システムで作るシステムであって、 既知のシステム・エネルギー拡散関数に非散乱光子に起
    因する未知のカウント数を乗ずるための手段；システム
    依存性コンプトン・エネルギー分布を提供するため、少
    なくとも、光子の物理的相互作用の確率から誘導される
    散乱光子の第１オーダー分布と既知のシステム・エネル
    ギー拡散関数とを繰り込むための手段；およびシステム
    依存性コンプトン・エネルギー分布に散乱光子のカウン
    トに対応する未知の数を乗ずるための手段から成る、試
    行関数を決定するための手段；散乱光子と非散乱光子に
    起因するすべてのカウントを含む患者のエネルギー・ス
    ペクトルを測定するための手段；非散乱光子に起因する
    未知のカウント数を決定するため、測定されたエネルギ
    ー・スペクトルに試行関数を適合させるための手段；お
    よび事実上コンプトン・フリーである映像を作るため、
    非散乱光子に起因するカウント数を使用するための手
    段；から成ることを特徴とするシステム。
  36. 【請求項３６】エネルギー・スペクトル測定の仕上げの
    数秒以内に、非散乱光子に起因する未知のカウント数を
    決定するための手段を含むことを特徴とする、請求項３
    ５のシステム。
  37. 【請求項３７】試行エネルギー関数を決定するための手
    段に、仁科−ケルビンの方程式を第１オーダー光子散乱
    確率エネルギー分布への変換；および第１オーダー散乱
    光子のエネルギー分布を得るために、変換された仁科−
    ケルビンの方程式を使用する手段；が含まれることを特
    徴とする、請求項３５のシステム。
  38. 【請求項３８】さらに高いオーダーの散乱光子のエネル
    ギー分布を得るため、第１オーダー散乱光子のエネルギ
    ー分布を繰り込むための手段を含むことを特徴とする、
    請求項３８のシステム。
  39. 【請求項３９】個別関数のセットを得るため、各マトリ
    ックスのコンプトン・エネルギー分布を統合させるため
    の手段を含むことを特徴とする、請求項３８のシステ
    ム。
  40. 【請求項４０】エネルギー関数を測定するための手段
    に、 システムの光子検出器に衝突する光子のＸ，Ｙ位置とエ
    ネルギーを決定するための手段；検出された光子のエネ
    ルギーを決定するための手段；衝突光子のカウントを集
    積するための、衝突位置に対応するＸ，Ｙピクセルを有
    する複数のマトリックス；前記各マトリックスは１個の
    エネルギー・ウインドウに対応し、前記複数のマトリッ
    クスはエネルギー範囲上に延伸する；およびマトリック
    ス中に集積されたカウントから測定されたエネルギー・
    スペクトルを作るための手段；が含まれることを特徴と
    する、請求項３９のシステム。
  41. 【請求項４１】個別関数を標準直交化するための手段が
    含まれることを特徴とする、請求項４０のシステム。
  42. 【請求項４２】直交化のためにグラハム−シュミットの
    手順を使用するための手段が含まれることを特徴とす
    る、請求項４１のシステム。
  43. 【請求項４３】測定されたエネルギー・スペクトルに試
    行関数を適合させるための手段に、最小自乗適合のため
    の手段が含まれることを特徴とする、請求項３５のシス
    テム。
  44. 【請求項４４】測定されたエネルギー・スペクトルに試
    行関数を適合させるための手段に、最大公算適合のため
    の手段が含まれることを特徴とする、請求項３５のシス
    テム。
  45. 【請求項４５】測定されたエネルギー・スペクトルに試
    行関数を適合させるための手段に、最小自乗適合と最大
    公算適合の組み合せのための手段が含まれることを特徴
    とする、請求項３５のシステム。
JP03139186A 1990-06-11 1991-06-11 コンプトン・フリー・ガンマ線映像 Expired - Fee Related JP3123024B2 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
IL94691A IL94691A0 (en) 1990-06-11 1990-06-11 Compton free gamma camera images
IL94691 1990-06-11

Publications (2)

Publication Number Publication Date
JPH075263A true JPH075263A (ja) 1995-01-10
JP3123024B2 JP3123024B2 (ja) 2001-01-09

Family

ID=11061295

Family Applications (1)

Application Number Title Priority Date Filing Date
JP03139186A Expired - Fee Related JP3123024B2 (ja) 1990-06-11 1991-06-11 コンプトン・フリー・ガンマ線映像

Country Status (7)

Country Link
US (1) US5293195A (ja)
JP (1) JP3123024B2 (ja)
CA (1) CA2042698A1 (ja)
DE (1) DE4118418A1 (ja)
FR (1) FR2663127B1 (ja)
GB (1) GB2246274B (ja)
IL (1) IL94691A0 (ja)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008039528A (ja) * 2006-08-03 2008-02-21 Fujifilm Ri Pharma Co Ltd 核医学データ処理装置、プログラム、記録媒体および方法
JP2008536138A (ja) * 2005-04-14 2008-09-04 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ Ｃｔにおけるエネルギー分布再構成
JP2008545135A (ja) * 2005-07-07 2008-12-11 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 多重散乱補正
JP2010517015A (ja) * 2007-01-17 2010-05-20 イノベイティブ アメリカン テクノロジー， インコーポレイテッド スペクトル解析のための高度パターン認識システム
JP2011089901A (ja) * 2009-10-22 2011-05-06 Sumitomo Heavy Ind Ltd 検出結果補正方法、その検出結果補正方法を用いた放射線検出装置、その検出結果補正方法を実行するためのプログラム、及びそのプログラムを記録する記録媒体
WO2014109374A1 (ja) * 2013-01-10 2014-07-17 株式会社 東芝 Ｘ線コンピュータ断層撮影装置、医用画像処理装置、および医用画像処理方法

Families Citing this family (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5371672A (en) * 1991-03-27 1994-12-06 Kabushiki Kaisha Toshiba Scintillation camera apparatus capable of quantitatively eliminating scattering signal components by setting multiple window and method for quantitatively eliminating scattering signal components
CA2070879A1 (en) * 1991-06-09 1992-12-10 Ge Medical Systems Israel Ltd. Gamma camera images having reduced artifacts
US5438202A (en) * 1992-05-28 1995-08-01 Elscint Ltd. Stabilized scatter free gamma camera images
US5434414A (en) * 1991-06-09 1995-07-18 Elscint Ltd. Reduction of image artifacts caused by unwanted photons
US5390225A (en) * 1992-06-29 1995-02-14 Siemens Medical Systems, Inc. Mapping slices of the human body into regions having a constant linear attenuation coefficient for correcting images acquired during a nuclear medicine study for attenuation artifacts
US5466939A (en) * 1992-09-24 1995-11-14 Shimadzu Corporation Nuclear medicine imaging apparatus
DE4233092A1 (de) * 1992-10-01 1994-04-07 Siemens Ag Verfahren zur Kompensation von Streuereignissen in einem Szintillationsbild
FR2705465B1 (fr) * 1993-05-18 1995-07-07 Inst Nat Sante Rech Med Procédé de décomposition d'images scintigraphiques en composantes d'absorption totale et diffusées.
US5561297A (en) * 1993-07-12 1996-10-01 Sopha Medical Systems, Inc. Scatter correcting gamma-ray camera
US5533091A (en) * 1995-04-28 1996-07-02 General Electric Company Noise suppression algorithm and system
WO1996041213A1 (en) * 1995-06-07 1996-12-19 Massachusetts Institute Of Technology X-ray detector and method for measuring energy of individual x-ray photons for improved imaging of subjects using reduced dose
US5633499A (en) * 1995-11-21 1997-05-27 Trionix Research Laboratory, Inc. Scatter elimination technique and apparatus in radionuclide emission and transmission imaging in a nuclear camera
US5689115A (en) * 1995-11-24 1997-11-18 Elscint Ltd. Advanced nuclear medicine system
US7085405B1 (en) * 1997-04-17 2006-08-01 Ge Medical Systems Israel, Ltd. Direct tomographic reconstruction
US6633626B2 (en) * 2002-02-01 2003-10-14 General Electric Company Methods and apparatus for correcting scatter
US7197172B1 (en) * 2003-07-01 2007-03-27 Analogic Corporation Decomposition of multi-energy scan projections using multi-step fitting
EP1946622A2 (en) * 2005-10-06 2008-07-23 Imaging Sciences International, Inc. Scatter correction
US8378310B2 (en) * 2009-02-11 2013-02-19 Prismatic Sensors Ab Image quality in photon counting-mode detector systems
JP6472336B2 (ja) 2014-06-18 2019-02-20 キヤノン株式会社 画像処理装置、画像処理方法およびプログラム
US10534096B1 (en) 2018-06-22 2020-01-14 Kromek Group, PLC Photon scatter imaging

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4780823A (en) * 1986-03-11 1988-10-25 Siemens Gammasonics, Inc. Multiple pixel area-weighted acquisition system for scintillation cameras
JPS6312985A (ja) * 1986-07-03 1988-01-20 Toshiba Corp 放射線エネルギ−弁別回路
US4839808A (en) * 1987-05-22 1989-06-13 The University Of Michigan Correction for compton scattering by analysis of energy spectra

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008536138A (ja) * 2005-04-14 2008-09-04 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ Ｃｔにおけるエネルギー分布再構成
JP2008545135A (ja) * 2005-07-07 2008-12-11 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 多重散乱補正
JP2008039528A (ja) * 2006-08-03 2008-02-21 Fujifilm Ri Pharma Co Ltd 核医学データ処理装置、プログラム、記録媒体および方法
JP2010517015A (ja) * 2007-01-17 2010-05-20 イノベイティブ アメリカン テクノロジー， インコーポレイテッド スペクトル解析のための高度パターン認識システム
JP2011089901A (ja) * 2009-10-22 2011-05-06 Sumitomo Heavy Ind Ltd 検出結果補正方法、その検出結果補正方法を用いた放射線検出装置、その検出結果補正方法を実行するためのプログラム、及びそのプログラムを記録する記録媒体
WO2014109374A1 (ja) * 2013-01-10 2014-07-17 株式会社 東芝 Ｘ線コンピュータ断層撮影装置、医用画像処理装置、および医用画像処理方法
JP2014147691A (ja) * 2013-01-10 2014-08-21 Toshiba Corp Ｘ線コンピュータ断層撮影装置、医用画像処理装置、および医用画像処理プログラム
US10001568B2 (en) 2013-01-10 2018-06-19 Toshiba Medical Systems Corporation X-ray computed tomography apparatus, medical image processing apparatus and medical image processing method

Also Published As

Publication number Publication date
JP3123024B2 (ja) 2001-01-09
GB9110628D0 (en) 1991-07-03
GB2246274B (en) 1994-12-07
FR2663127A1 (fr) 1991-12-13
FR2663127B1 (fr) 1993-11-05
IL94691A0 (en) 1991-04-15
GB2246274A (en) 1992-01-22
DE4118418A1 (de) 1991-12-12
CA2042698A1 (en) 1991-12-12
US5293195A (en) 1994-03-08

Similar Documents

Publication Publication Date Title
JPH075263A (ja) コンプトン・フリー・ガンマ線映像
Gagnon et al. Introduction to holospectral imaging in nuclear medicine for scatter subtraction
US7312455B2 (en) Method and system for scatter correction in a positron emission tomography system
US5438202A (en) Stabilized scatter free gamma camera images
US8705827B2 (en) Scatter correction methods
Beekman et al. Scatter compensation methods in 3D iterative SPECT reconstruction: a simulation study
US5434414A (en) Reduction of image artifacts caused by unwanted photons
US8058601B2 (en) Determining a multimodal pixon map for tomographic-image reconstruction
Ollinger Maximum-likelihood reconstruction of transmission images in emission computed tomography via the EM algorithm
EP3631763B1 (en) Method and devices for image reconstruction
Beekman et al. Selection of task-dependent diffusion filters for the post-processing of SPECT images
Chen et al. Sensor fusion in image reconstruction
Laurette et al. A three-dimensional ray-driven attenuation, scatter and geometric response correction technique for SPECT in inhomogeneous media
Selivanov et al. Cross-validation stopping rule for ML-EM reconstruction of dynamic PET series: effect on image quality and quantitative accuracy
US8014580B2 (en) Determining a pixon map for image reconstruction
Kole et al. Evaluation of the ordered subset convex algorithm for cone-beam CT
Kim et al. An effective post-filtering framework for 3-D PET image denoising based on noise and sensitivity characteristics
Qi et al. Scatter correction for positron emission mammography
Markiewicz et al. High accuracy multiple scatter modelling for 3D whole body PET
WO2021157403A1 (ja) 画像処理装置および方法、プログラム
Liacer et al. Statistically based image reconstruction for emission tomography
EP3433832B1 (en) Image reconstruction
Nunez et al. Maximum entropy and the concept of feasibility in tomographic image reconstruction
Williams et al. Elimination of loss of resolution at depth in single-photon nuclear images
Worstell et al. Diversity combining signal processing and NEC in list-mode PET

Legal Events

Date Code Title Description
LAPS Cancellation because of no payment of annual fees