International audiencePurpose:This work aims at investigating intensity corrected cone-beam x-ray computed tomography (CBCT) images for accurate dose calculation in adaptive intensity modulated proton therapy (IMPT) for prostate and head and neck (H&N) cancer. A deformable image registration (DIR)-based method and a scatter correction approach using the image data obtained from DIR as prior are characterized and compared on the basis of the same clinical patient cohort for the first time.Methods:Planning CT (pCT) and daily CBCT data (reconstructed images and measured projections) of four H&N and four prostate cancer patients have been considered in this study. A previously validated Morphons algorithm was used for DIR of the planning CT to the current CBCT image, yielding a so-called virtual CT (vCT). For the first time, this approach was translated from H&N to prostate cancer cases in the scope of proton therapy. The warped pCT images were also used as prior for scatter correction of the CBCT projections for both tumor sites. Single field uniform dose and IMPT (only for H&N cases) treatment plans have been generated with a research version of a commercial planning system. Dose calculations on vCT and scatter corrected CBCT (CBCT cor) were compared by means of the proton range and a gamma-index analysis. For the H&N cases, an additional diagnostic replanning CT (rpCT) acquired within three days of the CBCT served as additional reference. For the prostate patients, a comprehensive contour comparison of CBCT and vCT, using a trained physician’s delineation, was performed.Results:A high agreement of vCT and CBCT cor was found in terms of the proton range and gamma-index analysis. For all patients and indications between 95% and 100% of the proton dose profiles in beam’s eye view showed a range agreement of better than 3 mm. The pass rate in a (2%,2 mm) gamma-comparison was between 96% and 100%. For H&N patients, an equivalent agreement of vCT and CBCT cor to the reference rpCT was observed. However, for the prostate cases, an insufficient accuracy of the vCT contours retrieved from DIR was found, while the CBCT cor contours showed very high agreement to the contours delineated on the raw CBCT.Conclusions:For H&N patients, no considerable differences of vCT and CBCT cor were found. For prostate cases, despite the high dosimetric agreement, the DIR yields incorrect contours, probably due to the more pronounced anatomical changes in the abdomen and the reduced soft-tissue contrast in the CBCT. Using the vCT as prior, these inaccuracies can be overcome and images suitable for accurate delineation and dose calculation in CBCT-based adaptive IMPT can be retrieved from scatter correction of the CBCT projections

Purpose: To demonstrate the feasibility of proton dose calculation on scatter-corrected cone-beam computed tomographic (CBCT) images for the purpose of adaptive proton therapy. Methods: CBCT projection images were acquired from anthropomorphic phantoms and a prostate patient using an on-board imaging system of an Elekta infinity linear accelerator. Two previously introduced techniques were used to correct the scattered x-rays in the raw projection images: uniform scatter correction (CBCT us ) and a priori CT-based scatter correction (CBCT ap ). CBCT images were reconstructed using a standard FDK algorithm and GPU-based reconstruction toolkit. Soft tissue ROI-based HU shifting was used to improve HU accuracy of the uncorrected CBCT images and CBCT us , while no HU change was applied to the CBCT ap . The degree of equivalence of the corrected CBCT images with respect to the reference CT image (CT ref ) was evaluated by using angular profiles of water equivalent path length (WEPL) and passively scattered proton treatment plans. The CBCT ap was further evaluated in more realistic scenarios such as rectal filling and weight loss to assess the effect of mismatched prior information on the corrected images. Results: The uncorrected CBCT and CBCT us images demonstrated substantial WEPL discrepancies (7.3 ± 5.3 mm and 11.1 ± 6.6 mm, respectively) with respect to the CT ref , while the CBCT ap images showed substantially reduced WEPL errors (2.4 ± 2.0 mm). Similarly, the CBCT ap -based treatment plans demonstrated a high pass rate (96.0% ± 2.5% in 2 mm/2% criteria) in a 3D gamma analysis. Conclusions: A priori CT-based scatter correction technique was shown to be promising for adaptive proton therapy, as it achieved equivalent proton dose distributions and water equivalent path lengths compared to those of a reference CT in a selection of anthropomorphic phantoms. C 2015 American Association of Physicists in Medicine. [http://dx

To develop an online plan adaptation algorithm for intensity modulated proton therapy (IMPT) based on fast Monte Carlo dose calculation and cone beam CT (CBCT) imaging. A cohort of ten head and neck cancer patients with an average of six CBCT scans were studied. To adapt the treatment plan to the new patient geometry, contours were propagated to the CBCTs with a vector field (VF) calculated with deformable image registration between the CT and the CBCTs. Within the adaptive planning algorithm, beamlets were shifted following the VF at their distal falloff and raytraced in the CBCT to adjust their energies, creating a geometrically adapted plan. Four geometric adaptation modes were studied: unconstrained geometric shifts (Free), isocenter shift (Iso), a range shifter (RS), or isocenter shift and range shifter (Iso-RS). After evaluation of the geometrical adaptation, the weights of a selected subset of beamlets were automatically tuned using MC-generated influence matrices to fulfill the original plan requirements. All beamlet calculations were done with a fast Monte Carlo running on a GPU (graphics processing unit). Geometrical adaptation alone only worked with small anatomy changes. The weight-tuned adaptation worked for every fraction, with the Free and Iso modes performing similarly and being superior than the two range shifters modes. The mean V95 and V107 were 99.4  ±  0.9 and 6.4%  ±  4.7% in the Free mode with weight tuning. The calculation time per fraction was ~5 min, but further task parallelization could reduce it to ~1–2 min for delivery adaptation right after patient setup. An online adaptation algorithm was developed that significantly improved the treatment quality for inter-fractional geometry changes. Clinical implementation of the algorithm would allow delivery adaptation right before treatment and thus allow planning margin reductions for IMPT.