Abstract
Purpose:
To develop a projection quality-driven tube current modulation method in cone-beam computed tomography for image-guided radiotherapy based on the prior attenuation information obtained by the planning computed tomography and then evaluate its effect on a reduction in the imaging dose.
Materials and Methods:
The QCKV-1 phantom with different thicknesses (0-400 mm) of solid water upon it was used to simulate different attenuation (μ). Projections were acquired with a series of tube current–exposure time product (mAs) settings, and a 2-dimensional contrast to noise ratio was analyzed for each projection to create a lookup table of mAs versus 2-dimensional contrast to noise ratio, μ. Before a patient underwent computed tomography, the maximum attenuation
Results:
For the Catphan 503 phantom, the cone-beam computed tomography image obtained by the projection quality-driven tube current modulation method had a similar quality to that of conventional cone-beam computed tomography . However, the proposed method could reduce the imaging dose by 16% to 33% to achieve an equivalent contrast to noise ratio value. For the pelvis phantom, the structural similarity index was 0.992 with a dose reduction of 39.7% for the projection quality-driven tube current modulation method.
Conclusions:
The proposed method could reduce the additional dose to the patient while not degrading the image quality for cone-beam computed tomography. The projection quality-driven tube current modulation method could be especially beneficial to patients who undergo cone-beam computed tomography frequently during a treatment course.
Introduction
Kilovoltage cone-beam computed tomography (CBCT) installed on linear accelerators can provide detailed volumetric information on patient anatomy at the treatment position. It allows online verification of the 3-dimensional (3D) patient setup for image-guided radiotherapy (IGRT) 1 –3 and replanning for adaptive radiotherapy. 4 –6 As a result, CBCT is now used widely to guarantee highly accurate delivery of radiotherapy.
However, repeated use of CBCT in the course of radiotherapy produces a considerable radiation dose to patients, 7,8 which has aroused growing concern about a potential radiation risk. 9 –11 One of the main applicable methods to reduce the dose of CBCT is application of a low tube current–exposure time product (mAs) setting. Some studies 12,13 have investigated the image quality and imaging dose together as functions of mAs in a qualitative manner to choose the optimal setting. However, a low mAs scheme results in a reduced signal to noise ratio in the projection images due to fewer incident photons interacting with detectors. Accordingly, the reconstructed images are noisy and the image quality degraded. Moreover, if the patient is large (ie, when the X-rays are highly attenuated by an object), there will be streak artifacts due to “photon starvation” 14 at low mAs imaging, whereby the number of detected photons per X-ray is too low.
Cone-beam computed tomography systems for IGRT available in the clinic now adopt a constant mAs setting in a circular scan for a specific site. However, the patient’s anatomy is usually not cylindrical, especially in cases of imaging of the abdomen and/or pelvis. The attenuation of X-rays varies greatly at different scanning angles, which may result in higher mAs than that required in the short-X-ray penetration path or lower mAs in the long-X-ray penetration path. To reduce the patient dose without compromising image quality, automatic tube current modulation (ATCM) has been implemented on commercial systems of fan-beam computed tomography (FBCT). 15 –18 The first requirement for the operation of ATCM is the determination of patient attenuation, and this is primarily obtained from the scan projection radiograph (SPR). The methods of ATCM operation are different in CT scanners from different manufacturers. Scanners from Philips (Amsterdam, the Netherlands) and Siemens (Munich, Germany) use the concepts of reference image and reference mAs. However, scanners produced by Toshiba (Tokyo, Japan) and General Electric (Boston, Massachusetts) base the current modulation on a target noise setting. These ATCM systems modulate the X-ray tube current according to the patient’s anatomy along 2 directions. One modulation, named the “longitudinal” or “z” modulation, is based on different attenuation settings along the longitudinal axis of the patient. The other modulation, called the “angular” or “x-y” modulation, is used to modify the X-ray tube current during rotation around the patient. However, ATCM technology has yet to be applied to CBCT for IGRT.
The aim of the present study was to develop a projection quality-driven tube current modulation (PQD-TCM) method for CBCT. Although noise is the obvious parameter used as the image quality reference, the ability to detect low-contrast tissue in CBCT image depends on the contrast to noise ratio (CNR). Unlike the existing noise-based ATCM methods, we used CNR as the image quality reference to highlight the low-contrast visibility of CBCT. In addition, the patient attenuation could be exactly calculated with the planning CT (pCT) acquired before the use of CBCT, rather than the estimation from SPR. Performance of the proposed method was evaluated in terms of the imaging dose and image quality.
Materials and Methods
Phantoms
Four phantoms available commercially were used in the present study (Figure 1). The QCkV-1 phantom (Standard Imaging, Middleton, Wisconsin) is a dedicated kilovolt X-ray phantom used to analyze the image quality of a 2-dimensional (2D) projection (Figure 1A). It was used to measure and analyze the CNR for 2D projection images (CNR2D). The dosimetry phantom shown in Figure 1B is made of polymethyl methacrylate with a diameter of 32 cm, which meets the specification set by the International Electrotechnical Commission. 19 It was used to measure the computed tomography dose index (CTDI) of CBCT imaging.

