[論文レビュー] Fast non-coplanar beam orientation optimization based on group sparsity

Daniel O’Connor, Yevgen Voronenko+3 more|arXiv (Cornell University)|Oct 15, 2017

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

ひとこと要約

本稿では、グループスパarsityとFISTAアルゴリズムを用いた高速で非同一面性のビーム方位最適化手法を提案し、500〜800件の候補ビームから最適なビーム角度を効率的に選択可能にした。この手法は5分未満で優れた線量学的品質を達成し、従来のグループスパarsity手法が非同一面性ケースにおいてあまりに遅すぎて実用的でなかったのに対し、顕著に優れている。

ABSTRACT

The selection of beam orientations, which is a key step in radiation treatment planning, is particularly challenging for non-coplanar radiotherapy systems due to the large number of candidate beams. In this paper, we report progress on the group sparsity approach to beam orientation optimization, wherein beam angles are selected by solving a large scale fluence map optimization problem with an additional group sparsity penalty term that encourages most candidate beams to be inactive. The optimization problem is solved using an accelerated proximal gradient method, the Fast Iterative Shrinkage-Thresholding Algorithm (FISTA). We derive a closed-form expression for a relevant proximal operator which enables the application of FISTA. The proposed algorithm is used to create non-coplanar treatment plans for four cases (including head and neck, lung, and prostate cases), and the resulting plans are compared with clinical plans. The dosimetric quality of the group sparsity treatment plans is superior to that of the clinical plans. Moreover, the runtime for the group sparsity approach is typically about 5 minutes. Problems of this size could not be handled using the previous group sparsity method for beam orientation optimization, which was slow to solve much smaller coplanar cases. This work demonstrates for the first time that the group sparsity approach, when combined with an accelerated proximal gradient method such as FISTA, works effectively for non-coplanar cases with 500-800 candidate beams.

研究の動機と目的

• 非同一面性Intensity-Modulated Radiation Therapy (IMRT)におけるビーム方位選択の課題に対処すること。手動による選択は、膨大な数の候補ビームがあるため現実的ではない。
• ビーム数の最小化のための代替として凸なグループスパarsityペナルティを用い、組み合わせ的ビーム角度選択の計算的非実行可能性を克服すること。
• FISTAを用いた加速型近接勾配法の再定式化により、スケーラブルで包括的なビーム選択を可能とし、高速収束を実現すること。
• 非同一面性IMRTにおけるグループスパarsityに基づくビーム最適化の実現可能性と優位性を示すこと。これは、従来の実行時間制約により実現不可能であった。
• 従来のグループスパarsityアプローチが小規模な同一面性ケースですら数時間も要していたのに対し、臨床的に意味のある計画品質を大幅に短時間で達成すること。

提案手法

• ビーム方位最適化を、ビーム活性化におけるスパarsityを促進する$\ell_{2,1}$-ノルムグループスパarsityペナルティを備えた大規模なフルエンスマップ最適化問題として定式化する。
• 問題を解くためにFast Iterative Shrinkage-Thresholding Algorithm (FISTA)を適用し、$O(1/k^2)$の収束速度を活用して高速化を実現する。
• グループスパarsity項の閉形式近接作用素を導出することで、FISTAフレームワーク内での計算を効率化する。
• FISTAにおけるラインサーチ戦略を用い、ステップサイズを動的に調整し、グローバル収束保証を確保する。
• PTVおよびOARの線量制約を、重み付き最小二乗項を用いて目的関数に統合し、被照射領域の被覆と臓器遮断のバランスを取る。
• 各ビームのフルエンスをグループとして扱うブロック構造アプローチを実装し、近接作用素による効率的なグループワイドスレッショーティングを可能にする。

実験結果

リサーチクエスチョン

• RQ1 500〜800件の候補ビームを有する非同一面性IMRTにおいて、グループスパarsityに基づくビーム方位最適化が計算的に実行可能となるか？
• RQ2 閉形式近接作用素を有するFISTAを用いることで、従来の近接勾配法と比較して、ビーム方位最適化における実行時間が著しく短縮されるか？
• RQ3 グループスパarsityに基づくビーム選択は、非同一面性IMRTにおいて臨床的に選択されたビームと比較して線量学的に優れた計画を生成できるか？
• RQ4 FISTAの加速収束は、従来では非実行的とされた大規模な非同一面性ビーム選択問題を処理するのに十分か？
• RQ5 特に非同一面性ケースにおいて、本手法は従来のグループスパarsityアプローチと比較して、実行時間および計画品質で優れているか？

主な発見

• 提案されたFISTAベースのグループスパarsity手法は、頭頸部、肺、前立腺腫瘍を含む4例の症例において、臨床的計画を上回る治療計画品質を達成した。
• 500〜800件の候補ビームを有する問題におけるビーム方位最適化の平均実行時間は約5分であり、従来のグループスパarsity手法が小規模な同一面性ケースですら数時間も要していたのに対し、顕著な改善であった。
• 従来のグループスパarsity実装では計算ボトルネックにより非現実的だった非同一面性ビーム方位最適化問題を、本手法は正常に解消できた。
• FISTAと閉形式近接作用素を組み合わせたアプローチにより、非同一面性IMRTにおける包括的ビーム選択の実現可能性が示された。
• アルゴリズムは$O(1/k^2)$の収束を達成し、従来の近接勾配法の$O(1/k)$の収束速度を著しく上回った。これは大規模問題を処理する上で極めて重要であった。
• グループスパarsity項の近接作用素が正当化され、計算的に効率的であることが確認され、臨床的計画ワークフローにおけるリアルタイム応用が可能となった。