Kiyokazu Sato

Evaluation of patient DVH‐based QA metrics for prostate VMAT: correlation between accuracy of estimated 3D patient dose and magnitude of MLC misalignment

The purpose of this study was to evaluate the accuracy of commercially available software, using patient DVH‐based QA metrics, by investigating the correlation between estimated 3D patient dose and magnitude of MLC misalignments. We tested 3DVH software with an ArcCHECK. Two different calculating modes of ArcCHECK Planned Dose Perturbation (ACPDP) were used: “Normal Sensitivity” and “High Sensitivity”. Ten prostate cancer patients treated with hypofractionated VMAT (67.6 Gy/26 Fr) in our hospital were studied. For the baseline plan, we induced MLC errors (−0.75,−0.5,−0.25,0.25,0.5, and 0.75 mm for each single bank). We calculated the dose differences between the ACPDP dose with error and TPS dose with error using gamma passing rates and using DVH‐based QA metrics. The correlations between dose estimation error and MLC position error varied with each structure and metric. A comparison using 1%/1 mm gamma index showed that the larger was the MLC error‐induced, the worse were the gamma passing rates. Slopes of linear fit to dose estimation error versus MLC position error for mean dose and D95 to the PTV were 1.76 and 1.40% mm−1, respectively, for “Normal Sensitivity”, and −0.53 and 0.88% mm−1, respectively, for “High Sensitivity”, showing better accuracy for “High Sensitivity” than “Normal Sensitivity”. On the other hand, the slopes for mean dose to the rectum and bladder, V35 to the rectum and bladder and V55 to the rectum and bladder, were −1.00,−0.55,−2.56,−1.25,−3.53, and 1.85% mm−1, respectively, for “Normal Sensitivity”, and −2.89,−2.39,−4.54,−3.12,−6.24, and −4.11% mm−1, respectively, for “High Sensitivity”, showing significant better accuracy for “Normal Sensitivity” than “High Sensitivity”. Our results showed that 3DVH had some residual error for both sensitivities. Furthermore, we found that “Normal Sensitivity” might have better accuracy for the DVH metric for the PTV and that “High Sensitivity” might have better accuracy for DVH metrics for the rectum and bladder. We must be willing to tolerate this residual error in clinical care. PACS number: 87.55Qr

Comparison of DVH‐based plan verification methods for VMAT: ArcCHECK‐3DVH system and dynalog‐based dose reconstruction

Abstract The purpose of this study was comparing dose‐volume histogram (DVH)‐based plan verification methods for volumetric modulated arc therapy (VMAT) pretreatment QA. We evaluated two 3D dose reconstruction systems: ArcCHECK‐3DVH system (Sun Nuclear corp.) and Varian dynalog‐based dose reconstruction (DBDR) system, developed in‐house. Fifteen prostate cancer patients (67.6 Gy/26 Fr), four head and neck cancer patient (66 Gy/33 Fr), and four esophagus cancer patients (60 Gy/30 Fr) treated with VMAT were studied. First, ArcCHECK measurement was performed on all plans; simultaneously, the Varian dynalog data sets that contained the actual delivered parameters (leaf positions, gantry angles, and cumulative MUs) were acquired from the Linac control system. Thereafter, the delivered 3D patient dose was reconstructed by 3DVH software (two different calculating modes were used: High Sensitivity (3DVH‐HS) and Normal Sensitivity (3DVH‐NS)) and in‐house DBDR system. We evaluated the differences between the TPS‐calculated dose and the reconstructed dose using 3D gamma passing rates and DVH dose index analysis. The average 3D gamma passing rates (3%/3 mm) between the TPS‐calculated dose and the reconstructed dose were 99.1 ± 0.6%, 99.7 ± 0.3%, and 100.0 ± 0.1% for 3DVH–HS, 3DVH–NS, and DBDR, respectively. For the prostate cases, the average differences between the TPS‐calculated dose and reconstructed dose in the PTV mean dose were 1.52 ± 0.50%, −0.14 ± 0.55%, and −0.03 ± 0.07% for 3DVH–HS, 3DVH–NS, and DBDR, respectively. For the head and neck and esophagus cases, the dose difference to the TPS‐calculated dose caused by an effect of heterogeneity was more apparent under the 3DVH dose reconstruction than the DBDR. Although with some residual dose reconstruction errors, these dose reconstruction methods can be clinically used as effective tools for DVH‐based QA for VMAT delivery.

