Dosimetric and radiobiological impact of dose fractionation on respiratory motion induced IMRT delivery errors: a volumetric dose measurement study

Respiratory motion can introduce substantial dose errors during IMRT delivery. These errors are difficult to predict because of the nonsynchronous interplay between radiation beams and tissues. The present study investigates the impact of dose fractionation on respiratory motion induced dosimetric errors during IMRT delivery and their radiobiological implications by using measured 3D dose. We focused on IMRT delivery with dynamic multileaf collimation (DMLC-IMRT). IMRT plans using several beam arrangements were optimized for and delivered to a polystyrene phantom containing a simulated target and critical organs. The phantom was set in linear sinusoidal motion at a frequency of 15 cycles/min (0.25 Hz). The amplitude of the motion was +/- 0.75 cm in the longitudinal direction and +/- 0.25 cm in the lateral direction. Absolute doses were measured with a 0.125 cc ionization chamber while dose distributions were measured with transverse films spaced 6 mm apart. Measurements were performed for varying number of fractions with motion, with respiratory-gated motion, and without motion. A tumor control probability (TCP) model for an inhomogeneously irradiated tumor was used to calculate and compare TCPs for the measurements and the treatment plans. Equivalent uniform doses (EUD) were also computed. For individual fields, point measurements using an ionization chamber showed substantial dose deviations (-11.7% to 47.8%) for the moving phantom as compared to the stationary phantom. However, much smaller deviations (-1.7% to 3.5%) were observed for the composite dose of all fields. The dose distributions and DVHs of stationary and gated deliveries were in good agreement with those of treatment plans, while those of the nongated moving phantom showed substantial differences. Compared to the stationary phantom, the largest differences observed for the minimum and maximum target doses were -18.8% and +19.7%, respectively. Due to their random nature, these dose errors tended to average out over fractionated treatments. The results of five-fraction measurements showed significantly improved agreement between the moving and stationary phantom. The changes in TCP were less than 4.3% for a single fraction, and less than 2.3% for two or more fractions. Variation of average EUD per fraction was small (< 3.1 cGy for a fraction size of 200 cGy), even when the DVHs were noticeably different from that of the stationary tumor. In conclusion, IMRT treatment of sites affected by respiratory motion can introduce significant dose errors in individual field doses; however, these errors tend to cancel out between fields and average out over dose fractionation. 3D dose distributions, DVHs, TCPs, and EUDs for stationary and moving cases showed good agreement after two or more fractions, suggesting that tumors affected by respiration motion may be treated using IMRT without significant dosimetric and biological consequences.     

Improving the resolution of dynamic intensity modulated radiation therapy delivery by reducing the multileaf collimator sampling distance

The conformality of a dose distribution delivered by a multileaf collimator (MLC) for intensity modulated radiation therapy (IMRT) is limited in the direction perpendicular to leaf motion by the finite leaf width. Two methods of improving the resolution of IMRT intensity maps in this direction were investigated. In the first, the desired fluence distribution is considered to be sampled by the MLC, with the sampling distance being the center-to-center distance between the MLC leaves. The sampling distance is reduced below the leaf width by combining separate irradiations with a couch shift between them. This has been applied to static field therapy [Galvin et al., Int. J. Radiat. Oncol., Biol., Phys. 35, 89-94 (1996)], and was proposed for IMRT by Bortfeld et al. [Med. Phys. 27, 2494-2502 (2000)]. In the second method, two MLC component fluences, with leaf width L = 2deltay and offset by deltay, are combined to reproduce desired intensity bins with deltay width. The effect of MLC leaf sampling distance on dose resolution was quantified for both 1.0 and 0.5 cm MLC leaf widths, utilizing a high resolution bar-pattern fluence, an annular shaped fluence, and an intensity step-edge. Improvement in resolution was found for the 1.0 cm leaf width at a sampling distance of 0.5 cm, with only a small benefit for further reduction. For the 0.5 cm leaf width, a sampling distance of 0.25 cm resulted in a dose resolution that was nearly independent of direction.     

Real-time prediction of respiratory motion based on a local dynamic model in an augmented space

