검색
검색 팝업 닫기

Ex) Article Title, Author, Keywords

Article

Split Viewer

Original Article

Progress in Medical Physics 2021; 32(3): 59-69

Published online September 30, 2021

https://doi.org/10.14316/pmp.2021.32.3.59

Copyright © Korean Society of Medical Physics.

The Crucial Role of the Establishment of Computed Tomography Density Conversion Tables for Treating Brain or Head/Neck Tumors

Shu-Chin Yang1,2 , Su-Hua Lo3 , Li-Tsuen Shie1 , Sung-Wei Lee1 , Sheng-Yow Ho1,4,5

1Department of Radiation Oncology, Chi Mei Medical Center, Liouying, Tainan, 2Department of Medical Imaging and Radiological Sciences, Kaohsiung Medical University, Kaohsiung, 3Chiu Ho Medical System Co., LTD, Taipei, Taiwan, 4Graduate Institute of Medical Science, Chang Jung Christian University, 5Department of Radiation Oncology, Chi Mei Medical Center, Tainan, Taiwan

Correspondence to:Shu-Chin Yang (fr892072@gmail.com)
Tel: 886‑6‑6226999 (ext. 71025)
Fax: 886‑6‑6220226

Received: May 17, 2021; Revised: July 28, 2021; Accepted: August 6, 2021

Purpose: The relationship between computed tomography (CT) number and electron density (ED) has been investigated in previous studies. However, the role of these measures for guiding cancer treatment remains unclear.
Methods: The CT number was plotted against ED for different imaging protocols. The CT number was imported into ED tables for the Pinnacle treatment planning system (TPS) and was used to determine the effect on dose calculations. Conversion tables for radiation dose calculations were generated and subsequently monitored using a dosimeter to determine the effect of different CT scanning protocols and treatment sites. These tables were used to retrospectively recalculate the radiation therapy plans for 41 patients after an incorrect scanning protocol was inadvertently used. The gamma index was further used to assess the dose distribution, percentage dose difference (DD), and distance-to-agreement (DTA).
Results: For densities <1.1 g/cm3, the standard deviation of the CT number was ±0.6% and the greatest variation was noted for brain protocol conditions. For densities >1.1 g/cm3, the standard deviation of the CT number was ±21.2% and the greatest variation occurred for the tube voltage and head and neck (H&N) protocol conditions. These findings suggest that the factors most affecting the CT number are the tube voltage and treatment site (brain and H&N). Gamma index analyses for the 41 retrospective clinical cases, as well as brain metastases and H&N tumors, showed gamma passing rates >90% and <90% for the passing criterion of 2%/2 and 1%/1 mm, respectively.
Conclusions: The CT protocol should be carefully decided for TPS. The correct protocol should be used for the corresponding TPS based on the treatment site because this especially affects the dose distribution for brain metastases and H&N tumor recognition. Such steps could help reduce systematic errors.

KeywordsComputed tomography number, Electron density tables, Dose distribution, Gamma index, Imaging protocols

Computed tomography (CT) is used to determine the extent of a tumor and, subsequently, the treatment planning system (TPS). The CT number depends on various factors such as the tube voltage, field of view (FOV), scanner type, tube current, slice thickness, processing filters, and reconstruction algorithm including artifact reduction [1]. For TPS, tissue density measured using CT is converted to an electron density (ED) to calculate the dose [2]. Thus, the CT number must be accurate because it affects the dose calculation [3,4].

In most radiotherapy departments, the quality of CT images is ensured by daily, monthly, and yearly quality assurance mechanisms. The conversion of the CT number to ED is important because different CT scanning protocols use a unique conversion table to establish ED. The relationship between the CT number and ED has been evaluated in several studies [5]. Many radiotherapy departments use more than one scanning protocol to achieve better image quality assurance. The protocol used depends on the treatment site and the patient’s body characteristics. However, regardless of the protocol, most departments use only one conversion table. Thus, the dose calculations may be inaccurate [6].

The relationship between the CT number and ED has been the subject of several studies [5,7-11]. Previous studies have mainly used phantoms to illustrate the relationship between the CT number and ED. There is a concern that varying scanning parameters will change the CT number of the images and subsequently result in inaccurate dosimetric information in the TPS. In this study, the CT number is plotted against ED for different imaging protocols. All CT scanning protocols in clinical use in the departments of the present study were reviewed to create a specific conversion table for each scanning protocol. In addition, the conversion tables in 41 patients were retrospectively used and the effects on dose calculations were determined.

CT scanning was performed using a GE DiscoveryTM CT 590RT scanner (GE Medical Systems, Milwaukee, WI, USA) with a gantry aperture and a standard FOV of 80 and 65 cm, respectively. More than one scanning protocol was used to achieve better image quality, depending on the treatment site and patient’s body characteristics.

A cylindrical Virtual WaterTM phantom (Cheese phantom) was used for the different scanning protocols to plot the CT number–ED calibration curve. The phantom contains interchangeable plugs made of various tissue-equivalent materials. The device uses 12 plugs of different densities, each with a diameter and height of 2.7 and 7 cm, respectively (Fig. 1). The physical density (in gram per cubic centimeter) ranges from 0.3 (LN-300 lung) to 1.823 (T18336 cortical bone) for an ED relative to water of 1.107. The phantom was placed in the center of the CT gantry and carefully aligned using lasers. Scanning was performed using different imaging protocols and tube voltages (80–140 kVp), tube currents (300–400 mA), slice thicknesses (5 mm), and FOV (30, 50, and 65 cm). A 512×512 matrix was used to reconstruct the images. The CT scanning protocols for four treatment sites are presented in Table 1.

Table 1 Analyzed results of the computed tomography scanning parameters

Scan protocolParameter

Scan techniqueTube voltage (kVp)Tube current (mA)Slice thickness (mm)FOV (cm)
BrainAxial80, 120, and 140300530, 50, and 65
Head and neckHelical80, 120, and 140350530, 50, and 65
ChestHelical80, 120, and 140400530, 50, and 65
PelvicHelical80, 120, and 140400530, 50, and 65

FOV, field of view.


Figure 1.The cylindrical Cheese Virtual WaterTM phantom.

After image reconstruction, a circular region of interest (ROI) with a diameter of 1.5 cm was placed on each density plug and the mean CT number for each ROI was recorded. To minimize the effect of image artifacts and beam hardening, multiple CT scans of the phantom were acquired using different combinations of insertion positions and the resultant mean CT numbers were averages. The resultant CT number–ED conversion tables were compared for different tube voltages, tube currents, slice thicknesses, and FOVs as aforementioned.

The CT number–ED conversion tables were imported into Pinnacle TPS (Philips Radiation Oncology Systems, Fitchburg, WI, USA), which transforms the CT number as follows [12]:

CTHU1.024,TPSHU=0;1.024CTHU4.000;TPSHU=CTHU+1.000

The CT number–ED conversion tables for various tube voltages (80, 120, and 140 kVp) are shown in Fig. 2.

Figure 2.The computed tomography (CT) number–electron density conversion tables were split into two ranges: <1.1 and >1.1 g/cm3 densities. A straight-line fit was constructed for each range. The mean and standard deviation for each range and set of scanning conditions were used as a measure of the CT number relative to the CT scanning parameters.

A total of 48 CT number–ED conversion tables were produced and the relevant parameters were compared. The obtained curves were divided into two ranges based on Zurl et al. [11]: <1.1 and >1.1 g/cm3 of densities. A straight-line fit was constructed for each range. The mean and standard deviation for each range and set of scanning conditions were used as a measure of the CT number relative to the CT scanning parameters.

Then, the resulting new conversion tables were used to recalculate the dose distributions from the data of 41 patients, of which 9 exhibited brain metastases, 9 had head and neck cancer (H&N), 11 had lung cancer, and 12 had prostate cancer. Moreover, 123 plans were studied, including 27, 27, 33, and 36 for brain metastases, H&N tumors, lung tumors, and prostate tumors, respectively. The effect of the CT number on the dose was determined in terms of tube voltage and tumor position. These patients were retrospectively identified for inclusion in the analyses. The Institutional Review Board of the Chi Mei Medical Center, Liouying approved of the study and the approval number was 10811-L03. Informed consent was not required because this study had the lowest risk, and any possible risk to the patient did not exceed compared with those who did not participate in the study.

Of the patients, 12 and 29 received intensity-modulated radiation therapy (IMRT) and volumetric modulated arc therapy (VMAT), respectively. Each plan used a different CT number–ED conversion table for the specific tube voltage (80, 120, and 140 kVp), and the dose distribution was corrected for inhomogeneity. The remaining parameters, such as beam arrangements and monitor unit (MU), were the same.

