An assessment of the efficiency of methods for measurement of the computed tomography dose index (CTDI) for cone beam (CBCT) dosimetry by Monte Carlo simulation

Visualization and comparison of CT dose distribution between axial and helical acquisitions

BACKGROUND

It is challenging to assess the accuracy of volume CT Dose Index (CTDIvol ) when the axial scan modes corresponding to a helical scan protocol are not available. An alternative approach was proposed to directly measureC T D I v o l H $CTDI_{vol}^H$ using helical acquisitions and relatively small differences (< 20%) from CTDIvol were observed.

PURPOSE

To visually demonstrate the 3D dose distribution for both axial and helical CT acquisitions and quantitively compareC T D I v o l H $CTDI_{vol}^H$ and CTDIvol .

METHODS

3D dose distribution within the standard CTDI phantoms (16 and 32 cm diameter) from a single CT projection, Dp (x,y,z) was first generated using Monte Carlo simulation (GEANT4) with 9×108 photons per combination of tube voltage (80-140 kV), collimation width (1-8 cm), and z-axis location of the central ray of the x-ray beam, with a spatial resolution of 1 mm3 . These dose distributions from one single projection were analytically ensembled to simulate 3D dose volumes DA (x,y,z) and DH (x,y,z) for axial and helical scans, respectively, with different helical pitches (0.3-2) and scan lengths (100-150 mm). 2D planar dose distributions were obtained by integrating the inside 100 mm of the dose volumes. CTDIvol andC T D I v o l H $CTDI_{vol}^H\;$ were calculated using the planar dose data at corresponding pencil chamber locations and the percentage differences (PD) were reported.

RESULTS

High spatial resolution 3D CT dose volumes were generated and visualized. PDs betweenC T D I v o l H $CTDI_{vol}^H$ and CTDIvol had strong dependency on scan length and peripheral chamber locations, with subtle dependency on collimation width and pitch. PDs were mostly within the range of ± 3% for a scan length of 150 mm with four peripheral chamber locations.

CONCLUSIONS

With a scan length covering the entire phantom length,C T D I v o l H $CTDI_{vol}^H$ directly measured from helical scans can serve as an alternative to CTDIvol only if all four peripheral locations were measured.

Variations in size-specific effective dose with patient stature and beam width for kV cone beam CT imaging in radiotherapy

The facilities now available on linear accelerators for external beam radiotherapy enable radiation fields to be conformed to the shapes of tumours with a high level of precision. However, in order for the treatment delivered to take advantage of this, the patient must be positioned on the couch with the same degree of accuracy. Kilovoltage cone beam computed tomography systems are now incorporated into radiotherapy linear accelerators to allow imaging to be performed at the time of treatment, and image-guided radiation therapy is now standard in most radiotherapy departments throughout the world. However, because doses from imaging are much lower than therapy doses, less effort has been put into optimising radiological protection of imaging protocols. Standard imaging protocols supplied by the equipment vendor are often used with little adaptation to the stature of individual patients, and exposure factors and field sizes are frequently larger than necessary. In this study, the impact of using standard protocols for imaging anatomical phantoms of varying size from a library of 193 adult phantoms has been evaluated. Monte Carlo simulations were used to calculate doses for organs and tissues for each phantom, and results combined in terms of size-specific effective dose (SED). Values of SED from pelvic scans ranged from 11 mSv to 22 mSv for male phantoms and 8 mSv to 18 mSv for female phantoms, and for chest scans from 3.8 mSv to 7.6 mSv for male phantoms and 4.6 mSv to 9.5 mSv for female phantoms. Analysis of the results showed that if the same exposure parameters and field sizes are used, a person who is 5 cm shorter will receive a size SED that is 3%-10% greater, while a person who is 10 kg lighter will receive a dose that is 10%-14% greater compared with the average size.

A novel analytical method for computing dose from kilovoltage beams used in Image-Guided radiation therapy

PURPOSE

A modified convolution/superposition algorithm is proposed to compute dose from the kilovoltage beams used in IGRT. The algorithm uses material-specific energy deposition kernels instead of water-energy deposition kernels.

METHODS

Monte Carlo simulation was used to model the Elekta XVI unit and determine dose deposition characteristics of its kilovoltage beams. The dosimetric results were compared with ion chamber measurements. The dose from the kilovoltage beams was then computed using convolution/superposition along with material-specific energy deposition kernels and compared with Monte Carlo and measurements. The material-specific energy deposition kernels were previously generated using Monte Carlo.

RESULTS

The obtained gamma indices (using 2%/2mm criteria for 95% of points) were lower than 1 in almost all instances which indicates good agreement between simulated and measured depth doses and profiles. The comparisons of the algorithm with measurements in a homogeneous solid water slab phantom, and that with Monte Carlo in a head and neck CT dataset produced acceptable results. The calculated point doses were within 4.2% of measurements in the homogeneous phantom. Gamma analysis of the calculated vs. Monte Carlo simulations in the head and neck phantom resulted in 94% of points passing with a 2%/2mm criteria.

CONCLUSIONS

The proposed method offers sufficient accuracy in kilovoltage beams dose calculations and has the potential to supplement the conventional megavoltage convolution/superposition algorithms for dose calculations in low energy range.

