A predicted three-dimensional dose sequence based treatment planning optimization method for gynecologic IMRT.

In knowledge-based treatment planning (KBTP) for intensity-modulated radiation therapy (IMRT), the quality of the plan is dependent on the sophistication of the predicted dosimetric information and its application. In this paper, we propose a KBTP method that based on the effective and reasonable utilization of a three-dimensional (3D) dose prediction on planning optimization. We used an organs-at-risk (OARs) dose distribution prediction model to create a voxel-based dose sequence based optimization objective for OARs doses. This objective was used to reformulate a traditional fluence map optimization model, which involves a tolerable spatial re-assignment of the predicted dose distribution to the OAR voxels based on their current doses' positions at a sorted dose sequencing. The feasibility of this method was evaluated with ten gynecology (GYN) cancer IMRT cases by comparing its generated plan quality with the original clinical plan. Results showed feasible plan by proposed method, with comparable planning target volume (PTV) dose coverage and greater dose sparing of the OARs. Among ten GYN cases, the average V30 and V45 of rectum were decreased by 4%±4% (p = 0.02) and 4%±3% (p<0.01), respectively. V30 and V45 of bladder were decreased by 8%±2% (p<0.01) and 3%±2% (p<0.01), respectively. Our predicted dose sequence-based planning optimization method for GYN IMRT offered a flexible use of predicted 3D doses while ensuring the output plan consistency.

Radiation pneumonitis prediction after stereotactic body radiation therapy based on 3D dose distribution: dosiomics and/or deep learning-based radiomics features.

BACKGROUND: This study was designed to establish radiation pneumonitis (RP) prediction models using dosiomics and/or deep learning-based radiomics (DLR) features based on 3D dose distribution. METHODS: A total of 140 patients with non-small cell lung cancer who received stereotactic body radiation therapy (SBRT) were retrospectively included in this study. These patients were randomly divided into the training (n = 112) and test (n = 28) sets. Besides, 107 dosiomics features were extracted by Pyradiomics, and 1316 DLR features were extracted by ResNet50. Feature visualization was performed based on Spearman's correlation coefficients, and feature selection was performed based on the least absolute shrinkage and selection operator. Three different models were constructed based on random forest, including (1) a dosiomics model (a model constructed based on dosiomics features), (2) a DLR model (a model constructed based on DLR features), and (3) a hybrid model (a model constructed based on dosiomics and DLR features). Subsequently, the performance of these three models was compared with receiver operating characteristic curves. Finally, these dosiomics and DLR features were analyzed with Spearman's correlation coefficients. RESULTS: In the training set, the area under the curve (AUC) of the dosiomics, DLR, and hybrid models was 0.9986, 0.9992, and 0.9993, respectively; the accuracy of these three models was 0.9643, 0.9464, and 0.9642, respectively. In the test set, the AUC of these three models was 0.8462, 0.8750, and 0.9000, respectively; the accuracy of these three models was 0.8214, 0.7857, and 0.8571, respectively. The hybrid model based on dosiomics and DLR features outperformed other two models. Correlation analysis between dosiomics features and DLR features showed weak correlations. The dosiomics features that correlated DLR features with the Spearman's rho |ρ| ≥ 0.8 were all first-order features. CONCLUSION: The hybrid features based on dosiomics and DLR features from 3D dose distribution could improve the performance of RP prediction after SBRT.

Extending in aqua portal dosimetry with dose inhomogeneity conversion maps for accurate patient dose reconstruction in external beam radiotherapy.