Motion-adaptive radiotherapy aims to deliver ablative radiation dose to the tumor target with minimal normal tissue exposure, by accounting for real-time target movement. In practice, prediction is usually necessary to compensate for system latency induced by measurement, communication and control. This work focuses on predicting respiratory motion, which is most dominant for thoracic and abdominal tumors. We develop and investigate the use of a local dynamic model in an augmented space, motivated by the observation that respiratory movement exhibits a locally circular pattern in a plane augmented with a delayed axis. By including the angular velocity as part of the system state, the proposed dynamic model effectively captures the natural evolution of respiratory motion. The first-order extended Kalman filter is used to propagate and update the state estimate. The target location is predicted by evaluating the local dynamic model equations at the required prediction length. This method is complementary to existing work in that (1) the local circular motion model characterizes 'turning', overcoming the limitation of linear motion models; (2) it uses a natural state representation including the local angular velocity and updates the state estimate systematically, offering explicit physical interpretations; (3) it relies on a parametric model and is much less data-satiate than the typical adaptive semiparametric or nonparametric method. We tested the performance of the proposed method with ten RPM traces, using the normalized root mean squared difference between the predicted value and the retrospective observation as the error metric. Its performance was compared with predictors based on the linear model, the interacting multiple linear models and the kernel density estimator for various combinations of prediction lengths and observation rates. The local dynamic model based approach provides the best performance for short to medium prediction lengths under relatively low observation rate. Sensitivity analysis indicates its robustness toward the choice of parameters. Its simplicity, robustness and low computation cost makes the proposed local dynamic model an attractive tool for real-time prediction with system latencies below 0.4 s.     

Dosimetric consequences of tumour motion due to respiration for a scanned proton beam

A method for simulating spot-scanned delivery to a moving tumour was developed which uses patient-specific image and plan data. The magnitude of interplay effects was investigated for two patient cases under different fractionation and respiratory motion variation scenarios. The use of volumetric rescanning for motion mitigation was also investigated. For different beam arrangements, interplay effects lead to severely distorted dose distributions for a single fraction delivery. Baseline shift variations for single fraction delivery reduced the dose to the clinical target volume (CTV) by up to 14.1 Gy. Fractionated delivery significantly reduced interplay effects; however, local overdosage of 12.3% compared to the statically delivered dose remained for breathing period variations. Variations of the tumour baseline position and respiratory period were found to have the largest influence on target inhomogeneity; these effects were reduced with fractionation. Volumetric rescanning improved the dose homogeneity. For the CTV, underdosage was improved by up to 34% in the CTV and overdosage to the lung was reduced by 6%. Our results confirm that rescanning potentially increases the dose homogeneity; however, it might not sufficiently compensate motion-induced dose distortions. Other motion mitigation techniques may be required to additionally treat lung tumours with scanned proton beams.     

The effects of intra-fraction organ motion on the delivery of intensity-modulated field with a multileaf collimator

Intensity-modulated radiation therapy can be conveniently delivered with a multileaf collimator. With this method, the entire field is not delivered at once, but rather it is composed of many subfields defined by the leaf positions as a function of beam on time. At any given instant, only these subfields are delivered. During treatment, if the organ moves, part of the volume may move in or out of these subfields. Due to this interplay between organ motion and leaf motion the delivered dose may be different from what was planned. In this work, we present a method that calculates the effects of organ motion on delivered dose. The direction of organ motion may be parallel or perpendicular to the leaf motion, and the effect can be calculated for a single fraction or for multiple fractions. Three breast patients and four lung patients were included in this study,with the amplitude of the organ motion varying from +/- 3.5 mm to +/- 10 mm, and the period varying from 4 to 8 seconds. Calculations were made for these patients with and without organ motion, and results were examined in terms of isodose distribution and dose volume histograms. Each calculation was repeated ten times in order to estimate the statistical uncertainties. For selected patients, calculations were also made with conventional treatment technique. The effects of organ motion on conventional techniques were compared relative to that on IMRT techniques. For breast treatment, the effect of organ motion primarily broadened the penumbra at the posterior field edge. The dose in the rest of the treatment volume was not significantly affected. For lung treatment, the effect also broadened the penumbra and degraded the coverage of the planning target volume (PTV). However, the coverage of the clinical target volume (CTV) was not much affected, provided the PTV margin was adequate. The same effects were observed for both IMRT and conventional treatment techniques. For the IMRT technique, the standard deviations of ten samples of a 30-fraction calculation were very small for all patients, implying that over a typical treatment course of 30 fractions, the delivered dose was very close to the expected value. Hence, under typical clinical conditions, the effect of organ motion on delivered dose can be calculated without considering the interplay between the organ motion and the leaf motion. It can be calculated as the weighted average of the dose distribution without organ motion with the distribution of organ motion. Since the effects of organ motion on dose were comparable for both IMRT and conventional techniques, the PTV margin should remain the same for both techniques.     