A new cone-beam computed tomography dosimetry method providing optimal measurement positions: A Monte Carlo study

PURPOSE

The conventional weighted computed tomography dose index (CTDI_w) may not be suitable for cone-beam computed tomography (CBCT) dosimetry because a cross-sectional dose distribution is angularly inhomogeneous owing to partial angle irradiations. This study was conducted to develop a new dose metric (f(0)^CB_w) for CBCT dosimetry to determine a more accurate average dose in the central cross-sectional plane of a cylindrical phantom using Monte Carlo simulations.

METHODS

First, cross-sectional dose distributions of cylindrical polymethyl methacrylate phantoms over a wide range of phantom diameters (8-40 cm) were calculated for various CBCT scan protocols. Then, by obtaining linear least-squares fits of the full datasets of the cross-sectional dose distributions, the optimal radial positions, which represented measurement positions for the average phantom dose, were determined. Finally, the f(0)^CB_w method was developed by averaging point doses at the optimal radial positions of the phantoms. To demonstrate its validity, the relative differences between the average doses and each dose index value were estimated for the devised f(0)^CB_w, conventional CTDI_w, and Haba's CTDI_w methods, respectively.

RESULTS

The relative differences between the average doses and each dose index value were within 4.1%, 16.7%, and 11.9% for the devised, conventional CTDI_w, and Haba's CTDI_w methods, respectively.

CONCLUSIONS

The devised f(0)^CB_w value was calculated by averaging four "point doses" at 90° intervals and the optimal radial positions of the cylindrical phantom. The devised method can estimate the average dose more accurately than the previously developed CTDI_w methods for CBCT dosimetry.

Accuracy of weighted CTDI in estimating average dose delivered to CTDI phantoms: An experimental study

PURPOSE

The concept of the weighted computed tomography dose index ( CTDI w ) was proposed in 1995 to represent the average CTDI across an axial section of a cylindrical phantom. The purpose of this work was to experimentally re-examine the validity of the underlying assumptions behind CTDI w for modern MDCT systems.

METHODS

To enable experimental mapping of CTDI 100 in the axial plane, in-house 16 and 32 cm cylindrical phantoms were fabricated to allow the pencil chamber to reach any arbitrary axial location within the phantoms. The phantoms were scanned on a clinical MDCT with five beam collimation widths, three bowtie filters, and four kV levels. To evaluate the linearity and rotational invariance assumptions implicitly made when the weighting factors of 1/3 and 2/3 in the CTDI w formula were originally derived, CTDI 100 was measured at different radial and angular locations within the phantom for different collimation, bowtie, and kV combinations. The average CTDI ( CTDI avg ) across the axial plane was calculated from the experimental two-dimensional (2D) dose distribution and was compared with the traditional CTDI w .

RESULTS

For both phantoms under all scan conditions, the axial dose distributions were found to have significant angular dependence, potentially due to the x-ray attenuation by the patient couch or the head holder. The radial dose profiles were also found to significantly deviate from linearity in many cases due to the presence of the bowtie filter. When only the 12 o'clock peripheral CTDI 100 and the traditional weighting factors were used to calculate CTDI w , the average dose was overestimated in the 16 cm phantom by up to 8.4% at isocenter and up to 35.3% when the phantom was off-centered by 6 cm; in the 32 cm phantom at isocenter, the average dose was overestimated by up to 12.8%. Using an average of the four peripheral CTDI 100 measurements at the 12, 3, 6, and 9 o'clock locations reduced the error of CTDI w to within 1.2% in the 16 cm phantom. For the 32 cm phantom, even by using the average of the peripheral measurements, the traditional CTDI w underestimated the average dose by up to 4.3% due to aggressive drop-off of the CTDI 100 at the phantom periphery.

CONCLUSIONS

The linearity and rotational-invariance assumptions behind the traditional CTDI w formalism may not be valid for modern CT systems and thus CTDI w may not accurately represent the average dose or radiation output within a CTDI phantom. Utilizing data from all four peripheral locations always improves accuracy of CTDI w in representing the true average dose. For the large (32 cm) phantom, nonlinear models and more measurement points are needed if a more precise estimation of the average axial dose is required.

A Monte Carlo investigation of dose length product of cone beam computed tomography scans

