Technical Note: Range verification system using edge detection method for a scintillator and a CCD camera system.

PURPOSE: Three-dimensional irradiation with a scanned carbon-ion beam has been performed from 2011 at the authors' facility. The authors have developed the rotating-gantry equipped with the scanning irradiation system. The number of combinations of beam properties to measure for the commissioning is more than 7200, i.e., 201 energy steps, 3 intensities, and 12 gantry angles. To compress the commissioning time, quick and simple range verification system is required. In this work, the authors develop a quick range verification system using scintillator and charge-coupled device (CCD) camera and estimate the accuracy of the range verification. METHODS: A cylindrical plastic scintillator block and a CCD camera were installed on the black box. The optical spatial resolution of the system is 0.2 mm/pixel. The camera control system was connected and communicates with the measurement system that is part of the scanning system. The range was determined by image processing. Reference range for each energy beam was determined by a difference of Gaussian (DOG) method and the 80% of distal dose of the depth-dose distribution that were measured by a large parallel-plate ionization chamber. The authors compared a threshold method and a DOG method. RESULTS: The authors found that the edge detection method (i.e., the DOG method) is best for the range detection. The accuracy of range detection using this system is within 0.2 mm, and the reproducibility of the same energy measurement is within 0.1 mm without setup error. CONCLUSIONS: The results of this study demonstrate that the authors' range check system is capable of quick and easy range verification with sufficient accuracy.

[Evaluation of a Post-analysis Method for Cumulative Dose Distribution in Stereotactic Body Radiotherapy].

PURPOSE: The purpose of this study was to evaluate a post-analysis method for cumulative dose distribution in stereotactic body radiotherapy (SBRT) using volumetric modulated arc therapy (VMAT) . METHOD: VMAT is capable of acquiring respiratory signals derived from projection images and machine parameters based on machine logs during VMAT delivery. Dose distributions were reconstructed from the respiratory signals and machine parameters in the condition where respiratory signals were without division, divided into 4 and 10 phases. The dose distribution of each respiratory phase was calculated on the planned four-dimensional CT (4DCT). Summation of the dose distributions was carried out using deformable image registration (DIR), and cumulative dose distributions were compared with those of the corresponding plans. RESULTS AND DISCUSSION: Without division, dose differences between cumulative distribution and plan were not significant. In the condition where respiratory signals were divided, dose differences were observed over dose in cranial region and under dose in caudal region of planning target volume (PTV). Differences between 4 and 10 phases were not significant. CONCLUSION: The present method was feasible for evaluating cumulative dose distribution in VMAT-SBRT using 4DCT and DIR.

Experimental verification of gain drop due to general ion recombination for a carbon-ion pencil beam.

PURPOSE: Accurate dose measurement in radiotherapy is critically dependent on correction for gain drop, which is the difference of the measured current from the ideal saturation current due to general ion recombination. Although a correction method based on the Boag theory has been employed, the theory assumes that ionized charge density in an ionization chamber (IC) is spatially uniform throughout the irradiation volume. For particle pencil beam scanning, however, the charge density is not uniform, because the fluence distribution of a pencil beam is not uniform. The aim of this study was to verify the effect of the nonuniformity of ionized charge density on the gain drop due to general ion recombination. METHODS: The authors measured the saturation curve, namely, the applied voltage versus measured current, using a large plane-parallel IC and 24-channel parallel-plate IC with concentric electrodes. To verify the effect of the nonuniform ionized charge density on the measured saturation curve, the authors calculated the saturation curve using a method which takes into account the nonuniform ionized charge density and compared it with the measured saturation curves. RESULTS: Measurement values of the different saturation curves in the different channels of the concentric electrodes differed and were consistent with the calculated values. The saturation curves measured by the large plane-parallel IC were also consistent with the calculation results, including the estimation error of beam size and of setup misalignment. Although the impact of the nonuniform ionized charge density on the gain drop was clinically negligible with the conventional beam intensity, it was expected that the impact would increase with higher ionized charge density. CONCLUSIONS: For pencil beam scanning, the assumption of the conventional Boag theory is not valid. Furthermore, the nonuniform ionized charge density affects the prediction accuracy of gain drop when the ionized charge density is increased by a higher dose rate and/or lower beam size.

[Dose reconstruction using respiratory signals and machine parameters during treatment in stereotactic body radiotherapy].