Dosimetric verification of intensity modulated beams produced with dynamic multileaf collimation using an electronic portal imaging device

Dose distributions can often be significantly improved by modulating the two-dimensional intensity profile of the individual x-ray beams. One technique for delivering intensity modulated beams is dynamic multileaf collimation (DMLC). However, DMLC is complex and requires extensive quality assurance. In this paper a new method is presented for a pretreatment dosimetric verification of these intensity modulated beams utilizing a charge-coupled device camera based fluoroscopic electronic portal imaging device (EPID). In the absence of the patient, EPID images are acquired for all beams produced with DMLC. These images are then converted into two-dimensional dose distributions and compared with the calculated dose distributions. The calculations are performed with a pencil beam algorithm as implemented in a commercially available treatment planning system using the same absolute beam fluence profiles as used for calculation of the patient dose distribution. The method allows an overall verification of (i) the leaf trajectory calculation (including the models to incorporate collimator scatter and leaf transmission), (ii) the correct transfer of the leaf sequencing file to the treatment machine, and (iii) the mechanical and dosimetrical performance of the treatment unit. The method was tested for intensity modulated 10 and 25 MV photon beams; both model cases and real clinical cases were studied. Dose profiles measured with the EPID were also compared with ionization chamber measurements. In all cases both predictions and EPID measurements and EPID and ionization chamber measurements agreed within 2% (1 sigma). The study has demonstrated that the proposed method allows fast and accurate pretreatment verification of DMLC.     

The impact of respiratory motion and treatment technique on stereotactic body radiation therapy for liver cancer

Stereotactic body radiation therapy (SBRT), which delivers a much higher fractional dose than conventional treatment in only a few fractions, is an effective treatment for liver metastases. For patients who are treated under free-breathing conditions, however, respiration-induced tumor motion in the liver is a concern. Limited clinical information is available related to the impact of tumor motion and treatment technique on the dosimetric consequences. This study evaluated the dosimetric deviations between planned and delivered SBRT dose in the presence of tumor motion for three delivery techniques: three-dimensional conformal static beams (3DCRT), dynamic conformal arc (DARC), and intensity-modulated radiation therapy (IMRT). Five cases treated with SBRT for liver metastases were included in the study, with tumor motions ranging from 0.5 to 1.75 cm. For each case, three different treatment plans were developed using 3DCRT, DARC, and IMRT. The gantry/multileaf collimator (MLC) motion in the DARC plans and the MLC motion in the IMRT plans were synchronized to the patient's respiratory motion. Retrospectively sorted four-dimensional computed tomography image sets were used to determine patient-organ motion and to calculate the dose delivered during each respiratory phase. Deformable registration, using thin-plate-spline models, was performed to encode the tumor motion and deformation and to register the dose-per-phase to the reference phase images. The different dose distributions resulting from the different delivery techniques and motion ranges were compared to assess the effect of organ motion on dose delivery. Voxel dose variations occurred mostly in the high gradient regions, typically between the target volume and normal tissues, with a maximum variation up to 20%. The greatest CTV variation of all the plans was seen in the IMRT technique with the largest motion range (D99: -8.9%, D95: -8.3%, and D90: -6.3%). The greatest variation for all 3DCRT plans was less than 2% for D95. Dose variations for DARC fell between the 3DCRT and IMRT techniques. The dose volume histogram variations for normal organs were negligible. Therefore, the IMRT technique may be a preferable treatment choice in cases where the target volume and critical organs are in close proximity, or when normal organ protection is a high priority, provided that motion effect for the target volume can be managed.     

Dosimetric verification of a commercial 3D treatment planning system for conformal radiotherapy with a dynamic multileaf collimator