The 4 phantoms: (a) QCkV-1, (b) dosimetry, (c) Catphan 503, and (d) pelvis.
Catphan 503 (Phantom Laboratory, Salem, New York) was used to evaluate image quality in a quantitative manner. This phantom is cylindrical with a diameter of 20 cm and a length of 20 cm. To simulate the human body, the Catphan 503 phantom is inserted in an oval body annulus with a major axis of 35 cm and minor axis of 25 cm (Figure 1C).
The pelvis phantom (Phantom Laboratory) consists of soft tissue–equivalent materials and bone-equivalent materials (Figure 1D). Its absorption and scattering properties to X-rays are equivalent to those of human tissue. It was used to simulate a noncircular human body.
Cone-Beam Computed Tomography System
Cone-beam computed tomography was done using the X-ray volumetric imaging (XVI) of a Versa HD machine (Elekta, Stockholm, Sweden). The X-ray source uses a rotating anode X-ray tube (D604; Dunlee, Aurora, Illinois) with peak tube potential of 150 kV and maximum current of 500 mA. The detector is an indirect-detection, flat-panel imager with a spatial resolution of 1024 × 1024 arrays of 0.4 × 0.4 mm2 pixels (RID1640-A11; PerkinElmer, Waltham, Massachusetts). The source-to-axis distance and source-to-detector distance are 1000 and 1536 mm, respectively. Projections can be acquired with small, medium, and large field of view (FOV). The number of projections for a full 360° rotation is ≈660. X-ray volumetric imaging software uses a cone-beam reconstruction process based on the Feldkamp-Davis-Kress algorithm. A M20 collimator was applied for the body size phantoms. Projections were processed at high resolution to yield projections of dimension 512 × 512 with pixel size of 0.8 × 0.8 mm2. Cone-beam computed tomography images were reconstructed with a voxel size of 1.0 × 1.0 × 1.0 mm3.
Projection Quality-Driven Tube Current Modulation Method
For patients with noncylindrical anatomy, the attenuation of X-rays varies greatly at different scanning angles. A higher attenuation in the projection introduces more noise, which degrades the image quality of reconstructed CT. We proposed a PQD-TCM method to modulate the current based on the image quality of the projection.
The implementation procedure of the proposed PQD-TCM method is shown in Figure 2. It had 4 steps, as follows.