Background and purpose: In aqua dosimetry with electronic portal imaging devices (EPIDs) allows for dosimetric treatment verification in external beam radiotherapy by comparing EPID-reconstructed dose distributions (EPID_IA) with dose distributions calculated with the treatment planning system in water-equivalent geometries. The main drawback of the method is the inability to estimate the dose delivered to the patient. In this study, an extension to the method is presented to allow for patient dose reconstruction in the presence of inhomogeneities. Materials and methods: EPID_IA dose distributions were converted into patient dose distributions (EPID_IA_MC) by applying a 3D dose inhomogeneity conversion, defined as the ratio between patient and water-filled patient dose distributions computed using Monte Carlo calculations. EPID_IA_MC was evaluated against dose distributions calculated with a collapsed cone convolution superposition (CCCS) algorithm and with a GPU-based Monte Carlo dose calculation platform (GPUMCD) using non-transit EPID measurements of 25 plans. In vivo EPID measurements of 20 plans were also analyzed. Results: In the evaluation of EPID_IA_MC, the average γ-mean values (2% local/2mm, 50% isodose volume) were 0.70 ± 0.14 (1SD) and 0.66 ± 0.10 (1SD) against CCCS and GPUMCD, respectively. Percentage differences in median dose to the planning target volume were within 3.9% and 2.7%, respectively. The number of in vivo dosimetric alerts with EPID_IA_MC was comparable to EPID_IA. Conclusions: EPID_IA_MC accommodates accurate patient dose reconstruction for treatment disease sites with significant tissue inhomogeneities within a simple EPID-based direct dose back-projection algorithm, and helps to improve the clinical interpretation of both pre-treatment and in vivo dosimetry results.

Plastic scintillation dosimeter with a conical mirror for measuring 3D dose distribution.

PURPOSE: To test the measurement technique of the three-dimensional (3D) dose distribution measured image by capturing the scintillation light generated using a plastic scintillator and a scintillating screen. METHODS: Our imaging system constituted a column shaped plastic scintillator covered by a Gd2 O2 S:Tb scintillating screen, a conical mirror and a cooled CCD camera. The scintillator was irradiated with 6 MV photon beams. Meanwhile, the irradiated plan was prepared for the static field plans, two-field plan (2F plan) and the conformal arc plan (CA plan). The 2F plan contained 16 mm2 and 10 mm2 fields irradiated from gantry angles of 0° and 25°, respectively. The gantry was rotated counterclockwise from 45° to 315° for the CA plan. The field size was then obtained as 10 mm2 . A Monte Carlo simulation was performed in the experimental geometry to obtain the calculated 3D dose distribution as the reference data. Dose response was acquired by comparing between the reference and the measurement. The dose rate dependence was verified by irradiating the same MU value at different dose rates ranging from 100 to 600 MU/min. Deconvolution processing was applied to the measured images for the correction of light blurring. The measured 3D dose distribution was reconstructed from each measured image. Gamma analysis was performed to these 3D dose distributions. The gamma criteria were 3% for the dose difference, 2 mm for the distance-to-agreement and 10% for the threshold. RESULTS: Dose response for the scintillation light was linear. The variation in the light intensity for the dose rate ranging from 100 to 600 MU/min was less than 0.5%, while our system presents dose rate independence. For the 3D dose measurement, blurring of light through deconvolution processing worked well. The 3D gamma passing rate (3D GPR) for the 10 × 10 mm2 , 16 × 16 mm2 , and 20 × 20 mm2 fields were observed to be 99.3%, 98.8%, and 97.8%, respectively. Reproducibility of measurement was verified. The 3D GPR results for the 2F plan and the CA plan were 99.7% and 100%, respectively. CONCLUSIONS: We developed a plastic scintillation dosimeter and demonstrated that our system concept can act as a suitable technique for measuring the 3D dose distribution from the gamma results. In the future, we will attempt to measure the 4D dose distribution for clinical volumetric modulated arc radiation therapy (VMAT)-SBRTplans.

Three-dimensional dose-distribution measurement of therapeutic carbon-ion beams using a ZnS scintillator sheet.