Evaluation of deformable image registration between external beam radiotherapy and HDR brachytherapy for cervical cancer with a 3D‐printed deformable pelvis phantom

Purpose In this study, we developed a 3D‐printed deformable pelvis phantom for evaluating spatial DIR accuracy. We then evaluated the spatial DIR accuracies of various DIR settings for cervical cancer. Methods A deformable female pelvis phantom was created based on patient CT data using 3D printing. To create the deformable uterus phantom, we first 3D printed both a model of uterus and a model of the internal cavities of the vagina and uterus. We then made a mold using the 3D printed uterus phantom. Finally, urethane was poured into the mold with the model of the internal cavities in place, creating the deformable uterus phantom with a cavity into which an applicator could be inserted. To create the deformable bladder phantom, we first 3D printed models of the bladder and of the same bladder scaled down by 2 mm. We then made a mold using the larger bladder model. Finally, silicone was poured into the mold with the smaller bladder model in place to create the deformable bladder phantom with a wall thickness of 2 mm. To emulate the anatomical bladder, water was poured into the created bladder. We acquired phantom image without applicator for EBRT. Then, we inserted the applicator into the phantom to simulate BT. In this situation, we scanned the phantom again to obtain the phantom image for BT. We performed DIR using the two phantom images in two cases: Case A, with full bladder (170 ml) in both EBRT and BT images; and Case B with full bladder in the BT image and half‐full bladder (100 ml) in the EBRT image. DIR was evaluated using Dice similarity coefficients (DSCs) and 31 landmarks for the uterus and 25 landmarks for the bladder. A hybrid intensity and structure DIR algorithm implemented in RayStation with four DIR settings was evaluated. Results On visual inspection, reasonable agreement in shape of the uterus between the phantom and patient CT images was observed for both EBRT and BT, although some regional disagreements in shape of the bladder and rectum were apparent. The created phantom could reproduce the actual patients uterus deformation by the applicator. For both Case A and B, large variation was seen in landmark error among the four DIR parameters. In addition, although DSCs were comparable, moderate differences in landmark error existed between the two different DIR parameters selected from the four DIR parameters (i.e., DSC = 0.96, landmark error = 13.2 ± 5.7 mm vs. DSC = 0.97, landmark error = 9.7 ± 4.0 mm). This result suggests that landmark error evaluation might thus be more effective than DSC for evaluating DIR accuracy. Conclusions Our developed phantom enabled the evaluation of spatial DIR accuracy for the female pelvic region for the first time. Although the DSCs are high, the spatial errors can still be significant and our developed phantom facilitates their quantification. Our results showed that optimization is needed to identify suitable DIR settings. For determining suitable DIR settings, our method of evaluating spatial DIR accuracy using the 3D‐printed phantom may prove helpful.

Evaluation of the performance of deformable image registration between planning CT and CBCT images for the pelvic region: comparison between hybrid and intensity-based DIR

Abstract This study aimed to evaluate the performance of the hybrid deformable image registration (DIR) method in comparison with intensity-based DIR for pelvic cone-beam computed tomography (CBCT) images, using intensity and anatomical information. Ten prostate cancer patients treated with intensity-modulated radiation therapy (IMRT) were studied. Nine or ten CBCT scans were performed for each patient. First, rigid registration was performed between the planning CT and all CBCT images using gold fiducial markers, and then DIR was performed. The Dice similarity coefficient (DSC) and center of mass (COM) displacement were used to evaluate the quantitative DIR accuracy. The average DSCs for intensity-based DIR for the prostate, rectum, bladder, and seminal vesicles were 0.84 ± 0.05, 0.75 ± 0.05, 0.69 ± 0.07 and 0.65 ± 0.11, respectively, whereas those values for hybrid DIR were 0.98 ± 0.00, 0.97 ± 0.01, 0.98 ± 0.00 and 0.94 ± 0.03, respectively (P < 0.05). The average COM displacements for intensity-based DIR for the prostate, rectum, bladder, and seminal vesicles were 2.0 ± 1.5, 3.7 ± 1.4, 7.8 ± 2.2 and 3.6 ± 1.2 mm, whereas those values for hybrid DIR were 0.1 ± 0.0, 0.3 ± 0.2, 0.2 ± 0.1 and 0.6 ± 0.6 mm, respectively (P < 0.05). These results showed that the DSC for hybrid DIR had a higher DSC value and smaller COM displacement for all structures and all patients, compared with intensity-based DIR. Thus, the accumulative dose based on hybrid DIR might be trusted as a high-precision dose estimation method that takes into account organ movement during treatment radiotherapy.

