Target volume uncertainty and a method to visualize its effect on the target dose prescription

Dosimetric analysis of cervical cancer stage IIB patients treated with volumetric modulated arc therapy using plan uncertainty parameters module of Varian Eclipse treatment planning system

Introduction. The present study aims to investigate the dosimetric and radiobiological impact of patient setup errors (PSE) on the target and organs at risk (OAR) of the cervix carcinoma stage IIB patients treated with volumetric-modulated arc therapy (VMAT) delivery technique using plan uncertainty parameters module of Varian Eclipse treatment planning system and in-house developed DVH Analyzer program.Materials and Methods. A total of 976 VMAT plans were generated to simulate the PSE in the base plan that varies from -10 mm to 10 mm in a step size of 1 mm in x- (lateral), y- (craniocaudal), and z- (anteroposterior) directions. The different OAR and tumor (PTV) volumes were delineated in each case. Various plan quality metrics, such as conformity index (CI) and homogeneity index (HI), as well as radiobiological quantities, such as tumor control probability (TCP) and normal tissue control probability (NTCP), were calculated from the DVH bands generated from the cohort of treatment plans associated with each patient case, using an in-house developed 'DVH Analyzer' program. The extracted parameters were statistically analyzed and compared with the base plan's dosimetric parameters having no PSE.Results. The maximum variation of (i) 2.4%, 21.5%, 0.8%, 2.5% in D_2ccof bladder, rectum, small bowel and sigmoid colon respectively; (ii) 19.3% and 18.9% in D_maxof the left and right femoral heads (iii) 16.9% in D_95%of PTV (iv) 12.1% in NTCP of sigmoid colon were observed with change of PSE in all directions. TCP was found to be considerably affected for PSEs larger than 4 mm in x⁺, y⁺, z⁺directions and 7 mm in x^-, y^-and z^-directions, respectively.Conclusion. This study presents the effect of PSE on TCP and NTCP for the cervix carcinoma cases treated with VMAT technique and also recommends daily image guidance to mitigate the effects of PSE.

Bayesian characterization of uncertainty in intra-subject non-rigid registration

In settings where high-level inferences are made based on registered image data, the registration uncertainty can contain important information. In this article, we propose a Bayesian non-rigid registration framework where conventional dissimilarity and regularization energies can be included in the likelihood and the prior distribution on deformations respectively through the use of Boltzmann's distribution. The posterior distribution is characterized using Markov Chain Monte Carlo (MCMC) methods with the effect of the Boltzmann temperature hyper-parameters marginalized under broad uninformative hyper-prior distributions. The MCMC chain permits estimation of the most likely deformation as well as the associated uncertainty. On synthetic examples, we demonstrate the ability of the method to identify the maximum a posteriori estimate and the associated posterior uncertainty, and demonstrate that the posterior distribution can be non-Gaussian. Additionally, results from registering clinical data acquired during neurosurgery for resection of brain tumor are provided; we compare the method to single transformation results from a deterministic optimizer and introduce methods that summarize the high-dimensional uncertainty. At the site of resection, the registration uncertainty increases and the marginal distribution on deformations is shown to be multi-modal.

Intensity modulated radiation therapy with field rotation--a time-varying fractionation study

This paper proposes a novel mathematical approach to the beam selection problem in intensity modulated radiation therapy (IMRT) planning. The approach allows more beams to be used over the course of therapy while limiting the number of beams required in any one session. In the proposed field rotation method, several sets of beams are interchanged throughout the treatment to allow a wider selection of beam angles than would be possible with fixed beam orientations. The choice of beamlet intensities and the number of identical fractions for each set are determined by a mixed integer linear program that controls jointly for the distribution per fraction and the cumulative dose distribution delivered to targets and critical structures. Trials showed the method allowed substantial increases in the dose objective and/or sparing of normal tissues while maintaining cumulative and fraction size limits. Trials for a head and neck site showed gains of 25%-35% in the objective (average tumor dose) and for a thoracic site gains were 7%-13%, depending on how strict the fraction size limits were set. The objective did not rise for a prostate site significantly, but the tolerance limits on normal tissues could be strengthened with the use of multiple beam sets.

Visualization of a variety of possible dosimetric outcomes in radiation therapy using dose-volume histogram bands

PURPOSE

Dose-volume histograms (DVH) are the most common tool used in the appraisal of the quality of a clinical treatment plan. However, when delivery uncertainties are present, the DVH may not always accurately describe the dose distribution actually delivered to the patient. We present a method, based on DVH formalism, to visualize the variability in the expected dosimetric outcome of a treatment plan.

METHODS

For a case of chordoma of the cervical spine, we compared 2 intensity modulated proton therapy plans. Treatment plan A was optimized based on dosimetric objectives alone (ie, desired target coverage, normal tissue tolerance). Plan B was created employing a published probabilistic optimization method that considered the uncertainties in patient setup and proton range in tissue. Dose distributions and DVH for both plans were calculated for the nominal delivery scenario, as well as for scenarios representing deviations from the nominal setup, and a systematic error in the estimate of range in tissue. The histograms from various scenarios were combined to create DVH bands to illustrate possible deviations from the nominal plan for the expected magnitude of setup and range errors.

RESULTS

In the nominal scenario, the DVH from plan A showed superior dose coverage, higher dose homogeneity within the target, and improved sparing of the adjacent critical structure. However, when the dose distributions and DVH from plans A and B were recalculated for different error scenarios (eg, proton range underestimation by 3 mm), the plan quality, reflected by DVH, deteriorated significantly for plan A, while plan B was only minimally affected. In the DVH-band representation, plan A produced wider bands, reflecting its higher vulnerability to delivery errors, and uncertainty in the dosimetric outcome.

CONCLUSIONS

The results illustrate that comparison of DVH for the nominal scenario alone does not provide any information about the relative sensitivity of dosimetric outcome to delivery uncertainties. Thus, such comparison may be misleading and may result in the selection of an inferior plan for delivery to a patient. A better-informed decision can be made if additional information about possible dosimetric variability is presented; for example, in the form of DVH bands.

Using registration uncertainty visualization in a user study of a simple surgical task

We present a novel method to visualize registration uncertainty and a simple study to motivate the use of uncertainty visualization in computer-assisted surgery. Our visualization method resulted in a statistically significant reduction in the number of attempts required to localize a target, and a statistically significant reduction in the number of targets that our subjects failed to localize. Most notably, our work addresses the existence of uncertainty in guidance and offers a first step towards helping surgeons make informed decisions in the presence of imperfect data.