The dosimetric accuracy of a 3D treatment planning system (TPS) for conformal radiotherapy with a computer-assisted dynamic multileaf collimator (DMLC) was evaluated. The DMLC and the TPS have been developed for clinical applications where dynamic fields not greater than 10 x 10 cm2 and multiple non-coplanar arcs are required. Dosimetric verifications were performed by simulating conformal treatments of irregularly shaped targets using several arcs of irradiation with 6 MV x-rays and a spherical-shaped, tissue-simulating phantom. The accuracy of the delivered dose at the isocentre was verified using an ionization chamber placed in the centre of the phantom. Isodose distributions in the axial and sagittal planes passing through the centre of the phantom were measured using double-layer radiochromic films. Measured dose at the isocentre as well as isodose distributions were compared to those calculated by the TPS. The maximum percentage difference between measured and prescribed dose was less than 2.5% for all the simulated treatment plans. The mean (+/-SD) displacement between measured and calculated isodoses was, in the axial planes, 1.0 +/- 0.6 mm, 1.2 +/- 0.7 mm and 1.5 +/- 1.1 mm for 80%, 50% and 20% isodose curves, respectively, whereas in the sagittal planes it was 2.0 +/- 1.2 mm and 2.2 +/- 2 mm for 80% and 50% isodose curves, respectively. The results indicate that the accuracy of the 3D treatment planning system used with the DMLC is reasonably acceptable in clinical applications which require treatments with several non-coplanar arcs and small dynamic fields.     

Systematic errors in respiratory gating due to intrafraction deformations of the liver

This article shows the limitations of respiratory gating due to intrafraction deformations of the right liver lobe. The variability of organ shape and motion over tens of minutes was taken into account for this evaluation, which closes the gap between short-term analysis of a few regular cycles, as it is possible with 4DCT, and long-term analysis of interfraction motion. Time resolved MR volumes (4D MR sequences) were reconstructed for 12 volunteers and subsequent non-rigid registration provided estimates of the 3D trajectories of points within the liver over time. The full motion during free breathing and its distribution over the liver were quantified and respiratory gating was simulated to determine the gating accuracy for different gating signals, duty cycles, and different intervals between patient setup and treatment. Gating effectively compensated for the respiratory motion within short sequences (3 min), but deformations, mainly in the anterior inferior part (Couinaud segments IVb and V), led to systematic deviations from the setup position of more than 5 mm in 7 of 12 subjects after 20 min. We conclude that measurements over a few breathing cycles should not be used as a proof of accurate reproducibility of motion, not even within the same fraction, if it is longer than a few minutes. Although the diaphragm shows the largest magnitude of motion, it should not be used to assess the gating accuracy over the entire liver because the reproducibility is typically much more limited in inferior parts. Simple gating signals, such as the trajectory of skin motion, can detect the exhalation phase, but do not allow for an absolute localization of the complete liver over longer periods because the drift of these signals does not necessarily correlate with the internal drift.     

Real-time prediction of respiratory motion based on local regression methods

Recent developments in modulation techniques enable conformal delivery of radiation doses to small, localized target volumes. One of the challenges in using these techniques is real-time tracking and predicting target motion, which is necessary to accommodate system latencies. For image-guided-radiotherapy systems, it is also desirable to minimize sampling rates to reduce imaging dose. This study focuses on predicting respiratory motion, which can significantly affect lung tumours. Predicting respiratory motion in real-time is challenging, due to the complexity of breathing patterns and the many sources of variability. We propose a prediction method based on local regression. There are three major ingredients of this approach: (1) forming an augmented state space to capture system dynamics, (2) local regression in the augmented space to train the predictor from previous observation data using semi-periodicity of respiratory motion, (3) local weighting adjustment to incorporate fading temporal correlations. To evaluate prediction accuracy, we computed the root mean square error between predicted tumor motion and its observed location for ten patients. For comparison, we also investigated commonly used predictive methods, namely linear prediction, neural networks and Kalman filtering to the same data. The proposed method reduced the prediction error for all imaging rates and latency lengths, particularly for long prediction lengths.     

Dosimetric effect of respiration-gated beam on IMRT delivery

Intensity modulated radiation therapy (IMRT) with a dynamic multileaf collimator (DMLC) requires synchronization of DMLC leaf motion with dose delivery. A delay in DMLC communication is known to cause leaf lag and lead to dosimetric errors. The errors may be exacerbated by gated operation. The purpose of this study was to investigate the effect of leaf lag on the accuracy of doses delivered in gated IMRT. We first determined the effective leaf delay time by measuring the dose in a stationary phantom delivered by wedge-shaped fields. The wedge fields were generated by a DMLC at various dose rates. The so determined delay varied from 88.3 to 90.5 ms. The dosimetric effect of this delay on gated IMRT was studied by delivering wedge-shaped and clinical IMRT fields to moving and stationary phantoms at dose rates ranging from 100 to 600 MU/min, with and without gating. Respiratory motion was simulated by a linear sinusoidal motion of the phantom. An ionization chamber and films were employed for absolute dose and 2-D dose distribution measurements. Discrepancies between gated and nongated delivery to the stationary phantom were observed in both absolute dose and 2-D dose distribution measurements. These discrepancies increased monotonically with dose rate and frequency of beam interruptions, and could reach 3.7% of the total dose delivered to a 0.6 cm3 ion chamber. Isodose lines could be shifted by as much as 3 mm. The results are consistent with the explanation that beam hold-offs in gated delivery allowed the lagging leaves to catch up with the delivered monitor units each time that the beam was interrupted. Low dose rates, slow leaf speeds and low frequencies of beam interruptions reduce the effect of this delay-and-catch-up cycle. For gated IMRT it is therefore important to find a good balance between the conflicting requirements of rapid dose delivery and delivery accuracy.     