The gamma index was used to determine the effect on dose calculations for the TPS and was measured using dose distribution, percentage dose difference (DD), and distance-to-agreement (DTA). The distribution was displayed to provide a quantitative assessment of the quality of the calculation for the regions that passed and failed the acceptance criteria [13]. The gamma index produces a passing criterion assigned to each point, for which gamma ≤1 indicates pass and other values indicate failure. The percentage of passing points in the gamma distribution is referred to as passing criteria (or %GP). Each measured point was evaluated to determine if the percentage DD or DTA exceeded the established passing criterion (i.e., 2% and 2 mm, respectively) [14-16].

The difference in dose calculations for the gamma index (dose distributions, percentage DD, and DTA) for various ED tables was calculated. To compare the four treatment sites (brain metastases, H&N cancer, lung cancer, and prostate cancer), the MU calculated for the 140 kVp CT scan for optimum coverage of the planning target volume was used to calculate other CT number–ED conversion tables that used different tube voltage settings. The planning parameters, such as the beam arrangements, fields, and MU, were the same.

Each initial plan with a calibration table was used to generate a verification plan in PTW–VeriSoft® to generate the reference two-dimensional dose distribution. Local and global gamma index analyses were used. The local gamma index analysis normalizes the percent difference for every point to the expected dose at each point, whereas the global gamma index analysis normalizes the percent difference for every point to a globally used single value, which is usually the maximum planned dose [17]. Thus, the gamma passing rate calculated using the global gamma index is always higher than or equal to the local gamma index, where the same criteria and a lower dose threshold are used. Both the global and local gamma passing rates with various gamma criteria were 2%/2 mm, 1%/1 mm, and 0.5%/0.5 mm. These rates were calculated for the absolute doses using PTW–VeriSoft®, which linearly interpolates the calculated dose to a resolution <0.5 mm.

1. Computed tomography number-electron density conversion

Calculations using the Cheese phantom showed different effects for image protocols for 12 different density ranges. For densities <1.1 g/cm3, the standard deviation for the CT number was ±0.6% and the greatest variation was noted for brain protocol conditions. For densities >1.1 g/cm3, the standard deviation for the CT number was ±21.2% and the value varied most significantly with tube voltage and H&N protocol conditions (Fig. 3).

Figure 3.The computed tomography (CT) number–electron density conversion tables using different tube voltages. For densities >1.1 g/cm3, the standard deviation for the CT number was ±21.2%. For densities <1.1 g/cm3, the standard deviation for the CT number was ±0.6%.

This is probably due to the artifact in the plug that has different chemical compositions but the same physical density. The differences in CT numbers versus tube voltages are minimal in the density region from 0.3 (lung) to 1.0 (water), and no differences were noted in low-density materials (lung) of 4% (332 HU at 80 kVp vs. 320 HU at 140 kVp). Significant differences in high-density materials (cortical bone) were shown in 28% (2,856 HU at 80 kVp vs. 2,223 HU at 140 kVp) and a higher tube voltage that resulted in a lower CT number.

The CT number to relative ED calibration curves for the differences in FOV (30, 50, and 65 FOV) was compared. The differences in CT numbers versus FOV are maximum in the density region from 0.3 (lung) to 1.0 (water); the significant differences in CT number were 17% (371 HU at 30 FOV vs. 317 HU at 65 FOV) and 9% (2,332 HU at 30 FOV vs. 2,118 at 65 FOV) occurring for low-density materials (lung) and high-density materials (cortical bone), respectively.

Moreover, the CT number–ED tables showed that there is no significant relationship between tube current and slice thickness.

2. Global gamma passing rates for clinical cases

Fig. 4 shows the global gamma passing rates for various passing criteria for ED calibrations using 80 and 140 kVp and 120 and 140 kVp for various treatment sites as measured using a dose distribution dosimeter with PTW–VeriSoft® measurements.

Figure 4.Local gamma passing rates with various passing criteria calculated from the PTW–VeriSoft® dosimeter measurements are shown for (a) the 80 and 140 kVp plans and (b) the 120 and 140 kVp plans. Global gamma passing rates with various passing criteria calculated from the PTW–VeriSoft® dosimeter measurements are shown for (c) the 80 and 140 kVp plans and (d) the 120 s and 140 kVp plans. The dosimetric effects for the four treatment sites (brain metastases plans, head and neck cancer plans, lung cancer plans, and prostate cancer plans) were analyzed. The error bars above the histogram indicate the standard deviation.

The comparison of ED calibrations for 80 and 140 kVp for the four tumor sites showed the lowest gamma passing rates. For the passing criterion of 1%/1 mm, the mean gamma passing rates for the four treatment sites were 87.62%, 96.30%, 100%, and 99.62%, respectively (P<0.001). For the passing criterion of 2%/2 mm, the mean gamma passing rates for the treatment sites were 98.43%, 98.23%, 100%, and 99.85%, respectively (P=0.352). For 120 and 140 kVp using the passing criterion of 1%/1 mm, the mean gamma passing rates for the four treatment sites were 98.29%, 99.37%, 100%, and 100%, respectively (P=0.025). For the passing criterion of 2%/2 mm, the mean gamma passing rates for the four treatment sites were 99.88%, 99.79%, 100%, and 100%, respectively (P=0.649).

The global gamma passing rates for various passing criteria for ED calibrations with different treatment techniques (VMAT and IMAT) were compared, as measured using a dose distribution dosimeter with PTW–VeriSoft® measurements. For the passing criterion of 1%/1 mm, the mean gamma passing rates for the VMAT and IMRT were 97.25% and 98.70%, respectively (P=0.15). For the passing criterion of 2%/2 mm, the mean gamma passing rates for the VMAT and IMRT were 99.66% and 99.94%, respectively (P=0.21).

3. Local gamma passing rates for clinical cases

Fig. 4 shows the local gamma passing rate for various passing criteria for ED calibrations of 80 and 140 kVp and 120 and 140 kVp for various treatment sites measured using a dose distribution dosimeter with PTW–VeriSoft® measurements.

The ED calibrations for 80 and 140 kVp for brain and H&N plans were compared to determine the lowest gamma passing rate using a dose distribution dosimeter. For the passing criterion of 1%/1 mm, the mean gamma passing rates for the four treatment sites were 83.53%, 90.07%, 97.13%, and 99.00%, respectively (P=0.001). For the passing criterion of 2%/2 mm, the mean gamma passing rates for the four treatment sites were 98.11%, 98.49%, 98.75%, and 99.8%, respectively (P=0.248). Comparison of the results for 120 and 140 kVp showed that the mean gamma passing rates for the four tumor sites were 96.75%, 98.91%, 100%, and 100%, respectively, (P=0.005) at the passing criterion of 1%/1 mm. The mean gamma passing rates for the four treatment sites were 99.53%, 99.37%, 100%, and 100%, respectively, (P=0.452) for the passing criterion of 2%/2 mm.

This study determined that the factors most affecting the CT number were the tube voltage and treatment site (brain and H&N). Thus, the tube voltage (80–140 kVp) and treatment site (brain and H&N) were used to correct the dose distribution for inhomogeneity.

The sagittal view of the brain metastases treatment plan is parallel to the inferior orbitomeatal (IOM) line and across the frontal lobe (the top of the ROI), across the corpus callosum (the middle of the ROI), and across the cerebellum (the end of the ROI; Fig. 5). The sagittal view of the H&N tumor treatment plan system is parallel to the IOM line and across the maxillary sinus (the top of the ROI), across the mandible (the middle of the ROI), and across the C7 spinal disc (the end of the ROI). The gamma index was used to compare the dose distribution (Fig. 6).

Figure 5.Sagittal view of a brain metastases treatment plan (a) parallel to the IOM line and across the frontal lobe (the top of the ROI), (b) parallel to the IOM line and across the corpus callosum (the middle of the ROI), and (c) parallel to the IOM line and across the cerebellum (the end of the ROI). The local gamma passing rates for various passing criteria were compared for 1%/1 mm and show (d–f) ED calibrations of 80 and 140 kVp and (g–i) ED calibrations of 120 and 140 kVp. Passing criteria: gamma ≤1. Green 90.0%–100.0%, yellow 75.0%–90.0%, and red 0.0%–75.0%. IOM, inferior orbitomeatal; ROI, circular region of interest; ED, electron density.

Figure 6.Sagittal view of H&N tumor treatment plan system (a) parallel to the line and across the maxillary sinus (the top of the ROI), (b) parallel to the IOM line and across the mandible (the middle of the ROI), and (c) parallel to the IOM line and across the C7 (the end of the ROI). The local gamma passing rates for various passing criteria were compared for 1%/1 mm and show (d–f) ED calibrations of 80 and 120 kVp and (g–i) ED calibrations of 120 and 140 kVp. Passing criteria: gamma ≤1. Green 90.0%–100.0%, yellow 75.0%–90.0%, and red 0.0%–75.0%. IOM, inferior orbitomeatal; ROI, circular region of interest; ED, electron density.

