Streaking artifacts reduction in four-dimensional cone-beam computed tomography

Cone-beam computed tomography (CBCT) using an "on-board" x-ray imaging device integrated into a radiation therapy system has recently been made available for patient positioning, target localization, and adaptive treatment planning. One of the challenges for gantry mounted image-guided radiation therapy (IGRT) systems is the slow acquisition of projections for cone-beam CT (CBCT), which makes them sensitive to any patient motion during the scans. Aiming at motion artifact reduction, four-dimensional CBCT (4D CBCT) techniques have been introduced, where a surrogate for the target's motion profile is utilized to sort the cone-beam data by respiratory phase. However, due to the limited gantry rotation speed and limited readout speed of the on-board imager, fewer than 100 projections are available for the image reconstruction at each respiratory phase. Thus, severe undersampling streaking artifacts plague 4D CBCT images. In this paper, the authors propose a simple scheme to significantly reduce the streaking artifacts. In this method, a prior image is first reconstructed using all available projections without gating, in which static structures are well reconstructed while moving objects are blurred. The undersampling streaking artifacts from static structures are estimated from this prior image volume and then can be removed from the phase images using gated reconstruction. The proposed method was validated using numerical simulations, experimental phantom data, and patient data. The fidelity of stationary and moving objects is maintained, while large gains in streak artifact reduction are observed. Using this technique one can reconstruct 4D CBCT datasets using no more projections than are acquired in a 60 s scan. At the same time, a temporal gating window as narrow as 100 ms was utilized. Compared to the conventional 4D CBCT reconstruction, streaking artifacts were reduced by 60% to 70%.     

Dosimetric evaluation of radiation dose rate effect in respiratory gated intensity modulated radiation therapy

Purpose:

To investigate the dosimetric accuracy of the sliding window gated IMRT compared with the static treatment, using varying dose rates.

Materials and methods:

This study measured changes in output and diode array response with changing dose rate, verified the precision of the motion table, and measured changes in dose distribution accuracy with film and diodes at two depths with changing dose rate. During 4DCT (4 Dimensional Computed Tomography), the patient’s respiratory signals and target motion were recorded and imported to the XY4D simulation table of SUN NUCLEAR Corporation to simulate the patient’s respiration and tumour motion. A single field of each sliding window IMRT plan with 30º wedge and one for lung cancer were used in this study. Three irradiating conditions, static and moving target with and without gating, were applied to both plans.

Results:

The standard deviations of output, with the dose rates changing from 300–600 MU/min, were 0.065 cGy and 0.169 cGy for the ionisation chamber and diode, respectively. The verification of the motion table shows very good precision with 9.98 ± 0.02 cm (true value = 10.0 cm). The measurements by MapCheck show the gamma index of the planned absolute dose distribution in static and moving targets with gating, resulting in more than 96% passing for all dose rates. The absolute dose distribution measured by film for the static target was agreeable with the value of moving target with gating.

Conclusion:

The sliding window gated IMRT technique is able to deliver an accurate dose to a moving target with the dose rate of 300–600 MU/min that is suitable for clinical treatment.     

Respiratory correlated cone-beam computed tomography on an isocentric C-arm

A methodology for 3D image reconstruction from retrospectively gated cone-beam CT projection data has been developed. A mobile x-ray cone-beam device consisting of an isocentric C-arm equipped with a flat panel detector was used to image a moving phantom. Frames for reconstruction were retrospectively selected from complete datasets based on the known rotation of the C-arm and a signal from a respiratory monitor. Different sizes of gating windows were tested. A numerical criterion for blur on the reconstructed image was suggested. The criterion is based on minimization of an Ising energy function, similar to approaches used in image segmentation or restoration. It is shown that this criterion can be used for the determination of the optimal gating window size. Images reconstructed from the retrospectively gated projection sequences using the optimal gating window data showed a significant improvement compared to images reconstructed from the complete projection datasets.     

The effects of tumor motion on planning and delivery of respiratory-gated IMRT

