Computational analysis of interfractional anisotropic shape variations of the rectum in prostate cancer radiation therapy

Reducing systematic errors due to deformation of organs at risk in radiotherapy

PURPOSE

In radiotherapy (RT), the planning CT (pCT) is commonly used to plan the full RT-course. Due to organ deformation and motion, the organ shapes seen at the pCT will not be identical to their shapes during RT. Any difference between the pCT organ shape and the organ's mean shape during RT will cause systematic errors. We propose to use statistical shrinkage estimation to reduce this error using only the pCT and the population mean shape computed from training data.

METHODS

The method was evaluated for the rectum in a cohort of 37 prostate cancer patients that had a pCT and 7-10 treatment CTs with rectum delineations. Deformable registration was performed both within-patient and between patients, resulting in point-to-point correspondence between all rectum shapes, which enabled us to compute a population mean rectum. Shrinkage estimates were found by combining the pCTs linearly with the population mean. The method was trained and evaluated using leave-one-out cross validation. The shrinkage estimates and the patient mean shapes were compared geometrically using the Dice similarity index (DSI), Hausdorff distance (HD), and bidirectional local distance. Clinical dose/volume histograms, equivalent uniform dose (EUD) and minimum dose to the hottest 5% volume (D5%) were compared for the shrinkage estimate and the pCT.

RESULTS

The method resulted in moderate but statistically significant increase in similarity to the patient mean shape over the pCT. On average, the HD was reduced from 15.6 to 13.4 mm, while the DSI was increased from 0.74 to 0.78. Significant reduction in the bias of volume estimates was found in the DVH-range of 52.5-65 Gy, where the bias was reduced from -1.3 to -0.2 percentage points, but no significant improvement was found in EUD or D5%, CONCLUSIONS: The results suggest that shrinkage estimation can reduce systematic errors due to organ deformations in RT. The method has potential to increase the accuracy in RT of deformable organs and can improve motion modeling.

Interfractional robustness of scanning carbon ion radiotherapy for prostate cancer: An analysis based on dose distribution from daily in-room CT images

Purpose

We analyzed interfractional robustness of scanning carbon ion radiotherapy (CIRT) for prostate cancer based on the dose distribution using daily in‐room computed tomography (CT) images.

Materials and Methods

We analyzed 11 consecutive patients treated with scanning CIRT for localized prostate cancer in our hospital between December 2015 and January 2016. In‐room CT images were taken under treatment conditions in every treatment session. The dose distribution on each in‐room CT image was recalculated, while retaining the pencil beam arrangement of the initial treatment plan. Then, the dose–volume histogram (DVH) parameters including the percentage of the clinical target volume (CTV) with 95% and 90% of the prescribed dose area (V95% of CTV, V90% of CTV) and V80% of rectum were calculated. The acceptance criteria for the CTV and rectum were set at V95% of CTV ≥95%, V90% of CTV ≥98%, and V80% of rectum < 10 ml.

Results

V95% of CTV, V90% of CTV, and V80% of rectum for the reproduced plans were 98.8 ± 3.49%, 99.5 ± 2.15%, and 4.39 ± 3.96 ml, respectively. Acceptance of V95% of CTV, V90% of CTV, and V80% of rectum was obtained in 123 (94%), 125 (95%) and 117 sessions (89%), respectively. Acceptance of the mean dose of V95% of CTV, V90% of CTV, and V80% of rectum for each patient was obtained in 10 (91%), 10 (91%), and 11 patients (100%), respectively.

Conclusions

We demonstrated acceptable interfractional robustness based on the dose distribution in scanning CIRT for prostate cancer.

Analyses of integrated EPID images for on-treatment quality assurance to account for interfractional variations in volumetric modulated arc therapy

Purpose

To investigate the effects of interfractional variation, such as anatomical changes and setup errors, on dose delivery during treatment for prostate cancer (PC) and head and neck cancer (HNC) by courses of volumetric modulated arc therapy (VMAT) aided by on‐treatment electronic portal imaging device (EPID) images.

Methods

Seven patients with PC and 20 patients with HNC who had received VMAT participated in this study. After obtaining photon fluence at the position of the EPID for each treatment arc from on‐treatment integrated EPID images, we calculated the differences between the fluence for the first fraction and each subsequent fraction for each arc. The passing rates were investigated based on a tolerance level of 3% of the maximum fluence during the treatment courses and the correlations between the passing rates and anatomical changes.

Results

In PC, the median and lowest passing rates were 99.8% and 95.2%, respectively. No correlations between passing rates and interfractional variation were found. In HNC, the median passing rate of all fractions was 93.0%, and the lowest passing rate was 79.6% during the 35th fraction. Spearman’s correlation coefficients between the passing rates and changes in weight or neck volume were − 0.77 and − 0.74, respectively.

Conclusions

Analyses of the on‐treatment EPID images facilitates estimates of the interfractional anatomical variation in HNC patients during VMAT and thus improves assessments of the need for re‐planning or adaptive strategies and the timing thereof.

