The dose-volume histogram (DVH) is a clinically relevant criterion to evaluate the quality of a treatment plan. It is hence desirable to incorporate DVH constraints into treatment plan optimization for intensity modulated radiation therapy. Yet, the direct inclusion of the DVH constraints into a treatment plan optimization model typically leads to great computational difficulties due to the non-convex nature of these constraints. To overcome this critical limitation, we propose a new convex-moment-based optimization approach. Our main idea is to replace the non-convex DVH constraints by a set of convex moment constraints. In turn, the proposed approach is able to generate a Pareto-optimal plan whose DVHs are close to, or if possible even outperform, the desired DVHs. In particular, our experiment on a prostate cancer patient case demonstrates the effectiveness of this approach by employing two and three moment formulations to approximate the desired DVHs.

Stereotactic radiosurgery (SRS) is an effective technique to treat brain metastasis for which several inverse planning methods may be appropriate. We compare three different optimization models for segment duration optimization in SRS using Leksell Gamma Knife Icon (Elekta, Stockholm, Sweden). We investigate (1) a linear programming approach, (2) a piecewise quadratic penalty approach, and (3) an unconstrained convex moment-based penalty approach. We examine the performances of these approaches using anonymized data from 14 previously treated cases. In addition, we investigate the important modeling question of selecting weights for the objective functions where we use a simulated annealing algorithm to determine these weights for each model. The inverse plans obtained via optimization models are compared against each other and against the clinical plans. The three inverse planning models can all yield optimal treatment plans in a reasonable amount of time and the treatment plans obtained by these models meet or exceed clinical guidelines while displaying high conformity.

Purpose: To compare methods to incorporate the Dose Volume Histogram (DVH) curves into the treatment planning optimization. Method: The performance of three methods, namely, the conventional Mixed Integer Programming (MIP) model, a convex moment‐based constrained optimization approach, and an unconstrained convex moment‐based penalty approach, is compared using anonymized data of a prostate cancer patient. Three plans we generated using the corresponding optimization models. Four Organs at Risk (OARs) and one Tumor were involved in the treatment planning. The OARs and Tumor were discretized into total of 50,221 voxels. The number of beamlets was 943. We used commercially available optimization software Gurobi and Matlab to solve the models. Plan comparison was done by recording the model runtime followed by visual inspection of the resulting dose volume histograms. Conclusion: We demonstrate the effectiveness of the moment‐based approaches to replicate the set of prescribed DVH curves. The unconstrained convex moment‐based penalty approach is concluded to have the greatest potential to reduce the computational effort and holds a promise of substantial computational speed up.

No abstract