The Effect of Respiratory Motion in Breast Intensity-modulated Radiation Therapy: 3D-Printed Dynamic Phantom Study.

BACKGROUND/AIM: This study evaluated the dosimetric effect of respiratory motion in intensity-modulated radiation therapy (IMRT) for breast cancer using a three-dimensional (3D)-printed dynamic phantom. MATERIALS AND METHODS: Computed tomography (CT) data from breast cancer patients were used to create a 3D-printed breast phantom. Various types of treatment plans were generated using CT images acquired at the exhalation phase. Different infill densities in the 3D-printed phantom were tested to validate their effect on simulating the average human breast tissue density. Plans were delivered to the 3D-printed dynamic phantom in the exhalation position and free-breathing motion. Dosimetric verification was performed using Gafchromic EBT3 films. RESULTS: After changing the infill density to obtain Hounsfield Unit values similar to those of human breast tissue, a realistic patient-specific breast phantom was fabricated using a 3D printer at 80% infill density. The gamma passing rates of the dose distribution delivered in the exhalation phase and free-breathing motion were 92% or more. In addition, the dynamic phantom doses with free-breathing motion were directly compared with the static phantom dose in terms of sagittal dose profiles. Gamma passing rates of >93% and 90% were achieved at 3%/3 mm and 3%/2 mm, respectively. CONCLUSION: Despite a blurred dose distribution, the dose map of the film measurement with respiratory motion could be delivered without significantly increasing the dose heterogeneity of the tumor bed or loss of target coverage. Our findings demonstrate that the impact of respiration on breast IMRT for whole-breast irradiation was not significant, even in the tumor bed.

Application of error classification model using indices based on dose distribution for characteristics evaluation of multileaf collimator position errors.

This study aims to evaluate the specific characteristics of various multileaf collimator (MLC) position errors that are correlated with the indices using dose distribution. The dose distribution was investigated using the gamma, structural similarity, and dosiomics indices. Cases from the American Association of Physicists in Medicine Task Group 119 were planned, and systematic and random MLC position errors were simulated. The indices were obtained from distribution maps and statistically significant indices were selected. The final model was determined when all values of the area under the curve, accuracy, precision, sensitivity, and specificity were higher than 0.8 (p < 0.05). The dose-volume histogram (DVH) relative percentage difference between the error-free and error datasets was examined to investigate clinical relations. Seven multivariate predictive models were finalized. The common significant dosiomics indices (GLCM Energy and GLRLM_LRHGE) can characterize the MLC position error. In addition, the finalized logistic regression model for MLC position error prediction showed excellent performance with AUC > 0.9. Furthermore, the results of the DVH were related to dosiomics analysis in that it reflects the characteristics of the MLC position error. It was also shown that dosiomics analysis could provide important information on localized dose-distribution differences in addition to DVH information.

Development of Automated Delivery Quality Assurance Analysis Software for Helical Tomotherapy.

BACKGROUND: To develop a fully automated in-house gamma analysis software for the "Cheese" phantom-based delivery quality assurance (QA) of helical tomotherapy plans. METHODS: The developed in-house software was designed to automate several procedures, which need to be manually performed using commercial software packages. The region of interest for the analysis was automatically selected by cropping out film edges and thresholding dose values (>10% of the maximum dose). The film-measured dose was automatically aligned to the computed dose using an image registration algorithm. An optimal film scaling factor was determined to maximize the percentage of pixels passing gamma (gamma passing rate) between the measured and computed doses (3%/3 mm criteria). This gamma analysis was repeated by introducing setup uncertainties in the anterior-posterior direction. For 73 tomotherapy plans, the gamma analysis results using the developed software were compared to those analyzed by medical physicists using a commercial software package. RESULTS: The developed software successfully automated the gamma analysis for the tomotherapy delivery quality assurance. The gamma passing rate (GPR) calculated by the developed software was higher than that by the clinically used software by 3.0%, on average. While, for 1 of the 73 plans, the GPR by the manual gamma analysis was higher than 90% (pass/fail criteria), the gamma analysis using the developed software resulted in fail (GPR < 90%). CONCLUSIONS: The use of automated and standardized gamma analysis software can improve both the clinical efficiency and veracity of the analysis results. Furthermore, the gamma analyses with various film scaling factors and setup uncertainties will provide clinically useful information for further investigations.

Dosimetric comparison of robust angles in carbon-ion radiation therapy for prostate cancer.

