Search tips
Search criteria 


Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
PLoS One. 2013; 8(5): e64665.
Published online 2013 May 16. doi:  10.1371/journal.pone.0064665
PMCID: PMC3655997

Improvement of Internal Tumor Volumes of Non-Small Cell Lung Cancer Patients for Radiation Treatment Planning Using Interpolated Average CT in PET/CT

Nils Cordes, Editor


Respiratory motion causes uncertainties in tumor edges on either computed tomography (CT) or positron emission tomography (PET) images and causes misalignment when registering PET and CT images. This phenomenon may cause radiation oncologists to delineate tumor volume inaccurately in radiotherapy treatment planning. The purpose of this study was to analyze radiology applications using interpolated average CT (IACT) as attenuation correction (AC) to diminish the occurrence of this scenario. Thirteen non-small cell lung cancer patients were recruited for the present comparison study. Each patient had full-inspiration, full-expiration CT images and free breathing PET images by an integrated PET/CT scan. IACT for AC in PETIACT was used to reduce the PET/CT misalignment. The standardized uptake value (SUV) correction with a low radiation dose was applied, and its tumor volume delineation was compared to those from HCT/PETHCT. The misalignment between the PETIACT and IACT was reduced when compared to the difference between PETHCT and HCT. The range of tumor motion was from 4 to 17 mm in the patient cohort. For HCT and PETHCT, correction was from 72% to 91%, while for IACT and PETIACT, correction was from 73% to 93% (*p<0.0001). The maximum and minimum differences in SUVmax were 0.18% and 27.27% for PETHCT and PETIACT, respectively. The largest percentage differences in the tumor volumes between HCT/PET and IACT/PET were observed in tumors located in the lowest lobe of the lung. Internal tumor volume defined by functional information using IACT/PETIACT fusion images for lung cancer would reduce the inaccuracy of tumor delineation in radiation therapy planning.


PET/CT combines F18-FDG positron emission tomography (PET), and computed tomography (CT) images with both functional and anatomic information provide a more precise diagnostic reference for tumor volume, tumor locations and tumor staging. Therefore, PET/CT has been increasingly used for target volume delineation in radiotherapy treatment planning (RTP) to deliver the optimal radiation dose to tumors and to decrease the radiation dose to surrounding normal tissues [1][3]. The reduction of intra- and inter-observer variability in target volume delineation by contouring with PET/CT has also been reported in previous studies [4][5]. Moreover, PET is a useful imaging tool to differentiate between inflammation and malignance, such as lung atelectasis, mediastinal lymphadenopathy, and distant metastases [6]. The incorporation of PET information in RTP along with the CT-based gross tumor volume can improve the definition of tumor volume and has been extensively used in radiotherapy.

Respiration produces additional variation in imaging diagnosis and target contouring for radiotherapy in thoracic malignancies. CT can provide high spatial information and attenuation correction for PET in PET/CT. However, respiratory motion causes uncertainties in the tumor edges on either CT or PET images and causes misalignment when registering PET and CT images. Consequently, this phenomenon may influence oncologists when attempting to delineate tumor volume accurately in RTP [3], [7][9].

Many studies have improved the misalignment in PET/CT fusion image due to respiration by using gated (4D) imaging techniques [3], [10][12]. With gated technology aids, the reduction of motion artifacts and increases in the accuracy of tumor volume and localization delineation when compared to non-gated PET were achieved. Moreover, we previously proposed an interpolation method with interpolated average CT (IACT) for attenuation correction (AC) to reduce the PET/CT misalignment [12][13]. Using IACT imaging, standardized uptake value (SUV) can be corrected with a lower radiation dose compared to the use of gated imaging. In the present study, the comparison of tumor volumes for RTP is reported. We assess the differences in tumor volumes between 3D PET/CT and PET/IACT and evaluate the SUVmax changes in terms of tumor locations.

Materials and Methods

Patient selection

With IRB approval (DMR98-IRB-171-1) for the application of 4D PET/CT to tumor delineation in RTP, thirteen non-small cell lung cancer patients were recruited for this comparison study. All patients signed written, informed consent. There were 9 tumors in the upper lobe and 4 tumors in the lower lobe. The clinical characteristics are summarized in Table 1.