The local gamma passing rates for various passing criteria were compared for ED calibrations of 80 kVp and 140 kVp and 120 and 140 kVp for various brain metastases measured using a dose distribution dosimeter. The results are shown in Table 2. The measurement across the frontal lobe dosimeter for the passing criterion of 1%/1 mm showed mean gamma passing rates of 74.48±15.39 and 95.64±5.95, respectively (P<0.001). The measurement across the corpus callosum dosimeter for the passing criterion of 1%/1 mm showed mean gamma passing rates of 96.35±2.67 and 98.76±3.03, respectively (P=0.024). The measurement across the cerebellum dosimeter for the passing criterion of 1%/1 mm showed mean gamma passing rates of 79.32±6.48 and 95.50±3.77, respectively (P<0.001).

Table 2 The local gamma passing rates of the various passing criteria for various brain metastases measured

Treatment sitesPassing criterionLocal gamma passing rates (%)P-value

80 and 140 kVp120 and 140 kVp
Brain tumor
Across the frontal lobe2%/2 mm96.67±4.9699.53±0.670.009
1%/1 mm74.48±15.3995.64±5.95<0.001
0.5%/0.5 mm42.14±24.3664.31±21.420.007
Across the corpus callosum2%/2 mm99.70±0.4699.90±0.160.113
1%/1 mm96.35±2.6798.76±3.030.024
0.5%/0.5 mm81.65±7.4289.35±9.150.014
Across the cerebellum2%/2 mm98.14±1.7399.90±0.130.001
1%/1 mm79.32±6.4895.50±3.77<0.001
0.5%/0.5 mm39.40±13.5162.78±8.150.005

Values are presented as mean±standard deviation.


The local gamma passing rates for various passing criteria for ED calibrations of 80 and 140 kVp and 120 and 140 kVp for various H&N tumors measured using a dose distribution dosimeter are shown in Table 3. The measured value across the maxillary sinus dosimeter for the passing criterion of 1%/1 mm showed mean gamma passing rates of 88.69±13.31 and 98.96±1.56, respectively (P=0.019). The measured value across the mandible dosimeter for the passing criterion of 1%/1 mm showed mean gamma passing rates of 89.15±9.41 and 97.94±2.06, respectively (P=0.001). The measured value crossing the C7 spinal segment dosimeter for the passing criterion of 1%/1 mm showed mean gamma passing rates of 92.21±13.38 and 99.85±0.41, respectively (P=0.030).

Table 3 The local gamma passing rates of the various passing criteria for various head and neck tumors measured

Treatment sitePassing criterionLocal gamma passing rates (%)P-value

80 and 140 kVp120 and 140 kVp
H&N tumor
Across the maxillary sinus2%/2 mm98.50±2.2199.94±0.780.022
1%/1 mm88.69±13.3198.96±1.560.019
0.5%/0.5 mm65.60±18.2890.63±9.35<0.001
Across the mandible2%/2 mm97.59±2.4899.38±1.400.018
1%/1 mm89.15±9.4197.94±2.060.001
0.5%/0.5 mm69.85±17.6989.05±7.48<0.001
Across the C7 spinal segment2%/2 mm99.39±1.01100.00±0.000.022
1%/1 mm92.21±13.3899.85±0.410.030
0.5%/0.5 mm80.98±23.8795.50±9.220.031

Values are presented as mean±standard deviation.


The local gamma passing rates for various passing criteria for ED calibrations with different treatment techniques (VMAT and IMAT) were compared as measured using a dose distribution dosimeter with PTW–VeriSoft® measurements. For the passing criterion of 1%/1 mm, the mean gamma passing rates for the VMAT and IMRT were 96.03% and 97.24%, respectively (P=0.40). For the passing criterion of 2%/2 mm, the mean gamma passing rates for the VMAT and IMRT were 99.37% and 99.56%, respectively (P=0.51).

The CT number depends on the CT scanner image protocol conditions. In a previous study, the specifications and standards for a radiotherapy department at a regional hospital in southern Taiwan were assessed [18]. The study found that the smallest deviations in the relationship between the CT number and ED occurred for low ED materials and increased with increasing density. Thomas [19] noted that this deviation was <1% for materials with a low ED for different calibration techniques and manufacturers.

The CT number in this study was shown to be stable concerning the different acquisition parameters except for low tube voltage and materials with a higher ED, which resulted in deviations of about ±21.2%.

The results of a similar study assessing the relationship between the CT number and relative ED for various tube voltages reported a maximum difference in the mean dose of 1.1% [20]. Zurl et al. [11] used different image protocol conditions to assess variation in the HU of 20%, which resulted in an average dose error of approximately 1.5%. Moreover, Van Dyk et al. [21] and Gershkevitsh et al. [22] showed that a low tube voltage and an increased thickness of the anatomical structure being scanned resulted in a more accurate CT density conversion table that avoided calculation errors >2%. Kilby et al. [23] showed that radiotherapy at a reasonable dose level was possible for different CT imaging protocol conditions using an average CT density conversion table.

Taken together, these studies support that the factors most affecting the CT number are the tube voltage and treatment site. Given these findings, the gamma index analysis was used to compare the dose distribution in real clinical cases. Results showed that the gamma passing rates for ED calibrations of 80, 120, and 140 kVp with the passing criterion of 1%/1 mm were significant for the four treatment sites. Venselaar et al. [24] showed that the passing criterion of 2% or 2 mm in the TPS dose calculations in the central beam was acceptable. This study also noted a maximum difference of only gamma passing rate <90% for the passing criterion of 1%/1 mm for ED calibrations of 80 and 140 kVp at the treatment site (brain and H&N). Thus, the results were clinically insignificant.

du Plessis et al. [9] reported an average dose error of about 20% to 70% in the maxillary sinus region for a head model for different calibration techniques and manufacturers. The anatomical effects of homogeneities must be considered because the brain and H&N tumor sites are more likely to exhibit an uneven dose distribution. Thus, using the correct CT number–ED conversion table is essential.

Gamma passing rates were used to compare 80, 120, and 140 kVp for brain metastases and H&N tumors using the passing criteria of 1%/1 and 2%/2 mm. The results for both regions were statistically significant. The gamma passing rates were >90% and <90% for the passing criteria of 2%/2 and 1%/1 mm, respectively. Zurl et al. [11] used initial calculations for inhomogeneous features of tissue to show that the bone produces more statistically significant results than soft tissue. The skull almost entirely covers the brain tissue for the brain and H&N tumors at lower tube voltages. Thus, X-ray absorption is important when the X-ray beams traverse inhomogeneous areas. While there are very small differences for the passing criterion of 1%/1 mm for other treatment sites, the difference is greater for brain and H&N tumors.

In general, the CT number is most significant for high-density materials. For CT scanning parameter settings, adjusting the tube voltage has the most significant effect on the brain and H&N tumor imaging. Yang et al. [25] developed the Swiss Cheese model in 1990 to show that radiological errors can be caused by a series of mistakes. When a medical radiological technologist performs a CT simulation, the operating procedure checklist must be executed to ensure consistency in scanning parameters and to ensure that the correct tube voltage is used to prevent deviation in the TPS dose distribution [25].

The use of different CT scanning protocols resulted in a variation of up to 21.2% in the HU values. The local gamma passing rates for various passing criteria for ED calibrations of 80 and 140 kVp were calculated using real patient data, which revealed that brain metastases and H&N tumors showed a gamma passing rate of <90% for the passing criterion of 1%/1 mm. Regarding standardized scanning parameter settings, a specific conversion table needs to be created if the tube voltage in the CT scanner is adjusted. The correct protocol should be used for the corresponding TPS based on the treatment site, especially because this affects the dose distribution for brain metastases and H&N tumor recognition. Such steps could help reduce systematic errors.

The authors are greatly thankful to Mrs. Lo of Chiu Ho Medical System Co., LTD, Taiwan, for the knowledge on Pinnacle TPS. This study was supported by the Chi Mei Medical Center, Liouying, from January 2020 (grants CLFHR10914).


The data that support the findings of this study are available on request from the corresponding author.


