Evaluation of In-room Volumetric Imaging Doses for Image-guided Radiotherapy: A Multi-institutional Study

Aims

We investigated imaging dose and noise under clinical scan conditions at multiple institutions using a simple and unified method, and demonstrated the need for diagnostic reference levels in image-guided radiotherapy (IGRT).

Materials and Methods

Nine cone-beam and helical computed tomography (CT) scanners (Varian, Elekta, Accuray Inc., and BrainLAB) from seven institutions were investigated in this study. The weighted cone-beam dose index (CBDIw) was calculated for head and pelvic protocols using a 100 mm pencil chamber under the conditions used in actual clinical practice at each institution. Cone-beam CT image noise was evaluated using polymethylmethacrylate head and body phantoms with diameters of 16 and 32 cm, respectively.

Results

For head and pelvic protocols, CBDIw values ranged from 0.94-6.59 and 1.47-20.9 mGy, respectively. Similarly, standard deviation (SD) values ranged from 9.3-34.0 and 26.9-97.4 HU, respectively. The SD values tended to increase with decreasing imaging dose (r = -0.33 and -0.61 for the head and pelvic protocols, respectively).

Conclusions

Among the nine machines, the imaging dose for high imaging dose institutions was approximately 20 mGy to the pelvic phantom, and there was a 14-fold difference in dose compared with the other institutions. These results suggest the need to establish DRLs for IGRT to guide clinical decision-making.

Dose estimation for cone-beam computed tomography in image-guided radiation therapy using mesh-type reference computational phantoms and assuming head and neck cancer

This study aimed to estimate the additional dose the cone-beam computed tomography (CBCT) system integrated into the Varian TrueBeam linear accelerator delivers to a patient with head and neck cancer using mesh-type International Commission on Radiological Protection reference computational phantoms. In the first part, for use as a benchmark for the accuracy of the Monte Carlo geometry of CBCT, Particle and Heavy Ion Transport code System (PHITS) calculations were confirmed against measured lateral and depth dose profiles using a computed tomography dose profiler. After obtaining good agreement, organ dose calculations were performed by PHITS using mesh-type reference computational phantom (MRCP) and irradiating the neck region; the effective dose was calculated utilising absorbed organ doses and tissue weighting factors for male and female MRCP. Substantially, it has been found that the effective doses for male and female MRCP are 0.81 and 1.06 mSv, respectively. As this study aimed to assess the imaging dose from the CBCT system used in image-guided radiation therapy, it is required to take into account this dose in terms of both the target organ and surrounding tissues. Although the absorbed organ dose values and effective dose values obtained for both MRCP males and females were small, attention should be paid to the additional dose resulting from CBCT. This study can create awareness on the importance of doses arising from imaging techniques, especially CBCT.

Secondary cancer risk from modern external-beam radiotherapy of prostate cancer patients: Impact of fractionation and dose distribution

Modern radiotherapy (RT) uses altered fractionation, long beam-on time and image-guided procedure. This study aimed to compare secondary cancer risk (SCR) associated with primary field, scatter/leakage radiations and image-guided procedure in prostate treatment using intensity-modulated RT (IMRT), CyberKnife stereotactic body RT (CK-SBRT) in relative to 3-dimensional conformal RT (3D-CRT). Prostate plans were generated for 3D-CRT, IMRT (39 fractions of 2 Gy), and CK-SBRT (five fractions of 7.25 Gy). Excess absolute risk (EAR) was calculated for organs in the primary field using Schneider's mechanistic model and concept of organ equivalent dose (OED) to account for dose inhomogeneity. Doses from image-guided procedure and scatter/leakage radiations were determined by phantom measurements. The results showed that hypofractionation relative to conventional fractionation yielded lower SCR for organs in primary field (p ≤ 0.0001). SCR was further modulated by dose-volume distribution. For organs near the field edge, like the rectum and pelvic bone, CK-SBRT plan rendered better risk profiles than IMRT and 3D-CRT because of the absence of volume peak in high dose region (relative risk [RR]: 0.65, 0.22, respectively, p ≤ 0.0004). CK-SBRT and IMRT generated more scatter/leakage and imaging doses than 3D-CRT (p ≤ 0.0002). But primary field was the major contributor to SCR. EAR estimates (risk contributions, primary field: scatter/leakage radiations: imaging procedure) were 7.1 excess cases per 104 person-year (PY; 3.64:2.25:1) for CK-SBRT, 9.93 (7.32:2.33:1) for IMRT and 8.24 (15.99:2.35:1) for 3D-CRT (p ≤ 0.0002). We conclude that modern RT added more but small SCR from scatter/leakage and imaging doses. The primary field is a major contributor of risk which can be mitigated by the use of hypofractionation.