The objective of this study is to compare the plan robustness at various beam angles. Hence, the influence of the beam angles on robustness and linear energy transfer (LET) was evaluated in gantry-based carbon-ion radiation therapy (CIRT) for prostate cancer. 10 patients with prostate cancer were considered, and a total dose of 51.6 Gy (Relative biological effectiveness (RBE) was prescribed for the target volume in 12 fractions. Five beam field plans comprising two opposed fields with different angle pairs were characterized. Further, dose parameters were extracted, and the RBE-weighted dose and LET values for all angle pairs were compared. All plans considering the setup uncertainty satisfied the dose regimen. When a parallel beam pair was used for perturbed scenarios to take into account set-up uncertainty in the anterior direction, the LET clinical treatment volume (CTV) D 95% standard deviation was 1.5 times higher, and the standard deviation of RBE-weighted CTV D 95% was 7.9 times higher compared to an oblique pair. The oblique beam fields were superior in terms of dose sparing for the rectum compared to the dose distribution using two conventional lateral opposed fields for prostate cancer.

Novel framework for determining TPS-calculated doses corresponding to detector locations using 3D camera inin vivosurface dosimetry.

Purpose. To address the shortcomings of current procedures for evaluating the measured-to-planned dose agreement inin vivodosimetry (IVD), this study aimed to develop an accurate and efficient novel framework to identify the detector location placed on a patient's skin surface using a 3D camera and determine the planned dose at the same anatomical position corresponding to the detector location.Methods. Breast cancer treatment was simulated using an anthropomorphic adult female phantom (ATOM 702D; CIRS, Norfolk, VA, USA). An optically stimulated luminescent dosimeter was used for surface dose measurements (MyOSLchip, RadPro International GmbH, Germany) at six IVD points. Three-dimensional surface imaging (3DSI) of the phantom with the detector was performed in the treatment position using a 3D camera. The developed framework, iSMART, was designed to import 3DSI and treatment planning data for determining the position of the IVD detectors in the 3D treatment planning DICOM image. The clinical usefulness of iSMART was evaluated in terms of accuracy and efficiency, for comparison with the results obtained using cone-beam computed tomography (CBCT) image guidance.Results. The relative dose difference between the planned doses determined using iSMART and CBCT images displayed similar accuracies (within approximately ±2.0%) at all detector locations. The relative dose differences between the planned and measured doses at the six detector locations ranged from -4.8% to 3.1% for the CBCT images and -3.5% to 2.1% for iSMART. The total time required to read the planned doses at six detector locations averaged at 8.1 and 0.8 min for the CBCT images and iSMART, respectively.Conclusions. The proposed framework can improve the robustness of IVD analyses and aid in accurate and efficient evaluations of the measured-to-planned dose agreement.

Evaluating Mobius3D Dose Calculation Accuracy for Small-Field Flattening- Filter-Free Photon Beams.

Purpose: We aimed to investigate the dose calculation accuracy of Mobius3D for small-field flattening-filter-free x-rays, mainly utilized for stereotactic body radiation therapy (SBRT). The accuracy of beam modeling and multileaf collimator (MLC) modeling in Mobius3D, significantly affecting the dose calculation is investigated. Methods: The commissioning procedures of Mobius3D were performed for unflattened 6 MV and 10 MV x-ray beams of the linear accelerator, including beam model adjustment and dosimetric leaf gap (DLG) optimization. An experimental study with artificial plans was conducted to evaluate the accuracy of small-field modeling. The dose calculation accuracy of Mobius3D was also evaluated for retrospective SBRT plans with multiple targets. Results: Both studies evaluated the dose calculation accuracy through comparisons with the measured data. Relatively large differences were observed for off-axis distances over 5 cm and for small fields less than 1 cm field size. For the study with artificial plans, the maximum absolute error of 9.96% for unflattened 6 MV and 9.07% for unflattened 10 MV was observed when the field size was 1 cm. For the study with patient plans, the mean gamma passing rate with 3%/3 mm gamma criterion was 63.6% for unflattened 6 MV and 82.6% for unflattened 10 MV. The maximum of the average dose difference was -19.9% for unflattened 6MV and -10.1% for unflattened 10MV. Conclusions: The dose calculation accuracy uncertainties of Mobius3D for small-field flattening-filter-free photon beams were observed. The study results indicated that the beam and MLC modeling of Mobius3D must be improved for use in SBRT pretreatment QA in clinical practice.

Comprehensive clinical evaluation of TomoEQA for patient-specific pre-treatment quality assurance in helical tomotherapy.

BACKGROUND: Based on a previous study on the feasibility of TomoEQA, an exit detector-based patient-specific pre-treatment quality assurance (QA) method for helical tomotherapy, an in-depth clinical evaluation was conducted. METHODS: Data of one hundred patients were analyzed to evaluate the clinical usefulness of TomoEQA for patient-specific pre-treatment QA in comparison with the conventional phantom-based method. Additional investigations were also performed under unusual measurement conditions to validate the off-axis region. In addition to the clinical evaluation of TomoEQA, a statistical analysis was conducted to determine the plan parameters that affect the pass/failure results of pre-treatment QA. RESULTS: The average and standard deviations of the gamma passing rate and point dose error for TomoEQA were comparable to those of the conventional QA method. For TomoEQA, the average values of the gamma passing rate and point dose error were 96.32% (standard deviation (1 sigma) = 3.94; 95% confidence interval (CI), 95.55 to 97.09) and - 1.12% (standard deviation (1 sigma) = 1.04; CI, - 1.32 to - 0.92), respectively. For the conventional QA method, the average values of the gamma passing rate and point dose error were 95.95% (standard deviation (1 sigma) = 4.35; 95% confidence interval (CI), 95.10 to 96.80) and - 1.20% (standard deviation (1 sigma) = 1.61; CI, - 1.52 to - 0.88), respectively. Further experiments on the off-axis region demonstrated that TomoEQA can provide accurate results for 3D dose analysis, which is inherently difficult in the conventional QA method. Through a statistical analysis based on the results of TomoEQA, it was validated that the total fraction (Total Fx), monitor units, beam-on-time, leaf-of-time below 100 ms, and planning target volume diameter were statistically significant for the pass/failure of the pre-treatment QA results. CONCLUSIONS: TomoEQA is a clinically beneficial alternative to the conventional phantom-based QA method.

Feasibility of using dental putty-based custom molds for high-dose-rate brachytherapy of oral mucosal melanoma.

PURPOSE: To investigate the feasibility of using dental putty to fabricate custom molds for high-dose-rate (HDR) brachytherapy of oral mucosal melanoma. Practical guidelines of using dental putty for the custom mold fabrication are presented. METHODS: Dose measurements were performed using radiochromic films to understand the dosimetric characteristics of dental putty. The percentage depth dose (PDD) and profile curves were obtained using a single Co-60 source located on top of a cubic volume of the dental putty. Two custom molds were fabricated for a patient with mucosal melanoma lesions, which were spread throughout the right mandibular gingiva and soft palate regions. Furthermore, pretreatment dosimetry for both lesions was performed to evaluate the delivery quality of the resulting HDR brachytherapy plans and adjust the overall dose level. RESULTS: Prescribed doses for the two oral cavity regions were successfully delivered with a reasonable dose uniformity. Based on the measured single-source PDD curve, the maximum dose was observed at a depth of approximately 1 mm, which indicated steep dose gradients at depths ranging from 0 to 2 mm. Furthermore, a simulation with the measured single-source two-dimensional profile revealed that a source-to-source distance of 10 mm was appropriate to generate a uniform dose distribution at a source-to-surface distance of 5 mm. CONCLUSIONS: The use of a custom mold was found to be feasible for HDR brachytherapy while carefully considering the source-to-surface and source-to-source distances. Pretreatment delivery verification would be necessary when a uniform dose distribution is required.

Dosimetric outcomes of preoperative treatment planning with intraoperative optimization using stranded seeds in prostate brachytherapy.

This study aimed to evaluate the quality of low-dose-rate (LDR) prostate brachytherapy (BT) based on treatment-related dosimetric outcomes. Data of 100 patients treated using LDR BT with stranded seeds from November 2012 to November 2017 were collected. The prescription dose for the prostate was 145 Gy. The dose constraints for the preoperative plan were: V100% ≥ 95%, V150% ≤ 60%, V200% ≤ 20% for the prostate; V100% for rectum, ≤ 1 cc; and V200 Gy for urethra, 0.0 cc. Intraoperative real-time dose calculation and postoperative dose distribution analysis on days 0 and 30 were performed. Median dosimetric outcomes on days 0 and 30 respective were: V100% 92.28% and 92.23%, V200% 18.63% and 25.02%, and D90% 150.88 Gy and 151.46 Gy for the prostate; V100% for the rectum, 0.11 cc and 0.22 cc; and V200 Gy for the urethra, 0.00 cc and 0.00 cc, respectively. Twenty patients underwent additional seed implantation to compensate for insufficient dose coverage of the prostate. No loss or substantial migration of seeds or severe toxicity was reported. With stranded seed implantation and intraoperative optimization, appropriate dose delivery to the prostate without excessive dose to the organs at risk could be achieved.

Assessment of log-based fingerprinting system of Mobius3D with Elekta linear accelerators.

PURPOSE: The purpose of this study was to investigate the matching error that occurs when the Mobius3D fingerprinting system is applied in conjunction with an Elekta linear accelerator (LINAC) and to offer an acceptable and alternative method for circumventing this problem. MATERIAL AND METHODS: To avoid the multileaf collimator (MLC) conflicting error in the Mobius3D fingerprinting system, we developed an in-house program to move the MLC in the Digital Imaging and Communications in Medicine (DICOM) radiotherapy (RT)-Plan to pertinent positions, considering the relationship between log data and planned data. The re-delivered log files were calculated in the Mobius3D system, and the results were compared with those of corrected data (i.e., we analyzed a pair of re-collected log data and the previous DICOM RT-Plan data). The results were then evaluated by comparing several items, such as point dose errors, gamma index (GI) passing rates, and MLC root-mean-square (RMS) values. RESULTS: For the point dose error, the maximum difference found was below 2.0%. In the case of GI analysis of all plans, the maximum difference in the passing rates was below 1.4%. The statistical results obtained using a paired Student's t-test showed that there were no significant differences within the uncertainty. In the case of the RMS test, the maximum difference found was approximately 0.08 mm. CONCLUSIONS: Our results showed that all the mismatched log files were sufficiently acceptable within the uncertainty. We conclude that the matching error obtained when applying Mobius3D to an Elekta LINAC may be addressed using a simple modification of the fingerprinting system, and we expect that our study findings will help vendors resolve this issue in the near future.

Sensitivity of radio-photoluminescence glass dosimeters to accumulated doses.

BACKGROUND: This study investigated the effect of accumulated doses on radio-photoluminescence glass dosimeters (RPLGDs) from measurements involving mega-voltage photons. METHODS: Forty-five commercially available RPLGDs were irradiated to estimate their dose responses. Photon beams of 6, 10, and 15 MV were irradiated onto the RPLGDs inside a phantom, which were divided into five groups with different doses and energies. Groups 1 and 2 were irradiated at 1, 5, 10, 50, and 100 Gy in a sequential manner; Group 3 was irradiated 10 times with a dose of 10 Gy; and Groups 4 and 5 followed the same method as that of Group 3, but with doses of 50 Gy and 100 Gy, respectively. Each device was subjected to a measurement reading procedure each time irradiation. RESULTS: For the annealed Group 1, RPLGD exhibited a linearity response with variance within 5%. For the non-annealed Group 2, readings demonstrated hyperlinearity at 6 MV and 10 MV, and linearity at 15 MV. Following the 100 Gy irradiation, the readings for Group 2 were 118.7 ± 1.9%, 112.2 ± 2.7%, and 101.5 ± 2.3% at 6, 10, and 15 MV, respectively. For Groups 3, 4, and 5, the responsiveness of the RPLGDs gradually decreased as the number of repeated irradiations increased. The percentage readings for the 10th beam irradiation with respect to the readings for the primary beam irradiation were 84.6 ± 1.9%, 87.5 ± 2.4%, and 93.0 ± 3.0% at 6 MV, 10 MV, and 15 MV, respectively. CONCLUSIONS: The non-annealed RPLGD response to dose was hyperlinear for the 6 MV and 10 MV photon beams but not for the 15 MV photon beam. Additionally, the annealed RPLGD exhibited a fading phenomenon when the measurement was repeated several times and demonstrated a relatively large fading effect at low energies than at high energies.

TomoMQA: Automated analysis program for MVCT quality assurance of helical tomotherapy.

PURPOSE: In this study, we developed a simple but useful computer program, called TomoMQA, to offer an automated quality assurance for mega-voltage computed tomography (MVCT) images generated via helical tomotherapy. METHODS: TomoMQA is written in MATLAB and contains three steps for analysis: (a) open the DICOM dataset folder generated via helical tomotherapy (i.e., TomoTherapy® and Radixact™), (b) call the baseline data for the consistency test and click the "Analysis" button (or click the "Analysis" button without the baseline data and export the results as the baseline data), and (c) print an analyzed report. The overall procedure for the QA analysis included in TomoMQA is referred from the TG-148 recommendation. Here, the tolerances for MVCT QA were implemented from TG-148 recommended values as default; however, it can be modified by a user manually. RESULTS: To test the performance of the TomoMQA program, 15 MVCTs were prepared from five helical tomotherapy machines (1 of TomoTherapy® HD, 2 of TomoTherapy® HDA, and 2 of Radixact™) in 3 months and the QA procedures were performed using TomoMQA. From our results, the evaluation revealed that the developed program can successfully perform the MVCT QA analysis irrespective of the type of helical tomotherapy equipment. CONCLUSION: We successfully developed a new automated analysis program for MVCT QA of a helical tomotherapy platform, called TomoMQA. The developed program will be made freely downloadable from the TomoMQA-dedicated website.

Quantitative analysis of prompt gamma ray imaging during proton boron fusion therapy according to boron concentration.

The purpose of this study is to evaluate the prompt gamma ray imaging technique according to the clinical boron concentration range during proton boron fusion therapy (PBFT). To acquire a prompt gamma ray image from 32 projections, we simulated four head single photon emission computed tomography and a proton beam nozzle using a Monte Carlo simulation. We used modified ordered subset expectation maximization reconstruction algorithm with a graphic processing unit for fast image acquisition. Boron concentration was set as 20 to 100 µg at intervals of 20 µg. For quantitative analysis of the prompt gamma ray image, we acquired an image profile drawn through two boron uptake regions (BURs) and calculated the contrast value, signal-to-noise ratio (SNR), and difference between the physical target volume and volume of the prompt gamma ray image. The relative counts of prompt gamma rays were noticeably increased with increasing boron concentration. Although the intensities on the image profiles showed a similar tendency according to the boron concentration, the SNR and contrast value were improved with increasing boron concentration. This study suggests that a tumor monitoring technique using prompt gamma ray detection can be clinically applicable even if the boron concentration is relatively low.

Digital Imaging and Communications in Medicine (DICOM) information conversion procedure for SUV calculation of PET scanners with different DICOM header information.

PURPOSE: In nuclear medicine, the standardized uptake value (SUV) obtained using positron emission tomography with 2-deoxy-2-fluoro-D-glucose (FDG-PET) is widely used as a semi-quantitative diagnosis factor. We found that the header file of the Philips Allegro PET scanner using the Digital Imaging and Communications in Medicine (DICOM) standard was stored differently than with other scanners. Thus, the purpose of this study was to develop a DICOM header information conversion program to ensure compatibility between Allegro and other equipment. METHODS AND RESULTS: The NEMA IEC Body phantom was scanned using the Allegro PET scanner. We conducted measurements and performed calculations by using commercial software and the proposed self-developed program, respectively, to compare the SUVs by using conversion data. The program consists of three parts: an input part that can load data regardless of the number of DICOM images, and conversion and output parts that can be used to convert the DICOM header information and store it in the order of slices. The results of the calculation are in good agreement with the data measured at 12 circular regions of interest. The percent difference was lower than the 20%. CONCLUSION: In conclusion, this study suggested a simple and convenient method to solve the incompatibility through conversion of the DICOM header information. This study thus provides physicians more accurate information for diagnosis and treatment.

Prompt gamma ray imaging for verification of proton boron fusion therapy: A Monte Carlo study.

PURPOSE: The purpose of this study was to verify acquisition feasibility of a single photon emission computed tomography image using prompt gamma rays for proton boron fusion therapy (PBFT) and to confirm an enhanced therapeutic effect of PBFT by comparison with conventional proton therapy without use of boron. METHODS: Monte Carlo simulation was performed to acquire reconstructed image during PBFT. We acquired percentage depth dose (PDD) of the proton beams in a water phantom, energy spectrum of the prompt gamma rays, and tomographic images, including the boron uptake region (BUR; target). The prompt gamma ray image was reconstructed using maximum likelihood expectation maximisation (MLEM) with 64 projection raw data. To verify the reconstructed image, both an image profile and contrast analysis according to the iteration number were conducted. In addition, the physical distance between two BURs in the region of interest of each BUR was measured. RESULTS: The PDD of the proton beam from the water phantom including the BURs shows more efficient than that of conventional proton therapy on tumour region. A 719keV prompt gamma ray peak was clearly observed in the prompt gamma ray energy spectrum. The prompt gamma ray image was reconstructed successfully using 64 projections. Different image profiles including two BURs were acquired from the reconstructed image according to the iteration number. CONCLUSION: We confirmed successful acquisition of a prompt gamma ray image during PBFT. In addition, the quantitative image analysis results showed relatively good performance for further study.