The dose length product (DLP) provides a measurement related to energy imparted from a computed tomography (CT) scan. The DLP is based on the volume-averaged CT dose index (CTDI _vol), which is designed for fan beams. The aims of this study were to investigate the use of DLP for scans with wide beams used in cone beam CT (DLP _CBCT) in radiotherapy that would be analogous to the DLP of fan beam scans (DLP _CT), and to compare the efficiencies of DLP _CT and DLP _CBCT in reporting the total energy imparted in patients. A validated Monte Carlo model of a kV imaging system integrated into a Varian TrueBeam linac was employed. The DLP _CT was assessed by multiplying the CTDI _vol for a 20 mm fan beam by scan length, and the DLP _CBCT determined through multiplying the CTDI _vol, estimated for wide beams using a correction factor based on free-in-air measurements, by the beam width. Two scan protocols for head and body were investigated for tube potentials between 80 and 140 kV and a range of scan lengths/widths. Efficiency values were estimated by normalising the DLP _CT and DLP _CBCT with respect to the corresponding dose profile integrals (DPIs), which were evaluated within 900 mm long phantoms. The results show that the DLP _CBCT values were within 1% of those for DLP _CT of similar length performed on the same system, and the efficiencies decrease with tube potential. However, whereas DLP values for fan beams are approximately proportional to scan length, those for wide beams decrease by ∼2% between beam widths of 20 and 320 mm. As a result, while the DLP _CT efficiency is similar over all scan lengths, that for DLP _CBCT increases slightly with beam width. The DLP _CT and DLP _CBCT underestimated the total energy imparted by comparable amounts with efficiencies within the range of 80-81% and 80-83% for the head scans, and 71-76% and 70-77% for the body scans, respectively. The results indicate that the DLP _CBCT can be considered as an analogous dose index to the DLP _CT. It could, therefore, be used for quantification of doses from imaging in radiotherapy and provide a valuable tool to aid optimisation.

Dosimetric assessment of the exposure of radiotherapy patients due to cone-beam CT procedures

Cone-beam computed tomography (CBCT) is widely used for pre-treatment verification and patient setup in image-guided radiation therapy (IGRT). CBCT imaging is employed daily and several times per patient, resulting in potentially high cumulative imaging doses to healthy tissues that surround exposed target organs. Computed tomography dose index (CTDI) is the parameter used by CBCT equipment as indication of the radiation output to patients. This study aimed to increase the knowledge on the relation between CBCT organ doses and weighted CTDI (CTDI_W) for a thorax scanning protocol. A CBCT system was modelled using the Monte Carlo (MC) radiation transport program MCNPX2.7.0. Simulation results were validated against half-value layer (HVL), axial beam profile, patient skin dose (PSD) and CTDI measurements. For organ dose calculations, a male voxel phantom ("Golem") was implemented with the CBCT scanner computational model. After a successful MC model validation with measurements, a systematic comparison was performed between organ doses (and their distribution) and CTDI dosimetry concepts [CTDI_W and cumulative dose quantities f₁₀₀(150) and [Formula: see text]]. The results obtained show that CBCT organ doses vary between 1.2 ± 0.1 mGy and 3.3 ± 0.2 mGy for organs located within the primary beam. It was also verified that CTDI_W allows prediction of absorbed doses to tissues at distances of about 5 cm from the isocentre of the CBCT system, whereas f₁₀₀(150) allows prediction of organ doses at distances of about 10 cm from the isocentre, independently from its location. This study demonstrates that these dosimetric concepts are suitable methods that easily allow a good approximation of the additional CBCT imaging doses during a typical lung cancer IGRT treatment.

INVESTIGATION OF A PRACTICAL PATIENT DOSE INDEX FOR ASSESSMENT OF PATIENT ORGAN DOSE FROM CONE-BEAM COMPUTED TOMOGRAPHY IN RADIATION THERAPY USING A MONTE CARLO SIMULATION

The purpose of this study was to investigate a practical patient dose index for assessing the patient organ dose from a cone-beam computed tomography (CBCT) scan by comparing eight dose indices, i.e. CTDI100, CTDIIEC, CTDI∞, midpoint doses f(0)PMMA for a cylindrical polymethyl methacrylate (PMMA) phantom, f(0)Ap for an anthropomorphic phantom and f(0)Pat for a prostate cancer patient, as well as the conventional size specific dose estimations (SSDEconv) and modified SSDE (SSDEmod), with organ dose for the prostate (ODprost) obtained via Monte Carlo (MC) simulation. The ODprost was the reference dose used to find the practical dose index at the center of the pelvic region of a prostate cancer patient. The smallest error rate with respect to the ODprost of 19.3 mGy (reference) among eight dose indices was 5% for f(0)Pat. The practical patient dose index was the f(0)Pat, which showed the smallest error with respect to the reference dose.

On the improvement of CBCT image quality for soft tissue-based SRS localization

Purpose

We explore the optimal cone‐beam CT (CBCT) acquisition parameters to improve CBCT image quality to enhance intracranial stereotactic radiosurgery (SRS) localization and also assess the imaging dose levels associated with each imaging protocol.

Methods

Twenty‐six CBCT acquisition protocols were generated on an Edge^® linear accelerator (Varian Medical Systems, Palo Alto, CA) with different x‐ray tube current and potential settings, gantry rotation trajectories, and gantry rotation speeds. To assess image quality, images of the Catphan 504 phantom were analyzed to evaluate the following image quality metrics: uniformity, HU constancy, spatial resolution, low contrast detection, noise level, and contrast‐to‐noise ratio (CNR). To evaluate the imaging dose for each protocol, the cone‐beam dose index (CBDI) was measured. To validate the phantom results, further analysis was performed with an anthropomorphic head phantom as well as image data acquired for a clinical SRS patient.

Results