Evaluation of IGRT-Induced Imaging Doses and Secondary Cancer Risk for SBRT Early Lung Cancer Patients In Silico Study

Objectives:

This study performed dosimetry studies and secondary cancer risk assessments on using electronic portal imaging device (EPID) and cone beam computed tomography (CBCT) as image guided tools for the early lung cancer patients treated with SBRT.

Methods:

The imaging doses from MV-EPID and kV-CBCT of the Edge accelerator were retrospectively added to sixty-one SBRT treatment plans of early lung cancer patients. The MV-EPID imaging dose (6MV Photon beam) was calculated in Pinnacle TPS, and the kV-CBCT imaging dose was simulated and calculated by modeling of the kV energy beam in TPS using Pinnacle automatic modeling program. Three types of plans, namely Plan_EPID, Plan_CBCT and Plan_origin, were generated with incorporating doses of EPID, CBCT and no imaging, respectively, for analysis. The effects of imaging doses on dose-volume-histogram (DVH) and plan quality were analyzed, and the excess absolute risk (EAR) of secondary cancer for ipsilateral lung was evaluated.

Results:

The regions that received less than 50 cGy were significantly impacted by the imaging doses, while the isodose lines greater than 1000 cGy were barely changed. The DVH values of ipsilateral lung increased the most in Plan_EPID, followed by Plan_CBCT. Compared to Plan_origin on the average, the estimated EAR of ipsilateral lung in Plan_EPID increased by 3.43%, while the corresponding EAR increase in Plan_CBCT was much smaller (about 0.4%). Considering only the contribution of the imaging dose, the EAR values for the ipsilateral lung due to the MV-EPID dose in 5 years,10 years and 15 years were 1.49 cases, 2.09 cases and 2.88 cases per 10⁴PY respectively, and those due to the kV-CBCT dose were about 9 times lower, correspondingly.

Conclusions:

The imaging doses produced by MV-EPID and kV-CBCT had little effects on the target dose coverage. The secondary cancer risk caused by MV-EPID dose is more than 8.5 times that of kV-CBCT.

How low can you go? A CBCT dose reduction study

Purpose

Cone beam computed tomography (CBCT) is often used for patient setup based solely on bony anatomy. The goal of this work was to evaluate whether CBCT dose can be lowered to the level of kV image pair doses when used for bony anatomy‐based IGRT without compromising positioning accuracy.

Methods

An anthropomorphic phantom was CT scanned in the head, head and neck, chest, and pelvis regions and setup on the linear accelerator couch with the isocenter near the planned location. Cone beam computed tomographies were performed with the standard “full dose” protocol supplied by the linac vendor. With sequentially lowering the dose, three‐dimensional (3D) matching was performed for each without shifting the couch. The standard kV image pair protocol for each site was also used to image the phantoms. For all studies, six degrees of freedom was included in the 2D or 3D matching to the extent they could be employed. Imaging doses were determined in air at isocenter following the TG‐61 formalism.

Results