The purpose of this study is to investigate the effects of object motion on the planning and delivery of IMRT. Two phantoms containing objects were imaged using CT under a variety of motion conditions. The effects of object motion on axial CT acquisition with and without gating were assessed qualitatively and quantitatively. Measurements of effective slice width and position for the CT scans were made. Mutual information image fusion was adapted for use as a quantitative measure of object deformation in CT images. IMRT plans were generated on the CT scans of the moving and gated object images. These plans were delivered with motion, with and without gating, and the delivery error between the moving deliveries and a nonmoving delivery was assessed using a scalable vector-based index. Motion during CT acquisition produces motion artifact, object deformation, and object mispositioning, which can be substantially reduced with gating. Objects that vary in cross section in the direction of motion exhibit the most deformation in CT images. Mutual information provides a useful quantitative estimate of object deformation. The delivery of IMRT in the presence of target motion significantly alters the delivered dose distribution in relation to the planned distribution. The utilization of gating for IMRT treatment, including imaging, planning, and delivery, significantly reduces the errors introduced by object motion.     

Evaluation of the combined effects of target size, respiratory motion and background activity on 3D and 4D PET/CT images

Gated (4D) PET/CT has the potential to greatly improve the accuracy of radiotherapy at treatment sites where internal organ motion is significant. However, the best methodology for applying 4D-PET/CT to target definition is not currently well established. With the goal of better understanding how to best apply 4D information to radiotherapy, initial studies were performed to investigate the effect of target size, respiratory motion and target-to-background activity concentration ratio (TBR) on 3D (ungated) and 4D PET images. Using a PET/CT scanner with 4D or gating capability, a full 3D-PET scan corrected with a 3D attenuation map from 3D-CT scan and a respiratory gated (4D) PET scan corrected with corresponding attenuation maps from 4D-CT were performed by imaging spherical targets (0.5-26.5 mL) filled with (18)F-FDG in a dynamic thorax phantom and NEMA IEC body phantom at different TBRs (infinite, 8 and 4). To simulate respiratory motion, the phantoms were driven sinusoidally in the superior-inferior direction with amplitudes of 0, 1 and 2 cm and a period of 4.5 s. Recovery coefficients were determined on PET images. In addition, gating methods using different numbers of gating bins (1-20 bins) were evaluated with image noise and temporal resolution. For evaluation, volume recovery coefficient, signal-to-noise ratio and contrast-to-noise ratio were calculated as a function of the number of gating bins. Moreover, the optimum thresholds which give accurate moving target volumes were obtained for 3D and 4D images. The partial volume effect and signal loss in the 3D-PET images due to the limited PET resolution and the respiratory motion, respectively were measured. The results show that signal loss depends on both the amplitude and pattern of respiratory motion. However, the 4D-PET successfully recovers most of the loss induced by the respiratory motion. The 5-bin gating method gives the best temporal resolution with acceptable image noise. The results based on the 4D scan protocols can be used to improve the accuracy of determining the gross tumor volume for tumors in the lung and abdomen.     

Statistical analysis of target motion in gated lung stereotactic body radiation therapy

An external surrogate-based respiratory gating technique is a useful method to reduce target margins for the treatment of a moving lung tumor. The success of this technique relies on a good correlation between the motion of the external markers and the internal tumor as well as the repeatability of the respiratory motion. In gated lung stereotactic body radiation therapy (SBRT), the treatment time for each fraction could exceed 30 min due to large fractional dose. Tumor motion may experience pattern changes such as baseline shift during such extended treatment time. The purpose of this study is to analyze tumor motion traces in actual treatment situations and to evaluate the effect of the target baseline shift in gated lung SBRT treatment. Real-time motion data for both the external markers and tumors from 51 lung SBRT treatments with Cyberknife Synchrony technology were analyzed in this study. The treatment time is typically greater than 30 min. The baseline shift was calculated with a rolling average window equivalent to ～20 s and subtracted from that at the beginning. The magnitude of the baseline shift and its relationship with treatment time were investigated. Phase gating simulation was retrospectively performed on 12 carefully selected treatments with respiratory amplitude larger than 5 mm and regular phases. A customized gating window was defined for each individual treatment. It was found that the baseline shifts are specific to each patient and each fraction. Statistical analysis revealed that more than 69% treatments exhibited increased baseline shifts with the lapse of treatment time. The magnitude of the baseline shift could reach 5.3 mm during a 30 min treatment. Gating simulation showed that tumor excursion was caused mainly by the uncertainties in phase gating simulation and baseline shift, the latter being the primary factor. With a 5 mm gating window, 2 out of 12 treatments in the study group showed significant tumor excursion. Baseline shifts alone could cause up to 20% of tumor excursion outside the gating window. It is concluded that baseline shifts may increase with the treatment time and are more likely to act as a time-dependent systematic error. For phase-based gated lung SBRT, a baseline shift may be one of the major sources of targeting error during treatment.     