The Catphan data indicates that adjusting acquisition parameters had direct effects on the image noise level, low contrast detection, and CNR, but had minimal effects on uniformity, HU constancy, and spatial resolution. The noise level was reduced from 34.5 ± 0.3 to 18.5 ± 0.2 HU with a four‐fold reduction in gantry speed, and to 18.7 ± 0.2 HU with a four‐fold increase in tube current. Overall, the noise level was found to be proportional to inverse square root of imaging dose, and imaging dose was proportional to the product of total tube current‐time product and the cube of the x‐ray potential. Analysis of the anthropomorphic head phantom data and clinical SRS imaging data also indicates that noise is reduced with imaging dose increase.

Conclusions

Our results indicate that optimization of the imaging protocol, and thereby an increase in the imaging dose, is warranted for improved soft‐tissue visualization for intracranial SRS.

Initial evaluation of image performance of a 3-D x-ray system: phantom-based comparison of 3-D tomography with conventional computed tomography

Phantom-based initial performance assessment of a prototype three-dimensional (3-D) x-ray system and comparison of 3-D tomography with computed tomography (CT) were proposed. A 3-D image quality phantom was scanned with a prototype version of 3-D cone-beam CT imaging implemented on a twin robotic x-ray system using three trajectories (163 deg = table, 188 deg = upright, and 200 deg = side), six tube voltages (60, 70, 81, 90, 100, and 121 kV), and four detector doses (0.348, 0.696, 1.740, and [Formula: see text]). CT was obtained with a clinical protocol. Spatial resolution (line pairs/cm) and soft-tissue-contrast resolution were assessed by two independent readers. Radiation dose was assessed. Descriptive and analysis of variance (ANOVA) ([Formula: see text]) were performed. With 3-D tomography, a maximum of 16 lp/cm was visible and best soft-tissue-contrast resolution was 2 mm at 30 Hounsfield units (HU) for 160 projections. With CT, 10 lp/cm was visible and soft-tissue-contrast resolution was 4 mm at 20 HU. The upright trajectory yielded significantly better spatial resolution and soft tissue contrast, and the side trajectory yielded significantly higher soft tissue contrast than the table trajectory ([Formula: see text]). Radiation dose was higher in 3-D tomography (45 to 704 mGycm) than CT (44 mGycm). Three-dimensional tomography renders overall equal or higher spatial resolution and comparable soft tissue contrast to CT for medium- and high-dose protocols in the side and upright trajectories, but with higher radiation doses.

Evaluation of coefficients to derive organ and effective doses from cone-beam CT (CBCT) scans: a Monte Carlo study

Regular imaging is used throughout image guided radiation therapy to improve treatment delivery. In order for treatment procedures to be optimized, the doses delivered by imaging exposures should be taken into account. CT dosimetry methods based on the CT dose index (CTDI), measured with a 100 mm long pencil ionization chamber (CTDI₁₀₀) in standard phantoms, are not designed for cone-beam CT (CBCT) imaging systems used in radiotherapy, therefore a modified version has been proposed for CBCT by the International Electrotechnical Commission (CTDI_IEC). Monte Carlo simulations based on a Varian On-Board Imaging system were used to derive conversion coefficients that enable organ doses for ICRP reference phantoms to be determined from the CTDI_IEC for different scan protocols and different beam widths (80-320) mm. A dose-width product calculated by multiplying the CTDI_IEC by the width of the CBCT beam is proposed as a quantity that can be used for estimating effective dose. The variation in coefficients with CBCT beam width was studied. Coefficients to allow estimation of effective doses were derived, namely 0.0034 mSv (mGy cm)^-1 for the head, 0.0252 mSv (mGy cm)^-1 for the thorax, 0.0216 mSv (mGy cm)^-1 for the abdomen and 0.0150 mSv (mGy cm)^-1 for the pelvis, and these may be applicable more generally to other CBCT systems in radiotherapy. If data on effective doses are available, these can be used in making judgements on the contributions to patient dose from imaging, and thereby assist in optimization of the treatment regimes. The coefficients can also be employed in converting dosimetry data recorded in patient records into quantities relating directly to patient doses.

A Monte Carlo study of organ and effective doses of cone beam computed tomography (CBCT) scans in radiotherapy

Cone-beam CT (CBCT) scans utilised for image guided radiation therapy (IGRT) procedures have become an essential part of radiotherapy. The aim of this study was to assess organ and effective doses resulting from new CBCT scan protocols (head, thorax, and pelvis) released with a software upgrade of the kV on-board-imager (OBI) system. Organ and effective doses for protocols of the new software (V2.5) and a previous version (V1.6) were assessed using Monte Carlo (MC) simulations for the International Commission on Radiological Protection (ICRP) adult male and female reference computational phantoms. The number of projections and the mAs values were increased and the size of the scan field was extended in the new protocols. Influence of these changes on organ and effective doses of the scans was investigated. The OBI system was modelled in EGSnrc/BEAMnrc, and organ doses were estimated using EGSnrc/DOSXYZnrc. The MC model was benchmarked against experimental measurements. Organ doses resulting from the V2.5 protocols were higher than those of V1.6 for organs that were partially or fully inside the scans fields, and increased by (3-13)%, (10-77)%, and (13-21)% for the head, thorax, and pelvis protocols for both phantoms, respectively. As a result, effective doses rose by 14%, 17%, and 16% for the male phantom, and 13%, 18%, and 17% for the female phantom for the three scan protocols, respectively. The scan field extension for the V2.5 protocols contributed significantly in the dose increases, especially for organs that were partially irradiated such as the thyroid in head and thorax scans and colon in the pelvic scan. The contribution of the mAs values and projection numbers was minimal in the dose increases, up to 2.5%. The field size extension plays a major role in improving the treatment output by including more markers in the field of view to match between CBCT and CT images and hence setting up the patient precisely. Therefore, a trade-off between the risk and benefits of CBCT scans should be considered, and the dose increases should be monitored. Several recommendations have been made for optimisation of the patient dose involved for IGRT procedures.