Effects of averaging over motion and the resulting systematic errors in radiation therapy

The potential for systematic errors in radiotherapy of a breathing patient is considered using the statistical model of Bortfeld et al (2002 Phys. Med. Biol. 47 2203-20). It is shown that although averaging over 30 fractions does result in a narrow Gaussian distribution of errors, as predicted by the central limit theorem, the fact that one or a few samples of the breathing patient's motion distribution are used for treatment planning (in contrast to the many treatment fractions that are likely to be delivered) may result in a much larger error with a systematic component. The error distribution may be particularly large if a scan at breath-hold is used for planning.     

Effect of respiratory motion on the delivery of breast radiotherapy using SMLC intensity modulation

This study evaluates the effects of respiratory motion on breast radiotherapy delivered using segmented multileaf collimator (SMLC) intensity modulation. An anthropomorphic breast phantom was constructed of polystyrene plates between which radiographic films were inserted. The phantom was mounted on a moving platform to simulate one-dimensional sinusoidal oscillation with variable amplitude and frequency. The motion effect on two breast IMRT techniques, a beamlet-based plan created using the Corvus treatment planning system and an aperture-based plan, was evaluated via film comparison. Motion-induced differences in the treatment region are generally within +5%, with the exception of the posterior field edge and the apex of the breast in the Corvus IMRT plan. Considering the experimental uncertainty arising from the setup and film dosimetry, this result indicates that respiratory motion-induced dose variations are generally relatively insignificant. It appears that the anterior hot spots observed in the Corvus IMRT plan result from the high intensity fluence delivered to the "virtual bolus" area which must be created during the planning process in order to provide "flash" for the respiratory motion. The potential magnitude of such effects resulting from the interplay between fluence modulation and patient motion are unique to the individual planning system and planning technique, as well as the delivery equipment and technique. Such effects should be carefully investigated prior to the implementation of IMRT for breast radiotherapy.     

Dosimetric investigation of lung tumor motion compensation with a robotic respiratory tracking system: an experimental study

The benefits of using Synchrony Respiratory Tracking System (RTS) in conjunction with the CyberKnife robotic treatment device to treat a "breathing tumor" in an anthropomorphic, tissue-equivalent, thoracic phantom have been investigated. The following have been studied: (a) Synchrony's ability to allow the CyberKnife to deliver accurately a planned dose distribution to the free-breathing phantom and (b) the dosimetric implications when irregularities in the breathing cycle and phase differences between internal (tumor) and external (chest) motion exist in the course of one treatment fraction. The breathing phantom PULMONE (phantom used in lung motion experiments) has been used, which can imitate regular or irregular breathing patterns. The breathing traces from two patients with lung cancer have been selected as input. Both traces were irregular in amplitude, frequency, and base line. Patient B demonstrated a phase difference between internal and external motion, whereas patient A did not. The experiment was divided into three stages: In stage I-static, the treatment was delivered to the static phantom. In stage II-motion, the phantom was set to breathe, following the breathing trace of each of the two patients. Synchrony was switched off, so no motion compensation was made. In stage III-compensation, the phantom was set to breathe and Synchrony was switched on. A linear correspondence model was chosen to allow for phase differences between internal and external motion. Gafchromic EBT film was inserted in the phantom tumor to measure dose. To eradicate small errors in film alignment during readout, a gamma comparison with pass criteria of 3%/3 mm was selected. For a more quantitative approach, the percentage of pixels in each gamma map that exceeded the value of 1 (P1) was also used. For both breathing signals, the dose blurring caused by the respiratory motion of the tumor in stage II was degraded considerably compared with stage I (P1 = 15% for patient A and 8% for patient B). The motion compensation via the linear correspondence model was sufficient to provide a dose distribution that satisfied the set gamma criteria (P1=3% for patient A and 2% for patient B). Synchrony RTS has been found satisfactory in recovering the initial detail in dose distribution, for realistic breathing signals, even in the case where a phase delay between internal tumor motion and external chest displacement exists. For the signals applied here, a linear correspondence model provided an acceptable degree of motion compensation.     