Validation of target volume and position in respiratory gated CT planning and treatment

The capability of a commercial respiratory gating system based on video tracking of reflective markers to reduce motion-induced CT planning and treatment errors was evaluated. Spherical plastic shells (2.8-82 cm3), simulating the gross target volume (GTV), were placed in a water-filled body phantom that was moved sinusoidally along the longitudinal axis of the CT scanner and the accelerator for +/- 1 cm at 15-30 cycle/min. During gated CT imaging, the x-ray exposure was initiated by the gating system shortly before the end of expiration (so that the imaging time would be centered at the end of expiration); it was terminated by the scanner after completion of each slice. In nongated CT images, the target appeared distorted and often broken up. GTVs volume errors ranged 16%-110% in axial scans, and 7%-36% in spiral scans. In gated CT images, the spheres appeared 3 and 5 mm longer than their actual diameters (volume errors 2%-16%), at the respective respiration rates of 15 and 20 cycles/min. At 30 cycles/min the target appeared 1 cm longer, and volume error ranged 25%-53%. During treatment, gating kept the beam on for a duration equal to the CT acquisition time of 1 s/slice. The difference in positional errors between gated CT and portal films was 1 mm, regardless the size of residual motion errors. Because of the potential of suboptimal placement of the gating window between CT imaging and treatment, an extra 1.5-2.5 mm safety margin can be added regardless of the size of residual motion error. For respiratory rates > or = 30 cycles/min, the effectiveness of gating is limited by large residual motion in the 1 s CT acquisition time.     

Targeting accuracy of an image guided gating system for stereotactic body radiotherapy

Recently, a commercial system capable of x-ray image guided patient positioning and respiratory gated delivery has become available. Here we describe the operational principles of this system and investigate its geometric targeting accuracy under controlled conditions. The system tracks breathing via infrared (IR) detection of reflective markers located on the patient's abdomen. Localization kilovoltage (kV) x-rays are triggered from within the gated delivery window portion of the breathing trace and after positioning, the tumour will cross the linac isocentre during gated delivery. We tested geometric accuracy of this system by localizing and delivering gated fields to a moving phantom. Effects of phantom speed, gating window location, timing errors and phantom rotations on positioning and gating accuracy were investigated. The system delivered gated fields to both a moving and static phantom with equal accuracy. The position of the gating window affects accuracy only to the extent that an asymmetric breathing motion could affect dose distribution within its boundaries. Positioning errors were found to be less then 0.5 +/- 0.2 mm for phantom rotations up to 5 degrees. We found and corrected a synchronization error caused by a faulty x-ray duration setting and detected a 60 +/- 20 ms time delay in our linear accelerator.     

A technique for respiratory-gated radiotherapy treatment verification with an EPID in cine mode