Comparison of visual biofeedback system with a guiding waveform and abdomen-chest motion self-control system for respiratory motion management

Irregular breathing can influence the outcome of 4D computed tomography imaging and cause artifacts. Visual biofeedback systems associated with a patient-specific guiding waveform are known to reduce respiratory irregularities. In Japan, abdomen and chest motion self-control devices (Abches) (representing simpler visual coaching techniques without a guiding waveform) are used instead; however, no studies have compared these two systems to date. Here, we evaluate the effectiveness of respiratory coaching in reducing respiratory irregularities by comparing two respiratory management systems. We collected data from 11 healthy volunteers. Bar and wave models were used as visual biofeedback systems. Abches consisted of a respiratory indicator indicating the end of each expiration and inspiration motion. Respiratory variations were quantified as root mean squared error (RMSE) of displacement and period of breathing cycles. All coaching techniques improved respiratory variation, compared with free-breathing. Displacement RMSEs were 1.43 ± 0.84, 1.22 ± 1.13, 1.21 ± 0.86 and 0.98 ± 0.47 mm for free-breathing, Abches, bar model and wave model, respectively. Period RMSEs were 0.48 ± 0.42, 0.33 ± 0.31, 0.23 ± 0.18 and 0.17 ± 0.05 s for free-breathing, Abches, bar model and wave model, respectively. The average reduction in displacement and period RMSE compared with the wave model were 27% and 47%, respectively. For variation in both displacement and period, wave model was superior to the other techniques. Our results showed that visual biofeedback combined with a wave model could potentially provide clinical benefits in respiratory management, although all techniques were able to reduce respiratory irregularities.

TU-G-BRD-09: Evaluation of Patient DVH-Based QA Metrics for Prostate VMAT: Correlation Between Accuracy of Estimated 3D Patient Dose and MLC Position Error

Purpose: The purpose of this study was to evaluate the accuracy of patient DVH-based QA metrics and test the hypothesis: measurement-guided 3D dose reconstruction (MGDR) captures the induced dose error. Methods: We used 3DVH software with an ArcCHECK to estimate 3D patient dose by MGDR. Two different calculating modes of MGDR were used: “Normal Sensitivity” and “High Sensitivity”. Ten prostate cancer patients treated with hypo-fractionated VMAT (67.6 Gy/26 Fr) were studied. For the baseline plan, we induced MLC errors (−0.75, −0.5, −0.25, 0.25, 0.5 and 0.75 mm for each single bank), generated by in-house software. We compared the 3D patient dose estimated by 3DVH and that calculated by the treatment planning system. We evaluated the correlation between dose estimation error and MLC position error. Results: Slopes of linear fit to dose estimation error versus MLC position error for mean dose and D95 to the PTV were 1.76 and 1.40% mm-1, respectively, for “Normal Sensitivity” and −0.53 and −0.88% mm-1, respectively, for “High Sensitivity”, showing better accuracy for “High Sensitivity” than “Normal Sensitivity”. On the other hand, the slopes for mean dose to the rectum and bladder and V55 to the rectum and bladder were −1.00, −0.55, −3.53 and −1.85% mm-1, respectively, for “Normal Sensitivity” and −2.89, −2.39, −6.24 and −4.11% mm-1, respectively, for “High Sensitivity”, showing significant better accuracy for “Normal Sensitivity” than “High Sensitivity”. Conclusion: Our results showed that 3DVH had some residual error for both sensitivities, indicating the MGDR could capture a part of the induced error but not all of the induced error.Furthermore, we found that “Normal Sensitivity” might have better accuracy for the DVH metric for PTV and that “High Sensitivity” might have better accuracy for DVH metrics for the rectum and bladder. We must be willing to tolerate this residual error in clinical care.