Conceptualization: Shu-Chin Yang, Su-Hua Lo, Li-Tsuen Shie, Sung-Wei Lee, and Sheng-Yow Ho. Data curation: Shu-Chin Yang, Su-Hua Lo, and Li-Tsuen Shie. Formal analysis: Shu-Chin Yang and Su-Hua Lo. Funding acquisition: Shu-Chin Yang. Investigation: Shu-Chin Yang and Su-Hua Lo. Methodology: Shu-Chin Yang, Su-Hua Lo, Li-Tsuen Shie, Sung-Wei Lee, and Sheng-Yow Ho. Project administration: Shu-Chin Yang and Su-Hua Lo. Resources: Sung-Wei Lee. Software: Shu-Chin Yang and Su-Hua Lo. Supervision: Shu-Chin Yang. Validation: Shu-Chin Yang and Su-Hua Lo. Visualization: Shu-Chin Yang, Su-Hua Lo, Li-Tsuen Shie, Sung-Wei Lee, and Sheng-Yow Ho. Writing–original draft: Shu-Chin Yang and Su-Hua Lo. Writing–review & editing: Shu-Chin Yang, Su-Hua Lo, and Li-Tsuen Shie.

  1. Das IJ, Cheng CW, Cao M, Johnstone PA. Computed tomography imaging parameters for inhomogeneity correction in radiation treatment planning. J Med Phys. 2016;41:3-11.
    Pubmed KoreaMed CrossRef
  2. Geise RA, McCullough EC. The use of CT scanners in megavoltage photon-beam therapy planning. Radiology. 1977;124:133-141.
    Pubmed CrossRef
  3. Goodsitt MM, Christodoulou EG, Larson SC. Accuracies of the synthesized monochromatic CT numbers and effective atomic numbers obtained with a rapid kVp switching dual energy CT scanner. Med Phys. 2011;38:2222-2232.
    Pubmed CrossRef
  4. Lv P, Liu J, Zhang R, Jia Y, Gao J. Combined use of automatic tube voltage selection and current modulation with iterative reconstruction for CT evaluation of small hypervascular hepatocellular carcinomas: effect on lesion conspicuity and image quality. Korean J Radiol. 2015;16:531-540.
  5. Guan H, Yin FF, Kim JH. Accuracy of inhomogeneity correction in photon radiotherapy from CT scans with different settings. Phys Med Biol. 2002;47:N223-N231.
    Pubmed CrossRef
  6. Davis AT, Palmer AL, Nisbet A. Can CT scan protocols used for radiotherapy treatment planning be adjusted to optimize image quality and patient dose? A systematic review. Br J Radiol. 2017;90:20160406.
    Pubmed KoreaMed CrossRef
  7. Cropp RJ, Seslija P, Tso D, Thakur Y. Scanner and kVp dependence of measured CT numbers in the ACR CT phantom. J Appl Clin Med Phys. 2013;14:4417.
    Pubmed KoreaMed CrossRef
  8. Skrzyński W, Zielińska-Dabrowska S, Wachowicz M, Slusarczyk-Kacprzyk W, Kukołowicz PF, Bulski W. Computed tomography as a source of electron density information for radiation treatment planning. Strahlenther Onkol. 2010;186:327-333.
    Pubmed CrossRef
  9. du Plessis FC, Willemse CA, Lötter MG, Goedhals L. Comparison of the Batho, ETAR and Monte Carlo dose calculation methods in CT based patient models. Med Phys. 2001;28:582-589.
    Pubmed CrossRef
  10. Cozzi L, Fogliata A, Buffa F, Bieri S. Dosimetric impact of computed tomography calibration on a commercial treatment planning system for external radiation therapy. Radiother Oncol. 1998;48:335-338.
    Pubmed CrossRef
  11. Zurl B, Tiefling R, Winkler P, Kindl P, Kapp KS. Hounsfield units variations: impact on CT-density based conversion tables and their effects on dose distribution. Strahlenther Onkol. 2014;190:88-93.
    Pubmed CrossRef
  12. Görlitz E. Dosimetrische und verfahrenstechnische Untersuchungen zur Qualitätssicherung eines Bestrahlungsplanungsprogramms. [Thesis]. Hambur: Universität Hamburg. 2006.
  13. Low DA, Harms WB, Mutic S, Purdy JA. A technique for the quantitative evaluation of dose distributions. Med Phys. 1998;25:656-661.
    Pubmed CrossRef
  14. Cheng A, Gerber RL, Wong JW, Purdy JA; Harms WB Sr. Systematic verification of a three-dimensional electron beam dose calculation algorithm. Med Phys. 1996;23:685-693.
    Pubmed CrossRef
  15. Shiu AS, Tung S, Hogstrom KR, Wong JW, Gerber RL, Harms WB, et al. Verification data for electron beam dose algorithms. Med Phys. 1992;19:623-636.
    Pubmed CrossRef
  16. William BH, Daniel AL, James AP. A quantitative software analysis tool for verifying 3-d dose calculation programs. Int J Radiat Oncol Biol Phys. 1994;30(Suppl 1):187.
  17. Ezzell GA, Galvin JM, Low D, Palta JR, Rosen I, Sharpe MB, et al. Guidance document on delivery, treatment planning, and clinical implementation of IMRT: report of the IMRT Subcommittee of the AAPM Radiation Therapy Committee. Med Phys. 2003;30:2089-2115.
    Pubmed CrossRef
  18. Fraass B, Doppke K, Hunt M, Kutcher G, Starkschall G, Stern R, et al. American Association of Physicists in Medicine Radiation Therapy Committee Task Group 53: quality assurance for clinical radiotherapy treatment planning. Med Phys. 1998;25:1773-1829.
    Pubmed CrossRef
  19. Thomas SJ. Relative electron density calibration of CT scanners for radiotherapy treatment planning. Br J Radiol. 1999;72:781-786.
    Pubmed CrossRef
  20. Rhee DJ, Kim S, Jeong DH, Moon YM, Kim JK. Effects of the difference in tube voltage of the CT scanner on dose calculation. J Korean Phys Soc. 2015;67:123-128.
  21. Van Dyk J, Barnett RB, Cygler JE, Shragge PC. Commissioning and quality assurance of treatment planning computers. Int J Radiat Oncol Biol Phys. 1993;26:261-273.
    Pubmed CrossRef
  22. Gershkevitsh E, Schmidt R, Velez G, Miller D, Korf E, Yip F, et al. Dosimetric verification of radiotherapy treatment planning systems: results of IAEA pilot study. Radiother Oncol. 2008;89:338-346.
    Pubmed CrossRef
  23. 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. 2002;47:1485-1492.
    Pubmed CrossRef
  24. Venselaar J, Welleweerd H, Mijnheer B. Tolerances for the accuracy of photon beam dose calculations of treatment planning systems. Radiother Oncol. 2001;60:191-201.
    Pubmed CrossRef
  25. Yang SC, Chang YL, Chang LY, Ho SY. The effectiveness of skill learning assessment using direct observation of procedural skills (DOPS) in the training of post-graduate training year and internship students of radiological technologist. C J Radiologic Tech. 2018;42:89-96.

Article

Original Article

Progress in Medical Physics 2021; 32(3): 59-69

Published online September 30, 2021 https://doi.org/10.14316/pmp.2021.32.3.59

Copyright © Korean Society of Medical Physics.

The Crucial Role of the Establishment of Computed Tomography Density Conversion Tables for Treating Brain or Head/Neck Tumors

Shu-Chin Yang1,2 , Su-Hua Lo3 , Li-Tsuen Shie1 , Sung-Wei Lee1 , Sheng-Yow Ho1,4,5

1Department of Radiation Oncology, Chi Mei Medical Center, Liouying, Tainan, 2Department of Medical Imaging and Radiological Sciences, Kaohsiung Medical University, Kaohsiung, 3Chiu Ho Medical System Co., LTD, Taipei, Taiwan, 4Graduate Institute of Medical Science, Chang Jung Christian University, 5Department of Radiation Oncology, Chi Mei Medical Center, Tainan, Taiwan

Correspondence to:Shu-Chin Yang (fr892072@gmail.com)
Tel: 886‑6‑6226999 (ext. 71025)
Fax: 886‑6‑6220226

Received: May 17, 2021; Revised: July 28, 2021; Accepted: August 6, 2021

Abstract

Purpose: The relationship between computed tomography (CT) number and electron density (ED) has been investigated in previous studies. However, the role of these measures for guiding cancer treatment remains unclear.
Methods: The CT number was plotted against ED for different imaging protocols. The CT number was imported into ED tables for the Pinnacle treatment planning system (TPS) and was used to determine the effect on dose calculations. Conversion tables for radiation dose calculations were generated and subsequently monitored using a dosimeter to determine the effect of different CT scanning protocols and treatment sites. These tables were used to retrospectively recalculate the radiation therapy plans for 41 patients after an incorrect scanning protocol was inadvertently used. The gamma index was further used to assess the dose distribution, percentage dose difference (DD), and distance-to-agreement (DTA).
Results: For densities <1.1 g/cm3, the standard deviation of the CT number was ±0.6% and the greatest variation was noted for brain protocol conditions. For densities >1.1 g/cm3, the standard deviation of the CT number was ±21.2% and the greatest variation occurred for the tube voltage and head and neck (H&N) protocol conditions. These findings suggest that the factors most affecting the CT number are the tube voltage and treatment site (brain and H&N). Gamma index analyses for the 41 retrospective clinical cases, as well as brain metastases and H&N tumors, showed gamma passing rates >90% and <90% for the passing criterion of 2%/2 and 1%/1 mm, respectively.
Conclusions: The CT protocol should be carefully decided for TPS. The correct protocol should be used for the corresponding TPS based on the treatment site because this especially affects the dose distribution for brain metastases and H&N tumor recognition. Such steps could help reduce systematic errors.

Keywords: Computed tomography number, Electron density tables, Dose distribution, Gamma index, Imaging protocols

Introduction