Whole-body dose and energy measurements in radiotherapy by a combination of LiF:Mg,Cu,P and LiF:Mg,Ti

PURPOSE

Long-term survivors of cancer who were treated with radiotherapy are at risk of a radiation-induced tumor. Hence, it is important to model the out-of-field dose resulting from a cancer treatment. These models have to be verified with measurements, due to the small size, the high sensitivity to ionizing radiation and the tissue-equivalent composition, LiF thermoluminescence dosimeters (TLD) are well-suited for out-of-field dose measurements. However, the photon energy variation of the stray dose leads to systematic dose errors caused by the variation in response with radiation energy of the TLDs. We present a dosimeter which automatically corrects for the energy variation of the measured photons by combining LiF:Mg,Ti (TLD100) and LiF:Mg,Cu,P (TLD100H) chips.

METHODS

The response with radiation energy of TLD100 and TLD100H compared to ⁶⁰Co was taken from the literature. For the measurement, a TLD100H was placed on top of a TLD100 chip. The dose ratio between the TLD100 and TLD100H, combined with the ratio of the response curves was used to determine the mean energy. With the energy, the individual correction factors for TLD100 and TLD100H could be found. The accuracy in determining the in- and out-of-field dose for a nominal beam energy of 6MV using the double-TLD unit was evaluated by an end-to-end measurement. Furthermore, published Monte Carlo (M.C.) simulations of the mean photon energy for brachytherapy sources, stray radiation of a treatment machine and cone beam CT (CBCT) were compared to the measured mean energies. Finally, the photon energy distribution in an Alderson phantom was measured for different treatment techniques applied with a linear accelerator. Additionally, a treatment plan was measured with a cobalt machine combined with an MRI.

RESULTS

For external radiotherapy, the presented double-TLD unit showed a relative type A uncertainty in doses of -1%±2% at the two standard deviation level compared to an ionization chamber. The type A uncertainty in dose was in agreement with the theoretically calculated type B uncertainty. The measured energies for brachytherapy sources, stray radiation of a treatment machine and CBCT imaging were in agreement with M.C. simulations. A shift in energy with increasing distance to the isocenter was noticed for the various treatment plans measured with the Alderson phantom. The calculated type B uncertainties in energy were in line with the experimentally evaluated type A uncertainties.

CONCLUSION

The double-TLD unit is able to predict the photon energy of scatter radiation in external radiotherapy, X-ray imagine and brachytherapy sources. For external radiotherapy, the individual energy correction factors enabled a more accurate dose determination compared to conventional TLD measurements.

Technical note: No increase in effective dose from half compared to full rotation pelvis cone beam CT

PURPOSE

To image the abdomen of a patient with a gantry mounted imaging system of a linear accelerator, different cone beam computed tomography (CBCT) protocols are available. The whole-body dose of a full rotation abdomen CBCT and a half rotation CBCT was compared. In our clinic, both CBCT protocols are used in daily routine work.

METHODS

With an adult anthropomorphic Alderson phantom, the whole-body dose per CBCT scan was measured with thermoluminescence dosimeters. The half rotation CBCT was applied such that the gantry mounted X-ray source rotated around the right side of the phantom. The 183 measurement locations covered all ICRP recommended critical organs (except the gonads). The effective dose was calculated with the mean organ dose and the corresponding tissue weighting factors. A point-by-point dose comparison of both protocols was conducted.

RESULTS

The effective dose was 5.4 mSv ±5% and 5.0 mSv ±5% (estimated type B 1σ) for the full and the half rotation CBCT respectively. There was no significant difference (α = 0.05) in the effective dose within the precision of the measurement (1σ = 5%). The half rotation CBCT displayed an inhomogeneous dose distribution in a transversal phantom slice in contrast with the full rotation CBCT. In the imaging region, the mean dose was (20.5 ± 3.4) mGy and (19.2 ± 7.4) mGy (measured type A 1σ) for the full and the half rotation CBCT respectively.

CONCLUSION

The half compared to the full rotation CBCT displays a smaller field-of-view in a transversal slice and no significant difference in the effective dose. Hence, the full rotation CBCT is favorable compared to the half rotation CBCT. However, by using the half rotation protocol, critical volumes in the patient can be spared compared to the full rotation protocol.

Technical Note: Comparison of peripheral patient dose from MR-guided 60 Co therapy and 6 MV linear accelerator IGRT

PURPOSE

