[論文レビュー] A statistical reconstruction algorithm for positronium lifetime imaging using time-of-flight positron emission tomography

Hsin‐Hsiung Huang, Zheyuan Zhu+4 more|arXiv (Cornell University)|Jun 13, 2022

Medical Imaging Techniques and Applications被引用数 5

ひとこと要約

本論文は、時間飛行（TOF）PETにおけるポジトロニウム寿命イメージング（PLI）のための統計的再構成アルゴリズムを提案する。この手法は、取り込み量と寿命分布の同時尤度をモデル化することで、有限な時間分解能下でも高精度で定量的な寿命イメージの回復を可能にする。L-BFGS-B最適化と現実的なTOF分解能モデルを用いた最大尤度フレームワークにより、有限な時間分解能条件下でも、模擬TOF-PETデータから高忠実度の寿命イメージを再構成するのに成功した。

ABSTRACT

Positron emission tomography (PET) is an important modality for diagnosing diseases such as cancer and Alzheimer's disease, capable of revealing the uptake of radiolabeled molecules that target specific pathological markers of the diseases. Recently, positronium lifetime imaging (PLI) that adds to traditional PET the ability to explore properties of the tissue microenvironment beyond tracer uptake has been demonstrated with time-of-flight (TOF) PET and the use of non-pure positron emitters. However, achieving accurate reconstruction of lifetime images from data acquired by systems having a finite TOF resolution still presents a challenge. This paper focuses on the two-dimensional PLI, introducing a maximum likelihood estimation (MLE) method that employs an exponentially modified Gaussian (EMG) probability distribution that describes the positronium lifetime data produced by TOF PET. We evaluate the performance of our EMG-based MLE method against \st{traditional} approaches using exponential likelihood functions and penalized surrogate methods. Results from computer-simulated data reveal that the proposed EMG-MLE method can yield quantitatively accurate lifetime images. We also demonstrate that the proposed MLE formulation can be extended to handle PLI data containing multiple positron populations.

研究の動機と目的

• 実用的なTOF-PETシステムにおける有限な時間分解能を考慮したポジトロニウム寿命イメージング（PLI）のための包括的な統計的モデルの欠如に対処すること。
• リストモードTOF-PETデータから同時に取り込み量と寿命イメージを推定できる、共同最大尤度再構成フレームワークを開発すること。
• 有限な時間分解能および検出器応答の影響を統計モデルに組み込み、再構成過程での情報損失を回避すること。
• tracerの取り込み機構に依存しない正確な寿命分布の回復を可能にすることで、定量的PLIを実現すること。
• 現実的なシステムパラメータを用いたコンピュータシミュレーションによるTOF-PETデータを用いて、アルゴリズムの性能を検証すること。

提案手法

• 有限な時間分解能と直交ポジトロンウム（o-Ps）寿命の指数的減衰を考慮した2次元リストモードTOF-PETデータの統計的モデルを構築する。
• 検出器応答関数と時間飛行の不確実性を組み込んだ、取り込み量および寿命パラメータの共同対数尤度関数を導出する。
• 取り込み量および寿命イメージの両方に正の制約を課した上で、L-BFGS-B最適化アルゴリズムを用いて尤度関数を最大化する。
• 逆ラプラス変換問題を安定的に解けるように、ペナルティ付きサロゲート（PS）法を適用する。
• 寿命パラメータに関する尤度関数の解析的勾配を導出し、効率的な最適化を可能にする。
• 現実的なTOF分解能（200–600 ps FWHM）を想定したデータをシミュレートし、プロンプトカウント（PC）および事前設定時間（PT）の両取得モードで画像を再構成することで、モデルを検証する。

実験結果

リサーチクエスチョン

• RQ1 有限な時間分解能を正確に反映する統計的モデルを開発することで、TOF-PETシステムにおける信頼性の高いポジトロニウム寿命再構成が可能になるか？
• RQ2 取り込み量と寿命イメージの共同推定は、分離した再構成手法と比較して、PLIの精度と定量性をどのように向上させるか？
• RQ3 提案された最大尤度アルゴリズムは、現実的な検出器分解能を有する模擬TOF-PETデータから、真の寿命分布をどの程度正確に回復できるか？
• RQ4 時間分解能の影響を考慮した場合、PCモードとPTモードの両方でアルゴリズムの性能が一貫して維持されるか？
• RQ5 理想化された仮定（完全な時間分解能）に依存せずに、寿命成分を正確に分離できるか？

主な発見

• 提案されたアルゴリズムは、有限な時間分解能（200–600 ps FWHM）下でも、模擬TOF-PETデータから高定量的な精度で寿命イメージを再構成できた。
• 共同最大尤度推定フレームワークにより、同時に取り込み量と寿命分布の両方を回復でき、再構成の忠実度が向上した。
• アルゴリズムは異なる取得モードに対して頑健であり、PCおよびPTの両ケースで正のスケーリング要因を除いて同等のML推定値が得られた。
• 統計モデルに有限な時間分解能を組み込むことで、情報損失を防ぎ、正確な寿命成分の分離が可能になった。
• 正の制約を課したL-BFGS-B最適化により、安定した収束と物理的に妥当な画像再構成が実現された。
• 寿命パラメータに関する導出した勾配により、効率的かつ正確な最適化が可能となり、実用的イメージングシナリオにおける手法の実現可能性が示された。