Table 1
List of clinical characteristic for all patents.


FDG-PET/CT scans were obtained for tumor staging work-ups before cancer treatment. All patients had undergone the standard procedure of PET/CT (PET/CT-16 slice, Discovery STE, GE Medical System, Milwaukee, Wisconsin USA) scanning. Patients were injected with 370 MBq of 18F-FDG and rested during the pharmacokinetics uptake period. The original data included a series of helical CTs (HCT), two extreme phases CTs that considered the full-expiration and the full-inspiration CTs and a whole body PET. HCT images obtained with 120 kVp, variable mA (30–210 mA), 1.75[ratio]1 pitch, 8×3.75 mm x-ray collimation, and 0.5 gantry rotation time and gated CTs were acquired under the same conditions. PET data were acquired at the 3 min per 15 cm bed position; attenuation corrections of PETHCT and PETIACT utilized HCT and IACT, respectively.

Tumor motion and interpolated average CT

We previously proposed an IACT method from 4D-CT in comparison with 4D-Cine CT [12]. IACT is a robust, accurate low dose alternate to CACT and works well for a large range of respiratory motion amplitudes, as was reported in a simulation study [13]. The full-inspiration and full-expiration CT sets as two-extreme-phase images were used to generate the motion maps using the optical flow method (OFM), a deformable image registration algorithm [14]. The total motion range for each voxel in the forward motion map is equally spaced into 4 intervals, resulting in 3 sets of interpolated CT (ICT) image sets as the mid-phases from inspiration to expiration. (Figure 1). The 3 interpolated phases together with the two original phases, including one inhalation and expiration, compose a complete respiratory cycle. These 5 phases are averaged to generate the IACT for AC on PET data. The conclusion that the radiation dose using IACT could be reduced by 85% compared to that of 4D-CT was reported in previous study [12][13]. IACT serves as a low-dose alternative to 4D-CT. The OFM calculation was the following:

equation image

where n is the number of iterations, v(n) is the average velocity driven from the surrounding voxels, f(x, y, z, t) is the differentiable image intensity at position (x, y, z) at time t, and α is the weighting factor with an empirical value of 5. The given equations are applied to estimate the displacement for tumor motion between full-expiration and full-inspiration.

Figure 1
Illustrates the generation of IACT from full-inspiration and full-expiration images by OFM.

Tumor volume analysis

An experienced radiation oncologist manually delineated the internal tumor volume (TV) for all patients on HCT, IACT with fusion images of the individual PETHCT and PETIACT. Tumor volumes from HCTs (TVHCT) and from IACTs (TVIACT) consisted of the control group and the experiment group. Tumor volumes delineated with CT using CT/PET fusion information were compared with the percentage difference, which was calculated by the equation An external file that holds a picture, illustration, etc.
Object name is pone.0064665.e002.jpg. We also assessed the similarity between HCT/PETHCT and IACT/PETIACT fusion images using the ratio of the intersection to the union of TV. The correlation for compared TVs is defined as An external file that holds a picture, illustration, etc.
Object name is pone.0064665.e003.jpg, where A and B are the different tumor volumes from CT-based and PET-based images [15]. The metabolic rate of glucose, SUVmax from FDG-PET, was also applied to represent the physiology information within the tumor volume.


Figure 2 shows the PET/CT fusion for tumor contour delineation. The arrows indicate the mismatch observed in PETHCT/HCT fusion between PETHCT and HCT showing in the (A) transverse (B) coronal and (C) sagittal views. The correct image fusion with misalignment reduction on PETIACT and IACT represents the (D) transverse (E) coronal and (F) sagittal view. Table 2 shows the gross tumor volumes of HCT/PETHCT and IACT/PETIACT, along with estimation of the tumor motion, and the correction of HCT with PETHCT and IACT with PETIACT. The median (range) tumor volume for HCT/PET was 31 (4–169) mL, while for the corresponding IACT/PETIACT, the median tumor volume was 26 (3–149). The median tumor difference for HCT/PETHCT and IACT/PETIACT was 14% higher (range 5–24%). The range of tumor motion was from 4 to 17 mm. HCT and PETHCT correction was from 72% to 91%, while IACT and PETIACT was from 73% to 93% (*p<0.0001). Table 3 shows the SUVmax measurement for each tumor from PETHCT and PETIACT. The median SUVmax was about 9.65 (1.98–8.77) and 9.53 (1.97–18.59) for PETHCT and PETIACT, respectively. The maximum and minimum differences among SUVmax were 0.18% and 27.27%, respectively.

