- Open Access
A comparison of the different 3D CT scanning modes on the GTV delineation for the solitary pulmonary lesion
Radiation Oncology volume 9, Article number: 211 (2014)
To investigate the impacts of the different three-dimensional CT (3DCT) scanning modes on the GTV delineation for solitary pulmonary lesion (SPL) based on four-dimensional CT (4DCT), and to evaluate the feasibility of using the spiral CT scan in CT simulation.
Materials and methods
Twenty-one patients with SPL underwent axial CT scan, spiral CT scan and 4DCT simulation scan during free-breathing, respectively. The same clinical radiation oncologist delineated the gross tumor volume (GTV) under the same CT window setting. GTVA and GTVS were created from the axial and spiral images, respectively. ITVMIP was created from the maximum intensity projection (MIP) reconstructed 4D images. The target volumes and position between GTVA, GTVS and ITVMIP were compared. The matching index (MI) between GTVA and GTVS, and the correlation between MI and GTVS were evaluated.
ITVMIP was significantly larger than GTVA and GTVS (p s = 0.000). The ratios of ITVMIP to GTVA and GTVS were 1.57 ± 0.54 and 1.66 ± 0.61, respectively. There was no significant difference between GTVA and GTVS(p = 0.16). A comparison of the centroidal positions in x, y, and z directions for GTVA, GTVS and GTV4Dmip showed no significant difference (px = 0.17, py = 0.40, p z = 0.48). Additionally, there was no difference between distances from the centroidal positions of GTVA and GTVS to the origin of coordinates (p = 0.51). MI between GTVA and GTVS was 0.41 ± 0.24 (range 0–0.89), and it was positively correlated with the tumor volume (r = 0.64, p = 0.002).
There was no impact on the volume or centroidal position of GTV by the axial scan or spiral scan in 3DCT simulation for SPL. MI between GTVA and GTVS was small. A positively correlation was found between MI and GTVS. Relative to axial scanning mode, spiral CT scan was more timesaving and more efficient, it was feasible in 3DCT simulation for SPL.
The aim of radiation treatment is to increase the radiotherapy effects on the tumor while avoiding excessive toxicity [1, 2]. Computed tomography (CT) simulation and accurately determining the target margin were prerequisite for the successful oncology radiotherapy. Three-dimensional CT (3DCT) and four-dimensional CT (4DCT) simulation techniques are both in used, and 4DCT simulation is currently the leading simulation technology. It has enabled CT data acquisition to be correlated to the respiratory cycle(s), and allows a series of 3DCT data sets to be acquired at a number of points in a patient’s breathing cycle offering visualization of tumor motion on a patient by basis [3–5]. However, the 3DCT simulation is still widely used for 3D conformal radiotherapy (3D-CRT) at present. There were two scanning modes in 3DCT simulation, which were axial and spiral scanning mode. Because the speed of tube rotation in axial mode was slower than that in spiral mode, we often treat head and neck tumors that are not affected by respiratory motion with spiral scan . For thoracic and abdominal tumor, respiratory motion and heartbeat are important factors to determine the gross tumor volume (GTV). The slow speed scanning mode could collect more respiratory motion information, so we usually adopted axial scanning mode for simulation during free-breathing . But whether the axial scan and spiral scan that are usually adopted would influence the GTV of the lung tumor in motion during free-breathing was unclear. Basing on 4DCT technique, this study evaluates how the axial and spiral scanning modes affect the position, volume and spatial match relationships on the GTV delineation for the solitary pulmonary lesion (SPL) in 3DCT simulation.
Materials and methods
From March 2010 to September 2012, 21 patients with pathologically proven peripheral lung cancer or pulmonary metastasis who were treated with stereotactic body radiotherapy (SBRT) at Shandong Tumor Hospital were included in this study. Institutional Review Board of the hospital approval and informed consent were obtained for the present study. All patients had SPL with no adhesion between the tumor edge and pleura. And all of them could breathe freely and cooperate with CT scanning. The subjects consisted of 13 males and 8 females. The age of the patients ranged from 38 to 78 years, with a median age of 59 years. The lesion volumes ranged from 1.45 cm3 to 35.41 cm3, located in left lung in 9 cases and right lung in 12 cases, respectively.
CT simulation and image acquisition
The patients were immobilized in the supine position with their arms raised above head using vacuum bags. First, the axial scan was performed during free-breathing using Philip Brilliance big bore CT scanner (Philips Medical Systems, Inc., Cleveland, OH, USA). And then patients were conducted spiral scan and 4D scan consecutively, the image layer thickness was set to 3 mm. The axial scan cycle (scan + in-couch time) was 2.8 s. The pitch was 0.938 under spiral scanning. 4D technique was a kind of slow spiral scan mode, which would adjust the pitch according to respiratory rate in 4D scanning (pitch: 0.09 ~ 0.15) . The respiratory signal was recorded with the Real-Time Position Management (RPM, Varian Medical Systems, Palo Alto, CA, USA) by tracking the trajectory of the infrared markers placed on the patient’s abdomen . The respiratory signal was transferred to the CT scanning device, where the respiratory phase was matched with the acquired image(s). The 4D volume data were reconstructed into ten sequence images of different respiration phase (0%, 10%, 20% …90%) according to the respiration signal acquired by the RPM system, the 0% phase corresponded to the end-inhalation and the 50% phase corresponded to the end-exhalation. The ten sequence images were reconstructed with the same thickness of 3 mm.
Delineation, comparison and matching of GTV
The images obtained by three scanning modes were processed by Tumor Loc Software under the same CT window setting (W: 1600 C:-600). GTVA and GTVS were created from the two sets of CT images obtained by axial and spiral scanning modes, respectively (Figure 1A and B). The 3D coordinates of GTVA were Ax, Ay and Az. The coordinates of GTVS were Sx, Sy and Sz. calculating the distance from the centroidal position of GTVA and GTVS to the origin of coordinates, respectively. Tumor Loc software was used to retrospectively create the maximum intensity projection (MIP), and the ITVMIP was created from the MIP images  (Figure 1C). Using the Eclipse 8.6 treatment planning system, by which would conduct target registration on GTVA and GTVS, then calculate the matching index (MI) between GTVA and GTVS.MI = (GTVA∩GTVS)/(GTVA∪GTVS) (∩and ∪ were intersection and union, respectively) .
All statistical analyses were performed using the SPSS17.0 package (SPSS Inc, Chicago, IL, USA). The Wilcoxon test was used to compare the tumor volumes and the isocenter coordinates (x, y, z axial) gained by different scanning modes. A Pearson correlation analysis was utilized to study the relationship between MI and GTVS. Differences were considered to be significant if the p-value was less than 0.05.
The average ± standard deviation of GTVA and GTVS volumes were 9.94 ± 9.74 cm3 (range 1.72 ~ 36.56 cm3) and 10.14 ± 10.25 cm3 (range 1.45 ~ 35.41 cm3), respectively. There was no statistically significant difference between the volume of GTVA and GTVS (p = 0.16) (Table 1). The volume of ITVMIP was 14.33 ± 12.47 cm3 (rang 2.93 ~ 41.16 cm3). On pair wise analysis using Wilcoxon rank test, ITVMIP was significantly larger than GTVA and GTVS (p s = 0.000) (Table 1), and the ratios of ITVMIP to GTVA and GTVS were 1.57 ± 0.54 and 1.66 ± 0.61, respectively.
The center of tumor coordinates for GTVA, GTVS and ITVMIP were listed in Table 2 and the averages ± standard deviations were showed in Table 3. No differences were determined in x, y and z directions (px = 0.17, py = 0.40, p z = 0.48). The distances from the center of GTVA and GTVS to the origin of coordinate were 36.94 ± 24.64 cm and 36.88 ± 24.39 cm, respectively. There are no difference between the two distances (p = 0.51).
MI between GTVA and GTVS was 0.41 ± 0.24; the variation range was from 0 to 0.89. The correlation between the MI and GTVS was illustrated in Figure 2. There was a positive correlation between MI and GTVS. The correlation coefficients was (r = 0.64 p = 0.002).
SBRT is a safe and effective alternative therapy for patients with the inoperable lung cancer or metastatic tumor . Physiologic motions such as respiratory motion and heart beat would lead to tumor displacement and change in shape, so delineation of the GTV was a critical component in the radiotherapy process for thoracic tumor . Investigators have determined the target margin of thoracic tumor using the technology of active breathing control and 4DCT, which could reduce the interference on radiotherapy and increase the irradiation dose in target margin [13–16].
4DCT allowed for the clear delineation of the anatomical structures during free- breathing. This technique is widely used in evaluating the displacement of tumor and organ for thoracic and abdominal cancers. The ten sequence images reconstructed by 4D volume data could demonstrate movements of the tumor dynamically, and the ITVMIP represented the volume where tumor was present at any time within the respiratory cycle, which contained motion information of the tumor [9, 17–19]. The axial scanning mode was scanning alternates with the movement of scanner couch, that was, “scanning-in couch-scanning” pattern. As a result, the simulation time was long. Taking the ordinary 16 layers CT for example, the X -ray tube needs 1 s to rotate 360 degree and the scan cycle was 2.8 s by the axial mode, the scanning range was 24 mm (with 16 × 1.5 mm detector). The spiral mode was continuous scanning with CT scanner couch moving at a speed of 30 mm/s. There was no interval during the scanning process. It takes only 0.75 s to complete the scan rang of 24 mm in superior-inferior direction. In this study, we recorded and compared the scanning time for the same range by axial and spiral mode, respectively. The result revealed that axial mode took a longer time to complete the same scan than spiral mode, TA = 28.93 ± 4.63 s, TS = 8.81 ± 1.41 s (p = 0.000). Furthermore, the spiral mode could reduce motion artifacts, and increase the accuracy of tumor shape and position. Though the high speed of spiral mode, the scanning time for region of interest (ROI) was short for the small lesions and GTV contains little amount of respiratory motion information. The respiratory rate of the patients enrolled in the study was 16.00 ± 2.81breaths/min, and the respiratory cycle was 3.87 ± 0.73 s. The axial and spiral scans often complete ROI within 1 s. The key point was the ratio of the scanning time for ROI to the respiratory cycle. The higher the ratio was, the more motion information was included in the images. Among our subjects, there were 14 patients whose GTVA was larger than GTVS, while 7 smaller. There was no statistically significant difference in volume between the two scanning modes. We analyzed the tumor motion information included in GTVA and GTVS by comparing the volumes of ITVMIP with GTVA and GTVS, respectively. The ratio between ITVMIP and GTVA was 1.57 ± 0.54 (p = 0.000). It suggested that the tumor motion information included in GTVA was limited. The result was similar to the reports of Nakamura M [20, 21]. There was statistically significant difference in volume between ITVMIP and GTVS (p = 0.000). The ratio between ITVMIP and GTVS was 1.66 ± 0.61, which suggested that the motion information included in GTVS by spiral mode was limited too. The result of this study revealed that there was no difference in the volume of GTV by axial or spiral mode in 3DCT simulation. The tumor was moving when the axial and spiral scanning was performed during free-breathing. The GTV3D (GTVA and GTVS) created from the axial and spiral images included little respiration-induced motion information. It only represented the instant shape and volume that comes randomly by axial or spiral mode during the respiratory cycle.
ITVMIP was significantly larger than GTVA and GTVS. The differences in volume between GTV3D and ITVMIP were caused by the tumor motion in three dimensions. The motion information was included in ITVMIP. The motion amplitude in superior-inferior direction was larger than those in lateral and anterior-posterior directions. The superior-inferior motion of GTV centroid in lower lobe was larger than the in upper lobe [22, 23]. So the internal target volume (ITV) was generated from GTV3D by adding different margins in x, y and z directions. The compensated volume was still larger than ITVMIP created by 4D images.
The axial and spiral scan were all performed during free-breathing. The centroid coordinates of GTVA and GTVS were instant location during the respiratory cycle. There was no statistically significant difference in 3D coordinates of GTVA and GTVS, nor the distances from the centroidal positions to origin of coordinates. What’s more, there was no difference in the centroidal position between GTV3D and ITVMIP (px = 0.17, py = 0.40, p z = 0.48). This suggests that there was no obvious impact on for the tumors’ central position by the three scanning modes during free-breathing.
Although there was no difference in volume and the centroidal coordinates between the axial and spiral modes during free-breathing, there was still respiration-induced change in shape of tumor [24, 25]. We conducted target registration on GTVA and GTVS, MI range from 0 to 0.89. If the two target regions coincide completely, MI was 1, otherwise, that was 0. There was one case with no coincidence in our data, which suggested an obvious tumor displacement, and/or the tumor shapes varied greatly. Furthermore, there was a linear correlation between MI and GTVS (r = 0.64, p = 0.002). The small volume of GTV was one of the reasons for MI declination. The variation of MI was large in individuals. The low MI didn’t mean increase the margins, and different MI didn’t mean asymmetric margin either, because MI wasn’t determined by the displacement and the change in shape, it was affected by the tumor volume also.
In summary, our study found no impact on the volume and centroidal position of GTV by the axial or spiral scan in 3DCT simulation. There were no difference in the centroidal positions between 3DCT and 4DCT as well. The tumor motion information included in GTVA and GTVS was both limited. A linear correlation was found between MI and GTVS. Relative to axial scanning mode, spiral scan was more timesaving, reduced some motion artifacts and increased the accuracy of GTV due to its high scanning speed. So the spiral mode was a reasonable option to replace axial scan in 3DCT simulation for SPL.
Solitary pulmonary lesion
Three-dimensional conformal radiotherapy
Gross tumor volume
Gross tumor volume gained by axial scan
Gross tumor volume gained by spiral scan
Internal target volume gained by 4D maximum intensity projection
Maximum intensity projection
Stereotactic body radiotherapy
Real-Time Position Management.
Time of completing the scan by axial mode
Time of completing the scan by spiral mode
Region of interest
GTVA and GTVS
Internal target volume.
Nakayama H, Satoh H, Kurishima K, Ishikawa H, Tokuuye K: High-does conformal radiotherapy for patients with stage III non-small-cell lung carcinoma. Int J Radiat Oncol Biol Phys 2010, 78: 645-650. 10.1016/j.ijrobp.2009.08.037
Bradley JD, Nofal AN, El Naqa IM, Lu W, Liu J, Hubenschmidt J, Low DA, Drzymala RE, Khullar D: Comparison of helical, maximum intensityprojection (MIP), and averaged intensity (AI) 4D CT imaging for stereotactic body radiation therapy (SBRT) planning in lung cancer. Radio Oncol 2006, 81: 264-268. 10.1016/j.radonc.2006.10.009
Pan T, Lee TY, Rietzel E, Chen GT: 4D-CT imaging of a volume influenced by respiratory motion on multi-slice CT. Med Phys 2004, 31: 333-340. 10.1118/1.1639993
Vedam SS, Keall PJ, Kini VR, Mostafavi H, Shukla HP, Mohan R: Acquiring a four dimensional computed tomography dataset using an external respiratory signal. Phys Med Biol 2003, 48: 45-62. 10.1088/0031-9155/48/1/304
Ford EC, Mageras GS, Yorke E, Ling CC: Respiration-correlated spiral CT: a method of measuring respiratory-induced anatomic motion for radiation treatment planning. Med Phys 2003, 30: 88-97. 10.1118/1.1531177
Xianfu L, Bangxian T, Mi L, Yujun L, Yeqin Z, Jin H, Daiyuan M: Impact of CT scanning condition on the precision and accuracy of target position. China J Radiat Oncol 2009, 18: 142-145.
Wei H, Zheng F, Min F, Tonghai L, Heyi G, Baosheng L: Comparison of three CT scan methods used in precise radiotherapy of non-small cell lung cancer. China J Radiol Med Prot 2009, 29: 65-67.
Dunn L, Kron T, Tayor ML, Callahan J, Franich RD: A phantom for testing of 4D-CT for radiotherapy of small lesions. Med. Phys 2012, 39: 5372-5383. 10.1118/1.4742053
Cole AJ, O’Hare JM, McMahon SJ, McGarry CK, Butterworth KT, McAleese J, Jain S, Hounsell AR, Prise KM, Hanna GG, O’Sullivan JM: Investigating the potential impact of four-dimensional computed tomggraphy (4DCT) on toxicity, outcomes and dose escalation for radical lung cancer radiotherapy. Clin Oncol 2014, 26: 142-150. 10.1016/j.clon.2013.11.024
Underberg RW, Lagenwaard FJ, Slotman BJ, Cuijpers JP, Senan S: Use of maximum intensity projections (MIP) for target volume generation in 4DCT scans for lung cancer. Int J Radiat Oncol Biol Phy 2005, 63: 253-260. 10.1016/j.ijrobp.2005.05.045
Ezhil M, Vedam S, Balter P, Choi B, Mirkovic D, Starkschall G, Chang JY: Determination of patient-specific internal gross tumor volumes for lung cancer using four-dimensional computed tomography. Radia Oncol 2009, 4: 1-14. 10.1186/1748-717X-4-1
McGarry RC, Papies L, Williams R, Williams M, Whitford T, Timmerman RD: Stereotactic body radiation therapy of early-stage non-small-cell lung carcinoma: phase I study. Int J Radiat Oncol Biol Phys 2005, 63: 1010-1015. 10.1016/j.ijrobp.2005.03.073
Bosmans G, Van Baardwijk A, Dekker A, Ollers M, Boersma L, Minken A, Lambin P, De Ruysscher D: Intra-patient variability of rumor volume and tumor motion during conventionally fractionated radiotherapy for locally advanced on-small-cell lung cancer:a prospective clinical study. Int J Radiat Oncol Biol Phys 2006, 66: 748-753. 10.1016/j.ijrobp.2006.05.022
Yonghua Y, Yufen W, Shoufang G, Renben W, Wenwu L, Baosheng L: The influence on the position of peripheral lung cancer by the respiratory movement and the establishment of the mathematical model. China J Radiat Oncol 2004, 13: 83-85.
Cheung PC, Sixel KE, Tirona R, Along YC: Reproducibility of lung tumor position and reduction of lung mass within the planning target volume using active breathing control (ABC). Int J Radiat Oncol Biol Phys 2003, 57: 1437-1442. 10.1016/j.ijrobp.2003.08.006
Wang L, Hayes S, Paskalev K, Jin L, Buyyounouski MK, Ma CC, Feigenberg S: Dosimetric comparison of stereotactic body radiotherapy using 4D CT and multiphase CT images for treatment planning of lung cancer:eevaluation of the impact on daily dose coverage. Radiother Oncol 2009, 91: 314-324. 10.1016/j.radonc.2008.11.018
Low D: 4D imaging and 4D radiation therapy: a new era of therapy design and delivery. Front Radiat Ther Oncol 2011, 43: 99-117.
Li FX, Li JB, Zhang YJ, Liu TH, Tian SY, Xu M, Shang DP, Ma CS: Comparison of the planning target volume based on three dimensional CT and four dimensional CT images of non-small-cell lung cancer. Radio Oncol 2011, 99: 176-180. 10.1016/j.radonc.2011.03.015
Yu ZH, Lin SH, Balter P, Zhang L, Dong L: A comparison of tumor motion characteristics between early stage and locally advanced stage lung cancers. Radiother Oncol 2012, 104: 33-38. 10.1016/j.radonc.2012.04.010
Dongping S, Minghuan L, Jianbin L, Yong Y, Jinming Y, Jun D: The application of four-dimensional CT technique in the planning target volume of the solitary pulmonary lesion. China J Radiat Oncol 2011, 20: 417-419.
Nakamura M, Narita Y, Matsuo Y, Narabayashi M, Nakata M, Yano S, Miyabe Y, Matsugi K, Sawada A, Norihisa Y, Mizowaki T: Geometrical differences in target volume between slow CT and 4D CT imaging in stereotactic body radiotherapy for lung tumors in the upper and middle lobe. Med Phy 2008, 35: 4142-4148. 10.1118/1.2968096
Hof H, Rhein B, Haering P, Kopp-Schneider A, Debus J, Herfarth K: 4D CT based target volume definition in stereoactic radiotherapy of lung tumors:comparison with a convertional technique using individual margins. Radiother Oncol 2009, 93: 419-423. 10.1016/j.radonc.2009.08.040
Weiss E, Wijesooriya K, Dill SV, Keall PJ: Tumor and normal tissue motion in the thorax during respiration: analysis of volumetric and positional variations using 4DCT. Int J Radiat Oncol Biol Phys 2007, 67: 296-307. 10.1016/j.ijrobp.2006.09.009
Adamson J, Zhuang T, Yin FF: Contour based respiratory motion analysis for free breathing CT. Comput Biol Med 2011, 41: 908-915. 10.1016/j.compbiomed.2011.08.002
Hallman JL, Mori S, Sharp GC, Lu HM, Hong TS, Chen GT: A Four-dimensional computed tomography analysis of multiorgan abdominal motion. Int J Radiat Oncol Biol Phys 2012, 83: 435-441. 10.1016/j.ijrobp.2011.06.1970
This study was supported in part by Shandong Province Science and Technology Development Project (2012YD18089). Many colleagues helped to collect and analyze data. Their support and helps should be appreciated.
The authors declare that they have no competing interests.
D-PS delineated GTV in different series CT images and performed the statistical analysis. CL collected the CT images. YY conceived of the study, and helped to draft the manuscript. All authors read and approved the final manuscript.