A photon source model based on particle transport in a parameterized accelerator structure for Monte Carlo dose calculations

Purpose An accurate source model of a medical linear accelerator is essential for Monte Carlo (MC) dose calculations. This study aims to propose an analytical photon source model based on particle transport in parameterized accelerator structures, focusing on a more realistic determination of linac photon spectra compared to existing approaches. Methods We designed the primary and secondary photon sources based on the photons attenuated and scattered by a parameterized flattening filter. The primary photons were derived by attenuating bremsstrahlung photons based on the path length in the filter. Conversely, the secondary photons were derived from the decrement of the primary photons in the attenuation process. This design facilitates these sources to share the free parameters of the filter shape and be related to each other through the photon interaction in the filter. We introduced two other parameters of the primary photon source to describe the particle fluence in penumbral regions. All the parameters are optimized based on calculated dose curves in water using the pencil‐beam–based algorithm. To verify the modeling accuracy, we compared the proposed model with the phase space data (PSD) of the Varian TrueBeam 6 and 15 MV accelerators in terms of the beam characteristics and the dose distributions. The EGS5 Monte Carlo code was used to calculate the dose distributions associated with the optimized model and reference PSD in a homogeneous water phantom and a heterogeneous lung phantom. We calculated the percentage of points passing 1D and 2D gamma analysis with 1%/1 mm criteria for the dose curves and lateral dose distributions, respectively. Results The optimized model accurately reproduced the spectral curves of the reference PSD both on‐ and off‐axis. The depth dose and lateral dose profiles of the optimized model also showed good agreement with those of the reference PSD. The passing rates of the 1D gamma analysis with 1%/1 mm criteria between the model and PSD were 100% for 4 × 4, 10 × 10, and 20 × 20 cm2 fields at multiple depths. For the 2D dose distributions calculated in the heterogeneous lung phantom, the 2D gamma pass rate was 100% for 6 and 15 MV beams. The model optimization time was less than 4 min. Conclusion The proposed source model optimization process accurately produces photon fluence spectra from a linac using valid physical properties, without detailed knowledge of the geometry of the linac head, and with minimal optimization time.

Prognostic Value of Radiation Pneumonitis After Stereotactic Body Radiotherapy: Effect of Pulmonary Emphysema Quantitated Using CT Images

Background The aim of this study was to determine the prognostic factors of radiation pneumonitis (RP) after stereotactic body radiotherapy (SBRT). Patients and Methods A total of 50 patients (36 male and 14 female) were treated with SBRT for 42 primary lung cancers and 8 metastatic lung cancers. SBRT was performed with 48 Gy in 4 fractions to the isocenter or with 40 Gy in 4 fractions covering 95% of the planning target volume. Percentage of low attenuation area (%LAA) was defined as percentage of the lung area with attenuation of −860 Hounsfield units (HU) or lower (%LAA‐860) or of −960 HU or lower (%LAA‐960). The dosimetric parameter of V20 Gy, which means percentage volume of the lung receiving 20 Gy or more, was recalculated. RP was assessed using Common Terminology Criteria for Adverse Events version 4.0. Results The median follow‐up period was 39.0 months (range, 7.2‐94.5 months). RP of Grade 0, Grade 1, and Grade 2 to 3 was diagnosed in 11, 29, and 10 patients, respectively. Multivariate analyses (MVA) for Grade 1 showed that higher %LAA‐860 and higher %LAA‐960 were significantly associated with a lower rate of Grade 1 RP. MVA for Grade 2 to 3 showed that lower Brinkman index and lower lung V20 Gy were significantly associated with a lower rate of Grade 2 to 3 RP, and, in contrast, %LAA‐860 and %LAA‐960 had no association with Grade 2 to 3 RP. Conclusion This result suggests that high %LAA is associated with radiological changes (Grade 1) but that %LAA has no correlation with Grade 2 to 3 RP because symptomatic RP might also be affected by other factors. Micro‐Abstract The prognostic role of pulmonary emphysema for radiation pneumonitis (RP) after stereotactic body radiotherapy was investigated. It is true that patients with pulmonary emphysema showed a lower rate of abnormal shadow, but there was no association with Grade 2 to 3 RP. This is because patients with pulmonary emphysema had a low tolerance for symptomatic RP because of poor pulmonary function.