Computed tomography (CT) is used to determine the extent of a tumor and, subsequently, the treatment planning system (TPS). The CT number depends on various factors such as the tube voltage, field of view (FOV), scanner type, tube current, slice thickness, processing filters, and reconstruction algorithm including artifact reduction [1]. For TPS, tissue density measured using CT is converted to an electron density (ED) to calculate the dose [2]. Thus, the CT number must be accurate because it affects the dose calculation [3,4].

In most radiotherapy departments, the quality of CT images is ensured by daily, monthly, and yearly quality assurance mechanisms. The conversion of the CT number to ED is important because different CT scanning protocols use a unique conversion table to establish ED. The relationship between the CT number and ED has been evaluated in several studies [5]. Many radiotherapy departments use more than one scanning protocol to achieve better image quality assurance. The protocol used depends on the treatment site and the patient’s body characteristics. However, regardless of the protocol, most departments use only one conversion table. Thus, the dose calculations may be inaccurate [6].

The relationship between the CT number and ED has been the subject of several studies [5,7-11]. Previous studies have mainly used phantoms to illustrate the relationship between the CT number and ED. There is a concern that varying scanning parameters will change the CT number of the images and subsequently result in inaccurate dosimetric information in the TPS. In this study, the CT number is plotted against ED for different imaging protocols. All CT scanning protocols in clinical use in the departments of the present study were reviewed to create a specific conversion table for each scanning protocol. In addition, the conversion tables in 41 patients were retrospectively used and the effects on dose calculations were determined.

Materials|Methods

CT scanning was performed using a GE DiscoveryTM CT 590RT scanner (GE Medical Systems, Milwaukee, WI, USA) with a gantry aperture and a standard FOV of 80 and 65 cm, respectively. More than one scanning protocol was used to achieve better image quality, depending on the treatment site and patient’s body characteristics.

A cylindrical Virtual WaterTM phantom (Cheese phantom) was used for the different scanning protocols to plot the CT number–ED calibration curve. The phantom contains interchangeable plugs made of various tissue-equivalent materials. The device uses 12 plugs of different densities, each with a diameter and height of 2.7 and 7 cm, respectively (Fig. 1). The physical density (in gram per cubic centimeter) ranges from 0.3 (LN-300 lung) to 1.823 (T18336 cortical bone) for an ED relative to water of 1.107. The phantom was placed in the center of the CT gantry and carefully aligned using lasers. Scanning was performed using different imaging protocols and tube voltages (80–140 kVp), tube currents (300–400 mA), slice thicknesses (5 mm), and FOV (30, 50, and 65 cm). A 512×512 matrix was used to reconstruct the images. The CT scanning protocols for four treatment sites are presented in Table 1.

Table 1 . Analyzed results of the computed tomography scanning parameters.

Scan protocolParameter

Scan techniqueTube voltage (kVp)Tube current (mA)Slice thickness (mm)FOV (cm)
BrainAxial80, 120, and 140300530, 50, and 65
Head and neckHelical80, 120, and 140350530, 50, and 65
ChestHelical80, 120, and 140400530, 50, and 65
PelvicHelical80, 120, and 140400530, 50, and 65

FOV, field of view..



Figure 1. The cylindrical Cheese Virtual WaterTM phantom.

After image reconstruction, a circular region of interest (ROI) with a diameter of 1.5 cm was placed on each density plug and the mean CT number for each ROI was recorded. To minimize the effect of image artifacts and beam hardening, multiple CT scans of the phantom were acquired using different combinations of insertion positions and the resultant mean CT numbers were averages. The resultant CT number–ED conversion tables were compared for different tube voltages, tube currents, slice thicknesses, and FOVs as aforementioned.

The CT number–ED conversion tables were imported into Pinnacle TPS (Philips Radiation Oncology Systems, Fitchburg, WI, USA), which transforms the CT number as follows [12]:

CTHU1.024,TPSHU=0;1.024CTHU4.000;TPSHU=CTHU+1.000

The CT number–ED conversion tables for various tube voltages (80, 120, and 140 kVp) are shown in Fig. 2.

Figure 2. The computed tomography (CT) number–electron density conversion tables were split into two ranges: <1.1 and >1.1 g/cm3 densities. A straight-line fit was constructed for each range. The mean and standard deviation for each range and set of scanning conditions were used as a measure of the CT number relative to the CT scanning parameters.

A total of 48 CT number–ED conversion tables were produced and the relevant parameters were compared. The obtained curves were divided into two ranges based on Zurl et al. [11]: <1.1 and >1.1 g/cm3 of densities. A straight-line fit was constructed for each range. The mean and standard deviation for each range and set of scanning conditions were used as a measure of the CT number relative to the CT scanning parameters.

Then, the resulting new conversion tables were used to recalculate the dose distributions from the data of 41 patients, of which 9 exhibited brain metastases, 9 had head and neck cancer (H&N), 11 had lung cancer, and 12 had prostate cancer. Moreover, 123 plans were studied, including 27, 27, 33, and 36 for brain metastases, H&N tumors, lung tumors, and prostate tumors, respectively. The effect of the CT number on the dose was determined in terms of tube voltage and tumor position. These patients were retrospectively identified for inclusion in the analyses. The Institutional Review Board of the Chi Mei Medical Center, Liouying approved of the study and the approval number was 10811-L03. Informed consent was not required because this study had the lowest risk, and any possible risk to the patient did not exceed compared with those who did not participate in the study.

Of the patients, 12 and 29 received intensity-modulated radiation therapy (IMRT) and volumetric modulated arc therapy (VMAT), respectively. Each plan used a different CT number–ED conversion table for the specific tube voltage (80, 120, and 140 kVp), and the dose distribution was corrected for inhomogeneity. The remaining parameters, such as beam arrangements and monitor unit (MU), were the same.

The gamma index was used to determine the effect on dose calculations for the TPS and was measured using dose distribution, percentage dose difference (DD), and distance-to-agreement (DTA). The distribution was displayed to provide a quantitative assessment of the quality of the calculation for the regions that passed and failed the acceptance criteria [13]. The gamma index produces a passing criterion assigned to each point, for which gamma ≤1 indicates pass and other values indicate failure. The percentage of passing points in the gamma distribution is referred to as passing criteria (or %GP). Each measured point was evaluated to determine if the percentage DD or DTA exceeded the established passing criterion (i.e., 2% and 2 mm, respectively) [14-16].

The difference in dose calculations for the gamma index (dose distributions, percentage DD, and DTA) for various ED tables was calculated. To compare the four treatment sites (brain metastases, H&N cancer, lung cancer, and prostate cancer), the MU calculated for the 140 kVp CT scan for optimum coverage of the planning target volume was used to calculate other CT number–ED conversion tables that used different tube voltage settings. The planning parameters, such as the beam arrangements, fields, and MU, were the same.

Each initial plan with a calibration table was used to generate a verification plan in PTW–VeriSoft® to generate the reference two-dimensional dose distribution. Local and global gamma index analyses were used. The local gamma index analysis normalizes the percent difference for every point to the expected dose at each point, whereas the global gamma index analysis normalizes the percent difference for every point to a globally used single value, which is usually the maximum planned dose [17]. Thus, the gamma passing rate calculated using the global gamma index is always higher than or equal to the local gamma index, where the same criteria and a lower dose threshold are used. Both the global and local gamma passing rates with various gamma criteria were 2%/2 mm, 1%/1 mm, and 0.5%/0.5 mm. These rates were calculated for the absolute doses using PTW–VeriSoft®, which linearly interpolates the calculated dose to a resolution <0.5 mm.

Results

1. Computed tomography number-electron density conversion

Calculations using the Cheese phantom showed different effects for image protocols for 12 different density ranges. For densities <1.1 g/cm3, the standard deviation for the CT number was ±0.6% and the greatest variation was noted for brain protocol conditions. For densities >1.1 g/cm3, the standard deviation for the CT number was ±21.2% and the value varied most significantly with tube voltage and H&N protocol conditions (Fig. 3).

Figure 3. The computed tomography (CT) number–electron density conversion tables using different tube voltages. For densities >1.1 g/cm3, the standard deviation for the CT number was ±21.2%. For densities <1.1 g/cm3, the standard deviation for the CT number was ±0.6%.

This is probably due to the artifact in the plug that has different chemical compositions but the same physical density. The differences in CT numbers versus tube voltages are minimal in the density region from 0.3 (lung) to 1.0 (water), and no differences were noted in low-density materials (lung) of 4% (332 HU at 80 kVp vs. 320 HU at 140 kVp). Significant differences in high-density materials (cortical bone) were shown in 28% (2,856 HU at 80 kVp vs. 2,223 HU at 140 kVp) and a higher tube voltage that resulted in a lower CT number.

The CT number to relative ED calibration curves for the differences in FOV (30, 50, and 65 FOV) was compared. The differences in CT numbers versus FOV are maximum in the density region from 0.3 (lung) to 1.0 (water); the significant differences in CT number were 17% (371 HU at 30 FOV vs. 317 HU at 65 FOV) and 9% (2,332 HU at 30 FOV vs. 2,118 at 65 FOV) occurring for low-density materials (lung) and high-density materials (cortical bone), respectively.