The use of X-ray imaging in radiation therapy can give a substantial dose to the patient. A Cobalt machine combined with an magnetic resonance imaging (MRI) was recently introduced to clinical work. One positive aspect of using non-ionizing imaging devices is the reduction of the patient exposure. The purpose of this study was to quantify the difference in out-of-field dose to the patient between the image guided radiation therapy (IGRT) treatment applied with a linear accelerator with cone beam CT (CBCT) equipment and a Cobalt machine combined with an MRI.

METHODS

The treatment of a rhabdomyosarcoma in the prostate was planned and irradiated using different modalities and radiation therapy machines. The whole-body dose was measured for a 3D-conformal radiation therapy (3DCRT), an intensity-modulated radiation therapy (IMRT), and a volumetric-modulated arc therapy plan applied with a conventional linear accelerator operated at 6 MV beam energy. Additionally, the dose of an IMRT plan applied with a ⁶⁰ Co machine combined with an MRI was measured. Furthermore, the dose of one CBCT scan using the linear accelerator's on-board imaging system was determined. The 3D dose measurements were performed in an anthropomorphic phantom containing 168 slots for thermoluminescence dosimeters (TLDs). A combination of LiF:Mg,Ti (TLD100) and LiF:Mg,Cu,P (TLD100H) was used to accurately determine the in- and out-of-field dose. The plans were rescaled to different fractionation schemes (2 Gy, 3 Gy, and 5 Gy fraction dose) and the dose of one CBCT scan was additionally added to the treatment dose per fraction applied with the linear accelerator. The resulting absorbed doses per fraction of the two machines were compared.

RESULTS

In the target region, all measured treatment plans presented the same magnitude of dose, while the CBCT dose was a factor of 100 smaller. Close to the planned target volume (PTV), the dose from the ⁶⁰ Co machine was a factor of two higher compared with the 3DCRT + CBCT dose. Up to 45 cm from the PTV, the treatment applied with the ⁶⁰ Co-sources showed an increased out-of-field dose compared to the linear accelerator + CBCT IGRT treatments. Further away from the PTV in the region where leakage from the gantry head is dominating, the out-of-field dose of the Cobalt machine was smaller compared to the linear accelerator + CBCT.

CONCLUSION

The peripheral dose of the ⁶⁰ Co machine combined with an MRI is larger up to 45 cm from the PTV and further away, it is lower than the dose from a linear accelerator + CBCT treatment. The presented fractionation schemes had a marginal impact on the results.

EGS_cbct: Simulation of a fan beam CT and RMI phantom for measured HU verification

INTRODUCTION

A mathematical 3D model of an existing computed tomography (CT) scanner was created and used in the EGSnrc-based BEAMnrc and egs_cbct Monte Carlo codes. Simulated transmission dose profiles of a RMI-465 phantom were analysed to verify Hounsfield numbers against measured data obtained from the CT scanner.

METHODS AND MATERIALS

The modelled CT unit is based on the design of a Toshiba Aquilion 16 LB CT scanner. As a first step, BEAMnrc simulated the X-ray tube, filters, and secondary collimation to obtain phase space data of the X-ray beam. A bowtie filter was included to create a more uniform beam intensity and to remove the beam hardening effects. In a second step the Interactive Data Language (IDL) code was used to build an EGSPHANT file that contained the RMI phantom which was used in egs_cbct simulations. After simulation a series of profiles were sampled from the detector model and the Feldkamp-Davis-Kress (FDK) algorithm was used to reconstruct transversal images. The results were tested against measured data obtained from CT scans.

RESULTS

The egs_cbct code can be used for the simulation of a fan beam CT unit. The calculated bowtie filter ensured a uniform flux on the detectors. Good correlation between measured and simulated CT numbers was obtained.

CONCLUSIONS

In principle, Monte Carlo codes such as egs_cbct can model a fan beam CT unit. After reconstruction, the images contained Hounsfield values comparable to measured data.

Organ doses can be estimated from the computed tomography (CT) dose index for cone-beam CT on radiotherapy equipment

Cone beam computed tomography (CBCT) systems are fitted to radiotherapy linear accelerators and used for patient positioning prior to treatment by image guided radiotherapy (IGRT). Radiotherapists' and radiographers' knowledge of doses to organs from CBCT imaging is limited. The weighted CT dose index for a reference beam of width 20 mm (CTDIw,ref) is displayed on Varian CBCT imaging equipment known as an On-Board Imager (OBI) linked to the Truebeam linear accelerator. This has the potential to provide an indication of organ doses. This knowledge would be helpful for guidance of radiotherapy clinicians preparing treatments. Monte Carlo simulations of imaging protocols for head, thorax and pelvic scans have been performed using EGSnrc/BEAMnrc, EGSnrc/DOSXYZnrc, and ICRP reference computational male and female phantoms to derive the mean absorbed doses to organs and tissues, which have been compared with values for the CTDIw,ref displayed on the CBCT scanner console. Substantial variations in dose were observed between male and female phantoms. Nevertheless, the CTDIw,ref gave doses within  ±21% for the stomach and liver in thorax scans and 2  ×  CTDIw,ref can be used as a measure of doses to breast, lung and oesophagus. The CTDIw,ref could provide indications of doses to the brain for head scans, and the colon for pelvic scans. It is proposed that knowledge of the link between CTDIw for CBCT should be promoted and included in the training of radiotherapy staff.