Figure 2
Represents the PET/CT fusion for tumor contour delineation.
Table 2
Tumor volumes and correlation.
Table 3
SUVmax for tumor volumes and difference.


IACT used for attenuation in PET is able to resolve the CT/PET fusion misalignment of thoracic tumors caused by respiration. In this study, we investigated the potential tumor volume delineation using IACT/PETIACT for RTP in Table 2. Compared to tumor volume determined by HCT/PET, larger TVs were observed in all the subjects, and the percentage difference was from 5% to 24%. For PET, imaging is obtained during several breathing cycles and it represents a time-average map. HCT is obtained during a very short period, depending on the machine. The tumor appears in HCT and PET images without considering respiratory motion correction, increasing the uncertainty of the tumor boundaries [3]. Therefore, isotropic extension of the internal tumor volume is usually utilized for adequate coverage of tumors. Considering the tumor contour correction between CTs and PETs, IACT and PETIACT are better correlated than HCT and PETHCT. Because IACT/PETIACT aids in defining tumor volume for RTP, tumor volume delineation can be performed more accurately.

The thoracic tumor volume defined by its functional region using PET for RTP is still limited by respiratory motion, which often increases the real tumor size and reduces the SUV. Several studies suggest that gated PET/CT is a better solution for defining the physiological extent of moving tumors and to improve RTP for lung cancers [15][16]. In our study, the increase of the SUVmax in PETIACT compared to PETHCT was not obvious, as can be seen in Table 3. We found that the reason for this issue was because of the low electron density resulting from the IACT in the averaged ICTs. Because SUV for PET is based on the attenuation coefficient in terms of electron density, the increase in the SUVmax within the functional tumor volume was limited.

In addition, in Figure 3, we showed the percentage differences of the tumor volume and of SUVmax between delineations from HCT/PETHCT and IACT/PETIACT according to tumor location inside the lung. It is intuitively understood that tumors located in the lower lobe of the lung, which is closer to the diaphragm, move more than those in the upper lobe. The larger motion causes more motion blurring in both CT and PET, and therefore, it is possible to extract an inappropriate tumor volume using HCT/PET. Our results show that the percentage difference was more significant for the tumors located in the lower lobe for patents 11–13. A lager respiration motion was observed in the lower lobe of the lung.

Figure 3
Shows the percentage difference in the tumor volume and the SUVmax between delineations from HCT/PETHCT and IACT/PETIACT versus the tumor location inside the lung.

There are two limitations to the proposed method. IACT was averaged from the ICTs generated under the equal timing. The ICTs, as with real mid-phase CT, were only generated under the assumption that there was smooth breath from the patients. Second, the full-expiration CT and the full-inspiration CT acquisition can be performed on patients with normal lung function. For patients who are not able to hold their breath for imaging, the proposed interpolating method is limited. In addition, the lack of a ground truth in tumor volume delineation makes difficulty to evaluate the accuracy of the presented IACT method. However, IACT/PETIACT including the respiration information is still superior to the traditional 3D-PET fusion images for lung tumor delineation.


Our results suggest that tumor volume defined by PET using IACT/PETIACT fusion images for lung cancer would reduce the inaccuracy of tumor delineation compared to using HCT/PETHCT.

Funding Statement

This study was financially supported by China Medical University. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.


