[論文レビュー] Improved Gamma-Ray Point Source Quantification in Three Dimensions by Modeling Attenuation in the Scene

Mark S. Bandstra, Daniel Hellfeld+9 more|arXiv (Cornell University)|Nov 1, 2021

Medical Imaging Techniques and Applications参考文献 40被引用数 16

ひとこと要約

本稿では、LiDARで得られたボクセル化されたシーン幾何学を用いて中間物質による減衰をモデル化することで、ガンマ線源の3次元位置特定と定量化を行う手法を提示する。最大尤度推定とレイトレーシングを組み合わせて経路長に依存する減衰を計算することで、源の位置、活動量、減衰係数を同時に推定し、複雑な遮蔽を伴う現実世界の測定において、源パラメータを正確に再構成することを示した。

ABSTRACT

Using a series of detector measurements taken at different locations to localize a source of radiation is a well-studied problem. The source of radiation is sometimes constrained to a single point-like source, in which case the location of the point source can be found using techniques such as maximum likelihood. Recent advancements have shown the ability to locate point sources in 2-D and even 3-D but few have studied the effect of intervening material on the problem. In this work, we examine gamma-ray data taken from a freely moving system and develop voxelized 3-D models of the scene using data from its onboard light detection and ranging (LiDAR) unit. Ray casting is used to compute the distance each gamma ray travels through the scene material, which is then used to calculate attenuation assuming a single attenuation coefficient for solids within the geometry. Parameter estimation using maximum likelihood is performed to simultaneously find the attenuation coefficient, source activity, and source position that best match the data. Using a simulation, we validate the ability of this method to reconstruct the true location and activity of a source, along with the true attenuation coefficient of the structure it is inside, and then we apply the method to measured data with sources and find good agreement.

研究の動機と目的

• 中間物質による未モデル化の減衰が原因で生じる3次元環境下での不正確な源の定量と位置特定の課題に対処すること。
• 移動可能な検出器のデータを用いて、3次元シーン内で源の位置、活動量、減衰係数を同時に推定する手法を開発すること。
• シミュレーションおよび測定データの両方を用いて、本手法が真の源パラメータおよび減衰特性をどれだけ正確に回復できるかを検証すること。
• 移動型放射線検出システムにおけるポイント源尤度（PSL）フレームワークにシーンの減衰を組み込む可能性を示すこと。

提案手法

• 本手法は、搭載されたLiDAR点群データから生成されたボクセル化された3次元シーンモデルを用い、環境内の固体物質を表現する。
• 各検出器位置から候補となる源位置へとレイトレーシングを行い、シーン内における固体物質を通過する全経路長を計算する。
• 減衰は、すべての固体に対して一様な減衰係数（平均自由行程、λsolid）を用いてモデル化し、レイトレースされた経路に沿ってガンマ線の期待計数率を計算する。
• 最大尤度推定を用いて、測定された検出器計数に最もよく一致するよう、源の位置、源の活動量、λsolidを同時に最適化する。
• 最適化は候補となる源位置の離散グリッド上で実行され、不確実性は2σ信頼区間で評価される。
• 本手法は、シミュレーションデータを用いて検証され、鉄製輸送容器内に配置された137Csおよび133Ba源の実測データに応用された。

実験結果

リサーチクエスチョン

• RQ1 源の位置、活動量、減衰係数を同時に推定することで、複雑なシーンにおける3次元ガンマ線源位置特定の精度が向上するか？
• RQ2 中間物質による減衰が存在する状況下で、本手法が真の源パラメータをどの程度正確に回復できるか？
• RQ3 LiDARで得られた3次元シーンモデルは、直接の物質特性測定ができない状況下でも、どれほど正確な減衰モデル化を可能にするか？
• RQ4 実際の減衰物質が表面幾何形状のみで表現されない場合、本手法にどのような制限が生じるか？

主な発見

• 本手法は、実測データにおいて137Csおよび133Ba源の真の源位置と活動量を良好に再構成した。結果は真値とよく一致した。
• 137Cs源の場合、最適適合の減衰平均自由行程（λsolid）は2.58 mであった。これはボクセルサイズと有効壁厚さを考慮した場合の期待値と整合的であり、λsolid/∆ の比が約5.9（95%信頼区間4.5–9.8）の範囲内に収まった。
• 133Ba源の場合、最適適合のλsolid/∆ 比は9.2–58.7であったが、これは表記値を含まないが、複数のボクセルにまたがるより厚い減衰体積および内部の追加物質によって説明された。
• 本手法は、表面（例：コンテナの壁）による減衰が支配的である場合には良好に機能したが、内部のバルク物質が顕著に減衰に寄与し、かつLiDARで解像できない場合には性能が低下した。
• 処理時間は長く（単一CPUで100–210分）、しかし著者らはGPU加速および最適化の改善により、近似的にリアルタイム動作が可能になると指摘している。
• 未観測の内部体積をモデル化するためには、追加の文脈的情報（例：物質分類、セグメンテーション）の導入が必要であり、そのことが耐性向上に不可欠であると本研究は示唆している。