SU-F-I-24: Feasibility of Magnetic Susceptibility to Relative Electron Density Conversion Method for Radiation Therapy

PURPOSE The aim of this study is to develop radiation treatment planning using magnetic susceptibility obtained from quantitative susceptibility mapping (QSM) via MR imaging. This study demonstrates the feasibility of a method for generating a substitute for a CT image from an MRI. METHODS The head of a healthy volunteer was scanned using a CT scanner and a 3.0 T MRI scanner. The CT imaging was performed with a slice thickness of 2.5 mm at 80 and 120 kV (dual-energy scan). These CT images were converted to relative electron density (rED) using the CT-rED conversion table generated by a previous dual-energy CT scan. The CT-rED conversion table was generated using the conversion of the energy-subtracted CT number to rED via a single linear relationship. One T2 star-weighted 3D gradient echo-based sequence with four different echo times images was acquired using the MRI scanner. These T2 star-weighted images were used to estimate the phase data. To estimate the local field map, a Laplacian unwrapping of the phase and background field removal algorithm were implemented to process phase data. To generate a magnetic susceptibility map from the local field map, we used morphology enabled dipole inversion method. The rED map was resampled to the same resolution as magnetic susceptibility, and the magnetic susceptibility-rED conversion table was obtained via voxel-by-voxel mapping between the magnetic susceptibility and rED maps. RESULTS A correlation between magnetic susceptibility and rED is not observed through our method. CONCLUSION Our results show that the correlation between magnetic susceptibility and rED is not observed. As the next step, we assume that the voxel of the magnetic susceptibility map comprises two materials, such as water (0 ppm) and bone (-2.2 ppm) or water and marrow (0.81ppm). The elements of each voxel were estimated from the ratio of the two materials.

TU-AB-202-02: Deformable Image Registration Accuracy Between External Beam Radiotherapy and HDR Brachytherapy CT Images for Cervical Cancer Using a 3D-Printed Deformable Pelvis Phantom

PURPOSE Accurate deformable image registration (DIR) between external beam radiotherapy (EBRT) and HDR brachytherapy (BT) CT images in cervical cancer is challenging. DSC has been evaluated only on the basis of the consistency of the structure, and its use does not guarantee an anatomically reasonable deformation. We evaluate the DIR accuracy for cervical cancer with DSC and anatomical landmarks using a 3D-printed pelvis phantom. METHODS A 3D-printed, deformable female pelvis phantom was created on the basis of the patients CT image. Urethane and silicon were used as materials for creating the uterus and bladder, respectively, in the phantom. We performed DIR in two cases: case-A with a full bladder (170 ml) in both the EBRT and BT images and case-B with a full bladder in the BT image and a half bladder (100 ml) in the EBRT image. DIR was evaluated using DSCs and 70 uterus and bladder landmarks. A Hybrid intensity and structure DIR algorithm with two settings (RayStation) was used. RESULTS In the case-A, DSCs of the intensity-based DIR were 0.93 and 0.85 for the bladder and uterus, respectively, whereas those of hybrid-DIR were 0.98 and 0.96, respectively. The mean landmark error values of intensity-based DIR were 0.73±0.29 and 1.70±0.19 cm for the bladder and uterus, respectively, whereas those of Hybrid-DIR were 0.43±0.33 and 1.23±0.25 cm, respectively. In both cases, the Hybrid-DIR accuracy was better than the intensity-based DIR accuracy for both evaluation methods. However, for several bladder landmarks, the Hybrid-DIR landmark errors were larger than the corresponding intensity-based DIR errors (e.g., 2.26 vs 1.25 cm). CONCLUSION Our results demonstrate that Hybrid-DIR can perform with a better accuracy than the intensity-based DIR for both DSC and landmark errors; however, Hybrid-DIR shows a larger landmark error for some landmarks because the technique focuses on both the structure and intensity.