PURPOSE: Volumetric modulated arc therapy (VMAT) is a rotational intensity-modulated radiotherapy (IMRT) technique capable of acquiring projection images during treatment. The purpose of this study was to reconstruct the dose distribution from respiratory signals and machine parameters acquired during stereotactic body radiotherapy (SBRT). METHODS: The treatment plans created for VMAT-SBRT included the constraint of 1 mm/degree in multileaf collimator (MLC) for a moving phantom and three patients with lung tumors. The respiratory signals were derived from projection images acquired during VMAT delivery, while the machine parameters were derived from machine logs. The respiratory signals and machine parameters were then linked along with the gantry angle. With this data, the dose distribution of each respiratory phase was calculated on the planned four-dimensional CT (4D CT). The doses at the isocenter, the point of max dose and the centroid of the target were compared with those of the corresponding plans. RESULTS AND DISCUSSION: In the phantom study, the maximum dose difference between the plan and "in-treatment" results was -0.4% at the centroid of the target. In the patient study, the difference was -1.8 ± 0.4% at the centroid of the target. Dose differences of the evaluated points between 4 and 10 phases were not significant. CONCLUSION: The present method successfully reconstructed the dose distribution using the respiratory signals and machine parameters acquired during treatment. This is a feasible method for verifying the actual dose for a moving target.

Stereotactic body radiotherapy for small lung tumors in the University of Tokyo Hospital.

Our work on stereotactic body radiation therapy (SBRT) for primary and metastatic lung tumors will be described. The eligibility criteria for SBRT, our previous SBRT method, the definition of target volume, heterogeneity correction, the position adjustment using four-dimensional cone-beam computed tomography (4D CBCT) immediately before SBRT, volumetric modulated arc therapy (VMAT) method for SBRT, verifying of tumor position within internal target volume (ITV) using in-treatment 4D-CBCT during VMAT-SBRT, shortening of treatment time using flattening-filter-free (FFF) techniques, delivery of 4D dose calculation for lung-VMAT patients using in-treatment CBCT and LINAC log data with agility multileaf collimator, and SBRT method for centrally located lung tumors in our institution will be shown. In our institution, these efforts have been made with the goal of raising the local control rate and decreasing adverse effects after SBRT.

[Winston-lutz test and acquisition of flexmap using rotational irradiation].

PURPOSE: IGRT (image guided radiation therapy) is a useful technique for implementing precisely targeted radiation therapy. Quality assurance and quality control (QA/QC) medical linear accelerators with a portal imaging system (electronic portal imaging device: EPID) are the key to ensuring safe IGRT. The Winston-Lutz test (WLT) provides an evaluation of the MV isocenter, which is the intersection of radiation, collimator, and couch isocenters. A flexmap can indicate a displacement of EPID from the beam center axis as a function of gantry angles which can be removed from the images. The purpose of this study was to establish a novel method for simultaneously carrying out WLT and acquiring a flexmap using rotational irradiation. We also observed long-term changes in flexmaps over a period of five months. METHOD: We employed rotational irradiation with a rectangular field (30×30 mm). First, the displacement of EPID from the beam center axis, indicated by the ball bearing (BB) center, was evaluated using an in-house program. The location of the BB center was then modified according to WLT. Second, a second irradiation was used to acquire a flexmap. We performed this examination regularly and evaluated long-term changes in the flexmap. RESULTS AND DISCUSSION: It proved feasible to perform WLT and flexmap measurements using our proposed methods. The precision of WLT using rotational irradiation was 0.1 mm. In flexmap analysis, the maximum displacement from the mean value for each angle was 0.4 mm over five months. CONCLUSION: We have successfully established a novel method of simultaneously carrying out WLT and flexmap acquisition using rotational irradiation. Maximum displacement from the mean in each angle was 0.4 mm over five months.

Dose verification of volumetric modulated arc therapy (VMAT) by use of in-treatment linac parameters.

Linac parameters such as the multi-leaf collimator (MLC) position and jaw position, cumulative monitor units (MUs), and the corresponding gantry angle were recorded during the clinical delivery of volumetric modulated arc therapy for prostate, lung, and head/neck cancer patients. Then, linac parameters were converted into the beam-data format used in the treatment planning system, and the dose distribution was reconstructed. The dose-volume histogram and the dose difference (DD) were compared with the corresponding values in the treatment plan. A reproducible error of in-treatment linac parameters was observed when a sudden change of beam intensity or MLC/jaw speed occurred. The maximum cumulative MU error was more than 4 MU for lung cancer cases, and the maximum MLC position exceeded 5 mm for prostate and head/neck cancer patients. However, these errors were quickly compensated for at the next control point. All treatments analyzed in the present study were delivered within 0.4% accuracy at the planning target volume. The cumulative dose agreed with that of the plan within 3% of the prescribed dose. The 1% DD was 93.9, 99.9, and 93.4% of the prescription dose for prostate, lung, and head/neck cancer patients, respectively.