Evaluation of MVCT imaging dose levels during helical IGRT: comparison between ion chamber, TLD, and EBT3 films

The purpose of this investigation was to evaluate the dose on megavoltage CT (MVCT) images required for tomotherapy. As imaging possibilities are often used before each treatment and usually used several times before the session, we tried to evaluate the dose delivered during the procedure. For each scanning mode (fine, normal, and coarse), we first established the relative variation of these doses according to different technical parameters (explored length, patient setup). These dose variations measured with the TomoPhant, also known as Cheese phantom, showed the expected variations (due to the variation of scattered radiation) of 15% according to the explored length and ± 5% according to the phantom setup (due to the variation of the point of measurement in the bore). In order to estimate patient doses, an anthropomorphic phantom was used for thermoluminescent and film dosimetry. The degree of agreement between the two methods was very satisfactory (the differences correspond to 5 mGy per imaging session) for the three sites studied (head & neck, thorax, and abdomen). These measurements allowed us to estimate the delivered dose of between 1 cGy and 4 cGy according to the site and imaging mode. Finally, we attempted to investigate a way to calculate this delivered dose in our patients from the study conducted on a cylindrical phantom and by taking into account data from the initial kV-CT scan. The results we obtained were close to our measurements, with discrepancies below 5 mGy per MVCT.

Evaluation of cumulative dose for cone-beam computed tomography (CBCT) scans within phantoms made from different compositions using Monte Carlo simulations

Measurement of cumulative dose ƒ(0,150) with a small ionization chamber within standard polymethyl methacrylate (PMMA) CT head and body phantoms, 150 mm in length, is a possible practical method for cone-beam computed tomography (CBCT) dosimetry. This differs from evaluating cumulative dose under scatter equilibrium conditions within an infinitely long phantom ƒ(0,∞), which is proposed by AAPM TG-111 for CBCT dosimetry. The aim of this study was to investigate the feasibility of using ƒ(0,150) to estimate values for ƒ(0,∞) in long head and body phantoms made of PMMA, polyethylene (PE), and water, using beam qualities for tube potentials of 80-140 kV. The study also investigated the possibility of using 150 mm PE phantoms for assessment of ƒ(0,∞) within long PE phantoms, the ICRU/AAPM phantom. The influence of scan parameters, composition, and length of the phantoms was investigated. The capability of ƒ(0,150) to assess ƒ(0,∞) has been defined as the efficiency and assessed in terms of the ratios ε(ƒ(0,150) / ƒ(0,∞)). The efficiencies were calculated using Monte Carlo simulations for an On-Board Imager (OBI) system mounted on a TrueBeam linear accelerator. Head and body scanning protocols with beams of width 40-500 mm were used. Efficiencies ε(PMMA/PMMA) and ε(PE/PE) as a function of beam width exhibited three separate regions. For beam widths < 150 mm, ε(PMMA/PMMA) and ε(PE/PE) values were greater than 90% for the head and body phantoms. The efficiency values then fell rapidly with increasing beam width before levelling off at 74% for ε(PMMA/PMMA) and 69% for ε(PE/PE) for a 500 mm beam width. The quantities ε(PMMA/PE) and ε(PMMA/Water) varied with beam width in a different manner. Values at the centers of the phantoms for narrow beams were lower and increased to a steady state for ~100-150 mm wide beams, before declining with increasing the beam width, whereas values at the peripheries decreased steadily with beam width. Results for ε(PMMA/PMMA) were virtually independent of tube potential, but there was more variation for ε(PMMA/PE) and ε(PMMA/Water). ƒ(0,150) underestimated ƒ(0,∞) for beam widths used for CBCT scans, thus it is necessary to use long phantoms, or apply conversion factors (Cƒs) to measurements with standard PMMA CT phantoms. The efficiency values have been used to derive (Cƒs) to allow evaluation of ƒ(0,∞) from measurements of ƒ(0,150). The (Cƒs) only showed a weak dependence on scan parameters and scanner type, and so may be suitable for general application.

Investigation of practical approaches to evaluating cumulative dose for cone beam computed tomography (CBCT) from standard CT dosimetry measurements: a Monte Carlo study

