Reproducibility of target volumes generated using uncoached 4-dimensional CT scans for peripheral lung cancer
© van der Geld et al; licensee BioMed Central Ltd. 2006
Received: 19 September 2006
Accepted: 01 November 2006
Published: 01 November 2006
4-dimensional CT (4DCT) scans are increasingly used to account for mobility during radiotherapy planning. As variations in respiratory patterns can alter observed motion, with consequent changes in the generated target volumes, we evaluated the reproducibility of 4D target volumes generated during repeat uncoached quiet respiration.
A retrospective analysis was performed on two successive scans (4DCT1 and 4DCT2) generated at the same scanning session for 26 patients with peripheral lung cancer treated with stereotactic radiotherapy (SRT). The volume and position of planning target volumes (PTV4DCT1 and PTV4DCT2) contoured on both scans were compared, and a dosimetric analysis performed. A SRT plan optimized for each PTV was sequentially applied to the other PTV, and coverage by the 80% isodose was evaluated. Color intensity projections (CIP) were used to evaluate regions of underdosage.
No significant volumetric differences were observed between the two PTVs (t-Test p = 0.60). The average displacement of the center of mass between corresponding PTVs was 1.4 ± 1.0 mm, but differences in position were 2.0 mm or greater in 5 cases (19%). Coverage of both PTVs by the 80% prescription isodose exceeded 90% for all but one patient. For the latter, the prescription isodose covered only 82.5% of PTV4DCT1. CIP analysis revealed that the region of underdosage was an end-inspiratory position occupied by the tumor for only 10–20% of the respiratory cycle.
In nearly all patients with stage I lung cancer, the PTV derived from a single uncoached 4DCT achieves dosimetric coverage that is similar to that achieved using two such consecutive scans.
Organ motion in the thorax confounds imaging, radiotherapy treatment planning and treatment delivery for patients with lung cancer . When free breathing CT scans are used for planning, variable geometric errors can arise with respect to the position, shape, and the volume of the gross tumor volume (GTV) and mobile organs [2–4]. With the use of 4-dimensional CT (4DCT), or respiration-correlated CT scans, volumetric image data can be acquired at many different respiratory states, and the time-resolved 3D CT imaging allows organ motion to be characterized for treatment planning [5–8]. A common technique for 4DCT acquisition involves axial cine acquisition with simultaneous recording of an external respiratory signal with over-sampled CT images, which are then sorted into several bins based on the phase information obtained from the respiratory signal . 4DCT was found to be superior to multiple conventional CT scans for generating internal target volumes (ITV) for stereotactic radiotherapy (SRT) of stage I non-small cell lung cancer (NSCLC) .
However, the use of 4DCT scans does not remove all mobility-induced artefacts. A study using phantoms found that gated 4DCT acquisition reduced, but did not eliminate volumetric errors when imaging moving objects . Another study reported residual volumetric differences between a static spherical object scanned helically or with 4DCT while moving, to be within ± 5% . In patients, additional problems can arise if the respiratory pattern is not sufficiently reproducible. This results in missing images at a specific respiratory phase for some table indices, leading to gaps in resorted 4DCT volumes [10, 11]. Breathing patterns in patients with lung cancer can be complex and exhibit considerable diversity within the same patient . Routine breath coaching has been proposed as a method to improve regular respiratory motion during radiotherapy . However, delivery of fractionated SRT can take up to 40 minutes, which makes coached respiration not feasible in patients who often have severely compromised lung function. As respiration-gated delivery at our center is only used for stage III tumors, and not for SRT for small peripheral tumors, we derived an ITV that was based upon all mobility observed on two consecutive uncoached 4DCT scans.
The aim of the present study was to assess the reproducibility of the two consecutive PTVs based upon uncoached 4DCT scans that were generated at a single imaging session.
The successive planning 4DCT scans of 26 consecutive patients with stage I NSCLC who were treated with SRT were retrospectively analyzed.
4DCT scanning procedure
The 4DCT scanning approach used at our center has been reported in detail previously [6–8]. Briefly, patients are immobilized in the supine position with both arms positioned above the head on an adjustable arm support. The Real-Time Position Management system (RPM, Varian medical systems, Palo Alto, CA) uses a lightweight block containing two reflective markers placed on the upper abdomen to record the respiratory signal. 4DCT scans are performed during uncoached quiet respiration on a 16 slice CT scanner (Lightspeed 16 GE Medical Systems, Waukesha, WI). CT scanning commences when the RPM system shows a regular respiratory pattern, and multiple respiratory cycles are observed in order to determine the adequate cine scan duration parameter.
For thoracic scans, 8 contiguous slices of 2.5 mm are generated for a 2 cm total longitudinal coverage per gantry rotation with the scanner operated in axial cine mode. Other scanning parameters include 140 KV, 95 mA and tube rotation set closest to 1/10 of the average breathing cycle time to allow high temporal and spatial resolution. With the scanner couch in static mode, data is acquired for at least the duration of one full respiratory cycle with the addition of one second, after which the couch advances to the next position. Data acquisition ceases during the couch movement, and a full 4DCT scanning procedure of the thorax generally takes about 90 seconds. The radiation exposure from 4DCT acquisition is approximately 6 times the dose from a single conventional helical CT scan . As variations in breathing patterns during uncoached CT scans can lead to variations in the defined target volume, our routine imaging protocol for SRT planning requires that a second 4DCT scan (limited to the tumor region) be performed directly after the first. The additional information was considered crucial as doses of 60 Gy (to the 80% isodose) were delivered in 3 fractions with no intra-fractional imaging.
Advantage4D software (GE Medical Systems, Waukesha, WI) is used to retrospectively sort the acquired images into temporally coherent volumetric datasets. Each reconstructed image is assigned to one of 10 respiratory phases based on the temporal correlation between surface motion of the lightweight block and data acquisition. For each couch position, the reconstructed image with the respiratory phase nearest to the requested phase is selected to form the 3D dataset for a given phase of the breathing cycle. The user can define tolerance values for phase selection, and such phase tolerance is necessary in order to obtain complete volumetric information as images at different couch positions are reconstructed at different respiratory phase samples .
Generating Internal Target Volumes (ITVs)
The clinical SRT delivery for each patient was based on non-gated treatment delivery to an ITV that was derived from the summation of contouring on both 4DCT scans. ITVs were generated using maximum intensity projection (MIP) post-processing of data from each 4DCT scan in Advantage4D software. The MIP images reflect the maximum value of any given pixel intensity over all 10 phases of the 4DCT scan, and contoured MIP images have been shown to correspond well with ITVs derived from manual contouring of the corresponding GTVs in all 10 phases of a 4DCT .
A single end-expiratory phase of the 4DCT scan and both MIP datasets were imported into the BrainLab stereotactic planning system (Brainscan v 5.2, BrainLab AG, Heimstetten, Germany), and both MIP scans were co-registered with the single phase bin scan using the fusion software tool in Brainscan. If necessary, adjustments were performed to ensure a proper match of the spinal column. The ITV was derived from contouring the tumor on both MIP scans. A contour encompassing both ITV4DCT's was used to derive the ITV for the clinical plan. An isotropic ITV-PTV margin of 3 mm was added in order to account for residual patient setup inaccuracies on the Novalis Exactrac (Brainscan v 5.2, BrainLab AG, Heimstetten, Germany), which uses an online patient setup verification system. No separate margins were used to account for microscopic tumor extension .
Retrospective comparison of PTVs
Two separate PTVs, namely PTV4DCT1 and PTV4DCT2, were generated and the volumes and location were compared. For the latter, the displacement of the center of mass of both PTVs in the X (medio-lateral), Y (anterior-posterior), and Z (superior-inferior) direction was assessed using the Brainlab coordinate system.
A dosimetric analysis was performed in order to evaluate the clinical relevance of the variations between PTVs, if present. The clinical treatment plan, which consisted of between 8–12 non-coplanar treatment beams, was sequentially modified by adjusting the beam portals to both PTV4DCT1 and PTV4DCT2. Stereotactic treatment planning was normalized to the 100% isodose at the isocenter. As for the clinical treatment plan, PTV coverage by the 80% isodose was the goal of treatment planning. Coverage of the PTV4DCT2 was evaluated using the plan developed for PTV4DCT1, followed by the reverse.
Color intensity projections (CIP)
In order to obtain temporal information in cases with significant underdosage of the PTV based on a single 4DCT, a composite color intensity projection (CIP) image of the 4DCT data was generated . Briefly, the CIP technique generates color images that represent a projection of the intensity changes in the component images into a single image. Pixels whose intensity changes over the phases of the 4DCT appear in the CIP image as colored pixels, and the hue encodes the percentage of time that tissue is present in each location. CIP was implemented as a plug-in tool to ImageJ that was written in-house. ImageJ (version 1.31v) is a Java based image-processing package, and is free public domain software .
The Student paired t-test for comparison of volume of the PTV4DCT1 and PTV4DCT2 was performed using Excel software (Microsoft, Redmond, WA).
Volumetric and spatial comparison of PTV4DCT1 and PTV4DCT2
PTV4DCT1 – PTV4DCT2 (cm3)
PTV4DCT1 – PTV4DCT2 (%)
3D center of mass displacement (mm)
Mean ± SD
1.2 ± 1.1 cm 3
0.7 ± 7.6%
1.4 ± 1.0 mm
Percentage of each PTV encompassed by the 80% prescription isodose in a plan optimized for the other PTV.
Coverage of PTV4DCT2 with planning based on PTV4DCT1
Coverage of PTV4DCT1 with planning based on PTV4DCT2
Mean ± SD
97.8 ± 2.3%
97.0 ± 3.8%
Individualized margins are required to account for mobility of lung tumors, as no reliable correlation exists between tumor location and mobility . With the use of 4DCT, spatial and temporal information on tumor and normal organ mobility is generated by synchronously recording the respiratory waveform during imaging. The 4DCT process includes two steps, the first of which is the acquisition of several CT images per couch position through the volume of interest, which is followed by resorting of the acquired images into spatiotemporally coherent data sets.
Various approaches for 4DCT acquisition have been reported, one of which is a widely applied technique in clinical practice that uses the axial cine mode [6, 8, 9]. This approach has been shown to be fairly accurate in phantom studies, but residual motion artifacts arise due to factors such as partial volume projection effects, phase tolerance values for resorting, and irregularities in the breathing pattern of patients . A recent report suggested that the 4DCT system used in our patients fails to generate optimal data in approximately 30% of patients due to image sorting errors . Reasons cited for the latter include the fact that processing of respiratory traces may not accurately reflect the respiratory phase, and that the RPM-software excludes tracing from periods of irregular respiration and Advantage4D disregards the corresponding images. However, with the generation of ITVs using MIP datasets generated from all respiratory phases, image sorting inaccuracies become irrelevant, and remaining differences between 4DCT scans are therefore due to either lacking images of extreme tumor positions as in the described case, or due to breathing irregularities during the CT scanning session.
Although the absolute underdosage in the most caudal part of the PTV for this fraction of radiation appeared to be important, analysis using the CIP technique revealed that this region of the tumor was present in the area of underdosage for only 10–20% of the time. For the actual stereotactic treatment of this patient, the variation in PTVs was of no consequence as routine clinical plans were based upon ITVs derived from both 4DCT scans. Furthermore, we perform repeat 4DCT scans during the course of SRT .
Our results on the reproducibility of uncoached 4DCT scans appear to be even more reassuring in the context of conventionally fractionated 3D conformal radiotherapy. The larger GTV to PTV margins, and certainly also the much lower conformity index of a 3–4 field treatment plan , indicates that the variations between ITVs from uncoached 4DCT scans would not have a great impact. Furthermore, previous analysis indicates that incorporating the extremes of respiration-induced tumor movement during conventionally fractionated radiotherapy leads to an overestimation of the PTV [21–24].
Respiratory coaching during 4DCT imaging and treatment delivery has been advocated for patients undergoing gated radiotherapy. This approach is, however, not feasible in patients undergoing SRT as the delivery of a single fraction of non-gated SRT can take more than 30–40 minutes. Furthermore, it is doubtful if patients with medically inoperable stage I NSCLC who generally have severely impaired lung function can tolerate prolonged periods of coached breathing. We will study the variations in respiratory patterns of patients during the delivery of uncoached SRT, and interrupting treatment delivery when amplitude changes in RPM exceed those observed during the 4DCT scan may be one approach to ensure target reproducibility.
In the majority of patients with stage I lung cancer, the PTV derived from a single uncoached 4DCT achieves similar dosimetric coverage as is achieved using two consecutive 4DCT scans. This indicates that the intra-fractional reproducibility of uncoached 4DCT scans in this patient group is sufficiently reliable in clinical plans based upon target volumes that incorporate all mobility.
- Senan S, Chapet O, Lagerwaard FJ, Ten Haken RK: Defining target volumes for non-small cell lung carcinoma. Semin Radiat Oncol 2004, 14: 308-314. 10.1016/j.semradonc.2004.07.004View ArticlePubMedGoogle Scholar
- Balter JM, Ten Haken RK, Lawrence TS, Lam KL, Robertson JM: Uncertainties in CT-based radiation therapy treatment planning associated with patient breathing. Int J Radiat Oncol Biol Phys 1996, 36: 167-174. 10.1016/S0360-3016(96)00275-1View ArticlePubMedGoogle Scholar
- Shen S, Duan J, Fiveash JB, Brezovich IA, Plant BA, Spencer SA, Popple RA, Pareek PN, Bonner JA: Validation of target volume and position in respiratory gated CT planning and treatment. Med Phys 2003, 30: 3196-3205. 10.1118/1.1626121View ArticlePubMedGoogle Scholar
- Keall PJ, Starkschall G, Shukla H, Forster KM, Ortiz V, Stevens CW, Vedam SS, George R, Guerrero T, Mohan R: Acquiring 4D thoracic CT scans using a multislice helical method. Phys Med Biol 2004, 49: 2053-2067. 10.1088/0031-9155/49/10/015View ArticlePubMedGoogle Scholar
- 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/304View ArticlePubMedGoogle Scholar
- Underberg RWM, Lagerwaard FJ, Cuijpers JP, Slotman BJ, van Sörnsen de Koste JR, Senan S: Four-dimensional CT scans for treatment planning in stereotactic radiotherapy for stage I lung cancer. Int J Radiat Oncol Biol Phys 2004, 60: 1283-1290. 10.1016/j.ijrobp.2004.07.665View ArticlePubMedGoogle Scholar
- Underberg RWM, Lagerwaard FJ, Slotman BJ, Cuijpers JP, Senan S: Benefit of respiration-gated stereotactic radiotherapy for stage I lung cancer: an analysis of 4DCT datasets. Int J Radiat Oncol Biol Phys 2005, 62: 554-560. 10.1016/j.ijrobp.2005.01.032View ArticlePubMedGoogle Scholar
- van Sörnsen de Koste JR, Senan S, Kleynen CE, Slotman BJ, Lagerwaard FJ: Renal mobility during uncoached quiet respiration: An analysis of 4DCT scans. Int J Radiat Oncol Biol Phys 2006, 64: 799-803. 10.1016/j.ijrobp.2005.09.012View ArticlePubMedGoogle Scholar
- 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.1639993View ArticlePubMedGoogle Scholar
- Rietzel E, Pan T, Chen GT: Four-dimensional computed tomography: image formation and clinical protocol. Med Phys 2005, 32: 874-889. 10.1118/1.1869852View ArticlePubMedGoogle Scholar
- Keall P: 4-dimensional computed tomography imaging and treatment planning. Semin Radiat Oncol 2004, 14: 81-90. 10.1053/j.semradonc.2003.10.006View ArticlePubMedGoogle Scholar
- Ozhasoglu C, Murhpy MJ: Issues in respiratory motion compensation during external-beam radiotherapy. Int J Radiat Oncol Biol Phys 2002, 52: 1389-1399. 10.1016/S0360-3016(01)02789-4View ArticlePubMedGoogle Scholar
- Underberg RW, Lagerwaard 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 Phys 2005, 63: 253-260. 10.1016/j.ijrobp.2005.05.045View ArticlePubMedGoogle Scholar
- Potters L, Steinberg M, Rose C, Timmerman R, Ryu S, Hevezi JM, Welsh J, Mehta M, Larson DA, Janjan NA: American Society for Therapeutic Radiology and Oncology and American College of Radiology practice guideline for the performance of stereotactic body radiation therapy. Int J Radiat Oncol Biol Phys 2004, 60: 1026-1032. 10.1016/j.ijrobp.2004.07.701View ArticlePubMedGoogle Scholar
- Cover KS, Lagerwaard FJ, Senan S: Color intensity projections: A rapid approach for evaluating four-dimensional CT scans in treatment planning. Int J Radiat Oncol Biol Phys 2006, 64: 954-961. 10.1016/j.ijrobp.2005.10.006View ArticlePubMedGoogle Scholar
- van Sörnsen de Koste JR, Lagerwaard FJ, Nijssen-Visser MR, Graveland WJ, Senan S: Tumor location cannot predict the mobility of lung tumors: a 3D analysis of data generated from multiple CT scans. Int J Radiat Oncol Biol Phys 2003, 56: 348-354. 10.1016/S0360-3016(02)04467-XView ArticlePubMedGoogle Scholar
- Rietzel E, Chen GT: Improving retrospective sorting of 4D computed tomography data. Med Phys 2006, 33: 377-379. 10.1118/1.2150780View ArticlePubMedGoogle Scholar
- Underberg RW, Lagerwaard FJ, van Tinteren H, Cuijpers JP, Slotman BJ, Senan S: Time trends in target volumes for stage I non-small-cell lung cancer after stereotactic radiotherapy. Int J Radiat Oncol Biol Phys 2006, 64: 1221-1228. 10.1016/j.ijrobp.2005.09.045View ArticlePubMedGoogle Scholar
- van Sörnsen de Koste JR, Lagerwaard FJ, Schucchard-Schipper RH, Nijssen-Visser MR, Voet PW, Oei SS, Senan S: Dosimetric consequences of tumor mobility in radiotherapy of stage I non-small cell lung cancer – an analysis of data generated using 'slow' CT scans. Radioth Oncol 2001, 61: 93-99. 10.1016/S0167-8140(01)00373-5View ArticleGoogle Scholar
- Stroom JC, de Boer HC, Huizenga H, Visser AG: Inclusion of geometrical uncertainties in radiotherapy treatment planning by means of coverage probability. Int J Radiat Oncol Biol Phys 1999, 43: 905-919. 10.1016/S0360-3016(98)00468-4View ArticlePubMedGoogle Scholar
- Engelsman M, Damen EM, De Jaeger K, van Ingen KM, Mijnheer BJ: The effect of breathing and set-up errors on the cumulative dose to a lung tumor. Int J Radiat Oncol Biol Phys 2001, 60: 95-105.Google Scholar
- van Herk M, Witte M, van der Geer J, Schneider C, Lebesque JV: Biologic and physical fractionation effects of random geometric errors. Int J Radiat Oncol Biol Phys 2003, 57: 1460-1471. 10.1016/j.ijrobp.2003.08.026View ArticlePubMedGoogle Scholar
- Engelsman M, Sharp GC, Bortfeld T, Onimaru R, Shirato H: How much margin reduction is possible through gating or breath hold? Phys Med Biol 2005, 50: 477-490. 10.1088/0031-9155/50/3/006View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.