The accurate measurement of the 3D dose distribution of carbon-ion beams is essential for safe carbon-ion therapy. Although ionization chambers scanned in a water tank or air are conventionally used for this purpose, these measurement methods are time-consuming. We thus developed a rapid 3D dose-measurement tool that employs a silver-activated zinc sulfide (ZnS) scintillator with lower linear energy transfer (LET) dependence than gadolinium-based (Gd) scintillators; this tool enables the measurement of carbon-ion beams with small corrections. A ZnS scintillator sheet was placed vertical to the beam axis and installed in a shaded box. Scintillation images produced by incident carbon-ions were reflected with a mirror and captured with a charge-coupled device (CCD) camera. A 290 MeV/nucleon mono-energetic beam and spread-out Bragg peak (SOBP) carbon-ion passive beams were delivered at the Gunma University Heavy Ion Medical Center. A water tank was installed above the scintillator with the water level remotely adjusted to the measurement depth. Images were recorded at various water depths and stacked in the depth direction to create 3D scintillation images. Depth and lateral profiles were analyzed from the images. The ZnS-scintillator-measured depth profile agreed with the depth dose measured using an ionization chamber, outperforming the conventional Gd-based scintillator. Measurements were realized with smaller corrections for a carbon-ion beam with a higher LET than a proton. Lateral profiles at the entrance and the Bragg peak depths could be measured with this tool. The proposed method would make it possible to rapidly perform 3D dose-distribution measurements of carbon-ion beams with smaller quenching corrections.

Dosiomics improves prediction of locoregional recurrence for intensity modulated radiotherapy treated head and neck cancer cases.

OBJECTIVES: To investigate whether dosiomics can benefit to IMRT treated patient's locoregional recurrences (LR) prediction through a comparative study on prediction performance inspection between radiomics methods and that integrating dosiomics in head and neck cancer cases. MATERIALS AND METHODS: A cohort of 237 patients with head and neck cancer from four different institutions was obtained from The Cancer Imaging Archive and utilized to train and validate the radiomics-only prognostic model and integrate the dosiomics prognostic model. For radiomics, the radiomics features were initially extracted from images, including CTs and PETs, and selected on the basis of their concordance index (CI) values, then condensed via principle component analysis. Lastly, multivariate Cox proportional hazards regression models were constructed with class-imbalance adjustment as the LR prediction models by inputting those condensed features. For dosiomics integration model establishment, the initial features were similar, but with additional 3-dimensional dose distribution from radiation treatment plans. The CI and the Kaplan-Meier curves with log-rank analysis were used to assess and compare these models. RESULTS: Observed from the independent validation dataset, the CI of the model for dosiomics integration (0.66) was significantly different from that for radiomics (0.59) (Wilcoxon test, p=5.9×10-31). The integrated model successfully classified the patients into high- and low-risk groups (log-rank test, p=2.5×10-02), whereas the radiomics model was not able to provide such classification (log-rank test, p=0.37). CONCLUSION: Dosiomics can benefit in predicting the LR in IMRT-treated patients and should not be neglected for related investigations.

New calculation method for 3D dose distribution in tetrahedral-mesh phantoms in Geant4.

The tetrahedral-mesh (TM) geometry, which is a very promising geometry for computational human phantoms, has a limitation in 3D dose distribution calculation for medical applications. Even though Geant4 provides the read-out geometry for calculating 3D dose distribution in the TM geometry, this method significantly slows down the computation speed. In the present study, we developed a new method, called Moving Voxel-based Dose-Distribution Calculator (MVDDC), to rapidly calculate a 3D dose distribution in a TM geometry. To evaluate the performance of the MVDDC method, a simple TM cubic phantom and a human phantom were implemented in Geant4. Subsequently, the phantoms were irradiated with proton spot beams under various conditions, and the obtained results were compared with those of the read-out geometry method. The results show that there is no significant difference between the dose distributions calculated using the new method and the read-out geometry method. With respect to the computational performance, the speeds of simulations using the MVDDC were approximately 1.4-2.7 times faster than those of the simulations using the read-out geometry method.

Texture analysis of 3D dose distributions for predictive modelling of toxicity rates in radiotherapy.