A simulation study of irregular respiratory motion and its dosimetric impact on lung tumors

This study is aimed at providing a dosimetric evaluation of the irregular motion of lung tumors due to variations in patients' respiration. Twenty-three lung cancer patients are retrospectively enrolled in this study. The motion of the patient clinical target volume is simulated and two types of irregularities are defined: characteristic and uncharacteristic motions. Characteristic irregularities are representative of random fluctuations in the observed target motion. Uncharacteristic irregular motion is classified as systematic errors in determination of the target motion during the planning session. Respiratory traces from measurement of patient abdominal motion are also used for the target motion simulations. Characteristic irregular motion was observed to cause minimal changes in target dosimetry with the largest effect of 2.5% ± 0.9% (1σ) reduction in the minimum target dose (D(min)) observed for targets that move 2 cm on average and exhibiting 50% amplitude variations within a session. However, uncharacteristic irregular motion introduced more drastic changes in the clinical target volume (CTV) dose; 4.1% ± 1.7% reduction for 1 cm motion and 9.6% ± 1.7% drop for 2 cm. In simulations with patients' abdominal motion, corresponding changes in target dosimetry were observed to be negligible (<0.1%). Only uncharacteristic irregular motion was identified as a clinically significant source of dosimetric uncertainty.     

The impact of MLC transmitted radiation on EPID dosimetry for dynamic MLC beams

The purpose of this study was to experimentally quantify the change in response of an amorphous silicon (a-Si) electronic portal imaging device (EPID) to dynamic multileaf collimator (dMLC) beams with varying MLC-transmitted dose components and incorporate the response into a commercial treatment planning system (TPS) EPID prediction model. A combination of uniform intensity dMLC beams and static beams were designed to quantify the effect of MLC transmission on EPID response at the central axis of 10 x 10 cm2 beams, at off-axis positions using wide dMLC beam profiles, and at different field sizes. The EPID response to MLC transmitted radiation was 0.79 +/- 0.02 of the response to open beam radiation at the central axis of a 10 x 10 cm2 field. The EPID response to MLC transmitted radiation was further reduced relative to the open beam response with off-axis distance. The EPID response was more sensitive to field size changes for MLC transmitted radiation compared to open beam radiation by a factor of up to 1.17 at large field sizes. The results were used to create EPID response correction factors as a function of the fraction of MLC transmitted radiation, off-axis distance, and field size. Software was developed to apply the correction factors to each pixel in the TPS predicted EPID image. The corrected images agreed more closely with the measured EPID images in areas of intensity modulated fields with a large fraction of MLC transmission and, as a result the accuracy of portal dosimetry with a-Si EPIDs can be improved. Further investigation into the detector response function and the radiation source model are required to achieve improvements in accuracy for the general case.     

Fluoroscopic study of tumor motion due to breathing: facilitating precise radiation therapy for lung cancer patients

Target motion due to breathing is one of the major obstacles in dose escalation of radiation therapy to some tumors in the thoracoabdominal region. The development of beam gating or target motion tracking techniques provides a possibility to reduce normal tissue volume in a treatment field. Tumor motion monitoring in those techniques plays a crucial role, but has not yet been adequately explored. This paper reports our preliminary investigation on breath introduced tumor motion. Tumor locations and motion properties were determined from digitized fluoroscopic videos acquired during patient simulation. Image distortion due to irregularities in the imaging chain, such as the pincushion distortion, was corrected with a polynomial unwarping method. Temporal Fourier transformation of the fluoroscopic video was introduced to convert the motion information over time to a static view of a motion field, in which regions with different motion ranges can be directly measured. Patient breathing patterns vary from patient to patient and so does the kinematic behavior of individual tumors. In order to evaluate the feasibility for tracking internal target motion with nonionizing-radiation techniques, motion patterns between internal targets and external radio opaque markers placed on patient's chest during fluoroscopic video acquisition were compared. For some patients, significant motion phase discrepancies between an internal target and an external marker have been observed. Quantitative measurements are reported. These results will be useful in the design of a motion tracking or gated radiotherapy system.