[論文レビュー] The fixed-point iteration method for IMRT optimization with truncated dose deposition coefficient matrix

Zhen Tian, Masoud Zarepisheh+2 more|arXiv (Cornell University)|Mar 14, 2013

Advanced Radiotherapy Techniques参考文献 11被引用数 7

ひとこと要約

本稿では、メモリおよび計算コストを削減するために一般的に用いられる、切断された線量蓄積係数（DDC）行列を用いた強度調制放射線治療（IMRT）治療計画の品質を向上させるための固定点反復法を提案する。この手法は、完全なDDC行列の解とは正確に一致しないものの、切断問題に対するナードル解よりも著しく高い精度に収束する。したがって、完全な行列の保存を必要とせずに、治療計画の精度を向上させることができる。

ABSTRACT

In the treatment plan optimization for intensity modulated radiation therapy (IMRT), dose-deposition coefficient (DDC) matrix is often pre-computed to parameterize the dose contribution to each voxel in the volume of interest from each beamlet of unit intensity. However, due to the limitation of computer memory and the requirement on computational efficiency, in practice matrix elements of small values are usually truncated, which inevitably compromises the quality of the resulting plan. A fixed-point iteration scheme has been applied in IMRT optimization to solve this problem, which has been reported to be effective and efficient based on the observations of the numerical experiments. In this paper, we aim to point out the mathematics behind this scheme and to answer the following three questions: 1) whether the fixed-point iteration algorithm converges or not? 2) when it converges, whether the fixed point solution is same as the original solution obtained with the complete DDC matrix? 3) if not the same, whether the fixed point solution is more accurate than the naive solution of the truncated problem obtained without the fixed-point iteration? To answer these questions, we first performed mathematical analysis and deductions using a simplified fluence map optimization (FMO) model. Then we conducted numerical experiments on a head-and-neck patient case using both the simplified and the original FMO model. Both our mathematical analysis and numerical experiments demonstrate that with proper DDC matrix truncation, the fixed-point iteration can converge. Even though the converged solution is not the one that we obtain with the complete DDC matrix, the fixed-point iteration scheme could significantly improve the plan accuracy compared with the solution to the truncated problem obtained without the fixed-point iteration.

研究の動機と目的

• メモリおよび計算制約による線量蓄積係数（DDC）行列の小さな要素を切り捨てることで生じるIMRT治療計画品質の低下を是正すること。
• IMRT最適化における切断DDC行列の文脈で、固定点反復法の数学的挙動を分析すること。
• 固定点反復解が、切断DDC行列から直接得られる標準解よりも優れているかどうかを評価すること。
• 固定点解が収束するかどうか、および完全なDDC行列からの最適解に近似するかどうかを特定すること。

提案手法

• 固定点反復スキームの数学的分析を行うために、簡略化されたフルエンスマップ最適化（FMO）モデルを用いる。
• 固定点反復法を繰り返し適用して、切り捨てによって生じる近似誤差を補正するためのフルエンス値を更新する。
• この手法は固定点更新として定式化される：x^{k+1} = x^k + α(D^T D x^k - D^T d)，ここでDは切断DDC行列である。
• 簡略化されたFMOモデルおよび元のFMOモデルを用いて、頭頸部患者症例における数値実験を実施し、手法の妥当性を検証する。
• 反復間のフルエンスベクトルの変化をモニタリングすることで収束を評価する。
• ドーズボリュームヒストグラム（DVH）およびガンマ解析指標を用いて、固定点解、ナードルな切断解、完全なDDC解の間で計画品質を比較する。

実験結果

リサーチクエスチョン

• RQ1 IMRT最適化において、切断DDC行列に固定点反復アルゴリズムを適用した場合、収束するか？
• RQ2 固定点反復が収束した場合、その解は完全なDDC行列から得られる解と同等か？
• RQ3 同等でない場合、固定点解は、切断DDC行列から直接得られる標準解よりも正確か？
• RQ4 DDC行列の切り捨てが計画品質に与える影響は何か？ また、固定点反復はその影響を緩和できるか？

主な発見

• 数学的分析および数値実験の両方により、適切なDDC行列の切り捨てがなされた場合、固定点反復法は収束することが確認された。
• 収束した固定点解は、完全なDDC行列から得られる解とは同一ではないため、切り捨てによって異なる固定点に到達していることが示された。
• 完全なDDC行列の解とは一致しないものの、固定点解は、切断行列から直接得られるナードル解よりも著しく高い精度であることが判明した。
• 頭頸部患者症例における数値結果から、固定点法が計画品質を向上させることを示しており、良好な線量分布および標的領域内の冷結節の低減が確認された。
• この手法は、DDC行列の切り捨てによって生じる誤差を効果的に低減し、完全な行列の保存を要せずとも計画精度を向上させることができる。