4D registration and 4D verification of lung tumor position for stereotactic volumetric modulated arc therapy using respiratory-correlated cone-beam CT.

We propose a clinical workflow of stereotactic volumetric modulated arc therapy (VMAT) for a lung tumor from planning to tumor position verification using 4D planning computed tomography (CT) and 4D cone-beam CT (CBCT). A 4D CT scanner, an Anzai belt and a BodyFix were employed to obtain 10-phase respiratory-correlated CT data for a lung patient under constrained breathing conditions. A planning target volume (PTV) was defined by adding a 5-mm margin to an internal target volume created from 10 clinical target volumes, each of which was delineated on each of the 10-phase planning CT data. A single-arc VMAT plan was created with a D(95) prescription dose of 50 Gy in four fractions on the maximum exhalation phase CT images. The PTV contours were exported to a kilovoltage CBCT X-ray Volume Imaging (XVI) equipped with a linear accelerator (linac). Immediately before treatment, 10-phase 4D CBCT images were reconstructed leading to animated lung tumor imaging. Initial bone matching was performed between frame-averaged 4D planning CT and frame-averaged 4D CBCT datasets. Subsequently, the imported PTV contours and the animated moving tumor were simultaneously displayed on the XVI monitor, and a manual 4D registration was interactively performed on the monitor until the moving tumor was symmetrically positioned inside the PTV. A VMAT beam was delivered to the patient and during the delivery further 4D CBCT projection data were acquired to verify the tumor position. The entire process was repeated for each fraction. It was confirmed that the moving tumor was positioned inside the PTV during the VMAT delivery.

Four-dimensional measurement of the displacement of internal fiducial and skin markers during 320-multislice computed tomography scanning of breast cancer.

PURPOSE: To study the three-dimensional movement of internal tumor bed fiducial and breast skin markers, using 320-multislice computed tomography (CT); and to analyze intrafractional errors for breast cancer patients undergoing breast irradiation. METHODS AND MATERIALS: This study examined 280 markers on the skin of the breast (200 markers) and on the primary tumor bed (80 markers) of 20 patients treated by external-beam photon radiotherapy. Motion assessment was analyzed in 41 respiratory phases during 20 s of cine CT in the radiotherapy position. To assess intrafractional errors resulting from respiratory motion, four-dimensional CT scans were acquired for 20 patients. RESULTS: Motion in the anterior-posterior (A/P) and superior-inferior (S/I) directions showed a strong correlation (|r| > 0.7) with the respiratory curve for most markers (79% and 70%, respectively). The average marker displacements between maximum and minimum value during 20 s for the 200 breast skin metal markers were 1.1 ± 0.3 mm, 2.1 ± 0.6 mm, and 1.6 ± 0.4 mm in the left-right, A/P, and S/I directions, respectively. For the 80 tumor bed clips, displacements were 0.9 ± 0.2 mm in left-right, 1.7 ± 0.5 mm in A/P, and 1.1 ± 0.3 mm in S/I. There was no significant difference in the motion between breast quadrant regions or between the primary site and the other regions. CONCLUSIONS: Motion in primary breast tumors was evaluated with 320-multislice CT. Very little change was detected during individual radiation treatment fractions.

In-treatment 4D cone-beam CT with image-based respiratory phase recognition.

The use of respiration-correlated cone-beam computed tomography (4D-CBCT) appears to be crucial for implementing precise radiation therapy of lung cancer patients. The reconstruction of 4D-CBCT images requires a respiratory phase. In this paper, we propose a novel method based on an image-based phase recognition technique using normalized cross correlation (NCC). We constructed the respiratory phase by searching for a region in an adjacent projection that achieves the maximum correlation with a region in a reference projection along the cranio-caudal direction. The data on 12 lung cancer patients acquired just prior to treatment and on 3 lung cancer patients acquired during volumetric modulated arc therapy treatment were analyzed in the search for the effective area of cone-beam projection images for performing NCC with 12 combinations of registration area and segment size. The evaluation was done by a "recognition rate" defined as the ratio of the number of peak inhales detected with our method to that detected by eye (manual tracking). The average recognition rate of peak inhale with the most efficient area in the present method was 96.4%. The present method was feasible even when the diaphragm was outside the field of view. With the most efficient area, we reconstructed in-treatment 4D-CBCT by dividing the breathing signal into four phase bins; peak exhale, peak inhale, and two intermediate phases. With in-treatment 4D-CBCT images, it was possible to identify the tumor position and the tumor size in moments of inspiration and expiration, in contrast to in-treatment CBCT reconstructed with all projections.

Evaluation of heterogeneity dose distributions for Stereotactic Radiotherapy (SRT): comparison of commercially available Monte Carlo dose calculation with other algorithms.

