[论文解读] 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）模型，对定点迭代方案进行数学分析。
• 通过迭代方式使用截断DDC矩阵更新射野强度值，旨在校正截断引入的近似误差。
• 该方法被表述为定点更新：x^{k+1} = x^k + α(D^T D x^k - D^T d)，其中D为截断DDC矩阵。
• 在头颈部患者病例上，结合简化和原始FMO模型开展数值实验，以验证该方法。
• 通过监测迭代间射野强度向量的变化来评估收敛性。
• 利用剂量体积直方图（DVH）和伽马分析指标，比较定点解、朴素截断解与完整DDC解的计划质量。

实验结果

研究问题

• RQ1 在IMRT优化中，当应用于截断DDC矩阵时，定点迭代算法是否收敛？
• RQ2 当定点迭代收敛时，所得解是否等同于使用完整DDC矩阵所获得的解？
• RQ3 若不等价，定点解是否比直接从截断DDC矩阵导出的标准解更准确？
• RQ4 DDC矩阵截断对计划质量有何影响？定点迭代能否缓解此影响？

主要发现

• 定点迭代方法在合理截断DDC矩阵条件下可实现收敛，该结论得到数学分析与数值实验的共同验证。
• 收敛后的定点解与使用完整DDC矩阵所得解不完全相同，表明截断导致了不同的固定点。
• 尽管未完全匹配完整DDC解，定点解仍显著优于直接从截断矩阵导出的朴素解。
• 在头颈部患者病例的数值结果表明，定点方法提升了计划质量，表现为更优的剂量分布及靶区冷区减少。
• 该方法有效降低了DDC矩阵截断引入的误差，从而在不需完整矩阵存储的前提下提升了计划准确性。