Cone beam computed tomography dose was reduced by 81–98% of the standard CBCT protocol to nearly that of the standard kV image pair dose for each site. Relative to the standard CBCT shift values, translational shifts were within 0.3 and 1.6 mm for all sites, for the reduced dose CBCT and kV image pair, respectively. Rotational shifts were within 0.2 degree and 0.7 degrees for all sites, for the reduced dose CBCTs and kV image pair, respectively.

Conclusion

For bony anatomy‐based image guidance, CBCT dose can be reduced to a value similar to that of a kV image pair with similar or better patient positioning accuracy than kV image pair alignment. Where rotations are important to correct, CBCT will be superior to orthogonal kV imaging without significantly increased imaging dose. This is especially important for image guidance for pediatric patient treatments.

Radiation protection in radiological imaging: a survey of imaging modalities used in Japanese institutions for verifying applicator placements in high-dose-rate brachytherapy

Institutional imaging protocols for the verification of brachytherapy applicator placements were investigated in a survey study of domestic radiotherapy institutions. The survey form designed by a free on-line survey system was distributed via the mailing-list system of the Japanese Society for Radiation Oncology. Survey data of 75 institutions between August 2019 and October 2019 were collected. The imaging modalities used were dependent on resources available to the institutions. The displacement of a brachytherapy applicator results in significant dosimetric impact. It is essential to verify applicator placements using imaging modalities before treatment. Various imaging modalities used in institutions included a computed tomography (CT) scanner, an angiography X-ray system, a multi-purpose X-ray system and a radiotherapy simulator. The median total exposure time in overall treatment sessions was $\le$75 s for gynecological and prostate cancers. Some institutions used fluoroscopy to monitor the brachytherapy source movement. Institutional countermeasures for reducing unwanted imaging dose included minimizing the image area, changing the imaging orientation, reducing the imaging frequency and optimizing the imaging conditions. It is worth noting that half of the institutions did not confirm imaging dose regularly. This study reported on the usage of imaging modalities for brachytherapy in Japan. More caution should be applied with interstitial brachytherapy with many catheters that can lead to potentially substantial increments in imaging doses for monitoring the actual brachytherapy source using fluoroscopy. It is necessary to share imaging techniques, standardize imaging protocols and quality assurance/quality control among institutions, and imaging dose guidelines for optimization of imaging doses delivered in radiotherapy should be developed.

Evaluation of Dose Distribution and Normal Tissue Complication Probability of a Combined Dose of Cone-Beam Computed Tomography Imaging with Treatment in Prostate Intensity-Modulated Radiation Therapy

Purpose:

The purpose of this study is to evaluate the effects of cone-beam computed tomography (CBCT) on dose distribution and normal tissue complication probability (NTCP) by constructing a comprehensive dose evaluation system for prostate intensity-modulated radiation therapy (IMRT).

Methods:

A system that could combine CBCT and treatment doses with MATLAB was constructed. Twenty patients treated with prostate IMRT were studied. A mean dose of 78 Gy was prescribed to the prostate region, excluding the rectal volume from the target volume, with margins of 4 mm to the dorsal side of the prostate and 7 mm to the entire circumference. CBCT and treatment doses were combined, and the dose distribution and the NTCP of the rectum and bladder were evaluated.

Results:

The radiation dose delivered to 2% and 98% of the target volume increased by 0.90 and 0.74 Gy on average, respectively, in the half-fan mode and on average 0.76 and 0.72 Gy, respectively, in the full-fan mode. The homogeneity index remained constant. The percent volume of the rectum and bladder irradiated at each dose increased slightly, with a maximum increase of <1%. The rectal NTCP increased by approximately 0.07% from 0.46% to 0.53% with the addition of a CBCT dose, while the maximum NTCP in the bladder was approximately 0.02%.

Conclusions:

This study demonstrated a method to evaluate a combined dose of CBCT and a treatment dose using the constructed system. The combined dose distribution revealed increases of <1% volume in the rectal and bladder doses and approximately 0.07% in the rectal NTCP.