The implementation procedure of the projection quality-driven tube current modulation (PQD-TCM) method.
Step 1. Create a basic lookup table
The QCkV-1 phantom was placed in the center of the FOV with different thicknesses (0-400 mm) of solid water upon it to simulate different attenuation (μ) values. Projections (p) were acquired with a series of mAs settings, and the CNR2D was analyzed for each projection. Then, a lookup table relating mAs versus CNR2D, μ was created.
Step 2. Estimate the attenuation based on the prior pCT
Patients who underwent radiotherapy usually need to have a CT simulation to acquire their 3D pCT images for treatment planning. In the course of subsequent treatment, CBCT acquired at the treatment position is registered to the pCT for patient’s setup verification. Due to the same setup position of patient in the acquisition of pCT and CBCT, we can have an accurate attenuation (μ) estimation of the patient based on the pCT before CBCT scan. In this step, we first set the treatment center (isocenter) in pCT as the CBCT center. Then the 2D attenuation (μ) maps of each CBCT projection angle (θ) were acquired by projecting the pCT images. Finally, the maximum attenuation
Step 3. Calculate the mAs with the desired goal
This method allows the user to set a desired CNR2D value, which is used to estimate the tube current. Then, the mAs setting at θ is calculated with the lookup table using mAs versus CNR2D,
Step 4. Reconstruct 3D CBCT images using the projections acquired with the selected mAs
The CBCT system available could not change the mAs during a full circular scan. To simulate a series of scans with the selected mAs setting, reconstruction schemes with these projections were achieved by manually selecting the desirable projection frames in the XVI database.
Performance Evaluation
The image quality and imaging dose were measured and analyzed with and without the PQD-TCM method.
Analyses of image quality
The Catphan 503 phantom with the oval body annulus was used to evaluate the image quality in a quantitative manner. The CNR2D values of 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, and 0.8 were tested for the PQD-TCM method. Conventional CBCT using a set of constant mAs settings (ie, mAs = 2, 3, 4, 5.12, 6.4, 10, 12.8, and 16 per frame) were also studied for comparison. The indices for the quantitative comparison focused on 3 evaluation metrics including the CNR, the high-contrast resolution (HCR), and the noise level (noise).
Contrast to noise ratio was calculated using the CTP 404 module (Figure 3A). It has 8 embedded rods comprising different materials: air (2 rods), polymethylpentene, low-density polyethylene (LDPE), polystyrene (PS), acrylic, polyoxymethylene (Delrin), and polytetrafluoroethylene (Teflon). Regions of interest (ROIs) of size 4 × 4 mm2 within PS and LDPE inserts were used to measure the mean and standard deviation in Hounsfield units. The CNR was calculated using the following equation:

Three modules in Catphan 503 phantom: (A) CTP-404, (B) CTP-528, and (C) CTP-486.
where MeanPS and MeanLDPE are the mean voxel values in PS and LDPE, respectively, and
High-contrast resolution was measured using the CTP 528 module (Figure 3B and C). It was evaluated with the modulation transfer function (MTF) using a well-known method. 20 The evaluation value was 50% of the MTF value (MTF50%) from the MTF curves. Larger MTF50% value represents better HCR.
Image noise measurements were performed using the CTP 486 module (Figure 3D), which is uniform and water equivalent. The standard deviation (SD) in 5 ROIs of size 10 × 10 mm2 (including 1 at the center and 4 at the peripheral positions) was measured. Image noise was calculated as the mean of the 5 SD values.
Image quality was also evaluated using the pelvis phantom. The conventional CBCT used a default mAs setting in the clinic (ie, 16 mAs/frame), whereas the CNR2D value was selected to ensure a maximum mAs of 16 mAs/frame for the PQD-TCM method. The CBCT images obtained by the 2 methods were analyzed using the structural similarity (SSIM) index. The SSIM index is used to measure the similarity between 2 images
where
Measurement of imaging dose
For estimation of the radiation exposure of CBCT imaging, the volume CT dose index (
where
Results
Catphan 503
Figure 4A shows maximum attenuation of Catphan 503 phantom of each projection, and Figure 4B shows the corresponding mAs selection for CNR2D = 0.4 as an example. According to our method, the higher mAs should be selected at the projection angle with larger attenuation.

The maximum attenuation (A) and mAs selection for 2-dimensional contrast to noise ratio (CNR2D) = 0.4 (B) of each projection.
Figure 5 shows the mAs selection of each projection angle for the PQD-TCM method. The selected mAs increased with increases in the desired CNR2D setting. It also varied according to the attenuation of different projection angles at a fixed CNR2D setting. For some projection angles with a short penetrating path, a high mAs was not needed.

The mAs selection of each projection angle. (A), (B), (C), (D), (E), (F), (G) and (H) are for 2-dimensional contrast to noise ratio (CNR2D) = 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, and 0.8, respectively.
The results for the measured
Results Using the Catphan 503 Phantom.
Abbreviations: CNR, contrast to noise ratio;

Image quality as a function of imaging dose: (A) MTF50%; (B) contrast to noise ratio (CNR); (C) noise.

Cone-beam computed tomography (CBCT) images of Catphan 503: (A) Conventional method; (B) projection quality-driven tube current modulation (PQD-TCM) method.
Pelvis Phantom
The CBCT images for the pelvis phantom are shown in Figure 8. The outline of the phantom was not correct for conventional CBCT because of photon saturation through the thin edge. The overall SSIM index was 0.992 ± 0.001 for images acquired by the 2 methods. The

Cone-beam computed tomography (CBCT) images of the pelvis phantom: (A) and (B): Conventional method; (C) and (D): projection quality-driven tube current modulation (PQD-TCM) method.
Discussion
We have designed a PQD-TCM method in CBCT for IGRT. To the best of our knowledge, this task has not previously been reported. We compared the proposed method with the conventional method. This comparison reveals that the CBCT image obtained by the PQD-TCM method had a better low-contrast visibility, lower noise, and similar high-contrast visibility compared with that of conventional CBCT with the same imaging dose level. It could be beneficial to patients who undergo CBCT frequently.
Treatment plans in radiotherapy are designed to deliver a high radiation dose to the tumor target and minimize the radiation dose outside the treatment volume. This sharp dose distribution requires accurate positioning of the patient on the treatment couch. To ensure the precision of radiotherapy, CBCT is used frequently for position verification. If we assume that the average prescription dose is 200 cGy per fraction, the percentage dose from CBCT imaging to megavoltage treatment would be ≈ 1% for the body, which is relatively small. However, the “as low as reasonably achievable” principle is always recommended by the radiology community. The strategy presented here shows that there is a potential to reduce the CBCT dose significantly without sacrificing image quality.
Although the ATCM is available in FBCT now, it is different from the proposed PQD-TCM for CBCT in 3 ways. First, CBCT cannot modulate the current in the longitudinal direction because of its cone-shaped beam. Second, the tube current values per rotation are calculated prospectively using an SPR for FBCT, but the proposed method has the advantage of using the prior attenuation information obtained by the pCT. Third, the tube current was selected to achieve the standard or reference noise level for the ATCM systems of FBCT, whereas the PQD-TCM method calculates the tube current based on the CNR of 2D projections.
Tube current modulation technology is not available on current CBCT systems for IGRT, which use a constant current setting during a scan. To combine a group of projections acquired with the desired current at a certain gantry angle, we first scanned the phantoms with a series of current settings to create a database and then extracted the desired projections from the database to reconstruct 3D images. Although for the dose measurement, we measured the imaging dose of the projections with desired mAs setting separately and then summed them up manually. The methodology described here could aid the design of an imager that would provide tube current modulation technology.
In the present study, the CBCT system of Elekta was used, but it could be applied in other CBCT systems as well. The M20 collimator and F0 filter were tested in the present study as an example. However, analyses of other collimators and filters would be worthwhile to complement the results of the present study in the future. There are some differences between phantoms and human body. The heterogeneity of each person also does exist. Although the experiment results on the 2 phantoms show our method works well, more researches should be carried out before the clinical use.
The scanning FOV is much larger for cone-beam scanning, so it induces more scattering, which is well known to be the most important factor affecting the quality of CBCT images. The experiments demonstrate the proposed method can reduce the imaging dose to patients while not degrading the image quality, although the scatter wasn’t corrected. However, scatter may lead to inaccurate dose evaluation. Some researchers have developed different correction strategies. One study
21
estimated the scatter directly in each projection from pixel values near the edge of the detector behind the collimator leaves. Wu
Conclusion
This work describes a PQD-TCM framework that utilizes the prior information of PCT. The image quality and imaging dose were investigated. We showed that it is possible to reduce the additional dose to a patient while not degrading the image quality. These observations suggest that the proposed method could be applied to the frequently used CBCT in radiotherapy.
Footnotes
Abbreviations
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) declared the following potential conflicts of interest with respect to the research, authorship, and/or publication of this article: This work was supported by the National Natural Science Foundation of China (grant numbers 11605291;11475261) and the National Key Projects of Research and Development of China (grant number 2016YFC0904600; 2017YFC0107501).