Respiratory gating based on external surrogates is performed in many clinics. We have developed a new technique for treatment verification using an electronic portal imaging device (EPID) in cine mode for gated 3D conformal therapy. Implanted radiopaque fiducial markers inside or near the target are required for this technique. The markers are contoured on the planning CT set, enabling us to create digitally reconstructed radiographs (DRRs) for each treatment beam. During the treatment, a sequence of EPID images can be acquired without disrupting the treatment. Implanted markers are visualized in the images and their positions in the beam's eye view are calculated off-line and compared to the reference position by matching the field apertures in corresponding EPID and DRR images. The precision of the patient set-up, the placement of the beam-gating window, as well as the residual tumour motion can be assessed for each treatment fraction. This technique has been demonstrated with a case study patient, who had three markers implanted in his liver. For this patient, the intra-fractional variation of all marker positions in the gating window had a 95% range of 4.8 mm in the SI direction (the primary axis of motion). This was about the same (5 mm) as the residual motion considered in the planning process. The inter-fractional variation of the daily mean positions of the markers, which indicates the uncertainty in the set-up procedure, was within +8.3 mm/-4.5 mm (95% range) in the SI direction for this case.     

An evaluation of gating window size,delivery method,and composite field dosimetry of respiratory-gated IMRT

A respiratory gating system has been developed based on a commercial patient positioning system. The purpose of this study is to investigate the ability of the gating system to reproduce normal, nongated IMRT operation and to quantify the errors produced by delivering a nongated IMRT treatment onto a moving target. A moving phantom capable of simultaneous two-dimensional motion was built, and an analytical liver motion function was used to drive the phantom. Studies were performed to assess the effect of gating window size and choice of delivery method (segmented and dynamic multileaf collimation). Additionally, two multiple field IMRT cases were delivered to quantify the error in gated and nongated IMRT with motion. Dosimetric error between nonmoving and moving deliveries is related to gating window size. By reducing the window size, the error can be reduced. Delivery error can be reduced for both dynamic and segmented delivery with gating. For the implementation of dynamic IMRT delivery in this study, dynamic delivery was found to generate larger delivery errors than segmented delivery in most cases studied. For multiple field IMRT delivery, the largest errors were generated in regions where high field modulation was present parallel to the axis of motion. Gating was found to reduce these large errors to clinically acceptable levels.     

Extension of a data-driven gating technique to 3D, whole body PET studies

Respiratory gating can be used to separate a PET acquisition into a series of near motion-free bins. This is typically done using additional gating hardware; however, software-based methods can derive the respiratory signal from the acquired data itself. The aim of this work was to extend a data-driven respiratory gating method to acquire gated, 3D, whole body PET images of clinical patients. The existing method, previously demonstrated with 2D, single bed-position data, uses a spectral analysis to find regions in raw PET data which are subject to respiratory motion. The change in counts over time within these regions is then used to estimate the respiratory signal of the patient. In this work, the gating method was adapted to only accept lines of response from a reduced set of axial angles, and the respiratory frequency derived from the lung bed position was used to help identify the respiratory frequency in all other bed positions. As the respiratory signal does not identify the direction of motion, a registration-based technique was developed to align the direction for all bed positions. Data from 11 clinical FDG PET patients were acquired, and an optical respiratory monitor was used to provide a hardware-based signal for comparison. All data were gated using both the data-driven and hardware methods, and reconstructed. The centre of mass of manually defined regions on gated images was calculated, and the overall displacement was defined as the change in the centre of mass between the first and last gates. The mean displacement was 10.3 mm for the data-driven gated images and 9.1 mm for the hardware gated images. No significant difference was found between the two gating methods when comparing the displacement values. The adapted data-driven gating method was demonstrated to successfully produce respiratory gated, 3D, whole body, clinical PET acquisitions.     

Gated CT imaging using a free-breathing respiration signal from flow-volume spirometry