Characterization of the Megavoltage Cone-Beam Computed Tomography (MV-CBCT) System on HalcyonTM for IGRT: Image Quality Benchmark, Clinical Performance, and Organ Doses

Purpose: The Varian Halcyon includes an ultrafast 6 MV flattening filter free (FFF) cone-beam computed tomography (MV-CBCT). Although a kV-CBCT add-on is available, in the basic configuration MV is used for image guided radiotherapy (IGRT). We characterized the MV-CBCT imager in terms of reproducibility, linearity, field of view (FOV) dependence, detectability of soft-tissue, and the effect of metal implants. The performance of the MV-CBCT in the clinic, including resulting dose to organs, is also discussed herein.

Methods: A Gammex phantom was scanned using a Halcyon MV-CBCT and a 120 kVp Siemens Definition Edge CT. Mean and standard deviation of Hounsfield Units (HUs) for different electron density relative to water (ρeW) inserts were extracted. Doses to clinical patients due to MV-CBCT are calculated within Eclipse during treatment planning.

Results: A stable and near-linear HU-to-ρeW curve was obtained using the MV-CBCT. As the scan length increased from 10 to 28cm, the linearity of curve improved while the mean HUs decreased by 30%. All soft tissue inserts in the Gammex phantom were distinguishable. A crescent artifact affected HU measurements by up to 40 HUs. Soft-tissue contrast was sufficient for clinical online image-guidance in the low dose (5 MU) mode. Mean doses per fraction to organs-at-risk (OARs) were as high as 6 cGy for head and neck, 5 cGy for breast, and 4 cGy for pelvis patients. Metal rods did not affect HU values or introduce noticeable artifacts.

Conclusions: Halcyon's MV-CBCT has sufficient soft tissue contrast for IGRT and lacks metal-induced artifacts. Even though the absolute HU values vary with phantom size and scanning length, the HU-to-ρeW conversions are linear and stable day-to-day. In clinical cases, highest tissue doses from MV-CBCT ranged from 2-7cGy per fraction for various treatment sites, which could be significant for some organs at risk. Dose to out-of-treatment-field organs can be limited by reducing the scan length definition during planning and using the low dose mode. The high quality imaging mode did not provide material advantages over the low dose mode. Adequate IGRT was successfully delivered to multiple tumor sites using MV-CBCT.

Adaptive Imaging Versus Periodic Surveillance for Intrafraction Motion Management During Prostate Cancer Radiotherapy

Objective:

To evaluate the benefits of adaptive imaging with automatic correction compared to periodic surveillance strategies with either manual or automatic correction.

Methods:

Using Calypso trajectories from 54 patients with prostate cancer at 2 institutions, we simulated 5-field intensity-modulated radiation therapy and dual-arc volumetric-modulated arc therapy with periodic imaging at various frequencies and with continuous adaptive imaging, respectively. With manual/automatic correction, we assumed there was a 30/1 second delay after imaging to determine and apply couch shift. For adaptive imaging, real-time “dose-free” cine-MV images during beam delivery are used in conjunction with online-updated motion pattern information to estimate 3D displacement. Simultaneous MV-kV imaging is only used to confirm the estimated overthreshold motion and calculate couch shift, hence very low additional patient dose from kV imaging.

Results:

Without intrafraction intervention, the prostates could on average have moved out of a 3-mm margin for ∼20% of the beam-on time after setup imaging in current clinical situation. If the time interval from the setup imaging to beam-on can be reduced to only 30 seconds, the mean over-3 mm percentage can be reduced to ∼7%. For intensity-modulated radiation therapy simulation, with manual correction, 110 and 70 seconds imaging periods both reduced the mean over-3 mm time to ∼4%. Automatic correction could give another 1% to 2% improvement. However, with either manual or automatic correction, the maximum patient-specific over-3 mm time was still relatively high (from 6.4% to 12.6%) and those patients are actually clinically most important. In contrast, adaptive imaging with automatic intervention significantly reduced the mean percentage to 0.6% and the maximum to 2.7% and averagely only ∼1 kV image and ∼1 couch shift were needed per fraction. The results of volumetric-modulated arc therapy simulation show a similar trend to that of intensity-modulated radiation therapy.

