JP6235564B2 - 陽電子放出断層撮影における散乱された同時発生を用いた減衰マップ - Google Patents

陽電子放出断層撮影における散乱された同時発生を用いた減衰マップ Download PDF

Info

Publication number
JP6235564B2
JP6235564B2 JP2015509532A JP2015509532A JP6235564B2 JP 6235564 B2 JP6235564 B2 JP 6235564B2 JP 2015509532 A JP2015509532 A JP 2015509532A JP 2015509532 A JP2015509532 A JP 2015509532A JP 6235564 B2 JP6235564 B2 JP 6235564B2
Authority
JP
Japan
Prior art keywords
activity
map
events
event
attenuation
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.)
Expired - Fee Related
Application number
JP2015509532A
Other languages
English (en)
Other versions
JP2015519555A (ja
JP2015519555A5 (ja
Inventor
ベルカー,ヤンニック
シュルツ,フォルクマー
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.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips NV
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 Koninklijke Philips NV filed Critical Koninklijke Philips NV
Publication of JP2015519555A publication Critical patent/JP2015519555A/ja
Publication of JP2015519555A5 publication Critical patent/JP2015519555A5/ja
Application granted granted Critical
Publication of JP6235564B2 publication Critical patent/JP6235564B2/ja
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography
    • G06T11/005Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computed tomography [CT]
    • A61B6/037Emission tomography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5205Devices using data or image processing specially adapted for radiation diagnosis involving processing of raw data to produce diagnostic data
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5258Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01TMEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/16Measuring radiation intensity
    • G01T1/161Applications in the field of nuclear medicine, e.g. in vivo counting
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01TMEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/16Measuring radiation intensity
    • G01T1/161Applications in the field of nuclear medicine, e.g. in vivo counting
    • G01T1/164Scintigraphy
    • G01T1/1641Static instruments for imaging the distribution of radioactivity in one or two dimensions using one or several scintillating elements; Radio-isotope cameras
    • G01T1/1647Processing of scintigraphic data

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • Biomedical Technology (AREA)
  • General Health & Medical Sciences (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • Molecular Biology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Pathology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biophysics (AREA)
  • Public Health (AREA)
  • Radiology & Medical Imaging (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Nuclear Medicine (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Description

本発明は、一般的には核医学画像に関する。より特定的に、本発明は、陽電子放出断層撮影(Positron Emission Tomography:PET)における減衰の見積りに関連するアプリケーションに関する。しかしながら、本発明は、また、他の利用シナリオにおけるアプリケーションに使用され、必ずしも上述のアプリケーションに限定されないことが理解されよう。
陽電子放出断層撮影(PET)においては、放射性医薬品が患者に対して施される。放射性医薬品が崩壊において陽子を放出すると、陽子は近くの電子に衝突して、消滅イベントにおいてエネルギーに係る2つの光子を反対方向に放出する。光子は、およそ511keVのエネルギーである。放出された光子は、円周状に配置された検出器によって記録される。ほとんど同時に真の同時発生イベントが測定され、どこで光子が放出されたかを判断するために使用される。真の同時発生イベントに係る放出のパス(path)は、応答系統(line of response:LOR)と呼ばれている。飛行時間(time−of−flight:TOF)PETを用いて陽子のそれぞれの衝突の時間を使用することができ、放出された陽子の出所をポイントでより良く見積もるためにLORに沿った距離を計算する。放射性医薬品の崩壊によって放出された光子の量又は数量は、活動と呼ばれる。
放射性医薬品は、癌の成長といった特定の代謝性活動に係る場所を特定するために施される。放射性医薬品は、ターゲットの代謝性活動の領域に集中する。放射性医薬品がターゲット領域において崩壊して光子が放出されると、放出された光子に基づいてターゲット領域の活動が測定され、PETシステムによって画像が再構成される。再構成された画像は、放射性医薬品の測定された活動によってコントラスト(contrast)された、特定の代謝性活動の領域を示している。
しかしながら、いくらかの放出された光子は吸収されてしまう。減衰は、検出される以前に通らなければならない材料または物質に基づいて吸収される光子を含んでいる。光子は、吸収されるか、されないかのいずれかである。より密度の濃い材料は、より大きな減衰を有し、より多くの光子を吸収する。検出器に衝突する光子の測定においては、活動と減衰の両方が主要な変数である。多くの減衰を伴う大量の活動は、測定において、より少ない減衰を伴うより小さな活動と同様に見える。
減衰補正は、いくつかのシステムにおいてX線コンピュータトモグラフィ(CT)を通じて扱われている。いくつかのシステムは、組み合わせ又はハイブリッドシステムである。CT画像の輝度は減衰と強く相関しているので、CTシステムは減衰補正に対して役に立つ。例えば、骨は、輝度において非常に明るくて濃い。骨は、より大きな減衰を有している。しかしながら、複数のスタンドアロンのPETシステムが存在しており、磁気共鳴(MR)画像とPET画像の組み合わせについて、より大きな興味がある。CTを含まないものである。
スタンドアロンのPETシステムを用いた減衰補正は、最尤推定期待値最大化法(Maximum Likelihood Expectation Maximization:MLEM)及び活動と減衰の最尤再構成(Maximum−Likelihood of Activity and Attenuation:MLAA)といったアルゴリズムを含んでいる。これらの技術は、測定された陽子放出に係る真の同時発生を使用し、活動と減衰の組み合わせの反復を使用して、活動と減衰の両方を見積るユニークな解法に集中するように努めている。これらの技術は、時々、制限され追加の外部情報を必要とする。例えば、頭といった、円周状に対称な減衰が存在する場合、真の同時発生イベントは、活動と減衰との間の相違を識別するための情報を提供しない。
いくつかの放射された光子は、偏向され、または、コンプトン(Compton)散乱される。光子が散乱されると、偏向の角度に基づいて、光子はエネルギを失う。典型的なPETシステムにおいて、散乱された光子は、単純に除外される。例えば、511keV光子の測定された衝突は、第1の光子を確立する。450keVの第2の光子に係るほとんど同時発生の衝突は、対応する、しかし散乱された光子を示している。ほとんどのPETにおいて、2つのイベントは、さらなる処理から除外される。いくつかのPETシステムにおいては、エネルギが論理的な511keVに近い場合、LORの軌跡には不確定性が存在する。イベントが511keVに十分近い場合、LORは、感度を改善するために、再構成において使用される。他には、511keVからさらに離れたエネルギをもつ散乱された光子は、捨てられる。
以降は、散乱された同時発生に基づいた新規で改善された減衰マップを開示する。上記に参照された問題又は他のものを取り扱うものである。
一つの態様に従って、画像システムは、陽電子放出断層撮影(PET)スキャナと一つまたはそれ以上のプロセッサを含んでいる。PETスキャナは、真の同時発生イベントと散乱イベントを含むイベントデータを生成し、イベントデータは、応答系統(LOR)のそれぞれのエンドポイントと、それぞれのエンドポイントのエネルギを含んでいる。一つまたはそれ以上のプロセッサは、真の同時発生イベントに基づいて、複数の活動マップと減衰マップのペアを生成し、複数の活動マップと減衰マップのペアから散乱イベントに基づいて一つの活動マップと一つの減衰マップを選択するようにプログラムされている。
別の態様に従って、画像に係る方法は、真の同時発生イベントと散乱イベントを含むイベントデータを受け取るステップを含んでおり、イベントデータは、応答系統(LOR)のそれぞれのエンドポイントと、それぞれのエンドポイントのエネルギを含んでいる。真の同時発生イベントに基づいて、複数の活動マップと減衰マップのペアが生成される。散乱イベントに基づいて、複数の活動マップと減衰マップのペアから一つの活動マップと一つの減衰マップが選択される。
別の態様に従って、画像処理システムは、陽電子放出断層撮影(PET)イベントデータのリストと、一つまたはそれ以上のプロセッサを含んでいる。PETイベントデータのリストは、複数の散乱されていない同時発生イベントと散乱された同時発生イベントを含んでおり、かつ、イベントデータは、エンドポイントとそれぞれのイベントのエネルギを含んでいる。一つまたはそれ以上のプロセッサは、散乱イベントから減衰マップを生成するようにプログラムされている。一つまたはそれ以上のプロセッサは、さらに、同時発生イベントと生成された減衰マップに基づいて、画像を再構成するようにプログラムされている。
一つの利点は、PET情報を使用して見積られた減衰マップにある。
別の利点は、既存のPETの使用の容易さを含んでいる。
別の利点は、既存のハードウェアを用いて収集され得る情報を使用することにある。
さらなる利点が、以降の詳細な説明を読んで理解することで、当業者に対して正しく理解されよう。
本発明は、種々のコンポーネントとコンポーネントの構成、および、種々のステップとステップの構成の形式であってよい。図面は、好適な実施例を説明する目的のためだけのものであり、本発明を限定するものとして理解されるべきではない。
図1は、PETスキャナにおける対称的な画像シェルに係る一つの実施例を示している。 図2は、散乱イベントに基づいて減衰を見積るPETシステムに係る一つの実施例を示している。 図3は、種々の対称的なシェルの実施例について散乱イベントの分布を模式的に示している。 図4は、散乱イベントに基づいて減衰を見積る実施例を使用した一つの方法のフローチャートである。 図5は、散乱だけに基づいて減衰を生成する実施例を使用した一つの方法のフローチャートである。
図1を参照すると、PETスキャナにおける対称的な画像シェル10に係る一つの実施例が模式的に示されている。シェルは、検出器表面12の中、または、PETスキャナのボアの中に配置されている。画像シェル10の中心には、放射性医薬品14の点光源が置かれている。放射性医薬品が放射線を発し、放射線が材料に衝突すると、消滅イベントが発生する。消滅イベントは、511keVのガンマ光子を2つ生じ、光子は反対方向に放出される。真の同時発生イベントとは、光子が検出器表面12に向かって反対方向に妨げられることなく進行し、2つのエンドポイント16、18を確立するものである。真のイベントに対する応答系統は、エンドポイント間の直線である。
均一な密度と厚みのシェルが存在すると、減衰と散乱が生じる。いくつかの光子が吸収され、検出器への衝突回数または記録イベントの回数が減少する。活動は、光源によって光子が放射される際に測定される崩壊の量である。減衰は、イベントの記録された量を減じるように動作する。光子が検出器に到達しないからである。対称的な画像シェルに対して、記録された光子フラックス(flux)は、半径状に対称であり、A・e(−μ・d)に比例している。ここで、Aは点光源における活動であり、μは線形減衰係数であり、dはシェルの厚さである。不感時間効果(dead time effect)のために、真の同時発生レートは、シェルの半径と共にわずかに変化する。変化は、例えば、単調であり、かつ、統計的な変動のレベルを顕著に上回るものである。この項は、シェルの半径を含んでいないので、シェルを、PETシステムのボアの中のどこにでも配置することができる。そして、真(511keV)のイベントだけを用いて、活動と減衰の所定のペアを区別することができる。例えば、シェルの内側のあらゆる半径に対して、そして、任意のスケールパラメータα>0に対して、(A、μ、d)と(A、μ*α、d/α)の分布を区別することができない。より大まかに言えば、同一の減衰パスの積である∫μ(s)ds、sはLORに沿ったパラメータであって全ての測定されたLORに沿った(A、μ、d)としてのもの、を生じている全ての減衰分布は、認識することができない。この状況は、減衰自己クロストークと呼ばれる。さらに、(A、μ、d)と(A*e((αー1)*μ*d、μ*α、d)も認識可能ではない。後者は、活動と減衰との間のクロストークと呼ばれるもので、医療事例において理解され観察されてきており、活動または減衰のいずれかにおいて定数の既知でないオフセットについて解かれている。減衰自己クロストーク及び活動と減衰との間のクロストークは、散乱イベントに係る量、分布、および、特性を使用して解かれる。
散乱イベント22は、光子が、シェルといった材料に衝突して、エネルギ損失を伴って角度24で偏向された場合に生じる。オリジナルの放出パス18に沿ってエンドポイントに衝突する代わりに、光子は、偏向されて異なるエンドポイント26に衝突する。記録されたエンドポイント16、26は、応答系統(LOR)28を形成する。LOR角度30は、オリジナルの進行ライン20と記録されたエンドポイントとの間の角度として算出される。イベントは、散乱されたものとして知られている。変更された光子のエンドポイント26において記録されたエネルギが、508keVまたはそれ以下といった、511keVより小さいからである。
図2は、散乱イベントに基づいて減衰を見積るPETシステム36に係る一つの実施例を模式的に示している。システムは、PETスキャナ38を含んでいる。PETスキャナ38は、患者を受け入れる円環体形状の開口部40またはボアを有している。患者は、ベットまたはカウチといった患者サポート43によって支持されている。放射性医薬品が患者に施され、患者は、スキャナの開口部40またはボアの中に置かれる。開口部またはボアには、患者から放出される光子を検出する検出器42が並んでいる。PETスキャナは、同時発生プロセッサ44を含んでいる。それぞれの光子が検出器に衝突すると、同時発生プロセッサ44は、衝突の時間、エネルギ、そして場所をメモリの中に記録する。記録された衝突のリストは、メモリの中に、例えば、リストモードで、保管される。同時発生プロセッサは、散乱イベントと散乱されていないイベントがを特定するために、メモリに対して時間及びエネルギのゲーティングまたはフィルタリングを提供することができる。散乱イベントは、メモリの中に含まれている。例えば、190keVから509keVの間のエネルギをもったイベントを含んでおり、およそ511keVの別の衝突とペアにされたイベントは、散乱されたLORを含んでいる。エネルギ感度または検出器によるエネルギ損失は、ゲーティングパラメータにおいて調節される。別のゲーティングパラメータは、4.2°より大きいか等しいといった、LOR角度を含み得る。スキャナとメモリは、イベントデータの保管と検索のためにネットワーク48に接続されている。ストレージは、ソフトウェアストレージおよび画像ストレージを含んでよい。
システムは、ワークステーション50を含んでいる。ワークステーション50は、ネットワーク48に接続されている。ワークステーション50は、電子プロセッサまたは電子プロセッサデバイス52、少なくとも一つの入力デバイス54、および、マップ、画像、メニュー、そして、ユーザコントロールを表示するディスプレイ56を含んでいる。ワークステーション50は、コンピュータアレイ、アレイプロセッサ、デスクトップコンピュータ、ラップトップ、タブレット、モバイルコンピューティングデバイス、スマートフォン、および類似のもの、を含み得る。入力デバイスはキーボード、マウス、マイクロフォン、および類似のもの、であってよい。
活動と減衰マップ発生器58は、イベントデータを入力して、活動マップと減衰マップのペアを生成する。ペアは、イベントデータに基づいた類似の分布を含んでよい。活動と減衰マップ発生器は、真の散乱されていない同時発生ペアに対するイベントデータをフィルタして、マップのペアを生成する。アルゴリズムは、MLAA、MLEM、等を含み得る。
代替的な実施例において、活動と減衰マップ発生器58は、散乱だけから直接的に減衰マップを生成する。それぞれのLORに対して、クライン−ニシナの公式(Klein−Nishina formula)または、XCOMといった散乱クロスセクションデータベースを使用して、散乱された光子のエネルギに基づいてコンプトン角度が算出される。LORおよびコンプトン角度のエンドポイントから、可能な散乱ポイントに係る表面が計算される。エネルギが既知である限り、それぞれの可能な散乱ポイントは、LORエンドポイントを通る円弧のファミリ上に置かれる。エネルギ値におけるいくらかの不確実性を伴って、可能な散乱ポイントの表面は、可能な散乱ポイントのボリュームになっている。幾何学的な考察と他の先天的な情報からそれぞれの散乱ポイントの可能性が計算される。逆投射のステップにおいて、全てのLORに係る貢献の可能性が合計される。それぞれのLORは、再び、クラインーニシナの公式と幾何学的情報に基づいて、重み付けされ得る。任意的なフィルタリングは、減衰マップを生成するために、フィルタされた逆投射(FBP)におけるアプローチといったものをフォローすることができる。散乱イベントと同様に散乱されていないイベントを使用することも、また、考えられる。
シミュレータ60は、複数のマップペアを入力し、それぞれのマップペアを生成することができた散乱イベント(シミュレーションされたリストモードデータ)の分布をシミュレーションする。散乱イベントのシミュレーションされた分布は、比較器62によって、リストモードメモリからの散乱イベントに係る実際の分布と比較される。比較器は、代替的に、点光源としてのそれぞれのボクセルといったものにより、分布を組成コンポーネントに分解することができる。分布は、最も合う活動と減衰マップのペアを選択するために、最小二乗法に基づいて実際の分布と比較され得る。散乱イベントに基づいた活動と減衰マップのペアの選択は、類似の分布を含んでいる生成されたペア間のクロストークまたは自己クロストークを解決する。
再構成ユニット64は、選択された減衰マップとイベントとデータに基づいて、一つまたはそれ以上のマップを再構成する。散乱だけから減衰マップを直接に生成する実施例において、生成された減衰マップは選択された減衰マップである。イベントデータは、散乱イベントを取り除くようにフィルタされ、真の同時発生LORを再構成において入力する。選択された減衰マップは、再構成の最中に減衰補正を提供する。代替的に、活動マップは、患者の中の放射性医薬品または代謝性活動の存在を表わす一つまたはそれ以上の画像として再構成され得る。
同時発生プロセッサ44、活動と減衰マップ生成器58、シミュレータ60、比較器62、および、再構成ユニット64は、一つまたはそれ以上の電子データ処理デバイス、または、ネットワーク48によりワークステーション50と動作可能に接続されたネットワークベースのサーバーコンピュータ、または、類似のものによって、好適に具現化される。さらに、開示された活動と減衰マップ生成および選択技術は、電子データ処理デバイスによって読取り可能であり、電子データ処理デバイスによって実行可能なインストラクション(例えば、ソフトウェア)を保管している固定のストレージ媒体を使用して好適に実施され、開示された活動と減衰マップの生成および選択技術を実行する。
「コンピュータで読取り可能なストレージ媒体」は、ここにおいて使用されるように、コンピューティングデバイスのプロセッサによって実行可能なインストラクションを保管することができるあらゆる有形のストレージ媒体を含んでいる。コンピュータで読取り可能なストレージ媒体は、コンピュータで読取り可能な固定ストレージ媒体として参照されてよい。コンピュータで読取り可能なストレージ媒体は、また、固定のコンピュータで読取り可能な媒体として参照されてもよい。いくつかの実施例において、コンピュータで読取り可能なストレージ媒体は、また、コンピューティングデバイスのプロセッサによってアクセスし得るデータを保管することもできる。コンピュータで読取り可能なストレージ媒体の実施例は、これらに限定されるわけではないが、以下のものを含んでいる。フロッピー(登録商標)ディスク、磁気ハードディスクドライブ、半導体ハードドライブ、フラッシュメモリ、USBサムドライブ、ランダムアクセスメモリ(RAM)、読出し専用メモリ(ROM)、光ディスク、光磁気ディスク、および、プロセッサのレジスタファイル、である。光ディスクの実施例は、コンパクトディスク(CD)およびデジタル多目的ディスク(DVD)を含んでいる。例えば、CD−ROM、CD−RW、CD−R、DVD−ROM、DVD−RW、または、DVD−Rディスクである。コンピュータで読取り可能なストレージ媒体という用語は、また、ネットワークまたは通信リンクを介してコンピュータデバイスによってアクセスされ得る種々のタイプの記録媒体も参照する。例えば、データは、モデム、インターネット、または、ローカルエリアネットワークにわたり、リトリーブされてよい。コンピュータで読取り可能なストレージ媒体に対する参照は、おそらく複数のコンピュータで読取り可能なストレージ媒体として解釈されるべきである。プログラムに係る種々の実行可能なコンポーネントは、異なる場所に保管されてよい。コンピュータで読取り可能なストレージ媒体は、例えば、同一のコンピュータシステムの中の複数のコンピュータで読取り可能なストレージ媒体であってよい。コンピュータで読取り可能なストレージ媒体は、また、複数のコンピュータシステムまたはコンピューティングデバイスの中に分散されたコンピュータで読取り可能なストレージ媒体であってよい。
「メモリ」は、コンピュータで読取り可能なストレージ媒体の一つの実施例である。コンピュータメモリは、プロセッサに対して直接的にアクセス可能なあらゆるメモリである。コンピュータメモリの実施例は、これらに限定されるわけではないが、RAMメモリ、レジスタ、および、レジスタファイル、を含んでいる。「コンピュータメモリ」または「メモリ」に対する参照は、おそらく複数のメモリであるものとして解釈されるべきである。メモリは、例えば、同一のコンピュータシステムの中の複数のメモリであってよい。メモリは、また、複数のコンピュータシステムまたはコンピューティングデバイスの中に分散された複数のメモリであってよい。
「コンピュータストレージ」または「ストレージ」は、コンピュータで読取り可能なストレージ媒体の一つの実施例である。コンピュータストレージは、あらゆる不揮発性のコンピュータで読取り可能なストレージ媒体であってよい。コンピュータストレージの実施例は、これらに限定されるわけではないが、以下のものを含んでいる。ハードディスクドライブ、USBサムドライブ、フロッピー(登録商標)ディスクドライブ、スマートカード、DVD、CD−ROM、および、半導体ハードドライブ、である。いくつかの実施例において、コンピュータストレージは、また、コンピュータメモリであってよく、または、その逆も同様である。「コンピュータストレージ」または「ストレージ」に対する参照は、おそらく複数のストレージであるものとして解釈されるべきである。ストレージは、例えば、同一のコンピュータシステムの中の複数のストレージであってよい。ストレージは、また、複数のコンピュータシステムまたはコンピューティングデバイスの中に分散された複数のストレージであってよい。
「プロセッサ」は、ここにおいて使用されるように、プログラムまたはマシンで実行可能なインストラクションを実行することができる電子的コンポーネントを含んでいる。「プロセッサ」を含んでいるコンピューティングデバイスに対する参照は、おそらく一つまたはそれ以上のプロセッサまたはプロセッサコアを含んでいるものとして解釈されるべきである。プロセッサは、例えば、マルチコアプロセッサであってよい。プロセッサは、また、複数のコンピュータシステムの中の単一のコンピュータシステムまたは分散コンピュータシステムの集合を参照してよい。コンピューティングデバイスという用語は、それぞれがプロセッサを含んでいるコンピューティングデバイスの集合またはネットワークをおそらく参照するものと解釈されてもよい。多くのプログラムは、同一のコンピューティングデバイスまたは複数のコンピューティングデバイスにわたり分散されたものの中にある複数のプロセッサに、インスタンスを実行させる。
「ディスプレイ」または「ディスプレイ装置」は、ここにおいて使用されるように、画像またはデータを表示するように適合された出力デバイスまたはユーザインターフェイスを含んでいる。ディスプレイは、視覚、聴覚、及び/又は、触覚データを出力し得る。ディスプレイの実施例は、これらに限定されるわけではないが、以下のものを含んでいる。コンピュータモニタ、テレビ画面、タッチ画面、触覚電子ディスプレイ、点字スクリーン、ブラウン管(CRT)、ストレージ管、双安定ディスプレイ、電子ペーパー、ベクトルディスプレイ、フラットパネルディスプレイ、蛍光表示管(VF)、発光ダイオード(LED)ディスプレイ、エレクトロルミネッセントディスプレイ(ELD)、プラズマディスプレイパネル(PDP)、液晶ディスプレイ(LCD)、有機発光ダイオードディスプレイ(OLED)、プロジェクタ、および、ヘッドマウントディスプレイ、である。
図3は、種々の対称的なシェルの実施例について散乱イベントの分布を模式的に示している。それぞれのヒストグラムは、5cmから40cmまで5cmのステップの半径をもつ約10cmの均一な厚みであるシェルに対する散乱イベントの分布である。真の同時発生イベントに関して、活動と減衰マップは、シェルサイズ間の相違を識別しない。しかしながら、散乱イベントに関して、分布は、複数の可能性のある活動と減衰マップのペアの中で、シェルサイズを判断または区別するために使用することができる。
ヒストグラムの分布における相違は、散乱イベントがどのように自己クロストークを解決するための追加の情報を提供するかを説明している。類似の代替的な構成は、点光源の活動および全てのLORに沿った経路積分といった減衰の両方を変化させることによって、活動と減衰のクロストークをどのように解決するかを説明することができる。第1のヒストグラムは、小さなシェル(r=5cm)に対する散乱イベントの分布を含んでいる。分布は、120秒といった同一の取得時間に対して収集された複数のLOR角度にわたり記録されたエネルギレベルに基づいて、示されている。小さなシェルは、LOR角度の増加につれてより緩やかになるスロープを伴う、より幅広い分布を示している。シェルサイズが増加すると、分布は、LOR角度の増加につれてより険しくなるスロープを伴って、幅がより狭くなる。散乱イベントの実際の分布に基づいて、分布は、真の同時発生イベントに基づき生成された活動と減衰マップのペア間のクロストークおよび自己クロストークの曖昧性を解決するための情報を提供する。
図4は、散乱イベントに基づいて減衰を見積る実施例を使用した一つの方法のフローチャートである。ステップ70において、イベントデータが取得される。患者は、放射性医薬品を施され、患者サポート上でスキャナのボアまたは開口部の中に配置される。放射性医薬品が崩壊する際に、陽子が放出され、そして、陽子は、消滅イベントを引き起こして、光子ペアが放出される。放出された光子は、衝突の際に、検出器によって記録される。検出器が衝突を検出すると、イベントペアとそれぞれのイベントのエネルギが、リストモードのメモリの中に、同時発生プロセッサによって記録される。同時発生プロセッサは、リストモードのメモリの中に、散乱イベントと真の同時発生イベントを含んでいる。
活動と減衰マップ生成器は、リストモードのメモリからリストモードデータを受け取り、ステップ72において、複数の活動マップと減衰マップのペアを生成する。ペアは、MLAA、MLEM、等といった、アルゴリズムを使用して生成され得る。アルゴリズムは、リストモードのメモリからフィルタされた真の同時発生イベントを使用する。
ステップ74において、シミュレータは、複数の活動と減衰マップのペアに基づいて、散乱イベントに係る一つまたはそれ以上の分布をシミュレーションする(例えば、モンテカルロシミュレーション)。存在する場合、飛行時間(TOF)は追加の情報を提供することができる。情報は、散乱イベントの分布からの情報を補足することができ、または、シミュレーションの量を削減するために使用される。シミュレーションされた散乱イベントの分布は、リストモードデータからフィルタされた散乱イベントの実際の分布と比較またはフィットされ、ステップ76において、最も合う活動マップと減衰マップが選択される。選択プロセスは、クロストークおよび自己クロストークの問題を解決し、生成された活動と減衰マップのペアに対する収束性の解決策を提供する。
ステップ78において、再構成ユニットによって一つまたはそれ以上の画像が再構成される。再構成ユニットは、一つまたはそれ以上の画像を再構成するために、フィルタされた真の同時発生イベントと選択された減衰マップを使用する。代替的に、再構成ユニットは、選択された活動と減衰マップのペアに係る活動マップから再構成することができる。一つまたはそれ以上の画像は、ワークステーションのディスプレイ56または他のディスプレイ上に表示される。
図5を参照すると、別の実施例においては、散乱だけに基づいて減衰マップが生成される。ステップ80において、それぞれのLORに対して、クライン−ニシナの公式、または、XCOMといった散乱されたクロスセクションのデータベースを使用して、散乱された光子のエネルギからコンプトン角度が計算される。ステップ82において、LORのエンドポイントとコンプトン角度から、可能な散乱ポイントに係る表面が計算される。エネルギが既知である限り、それぞれの可能な散乱ポイントは、LORエンドポイントを通る円弧のファミリ上に置かれる。エネルギ値におけるいくらかの不確実性を伴って、可能な散乱ポイントの表面は、可能な散乱ポイントのボリュームになっている。ステップ84において、幾何学的な考察と他の先天的な情報からそれぞれの散乱ポイントの可能性が計算される。それぞれのLORは、再び、クラインーニシナの公式と幾何学的情報に基づいて、重み付けされ得る。ステップ86において、フィルタされた逆投射(FBP)におけるアプローチといった、それぞれのLORに対する可能な分布に対して、フィルタリングと逆投射が適用されて、減衰マップと画像が生成される。散乱と同様に、散乱されていないイベントの使用も考えられる。
全てのイベントからの減衰補正された画像再構成は、散乱補正及び類似のものを含んでいる、真のもの及び光電ピーク散乱からの標準の画像再構成を含み得るが、追加の散乱イベントを合体してもよい。
ここにおいて提示された説明的な実施例に関して、所定の構成及び/又は機能の特徴が、定められたエレメント及び/又はコンポーネントの中に組み入れられているものとして記述されていることが正しく理解されるべきである。しかしながら、これらの特徴は、また、同一または類似の利点に対して、好適な他のエレメント及び/又はコンポーネントの中に同様に組み込まれているものと考えられる。所望のアプリケーションに対して適切な他の代替的な実施例を達成するために、典型的な実施例に係る異なる態様が、好適であるとして選択的に使用され得ることも正しく理解されるべきである。態様に係るそれぞれの利点をそれによって実現する他の代替的な実施例は、その中に組み込まれている。
ここにおいて説明された所定のエレメントまたはコンポーネントは、ハードウェア、ソフトウェア、ファームウェア、または、それらの組み合わせを介して、その機能性を好適に実施させ得ることも、正しく理解されるべきである。加えて、一緒に組み込まれているものとしてここにおいて説明された所定のエレメントは、好適な環境下において、スタンドアロンなエレメントであるか、そうでなければ、分割されることが正しく理解されるべきである。同様に、一つの所定のエレメントによって実行されるものとして説明された複数の所定の機能性は、個々の機能性を個別に実行するように動作する複数の別個のエレメントによって実行されてよい。または、所定の個々の機能性は、分割され、協調している複数の別個のエレメントによって実行されてよい。代替的に、そうでなければ、ここにおいて、お互いに別個のものとして説明及び/又は示される、いくつかのエレメントまたはコンポーネントは、物理的または機能的に好適に組み合わされてよい。
つまり、本発明の仕様は、好ましい実施例に関して明らかにされてきた。本発明の明細書を読んで理解すれば、他人に対して変形および代替が生じるのが明らかである。全てのこうした変形および代替は、それが添付の特許請求の範囲内にあるか、それらの均等物である限り、本発明に含まれるものとして理解されるべきである。つまり、上記に開示された種々のもの、他の特徴と機能性、または、それらの代替は、多くの他の異なるシステムまたはアプリケーションの中に望ましくは結合されてよいことが正しく理解されるべきである。そして、また、種々の現在では予見できない、もしくは、予期しない代替、変形、変化、または、改善が、当業者によって後に行われ、同様に、以降の特許請求の範囲によって包含されるように意図されていることも、正しく理解されるべきである。

Claims (15)

  1. 真の同時発生イベントと散乱イベントを含んでいるイベントデータを生成する陽電子放出断層撮影(PET)スキャナであり、前記イベントデータは、応答系統(LOR)のそれぞれのエンドポイントとそれぞれのエンドポイントのエネルギを含む、陽電子放出断層撮影スキャナと、
    一つまたはそれ以上のプロセッサであり、
    前記真の同時発生イベントに基づいて、複数の活動マップと減衰マップのペアを生成し、かつ、
    前記散乱イベントに基づいて、前記複数の活動と減衰マップのペアから一つの活動マップと一つの減衰マップを選択する、
    ようにプログラムされているプロセッサと、
    を含む、画像システム。
  2. 前記一つまたはそれ以上のプロセッサは、さらに、
    前記活動と減衰マップのペアのそれぞれに基づいて、少なくとも前記散乱イベントの分布をシミュレーションし、かつ、
    散乱イベントの前記シミュレーションされた分布と記録された散乱イベントから構成された散乱イベントの記録された分布との間の比較に基づいて、前記活動と減衰マップのペアを選択する、ようにプログラムされている、
    請求項1に記載の画像システム。
  3. 前記一つまたはそれ以上のプロセッサは、さらに、
    前記エンドポイントの位置に基づいて、散乱イベントの分布を減衰コンポーネントへと分割する、ようにプログラムされている、
    請求項1または2に記載の画像システム。
  4. 前記一つまたはそれ以上のプロセッサは、さらに、
    前記真の同時発生イベントと前記選択された減衰マップ、または、前記選択された活動マップ、のうち少なくとも一つに基づいて、画像を再構成する、ようにプログラムされている、
    請求項1乃至3いずれか一項に記載の画像システム。
  5. 前記複数の活動マップと減衰マップのペアは、類似の分布である、
    請求項1乃至4いずれか一項に記載の画像システム。
  6. 前記分割は、活動の点光源としてのボクセルを含んでいる、
    請求項3に記載の画像システム。
  7. 前記散乱イベントの分布は、計算された偏向角度を含んでいる、
    請求項2乃至6いずれか一項に記載の画像システム。
  8. 前記散乱イベントの分布は、偏向角度のオフセットの発生頻度を含んでいる、
    請求項2乃至7いずれか一項に記載の画像システム。
  9. 散乱イベントは、190keVと509keVとの間のエネルギを含んでおり、応答系統(LOR)の偏向は、4.2°以上である、
    請求項1乃至8いずれか一項に記載の画像システム。
  10. 真の同時発生イベントと散乱イベントを含んでいるイベントデータを受け取るステップであり、前記イベントデータは、応答系統(LOR)のそれぞれのエンドポイントとそれぞれのエンドポイントのエネルギを含む、ステップと、
    前記真の同時発生イベントに基づいて、複数の活動マップと減衰マップのペアを生成するステップと、
    前記散乱イベントに基づいて、前記複数の活動と減衰マップのペアから一つの活動マップと一つの減衰マップを選択するステップと、
    を含む、画像表示の方法。
  11. 前記方法は、さらに、
    前記活動と減衰マップのペアのそれぞれに基づいて、少なくとも前記散乱イベントの分布をシミュレーションするステップと、
    散乱イベントの前記シミュレーションされた分布と記録された散乱イベントから構成された散乱イベントの記録された分布との間の比較に基づいて、前記活動と減衰マップのペアを選択するステップと、を含む、
    請求項10に記載の方法。
  12. 前記方法は、さらに、
    前記エンドポイントの位置に基づいて、散乱イベントの分布を減衰コンポーネントへと分割するステップと、を含む、
    請求項10または11に記載の方法。
  13. 前記方法は、さらに、
    前記真の同時発生イベントと前記選択された減衰マップに基づいて、画像を再構成するステップと、を含む、
    請求項10乃至12いずれか一項に記載の方法。
  14. 散乱イベントの分布は、計算された偏向角度を含んでいる、
    請求項10乃至13いずれか一項に記載の方法。
  15. 請求項10乃至14いずれか一項に記載の方法を実行するように、一つまたはそれ以上の電子データ処理デバイスをコントロールするソフトウェア、を運んでいる固定のコンピュータで読取り可能なストレージ媒体。
JP2015509532A 2012-05-04 2013-04-23 陽電子放出断層撮影における散乱された同時発生を用いた減衰マップ Expired - Fee Related JP6235564B2 (ja)

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
US201261642730P 2012-05-04 2012-05-04
US61/642,730 2012-05-04
US201261661455P 2012-06-19 2012-06-19
US61/661,455 2012-06-19
PCT/IB2013/053181 WO2013164731A1 (en) 2012-05-04 2013-04-23 Attenuation map with scattered coincidences in positron emission tomography

Publications (3)

Publication Number Publication Date
JP2015519555A JP2015519555A (ja) 2015-07-09
JP2015519555A5 JP2015519555A5 (ja) 2017-06-22
JP6235564B2 true JP6235564B2 (ja) 2017-11-22

Family

ID=48539336

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2015509532A Expired - Fee Related JP6235564B2 (ja) 2012-05-04 2013-04-23 陽電子放出断層撮影における散乱された同時発生を用いた減衰マップ

Country Status (6)

Country Link
US (1) US9734600B2 (ja)
EP (1) EP2844148B1 (ja)
JP (1) JP6235564B2 (ja)
CN (1) CN104271045B (ja)
RU (1) RU2014148796A (ja)
WO (1) WO2013164731A1 (ja)

Families Citing this family (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2620862C2 (ru) * 2012-05-21 2017-05-30 Конинклейке Филипс Н.В. Быстрая оценка рассеяния при реконструкции посредством позитронно-эмиссионной томографии
WO2014074666A1 (en) * 2012-11-07 2014-05-15 Massachusetts Institute Of Technology Inter-detector scatter enhanced emission tomography
US11058371B2 (en) * 2013-10-30 2021-07-13 The Board Of Trustees Of The Leland Stanford Junior University Simultaneous attenuation and activity reconstruction for positron emission tomography
CN107250844B (zh) 2014-12-12 2020-06-09 爱丁堡大学董事会 用于处理正电子发射断层摄影数据的方法和装置
CN107635469B (zh) * 2015-05-19 2021-04-20 皇家飞利浦有限公司 基于pet系统中的散射符合的衰减图的估计
WO2016197127A1 (en) * 2015-06-04 2016-12-08 Rensselaer Polytechnic Institute Attenuation map reconstruction from tof pet data
CN105193442B (zh) * 2015-09-17 2018-12-25 沈阳东软医疗系统有限公司 一种pet扫描时间的确定方法和装置
JP6516019B2 (ja) * 2015-11-19 2019-05-22 株式会社島津製作所 放射線断層撮影装置
US11354832B2 (en) * 2017-05-01 2022-06-07 Koninklijke Philips N.V. Generation of accurate hybrid datasets for quantitative molecular imaging
WO2018202878A1 (en) * 2017-05-04 2018-11-08 Koninklijke Philips N.V. Positron emission tomography (pet) timing calibration using coincidences involving high-energy cascade gamma from prompt-gamma positron emitters
JP6761610B2 (ja) * 2017-05-29 2020-09-30 株式会社島津製作所 吸収係数画像推定方法、吸収係数画像推定プログラム並びにそれを搭載したポジトロンct装置
US10438378B2 (en) 2017-08-25 2019-10-08 Uih America, Inc. System and method for determining an activity map and an attenuation map
JP6932250B2 (ja) * 2017-09-14 2021-09-08 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. 陽電子放出断層撮影における狭いエネルギー窓カウントから再構成された放出画像推定を使用した散乱補正
CN111588399B (zh) * 2020-05-27 2023-05-16 上海联影医疗科技股份有限公司 医学成像设备状态监控的方法、设备和计算机设备
WO2022015834A1 (en) * 2020-07-14 2022-01-20 Yale University Ct-free attenuation correction for spect using deep learning with imaging and non-imaging information
US20220130079A1 (en) * 2020-10-23 2022-04-28 Siemens Medical Solutions Usa, Inc. Systems and methods for simultaneous attenuation correction, scatter correction, and de-noising of low-dose pet images with a neural network

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6310968B1 (en) * 1998-11-24 2001-10-30 Picker International, Inc. Source-assisted attenuation correction for emission computed tomography
JP4659962B2 (ja) * 2000-10-04 2011-03-30 株式会社東芝 核医学診断装置
JP3925278B2 (ja) * 2002-04-04 2007-06-06 株式会社島津製作所 核医学イメージング装置
US7414246B2 (en) * 2006-01-03 2008-08-19 Koninklijke Philips Electronics N.V. Achieving accurate time-of-flight calibrations with a stationary coincidence point source
RU2010123369A (ru) * 2007-11-09 2011-12-20 Конинклейке Филипс Электроникс Н.В. (Nl) Устройство и способ создания карты ослаблений
JP4836094B2 (ja) * 2008-06-18 2011-12-14 株式会社日立製作所 核医学撮像装置およびその初期散乱位置判定方法
US8218848B2 (en) 2008-07-23 2012-07-10 Siemens Aktiengesellschaft System and method for the generation of attenuation correction maps from MR images
US9453922B2 (en) * 2008-09-09 2016-09-27 Multi-Magnetics Incorporated System and method for correcting attenuation in hybrid medical imaging
US8384036B2 (en) * 2008-09-11 2013-02-26 Siemens Medical Solutions Usa, Inc. Positron emission tomography (PET) imaging using scattered and unscattered photons
US8017914B2 (en) 2008-11-13 2011-09-13 General Electric Company Method and system for scatter correction
EP2433261B1 (en) 2009-05-20 2014-03-05 Koninklijke Philips N.V. Continuous time-of-flight scatter simulation method
US8299438B2 (en) * 2009-07-16 2012-10-30 Siemens Medical Solutions Usa, Inc. Model based estimation of a complete or partial positron emission tomography attenuation map using maximum likelihood expectation maximization
US8620053B2 (en) 2009-11-04 2013-12-31 Siemens Medical Solutions Usa, Inc. Completion of truncated attenuation maps using maximum likelihood estimation of attenuation and activity (MLAA)
WO2014205423A1 (en) * 2013-06-21 2014-12-24 O.N.Diagnostics, LLC Quantitative phantomless calibration of computed tomography scans

Also Published As

Publication number Publication date
EP2844148A1 (en) 2015-03-11
JP2015519555A (ja) 2015-07-09
CN104271045A (zh) 2015-01-07
US9734600B2 (en) 2017-08-15
CN104271045B (zh) 2018-10-02
US20150098640A1 (en) 2015-04-09
EP2844148B1 (en) 2019-10-09
WO2013164731A1 (en) 2013-11-07
RU2014148796A (ru) 2016-06-27

Similar Documents

Publication Publication Date Title
JP6235564B2 (ja) 陽電子放出断層撮影における散乱された同時発生を用いた減衰マップ
RU2620862C2 (ru) Быстрая оценка рассеяния при реконструкции посредством позитронно-эмиссионной томографии
Elschot et al. Quantitative monte carlo–based 90y spect reconstruction
EP2984631B1 (en) Method for modeling and accounting for cascade gammas in images
CN106659452B (zh) 在定量单光子发射计算机断层扫描中利用多个光电峰的重构
WO2015125086A1 (en) In-reconstruction filtering for positron emission tomography(pet) list mode iterative reconstruction
CN110023999A (zh) 发射和透射断层摄影中的交互式目标超快重建
Siman et al. Effects of image noise, respiratory motion, and motion compensation on 3D activity quantification in count-limited PET images
US10772580B2 (en) Multiple emission energies in single photon emission computed tomography
D′ Ambrosio et al. Attenuation correction for small animal PET images: a comparison of two methods
Walrand et al. The origin and reduction of spurious extrahepatic counts observed in 90Y non-TOF PET imaging post radioembolization
US9592408B2 (en) Dose-volume kernel generation
Nuyts Nuclear medicine technology and techniques
Schmidtlein et al. Using an external gating signal to estimate noise in PET with an emphasis on tracer avid tumors
Wei et al. Evaluation of quantitative, efficient image reconstruction for VersaPET, a compact PET system
Pryma Nuclear medicine: practical physics, artifacts, and pitfalls
Martins et al. A direct time-of-flight reconstruction for whole-body single-bed RPC-PET: Results from lesion and anthropomorphic simulated data
Theodorakis et al. Prospective PET image quality gain calculation method by optimizing detector parameters
Moreau et al. Assessment of a fully 3D Monte Carlo reconstruction method for preclinical PET with iodine-124
Juma Lesion Synthesis Toolbox: Development and validation of dedicated software for synthesis of lesions in raw PET and CT patient images
Gabrani-Juma Lesion Synthesis Toolbox: Development and validation of dedicated software for synthesis of realistic lesions in raw PET and CT patient images
Anderl Selective internal radiotherapy-determination of the 90Y biodistribution using computer simulation of PET
Kalaitzidis From Monte Carlo PET Simulations to Reconstructed Images: Modelling and Optimisation for 68Ga Theragnostics
Hicham Single Photon Emission Computed Tomography (SPECT) imaging systems optimization in nuclear medicine: Monte Carlo study
Mille A study of shape-dependent partial volume correction in PET imaging using ellipsoidal phantoms fabricated via rapid prototyping

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20160420

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20170515

A871 Explanation of circumstances concerning accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A871

Effective date: 20170515

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20170516

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20170614

AA91 Notification of revocation by ex officio

Free format text: JAPANESE INTERMEDIATE CODE: A971091

Effective date: 20170808

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20170823

TRDD Decision of grant or rejection written
A975 Report on accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A971005

Effective date: 20170829

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20171003

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20171026

R150 Certificate of patent or registration of utility model

Ref document number: 6235564

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

LAPS Cancellation because of no payment of annual fees