Respiration-induced tumor motion is known to cause artifacts on free-breathing spiral CT images used in treatment planning. This leads to inaccurate delineation of target volumes on planning CT images. Flow-volume spirometry has been used previously for breath-holds during CT scans and radiation treatments using the active breathing control (ABC) system. We have developed a prototype by extending the flow-volume spirometer device to obtain gated CT scans using a PQ 5000 single-slice CT scanner. To test our prototype, we designed motion phantoms to compare image quality obtained with and without gated CT scan acquisition. Spiral and axial (nongated and gated) CT scans were obtained of phantoms with motion periods of 3-5 s and amplitudes of 0.5-2 cm. Errors observed in the volume estimate of these structures were as much as 30% with moving phantoms during CT simulation. Application of motion-gated CT with active breathing control reduced these errors to within 5%. Motion-gated CT was then implemented in patients and the results are presented for two clinical cases: lung and abdomen. In each case, gated scans were acquired at end-inhalation, end-exhalation in addition to a conventional free-breathing (nongated) scan. The gated CT scans revealed reduced artifacts compared with the conventional free-breathing scan. Differences of up to 20% in the volume of the structures were observed between gated and free-breathing scans. A comparison of the overlap of structures between the gated and free-breathing scans revealed misalignment of the structures. These results demonstrate the ability of flow-volume spirometry to reduce errors in target volumes via gating during CT imaging.     

Linac-integrated 4D cone beam CT: first experimental results

A new online imaging approach, linac-integrated cone beam CT (CBCT), has been developed over the past few years. It has the advantage that a patient can be examined in their treatment position directly before or during a radiotherapy treatment. Unfortunately, respiratory organ motion, one of the largest intrafractional organ motions, often leads to artefacts in the reconstructed 3D images. One way to take this into account is to register the breathing phase during image acquisition for a phase-correlated image reconstruction. Therefore, the main focus of this work is to present a system which has the potential to investigate the correlation between internal (movement of the diaphragm) and external (data of a respiratory gating system) information about breathing phase and amplitude using an inline CBCT scanner. This also includes a feasibility study about using the acquired information for a respiratory-correlated 4D CBCT reconstruction. First, a moving lung phantom was used to develop and to specify the required methods which are based on an image reconstruction using only projections belonging to a certain moving phase. For that purpose, the corresponding phase has to be detected for each projection. In the case of the phantom, an electrical signal allows one to track the movement in real time. The number of projections available for the image reconstruction depends on the breathing phase and the size of the position range from which projections should be used for the reconstruction. The narrower this range is, the better the inner structures can be located, but also the noise of the images increases due to the limited number of projections. This correlation has also been analysed. In a second step, the methods were clinically applied using data sets of patients with lung tumours. In this case, the breathing phase was detected by an external gating system (AZ-733V, Anzai Medical Co.) based on a pressure sensor attached to the patient's abdominal region with a fixation belt. The comparison of the reconstructed 4D CBCT images and the corresponding 4D CT images used for the treatment planning provides the required information for the calculation of possible setup errors. So, a repositioning of the patient is feasible even though the patient moves due to respiration. In addition to the external signal, the position of the diaphragm in the cranial-caudal direction could be extracted from each projection. Both independent sources of information show a very good agreement of the phase and even the amplitude of the movement and the external signal respectively. This suggests the usability of such a system for a gated dose delivery approach. However, more studies involving patients with different incidences have to be carried out to confirm these first results.     

Automatic segmentation of phase-correlated CT scans through nonrigid image registration using geometrically regularized free-form deformation

Conventional radiotherapy is planned using free-breathing computed tomography (CT), ignoring the motion and deformation of the anatomy from respiration. New breath-hold-synchronized, gated, and four-dimensional (4D) CT acquisition strategies are enabling radiotherapy planning utilizing a set of CT scans belonging to different phases of the breathing cycle. Such 4D treatment planning relies on the availability of tumor and organ contours in all phases. The current practice of manual segmentation is impractical for 4D CT, because it is time consuming and tedious. A viable solution is registration-based segmentation, through which contours provided by an expert for a particular phase are propagated to all other phases while accounting for phase-to-phase motion and anatomical deformation. Deformable image registration is central to this task, and a free-form deformation-based nonrigid image registration algorithm will be presented. Compared with the original algorithm, this version uses novel, computationally simpler geometric constraints to preserve the topology of the dense control-point grid used to represent free-form deformation and prevent tissue fold-over. Using mean squared difference as an image similarity criterion, the inhale phase is registered to the exhale phase of lung CT scans of five patients and of characteristically low-contrast abdominal CT scans of four patients. In addition, using expert contours for the inhale phase, the corresponding contours were automatically generated for the exhale phase. The accuracy of the segmentation (and hence deformable image registration) was judged by comparing automatically segmented contours with expert contours traced directly in the exhale phase scan using three metrics: volume overlap index, root mean square distance, and Hausdorff distance. The accuracy of the segmentation (in terms of radial distance mismatch) was approximately 2 mm in the thorax and 3 mm in the abdomen, which compares favorably to the accuracies reported elsewhere. Unlike most prior work, segmentation of the tumor is also presented. The clinical implementation of 4D treatment planning is critically dependent on automatic segmentation, for which is offered one of the most accurate algorithms yet presented.     