Conclusions:

Adaptive continuous monitoring with automatic motion compensation is more beneficial than periodic imaging surveillance at similar or even less imaging 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.

Quantification of the kV X-ray imaging dose during real-time tumor tracking and from three- and four-dimensional cone-beam computed tomography in lung cancer patients using a Monte Carlo simulation

Knowledge of the imaging doses delivered to patients and accurate dosimetry of the radiation to organs from various imaging procedures is becoming increasingly important for clinicians. The purposes of this study were to calculate imaging doses delivered to the organs of lung cancer patients during real-time tumor tracking (RTTT) with three-dimensional (3D), and four-dimensional (4D) cone-beam computed tomography (CBCT), using Monte Carlo techniques to simulate kV X-ray dose distributions delivered using the Vero4DRT. Imaging doses from RTTT, 3D-CBCT and 4D-CBCT were calculated with the planning CT images for nine lung cancer patients who underwent stereotactic body radiotherapy (SBRT) with RTTT. With RTTT, imaging doses from correlation modeling and from monitoring of imaging during beam delivery were calculated. With CBCT, doses from 3D-CBCT and 4D-CBCT were also simulated. The doses covering 2-cc volumes (D2cc) in correlation modeling were up to 9.3 cGy for soft tissues and 48.4 cGy for bone. The values from correlation modeling and monitoring were up to 11.0 cGy for soft tissues and 59.8 cGy for bone. Imaging doses in correlation modeling were larger with RTTT. On a single 4D-CBCT, the skin and bone D2cc values were in the ranges of 7.4-10.5 cGy and 33.5-58.1 cGy, respectively. The D2cc from 4D-CBCT was approximately double that from 3D-CBCT. Clinicians should Figure that the imaging dose increases the cumulative doses to organs.

A Projection Quality-Driven Tube Current Modulation Method in Cone-Beam CT for IGRT: Proof of Concept

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 μmaxθ within the 95% range of each projection angle (θ) was estimated according to the planning computed tomography images. Then, a desired 2-dimensional contrast to noise ratio value was selected, and the mAs setting at θ was calculated with the lookup table of mAs versus 2-dimensional contrast to noise ratio,μmaxθ. Three-dimensional cone-beam computed tomography images were reconstructed using the projections acquired with the selected mAs. The imaging dose was evaluated with a polymethyl methacrylate dosimetry phantom in terms of volume computed tomography dose index. Image quality was analyzed using a Catphan 503 phantom with an oval body annulus and a pelvis phantom.

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.

A motion-compensated image filter for low-dose fluoroscopy in a real-time tumor-tracking radiotherapy system

In the real-time tumor-tracking radiotherapy system, a surrogate fiducial marker inserted in or near the tumor is detected by fluoroscopy to realize respiratory-gated radiotherapy. The imaging dose caused by fluoroscopy should be minimized. In this work, an image processing technique is proposed for tracing a moving marker in low-dose imaging. The proposed tracking technique is a combination of a motion-compensated recursive filter and template pattern matching. The proposed image filter can reduce motion artifacts resulting from the recursive process based on the determination of the region of interest for the next frame according to the current marker position in the fluoroscopic images. The effectiveness of the proposed technique and the expected clinical benefit were examined by phantom experimental studies with actual tumor trajectories generated from clinical patient data. It was demonstrated that the marker motion could be traced in low-dose imaging by applying the proposed algorithm with acceptable registration error and high pattern recognition score in all trajectories, although some trajectories were not able to be tracked with the conventional spatial filters or without image filters. The positional accuracy is expected to be kept within ±2 mm. The total computation time required to determine the marker position is a few milliseconds. The proposed image processing technique is applicable for imaging dose reduction.