Comparison of volumetric-modulated arc therapy and intensity-modulated radiation therapy prostate cancer plans accounting for cold spots

This study compared dosimetric indices of volumetric-modulated arc therapy (VMAT) with intensity-modulated radiation therapy (IMRT) accounting for cold spots in prostate cancer plans. IMRT plans were retrospectively generated from 30 prostate cancer patients with ten cases for each risk group, who received VMAT plans. The mean, maximum, and minimum doses, and conformity and homogeneity indexes were evaluated for planning target volume (PTV) and the mean dose and V20-V70 for organs at risk (OAR) including the rectum, bladder, right and left femoral heads, and rectum overlapped with PTV (ROP) regions. The numbers and volume percentages of cold spots within PTVs and ROP regions were measured using in-house software. Three-dimensional probabilistic distributions of the probability and distributions of cold spots were generated using a centroid matching technique for visualization and analysis. There was a statistically better dose conformity in the PTV, rectum, and bladder dose-sparing in VMAT compared to that in the IMRT plans, whereas VMAT had statistically worse target dose homogeneity, and right and left femoral head dose-sparing than those of the IMRT plans. The average volume percentage of cold spots per PTV for the VMAT was 4.37 ± 2.68%, which was smaller than the 5.72 ± 1.84% observed for IMRT plans (P = 0.007). The volume percentage of cold spots per ROP for the VMAT did not significantly differ from those for the IMRT plans. Compared with IMRT, the VMAT plans achieved better PTV dose conformity, OAR dose-sparing, and smaller cold spots in the treatment of prostate cancer.

Rotationally resliced 3D prostate segmentation of MR images using Bhattacharyya similarity and active band theory

PURPOSE

In this article, we propose a novel, semi-automatic segmentation method to process 3D MR images of the prostate using the Bhattacharyya coefficient and active band theory with the goal of providing technical support for computer-aided diagnosis and surgery of the prostate.

METHODS

Our method consecutively segments a stack of rotationally resectioned 2D slices of a prostate MR image by assessing the similarity of the shape and intensity distribution in neighboring slices. 2D segmentation is first performed on an initial slice by manually selecting several points on the prostate boundary, after which the segmentation results are propagated consecutively to neighboring slices. A framework of iterative graph cuts is used to optimize the energy function, which contains a global term for the Bhattacharyya coefficient with the help of an auxiliary function. Our method does not require previously segmented data for training or for building statistical models, and manual intervention can be applied flexibly and intuitively, indicating the potential utility of this method in the clinic.

RESULTS

We tested our method on 3D T2-weighted MR images from the ISBI dataset and PROMISE12 dataset of 129 patients, and the Dice similarity coefficients were 90.34 ± 2.21% and 89.32 ± 3.08%, respectively. The comparison was performed with several state-of-the-art methods, and the results demonstrate that the proposed method is robust and accurate, achieving similar or higher accuracy than other methods without requiring training.

CONCLUSION

The proposed algorithm for segmenting 3D MR images of the prostate is accurate, robust, and readily applicable to a clinical environment for computer-aided surgery or diagnosis.

MRI-guided localization of the dominant intraprostatic lesion and dose analysis of volumetric modulated arc therapy planning for prostate cancer

PURPOSE

Primary radiation therapy is a curative treatment option for prostate cancer. The aim of this study was to evaluate the detection of the dominant intraprostatic lesion (DIL) with magnetic resonance imaging (MRI) for radiotherapy treatment planning, the comparison with transrectal ultrasound (TRUS)-guided biopsies and the examination of the dose distribution in relation to the DIL location.

MATERIALS AND METHODS

In all, 54 patients with treatment planning MRI for primary radiotherapy of prostate cancer from 03/2015 to 03/2017 at the Universitätsklinikum Würzburg were identified. The localization of the DIL was based on MRI with T2- and diffusion-weighted imaging. After registration of the MR image sets within Pinnacle³ (Philips Radiation Oncology Systems, Fitchburg, WI, USA), the dose distribution was analyzed. The location of the DIL was compared to the pathology reports in a side-based manner.

RESULTS

The DIL mean dose (Dmean) was 77.51 ± 0.77 Gy and in 50/51 cases within the tolerance range or exceeded the prescribed dose. There was a significant difference in Dmean between ventral (n = 21) and dorsal (n = 30) DIL (77.87 ± 0.67 vs. 77.26 ± 0.77 Gy; p = 0.005). MRI-guided localization showed an accuracy and sensitivity of up to 78.8% and 82.1% for inclusion of secondary lesions, respectively.

CONCLUSION

Up to 82.1% of histologically verified intraprostatic lesions were identified in the context of MRI-guided radiotherapy treatment planning. As expected, dorsal DIL tend to be minimally underdosed in comparison to ventral DIL. Adequate dose coverage was achieved in over 98% of patients.