An automatic contour propagation method to follow parotid gland deformation during head-and-neck cancer tomotherapy

We developed an efficient technique to auto-propagate parotid gland contours from planning kVCT to daily MVCT images of head-and-neck cancer patients treated with helical tomotherapy. The method deformed a 3D surface mesh constructed from manual kVCT contours by B-spline free-form deformation to generate optimal and smooth contours. Deformation was calculated by elastic image registration between kVCT and MVCT images. Data from ten head-and-neck cancer patients were considered and manual contours by three observers were included in both kVCT and MVCT images. A preliminary inter-observer variability analysis demonstrated the importance of contour propagation in tomotherapy application: a high variability was reported in MVCT parotid volume estimation (p = 0.0176, ANOVA test) and a larger uncertainty of MVCT contouring compared with kVCT was demonstrated by DICE and volume variability indices (Wilcoxon signed rank test, p < 10(-4) for both indices). The performance analysis of our method showed no significant differences between automatic and manual contours in terms of volumes (p > 0.05, in a multiple comparison Tukey test), center-of-mass distances (p = 0.3043, ANOVA test), DICE values (p = 0.1672, Wilcoxon signed rank test) and average and maximum symmetric distances (p = 0.2043, p = 0.8228 Wilcoxon signed rank tests). Results suggested that our contour propagation method could successfully substitute human contouring on MVCT images.     

Quantification of accuracy of the automated nonlinear image matching and anatomical labeling (ANIMAL) nonlinear registration algorithm for 4D CT images of lung

The performance of the ANIMAL (Automated Nonlinear Image Matching and Anatomical Labeling) nonlinear registration algorithm for registration of thoracic 4D CT images was investigated. The algorithm was modified to minimize the incidence of deformation vector discontinuities that occur during the registration of lung images. Registrations were performed between the inhale and exhale phases for five patients. The registration accuracy was quantified by the cross-correlation of transformed and target images and distance to agreement (DTA) measured based on anatomical landmarks and triangulated surfaces constructed from manual contours. On average, the vector DTA between transformed and target landmarks was 1.6 mm. Comparing transformed and target 3D triangulated surfaces derived from planning contours, the average target volume (GTV) center-of-mass shift was 2.0 mm and the 3D DTA was 1.6 mm. An average DTA of 1.8 mm was obtained for all planning structures. All DTA metrics were comparable to inter observer uncertainties established for landmark identification and manual contouring.     

Feasibility of four-dimensional conformal planning for robotic radiosurgery

