Use of combined maximum and minimum intensity projections to determine internal target volume in 4-dimensional CT scans for hepatic malignancies
Radiation Oncologyvolume 7, Article number: 11 (2012)
To evaluate the accuracy of the combined maximum and minimum intensity projection-based internal target volume (ITV) delineation in 4-dimensional (4D) CT scans for liver malignancies.
4D CT with synchronized IV contrast data were acquired from 15 liver cancer patients (4 hepatocellular carcinomas; 11 hepatic metastases). We used five approaches to determine ITVs: (1). ITVAllPhases: contouring gross tumor volume (GTV) on each of 10 respiratory phases of 4D CT data set and combining these GTVs; (2). ITV2Phase: contouring GTV on CT of the peak inhale phase (0% phase) and the peak exhale phase (50%) and then combining the two; (3). ITVMIP: contouring GTV on MIP with modifications based on physician's visual verification of contours in each respiratory phase; (4). ITVMinIP: contouring GTV on MinIP with modification by physician; (5). ITV2M: combining ITVMIP and ITVMinIP. ITVAllPhases was taken as the reference ITV, and the metrics used for comparison were: matching index (MI), under- and over-estimated volume (Vunder and Vover).
4D CT images were successfully acquired from 15 patients and tumor margins were clearly discernable in all patients. There were 9 cases of low density and 6, mixed on CT images. After comparisons of metrics, the tool of ITV2M was the most appropriate to contour ITV for liver malignancies with the highest MI of 0.93 ± 0.04 and the lowest proportion of Vunder (0.07 ± 0.04). Moreover, tumor volume, target motion three-dimensionally and ratio of tumor vertical diameter over tumor motion magnitude in cranio-caudal direction did not significantly influence the values of MI and proportion of Vunder.
The tool of ITV2M is recommended as a reliable method for generating ITVs from 4D CT data sets in liver cancer.
Primary and metastatic hepatic malignancies are commonly treated by surgery, but radiation therapy is also one of options as non-surgical modalities. It has been demonstrated that radiation therapy is feasible and the outcomes are promising [1, 2]. However, due to respiration liver motion up to 3 cm  is one of obstacles to accurately localize the target. Moreover, respiratory-induced tumor motion is known to be anisotropic, thus individual determination of internal margin around gross tumor volume (GTV) is crucial to form an internal target volume (ITV), which can avoid both inadequate tumor coverage and unnecessary liver parenchymal irradiation for individual patient. Four-dimensional computed tomography (4D CT) is one of appropriate approaches to estimate and determine ITV for tumor with respiratory motion [4, 5].
4D CT has been widely used in lung cancer to determine ITV [6, 7]. Ideally, ITV should be delineated by manually contouring GTV in all 10 breath phases of a 4D scan image sets to form ITV, which is the most accurate tool to determine ITV, but it is a time-consuming and labor-intensive task. To reduce the workload of contouring multiple GTVs, one solution is to contour only two extreme phases at end-inhalation and end-exhalation and then to sum of the two becoming ITV [8, 9]; and the other is to use the post-processing tools of maximum intensity projection (MIP) and minimum intensity projection (MinIP) from 4D CT data sets to generate ITV. MIP-based ITV delineation is performed on a single 3-D CT data set, where each pixel in this set represents the brightest object encountered by corresponding voxels in all volumetric 4D CT data sets, for instance, MIP-based ITV delineation for lung cancer, which was recommended as a reliable tool and a good first estimation [10–12]. Conversely, MinIP-based ITV is on the CT set, where each pixel represents the lowest data value in the volumetric data .
MIP and MinIP methods seem not suitable for liver cancer because most tumors in the liver have similar attenuation to the normal liver parenchyma and therefore are not easily discernable. Contrast-enhanced CT scan has been routinely used for radiation oncologists to differentiate the tumor from normal tissues. It should be noted that regardless of with contrast enhancement or not, most liver tumors present inhomogeneous density, either because of the inherent nature of the tumor, such as the routes of blood supply, vascular volume and permeability, or because of areas of fluid, hemorrhage, and necrosis within tumors, or because of secondary change due to treatments, for example, iodine deposition as a result of transcatheter arterial chemoembolization (TACE). Recently, Beddar et al described a simple method for 4D CT acquisition by using synchronized intravenous contrast injection to improve the accuracy of liver tumor delineation. By this way most liver metastases and cholangiocarcinomas can be identified on image of portal venous phase, while HCC, most visible in the delayed phase .
In theory, for tumors with homogeneous hyperdensity or hypodensity comparing to the surrounding normal liver, MIP or MinIP projections should accordingly reflect the tumor trajectory across all time-resolved data sets. Visualization of tumor with mixed-density means the tumor border should be discernable; regardless which part of it is more hyperdense or hypodense than the adjacent liver parenchyma. Using MIP technique or MinIP technique only definitely misses the spatial information of the moving liver tumor. Therefore, our hypothesis is to combine MIP and MinIP, which may fit for the situation of mixed tumor density. Thus, the purpose of this study is to evaluate the feasibility and accuracy of the MIP and/or MinIP-based ITV delineation in 4DCT scans for liver tumors.
Patients who met the following criteria were qualified to this study: (1). Patient suffered from hepatic malignancies, primary or metastatic, and were planned to receive irradiation; (2). The margins of hepatic lesion were clear on intravenous contrast enhanced CT; (3). Patient's breath was regular after a training session; (4). Patient did not have the history of allege to contrast; and (5). Informed consent was obtained.
4D CT image acquisition
In order to enhance the visibility of tumors on 4D CT, 100 mL of contrast at a concentration of 300 mg I/mL was injected synchronously with 4D CT image acquisition. All patients were imaged during the portal venous phase. A time delay was programmed within 4D CT image acquisition protocol so that the start of contrast injection is initiated simultaneously with the start of the scanner's timer countdown. For those with liver metastases, the liver was scanned with a 45 s time delay, while for hepatocellular carcinoma (HCC) patients the time delay was 65 s. This CT scan protocol was proposed by Beddar et al .
All patients were immobilized using customized vacuum lock in supine position with arms placed on their forehead. A16-slice Brilliance Big Bore CT scanner (Philips Co.) was used for 4D CT image acquisition. The patient respiration was tracked using Real-Time Position Management (RPM) System (Varian Medical Systems, Palo Alto, CA). The region of interest usually comprises the area from 2 to 3 cm above diaphragm to iliac crest. A 4D CT scan is performed in cine mode with at least one complete breath cycle for each couch position. After scanning 4-dimensional images were binned based on the respiratory traces to become a complete image set, which covered each of 10 breathing phases. MIP and minimum intensity projections (MinIP) were then generated from the raw data set of 4D CT scans.
All reconstructed CT series were transferred to MIM software (Version 5.1, MIMvista Corp., Cleveland, OH). Window/level was adjusted to optimize the visual contrast between the tumor and normal parenchyma regions. All the contours were drawn by a single radiation oncologist (JL) and verified by a senior radiation oncologist (JDZ). We used 5 approaches to determine ITV, which were: (1). ITVAllPhases: contouring GTV on each of 10 respiratory phases of 4D CT data set and combining these GTVs; (2). ITV2Phase: contouring GTV on CT of the peak inhale phase (0% phase) and the peak exhale phases (50%) and then combining the two; (3). ITVMIP: contouring GTV on MIP with modifications based on physician's visual verification of contours in each respiratory phase; (4). ITVMinIP: contouring GTV on MinIP with modification by physician; (5). ITV2M: combining ITVMIP and ITVMinIP. Figure 1 illustrates different approaches in the determination of ITV for patient 13.
Evaluation of target motion
We used MIM software (version 5.1) to measure cranio-caudal, left-right and anterior-posterior movement of the tumor. Tumor motion was also expressed as a 3D vector, which is the quadratic mean of the motions in 3 orthogonal directions.
Data analysis and statistics
ITV2Phase, ITVMIP, ITVMinIP and ITV2M, as the tested ITVs, would be compared with the reference of ITVAllPhases, respectively. The metrics used for comparison were: (1) Matching index (MI), which was the percentage of overlapped volume in 3-dimensions between 2 ITVs. When 2 ITVs were totally overlapped, MI was 1.00, whereas when 2 ITVs were totally separate, MI was 0; and (2) under- and over-estimated volume (Vunder and Vover). A tested ITV (V test) was compared to the standard volume, ITVAllPhases. The formulas were: Vunder = VAllPhases\Vtest, and Vover = Vtest \ VAllPhases.
The definitions and calculations of those metrics were referred to Ezhil .
Paired sample t-test and Independent-Samples T test (SPSS v.13, SPSS Inc., Chicago, IL) were applied to compare the differences with p value of < 0.05 considered significant.
From August 2010 to February 2011, 15 eligible patients with liver cancer were enrolled in this study and underwent 4D CT simulation for irradiation treatment planning in our institution. Of these patients 11 had metastatic liver cancers and 4, HCC with mean lesion volume of 152 cm3 (range, 2 cm3 - 932 cm3). 4D CT images were successfully obtained from 15 patients and tumor margins were clearly discernable in all patients. There were 9 cases of low density and 6, mixed on CT images (Table 1).
The cranio-caudal motion of the target was predominant with a mean distance of 8.0 mm ± 3.3 mm, while the left-right and anterior-posterior motions were much less with mean values of 1.6 mm ± 0.9 mm and 3.2 mm ± 2.2 mm, respectively (Table 2).
Comparison of ITVs countered by different approaches
(1). ITVs volume
Table 3 shows all volumes of ITV. The mean volume of ITV2M was closest to that of ITVAllPhases, and then followed by ITV2Phase. The volume difference between ITV2M and ITV2Phase was statistically significant with a p value of 0.04. Taking ITVAllPhases as the reference, the mean ratios of the tested ITVs to ITVAllPhases were 88.9% ± 5.7%, 82.7% ± 12.6%, 82.5% ± 10.8% and 94.0% ± 3.6%, respectively for ITV2phase, ITVMinIP, ITVMIP and ITV2M.
As shown in Table 4, ITV2M was closest matched with ITVAllPhases with mean MI of 0.93 ± 0.04, and mean MIs were 0.89 ± 0.06, 0.82 ± 0.12 and 0.82 ± 0.10, respectively for ITV2Phase, ITVMinIP and ITVMIP. The differences of MI were statistically significant between ITV2M and ITV2Phase (p = 0.004), between ITV2M and ITVMinIP (p = 0.003), and between ITV2M and ITVMIP (p = 0.000). All the other comparisons between ITVs were not significant.
We further analyze the tumor characteristics, which would impact MI.
1). MI and tumor volume: As lesion of > 5 cm is not a candidate for stereotactic body radiation therapy in our practice, the patients were split into two groups: ≤ 65.4 cm3 (equivalent to the volume of a sphere with 5 cm in diameter) vs. > 65.4 cm3. MIs of ITV2M were 0.94 ± 0.02 for the former and 0.92 ± 0.04 for the latter (p = 0.260). Nevertheless, there was significant difference between tumor size and MIs of ITV2Phase, MI being 0.92 ± 0.03 for tumor of ≤ 65.4 cm3 and 0.85 ± 0.06 for tumor of > 65.4 cm3 (p = 0.015). For ITV2M as tumor volume increased MI did not change much with no correlation between them (R = 0.364, p = 0.182). However, for ITV2Phase, when tumor volume increased MI was significantly enhanced with positive correlation (R = 0.527, p = 0.044) (Figure 2).
2). MI and target motion three-dimensionally: It is recommended that for tumor motion of > 10 mm, we need to reduce the movement by breath control devices, such as Active Breathing Coordinator or RPM gating system. The patients were divided into two groups: 3D vector of ≤ 10 mm vs. > 10 mm. Mean MIs of ITV2M were 0.94 ± 0.01 and 0.92 ± 0.04 (p = 0.542), and mean of MIs of ITV2Phase, 0.88 ± 0.04 and 0.89 ± 0.06 (p = 0.756), respectively for the former and the latter. There was no significant correlation between the magnitude of target motion and MI.
3). MI and ratio of tumor vertical diameter over tumor motion magnitude in cranio-caudal direction: As shown in Figure 3, for ITV2M when vertical diameter over tumor motion magnitude in cranio-caudal direction increased MI was increased slightly, but with no significant correlation (R = 0.352, p = 0.198). However, for ITV2Phase there was a positive correlation between them (R = 0.535, p = 0.040).
4). MI and tumor density: For the tumor of low density or mixed density, there was no significant difference in MIs between them for both ITV2M and ITV2Phase, MIs of ITV2M being 0.93 ± 0.02 and 0.92 ± 0.05 (p = 0.676), and MIs of ITV2Phase, 0.88 ± 0.07 and 0.90 ± 0.05 (p = 0.702), respectively for low density tumor and the mixed. However, for low density tumor MI of ITVMinIP was better than that for mixed density tumor with MI of 0.90 ± 0.05 and 0.71 ± 0.11 (p = 0.001), respectively.
We also noticed that for low density tumors, which located within liver parenchyma and were not closed to adjacent organs, such as 5th, 6th, 8th and 12th patient (Table 1), when using ITVMinIP mean MI was 0.92 ± 0.04, while using ITV2M it was 0.93 ± 0.03 with no significant difference between ITVMinIP and ITV2M (p = 0.114).
(3). Proportion of Vunder
Figure 4 illustrates the proportions of Vunder in 15 patients. Compared to ITVAllPhases, the proportional Vunder of ITV2M was the lowest (0.07 ± 0.04) with the maximum of 0.17 among ITV2Phase, ITVMinIP, ITVMIP and ITV2M. While proportional Vunder of ITV2Phase was 0.11 ± 0.06 with the maximum of 0.26. The mean proportion of Vunder for ITV2M were significantly less than that for ITV2Phase (p = 0.001). However, ITVMinIP and ITVMIP underestimated larger volumes, the proportions of Vunder being 0.18 ± 0.12 and 0.18 ± 0.11, respectively.
The analyses of tumor characteristics, which would impact proportion of Vunder were as follows.
1). Proportion of V under and tumor volume: There were no correlations between the diameter of GTV and the proportion of Vunder for ITV2M, no matter the diameter over than or less than 5 cm. The proportions of Vunder were respectively 0.06 ± 0.02 and 0.08 ± 0.05 (p = 0.244). However, there was significant difference in proportion of Vunder between tumor size of ≤ 5 cm and > 5 cm with proportions of ITV2Phase being 0.08 ± 0.03 and 0.15 ± 0.07 (p = 0.018), respectively. For ITV2M as tumor size increased the proportions of Vunder did not change significantly with no correlations between them (R = -0. 340, p = 0.215). In contrast, for ITV2Phase there was negative correlation between GTV volume and the proportion of Vunder (R = -0.528, p = 0.043) (Figure 5).
2). Proportion of V under and target motion three-dimensionally: There was no strong correlation between 3 victor tumor motion of ≤ 10 mm and > 10 mm, and the proportions of Vunder were 0.06 ± 0.01 and 0.07 ± 0.04 (p = 0.480) for ITV2M, and 0.12 ± 0.04 and 0.10 ± 0.07 (p = 0.758) for ITV2phase.
3). Proportion of V under and ratio of tumor vertical diameter over tumor motion: For ITV2M there was no strong correlation between them (R = -0.344, p = 0.210). However, for ITV2phase as ratio of tumor vertical diameter over tumor motion increased the proportions of Vunder decreased significantly (R = -0.523, p = 0.040) (Figure 6).
4). Proportion of V under and tumor density: For both low density and mixed density tumors there was no significant difference in the underestimations, regardless ITV2M or ITV2Phase, For ITV2M the proportions of Vunder were 0.06 ± 0.03 and 0.07 ± 0.05 (p = 0.723), respectively for low density tumor and the mixed, and for ITV2Phase they were 0.12 ± 0.07 and 0.11 ± 0.05 (p = 0.680). However, when using ITVMinIP the underestimated proportion was 0.10 ± 0.06 for 9 low density tumors, but it was 0.29 ± 0.11 for 6 mixed density tumors (p = 0.001).
As for MI in 5th, 6th, 8th and 12th patient, when ITVMinIP was used the mean proportion of Vunder was 0.08 ± 0.04, while it was 0.07 ± 0.03 when ITV2M was used (p = 0.094).
Compared to ITVAllPhases, the percentages of Vover were all less than 1% for ITV2phase, ITVMinIP, ITVMIP and ITV2M.
MIP method as an image post-processing is based on more complex algorithms and can be used for generating three-dimensional vascular reconstructions [14, 15]. In lung cancer, MIPs are believed to be a reliable tool to generate ITVs from 4D CT data sets , however, it is mandatory to modify each individual MIP to improve ITV delineation for tumors adjacent to the thoracic wall, mediastinum, heart, or diaphragm [11, 12]. In the current study, we also had to verify ITVs contoured on MIP and MinIP CT by overlaying it on a movie loop displaying 4D CT data and then editing it, especially for those closed to adjacent organs. Mancosu had recently proposed a semiautomatic technique, which allowed for inclusion of the residual part of ITV covered by liver and spleen cupola when using MIP algorithm. It was validated on phantom and selected patients, which revealed this possibility when lesion located near liver and spleen cupola by performing only the contours on MIP series . Thus, the dedicated software needs to be developed to exclude diaphragm and chest wall in some breathing phases using 4D CT for better tumor MIP/MinIP images.
Theoretically, for tumors with homogeneous hyperdensity or hypodensity compared to the surrounding normal liver, MIP or MinIP projections should accordingly reflect the tumor trajectory across all time-resolved data sets. In patients when the lesions were homogenous low CT density, located intrahepatic, not adjacent to perihepatic organs and also small size (5th, 6th, 8th and 12th patient), using MinIP was also an appropriate tool with good MI of 0.92 ± 0.04 and low proportion of Vunder of 0.08 ± 0.04. Thus, ITVMinIP method can be used for fast contouring of liver tumor with homogeneous low CT density including respiratory motion.
However, a number of liver cancers present inhomogeneous density, results from this study showed combined MIP and MinIP fit excellently in this situation. In current study we compared MIs and proportions of Vunder resulting from ITV2Phase, ITVMinIP, ITVMIP and ITV2M contoured by 4 approaches, and found that the closest to ITVAllPhases was the combined one (ITV2M) of ITVMinIP and ITVMIP, which were contoured on MinIP and MIP series of 4D CT set because it resulted in the highest MI and lowest proportion of Vunder. Moreover, the size of tumor and the ratio of tumor vertical diameter over cranio-caudal movement did not have influence on MI or proportion of Vunder when using ITV2M. In other words, no matter how big the tumor was, and the tumor vertical diameter over cranio-caudal movement was small or big, the tool of ITV2M could always result in the best outcome.
For moving target it was also a practice to sum 2 GTVs to generate ITV, one contoured on CT image set acquired after end-inhale and holding breath and the other contoured on CT acquired after end-exhale and holding breath when 4D CT was not available. The similar method was also used for 4D CT, i.e., to contour only two extreme phases at end-inhalation and end-exhalation. However, possible hysteresis effects would be neglected as occurred in lung cancer [17, 18]. Seppenwoolde and colleagues  reported that when the tumor was small and had a large range of motion, the separation between the positions of the images of inspiration and expiration phases was relatively obvious and the information of the intermediate breathing might not be comprehensive. Besides, the combined images of the two time phases might omit the lag of the tumor. The phenomenon was caused by the time difference among the document recorded by the computer, the transition of respiratory cycle and the transition between inspiration and expiration in a respiratory cycle. Xi  reported the feasibility of using limited 4D CT images for treatment planning for liver radiotherapy. As recognized by the authors, deriving ITV by two extreme phases was reasonably safe only for low and medium tumor motion amplitude (< 1.6 cm). The tumor motion in cranio-caudal direction between Xi's study and ours' were comparable, but the tumor size of our data (152.20 cm3 ± 242.85 cm3) was more diverse than Xi's (70.36 cm3 ± 66.23 cm3). Xi did not investigate the influence of tumor size on ITV determination using 4D CT data. Our study did find a smaller volume of ITV2Phase (186.8 cm3) than that of ITVAllPhases (200.8 cm3) with significant difference (p = 0.004). And also we found that MI and proportion of Vunder were influenced by the tumor volume and the ratio of tumor size over tumor motion magnitude significantly. Whereas, ITV2M was the closest to ITVAllPhases, and MI and proportion of Vunder were not influenced by tumor volume and the ratio of tumor size over tumor motion magnitude when using tool of ITV2M.
To reduce the workload of contouring multiple GTVs in 4D CT data sets, contouring only two extreme phases is appropriate only when tumor volume is big and tumor motion magnitude is relatively small. For hepatic malignancies with inhomogeneous density we found that the method of using ITV2M was a more reliable and appropriate tool for generating ITVs from 4D CT data sets, compared to the others.
Swaminath A, Dawson LA: Emerging role of radiotherapy in the management of liver metastases. Cancer J 2010, 16: 150-155. 10.1097/PPO.0b013e3181d7e8b3
Ren ZG, Zhao JD, Gu K, Chen Z, Lin JH, Xu ZY, Hu WG, Zhou ZH, Liu LM, Jiang GL: Three-dimensional conformal radiation therapy and intensity-modulated radiation therapy combined with transcatheter arterial chemoembolization for locally advanced hepatocellular carcinoma: an irradiation dose escalation study. Int J Radiat Oncol Biol Phys 2011, 79: 496-502. 10.1016/j.ijrobp.2009.10.070
Langen KM, Jones DT: Organ motion and its management. Int J Radiat Oncol Biol Phys 2001, 50: 265-278. 10.1016/S0360-3016(01)01453-5
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
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/015
Liu HH, Balter P, Tutt T, Choi B, Zhang J, Wang C, Chi M, Luo D, Pan T, Hunjan S, Starkschall G, Rosen I, Prado K, Liao Z, Chang J, Komaki R, Cox JD, Mohan R, Dong L: 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 2007, 68: 531-540. 10.1016/j.ijrobp.2006.12.066
Jin JY, Ajlouni M, Chen Q, Yin FF, Movsas B: A technique of using gated-CT images to determine internal target volume (ITV) for fractionated stereotactic lung radiotherapy. Radiother Oncol 2006, 78: 177-184. 10.1016/j.radonc.2005.11.012
Guckenberger M, Sweeney RA, Wilbert J, Krieger T, Richter A, Baier K, Mueller G, Sauer O, Flentje M: Image-guided radiotherapy for liver cancer using respiratory-correlated computed tomography and cone-beam computed tomography. Int J Radiat Oncol Biol Phys 2008, 71: 297-304. 10.1016/j.ijrobp.2008.01.005
Xi M, Liu MZ, Zhang L, Li QQ, Huang XY, Liu H, Hu YH: How many sets of 4DCT images are sufficient to determine internal target volume for liver radiotherapy. Radiother Oncol 2009, 92: 255-259. 10.1016/j.radonc.2009.05.007
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.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. Radiat Oncol 2009, 4: 4. 10.1186/1748-717X-4-4
Rietzel E, Liu AK, Chen GT, Choi NC: Maximum-intensity volumes for fast contouring of lung tumors including respiratory motion in 4DCT planning. Int J Radiat Oncol Biol Phys 2008, 71: 1245-1252. 10.1016/j.ijrobp.2008.03.030
Beddar AS, Briere TM, Balter P, Pan T, Tolani N, Ng C, Szklaruk J, Krishnan S: 4D-CT imaging with synchronized intravenous contrast injection to improve delineation of liver tumors for treatment planning. Radiother Oncol 2008, 87: 445-448. 10.1016/j.radonc.2007.12.009
Sato Y, Shiraga N, Nakajima S, Tamura S, Kikinis R: Local maximum intensity projection (LMIP): a new rendering method for vascular visualization. J Comput Assist Tomogr 1998, 22: 912-917. 10.1097/00004728-199811000-00014
Lin W, Haacke EM, Masaryk TJ, Smith AS: Automated local maximum-intensity projection with three-dimensional vessel tracking. J Magn Reson Imaging 1992, 2: 519-526. 10.1002/jmri.1880020509
Mancosu P, Sghedoni R, Bettinardi V, Aquilina MA, Navarria P, Cattaneo GM, Di MN, Cozzi L, Scorsetti M: Semiautomatic technique for defining the internal gross tumor volume of lung tumors close to liver/spleen cupola by 4D-CT. Med Phys 2010, 37: 4572-4576. 10.1118/1.3471379
Mageras GS, Pevsner A, Yorke ED, Rosenzweig KE, Ford EC, Hertanto A, Larson SM, Lovelock DM, Erdi YE, Nehmeh SA, Humm JL, Ling CC: Measurement of lung tumor motion using respiration-correlated CT. Int J Radiat Oncol Biol Phys 2004, 60: 933-941.
Seppenwoolde Y, Shirato H, Kitamura K, Shimizu S, van HM, Lebesque JV, Miyasaka K: Precise and real-time measurement of 3D tumor motion in lung due to breathing and heartbeat, measured during radiotherapy. Int J Radiat Oncol Biol Phys 2002, 53: 822-834. 10.1016/S0360-3016(02)02803-1
This study was funded by National Science Foundation of China (Grant No. 30800279).
The authors declare that they have no competing interests.
ZJD and JGL designed the study. LJ, WJZ and ZJD did the study and wrote the manuscript. JGL was responsible for manuscript revision and submission. WJZ and XZY were involved in 4D CT simulation and data analysis. All authors read and approved the final version of the manuscript.