Moreover, the CT number–ED tables showed that there is no significant relationship between tube current and slice thickness.

2. Global gamma passing rates for clinical cases

Fig. 4 shows the global gamma passing rates for various passing criteria for ED calibrations using 80 and 140 kVp and 120 and 140 kVp for various treatment sites as measured using a dose distribution dosimeter with PTW–VeriSoft® measurements.

Figure 4. Local gamma passing rates with various passing criteria calculated from the PTW–VeriSoft® dosimeter measurements are shown for (a) the 80 and 140 kVp plans and (b) the 120 and 140 kVp plans. Global gamma passing rates with various passing criteria calculated from the PTW–VeriSoft® dosimeter measurements are shown for (c) the 80 and 140 kVp plans and (d) the 120 s and 140 kVp plans. The dosimetric effects for the four treatment sites (brain metastases plans, head and neck cancer plans, lung cancer plans, and prostate cancer plans) were analyzed. The error bars above the histogram indicate the standard deviation.

The comparison of ED calibrations for 80 and 140 kVp for the four tumor sites showed the lowest gamma passing rates. For the passing criterion of 1%/1 mm, the mean gamma passing rates for the four treatment sites were 87.62%, 96.30%, 100%, and 99.62%, respectively (P<0.001). For the passing criterion of 2%/2 mm, the mean gamma passing rates for the treatment sites were 98.43%, 98.23%, 100%, and 99.85%, respectively (P=0.352). For 120 and 140 kVp using the passing criterion of 1%/1 mm, the mean gamma passing rates for the four treatment sites were 98.29%, 99.37%, 100%, and 100%, respectively (P=0.025). For the passing criterion of 2%/2 mm, the mean gamma passing rates for the four treatment sites were 99.88%, 99.79%, 100%, and 100%, respectively (P=0.649).

The global gamma passing rates for various passing criteria for ED calibrations with different treatment techniques (VMAT and IMAT) were compared, as measured using a dose distribution dosimeter with PTW–VeriSoft® measurements. For the passing criterion of 1%/1 mm, the mean gamma passing rates for the VMAT and IMRT were 97.25% and 98.70%, respectively (P=0.15). For the passing criterion of 2%/2 mm, the mean gamma passing rates for the VMAT and IMRT were 99.66% and 99.94%, respectively (P=0.21).

3. Local gamma passing rates for clinical cases

Fig. 4 shows the local gamma passing rate for various passing criteria for ED calibrations of 80 and 140 kVp and 120 and 140 kVp for various treatment sites measured using a dose distribution dosimeter with PTW–VeriSoft® measurements.

The ED calibrations for 80 and 140 kVp for brain and H&N plans were compared to determine the lowest gamma passing rate using a dose distribution dosimeter. For the passing criterion of 1%/1 mm, the mean gamma passing rates for the four treatment sites were 83.53%, 90.07%, 97.13%, and 99.00%, respectively (P=0.001). For the passing criterion of 2%/2 mm, the mean gamma passing rates for the four treatment sites were 98.11%, 98.49%, 98.75%, and 99.8%, respectively (P=0.248). Comparison of the results for 120 and 140 kVp showed that the mean gamma passing rates for the four tumor sites were 96.75%, 98.91%, 100%, and 100%, respectively, (P=0.005) at the passing criterion of 1%/1 mm. The mean gamma passing rates for the four treatment sites were 99.53%, 99.37%, 100%, and 100%, respectively, (P=0.452) for the passing criterion of 2%/2 mm.

This study determined that the factors most affecting the CT number were the tube voltage and treatment site (brain and H&N). Thus, the tube voltage (80–140 kVp) and treatment site (brain and H&N) were used to correct the dose distribution for inhomogeneity.

The sagittal view of the brain metastases treatment plan is parallel to the inferior orbitomeatal (IOM) line and across the frontal lobe (the top of the ROI), across the corpus callosum (the middle of the ROI), and across the cerebellum (the end of the ROI; Fig. 5). The sagittal view of the H&N tumor treatment plan system is parallel to the IOM line and across the maxillary sinus (the top of the ROI), across the mandible (the middle of the ROI), and across the C7 spinal disc (the end of the ROI). The gamma index was used to compare the dose distribution (Fig. 6).

Figure 5. Sagittal view of a brain metastases treatment plan (a) parallel to the IOM line and across the frontal lobe (the top of the ROI), (b) parallel to the IOM line and across the corpus callosum (the middle of the ROI), and (c) parallel to the IOM line and across the cerebellum (the end of the ROI). The local gamma passing rates for various passing criteria were compared for 1%/1 mm and show (d–f) ED calibrations of 80 and 140 kVp and (g–i) ED calibrations of 120 and 140 kVp. Passing criteria: gamma ≤1. Green 90.0%–100.0%, yellow 75.0%–90.0%, and red 0.0%–75.0%. IOM, inferior orbitomeatal; ROI, circular region of interest; ED, electron density.

Figure 6. Sagittal view of H&N tumor treatment plan system (a) parallel to the line and across the maxillary sinus (the top of the ROI), (b) parallel to the IOM line and across the mandible (the middle of the ROI), and (c) parallel to the IOM line and across the C7 (the end of the ROI). The local gamma passing rates for various passing criteria were compared for 1%/1 mm and show (d–f) ED calibrations of 80 and 120 kVp and (g–i) ED calibrations of 120 and 140 kVp. Passing criteria: gamma ≤1. Green 90.0%–100.0%, yellow 75.0%–90.0%, and red 0.0%–75.0%. IOM, inferior orbitomeatal; ROI, circular region of interest; ED, electron density.

The local gamma passing rates for various passing criteria were compared for ED calibrations of 80 kVp and 140 kVp and 120 and 140 kVp for various brain metastases measured using a dose distribution dosimeter. The results are shown in Table 2. The measurement across the frontal lobe dosimeter for the passing criterion of 1%/1 mm showed mean gamma passing rates of 74.48±15.39 and 95.64±5.95, respectively (P<0.001). The measurement across the corpus callosum dosimeter for the passing criterion of 1%/1 mm showed mean gamma passing rates of 96.35±2.67 and 98.76±3.03, respectively (P=0.024). The measurement across the cerebellum dosimeter for the passing criterion of 1%/1 mm showed mean gamma passing rates of 79.32±6.48 and 95.50±3.77, respectively (P<0.001).

Table 2 . The local gamma passing rates of the various passing criteria for various brain metastases measured.

Treatment sitesPassing criterionLocal gamma passing rates (%)P-value

80 and 140 kVp120 and 140 kVp
Brain tumor
Across the frontal lobe2%/2 mm96.67±4.9699.53±0.670.009
1%/1 mm74.48±15.3995.64±5.95<0.001
0.5%/0.5 mm42.14±24.3664.31±21.420.007
Across the corpus callosum2%/2 mm99.70±0.4699.90±0.160.113
1%/1 mm96.35±2.6798.76±3.030.024
0.5%/0.5 mm81.65±7.4289.35±9.150.014
Across the cerebellum2%/2 mm98.14±1.7399.90±0.130.001
1%/1 mm79.32±6.4895.50±3.77<0.001
0.5%/0.5 mm39.40±13.5162.78±8.150.005

Values are presented as mean±standard deviation..



The local gamma passing rates for various passing criteria for ED calibrations of 80 and 140 kVp and 120 and 140 kVp for various H&N tumors measured using a dose distribution dosimeter are shown in Table 3. The measured value across the maxillary sinus dosimeter for the passing criterion of 1%/1 mm showed mean gamma passing rates of 88.69±13.31 and 98.96±1.56, respectively (P=0.019). The measured value across the mandible dosimeter for the passing criterion of 1%/1 mm showed mean gamma passing rates of 89.15±9.41 and 97.94±2.06, respectively (P=0.001). The measured value crossing the C7 spinal segment dosimeter for the passing criterion of 1%/1 mm showed mean gamma passing rates of 92.21±13.38 and 99.85±0.41, respectively (P=0.030).

Table 3 . The local gamma passing rates of the various passing criteria for various head and neck tumors measured.

Treatment sitePassing criterionLocal gamma passing rates (%)P-value

80 and 140 kVp120 and 140 kVp
H&N tumor
Across the maxillary sinus2%/2 mm98.50±2.2199.94±0.780.022
1%/1 mm88.69±13.3198.96±1.560.019
0.5%/0.5 mm65.60±18.2890.63±9.35<0.001
Across the mandible2%/2 mm97.59±2.4899.38±1.400.018
1%/1 mm89.15±9.4197.94±2.060.001
0.5%/0.5 mm69.85±17.6989.05±7.48<0.001
Across the C7 spinal segment2%/2 mm99.39±1.01100.00±0.000.022
1%/1 mm92.21±13.3899.85±0.410.030
0.5%/0.5 mm80.98±23.8795.50±9.220.031

Values are presented as mean±standard deviation..