Organ motion can have a severe impact on the dose delivered by radiation therapy, and different procedures have been developed to address its effects. Conventional techniques include breath hold methods and gating. A different approach is the compensation for target motion by moving the treatment beams synchronously. Practical results have been reported for robot based radiosurgery, where a linear accelerator mounted on a robotic arm delivers the dose. However, not all organs move in the same way, which results in a relative motion of the beams with respect to the body and the tissues in the proximity of the tumor. This relative motion can severely effect the dose delivered to critical structures. We propose a method to incorporate motion in the treatment planning for robotic radiosurgery to avoid potential overdosing of organs surrounding the target. The method takes into account the motion of all considered volumes, which is discretized for dose calculations. Similarly, the beam motion is taken into account and the aggregated dose coefficient over all discrete steps is used for planning. We simulated the treatment of a moving target with three different planning methods. First, we computed beam weights based on a 3D planning situation and simulated treatment with organ motion and the beams moving synchronously to the target. Second, beam weights were computed by the 4D planning method incorporating the organ and beam motion and treatment was simulated for beams moving synchronously to the target. Third, the beam weights were determined by the 4D planning method with the beams fixed during planning and simulation. For comparison we also give results for the 3D treatment plan if there was no organ motion and when the plan is delivered by fixed beams in the presence of organ motion. The results indicate that the new 4D method is preferable and can further improve the overall conformality of motion compensated robotic radiosurgery.     

Automatic re-contouring in 4D radiotherapy

Delineating regions of interest (ROIs) on each phase of four-dimensional (4D) computed tomography (CT) images is an essential step for 4D radiotherapy. The requirement of manual phase-by-phase contouring prohibits the routine use of 4D radiotherapy. This paper develops an automatic re-contouring algorithm that combines techniques of deformable registration and surface construction. ROIs are manually contoured slice-by-slice in the reference phase image. A reference surface is constructed based on these reference contours using a triangulated surface construction technique. The deformable registration technique provides the voxel-to-voxel mapping between the reference phase and the test phase. The vertices of the reference surface are displaced in accordance with the deformation map, resulting in a deformed surface. The new contours are reconstructed by cutting the deformed surface slice-by-slice along the transversal, sagittal or coronal direction. Since both the inputs and outputs of our automatic re-contouring algorithm are contours, it is relatively easy to cope with any treatment planning system. We tested our automatic re-contouring algorithm using a deformable phantom and 4D CT images of six lung cancer patients. The proposed algorithm is validated by visual inspections and quantitative comparisons of the automatic re-contours with both the gold standard segmentations and the manual contours. Based on the automatic delineated ROIs, changes of tumour and sensitive structures during respiration are quantitatively analysed. This algorithm could also be used to re-contour daily images for treatment evaluation and adaptive radiotherapy.     

Residual motion of lung tumors in end-of-inhale respiratory gated radiotherapy based on external surrogates

It has been noted that some lung tumors exhibit large periodic motion due to respiration. To limit the amount of dose to healthy lung tissues, many clinics have begun gating radiotherapy treatment using externally placed surrogates. It has been observed by several institutions that the end-of-exhale (EOE) tumor position is more reproducible than other phases of the breathing cycle, so the gating window is often set there. From a treatment planning perspective, end-of-inhale (EOI) phase might be preferred for gating because the expanded lungs will further decrease the healthy tissue within the treatment field. We simulate gated treatment at the EOI phase, using a set of recently measured internal/external anatomy patient data. This paper attempts to answer three questions: (1) How much is the tumor residual motion when we use an external surrogate gating window at EOI? (2) How could we reduce the residual motion in the EOI gating window? (3) Is there a preference for amplitude- versus phase-based gating at EOI? We found that under free breathing conditions the residual motion of the tumors is much larger for EOI phase than for EOE phase. The mean values of residual motion at EOI were found to be 2.2 and 2.7 mm for amplitude- and phase-based gating, respectively, and, at EOE, 1.0 and 1.2 mm for amplitude- and phase-based gating, respectively. However, we note that the residual motion in the EOI gating window is correlated well with the reproducibility of the external surface position in the EOI phase. Using the results of a published breath-coaching study, we deduce that the residual motion of a lung tumor at EOI would approach that at EOE, with the same duty cycle (30%), under breath-coaching conditions. Additionally, we found that under these same conditions, phase-based gating approaches the same residual motion as amplitude-based gating, going from a 28% difference to 11%, for the patient with the largest difference between the two gating modalities. We conclude that it is feasible to achieve the same reproducibility of tumor location at EOI as at EOE if breath coaching is implemented, enabling us to reap the benefits of the dosimetric advantage of EOI gating.