BACKGROUND: The purpose of this study was to compare dose distributions from three different algorithms with the x-ray Voxel Monte Carlo (XVMC) calculations, in actual computed tomography (CT) scans for use in stereotactic radiotherapy (SRT) of small lung cancers. METHODS: Slow CT scan of 20 patients was performed and the internal target volume (ITV) was delineated on Pinnacle3. All plans were first calculated with a scatter homogeneous mode (SHM) which is compatible with Clarkson algorithm using Pinnacle3 treatment planning system (TPS). The planned dose was 48 Gy in 4 fractions. In a second step, the CT images, structures and beam data were exported to other treatment planning systems (TPSs). Collapsed cone convolution (CCC) from Pinnacle3, superposition (SP) from XiO, and XVMC from Monaco were used for recalculating. The dose distributions and the Dose Volume Histograms (DVHs) were compared with each other. RESULTS: The phantom test revealed that all algorithms could reproduce the measured data within 1% except for the SHM with inhomogeneous phantom. For the patient study, the SHM greatly overestimated the isocenter (IC) doses and the minimal dose received by 95% of the PTV (PTV95) compared to XVMC. The differences in mean doses were 2.96 Gy (6.17%) for IC and 5.02 Gy (11.18%) for PTV95. The DVH's and dose distributions with CCC and SP were in agreement with those obtained by XVMC. The average differences in IC doses between CCC and XVMC, and SP and XVMC were -1.14% (p = 0.17), and -2.67% (p = 0.0036), respectively. CONCLUSIONS: Our work clearly confirms that the actual practice of relying solely on a Clarkson algorithm may be inappropriate for SRT planning. Meanwhile, CCC and SP were close to XVMC simulations and actual dose distributions obtained in lung SRT.

4D-CBCT reconstruction using MV portal imaging during volumetric modulated arc therapy.

BACKGROUND: Recording target motion during treatment is important for verifying the irradiated region. Recently, cone-beam computed tomography (CBCT) reconstruction from portal images acquired during volumetric modulated arc therapy (VMAT), known as VMAT-CBCT, has been investigated. In this study, we developed a four-dimensional (4D) version of the VMAT-CBCT. MATERIALS AND METHODS: The MV portal images were sequentially acquired from an electronic portal imaging device. The flex, background, monitor unit, field size, and multi-leaf collimator masking corrections were considered during image reconstruction. A 4D VMAT-CBCT requires a respiratory signal during image acquisition. An image-based phase recognition (IBPR) method was performed using normalised cross correlation to extract a respiratory signal from the series of portal images. RESULTS: Our original IBPR method enabled us to reconstruct 4D VMAT-CBCT with no external devices. We confirmed that 4D VMAT-CBCT was feasible for two patients and in good agreement with in-treatment 4D kV-CBCT. CONCLUSION: The visibility of the anatomy in 4D VMAT-CBCT reconstruction for lung cancer patients has the potential of using 4D VMAT-CBCT as a tool for verifying relative positions of tumour for each respiratory phase.

Field-size effect of physical doses in carbon-ion scanning using range shifter plates.

A field-size effect of physical doses was studied in scanning irradiation with carbon ions. For the target volumes of 60 x 60 x 80, 40 x 40 x 80, and 20 x 20 x 80 mm3, the doses along the beam axis within the spread-out Bragg peaks reduced to 99.4%, 98.2%, and 96.0% of the dose for the target of 80 x 80 x 80 mm3, respectively. The present study revealed that the observed reductions can be compensated for by adopting the three-Gaussian form of lateral dose distributions for the pencil beam model used in the treatment planning system. The parameters describing the form were determined through the irradiation experiments making flat concentric squared frames with a scanned carbon beam. Since utilizing the three-Gaussian model in the dose optimization loop is at present time consuming, the correction for the field-size effect should be implemented as a "predicted-dose scaling factor." The validity of this correction method was confirmed through the irradiation of a target of 20 x 20 x 80 mm3.

Delivery verification using 3D dose reconstruction based on fluorescence measurement in a carbon beam scanning irradiation system.

The authors have developed a method to reconstruct the 3D dose distribution in a particle beam scanning irradiation system. In this scheme, 3D dose distribution is reconstructed by using the measured images of fluence distribution, which are taken for each iso-energy slice (i.e., the unit of the depth scanning). A fluorescent screen with a CCD camera is used to measure the fluence distribution. This system was installed at the HIMAC experimental port and tested by using carbon ion beams. Since a maximum difference between the reconstructed dose and the ionization chamber measurement was around 5% in the target volume, this system can be useful for quick verification of 3D dose distribution.