Robust optimization for intensity-modulated proton therapy with soft spot sensitivity regularization.
Proton dose distribution is sensitive to uncertainties in range estimation and patient positioning. Currently, the proton robustness is managed by worst-case scenario optimization methods, which are computationally inefficient. To overcome these challenges, we develop a novel intensity-modulated proton therapy (IMPT) optimization method that integrates dose fidelity with a sensitivity term that describes dose perturbation as the result of range and positioning uncertainties.
In the integrated optimization framework, the optimization cost function is formulated to include two terms: a dose fidelity term and a robustness term penalizing the inner product of the scanning spot sensitivity and intensity. The sensitivity of an IMPT scanning spot to perturbations is defined as the dose distribution variation induced by range and positioning errors. To evaluate the sensitivity, the spatial gradient of the dose distribution of a specific spot is first calculated. The spot sensitivity is then determined by the total absolute value of the directional gradients of all affected voxels. The fast iterative shrinkage-thresholding algorithm is used to solve the optimization problem. This method was tested on three skull base tumor (SBT) patients and three bilateral head-and-neck (H&N) patients. The proposed sensitivity-regularized method (SenR) was implemented on both clinic target volume (CTV) and planning target volume (PTV). They were compared with conventional PTV-based optimization method (Conv) and CTV-based voxel-wise worst-case scenario optimization approach (WC).
Under the nominal condition without uncertainties, the three methods achieved similar CTV dose coverage, while the CTV-based SenR approach better spared organs at risks (OARs) compared with the WC approach, with an average reduction of [Dmean, Dmax] of [4.72, 3.38] GyRBE for the SBT cases and [2.54, 3.33] GyRBE for the H&N cases. The OAR sparing of the PTV-based SenR method was comparable with the WC method. The WC method, and SenR approaches all improved the plan robustness from the conventional PTV-based method. On average, under range uncertainties, the lowest [D95%, V95%, V100%] of CTV were increased from [93.75%, 88.47%, 47.37%] in the Conv method, to [99.28%, 99.51%, 86.64%] in the WC method, [97.71%, 97.85%, 81.65%] in the SenR-CTV method and [98.77%, 99.30%, 85.12%] in the SenR-PTV method, respectively. Under setup uncertainties, the average lowest [D95%, V95%, V100%] of CTV were increased from [95.35%, 94.92%, 65.12%] in the Conv method, to [99.43%, 99.63%, 87.12%] in the WC method, [96.97%, 97.13%, 77.86%] in the SenR-CTV method, and [98.21%, 98.34%, 83.88%] in the SenR-PTV method, respectively. The runtime of the SenR optimization is eight times shorter than that of the voxel-wise worst-case method.
We developed a novel computationally efficient robust optimization method for IMPT. The robustness is calculated as the spot sensitivity to both range and shift perturbations. The dose fidelity term is then regularized by the sensitivity term for the flexibility and trade-off between the dosimetry and the robustness. In the stress test, SenR is more resilient to unexpected uncertainties. These advantages in combination with its fast computation time make it a viable candidate for clinical IMPT planning.
Gu W
，Ruan D
，O'Connor D
，Zou W
，Dong L
，Tsai MY
，Jia X
，Sheng K
... -
《-》
Exploratory Study of 4D versus 3D Robust Optimization in Intensity Modulated Proton Therapy for Lung Cancer.
The purpose of this study was to compare the impact of uncertainties and interplay on 3-dimensional (3D) and 4D robustly optimized intensity modulated proton therapy (IMPT) plans for lung cancer in an exploratory methodology study.
IMPT plans were created for 11 nonrandomly selected non-small cell lung cancer (NSCLC) cases: 3D robustly optimized plans on average CTs with internal gross tumor volume density overridden to irradiate internal target volume, and 4D robustly optimized plans on 4D computed tomography (CT) to irradiate clinical target volume (CTV). Regular fractionation (66 Gy [relative biological effectiveness; RBE] in 33 fractions) was considered. In 4D optimization, the CTV of individual phases received nonuniform doses to achieve a uniform cumulative dose. The root-mean-square dose-volume histograms (RVH) measured the sensitivity of the dose to uncertainties, and the areas under the RVH curve (AUCs) were used to evaluate plan robustness. Dose evaluation software modeled time-dependent spot delivery to incorporate interplay effect with randomized starting phases of each field per fraction. Dose-volume histogram (DVH) indices comparing CTV coverage, homogeneity, and normal tissue sparing were evaluated using Wilcoxon signed rank test.
4D robust optimization plans led to smaller AUC for CTV (14.26 vs 18.61, respectively; P=.001), better CTV coverage (Gy [RBE]) (D95% CTV: 60.6 vs 55.2, respectively; P=.001), and better CTV homogeneity (D5%-D95% CTV: 10.3 vs 17.7, respectively; P=.002) in the face of uncertainties. With interplay effect considered, 4D robust optimization produced plans with better target coverage (D95% CTV: 64.5 vs 63.8, respectively; P=.0068), comparable target homogeneity, and comparable normal tissue protection. The benefits from 4D robust optimization were most obvious for the 2 typical stage III lung cancer patients.
Our exploratory methodology study showed that, compared to 3D robust optimization, 4D robust optimization produced significantly more robust and interplay-effect-resistant plans for targets with comparable dose distributions for normal tissues. A further study with a larger and more realistic patient population is warranted to generalize the conclusions.
Liu W
，Schild SE
，Chang JY
，Liao Z
，Chang YH
，Wen Z
，Shen J
，Stoker JB
，Ding X
，Hu Y
，Sahoo N
，Herman MG
，Vargas C
，Keole S
，Wong W
，Bues M
... -
《-》
ResearchGate
(全网免费下载)
钛学术
(全网免费下载)
