Abstract
This study aimed to evaluate the effect of image distortion and Hounsfield unit (HU) variation due to the extended field-of-view (eFOV) of the large-bore (LB) computed tomography (CT) on dose distribution. Both home-made inhomogeneity and breast phantoms were scanned at the geometric center position and four different offset positions. We also performed dose optimizations based on different breast phantom CT sets for evaluating the effects of image artifacts on the intensity-modulated radiation techniques. The volume changes were 0.0% to 0.5% in the air, −0.5% to 3.0% in the water, and 4.0% to 5.0% in the high-density material of the inhomogeneity phantom. Both phantoms scanning results indicate that more distortions occurred in the eFOV area due to the biased scanning center. The gamma index differences ranged from 0.87% to 4.87% for the FIF plan and from 0.52% to 6.26% for the VMAT plan. This resulted in decrease of the minimum (7.3–13.1%), maximum (−0.8–2.2%), and mean doses (−0.2–4.4%). We recommend that it should be evaluated whether the applied CT would have an appropriate eFOV range for clinical radiation treatment planning for patients.
Similar content being viewed by others
Introduction
In modern radiation therapy simulation, computed tomography (CT) is essential for delineating critical organs and targets, combined with other high-resolution imaging modalities, such as magnetic resonance imaging (MRI) and positron emission tomography (PET). In particular, the use of CT images enables clinicians to accommodate various electron densities of the tissues for calculating dose distributions with heterogeneity correction during dose optimization1. CT image calibration, which stands for electron density vs. Hounsfield Unit (HU), must be performed prior to their clinical application in the radiation treatment planning2,3. Generally, a large bore (LB) size allows for applying flexible set-up when a patient is scanned with sophisticated immobilization devices, such as breast tilting boards. The commercial LB models of various CT simulators provide larger bore diameters, ranging from 80 to 90 cm, than those of the diagnostic CT scanners (typically 65–70 cm)4,5. The LB model CTs have an extended range of HU value (−32738 to 32767) as compared to the HU range (1024 to 8191) used in the scanned images depending on the manufacturers. It also provides extended field-of-view (eFOV), larger than the scan field-of-view (sFOV) from the reconstruction algorithm6.
If the treated lesions are located laterally with respect to the center of the body, such as breast cancer lesions, the treatment area could be close to the edge of or outside the sFOV. In such cases, the eFOV algorithm could make image distortion and HU errors due to the partial sampling of the scanning data between the diameter of the sFOV and that of the eFOV7.
The image reconstruction algorithm provided by most CT scanners uses either a Filtered Back Projection (FBP) or a repetitive reconstruction method. As the FBP algorithm is a method of reconstructing scanning data within the range of the sFOV, incomplete data provided by the eFOV could result in image distortion or undesirable HU variation8,9.
It could lead to inaccurate and imprecise dose calculation based on the incorrectly reconstructed CT set in the treatment planning10,11.
In this study, we created standardized radiation treatment plans (RTPs) using tangential irradiation techniques, such as Field-in-Field (FIF) and Volumetric Modulated Arc Therapy (VMAT), based on the different CT image sets.
We analyzed the effect of image distortion and HU variation due to the eFOV on the calculated dose distribution for patients with breast cancer.
Results
Volume deformation and HU value variation
Table 1 shows the differences between the contoured volumes of air, water, and high-density material at the different scanned center positions, which are center, Off-20 cm, and Off-30 cm. For the reference position of scanning at the CT bore center, volume changes showed differences of 0.0–0.5% in the air, −0.5–3.0% in the water, and 4.0–5.0% in the high-density material. There were greater differences in the high-density material than in the air and water and in the Off-30 cm than in the Off-20 cm. These results indicate that the reconstructed volumes could be affected by its locations and densities in the vicinity of the range between an sFOV of 75 cm and an eFOV of 85 cm diameters. For the breast phantom, the variation of contoured volumes of the lung, heart, and left breast are listed in Table 2. Depending on the off displacements from the geometrical center, the volumes changed in the range of −0.3–12.5%, 0.1–0.5%, and 0.5–12.8% for the lung, heart, and left breast PTV, respectively. The results from both phantom scans indicate that many distortions may have occurred near the predicted eFOV boundary of Off-25 cm or 30 cm.
As depicted in Fig. 1, more severe variations of the HU profiles were found in the Off-30 cm CT set as compared to the center CT set, mainly in the high-density material. Table 3 shows that the HU values tended to increase in the lung and heart, while those of the left breast tended to fluctuate at the Off-25 cm and Off-30 cm. These results indicate that the HU values tend to increase with the increasing of the offset distance. However, if the offset distance is getting closer to the eFOV beyond the sFOV, the HU could be decreased.
Effects on dose distributions of breast phantom plans
In the comparison of FIF plans based on the different offset-CT, the maximum doses were approximately 2% higher in the Off-25 cm and Off-30 cm than that in the center. Other than Off-30 cm, other offset plans illustrated very similar line profiles of dose distribution in the breast PTV area (Fig. 2). This finding reflects similar results from reconstruction errors of HU variations and image distortions. As shown in Table 4, in comparison to the center scanning image, GI differences ranged from 0.87 to 4.87% for the FIF plans and from 0.52 to 6.26% for the VMAT plan. The trend seen in these results also shows that the differences in GI tend to increase with increasing of the offset distances. As depicted in Figs. 3–5, the differential dose-volume histograms (DDVHs) show the effect of the image distortion on the dose distributions in the lung, heart and left breast in VMAT plans. There are many discrepancies in the range of 1000 to 5000 cGy in the DDVH of the lung (Fig. 3), DDVHs of the heart show nearly identical dose distributions (Fig. 4). A possible explanation could be the influence of the deformed lung portion in the high dose region close to the breast PTV. Table 5 and Fig. 5 show the dose variations according to the breast volume changes and a differential dose-volume comparison of breast PTV. As seen in Table 5, as offset distances increase, the volumes in the breast tend to decrease. This decrease in breast volumes reduces the minimum (7.3–13.1%), maximum (−0.8–2.2%), and mean doses (−0.2–4.4%) in the offset-CT plans.
Discussion
The use of these sophisticated positioning devices has intensified the need for CT simulators to operate with bore sizes larger than those of diagnostic CT scanners. It is critical to reconstruct accurately any portion within a geometrical CT bore range for high precision radiation treatment planning. The reconstructed CT images by eFOV were found to contain significant artifacts and image distortion, resulting in unwanted dose discrepancies during dose planning, with wide bore CT simulators (80–82 cm diameter) in many previous studies12,13. Wu et al. showed similar results regarding the dosimetric impacts of image artifacts from wide-bore CT, but there were some differences between 3D conformal radiation therapy (3D-CRT) and VMAT plans. Those results from Wu et al. demonstrated that the percent differences (within −0.5%) of VMAT were less than those (2–3%) of 3D-CRT plans10. In contrast, our study showed that VMAT plans had greater differences (0.5–6.3%) that those of 3D FIF plans (0.9–4.9%) based on breast phantom, which is very close to extended FOV limited range. These findings indicate that the dose deformation could be more influenced by the deformed breast PTV than by the HU variation. Recently, Cheung et al. also evaluated the impact of eFOV on the CT values and dosimetric accuracy by comparing the inserted CT calibration phantom with the mixed body phantoms. Their phantom studies demonstrated relatively small differences (less than 50 HU) for the inserts, except for the adipose, breast, and dense bone insets. Although the CT value accuracy should be within ± 20 HU from the manufacturer suggested values according to the IAEA guidelines14, many researchers showed that the HU variation of approximately ± 100 HU did not induce significant dose differences15. The artifacts caused by the inevitably biased position, such as in the breast, would result in severe volume changes that could significantly affect dose optimization, particularly in VMAT plans for breast cancer. These studies show that complex treatment techniques, such as intensity-modulated radiation therapy (IMRT), could be very vulnerable to shape and volume changes of the structures involving the dose-volume optimization14,15. This study has some limitations. The major limitation is that we used a home-made incomplete phantom for the planning study. For a more practical evaluation of the dose deformation due to HU variation and volume deformation close to the eFOV boundary, actual dose measurement using dedicated anthropomorphic dosimetric phantoms would be more beneficial for accurate analysis. In addition, while we found dose discrepancies around the reconstructed images close to the eFOV boundary even in the LB-CT simulator, establishing clinical guidance to avoid these errors could be difficult due to the influence of the patients’ physical conditions or the presence of additional immobilization devices in various clinical situations.
We investigated the use of the LB-CT, which has a 90-cm diameter, to evaluate the effect of the 85 cm-eFOV mode using MUSCOT algorithm with home-made inhomogeneity and humanoid breast phantoms. While image distortions did not equate to significant volume and HU variations within offset-25 cm, this offset still could be attributed to some dose distribution discrepancies in this study. The optimization of dose fluence based on segmented volumes on registered CT images for IMRT treatment could render these planning techniques more vulnerable to variations created by eFOV-based reconstructions. These variations indicate that precaution is needed when the dose calculation involves portions in the predicted eFOV area, particularly in IMRT techniques, such as VMAT and FIF. Because IMRT techniques optimize dose fluence based on segmented volumes on registered CT images, it would be more crucial to use CT images by eFOV-based reconstructions. We recommend that it should be evaluated whether the applied CT would have an appropriate range of eFOV prior to clinical application for radiation treatment planning in patients.
Methods
The LB CT (Aquilion LB, Toshiba, Japan) used in the experiment has a 70-cm sFOV, 85-cm eFOV, and 90-cm bore size. The scanning was performed at 120 kVp, 40 mA, 0.5 second rotation time, 300-mm scan range, and the CT images were reconstructed by 85-cm eFOV mode using multislice cone-beam tomography algorithm (MUSCOT)16. The inhomogeneous phantom array and a humanoid breast phantom were used in the experiments (Fig. 6).
Two home-made phantoms for CT image acquisitions
The inhomogeneous phantom array was composed of air, water, and high-density materials using 190-cm syringes (Bayer, Germany) on an acrylic plate. The high-density material was made by injection with a mixture of contrast agent (Ultravist, Bayer, Germany) and distilled water for imitating bonny density in the contrast syringes. The concentration of the injected high-density material was set to the HU value (800–1020) shown on the CT image (Fig. 7a). The home-made human phantom, which was to mimic breast and surrounding tissues, was fabricated to simulate the radiation treatment planning of breast cancer. In order to adapt the similar HU values to generic patient CT images (lung, −740 to −870; heart, 30 to 65; breast, 6 to 40) in the humanoid phantom, the lung material was chosen to cork, while the heart and breast portions were composed of thermoplastic material (Fig. 7b).
The CT scans were performed using two phantoms with different center positions, which were at the geometrical center of the bore, and various centers shifting: Off-20 cm and Off-30 cm for the inhomogeneity array phantom and Off-15 cm, Off-20 cm, Off-25 cm, and Off-30 cm for the breast phantom (Fig. 8).
Treatment planning in the breast phantom
The acquired CT images were transferred to the treatment planning system (Eclipse v.13.6, Varian, USA) for image registration. Based on the registered images, after organs at risk and planning target volume (PTV) delineations, the segmented volumes and HU line profiles were measured on the air, water, and high-density material of the inhomogeneous phantom array and humanoid breast phantom. The breast PTV was delineated on the breast phantom with a minus margin of 5 mm for dose-volume optimization (DVO) of VMAT. The dose prescription was 4256 cGy/16 fractions for both the FIF and VMAT plans. The dose calculation algorithm used for the FIF and VMAT planning was analytical anisotropic algorithm (AAA).
The relative priority settings for PTV, heart, and left lung for DVO were 990, 300, and 500, respectively. The applied parameters were same for the VMAT plans of the different offset CT images.
The FIF and VMAT treatment plans were made to investigate the impact of HU variation and image distortion in the vicinity of the biased range of the reconstructed CT images through eFOV on the treatment planning results. The FIF treatment plans were made using two main tangential oblique fields with several subfields to improve dose homogeneity. The VMAT plans were also performed along the two tangentially 180° (gantry angle from 330° to 150°) clock-wise rotating volumetric modulated arc fields to cover the whole range of eFOV based on the reconstructed portion of the breast area (Fig. 9).
The dose line profile, maximum dose, and gamma index (GI) were used for analysis of the FIF plans17. The dose difference of 2% and distance to agreement of 2% were determined for GI calculations for analyzing the dose distributions in the FIF and VMAT plans. The commercial software for dosimetry (OmniPro-I’mrt v.1.7.0021, IBA, Germany) was used for the quantitative analysis. The DDVHs were also compared to evaluate the effect of HU variation and volume deformation of the reconstructed images on the dose distribution in breast PTV.
Ethics approval and informed consent
Ethics approval and formal consent were not required as this was a phantom study.
Data availability
The datasets used and analyzed during the current study are available upon reasonable request. Please contact the authors for data requests.
References
Seco, J. & Evans, P. M. Assessing the effect of electron density in photon dose calculations. Med. Phys. 33, 540–552 (2006).
Skrzynski, W. et al. Computed tomography as a source of electron density information for radiation treatment planning. Strahlenther Onkol. 186, 327–333 (2010).
Kilby, W., Sage, J. & Rabett, V. Tolerance levels for quality assurance of electron density values generated from CT in radiotherapy treatment planning. Phys Med. Biol. 47, 1485–1492 (2002).
Garcia-Ramirez, J. L., Mutic, S., Dempsey, J. F., Low, D. A. & Purdy, J. A. Performance evaluation of an 85-cm-bore X-ray computed tomography scanner designed for radiation oncology and comparison with current diagnostic CT scanners. Int J Radiat Oncol Biol Phys. 52, 1123–1131 (2002).
Mutic, S. et al. Quality assurance for computed-tomography simulators and the computed-tomography-simulation process: report of the AAPM Radiation Therapy Committee Task Group No. 66. Med. Phys. 30, 2762–2792 (2003).
The ImPACT CT scanner evaluation group. Comparative Specifications-Wide bore CT scanners. Centre for Evidence-based Purchasing (CEP) 08029. NHS PASA (2009).
Ebert, M. A., Kenny, J. & Greer, P. B. Experience converting an RT department to full CT simulation: technical issues identified during commissioning of a wide-bore scanner. J Med. Imaging Radiat Oncol. 53, 325–330 (2009).
Hsieh, J. et al. A novel reconstruction algorithm to extend the CT scan field-of-view. Med. Phys. 31, 2385–2391 (2004).
Beeksma, B., Truant, D., Holloway, L. & Arumugam, S. An assessment of image distortion and CT number accuracy within a wide-bore CT extended field of view. Australas Phys Eng Sci Med. 38, 255–261 (2015).
Wu, V., Podgorsak, M. B., Tran, T. A., Malhotra, H. K. & Wang, I. Z. Dosimetric impact of image artifact from a wide-bore CT scanner in radiotherapy treatment planning. Med. Phys. 38, 4451–4463 (2011).
Guan, H., Yin, F. F. & Kim, J. H. Accuracy of inhomogeneity correction in photon radiotherapy from CT scans with different settings. Phys Med. Biol. 47, N223–231 (2002).
Ohnesorge, B., Flohr, T., Schwarz, K., Heiken, J. P. & Bae, K. T. Efficient correction for CT image artifacts caused by objects extending outside the scan field of view. Med. Phys. 27, 39–46 (2000).
Zamyatin, A. A. & Nakanishi, S. Extension of the reconstruction field of view and truncation correction using sinogram decomposition. Med. Phys. 34, 1593–1604 (2007).
McLean ID et al. IAEA Human Health Series No. 19, Quality assurance programme for computed tomography: diagnostic and therapy applications (2012).
Cheung, J. P., Shugard, E., Mistry, N., Pouliot, J. & Chen, J. Evaluating the impact of extended field-of-view CT reconstructions on CT values and dosimetric accuracy for radiation therapy. Med. Phys. 46, 892–901 (2019).
Seeram, E. Computed Tomography: Physical Principles and Recent Technical Advances. J Med. Imaging Radiat Sci. 41, 87–109 (2010).
Chen, M., Lu, W., Chen, Q., Ruchala, K. & Olivera, G. Efficient gamma index calculation using fast Euclidean distance transform. Phys Med. Biol. 54, 2037–2047 (2009).
Acknowledgements
We are very grateful to all collaborative Radiation Oncology department members for their assistance in this study.
Author information
Authors and Affiliations
Contributions
Y.B. and J.L. contributed equally to this work. J.L. and S.H. conceived the idea of this study. J.L. designed the study, and Y.B. performed the experiments and analyzed the data with J.L. All authors contributed to organizing this manuscript. S.H. and J.L. reviewed the final version of the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Bae, YK., Lee, JW. & Hong, S. Effects of image distortion and Hounsfield unit variations on radiation treatment plans: An extended field-of-view reconstruction in a large bore CT scanner. Sci Rep 10, 473 (2020). https://doi.org/10.1038/s41598-020-57422-y
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-020-57422-y
This article is cited by
-
Validation of deep learning-based CT image reconstruction for treatment planning
Scientific Reports (2023)
-
Assessing radiomics feature stability with simulated CT acquisitions
Scientific Reports (2022)
-
One step acid activation of bentonite derived adsorbent for the effective remediation of the new generation of industrial pesticides
Scientific Reports (2020)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.