- Open Access
The use of positron emission tomography/computed tomography imaging in radiation therapy: a phantom study for setting internal target volume of biological target volume
Radiation Oncology volume 10, Article number: 1 (2015)
Fluorodeoxyglucose (18F-FDG) positron emission tomography/computed tomography (PET/CT) is an important method for detecting tumours, planning radiotherapy treatment, and evaluating treatment responses. However, using the standardized uptake value (SUV) threshold with PET imaging may be suitable not to determine gross tumour volume but to determine biological target volume (BTV). The aim of this study was to extract internal target volume of BTV from PET images.
Three spherical densities of 18F-FDG were employed in a phantom with an air or water background with repetitive motion amplitudes of 0–30 mm. The PET data were reconstructed with attenuation correction (AC) based on CT images obtained by slow CT scanning (SCS) or helical CT scanning (HCS). The errors in measured SUVmax and volumes calculated using SUV threshold values based on SUVmax (THmax) in experiments performed with varying extents of respiratory motion and AC were analysed.
A partial volume effect (PVE) was not observed in spheres with diameters of ≥ 28 mm. When calculating SUVmax and THmax, using SCS for AC yielded smaller variance than using HCS (p < 0.05). For spheres of 37- and 28-mm diameters in the phantom with either an air or water background, significant differences were observed when mean THmax of 30-, 20-, or 10-mm amplitude were compared with the stationary conditions (p < 0.05). The average THmax values for 37-mm and 28-mm spheres with an air background were 0.362 and 0.352 in non-motion, respectively, and the mean THmax values for 37-mm and 28-mm spheres with a water background were 0.404 and 0.387 in non-motion and 0.244 and 0.263 in motion, respectively. When the phantom background was air, regardless of sphere concentration or size, THmax was dependent only on motion amplitude.
We found that there was no PVE for spheres with ≥ 28-mm diameters, and differences between SUVmax and THmax were reduced by using SCS for AC. In the head-and-neck and the abdomen, the standard values of THmax were 0.25 and 0.40 with and without respiratory movement, respectively. In the lungs, the value of THmax became the approximate expression depending on motion amplitude.
Positron emission tomography/computed tomography (PET/CT) has become an important and commonly used imaging method to detect tumours, plan radiotherapy treatment, and evaluate responses to therapy, as it allows for simultaneous functional and anatomical imaging. During treatment planning for radiation therapy, gross tumour volume (GTV) is calculated based on CT images. Moreover, it has been reported that when planning treatment for tumours of the head-and-neck, inter- and intra-observer differences in GTV may be reduced by using PET images .
There are many issues that arise when determining GTV with only PET imaging, such as accounting for respiratory movement and distinguishing inflammation. PET imaging also involves limited spatial resolution and the partial volume effect (PVE). Moreover, it has been reported that extraction volume differs according to the standardized uptake values (SUV) of normal tissue and tumour . Thus, it is not suitable to automatically extract GTV using the SUV from PET images. Realistically, it is preferable that the radiotherapist uses the PET images as a reference .
Given that the PET image was a function image, the extracted volume is thought to be biological target volume (BTV). BTV often indicates the part with especially high tumour revitalization in GTV . It became possible to irradiate GTV and BTV by a separate dose by the development of the radiation therapy machine and the plan device. Therefore, BTV is being used for radiation therapy .
Several research reports suggest extracting BTV from PET images [6-12]. However, a concrete extraction condition that tempers with various situations has not been established to date. The purpose of this study was to obtain conditions to automatically extract internal target volume of BTV (ITV-b) from phantom experiments using SUV threshold values based on SUVmax (THmax). The items that should be considered for extracting ITV-b using THmax are as follows: (1) tumour size (PVE), (2) motion amplitude, (3) tumour site (lung or head-and-neck, abdomen), and (4) attenuation correction (AC). We examined tumour sizes for which a PVE was not observed [13,14] in addition to changes in SUVmax and THmax according to the degree of motion amplitude [15-17], the position of the tumour , and the method of attenuation correction [19-21].
The National Electrical Manufacturers Association (NEMA) 2001 International Electrotechnical Commission (IEC) Phantom
In this study, we imaged a PET image quality phantom (NEMA/IEC phantom; Data Spectrum Corp.). This phantom contained six hollow spheres with internal diameters of 10, 13, 17, 22, 28, and 37 mm, inside a simulated body cavity that may be filled with water. The phantom background was filled with either water or air. The phantom background was filled with water to simulate a human head-and-neck, abdominal and pelvic tumour. In this series of experiments, all spheres were filled with 2-[fluorine-18] fluoro-2-deoxy-D-glucose (18F-FDG) solution that was measured in a dose calibrator and diluted with water to obtain a radioisotope concentration (RC) of 16–64 kBq/mL. The total weight of the filled phantom was 12.4 kg. The phantom background was filled with air to simulate a human lung tumour. In this series of experiments, all spheres were filled with 18F-FDG solution diluted with water to obtain an RC of 64–256 kBq/mL. The total weight of the phantom in this situation was 2.4 kg.
The phantom background was filled with water or air, yielding sphere-to-background (S/B) ratios of approximately 4, 8, and 16 to 1 or null. To achieve the same S/B ratio, each given dose was adjusted by the difference of the background. The RC values of the spheres and background and the S/B ratios were chosen to simulate the ranges obtained from clinical conditions based on the tumour concentration of 18F-FDG and the SUV. The SUV is commonly used as a relative measure of FDG uptake. The basic formula for SUV was as follows:
where C was the RC (kBq/mL) measured by the PET scanner within the region of interest, A was the amount of injected activity (kBq), and W was the weight of the patient (g).
A motion system (QRP-01; Qualita, Ageo City, Saitama Prefecture, Japan) was used in this study to simulate respiratory motion. The NEMA/IEC phantom was placed on a moving table. The moving table was oscillated with a displacement of 10, 20, and 30 mm in the cranio-caudal direction at a frequency of 15 cycles/min. These parameters were selected to simulate displacements and respiratory cycles typically observed in normal respiratory motion.
CT transmission (for AC) and PET scans were performed using a PET/CT scanner (Discovery PET/CT 600 Motion; General Electric Medical Systems [GEMS], Milwaukee, WI, USA) with a stationary or moving phantom. Scans of the moving phantom were acquired without gating. Slow CT scanning (SCS; 120 kV, 45 mA, and 4.0 s/rotation) or the helical CT scanning (HCS; 120 kV, 360 mA, 0.5 s/rotation, and pitch 0.938) imaging data, followed by PET emission data, were acquired for 300 s from one field of view. Both PET scan types (stationary and motion) were acquired in 3-dimensional (3D) mode. The 3D data were first Fourier rebinned, and then all emission scans were reconstructed into a 256 matrix using ordered-subset expectation maximization algorithms. The numbers of iterations and subsets were 2 and 32, respectively. The PET-reconstructed slice thickness was 3.75 mm.
The NEMA/IEC phantom, which contained six spheres ranging from 10 to 37 mm in diameter, was filled with 18F-FDG. PET/CT data acquired with the stationary phantom was used as a standard (non-motion). Respiratory motion was simulated by a motor-driven plastic platform that moved the phantom with a displacement of 10, 20, and 30 mm in the cranio-caudal direction at a frequency of 15 cycles/min. A 4-s period of oscillation was selected based on the average respiratory period measured in 50 patients with lung tumours at our centre using a respiratory gating system (Real-time Positioning Management System; Varian, Palo Alto, CA, USA). With the phantom in motion, PET data were acquired in the static mode, wherein data were collected for a definite time period according to the width of the detector’s body axis instead of continuous movement like in CT. Acquisitions using three types of S/B ratios, four different motion amplitudes, two different background media, and two types of CT scans were repeated three times; for a total of 144 image acquisitions performed.
Image analysis was performed using PET VCAR™ software and the Advantage Workstation (GE Healthcare, Milwaukee, WI, USA). This software is designed to measure SUV and extracted volume for any PET-defined region with detected activity. SUVmax values were measured in spheres with 144 image acquisitions. The SUV threshold values based on SUVmax (THmax) were calculated so that the volume extracted by the THmax matched the ideal volume of the spherical tumours in the phantom, taking into consideration the motion amplitude. Thus, an ideal phantom volume becomes a total of the volume of spheres and the volume of the column corresponding to the amount of the amplitude (Figure 1).
Evaluation of PVE
To correctly assess the influence of respiratory movement, it is necessary to use a volume that will not be affected by the PVE. Therefore, the PVE was evaluated using the data from scans with the stationary phantom, with the recovery coefficient (RCo) calculated as follows:
where Ci was the maximum radioactivity of each sphere at each S/B ratio, and Cr was the maximum radioactivity of the 37-mm sphere at each S/B ratio in non-motion, both measured in scans with no phantom motion. Evaluation was performed using two types of background (air and water) and attenuation corrections (SCS and HCS). An analysis of variance (ANOVA) was used to check for significant differences in the sphere diameter. Next, multiple comparisons between six types of spheres of (10–37 mm) were performed using Dunnett’s method and the size without the influence of PVE was then examined.
Analysis of SUVmax and THmax
Effects of sphere size, motion amplitude, background medium, and attenuation correction on measured SUVmax and THmax were examined. The spheres with 37-mm and 28-mm diameters were used for this evaluation. The phantom, with backgrounds of both air and water, were scanned with the SCS and HCS, as described above. During scanning, the phantom was moved with amplitudes of 0, 10, 20, or 30 mm. To evaluate the influence of all combined S/B ratios on SUVmax, the ratios of SUVmax (SUVr) were calculated using the average SUVmax for the three S/B ratios in scans performed with no phantom motion.
The SUVr was calculated as follows:
where Cj was the average of SUVmax values obtained for all S/B ratios with a particular sphere size, and Ca was the average of SUVmax values obtained for all S/B ratios with the same sphere size and no phantom motion. The variances in SUVr were compared between the attenuation corrections using the F-test. The differences between the mean SUVr obtained with different motion amplitudes were evaluated using ANOVA with Dunnett’s multiple comparisons.
The variances of THmax were compared with the attenuation corrections using the F-test. The variations of the mean THmax with different amplitudes were evaluated by ANOVA and Dunnett’s multiple comparisons.
Partial volume effects
Table 1 shows the RCo values measured in our study. A significant difference was observed between RCo values measured in spheres of different sizes with phantoms using either air or water as background, scanned with SCS or HCS. When evaluating all RCo values measured using each of the CT methods and both of the background media used, no PVE influence on RCo was observed in spheres with diameters of 28 mm or larger. Conversely, we demonstrated that the spheres with diameters of 22 mm or less might be affected by PVE. Therefore, we only analysed 28- and 37-mm spheres in subsequent experiments.
The SUVr measured in the phantom with an air background and spheres sized 37 and 28 mm are presented as box plots in Figure 2a and b, respectively. The variance of SUVr using the SCS CT technique was significantly smaller than that of SUVr obtained using the HCS approach. Figure 3a and b show box plots of the SUVr of spheres of 37 and 28-mm diameters, respectively, measured in a phantom with water background. The variance in SUVr calculated using the SCS CT technique was significantly smaller than that of SUVr obtained by HCS. Considering SUVr calculated in phantoms with air and water backgrounds together, the SCS CT technique was found to be a more suitable approach for AC, compared with HCS. Consequently, we only analysed SUVr calculated using the SCS CT technique in subsequent experiments.
The effect of motion amplitude on SUVr was then investigated. There was a significant difference between SUVr values obtained with different amplitudes. Conversely, the 10-mm motion amplitude did not affect the SUVr. With motion amplitudes of 20 mm and above, the decrease in SUVr was found to be relative to the amplitude of phantom motion. The largest observed SUVr difference, compared with the SUVr obtained under stationary conditions, was 30.9%, calculated for the 28-mm sphere in the phantom with water background at 30-mm motion amplitude (Figure 3b).
THmax values obtained in the phantom with an air background using 37- and 28-mm spheres are presented as box plots in Figure 4a and b, respectively. The variance of THmax values obtained using the SCS technique was significantly smaller than that obtained using the HCS approach. Figure 5a and b present box plots of the THmax calculated for 37- and 28-mm spheres, respectively, using measurements obtained in a phantom with a water background. For the 37- and 28-mm spheres, the variance of THmax did not differ significantly between SCS and HCS. When we consider the results of THmax for the phantoms with a water or air background together, the use of the SCS method for AC of PET data resulted in less variance, compared with the HCS approach. Therefore, in subsequent experiments we only analysed the THmax values derived from PET data with AC utilizing the SCS method.
The effect of motion amplitude on THmax was then investigated. We observed a significant difference between THmax values obtained at different amplitudes. For spheres of 37- and 28-mm diameters in the phantom with an air background, significant differences were compared with the stationary conditions (Figure 4a and b). The average THmax values for 37-mm and 28-mm spheres were 0.362 and 0.352 in non-motion, respectively. For 37- and 28-mm spheres with the water background phantom, significant differences were individually compared with the stationary conditions (Figure 5a and b). The average THmax value for 37-mm and 28-mm spheres were 0.404 and 0.387 in non-motion, respectively. Moreover, the average THmax value for 37-mm and 28-mm spheres were 0.244 and 0.263 in motion, respectively. When the phantom background was filled with air, regardless of the sphere concentration or size, THmax was dependent only on the motion amplitude. Therefore, it may be possible to approximate THmax using the following formula (Figure 6):
where y is the approximate value of THmax, and x is the motion amplitude.
The largest THmax deviation from the stationary conditions was 82.0%, calculated for the 37-mm diameter sphere with an air background phantom and 30-mm motion amplitude (Figure 4a). When the phantom background was water, no relation between the decrease in THmax and amplitude was observed. At amplitudes greater than 10 mm, THmax was reduced in comparison to the stationary conditions (Figure 5a and b).
In this study, we aimed to extract ITV-b from PET images using calculations from phantom studies. The results of the current study have demonstrated that tumour size, motion amplitude, tumour site (lung or head-and-neck, abdomen), and AC can have a significant impact on PET/CT imaging outcomes, including SUVmax and THmax. These four factors are individually discussed below.
In considering tumour size, we found no influence of PVE in the experiments with 28- and 37-mm spheres. In spheres with a diameter of 22 mm or less, reduced SUV estimates can result from PVE. Because the underestimation of SUV is exacerbated by respiratory movement, one should be careful when evaluating targets smaller than 28 mm by estimating the SUV. In cases with volumes larger than spheres of 37-mm diameter, the potential influence of PVE remains to be determined.
When evaluating the effect of target motion amplitude, we observed no SUVr underestimation for 10-mm respiratory motion. In the case of amplitudes of 20 mm or more, SUVr decreased in an amplitude-dependent manner. Moreover, differences in background affected calculated volumes. When the background was air, a significant amplitude-dependent difference in THmax was observed. When the background was water, there was no significant effect of amplitude on THmax. With amplitudes greater than 10 mm, the THmax was reduced compared with the stationary condition. THmax was found to be more strongly affected by motion than SUVmax. Because SUVmax is the value of a certain point in an ideal phantom volume, the impact of movement on it is less severe than the impact that movement has on THmax.
When considering the differences between simulated tumour sites (lung or head-and-neck, abdomen), head-and-neck and abdominal lesions were less susceptible to error caused by AC because the attenuation factors were uniform throughout water. In contrast, lung lesions may be expected to be more susceptible to errors arising from AC differences because of large density differences between sphere phantoms and air. Using air as the phantom background, regardless of RC or size, THmax was dependent only on amplitude. Therefore, in this context, it was possible to approximate THmax using the equation derived from the cubic equation fit to the data (Eq. 4). Using water as the background for the phantom, THmax was dependent on RC, size, and amplitude. When the average threshold of width in the value is selected, the error margin is included in the extraction capacity. However, we believe that this will become the standard for determining ITV-b.
When we consider differences between ACs, the variance of SUVr and THmax with SCS was not significantly different than that of SUVr and THmax with HCS with water as the phantom background. However, when the phantom background was air, the variance of SUVr and THmax with SCS was significantly smaller than that of SUVr and THmax with HCS. In evaluating the response to therapy, or for planning radiotherapy treatment using PET images, SUV and extracted volumes must not change between scans. Our data suggest that the SCS method provides superior AC in comparison to HCS.
Past phantom experiments indicate variation in the value of THmax, although they have been reported 40-50% in the non-motion and 25-35% in the moving phantoms [6-9]. THmax for tumours in different areas has also been reported. They are 15-50% for lung tumours [22,23], 40-50% for head-and-neck tumours , and 40% pelvic tumours . In summary, when ITV-b is extracted, the approximate expression based on motion amplitude can be used in the lung field. The approximate expression can be used by understanding the amount of the breath movement by the fluoroscope image. Because THmax has a width with the same amplitude in the head-and-neck and the abdomen, it is difficult to compute a concrete threshold value. However, we believe that extracting BTV using the mean value of THmax may aid in treatment planning. As for the value of THmax, 0.40 is the standard value in head-and-neck and lower abdominal areas without respiratory movement. Moreover, 0.25 is the standard value in upper abdominal areas with respiratory movement. However, these threshold values resemble those reported previously and become the first choice, it will be necessary to check the extracted volume.
In this study, we demonstrate concrete conditions to set ITV-b. There was no PVE with tumour sizes ≥ 28 mm, and the variation between SUVmax and THmax may be reduced by using SCS for AC. There was no decrease in respiratory movement from approximately 10 mm for SUVmax. As for the value of THmax, 0.40 became the standard value in areas without respiratory movement in the head-and-neck and the lower abdominal area, and 0.25 became the standard value in upper abdominal areas with respiratory movement. The approximate expression based on the motion amplitude may be adjusted in the lung area.
Positron emission tomography
Gross tumour volume
Partial volume effect
Standardized uptake value
Biological target volume
Internal target volume
- SUVmax :
Maximum standardized uptake value
Ratios of SUVmax
- THmax :
Threshold values based on SUVmax
National Electrical Manufacturers Association
International Electrotechnical Commission
Slow computed tomography scanning
Helical computed tomography scanning
Analysis of variance
Yu H, Caldwell C, Mah K, Poon I, Balogh J, MacKenzie R. Automated radiation targeting in head-and-neck cancer using region-based texture analysis of PET and CT images. Int J Radiat Oncol Biol Phys. 2009;75:618–25.
Le Maitre A, Hatt M, Pradier O, Cheze-le Rest C, Visvikis D. Impact of the accuracy of automatic tumour functional volume delineation on radiotherapy treatment planning. Phys Med Biol. 2012;57:5381–97.
MacManus M, Nestle U, Rosenzweig KE, Carrio I, Messa C, Belohlavek O, et al. Use of PET and PET/CT for radiation therapy planning: IAEA expert report 2006–2007. Radiother Oncol. 2009;91:85–94.
Jingu K, Takeda K, Metoki T, Ogawa Y, Ariga H, Yamada S, et al. The difference between gross tumor volume and biological target volume in postoperative local recurrent rectal cancer; a preliminary study of PET/CT radiation therapy planning. J Jpn Soc Ther Radiol Oncol. 2008;20:127–31.
Devic S, Tomic N, Faria S, Dean G, Lisbona R, Parker W, et al. Impact of 18FDG-PET/CT on biological target volume (BTV) definition for treatment planning for non-small cell lung cancer patients. Nucl Instrum Methods Phys Res Section A. 2007;571:89–92.
Erdi YE, Mawlawi O, Larson SM, Imbriaco M, Yeung H, Finn R, et al. Segmentation of lung lesion volume by adaptive positron emission tomography image thresholding. Cancer. 1997;80:2505–9.
Nagel CC, Bosmans G, Dekker AL, Ollers MC, De Ruysscher DK, Lambin P, et al. Phased attenuation correction in respiration correlated computed tomography/positron emitted tomography. Med Phys. 2006;33:1840–7.
Okubo M, Nishimura Y, Nakamatsu K, Okumura M, Shibata T, Kanamori S, et al. Static and moving phantom studies for radiation treatment planning in a positron emission tomography and computed tomography (PET/CT) system. Ann Nucl Med. 2008;22:579–86.
Yaremko B, Riauka T, Robinson D, Murray B, Alexander A, McEwan A, et al. Thresholding in PET images of static and moving targets. Phys Med Biol. 2005;50:5969–82.
Sakaguchi Y, Mitsumoto T, Zhang T, Mitsumoto K, Tachiya Y, Ohya N, et al. Importance of gated CT acquisition for the quantitative improvement of the gated PET/CT in moving phantom. Ann Nucl Med. 2010;24:507–14.
Toya R, Murakami R, Tashiro K, Yoshida M, Sakamoto F, Kawanaka K, et al. FDG-PET/CT-based gross tumor volume contouring for radiation therapy planning: an experimental phantom study. J Radiat Res. 2012;53:338–41.
Niyazi M, Landrock S, Elsner A, Manapov F, Hacker M, Belka C, et al. Automated biological target volume delineation for radiotherapy treatment planning using FDG-PET/CT. Radiat Oncol. 2013;8:180.
Liu C, Pierce 2nd LA, Alessio AM, Kinahan PE. The impact of respiratory motion on tumor quantification and delineation in static PET/CT imaging. Phys Med Biol. 2009;54:7345–62.
Geramifar P, Zafarghandi MS, Ghafarian P, Rahmim A, Ay MR. Respiratory-induced errors in tumor quantification and delineation in CT attenuation-corrected PET images: effects of tumor size, tumor location, and respiratory trace: a simulation study using the 4D XCAT phantom. Mol Imaging Biol. 2013;15:655–65.
Nehmeh SA, Erdi YE. Respiratory motion in positron emission tomography/computed tomography: a review. Semin Nucl Med. 2008;38:167–76.
Kasuya T, Tateishi U, Suzuki K, Daisaki H, Nishiyama Y, Hata M, et al. Role of respiratory-gated PET/CT for pancreatic tumors: a preliminary result. Eur J Radiol. 2013;82:69–74.
Callahan J, Kron T, Schneider-Kolsky M, Hicks RJ. The clinical significance and management of lesion motion due to respiration during PET/CT scanning. Cancer Imaging. 2011;11:224–36.
Park SJ, Ionascu D, Killoran J, Mamede M, Gerbaudo VH, Chin L, et al. Evaluation of the combined effects of target size, respiratory motion and background activity on 3D and 4D PET/CT images. Phys Med Biol. 2008;53:3661–79.
Didierlaurent D, Ribes S, Batatia H, Jaudet C, Dierickx LO, Zerdoud S, et al. The retrospective binning method improves the consistency of phase binning in respiratory-gated PET/CT. Phys Med Biol. 2012;57:7829–41.
Lamare F, Ledesma Carbayo MJ, Cresson T, Kontaxakis G, Santos A, Le Rest CC, et al. List-mode-based reconstruction for respiratory motion correction in PET using non-rigid body transformations. Phys Med Biol. 2007;52:5187–204.
Pevsner A, Nehmeh SA, Humm JL, Mageras GS, Erdi YE. Effect of motion on tracer activity determination in CT attenuation corrected PET images: a lung phantom study. Med Phys. 2005;32:2358–62.
Yu J, Li X, Xing L, Mu D, Fu Z, Sun X, et al. Comparison of tumor volumes as determined by pathologic examination and FDG-PET/CT images of non-small-cell lung cancer: a pilot study. Int J Radiat Oncol Biol Phys. 2009;75:1468–74.
Biehl KJ, Kong FM, Dehdashti F, Jin JY, Mutic S, Naqa EI, et al. 18F-FDG PET definition of gross tumor volume for radiotherapy of non-small cell lung cancer: is a single standardized uptake value threshold approach appropriate? J Nucl Med. 2005;61:1385–92.
Paulino AC, Koshy M, Howell R, Schuster D, Davis LW. Comparison of CT- and FDG-PET- defined gross tumor volume in intensity-modulated radiotherapy for head-and-neck cancer. Int J Radiat Oncol Biol Phys. 2009;75:1468–74.
Bassi MC, Turri L, Sacchetti G, Loi G, Cannillo B, La Mattina P, et al. FDG-PET/CT imaging for staging and target volume delineation in preoperative conformal radiotherapy of rectal cancer. Int J Radiat Oncol Biol Phys. 2008;70:1423–6.
This study was partly supported by the Foundation for Promotion of Hokkoku Cancer Research.
The authors declare that they have no competing interests.
WK drafted the manuscript and performed all mathematical analyses. SY participated in the sequence alignment. KY, KN, and KK participated in the design of the study and performed the statistical analysis. AT conceived of the study, participated in its design and coordination, and helped to draft the manuscript. All authors read and approved the final manuscript.
About this article
Cite this article
Kawakami, W., Takemura, A., Yokoyama, K. et al. The use of positron emission tomography/computed tomography imaging in radiation therapy: a phantom study for setting internal target volume of biological target volume. Radiat Oncol 10, 1 (2015) doi:10.1186/s13014-014-0315-2
- Respiratory motion
- Attenuation correction