A function called Gx(L) was introduced by the International Commission on Radiation Units and Measurements (ICRU) Report-87 to facilitate measurement of cumulative dose for CT scans within long phantoms as recommended by the American Association of Physicists in Medicine (AAPM) TG-111. The Gx(L) function is equal to the ratio of the cumulative dose at the middle of a CT scan to the volume weighted CTDI (CTDIvol), and was investigated for conventional multi-slice CT scanners operating with a moving table. As the stationary table mode, which is the basis for cone beam CT (CBCT) scans, differs from that used for conventional CT scans, the aim of this study was to investigate the extension of the Gx(L) function to CBCT scans. An On-Board Imager (OBI) system integrated with a TrueBeam linac was simulated with Monte Carlo EGSnrc/BEAMnrc, and the absorbed dose was calculated within PMMA, polyethylene (PE), and water head and body phantoms using EGSnrc/DOSXYZnrc, where the body PE body phantom emulated the ICRU/AAPM phantom. Beams of width 40-500 mm and beam qualities at tube potentials of 80-140 kV were studied. Application of a modified function of beam width (W) termed Gx(W), for which the cumulative dose for CBCT scans f (0) is normalized to the weighted CTDI (CTDIw) for a reference beam of width 40 mm, was investigated as a possible option. However, differences were found in Gx(W) with tube potential, especially for body phantoms, and these were considered to be due to differences in geometry between wide beams used for CBCT scans and those for conventional CT. Therefore, a modified function Gx(W)100 has been proposed, taking the form of values of f (0) at each position in a long phantom, normalized with respect to dose indices f 100(150)x measured with a 100 mm pencil ionization chamber within standard 150 mm PMMA phantoms, using the same scanning parameters, beam widths and positions within the phantom. f 100(150)x averages the dose resulting from a CBCT scan over the 100 mm length. Like the Gx(L) function, the Gx(W)100 function showed only a weak dependency on tube potential at most positions for the phantoms studied. The results were fitted to polynomial equations from which f (0) within the longer PMMA, PE, or water phantoms can be evaluated from measurements of f 100(150)x. Comparisons with other studies, suggest that these functions may be suitable for application to any CT or CBCT scan acquired with stationary table mode.

Radiological Protection in Cone Beam Computed Tomography (CBCT). ICRP Publication 129

The objective of this publication is to provide guidance on radiological protection in the new technology of cone beam computed tomography (CBCT). Publications 87 and 102 dealt with patient dose management in computed tomography (CT) and multi-detector CT. The new applications of CBCT and the associated radiological protection issues are substantially different from those of conventional CT. The perception that CBCT involves lower doses was only true in initial applications. CBCT is now used widely by specialists who have little or no training in radiological protection. This publication provides recommendations on radiation dose management directed at different stakeholders, and covers principles of radiological protection, training, and quality assurance aspects. Advice on appropriate use of CBCT needs to be made widely available. Advice on optimisation of protection when using CBCT equipment needs to be strengthened, particularly with respect to the use of newer features of the equipment. Manufacturers should standardise radiation dose displays on CBCT equipment to assist users in optimisation of protection and comparisons of performance. Additional challenges to radiological protection are introduced when CBCT-capable equipment is used for both fluoroscopy and tomography during the same procedure. Standardised methods need to be established for tracking and reporting of patient radiation doses from these procedures. The recommendations provided in this publication may evolve in the future as CBCT equipment and applications evolve. As with previous ICRP publications, the Commission hopes that imaging professionals, medical physicists, and manufacturers will use the guidelines and recommendations provided in this publication for implementation of the Commission's principle of optimisation of protection of patients and medical workers, with the objective of keeping exposures as low as reasonably achievable, taking into account economic and societal factors, and consistent with achieving the necessary medical outcomes.

A Monte Carlo investigation of cumulative dose measurements for cone beam computed tomography (CBCT) dosimetry

Many studies have shown that the computed tomography dose index (CTDI100) which is considered as a main dose descriptor for CT dosimetry fails to provide a realistic reflection of the dose involved in cone beam computed tomography (CBCT) scans. Several practical approaches have been proposed to overcome drawbacks of the CTDI100. One of these is the cumulative dose concept. The purpose of this study was to investigate four different approaches based on the cumulative dose concept: the cumulative dose (1) f(0,150) and (2) f(0,∞) with a small ionization chamber 20 mm long, and the cumulative dose (3) f100(150) and (4) f100(∞) with a standard 100 mm pencil ionization chamber. The study also aimed to investigate the influence of using the 20 and 100 mm chambers and the standard and the infinitely long phantoms on cumulative dose measurements. Monte Carlo EGSnrc/BEAMnrc and EGSnrc/DOSXYZnrc codes were used to simulate a kV imaging system integrated with a TrueBeam linear accelerator and to calculate doses within cylindrical head and body PMMA phantoms with diameters of 16 cm and 32 cm, respectively, and lengths of 150, 600, 900 mm. f(0,150) and f100(150) approaches were studied within the standard PMMA phantoms (150 mm), while the other approaches f(0,∞) and f100(∞) were within infinitely long head (600 mm) and body (900 mm) phantoms. CTDI∞ values were used as a standard to compare the dose values for the approaches studied at the centre and periphery of the phantoms and for the weighted values. Four scanning protocols and beams of width 20-300 mm were used. It has been shown that the f(0,∞) approach gave the highest dose values which were comparable to CTDI∞ values for wide beams. The differences between the weighted dose values obtained with the 20 and 100 mm chambers were significant for the beam widths <120 mm, but these differences declined with increasing beam widths to be within 4%. The weighted dose values calculated within the infinitely long phantoms with both the chambers for the beam widths ≤140 were within 3% of those within the standard phantoms, but the differences rose to be within 15% at wider beams. By comparing the approaches studied in this investigation with other methodologies taking into account the efficiency of the approach as a dose descriptor and the simplicity of the implementation in the clinical environment, the f(0,150) method may be the best for CBCT dosimetry combined with the use of correction factors.