The local gamma passing rates for various passing criteria for ED calibrations with different treatment techniques (VMAT and IMAT) were compared as measured using a dose distribution dosimeter with PTW–VeriSoft® measurements. For the passing criterion of 1%/1 mm, the mean gamma passing rates for the VMAT and IMRT were 96.03% and 97.24%, respectively (P=0.40). For the passing criterion of 2%/2 mm, the mean gamma passing rates for the VMAT and IMRT were 99.37% and 99.56%, respectively (P=0.51).

Discussion

The CT number depends on the CT scanner image protocol conditions. In a previous study, the specifications and standards for a radiotherapy department at a regional hospital in southern Taiwan were assessed [18]. The study found that the smallest deviations in the relationship between the CT number and ED occurred for low ED materials and increased with increasing density. Thomas [19] noted that this deviation was <1% for materials with a low ED for different calibration techniques and manufacturers.

The CT number in this study was shown to be stable concerning the different acquisition parameters except for low tube voltage and materials with a higher ED, which resulted in deviations of about ±21.2%.

The results of a similar study assessing the relationship between the CT number and relative ED for various tube voltages reported a maximum difference in the mean dose of 1.1% [20]. Zurl et al. [11] used different image protocol conditions to assess variation in the HU of 20%, which resulted in an average dose error of approximately 1.5%. Moreover, Van Dyk et al. [21] and Gershkevitsh et al. [22] showed that a low tube voltage and an increased thickness of the anatomical structure being scanned resulted in a more accurate CT density conversion table that avoided calculation errors >2%. Kilby et al. [23] showed that radiotherapy at a reasonable dose level was possible for different CT imaging protocol conditions using an average CT density conversion table.

Taken together, these studies support that the factors most affecting the CT number are the tube voltage and treatment site. Given these findings, the gamma index analysis was used to compare the dose distribution in real clinical cases. Results showed that the gamma passing rates for ED calibrations of 80, 120, and 140 kVp with the passing criterion of 1%/1 mm were significant for the four treatment sites. Venselaar et al. [24] showed that the passing criterion of 2% or 2 mm in the TPS dose calculations in the central beam was acceptable. This study also noted a maximum difference of only gamma passing rate <90% for the passing criterion of 1%/1 mm for ED calibrations of 80 and 140 kVp at the treatment site (brain and H&N). Thus, the results were clinically insignificant.

du Plessis et al. [9] reported an average dose error of about 20% to 70% in the maxillary sinus region for a head model for different calibration techniques and manufacturers. The anatomical effects of homogeneities must be considered because the brain and H&N tumor sites are more likely to exhibit an uneven dose distribution. Thus, using the correct CT number–ED conversion table is essential.

Gamma passing rates were used to compare 80, 120, and 140 kVp for brain metastases and H&N tumors using the passing criteria of 1%/1 and 2%/2 mm. The results for both regions were statistically significant. The gamma passing rates were >90% and <90% for the passing criteria of 2%/2 and 1%/1 mm, respectively. Zurl et al. [11] used initial calculations for inhomogeneous features of tissue to show that the bone produces more statistically significant results than soft tissue. The skull almost entirely covers the brain tissue for the brain and H&N tumors at lower tube voltages. Thus, X-ray absorption is important when the X-ray beams traverse inhomogeneous areas. While there are very small differences for the passing criterion of 1%/1 mm for other treatment sites, the difference is greater for brain and H&N tumors.

In general, the CT number is most significant for high-density materials. For CT scanning parameter settings, adjusting the tube voltage has the most significant effect on the brain and H&N tumor imaging. Yang et al. [25] developed the Swiss Cheese model in 1990 to show that radiological errors can be caused by a series of mistakes. When a medical radiological technologist performs a CT simulation, the operating procedure checklist must be executed to ensure consistency in scanning parameters and to ensure that the correct tube voltage is used to prevent deviation in the TPS dose distribution [25].

Conclusions

The use of different CT scanning protocols resulted in a variation of up to 21.2% in the HU values. The local gamma passing rates for various passing criteria for ED calibrations of 80 and 140 kVp were calculated using real patient data, which revealed that brain metastases and H&N tumors showed a gamma passing rate of <90% for the passing criterion of 1%/1 mm. Regarding standardized scanning parameter settings, a specific conversion table needs to be created if the tube voltage in the CT scanner is adjusted. The correct protocol should be used for the corresponding TPS based on the treatment site, especially because this affects the dose distribution for brain metastases and H&N tumor recognition. Such steps could help reduce systematic errors.

Acknowledgements

The authors are greatly thankful to Mrs. Lo of Chiu Ho Medical System Co., LTD, Taiwan, for the knowledge on Pinnacle TPS. This study was supported by the Chi Mei Medical Center, Liouying, from January 2020 (grants CLFHR10914).

Conflicts of Interest


The authors have nothing to disclose.

Availability of Data and Materials


The data that support the findings of this study are available on request from the corresponding author.

Author Contributions


Conceptualization: Shu-Chin Yang, Su-Hua Lo, Li-Tsuen Shie, Sung-Wei Lee, and Sheng-Yow Ho. Data curation: Shu-Chin Yang, Su-Hua Lo, and Li-Tsuen Shie. Formal analysis: Shu-Chin Yang and Su-Hua Lo. Funding acquisition: Shu-Chin Yang. Investigation: Shu-Chin Yang and Su-Hua Lo. Methodology: Shu-Chin Yang, Su-Hua Lo, Li-Tsuen Shie, Sung-Wei Lee, and Sheng-Yow Ho. Project administration: Shu-Chin Yang and Su-Hua Lo. Resources: Sung-Wei Lee. Software: Shu-Chin Yang and Su-Hua Lo. Supervision: Shu-Chin Yang. Validation: Shu-Chin Yang and Su-Hua Lo. Visualization: Shu-Chin Yang, Su-Hua Lo, Li-Tsuen Shie, Sung-Wei Lee, and Sheng-Yow Ho. Writing–original draft: Shu-Chin Yang and Su-Hua Lo. Writing–review & editing: Shu-Chin Yang, Su-Hua Lo, and Li-Tsuen Shie.

Fig 1.

Figure 1.The cylindrical Cheese Virtual WaterTM phantom.
Progress in Medical Physics 2021; 32: 59-69https://doi.org/10.14316/pmp.2021.32.3.59

Fig 2.

Figure 2.The computed tomography (CT) number–electron density conversion tables were split into two ranges: <1.1 and >1.1 g/cm3 densities. A straight-line fit was constructed for each range. The mean and standard deviation for each range and set of scanning conditions were used as a measure of the CT number relative to the CT scanning parameters.
Progress in Medical Physics 2021; 32: 59-69https://doi.org/10.14316/pmp.2021.32.3.59

Fig 3.

Figure 3.The computed tomography (CT) number–electron density conversion tables using different tube voltages. For densities >1.1 g/cm3, the standard deviation for the CT number was ±21.2%. For densities <1.1 g/cm3, the standard deviation for the CT number was ±0.6%.
Progress in Medical Physics 2021; 32: 59-69https://doi.org/10.14316/pmp.2021.32.3.59

Fig 4.

Figure 4.Local gamma passing rates with various passing criteria calculated from the PTW–VeriSoft® dosimeter measurements are shown for (a) the 80 and 140 kVp plans and (b) the 120 and 140 kVp plans. Global gamma passing rates with various passing criteria calculated from the PTW–VeriSoft® dosimeter measurements are shown for (c) the 80 and 140 kVp plans and (d) the 120 s and 140 kVp plans. The dosimetric effects for the four treatment sites (brain metastases plans, head and neck cancer plans, lung cancer plans, and prostate cancer plans) were analyzed. The error bars above the histogram indicate the standard deviation.
Progress in Medical Physics 2021; 32: 59-69https://doi.org/10.14316/pmp.2021.32.3.59

Fig 5.

Figure 5.Sagittal view of a brain metastases treatment plan (a) parallel to the IOM line and across the frontal lobe (the top of the ROI), (b) parallel to the IOM line and across the corpus callosum (the middle of the ROI), and (c) parallel to the IOM line and across the cerebellum (the end of the ROI). The local gamma passing rates for various passing criteria were compared for 1%/1 mm and show (d–f) ED calibrations of 80 and 140 kVp and (g–i) ED calibrations of 120 and 140 kVp. Passing criteria: gamma ≤1. Green 90.0%–100.0%, yellow 75.0%–90.0%, and red 0.0%–75.0%. IOM, inferior orbitomeatal; ROI, circular region of interest; ED, electron density.
Progress in Medical Physics 2021; 32: 59-69https://doi.org/10.14316/pmp.2021.32.3.59

Fig 6.

