Use of kilovoltage X-ray volume imaging in patient dose calculation for head-and-neck and partial brain radiation therapy
Radiation Oncologyvolume 5, Article number: 29 (2010)
To evaluate the accuracy of using kilovoltage x-ray cone-beam computed tomography (kV-CBCT) imaging for in vivo dose calculations.
A Region-of-Interest (ROI) CT number mapping method was developed to generate the cone-beam CT number vs. relative electron density calibration curve for 3D dose calculations. The stability of the results was validated for three consecutive months. The method was evaluated on three brain tumors and three head-and-neck tumor cases. For each patient, kV-CBCT images were acquired on the first treatment day and two-week intervals on the Elekta XVI system. The delivered dose distributions were calculated by applying the patients' treatment plans to the kV-CBCT images. The resulting dose distributions and dose volume histograms (DVHs) of the tumor and critical structures were compared to the original treatment plan.
The kV-CBCT electron density calibration was stable within 1.5% over a three-month period. The DVH and dose distribution comparison based on the planning CT and the initial kV-CBCT showed good agreements for majority of cases. The doses calculated from the planning CT and kV-CBCT were compared on planes perpendicular to the beam axes and passing through the isocenter. Using γ analysis with a criterion of 2 mm/2% and a threshold of 10%, more than 99.5% of the points on the iso-planes exhibited γ <1. For one patient, kV-CBCT images detected 5.8% dose variation in the right parotid due to tumor shrinkage and patient weight loss.
ROI mapping method is an effective method for the creation of kV-CBCT electron density calibration curves for head-and-neck and brain tumor patients. Dose variations as monitored using kV-CBCT imaging suggest that some patients can benefit from adaptive treatment plan re-optimization.
Patients with head-and-neck and definitive brain tumor are routinely treated with intensity-modulated radiotherapy (IMRT) to enable delivery of highly conformal dose distribution to the tumor while sparing surrounding critical structures. Precise target localization is important for such treatments [1–3]. Ideally, the cumulative dose delivered over the whole treatment course should match the total planned dose. However, many uncertainties can be incurred due to patient set-up, anatomic changes and the organ motions during the course of treatment. Barker JL Jr. et al. reported that relative median loss in gross tumor volume was 69.5% and measurable anatomic changes were found throughout the fractionated radiotherapy in head-and-neck patients. As a result of these changes the actual delivered dose deviates from the original planned dose distribution, potentially affecting the tumor control and the normal tissue complication rates.
Cone-beam computed tomography (CBCT) systems mounted on the linear accelerator has become available for image-guided radiotherapy (IGRT). Currently, there are two types of commercially available CBCT imaging systems: (1) the kV-CBCT system, which includes the Varian On-Board-Imaging (OBI) (Varian Medical Systems, Palo Alto, CA) and the Elekta XVI Synergy system (Elekta, Stockholm, Sweden); and (2) the Siemens MVision system (Siemens Medical Solutions, Malvern, PA) [5–7]. In our hospital, we commissioned an Elekta Synergy™ accelerator with on-board kV-CBCT in 2006. The main clinical application of CBCT is to improve the geometric accuracy of target localization in radiation therapy, where the volumetric images of patient acquired immediately before the treatment are registered to the reference planning CT images to correct the patient setup error [8, 9]. KV-CBCT imaging has shown enough soft tissue contrast and spatial resolution for soft-tissue based setup, but the image quality is affected by the acquisition parameters. In principle the kV-CBCT data set can be used to calculate the dose distribution, which means that the planned dose distribution can be evaluated and verified on every treatment day . In order to use CBCT images for dose calculation, the image pixel values need to be converted from dimensionless CT numbers to either electron or physical density. Methods for calibrating conventional fan-beam CT to electron density have been widely used in clinical dose calculation. However, compared to conventional CT scanners, kV-CBCT images have increased artifacts and reduced contrast due to photon scatter. As a result, the calibration of kV-CBCT images for dose calculation is an active area of research [12–15].
The purpose of this study is to assess the feasibility of using a mapping method to calibrate the kV-CBCT images for dose calculation in head-and-neck and definitive brain tumor radiation treatments. By monitoring the dose that patient receives from each fraction, physicians will be able to track the dose distribution during the course of radiation therapy and modify the treatment plan as needed based on the actual dose delivered.
KV CBCT data acquisition
The kV CBCT images were acquired on a linear accelerator equipped with an integrated kV X-ray volumetric imaging system (Elekta, Synergy S, XVI, Crawley, UK). For imaging the head-and-neck and brain tumor patients, we used the following parameters: 100 kVp, S20 collimator and F0 filter, total 65 mAs and a high-resolution reconstruction (512 × 512). A total of about 650 projections were acquired in a full rotation. The CBCT images were reconstructed with slice thickness of 2.5 mm and then transferred to the treatment planning system (TPS, Philips Pinnacle3 V8.0d, Fitchburg, WI, USA) for image registration and dose calculations.
A phantom, Catphan-600 module CTP503 (Phantom Laboratory, NY) was used to evaluate the stability and uniformity of the CBCT numbers. The phantom has seven embedded rods made of different materials: air, PMP, LDPE, polystyrene, acrylic, Delrin and Teflon. Their electron densities relative to water range from 0.00 to 2.16. The CBCT image of the phantom was acquired every month for three consecutive months, and the CBCT numbers were obtained from the TPS and the relative electron densities were recorded accordingly. We also evaluated the maximal fluctuation in CBCT numbers on the image uniformity module part of the phantom.
Calibration of relative electron density
For the dose calculation in a treatment planning system, the relative electron density or physical density of each voxel of the CT images is required for inhomogeneity corrections . In this study, calibration of conventional CT (AcQsim CT Simulator, Philips Medical System, Cleveland, OH) number to physical density was performed on a CT phantom (CIRS model 062, Norfolk, VA). However, each individual patient's CBCT scan has a different scatter component that affects the HU measured in the image. A significant dose error was observed if we directly applied the calibration method for conventional CT to kV-CBCT.
We used a Region-of-interest (ROI) CT number mapping method similar to Richter's report to generate the CT number to physical density conversion curve for the dose calculation . This process was applied to a single patient with head and neck cancer. The CBCT images for this patient were acquired at the same day of planning CT, so that there was no visible change in patient anatomy between the two images. A brief description of the calibration progress is as follows: (1) register the planning CT images and kV-CBCT images in the ADAC Pinnacle treatment planning system; (2) map the regions of interests (ROIs) from conventional CT dataset to the CBCT dataset, and record the mean CBCT number values of these ROIs, and (3) Generate the kV-CBCT numbers to physical density calibration curve based on the density values measured on the conventional CT.
Three head-and-neck cases and three brain patients with different tumor sites treated on Elekta Synergy were selected for retrospective evaluation of the accuracy of CBCT-based dose calculations. The head-and-neck cases included two natural killer/T-cell (NK/T) lymphoma cases and one nasopharyngeal carcinoma (NPC) case. For all patients, conventional CT was acquired with slice thickness of 5 mm and the target and critical structures were delineated by the attending physicians. IMRT plans were designed according to the physician's prescriptions with beams of 6 MV. The beam angles were 0, 50, 110, 250 and 310 degree for NK/T cases, and 0, 45, 90, 120, 160, 200, 240, 280 and 320 degree for NPC cases. For brain cases, gantry angles were 0, 60, 230, 300 and 45 degree with a 90-degree couch kick. Two sets of CBCT images were acquired, one on the first treatment day and the other on two weeks later. The patients were initially set up to the skin markers then followed with a CBCT scanning. The CBCT images were acquired according to the appropriate protocol and then reconstructed with slice thickness of 2.5 mm. All images were transferred to treatment planning system for analysis. For each case, the CBCT images were first registered to the conventional kV-CT images using an automatic registration method based on normalized mutual information algorithm, and then manual adjustments were performed to achieve the optimal match. A second set of CBCT images were acquired two weeks later and registered to the reference planning CT with the same method. The contours were mapped from CT to CBCT images with slight changes if them were beyond the skin.
For dose calculation, all the parameters (e.g., iso-center location, beam angles, MLC shapes, and monitor units) in the initial treatment plan were applied to the kV-CBCT images, and then the dose distribution was recalculated based on the new calibration curve. The dose calculation was performed in the Pinnacle treatment planning system using the collapsed cone superposition convolution algorithm with an isotropic 2 mm dose grid resolution. The contours delineated on the conventional CT were also mapped onto the kV-CBCT image data sets. Finally, the initial dose distribution matrix calculated on the planning CT was imported in the treatment planning system and displayed on the kV-CBCT dataset using scripts developed in-house. Dose volume histogram and the dose to tumor and normal structures were compared on the two image data sets. The differences in the dose distributions of the two plans were analyzed using γ analysis along planes through the isocenter perpendicular to each beam axis using commercial software (MapCheck, Version4.0, Sun Nuclear, Melbourne, FL).
The stability of kV-CBCT numbers
Because most patients complete their treatment courses within five weeks, we consider the three-month length of the stability test to be adequate. The maximal difference in CBCT numbers was 21, with a maximum standard error of less than 1.5%. The stability of kV-CBCT number and electron density indicates acceptable overall performance of the kV-CBCT system. The kV-CBCT images of the uniformity section of the phantom shows the maximal fluctuation of the CBCT numbers is ± 35 Hounsfield unit (HU), which translates to a fluctuation of approximately 1% in electron density values.
Conversion of kV-CBCT numbers to relative electron density
A total of 13 different ROIs were used in generating the conversion curve, which include air, skin, muscle, brain stem, spinal cord, parotid gland, outer bone, inner bone, tooth and other structures. Table 1 shows the CBCT numbers and their corresponding physical density values. The calibration curves, as shown in Figure 1, were implemented in treatment planning system for the dose calculations. Large discrepancies were noted from these two curves. In particular, some discontinuous steps were observed on the calibration curve of kV-CBCT images.
Only minimal changes and deformations were observed in the anatomical structures on the patients' first CBCT images as compared with the reference planning CT. The DVHs of one NK/T lymphoma case (patient1), one NPC case (patient2) and one brain tumor case (patient3) are shown in Figure 2 as an example. The solid lines represented the DVHs based on conventional CT images and the dash lines were based on the dose calculated from the KV-CBCT images. Figure 3(a)-(c) are the dose distributions on the transverse planes of the three patients. The left images represent dose distributions based on the kV-CBCT and the right had images represent the dose on kV-CT images. There is no significant dose difference between the conventional CT images and kV-CBCT images.
For all clinical cases, the dose comparison was performed at a plane through the isocenter for each individual beam. Good agreement was found between the conventional CT and the first kV-CBCT based dose calculations. Using the γ analysis with a criterion of 2 mm and 2% and a threshold of 10%, more than 99.5% of the points at the iso-plane have the γ value less than 1.0. Table 2 shows the distance to agreement (DTA) and gamma analysis results of the three cases. For most of the beams, the pass rate for distance to agreement were better than 96% except one beam which has the data of 94.5%.
Table 3 shows the dose to the tumor and some normal structures of the three patients in the planning CT data sets and the first CBCT data sets. The differences of the dose to tumor and some normal tissues were within 1% and 3.2%, respectively. The difference of maximal dose in tumor is 0.49% and in normal structures are 3.15%.
Five out of six patients didn't show significant anatomy changes and setup variations between the first CBCT images and the second CBCT images. But for one NK/T patient (patient5), a slight anatomical change in the patient's skin contour and air cavity was found in the second CBCT images compared to the conventional CT images, as shown in Figure 4. The dose comparisons of the reference kVCT, the 1st and 2nd kV CBCTs for that patient are listed in Table 4. On the first treatment day, the dose difference in gross tumor volume (GTV), clinical tumor volume (CTV) and planning tumor volume (PTV) between reference CT and cone beam CT (kV-CBCT1) were 0.98%, 0.54%, 0.54%, respectively. The maximal dose difference was found on the spinal cord (-1.87%). For the second cone beam CT (kV-CBCT2) acquired two weeks later, the maximal dose difference of spinal cord increased to 3.77%, and the maximal dose difference was found in the right parotid (5.81%). While for tumor and other structures, the dose agreement was still within 1.0%.
On-board CBCT volumetric imaging can improve the accuracy of radiation therapy in two aspects, namely target localization and delivered dose verification [12, 17, 18]. By acquiring 3D CBCT images with patient on the treatment couch just before the treatment delivery, patient setup error can be corrected and the accuracy of target positioning localization accuracy can be improved. CBCT image data sets obtained throughout the treatment course can be used for dose calculation, hence providing a clinical quality assurance tool for radiotherapy. However, the CBCT image quality is susceptible to many factors, such as scattering, beam hardening effects and organ motion, etc . Morin et al studied dose calculation based on MV CBCT images and reported that the MV CBCT could be used to estimate the dose variation due to the anatomical changes in the head-and-neck region . In this study, we investigated the feasibility and the accuracy of using kV CBCT images for direct dose calculation in head-and-neck and brain tumor radiotherapy with a simple and effective method.
The relative electron density can directly affect the dose calculation accuracy when inhomogeneity correction is involved. Unlike conventional kV CT or MV CBCT, kV-CBCT has a larger scatter radiation component and the image quality suffers from the beam hardening effect [21–23]. It has been reported that the effect of scatter radiation can be partly corrected or reduced by calibrating the kV-CBCT system [24, 25]. The conventional CT number to the relative electron density conversion was performed with a CT number calibration phantom embedded with different types of tissue-equivalent inserts. However if such method was used directly for kV-CBCT, a dose calculation error can be introduced. Based on our tests, if we used the calibration curve generated by the phantom directly, the dose difference between first CBCT and planning CT would be more than 5%, which agrees with the reports from Yang et al and Tucking et al [12, 26].
In this study, we used the ROI mapping method to generate calibration curve for kV CBCT image-based dose calculation . Obviously, accurate image registration is needed for this method. The registration of different image modalities is widely used in radiotherapy for delineating the region of interests. As the registration algorithms in the commercial treatment planning systems generally use rigid body transformations, we selected the head-and-neck and brain tumor cases for our study, where this assumption was generally valid.
The calibration curve for kV-CBCT is different from that for conventional CT. The conventional CT number is zero for the air outside the patient skin in the planning system; however, the CBCT number in such situation is much greater than zero. The mean CBCT value in the air around the skin is 1379 for the selected case, similar to other report . The steps in the kV CBCT electron density conversion curve is mainly caused by scatter and beam hardening effects.
A good agreement of the calculated doses to the tumor and normal structures was found between the conventional CT and the first kV-CBCT images because there were virtually no anatomical changes between these images. The maximal dose deviation was found in the eye mainly due to the residual registration error and contour deviations, as the slice thickness was 5 mm for the conventional CT and 2.5 mm for CBCT images. The structures near the skin showed larger differences. The DTA and γ index analysis results also showed the good agreement between kV-CBCT based and conventional CT based dose calculation. Richter et al. used the same method and reported the dose difference between the planning CT and CBCT was 1.36% ± 1.96% in head patients with three-dimensional conformal plans. Our data showed the difference was within 1% of the target, which was consistent with their result. Our results demonstrated that the mapping method for CBCT correction is accurate both for three-dimensional conformal plans and IMRT plans in head and brain cases.
Furthermore, we generated the density conversion table based on one patient and applied the same table to the other patients. There were only a small discrepancy between the doses calculated by using kV CBCT and conventional CT in all 6 cases with different tumor locations. This result suggests that, for head-and-neck and brain patients, variations in the scatter effect in imaging different tumor sites is relatively small from patient to patient, and therefore it is reasonable to use the same electron density conversion curve for kV CBCT based dose calculation. Compared to the patient group based conversion table in report of Richter et al. or CT-based HU mapping method in Mathilda et al., this specific case mapping method is less complex to develop and implement, but it is limited to the preset scanning parameters.
Overall, our study showed good accuracy in CBCT based dose calculation. However, it is not recommended to replace the conventional planning CT by kV CBCT for the purpose of treatment planning as the inferior image quality of kV CBCT may affect the accuracy of target and normal structures delineation.
The kV CBCT can also be used to evaluate the dose to tumor or the normal structures. In this study, one NK/T patient had slight changes in anatomy after two weeks' treatment, dose variations were found in the spinal cord and the right parotid gland. These results suggested that even within a relatively short period such as 2 weeks, dose verification based on CBCT or CT will be necessary for certain patients to account for dosimetric effects due to patient anatomical changes.
Anatomic changes for head-and-neck patients, including nodal mass shrinkage and patient weight loss during the course of radiation therapy, can occur [28, 29]. For these cases, repeat CT imaging and re-planning may be essential to ensure the adequate dose delivered to the tumor and proper sparing of the surrounding sensitive structures.
Technically, the 26 cm field of view for the S20 collimator may limit the use of kV CBCT for dose calculation of patients with beams going through their shoulders. However for most head-and-neck patients, the FOV is sufficiently large to evaluate the dose to PTV, brain stem, spinal cord, eyes and parotid glands. For those patients who receive thoracic or pelvic treatment, S20 is not large enough to encompass all the structures and skin. Dose verification for other sites is part of our future research.
ROI mapping method is a feasible method to overcome the effects of scatter for generating the kV CBCT relative electron density calibration curve for head-and-neck cancer and brain tumor patients. Dose variations as monitored using kV CBCT imaging were observed in a relatively short period of two weeks, which suggests potential benefits of adaptive treatment plan re-optimization for certain head-and-neck and brain tumor patients.
Hong TS, Tomé WA, Chappell RJ, Chinnaiyan P, Mehta MP, Harari PM: The impact of daily setup variations on head-and-neck intensity-modulated radiation therapy. Int J Radiat Oncol Biol Phys 2005, 61: 779-88. 10.1016/j.ijrobp.2004.07.696
Xing L, Lin Z, Donaldson SS, Le QT, Tate D, Goffinet DR, Wolden S, Ma L, Boyer AL: Dosimetric effects of patient displacement and collimator and gantry angle misalignment on intensity modulataed radiation therapy. Radiother Oncol 2000, 56: 97-108. 10.1016/S0167-8140(00)00192-4
Han C, Chen YJ, Liu A, Schultheiss TE, Wong JY: Actual dose variation of parotid glands and spinal cord for nasopharyngeal cancer patients during radiotherapy. Int J Radiat Oncol Biol Phys 2008, 70: 1256-62. 10.1016/j.ijrobp.2007.10.067
Barker JL Jr, Garden AS, Ang KK, O'Daniel JC, Wang H, Court LE, Morrison WH, Rosenthal DI, Chao KS, Tucker SL, Mohan R, Dong L: Quantification of volumetric and geometric changes occurring during fractionated radiotherapy for head-and-neck cancer using an integrated CT/linear accelerator system. Int J Radiat Oncol Biol Phys 2004, 59: 960-70. 10.1016/j.ijrobp.2003.12.024
Kim GY, Pawlicki T, Le QT, Luxton G: Linac-based on-board imaging feasibility and the dosimetric consequences of head roll in head-and-neck IMRT plans. Med Dosim 2008,33(1):93-9. 10.1016/j.meddos.2007.05.004
Oelfke U, Tücking T, Nill S, Seeber A, Hesse B, Huber P, Thilmann C: Linac-integrated kV-cone beam CT: technical features and first applications. Med Dosim 2006, 31: 62-70. 10.1016/j.meddos.2005.12.008
Morin O, Gillis A, Chen J, Aubin M, Bucci MK, Roach M, Pouliot J: Megavoltage cone-beam CT: system description and clinical applications. Med Dosim 2006, 31: 51-61. 10.1016/j.meddos.2005.12.009
Chang J, Mageras GS, Yorke E, De Arruda F, Sillanpaa J, Rosenzweig KE, Hertanto A, Pham H, Seppi E, Pevsner A, Ling CC, Amols H: Observation of interfractional variations in lung tumor position using respiratory gated and ungated megavoltage cone-beam computed tomography. Int J Radiat Oncol Biol Phys 2007, 67: 1548-58. 10.1016/j.ijrobp.2006.11.055
Fatunase T, Wang Z, Yoo S, Hubbs JL, Prosnitz RG, Yin FF, Marks LB: Assessment of the Residual Error in Soft Tissue Setup in Patients Undergoing Partial Breast Irradiation: Results of a Prospective Study Using Cone-Beam Computed Tomography. Int J Radiat Oncol Biol Phys 2008, 70: 1025-34. 10.1016/j.ijrobp.2007.07.2344
Varadhan R, Hui SK, Way S, Nisi K: Assessing prostate, bladder and rectal doses during image guided radiation therapy--need for plan adaptation? J Appl Clin Med Phys 2009,10(3):2883. 10.1120/jacmp.v10i3.2883
Thomas SJ: Relative electron density calibration of CT scanners for radiotherapy treatment planning. Br J Radiol 1999, 72: 781-6.
Yang Y, Schreibmann E, Li T, Wang C, Xing L: Evaluation of on-board kV cone -beam CT (CBCT)-based dose calculation. Phys Med Biol 2007,52(3):685-705. 10.1088/0031-9155/52/3/011
van Zijtveld M, Dirkx M, Heijmen B: Correction of conebeam CT values using a planning CT for derivation of the "dose of the day". Radiother Oncol 2007,85(2):195-200. 10.1016/j.radonc.2007.08.010
Richter A, Hu Q, Steglich D, Baier K, Wilbert J, Guckenberger M, Flentje M: Investigation of the usability of conebeam CT data sets for dose calculation. Radiat Oncol 2008, 3: 42. 10.1186/1748-717X-3-42
Guan H, Dong H: Dose calculation accuracy using cone-beam CT (CBCT) for pelvic adaptive radiotherapy. Phys Med Biol 2009,54(20):6239-50. 10.1088/0031-9155/54/20/013
Low DA, Harms WB, Mutic S, Purdy JA: A technique for the quantitative evaluation of dose distributions. Med Phys 1998, 25: 656-61. 10.1118/1.598248
Thilmann C, Nill S, Tücking T, Höss A, Hesse B, Dietrich L, Bendl R, Rhein B, Häring P, Thieke C, Oelfke U, Debus J, Huber P: Correction of patient positioning errors based on in-line cone beam CTs: clinical implementation and first experiences. Radiat Oncol 2006, 1: 16. 10.1186/1748-717X-1-16
Grills IS, Hugo G, Kestin LL, Galerani AP, Chao KK, Wloch J, Yan D: Image-Guided Radiotherapy Via Daily Online Cone-Beam CT Substantially Reduces Margin Requirements for Stereotactic Lung Radiotherapy. Int J Radiat Oncol Biol Phys 2008, 70: 1045-56. 10.1016/j.ijrobp.2007.07.2352
Jaffray DA, Battista JJ, Fenster A, Munro P: X-ray scatter in megavoltage transmission radiography: physical characteristics and influence on image quality. Med Phys 1994, 21: 45-60. 10.1118/1.597255
Morin O, Chen J, Aubin M, Gillis A, Aubry JF, Bose S, Chen H, Descovich M, Xia P, Pouliot J: Dose calculation using megavoltage cone-beam CT. Int J Radiat Oncol Biol Phys 2007, 67: 1201-10. 10.1016/j.ijrobp.2006.10.048
Ding GX, Duggan DM, Coffey CW: Characteristics of kilovoltage x-ray beams used for cone-beam computed tomography in radiation therapy. Phys Med Biol 2007, 52: 1595-615. 10.1088/0031-9155/52/6/004
Siewerdsen JH, Jaffray DA: Cone-beam computed tomography with a flat-panel imager: magnitude and effects of x-ray scatter. Med Phys 2001, 28: 220-31. 10.1118/1.1339879
Reitz I, Hesse BM, Nill S, Tücking T, Oelfke U: Enhancement of image quality with a fast iterative scatter and beam hardening correction method for kV CBCT. Z Med Phys 2009,19(3):158-72.
Létourneau D, Wong R, Moseley D, Sharpe MB, Ansell S, Gospodarowicz M, Jaffray DA: Online planning and delivery technique for radiotherapy of spinal metastases using cone-beam CT: image quality and system performance. Int J Radiat Oncol Biol Phys 2007, 67: 1229-37.
Siewerdsen JH, Jaffray DA: Optimization of x-ray imaging geometry (with specific application to flat-panel cone-beam computed tomography). Med Phys 2000, 27: 1903-14. 10.1118/1.1286590
Tucking T, Nill S, Oelfke U: Dose calculation on linac integrated KV-cone beam CT. Radiotherapy and oncology 2006,81(Supplement 1):26.
Brock KK: Image registration in intensity-modulated, image-guided and stereotactic body radiation therapy. Front Radiat Ther Oncol 2007, 40: 94-115. full_text
Beaver ME, Matheny KE, Roberts DB, Myers JN: Predictors of weight loss during radiation therapy. Otolaryngol Head Neck Sur 2001, 125: 645-8. 10.1067/mhn.2001.120428
Hansen EK, Bucci MK, Quivey JM, Weinberg V, Xia P: Repeat CT imaging and replanning during the course of IMRT for head-and-neck cancer. Int J Radiat Oncol Biol Phys 2006, 64: 355-62. 10.1016/j.ijrobp.2005.07.957
The authors thank the Drs. Lijun Ma and Andrew Huwang for helpful discussions and editing of the paper.
Declaration of competing interests
The authors declare that they have no competing interests.
Each author has participated sufficiently in the work to take public responsibility for appropriate portions of the content. JY, ZZ designed the study. WH, JW performed the study and analysis. XM provided the patients' images. The manuscript was written by WH, all other authors helped and finally approved the final manuscript.