BACKGROUND AND PURPOSE: To explore the use of texture analysis (TA) features of patients' 3D dose distributions to improve prediction modelling of treatment complication rates in prostate cancer radiotherapy. MATERIAL AND METHODS: Late toxicity scores, dose distributions, and non-treatment related (NTR) predictors for late toxicity, such as age and baseline symptoms, of 351 patients of the hypofractionation arm of the HYPRO randomized trial were used in this study. Apart from DVH parameters, also TA features of rectum and bladder 3D dose distributions were used for predictive modelling of gastrointestinal (GI) and genitourinary (GU) toxicities. Logistic Normal Tissue Complication Probability (NTCP) models were derived, using only NTR parameters, NTR + DVH, NTR + TA, and NTR + DVH + TA. RESULTS: For rectal bleeding, the area under the curve (AUC) for using only NTR parameters was 0.58, which increased to 0.68, and 0.73, when adding DVH or TA parameters respectively. For faecal incontinence, the AUC went up from 0.63 (NTR only), to 0.68 (+DVH) and 0.73 (+TA). For nocturia, adding TA features resulted in an AUC increase from 0.64 to 0.66, while no improvement was seen when including DVH parameters in the modelling. For urinary incontinence, the AUC improved from 0.68 to 0.71 (+DVH) and 0.73 (+TA). For GI, model improvements resulting from adding TA parameters to NTR instead of DVH were statistically significant (p < 0.04). CONCLUSION: Inclusion of 3D dosimetric texture analysis features in predictive modelling of GI and GU toxicity rates in prostate cancer radiotherapy improved prediction performance, which was statistically significant for GI.

[Fricke Gel Dosimeters].

Fricke gel dosimeters are based on the aqueous ferrous sulfate solution that has been used as a reliable chemical dosimeter for more than 90 years. The amount of radiation-induced oxidation of ferrous ions to ferric ions could be evaluated three-dimensionally via optical computed tomography (OCT) or magnetic resonance imaging (MRI). Three-dimensional dosimetry is expected to be a useful tool, in particular, for the verification of complex radiation dose planning in advanced radiation cancer therapy. Compared with other 3D gel dosimeters, there are several problems such as retention of dose distribution or dose sensitivity; however, they can easily be prepared in a chemistry laboratory. In addition, a unique gel dosimeter was also reported for heavy-ion beam irradiation. In this review, recent papers concerning the Fricke gel dosimeter and its application in 3D dosimetry are summarized.

[Polymer Gel Dosimeter].

With rapid advances being made in radiotherapy treatment, three-dimensional (3D) dose measurement techniques of great precision are required more than ever before. It is expected that 3D polymer gel dosimeters will satisfy clinical needs for an effective detector that can measure the complex 3D dose distributions. Polymer gel dosimeters are devices that utilize the radiation-induced polymerization reactions of vinyl monomers in a gel to store information about radiation dose. The 3D absorbed dose distribution can be deduced from the resulting polymer distribution using several imaging modalities, such as MRI, X-ray and optical CTs. In this article, the fundamental characteristics of polymer gel dosimeter are reviewed and some challenging keys are also suggested for the widely spread in clinical use.

[Evaluation of Radiochromic Gel Dosimeter and Optical CT].

Three-dimensional (3D) radiation dosimeter has received growing interest with the implementation of highly conformal radiotherapy treatments. Optical Computed tomography (CT) scanners have been developed for 3D quantitative measurement of optical attenuation. Malachite green and Leuco Crystal Violet (LCV) dye radiochromic gels were recently developed as radiation dosimeters in combination with optical CT scanners. Malachite green and LCV are colorless dyes and they become green or violet after irradiation, respectively. The dyes mixed with gelatin and surfactant became hardened and can be used as a 3D gel dosimeter. The radiochromic gels dosimeter with optical CT scanners provides many of desired features such as: low cost, easy to fabricate, low toxicity and fast readout. In this article, the fundamental characteristics of radiochromic gel dosimeter with optical CT scanners are reviewed.

Adaptation, Commissioning, and Evaluation of a 3D Treatment Planning System for High-Resolution Small-Animal Irradiation.