1. Ashamalla H, Rafla S, Parikh K, Mokhtar B, Goswami G, et al. (2005) The contribution of integrated PET/CT to the evolving definition of treatment volumes in radiation treatment planning in lung cancer. Int. J. Radiat. Oncol. Biol. Phys. 63(4): 1016–23 [PubMed]
2. Kao CH, Hsieh TC, Yu CY, Yen KY, Yang SN, et al. (2010) 18F-FDG PET/CT-based gross tumor volume definition for radiotherapy in head and neck cancer: a correlation study between suitable uptake value threshold and tumor parameters. Radiat. Oncol. 5: 76 [PMC free article] [PubMed]
3. Aristophanous M, Berbeco RI, Killoran JH, Yap JT, Sher DJ, et al. (2012) Clinical utility of 4D FDG-PET/CT scans in radiation treatment planning. Int. J. Radiat. Oncol. Biol. Phys. 82(1): 99–105 [PubMed]
4. Steenbakkers RJ, Duppen JC, Fitton I, Deurloo KE, Zijp LJ, et al. (2006) Reduction of observer variation using matched CT-PET for lung cancer delineation: a three dimensional analysis. Int. J. Radiat. Oncol. Biol. Phys. 64: 435–48 [PubMed]
5. van Baardwijk A, Bosmans G, Boersma L, Buijsen J, Wanders S, et al. (2007) PET-CT–based auto-contouring in non–small-cell lung cancer correlates with pathology and reduces interobserver variability in the delineation of the primary tumour and involved nodal volumes. Int. J. Radiat. Oncol. Biol. Phys. 68: 771–8 [PubMed]
6. Bradley J, Thorstad WL, Mutic S, Miller TR, Dehdashti F, et al. (2004) Impact of FDG-PET on radiation therapy volume delineation in non-small cell lung cancer. Int. J. Radiat. Oncol. Biol. Phys. 59: 78–86 [PubMed]
7. Ekberg L, Holmberg O, Wittgren L, Bjelkengren G, Landberg T (1998) What margins should be added to the clinical target volume in radiotherapy treatment planning for lung cancer? Radiother. Oncol. 48(1): 71–7 [PubMed]
8. Liu HH, Balter P, Tutt T, Choi B, Zhang J, et al. (2007) Assessing respiration-induced tumor motion and internal target volume using four-dimensional computed tomography for radiotherapy of lung cancer. Int. J. Radiat. Oncol. Biol. Phys. 68(2): 531–40 [PubMed]
9. Britton KR, Starkschall G, Tucker SL, Pan T, Nelson C, et al. (2007) Assessment of gross tumor volume regression and motion changes during radiotherapy for non-small-cell lung cancer as measured by four-dimensional computed tomography. Int. J. Radiat. Oncol. Biol. Phys. 68(4): 1036–46 [PubMed]
10. Nehmeh SA, Erdi YE, Ling CC, Rosenzweig KE, Schoder H, et al. (2002) Effect of respiratory gating on quantifying PET images of lung cancer. J. Nucl. Med. 43(7): 876–81 [PubMed]
11. Pan T, Mawlawi O, Nehmeh SA, Erdi YE, Luo D, et al. (2005) Attenuation correction of PET images with respiration-averaged CT images in PET/CT. J. Nucl. Med. 46: 1481–87 [PubMed]
12. Huang T, Mok G, Wang S, Wu T, Zhang G (2011) Attenuation correction of PET images with interpolated average CT for thoracic tumors. Physics in Medicine and Biology. 56: 2559–2567 [PubMed]
13. Greta S. P Mok, Tao Sun, Tzung-Chi Huang, Mang I Vai (2013) Interpolated Average CT for Attenuation Correction in PET – A Simulation Study. IEEE Transactions on Biomedical Engineering. DOI:10.1109/TBME.2013.2245132 [PubMed]
14. Horn B K P, Schunck B G (1981) Determining optical flow. Artif. Intell. 17: 185–203
15. Studholme C, Hill D L G, Hawkes DJ (1999) An overlap invariant entropy measure of 3D medical image alignment. Pattern Recognit. 32: 71–86
16. Hatt M, Maitre AL, Wallach D (2012) Comparison of different methods of incorporating respiratory motion for lung cancer tumor volume delineation on PET images: a simulation study. Phys. Med. Biol. 57: 7409–7430 [PubMed]

Articles from PLoS ONE are provided here courtesy of Public Library of Science