Figure 6.Sagittal view of H&N tumor treatment plan system (a) parallel to the line and across the maxillary sinus (the top of the ROI), (b) parallel to the IOM line and across the mandible (the middle of the ROI), and (c) parallel to the IOM line and across the C7 (the end of the ROI). The local gamma passing rates for various passing criteria were compared for 1%/1 mm and show (d–f) ED calibrations of 80 and 120 kVp and (g–i) ED calibrations of 120 and 140 kVp. Passing criteria: gamma ≤1. Green 90.0%–100.0%, yellow 75.0%–90.0%, and red 0.0%–75.0%. IOM, inferior orbitomeatal; ROI, circular region of interest; ED, electron density.
Progress in Medical Physics 2021; 32: 59-69https://doi.org/10.14316/pmp.2021.32.3.59

Table 1 Analyzed results of the computed tomography scanning parameters

Scan protocolParameter

Scan techniqueTube voltage (kVp)Tube current (mA)Slice thickness (mm)FOV (cm)
BrainAxial80, 120, and 140300530, 50, and 65
Head and neckHelical80, 120, and 140350530, 50, and 65
ChestHelical80, 120, and 140400530, 50, and 65
PelvicHelical80, 120, and 140400530, 50, and 65

FOV, field of view.


Table 2 The local gamma passing rates of the various passing criteria for various brain metastases measured

Treatment sitesPassing criterionLocal gamma passing rates (%)P-value

80 and 140 kVp120 and 140 kVp
Brain tumor
Across the frontal lobe2%/2 mm96.67±4.9699.53±0.670.009
1%/1 mm74.48±15.3995.64±5.95<0.001
0.5%/0.5 mm42.14±24.3664.31±21.420.007
Across the corpus callosum2%/2 mm99.70±0.4699.90±0.160.113
1%/1 mm96.35±2.6798.76±3.030.024
0.5%/0.5 mm81.65±7.4289.35±9.150.014
Across the cerebellum2%/2 mm98.14±1.7399.90±0.130.001
1%/1 mm79.32±6.4895.50±3.77<0.001
0.5%/0.5 mm39.40±13.5162.78±8.150.005

Values are presented as mean±standard deviation.


Table 3 The local gamma passing rates of the various passing criteria for various head and neck tumors measured

Treatment sitePassing criterionLocal gamma passing rates (%)P-value

80 and 140 kVp120 and 140 kVp
H&N tumor
Across the maxillary sinus2%/2 mm98.50±2.2199.94±0.780.022
1%/1 mm88.69±13.3198.96±1.560.019
0.5%/0.5 mm65.60±18.2890.63±9.35<0.001
Across the mandible2%/2 mm97.59±2.4899.38±1.400.018
1%/1 mm89.15±9.4197.94±2.060.001
0.5%/0.5 mm69.85±17.6989.05±7.48<0.001
Across the C7 spinal segment2%/2 mm99.39±1.01100.00±0.000.022
1%/1 mm92.21±13.3899.85±0.410.030
0.5%/0.5 mm80.98±23.8795.50±9.220.031

Values are presented as mean±standard deviation.


References

  1. Das IJ, Cheng CW, Cao M, Johnstone PA. Computed tomography imaging parameters for inhomogeneity correction in radiation treatment planning. J Med Phys. 2016;41:3-11.
    Pubmed KoreaMed CrossRef
  2. Geise RA, McCullough EC. The use of CT scanners in megavoltage photon-beam therapy planning. Radiology. 1977;124:133-141.
    Pubmed CrossRef
  3. Goodsitt MM, Christodoulou EG, Larson SC. Accuracies of the synthesized monochromatic CT numbers and effective atomic numbers obtained with a rapid kVp switching dual energy CT scanner. Med Phys. 2011;38:2222-2232.
    Pubmed CrossRef
  4. Lv P, Liu J, Zhang R, Jia Y, Gao J. Combined use of automatic tube voltage selection and current modulation with iterative reconstruction for CT evaluation of small hypervascular hepatocellular carcinomas: effect on lesion conspicuity and image quality. Korean J Radiol. 2015;16:531-540.
  5. Guan H, Yin FF, Kim JH. Accuracy of inhomogeneity correction in photon radiotherapy from CT scans with different settings. Phys Med Biol. 2002;47:N223-N231.
    Pubmed CrossRef
  6. Davis AT, Palmer AL, Nisbet A. Can CT scan protocols used for radiotherapy treatment planning be adjusted to optimize image quality and patient dose? A systematic review. Br J Radiol. 2017;90:20160406.
    Pubmed KoreaMed CrossRef
  7. Cropp RJ, Seslija P, Tso D, Thakur Y. Scanner and kVp dependence of measured CT numbers in the ACR CT phantom. J Appl Clin Med Phys. 2013;14:4417.
    Pubmed KoreaMed CrossRef
  8. Skrzyński W, Zielińska-Dabrowska S, Wachowicz M, Slusarczyk-Kacprzyk W, Kukołowicz PF, Bulski W. Computed tomography as a source of electron density information for radiation treatment planning. Strahlenther Onkol. 2010;186:327-333.
    Pubmed CrossRef
  9. du Plessis FC, Willemse CA, Lötter MG, Goedhals L. Comparison of the Batho, ETAR and Monte Carlo dose calculation methods in CT based patient models. Med Phys. 2001;28:582-589.
    Pubmed CrossRef
  10. Cozzi L, Fogliata A, Buffa F, Bieri S. Dosimetric impact of computed tomography calibration on a commercial treatment planning system for external radiation therapy. Radiother Oncol. 1998;48:335-338.
    Pubmed CrossRef
  11. Zurl B, Tiefling R, Winkler P, Kindl P, Kapp KS. Hounsfield units variations: impact on CT-density based conversion tables and their effects on dose distribution. Strahlenther Onkol. 2014;190:88-93.
    Pubmed CrossRef
  12. Görlitz E. Dosimetrische und verfahrenstechnische Untersuchungen zur Qualitätssicherung eines Bestrahlungsplanungsprogramms. [Thesis]. Hambur: Universität Hamburg. 2006.
  13. Low DA, Harms WB, Mutic S, Purdy JA. A technique for the quantitative evaluation of dose distributions. Med Phys. 1998;25:656-661.
    Pubmed CrossRef
  14. Cheng A, Gerber RL, Wong JW, Purdy JA; Harms WB Sr. Systematic verification of a three-dimensional electron beam dose calculation algorithm. Med Phys. 1996;23:685-693.
    Pubmed CrossRef
  15. Shiu AS, Tung S, Hogstrom KR, Wong JW, Gerber RL, Harms WB, et al. Verification data for electron beam dose algorithms. Med Phys. 1992;19:623-636.
    Pubmed CrossRef
  16. William BH, Daniel AL, James AP. A quantitative software analysis tool for verifying 3-d dose calculation programs. Int J Radiat Oncol Biol Phys. 1994;30(Suppl 1):187.
  17. Ezzell GA, Galvin JM, Low D, Palta JR, Rosen I, Sharpe MB, et al. Guidance document on delivery, treatment planning, and clinical implementation of IMRT: report of the IMRT Subcommittee of the AAPM Radiation Therapy Committee. Med Phys. 2003;30:2089-2115.
    Pubmed CrossRef
  18. Fraass B, Doppke K, Hunt M, Kutcher G, Starkschall G, Stern R, et al. American Association of Physicists in Medicine Radiation Therapy Committee Task Group 53: quality assurance for clinical radiotherapy treatment planning. Med Phys. 1998;25:1773-1829.
    Pubmed CrossRef
  19. Thomas SJ. Relative electron density calibration of CT scanners for radiotherapy treatment planning. Br J Radiol. 1999;72:781-786.
    Pubmed CrossRef
  20. Rhee DJ, Kim S, Jeong DH, Moon YM, Kim JK. Effects of the difference in tube voltage of the CT scanner on dose calculation. J Korean Phys Soc. 2015;67:123-128.
  21. Van Dyk J, Barnett RB, Cygler JE, Shragge PC. Commissioning and quality assurance of treatment planning computers. Int J Radiat Oncol Biol Phys. 1993;26:261-273.
    Pubmed CrossRef
  22. Gershkevitsh E, Schmidt R, Velez G, Miller D, Korf E, Yip F, et al. Dosimetric verification of radiotherapy treatment planning systems: results of IAEA pilot study. Radiother Oncol. 2008;89:338-346.
    Pubmed CrossRef
  23. 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. 2002;47:1485-1492.
    Pubmed CrossRef
  24. Venselaar J, Welleweerd H, Mijnheer B. Tolerances for the accuracy of photon beam dose calculations of treatment planning systems. Radiother Oncol. 2001;60:191-201.
    Pubmed CrossRef
  25. Yang SC, Chang YL, Chang LY, Ho SY. The effectiveness of skill learning assessment using direct observation of procedural skills (DOPS) in the training of post-graduate training year and internship students of radiological technologist. C J Radiologic Tech. 2018;42:89-96.
Korean Society of Medical Physics

Vol.35 No.3
September 2024

pISSN 2508-4445
eISSN 2508-4453
Formerly ISSN 1226-5829

Frequency: Quarterly

Current Issue   |   Archives

Stats or Metrics

Share this article on :

  • line