Although spatially precise systems are now available for small-animal irradiations, there are currently limited software tools available for treatment planning for such irradiations. We report on the adaptation, commissioning, and evaluation of a 3-dimensional treatment planning system for use with a small-animal irradiation system. The 225-kV X-ray beam of the X-RAD 225Cx microirradiator (Precision X-Ray) was commissioned using both ion-chamber and radiochromic film for 10 different collimators ranging in field size from 1 mm in diameter to 40 × 40 mm(2) A clinical 3-dimensional treatment planning system (Metropolis) developed at our institution was adapted to small-animal irradiation by making it compatible with the dimensions of mice and rats, modeling the microirradiator beam orientations and collimators, and incorporating the measured beam data for dose calculation. Dose calculations in Metropolis were verified by comparison with measurements in phantoms. Treatment plans for irradiation of a tumor-bearing mouse were generated with both the Metropolis and the vendor-supplied software. The calculated beam-on times and the plan evaluation tools were compared. The dose rate at the central axis ranges from 74 to 365 cGy/min depending on the collimator size. Doses calculated with Metropolis agreed with phantom measurements within 3% for all collimators. The beam-on times calculated by Metropolis and the vendor-supplied software agreed within 1% at the isocenter. The modified 3-dimensional treatment planning system provides better visualization of the relationship between the X-ray beams and the small-animal anatomy as well as more complete dosimetric information on target tissues and organs at risk. It thereby enhances the potential of image-guided microirradiator systems for evaluation of dose-response relationships and for preclinical experimentation generally.

OAR predicted dose distribution and gEUD based treatment planning optimization for IMRT / 中华放射医学与防护杂志

Objective To propose a treatment planning optimization algorithm which can make full use of OAR dose distribution prediction meanwhile improving the output planning quality as much as possible.Methods We had reformulated an FMO function under the guidance of dose distribution prediction and also integrated equivalent uniform dose (gEUD) based on the consideration of prediction uncertainty,for providing optimal solution.Performance of the method was evaluated by comparing the optimized IMRT plan quality of 8 cervical cancers in the term of DVH curves,dose distribution and dosimetric endpoints with the original ones.Results The proposed method had a feasible,fast solution.Compared with original plan,its output plan had better plan quality in better dose homogeneity,less hot spot and further dose sparing for OARs.V30,V45 of rectum was decreased by (6.60±3.53)％ and (17.03±7.44)％,respectively,with the statistically significant difference (t=-4.954,-6.055,P＜0.05).V30,V45 of bladder was decreased by (14.74 ± 5.61) ％ and (14.99 ± 4.53) ％,respectively,with the statistically significant difference (t=-6.945,-8.759,P＜0.05).Conclusions We have successfully developed a predicted dose distribution and equivalent uniform dose-based planning optimization method,which is able to make good use of 3D dose prediction and ensure the output plan quality for intensity modulated radiation therapy.

3-Dimensional Dosimetry of Small Field Photon Beam / 의학물리

A polymer gel dosimeter was fabricated. A 3-dimensional dosimetry experiment was performed in the small field of the photon of the cyberknife. The dosimeter was installed in a head and neck phantom. It was manufactured from the acrylic and it was used in dosimetry. By using the head and neck CT protocol of the CyberKnife system, CT images of the head and neck phantom were obtained and delivered to the treatment planning system. The irradiation to the dosimeter in the treatment planning was performed, and then, the image was obtained by using 3.0T magnetic resonance imaging (MRI) after 24 hours. The dose distribution of the phantom was analyzed by using MATLAB. The results of this measurement were compared to the results of calculation in the treatment planning. In the isodose curve on the axial direction, the dose distribution coincided with the high dose area, 0.76mm difference on 80%, rather than the low dose area, 1.29 mm difference on 40%. In this research, the fact that the polymer gel dosimeter and MRI can be applied for analyzing a small field in a 3 dimensional dosimetry was confirmed. Moreover, the feasibility of using these for the therapeutic radiation quality control was also confirmed.