[論文レビュー] Proton Tracking Algorithm in a Pixel-Based Range Telescope for Proton Computed Tomography

Helge Egil Seime Pettersen, Ilker Meriç+4 more|arXiv (Cornell University)|Jun 17, 2020

Radiation Therapy and Dosimetry被引用数 6

ひとこと要約

本論文は、陽子コンputedトモグラフィー（CT）向けに設計されたピクセル型レンジテレスコープ（デジタルトラッキングカルシメーター、DTC）において、陽子軌道を再構成するためのトラックフォローポンティングアルゴリズムを提示する。このアルゴリズムは、前層のピクセルからトラックを成長させることで、回転偏差を最小化し、数100本の同時陽子トラックを高い効率で再構成に成功しており、1秒あたり500万～2500万個の陽子ビーム強度を実現する。

ABSTRACT

Proton Computed Tomography (CT) is a prototype imaging modality for the reconstruction of the Relative Stopping Power of a patient, for more accurate calculations of the dose distributions in proton therapy dose planning. The measurements needed for the reconstruction of a proton CT image are: i) each initial proton vector incident on the imaged object, ii) each proton vector incident on the front face of the detector and iii) the stopping depth of each proton in the detector. In this study, a track reconstruction algorithm is adapted for a planned pixel-based particle-tracking range telescope for proton CT, called the Digital Tracking Calorimeter (DTC). The algorithm is based on the track-following scheme, in which a growing track searches for deeper-laying activated pixels, while minimizing the accumulated angular change. The algorithm is applied to Monte Carlo-simulated output of the DTC, showing that the DTC is able to reconstruct the tracks and find the depths of up to several hundred simultaneous proton tracks.

研究の動機と目的

• 臨床用陽子CTを可能にするために、ピクセル型レンジテレスコープ向けに、陽子軌道を正確に再構成できる耐障害性の高い陽子トラッキングアルゴリズムの開発。
• 複数のクーロン散乱および非弾性核反応の影響下でも、複数の同時陽子トラックを再構成する課題に対処すること。
• 偽トラック再構成を最小限に抑え、トラッキング効率を向上させることで、陽子CTシステムにおけるビーム強度を最大化すること。
• 特にレンジおよび角度再構成に特化した陽子CTのニーズに適合する、実証済みの高エネルギー物理学のトラッキング技術の適応。
• 散乱や相互作用の影響を含む現実的なビーム条件下でのアルゴリズム性能の評価。

提案手法

• アルゴリズムは、トラックを前部から成長させるトラックフォローポンティング方式を採用し、累積的な回転偏差が最小となるような、より奥のレイヤーの励起ピクセルを探索する。
• 誤再構成を低減するため、現在のトラック方向からの回転偏差が最小となるピクセルクラスタを優先するグリーディー探索戦略を採用する。
• 再構成は、最初の活性ピクセルレイヤーから開始し、層ごとに逐次実行され、有効なピクセル接続がこれ以上見つからなくなるまで継続する。
• アルゴリズムはC++/ROOT6で実装され、単一CPU環境での最適化が行われており、並列処理およびさらなる最適化が可能である。
• 偽トラックの排除は、後処理フィルタを用いて実現され、特に水等価パス長（WEPL）値の3σカットを用いて偽陽性を低減する。
• 実際のビーム条件（散乱や核反応を含む）下での性能評価には、DTCプロトタイプのモンテカルロシミュレーションが用いられる。

実験結果

リサーチクエスチョン

• RQ1 ピクセル型レンジテレスコープにおいて、現実的な散乱および相互作用条件下でも、トラックフォローポンティングアルゴリズムが複数の同時陽子トラックを効果的に再構成できるか？
• RQ2 トラックの混同や偽陽性が顕著に増加する前に、アルゴリズムが処理できる陽子ビーム強度の上限は何か？
• RQ3 非弾性核散乱および複数回のクーロン散乱は、トラック再構成の正確性と偽陽性率にどのように影響を与えるか？
• RQ4 臨床的陽子CTの撮影時間に適したビーム強度において、アルゴリズムは高い効率と低い偽陽性率を維持できるか？
• RQ5 再構成のシード選択（前層 vs. 遠位層）が、トラック再構成性能および計算負荷にどのように影響を与えるか？

主な発見

• アルゴリズムは、DTCプロトタイプからのシミュレートデータにおいて、数100本の同時陽子トラックを成功裏に再構成した。
• 低陽子密度（np < 5）では、偽トラックの主な原因は非弾性核散乱であり、偽再構成の約1%～2%を占める。
• 高密度（np > 50）では、複数回のクーロン散乱が混乱の主な要因となり、再構成品質が著しく低下する。
• アルゴリズムは、ビームスポットサイズや許容可能な偽陽性率に応じて、最大1秒あたり500万～2500万個のビーム強度を達成した。
• 現在の再構成時間は1陽子あたり26 msであり、臨床用途には不十分であり、さらなる最適化および並列処理が必要である。
• WEPL値の3σフィルタは、偽トラックを効果的に低減することができ、高レート条件下での再構成忠実度向上の有効な道筋を示している。