[论文解读] Feasibility-Seeking and Superiorization Algorithms Applied to Inverse Treatment Planning in Radiation Therapy

R. Davidi, Y. Censor+3 more|arXiv (Cornell University)|Feb 6, 2014

Spectroscopy Techniques in Biomedical and Chemical Research被引用 4

一句话总结

本文提出使用优越化方法（SM）改进调强放射治疗（IMRT）中的逆向治疗计划，以计算成本较高的最小化方法替代具有抗扰动能力的可行性搜索算法，从而降低射野强度的总变差（TV）。SM在比标准可行性搜索方法更少的迭代次数内实现了临床可接受的剂量分布，且在前列腺癌病例研究中，符合RTOG 0815标准，提高了适形性并减少了热点区域。

ABSTRACT

We apply the recently proposed superiorization methodology (SM) to the inverse planning problem in radiation therapy. The inverse planning problem is represented here as a constrained minimization problem of the total variation (TV) of the intensity vector over a large system of linear two-sided inequalities. The SM can be viewed conceptually as lying between feasibility-seeking for the constraints and full-fledged constrained minimization of the objective function subject to these constraints. It is based on the discovery that many feasibility-seeking algorithms (of the projection methods variety) are perturbation-resilient, and can be proactively steered toward a feasible solution of the constraints with a reduced, thus superiorized, but not necessarily minimal, objective function value.

研究动机与目标

• 解决当前IMRT逆向计划方法中使用完整约束最小化所导致的高计算成本问题。
• 研究优越化是否能在计算负担减少的情况下生成临床可接受的治疗计划。
• 评估SM是否能在保持约束可行性的同时，提高解的质量（即降低目标函数值）。
• 比较SM与标准可行性搜索算法在收敛速度和剂量分布质量方面的表现。
• 证明使用TV-优越化生成可交付的、分段常数的强度图在IMRT中的可行性。

提出的方法

• 将逆向治疗计划问题表述为在大规模双侧线性不等式系统上对总变差（TV）进行约束最小化。
• 将优越化方法（SM）应用于不等式系统的代数重建技术（ART），该方法已知具有有界扰动鲁棒性。
• SM在可行性迭代过程中引入受控扰动，以引导解向更低的TV值方向发展，而无需进行完整最小化。
• 该方法使用基于投影的算法，交替执行可行性搜索步骤和通过TV最小化降低目标函数值的步骤。
• 算法以零或均匀的射野强度初始化，并在强制执行剂量约束的同时迭代更新解。
• 最终解通过剂量体积直方图（DVH）分析以及符合RTOG 0815临床接受标准进行评估。

实验结果

研究问题

• RQ1 优越化是否能比标准可行性搜索算法更快生成临床可接受的IMRT计划？
• RQ2 优越化是否能在保持剂量约束的前提下降低射野强度的总变差？
• RQ3 SM是否能在不增加计算成本的情况下改善计划质量（如减少热点、提高适形性）？
• RQ4 初始射野强度的选择如何影响优越化算法的收敛速度和计划质量？
• RQ5 当标准可行性搜索方法在相同迭代限制内失败时，SM是否足够稳健以达到临床目标？

主要发现

• 在实验1（零初始化）中，TV-优越化算法在12次迭代内达到RTOG 0815标准，而可行性搜索算法未能达到最低PTV剂量（56.13 Gy vs. 所需的75.24 Gy）。
• 在实验2（均匀初始化）中，TV-优越化算法在7次迭代内达到可接受计划，PTV最低剂量为77.80 Gy，而未优越化的版本为76.15 Gy。
• 在实验1中，受照剂量超过60 Gy的直肠体积减少至34.50%（TV-优越化）对比8.50%（未优越化），表明对危及器官（OAR）的保护更好。
• 在实验2中，优越化方案的直肠体积超过60 Gy为36.90%，而未优越化方案为40.50%，表明OAR保护有所改善。
• 在实验1中，直肠最大剂量分别为82.64 Gy（TV-优越化）和82.71 Gy（未优越化），表明OAR保护水平相当。
• DVH曲线显示，优越化解在更少的迭代次数内实现了更好的适形性和均匀性，尤其在PTV覆盖和OAR保护方面表现更优。