Dosimetric impact of motion in free-breathing and gated lung radiotherapy: a 4D Monte Carlo study of intrafraction and interfraction effects

Current Research Status of Respiratory Motion for Thorax and Abdominal Treatment: A Systematic Review

Malignant tumors have become one of the serious public health problems in human safety and health, among which the chest and abdomen diseases account for the largest proportion. Early diagnosis and treatment can effectively improve the survival rate of patients. However, respiratory motion in the chest and abdomen can lead to uncertainty in the shape, volume, and location of the tumor, making treatment of the chest and abdomen difficult. Therefore, compensation for respiratory motion is very important in clinical treatment. The purpose of this review was to discuss the research and development of respiratory movement monitoring and prediction in thoracic and abdominal surgery, as well as introduce the current research status. The integration of modern respiratory motion compensation technology with advanced sensor detection technology, medical-image-guided therapy, and artificial intelligence technology is discussed and analyzed. The future research direction of intraoperative thoracic and abdominal respiratory motion compensation should be non-invasive, non-contact, use a low dose, and involve intelligent development. The complexity of the surgical environment, the constraints on the accuracy of existing image guidance devices, and the latency of data transmission are all present technical challenges.

Anthropomorphic phantom for deformable lung and liver CT and MR imaging for radiotherapy

In this study, a functioning and ventilated anthropomorphic phantom was further enhanced for the purpose of CT and MR imaging of the lung and liver. A deformable lung, including respiratory tract was 3D printed. Within the lung's inner structures is a solid region shaped from a patient's lung tumour and six nitro-glycerine capsules as reference landmarks. The full internal mesh was coated, and the tumour filled, with polyorganosiloxane based gel. A moulded liver was created with an external casing of silicon filled with polyorganosiloxane gel and flexible plastic internal structures. The liver, fitted to the inferior portion of the right lung, moves along with the lung's ventilation. In the contralateral side, a cavity is designed to host a dosimeter, whose motion is correlated to the lung pressure. A 4DCT of the phantom was performed along with static and 4D T1 weighted MR images. The CT Hounsfield units (HU) for the flexible 3D printed material were -600-100 HU (lung and liver structures), for the polyorganosiloxane gel 30-120 HU (lung coating and liver filling) and for the silicon 650-800 HU (liver casing). The MR image intensity units were 0-40, 210-280 and 80-130, respectively. The maximum range of motion in the 4D imaging for the superior lung was 1-3.5 mm and 3.5-8 mm in the inferior portion. The liver motion was 5.5-8.0 mm at the tip and 5.7-10.0 mm at the dome. No measurable drift in motion was observed over a 2 h session and motion was reproducible over three different sessions for sin²(t), sin⁴(t) and a patient-like breathing curve with the interquartile range of amplitudes for all breathing cycles within 0.5 mm. The addition of features within the lung and of a deformable liver will allow the phantom to be used for imaging studies such as validation of 4DMRI and pseudo CT methods.

4D Monte Carlo dose calculations for pre-treatment quality assurance of VMAT SBRT: a phantom-based feasibility study

Volumetric arc therapy (VMAT) for lung stereotactic body radiotherapy (SBRT) is challenging due to both breathing-induced motion and the dynamic components of the linear accelerator. In this study, a 4D Monte Carlo (4DMC) dose calculation method for VMAT SBRT is proposed and the feasibility of the method is evaluated. A rigidly-moving lung phantom was imaged using four dimensional computed tomography (4DCT). VMAT SBRT plans were generated on the average intensity projection dataset using the internal target volume (ITV) strategy (ITV-plan) and a single phase to simulate a dynamic treatment-couch tracking technique (TRACKING-plan). 4DMC simulations were performed and compared to 3D Monte Carlo (3DMC) and 3D- and 4D- calculations in the treatment planning system using the adaptive convolution (AC) algorithm. Dose metrics calculated for the ITV-plan showed an overestimation with 3D adaptive convolution (3DAC) for D[Formula: see text] (GTV) by 3.5% and by 2.0% for 3DMC, both compared to 4DMC. The TRACKING-plan D[Formula: see text] (GTV) calculated with the 3DAC method overestimated by 2.0% compared with 4DMC. Deviations between the calculation methods for D _mean (Lung) and D[Formula: see text] (PTV) were minimal. For both plans, measurements were taken with EBT3 film inside the phantom tumour. EBT3 film profiles showed good agreement with 4DMC for the TRACKING-plan giving a gamma pass rate of 97.2% for 3%/3 mm global and for 3DAC compared with measured, 95.8%. Whereas for the ITV-plan, the 3D profiles varied from film in the ITV periphery region with a pass rates of 50% and 48.6% for 3DAC and 3DMC, respectively. 4DMC agreed more closely to measurements for this plan with a pass rate of 95.8%. We have proposed an accurate method to perform 4D dose calculations for pre-treatment quality assurance of VMAT SBRT. The method was compared to experimental measurements and for both plans, 4DMC dose agreed with measurements more closely than other evaluated dose calculation methods. This study has demonstrated the feasibility of this 4DMC method.

GIFTed Demons: deformable image registration with local structure-preserving regularization using supervoxels for liver applications

Deformable image registration, a key component of motion correction in medical imaging, needs to be efficient and provides plausible spatial transformations that reliably approximate biological aspects of complex human organ motion. Standard approaches, such as Demons registration, mostly use Gaussian regularization for organ motion, which, though computationally efficient, rule out their application to intrinsically more complex organ motions, such as sliding interfaces. We propose regularization of motion based on supervoxels, which provides an integrated discontinuity preserving prior for motions, such as sliding. More precisely, we replace Gaussian smoothing by fast, structure-preserving, guided filtering to provide efficient, locally adaptive regularization of the estimated displacement field. We illustrate the approach by applying it to estimate sliding motions at lung and liver interfaces on challenging four-dimensional computed tomography (CT) and dynamic contrast-enhanced magnetic resonance imaging datasets. The results show that guided filter-based regularization improves the accuracy of lung and liver motion correction as compared to Gaussian smoothing. Furthermore, our framework achieves state-of-the-art results on a publicly available CT liver dataset.

Internal Motion Estimation by Internal-external Motion Modeling for Lung Cancer Radiotherapy

The aim of this study is to develop an internal-external correlation model for internal motion estimation for lung cancer radiotherapy. Deformation vector fields that characterize the internal-external motion are obtained by respectively registering the internal organ meshes and external surface meshes from the 4DCT images via a recently developed local topology preserved non-rigid point matching algorithm. A composite matrix is constructed by combing the estimated internal phasic DVFs with external phasic and directional DVFs. Principle component analysis is then applied to the composite matrix to extract principal motion characteristics, and generate model parameters to correlate the internal-external motion. The proposed model is evaluated on a 4D NURBS-based cardiac-torso (NCAT) synthetic phantom and 4DCT images from five lung cancer patients. For tumor tracking, the center of mass errors of the tracked tumor are 0.8(±0.5)mm/0.8(±0.4)mm for synthetic data, and 1.3(±1.0)mm/1.2(±1.2)mm for patient data in the intra-fraction/inter-fraction tracking, respectively. For lung tracking, the percent errors of the tracked contours are 0.06(±0.02)/0.07(±0.03) for synthetic data, and 0.06(±0.02)/0.06(±0.02) for patient data in the intra-fraction/inter-fraction tracking, respectively. The extensive validations have demonstrated the effectiveness and reliability of the proposed model in motion tracking for both the tumor and the lung in lung cancer radiotherapy.

Effects of collimator angle, couch angle, and starting phase on motion-tracking dynamic conformal arc therapy (4D DCAT)

Purpose

The aim of this study was to find an optimized configuration of collimator angle, couch angle, and starting tracking phase to improve the delivery performance in terms of MLC position errors, maximal MLC leaf speed, and total beam‐on time of DCAT plans with motion tracking (4D DCAT).

Method and materials

Nontracking conformal arc plans were first created based on a single phase (maximal exhalation phase) of a respiratory motion phantom with a spherical target. An ideal model was used to simulate the target motion in superior‐inferior (SI), anterior‐posterior (AP), and left‐right (LR) dimensions. The motion was decomposed to the MLC leaf position coordinates for motion compensation and generating 4D DCAT plans. The plans were studied with collimator angle ranged from 0° to 90°; couch angle ranged from 350°(−10°) to 10°; and starting tracking phases at maximal inhalation (θ=π/2) and exhalation (θ=0) phases. Plan performance score (PPS) evaluates the plan complexity including the variability in MLC leaf positions, degree of irregularity in field shape and area. PPS ranges from 0 to 1, where low PPS indicates a plan with high complexity. The 4D DCAT plans with the maximal and the minimal PPS were selected and delivered on a Varian TrueBeam linear accelerator. Gafchromic‐EBT3 dosimetry films were used to measure the dose delivered to the target in the phantom. Gamma analysis for film measurements with 90% passing rate threshold using 3%/3 mm criteria and trajectory log files were analyzed for plan delivery accuracy evaluation.

Results

The maximal PPS of all the plans was 0.554, achieved with collimator angle at 87°, couch angle at 350°, and starting phase at maximal inhalation (θ=π/2). The maximal MLC leaf speed, MLC leaf errors, total leaf travel distance, and beam‐on time were 20 mm/s, 0.39 ± 0.16 mm, 1385 cm, and 157 s, respectively. The starting phase, whether at maximal inhalation or exhalation had a relatively small contribution to PPS (0.01 ± 0.05).

Conclusions

By selecting collimator angle, couch angle, and starting tracking phase, 4D DCAT plans with the maximal PPS demonstrated less MLC leaf position errors, lower maximal MLC leaf speed, and shorter beam‐on time which improved the performance of 4D motion‐tracking DCAT delivery.

Evaluation of Dose Distribution in Intensity Modulated Radiosurgery for Lung Cancer under Condition of Respiratory Motion

The dose of a real tumor target volume and surrounding organs at risk (OARs) under the effect of respiratory motion was calculated for a lung tumor plan, based on the target volume covering the whole tumor motion range for intensity modulated radiosurgery (IMRS). Two types of IMRS plans based on simulated respiratory motion were designed using humanoid and dynamic phantoms. Delivery quality assurance (DQA) was performed using ArcCHECK and MapCHECK2 for several moving conditions of the tumor and the real dose inside the humanoid phantom was evaluated using the 3DVH program. This evaluated dose in the tumor target and OAR using the 3DVH program was higher than the calculated dose in the plan, and a greater difference was seen for the RapidArc treatment than for the standard intensity modulated radiation therapy (IMRT) with fixed gantry angle beams. The results of this study show that for IMRS plans based on target volume, including the whole tumor motion range, tighter constraints of the OAR should be considered in the optimization process. The method devised in this study can be applied effectively to analyze the dose distribution in the real volume of tumor target and OARs in IMRT plans targeting the whole tumor motion range.

3D delivered dose assessment using a 4DCT-based motion model

PURPOSE

The purpose of this work is to develop a clinically feasible method of calculating actual delivered dose distributions for patients who have significant respiratory motion during the course of stereotactic body radiation therapy (SBRT).

METHODS

A novel approach was proposed to calculate the actual delivered dose distribution for SBRT lung treatment. This approach can be specified in three steps. (1) At the treatment planning stage, a patient-specific motion model is created from planning 4DCT data. This model assumes that the displacement vector field (DVF) of any respiratory motion deformation can be described as a linear combination of some basis DVFs. (2) During the treatment procedure, 2D time-varying projection images (either kV or MV projections) are acquired, from which time-varying "fluoroscopic" 3D images of the patient are reconstructed using the motion model. The DVF of each timepoint in the time-varying reconstruction is an optimized linear combination of basis DVFs such that the 2D projection of the 3D volume at this timepoint matches the projection image. (3) 3D dose distribution is computed for each timepoint in the set of 3D reconstructed fluoroscopic images, from which the total effective 3D delivered dose is calculated by accumulating deformed dose distributions. This approach was first validated using two modified digital extended cardio-torso (XCAT) phantoms with lung tumors and different respiratory motions. The estimated doses were compared to the dose that would be calculated for routine 4DCT-based planning and to the actual delivered dose that was calculated using "ground truth" XCAT phantoms at all timepoints. The approach was also tested using one set of patient data, which demonstrated the application of our method in a clinical scenario.

RESULTS

For the first XCAT phantom that has a mostly regular breathing pattern, the errors in 95% volume dose (D95) are 0.11% and 0.83%, respectively for 3D fluoroscopic images reconstructed from kV and MV projections compared to the ground truth, which is clinically comparable to 4DCT (0.093%). For the second XCAT phantom that has an irregular breathing pattern, the errors are 0.81% and 1.75% for kV and MV reconstructions, both of which are better than that of 4DCT (4.01%). In the case of real patient, although it is impossible to obtain the actual delivered dose, the dose estimation is clinically reasonable and demonstrates differences between 4DCT and MV reconstruction-based dose estimates.

CONCLUSIONS

With the availability of kV or MV projection images, the proposed approach is able to assess delivered doses for all respiratory phases during treatment. Compared to the planning dose based on 4DCT, the dose estimation using reconstructed 3D fluoroscopic images was as good as 4DCT for regular respiratory pattern and was a better dose estimation for the irregular respiratory pattern.

Robustness of sweeping-window arc therapy treatment sequences against intrafractional tumor motion

PURPOSE

Due to the potentially periodic collimator dynamic in volumetric modulated arc therapy (VMAT) dose deliveries with the sweeping-window arc therapy (SWAT) technique, additional manifestations of dosimetric deviations in the presence of intrafractional motion may occur. With a fast multileaf collimator (MLC), and a flattening filter free dose delivery, treatment times close to 60 s per fraction are clinical reality. For these treatment sequences, the human breathing period can be close to the collimator sweeping period. Compared to a random arrangement of the segments, this will cause a further degradation of the dose homogeneity.

METHODS

Fifty VMAT sequences of potentially moving target volumes were delivered on a two dimensional ionization chamber array. In order to detect interplay effects along all three coordinate axes, time resolved measurements were performed twice--with the detector aligned in vertical (V) or horizontal (H) orientation. All dose matrices were then moved within a simulation software by a time-dependent motion vector. The minimum relative equivalent uniform dose EUDr,m for all breathing starting phases was determined for each amplitude and period. Furthermore, an estimation of periods with minimum EUD was performed. Additionally, LINAC logfiles were recorded during plan delivery. The MLC, jaw, gantry angle, and monitor unit settings were continuously saved and used to calculate the correlation coefficient between the target motion and the dose weighed collimator motion component for each direction (CC, LR, AP) separately.

RESULTS

The resulting EUDr,m were EUDr,m(CCV) = (98.3 ± 0.6)%, EUDr,m(CCH) = (98.6 ± 0.5)%, EUDr,m(APV) = (97.7 ± 0.9)%, and EUDr,m(LRH) = (97.8 ± 0.9)%. The overall minimum relative EUD observed for 360(∘) arc midventilation treatments was 94.6%. The treatment plan with the shortest period and a minimum relative EUD of less than 97% was found at T = 6.1 s. For a partial 120(∘) arc, an EUDr,m = 92.0% was found. In all cases, a correlation coefficient above 0.5 corresponded to a minimum in EUD.

CONCLUSIONS

With the advent of fast VMAT delivery techniques, nonrobust treatment sequences for human breathing patterns can be generated. These sequences are characterized by a large correlation coefficient between a target motion component and the corresponding collimator dynamic. By iteratively decreasing the maximum allowed dose rate, a low correlation coefficient and consequentially a robust treatment sequence are ensured.

Markerless tumor tracking using short kilovoltage imaging arcs for lung image-guided radiotherapy

The ability to monitor tumor motion without implanted markers is clinically advantageous for lung image-guided radiotherapy (IGRT). Existing markerless tracking methods often suffer from overlapping structures and low visibility of tumors on kV projection images. We introduce the short arc tumor tracking (SATT) method to overcome these issues. The proposed method utilizes multiple kV projection images selected from a nine-degree imaging arc to improve tumor localization, and respiratory-correlated 4D cone-beam CT (CBCT) prior knowledge to minimize the effects of overlapping anatomies. The 3D tumor position is solved as an optimization problem with prior knowledge incorporated via regularization. We retrospectively validated SATT on 11 clinical scans from four patients with central tumors. These patients represent challenging scenarios for markerless tumor tracking due to the inferior adjacent contrast. The 3D trajectories of implanted fiducial markers were used as the ground truth for tracking accuracy evaluation. In all cases, the tumors were successfully tracked at all gantry angles. Compared to standard pre-treatment CBCT guidance alone, trajectory errors were significantly smaller with tracking in all cases, and the improvements were the most prominent in the superior-inferior direction. The mean 3D tracking error ranged from 2.2-9.9 mm, which was 0.4-2.6 mm smaller compared to pre-treatment CBCT. In conclusion, we were able to directly track tumors with inferior visibility on kV projection images using SATT. Tumor localization accuracies are significantly better with tracking compared to the current standard of care of lung IGRT. Future work involves the prospective evaluation and clinical implementation of SATT.

The impact of low-Z and high-Z metal implants in IMRT: a Monte Carlo study of dose inaccuracies in commercial dose algorithms

PURPOSE

The aim of the study was to evaluate the dosimetric impact of low-Z and high-Z metallic implants on IMRT plans.

METHODS

Computed tomography (CT) scans of three patients were analyzed to study effects due to the presence of Titanium (low-Z), Platinum and Gold (high-Z) inserts. To eliminate artifacts in CT images, a sinogram-based metal artifact reduction algorithm was applied. IMRT dose calculations were performed on both the uncorrected and corrected images using a commercial planning system (convolution/superposition algorithm) and an in-house Monte Carlo platform. Dose differences between uncorrected and corrected datasets were computed and analyzed using gamma index (Pγ<1) and setting 2 mm and 2% as distance to agreement and dose difference criteria, respectively. Beam specific depth dose profiles across the metal were also examined.

RESULTS

Dose discrepancies between corrected and uncorrected datasets were not significant for low-Z material. High-Z materials caused under-dosage of 20%-25% in the region surrounding the metal and over dosage of 10%-15% downstream of the hardware. Gamma index test yielded Pγ<1>99% for all low-Z cases; while for high-Z cases it returned 91% < Pγ<1< 99%. Analysis of the depth dose curve of a single beam for low-Z cases revealed that, although the dose attenuation is altered inside the metal, it does not differ downstream of the insert. However, for high-Z metal implants the dose is increased up to 10%-12% around the insert. In addition, Monte Carlo method was more sensitive to the presence of metal inserts than superposition/convolution algorithm.

CONCLUSIONS

The reduction in terms of dose of metal artifacts in CT images is relevant for high-Z implants. In this case, dose distribution should be calculated using Monte Carlo algorithms, given their superior accuracy in dose modeling in and around the metal. In addition, the knowledge of the composition of metal inserts improves the accuracy of the Monte Carlo dose calculation significantly.

A mass-conserving 4D XCAT phantom for dose calculation and accumulation

PURPOSE

The XCAT phantom is a realistic 4D digital torso phantom that is widely used in imaging and therapy research. However, lung mass is not conserved between respiratory phases of the phantom, making detailed dosimetric simulations and dose accumulation unphysical. A framework is developed to correct this issue by enforcing local mass conservation in the XCAT lung. Dose calculations are performed to assess the implications of neglecting mass conservation, and to demonstrate an application of the phantom to calculate the accumulated delivered dose in an irregularly breathing patient.

METHODS

A displacement vector field (DVF) between each respiratory state and a reference image is generated from the XCAT motion model and its divergence is calculated and used to correct the lung density. A series of phantoms with regular and irregular breathing (based on patient data) are generated and modified to conserve mass. Monte Carlo methods are used to simulate conventional and SBRT treatment delivery. The calculated dose is deformed and accumulated using the DVF. Results from the mass-conserving and original XCAT are compared. A 4DCT is simulated for the irregularly breathing patient, and a 4DCT-based dose estimate is compared with the accumulated delivered dose.

RESULTS

The presented framework successfully conserves mass in the XCAT lung. The spatial distribution of the lung dose was qualitatively changed by the use of a mass conservation in the XCAT; however, the corresponding DVH did not change significantly. The comparison of the delivered dose with the 4DCT-based prediction shows similar lung metric results, however dose differences of 10% can be seen in some spatial regions.

CONCLUSIONS

The XCAT phantom has been successfully modified so that it conserves lung mass during respiration, enabling it to be used as a tool to perform dose accumulation studies in the lung without relying on deformable image registration. Neglecting mass conservation can result in erroneous spatial distributions of the dose in the lung. Using this tool to simulate patient treatments reveals differences between the planned dose and the calculated delivered dose for the full treatment. The software is freely available from the authors.

Experimental evaluations of the accuracy of 3D and 4D planning in robotic tracking stereotactic body radiotherapy for lung cancers

PURPOSE

Due to the complexity of 4D target tracking radiotherapy, the accuracy of this treatment strategy should be experimentally validated against established standard 3D technique. This work compared the accuracy of 3D and 4D dose calculations in respiration tracking stereotactic body radiotherapy (SBRT).

METHODS

Using the 4D planning module of the CyberKnife treatment planning system, treatment plans for a moving target and a static off-target cord structure were created on different four-dimensional computed tomography (4D-CT) datasets of a thorax phantom moving in different ranges. The 4D planning system used B-splines deformable image registrations (DIR) to accumulate dose distributions calculated on different breathing geometries, each corresponding to a static 3D-CT image of the 4D-CT dataset, onto a reference image to compose a 4D dose distribution. For each motion, 4D optimization was performed to generate a 4D treatment plan of the moving target. For comparison with standard 3D planning, each 4D plan was copied to the reference end-exhale images and a standard 3D dose calculation was followed. Treatment plans of the off-target structure were first obtained by standard 3D optimization on the end-exhale images. Subsequently, they were applied to recalculate the 4D dose distributions using DIRs. All dose distributions that were initially obtained using the ray-tracing algorithm with equivalent path-length heterogeneity correction (3D EPL and 4D EPL) were recalculated by a Monte Carlo algorithm (3D MC and 4D MC) to further investigate the effects of dose calculation algorithms. The calculated 3D EPL, 3D MC, 4D EPL, and 4D MC dose distributions were compared to measurements by Gafchromic EBT2 films in the axial and coronal planes of the moving target object, and the coronal plane for the static off-target object based on the γ metric at 5%/3mm criteria (γ5%/3mm). Treatment plans were considered acceptable if the percentage of pixels passing γ5%/3mm (Pγ<1) ≥ 90%.

RESULTS

The averaged Pγ<1 values of the 3D EPL, 3D MC, 4D EPL, and 4D MC dose calculation methods for the moving target plans are 95%, 95%, 94%, and 95% for reproducible motion, and 95%, 96%, 94%, and 93% for nonreproducible motion during actual treatment delivery. The overall measured target dose distributions are in better agreement with the 3DMC dose distributions than the 4DMC dose distributions. Conversely, measured dose distributions agree much better with the 4D EPL/MC than the 3D EPL/MC dose distributions in the static off-target structure, resulting in higher Pγ<1 values with 4D EPL/MC (91%) vs 3D EPL (24%) and 3D MC (25%). Systematic changes of target motion reduced the averaged Pγ<1 to 47% and 53% for 4D EPL and 4D MC dose calculations, and 22% for 3D EPL/MC dose calculations in the off-target films.

CONCLUSIONS

In robotic tracking SBRT, 4D treatment planning was found to yield better prediction of the dose distributions in the off-target structure, but not necessarily in the moving target, compared to standard 3D treatment planning, for reproducible and nonreproducible target motion. It is important to ensure on a patient-by-patient basis that the cumulative uncertainty associated with the 4D-CT artifacts, deformable image registration, and motion variability is significantly smaller than the cumulative uncertainty occurred in standard 3D planning in order to make 4D planning a justified option.

Visualisation of respiratory tumour motion and co-moving isodose lines in the context of respiratory gating, IMRT and flattening-filter-free beams

Respiratory motion during percutaneous radiotherapy can be considered based on respiration-correlated computed tomography (4DCT). However, most treatment planning systems perform the dose calculation based on a single primary CT data set, even though cine mode displays may allow for a visualisation of the complete breathing cycle. This might create the mistaken impression that the dose distribution were independent of tumour motion. We present a movie visualisation technique with the aim to direct attention to the fact that the dose distribution migrates to some degree with the tumour and discuss consequences for gated treatment, IMRT plans and flattening-filter-free beams. This is a feasibility test for a visualisation of tumour and isodose motion. Ten respiratory phases are distinguished on the CT, and the dose distribution from a stationary IMRT plan is calculated on each phase, to be integrated into a movie of tumour and dose motion during breathing. For one example patient out of the sample of five lesions, the plan is compared with a gated treatment plan with respect to tumour coverage and lung sparing. The interplay-effect for small segments in the IMRT plan is estimated. While the high dose rate, together with the cone-shaped beam profile, makes the use of flattening-filter-free beams more problematic for conformal and IMRT treatment, it can be the option of choice if gated treatment is preferred. The different effects of respiratory motion, dose build-up and beam properties (segments and flatness) for gated vs. un-gated treatment can best be considered if planning is performed on the full 4DCT data set, which may be an incentive for future developments of treatment planning systems.

Dosimetric assessment of a novel metal artifact reduction method in CT images

The aim of this study was to assess the ability of metal artifact reduction (MAR) algorithm in restoring the CT image quality while correcting the tissue density information for the accurate estimation of the absorbed dose. A phantom filled with titanium (low-Z metal) and Cerrobend (high-Z metal) inserts was used for this purpose. The MAR algorithm was applied to phantom's CT dataset. Static intensity-modulated radiation therapy (IMRT) plans, including five beam angles, were designed and optimized on the uncorrected images to deliver 10 Gy on the simulated target. Monte Carlo dose calculation was computed on uncorrected, corrected, and ground truth image datasets. It was firstly verified that MAR methodology was able to correct HU errors due to the metal presence. In the worst situation (high-Z phantom), the image difference, uncorrected ground truth and corrected ground truth, went from -4.4 ± 118.8 HU to 0.4 ± 10.8 HU, respectively. Secondly, it was observed that the impact of dose errors estimation depends on the atomic number of the metal: low-Z inserts do not produce significant dose inaccuracies, while high-Z implants substantially influence the computation of the absorbed dose. In this latter case, dose errors in the PTV region were up to 23.56% (9.72% mean value) when comparing the uncorrected vs. the ground truth dataset. After MAR correction, errors dropped to 0.11% (0.10% mean value). In conclusion, it was assessed that the new MAR algorithm is able to restore image quality without distorting mass density information, thus producing a more accurate dose estimation.

Numerical solutions of the γ-index in two and three dimensions

The γ-index is used routinely to establish correspondence between two dose distributions. The definition of the γ-index can be written with a single equation but solving this equation at millions of points is computationally expensive, especially in three dimensions. Our goal is to extend the vector-equation method in Bakai et al (2003 Phys. Med. Biol.48 3543-53) to higher order for better accuracy and, as important, to determine the magnitude of accuracy in a higher order solution. We construct a numerical framework for calculating the γ-index in two and three dimensions and present an efficient method for calculating the γ-index with zeroth-, first- and second-order methods using tricubic spline interpolation. For an intensity-modulated radiation therapy example with 1.78 × 10⁶ voxels, the zeroth-order, first-order, first-order iterations and semi-second-order methods calculate the three-dimensional γ-index in 1.5, 4.7, 34.7 and 35.6 s with 36.7%, 1.1%, 0.2% and 0.8% accuracy, respectively. The accuracy of linear interpolation with this example is 1.0%. We present efficient numerical methods for calculating the three-dimensional γ-index with tricubic spline interpolation. The first-order method with iterations is the most accurate and fastest choice of the numerical methods if the dose distributions may have large second-order gradients. Furthermore, the difference between iterations can be used to determine the accuracy of the method.

Cardiac-sparing whole lung IMRT in children with lung metastasis

PURPOSE

To demonstrate the dosimetric advantages of cardiac-sparing (CS) intensity modulated radiation therapy (IMRT) in children undergoing whole lung irradiation (WLI).

METHODS AND MATERIALS

Chest CT scans of 22 children who underwent simulation with 3-dimensional (n=10) or 4-dimensional (n=12) techniques were used for this study. Treatment planning was performed using standard anteroposterior-posteroanterior (S-RT) technique and CS-IMRT. Left and right flank fields were added to WLI fields to determine whether CS-IMRT offered any added protection to normal tissues at the junction between these fields. The radiation dose to the lung PTV, cardiac structures, liver, and thyroid were analyzed and compared.

RESULTS

CS-IMRT had 4 significant advantages over S-RT: (1) superior cardiac protection (2) superior 4-dimensional lung planning target volume coverage, (3) superior dose uniformity in the lungs with fewer hot spots, and (4) significantly lower dose to the heart when flank RT is administered after WLI.

CONCLUSIONS

The use of CS-IMRT and 4-dimensional treatment planning has the potential to improve tumor control rates and reduce cardiac toxicity in children receiving WLI.

Quality assurance for nonradiographic radiotherapy localization and positioning systems: report of Task Group 147

New technologies continue to be developed to improve the practice of radiation therapy. As several of these technologies have been implemented clinically, the Therapy Committee and the Quality Assurance and Outcomes Improvement Subcommittee of the American Association of Physicists in Medicine commissioned Task Group 147 to review the current nonradiographic technologies used for localization and tracking in radiotherapy. The specific charge of this task group was to make recommendations about the use of nonradiographic methods of localization, specifically; radiofrequency, infrared, laser, and video based patient localization and monitoring systems. The charge of this task group was to review the current use of these technologies and to write quality assurance guidelines for the use of these technologies in the clinical setting. Recommendations include testing of equipment for initial installation as well as ongoing quality assurance. As the equipment included in this task group continues to evolve, both in the type and sophistication of technology and in level of integration with treatment devices, some of the details of how one would conduct such testing will also continue to evolve. This task group, therefore, is focused on providing recommendations on the use of this equipment rather than on the equipment itself, and should be adaptable to each user's situation in helping develop a comprehensive quality assurance program.

Respiration-correlated image guidance is the most important radiotherapy motion management strategy for most lung cancer patients

PURPOSE

The purpose of this study was to quantify the effects of four-dimensional computed tomography (4DCT), 4D image guidance (4D-IG), and beam gating on calculated treatment field margins in a lung cancer patient population.

MATERIALS AND METHODS

Images were acquired from 46 lung cancer patients participating in four separate protocols at three institutions in Europe and the United States. Seven patients were imaged using fluoroscopy, and 39 patients were imaged using 4DCT. The magnitude of respiratory tumor motion was measured. The required treatment field margins were calculated using a statistical recipe (van Herk M, et al. Int J Radiat Oncol Biol Phys 2000;474:1121-1135), with magnitudes of all uncertainties, except respiratory peak-to-peak displacement, the same for all patients, taken from literature. Required margins for respiratory motion management were calculated using the residual respiratory tumor motion for each patient for various motion management strategies. Margin reductions for respiration management were calculated using 4DCT, 4D-IG, and gated beam delivery.

RESULTS

The median tumor motion magnitude was 4.4 mm for the 46 patients (range 0-29.3 mm). This value corresponded to required treatment field margins of 13.7 to 36.3 mm (median 14.4 mm). The use of 4DCT, 4D-IG, and beam gating required margins that were reduced by 0 to 13.9 mm (median 0.5 mm), 3 to 5.2 mm (median 5.1 mm), and 0 to 7 mm (median 0.2 mm), respectively, to a total of 8.5 to 12.4 mm (median 8.6 mm).

CONCLUSION

A respiratory management strategy for lung cancer radiotherapy including planning on 4DCT scans and daily image guidance provides a potential reduction of 37% to 47% in treatment field margins. The 4D image guidance strategy was the most effective strategy for >85% of the patients.

Towards accurate dose accumulation for Step-&-Shoot IMRT: Impact of weighting schemes and temporal image resolution on the estimation of dosimetric motion effects

PURPOSE

Breathing-induced motion effects on dose distributions in radiotherapy can be analyzed using 4D CT image sequences and registration-based dose accumulation techniques. Often simplifying assumptions are made during accumulation. In this paper, we study the dosimetric impact of two aspects which may be especially critical for IMRT treatment: the weighting scheme for the dose contributions of IMRT segments at different breathing phases and the temporal resolution of 4D CT images applied for dose accumulation.

METHODS

Based on a continuous problem formulation a patient- and plan-specific scheme for weighting segment dose contributions at different breathing phases is derived for use in step-&-shoot IMRT dose accumulation. Using 4D CT data sets and treatment plans for 5 lung tumor patients, dosimetric motion effects as estimated by the derived scheme are compared to effects resulting from a common equal weighting approach. Effects of reducing the temporal image resolution are evaluated for the same patients and both weighting schemes.

RESULTS

The equal weighting approach underestimates dosimetric motion effects when considering single treatment fractions. Especially interplay effects (relative misplacement of segments due to respiratory tumor motion) for IMRT segments with only a few monitor units are insufficiently represented (local point differences >25% of the prescribed dose for larger tumor motion). The effects, however, tend to be averaged out over the entire treatment course. Regarding temporal image resolution, estimated motion effects in terms of measures of the CTV dose coverage are barely affected (in comparison to the full resolution) when using only half of the original resolution and equal weighting. In contrast, occurence and impact of interplay effects are poorly captured for some cases (large tumor motion, undersized PTV margin) for a resolution of 10/14 phases and the more accurate patient- and plan-specific dose accumulation scheme.

CONCLUSIONS

Radiobiological consequences of reported single fraction local point differences >25% of the prescribed dose are widely unclear and should be subject to future investigation. Meanwhile, if aiming at accurate and reliable estimation of dosimetric motion effects, precise weighting schemes such as the presented patient- and plan-specific scheme for step-&-shoot IMRT and full available temporal 4D CT image resolution should be applied for IMRT dose accumulation.

Comparison of helical and average computed tomography for stereotactic body radiation treatment planning and normal tissue contouring in lung cancer

AIMS

To compare average computed tomography (CT(AVE)) datasets with free breathing helical computed tomography (CT(HEL)) for contouring organs at risk (OARs) and radiation treatment planning in patients receiving stereotactic body radiation therapy (SBRT) in the lung.

MATERIALS AND METHODS

Ten SBRT patients with peripheral stage I non-small cell lung cancer underwent a CT(HEL) and a four-dimensional computed tomography scan in the treatment position. CT(AVE) datasets were generated from the four-dimensional computed tomography scan. The following OARs were delineated on the CT(HEL) and CT(AVE) datasets of each patient: lung minus internal target volume, trachea/main bronchus, heart, oesophagus and spinal cord. Volumes and geometric centres of the OARs, as well as the dosimetric impact of planning with these different datasets, were examined.

RESULTS

There were no statistical differences in the OAR geometric centre coordinates nor in the OAR volumes between the CT(HEL) and CT(AVE) datasets, except that CT(AVE)-defined trachea/main bronchus and lung minus internal target volume mean volumes were larger than those defined on the CT(HEL) (46, 43 cm³ and 3516, 3378 cm³, respectively, P<0.05). Despite this, there were no significant differences in the mean and maximum doses to the OAR contours when using the CT(HEL) or CT(AVE) for planning (<4% average change in the maximum and mean doses for all OARs, P>0.05). There were also no significant differences in the locations of the hotspots within OARs among the CT(HEL) or CT(AVE) datasets (P>0.05).

CONCLUSIONS

CT(AVE) datasets may be used in place of CT(HEL) for OAR contouring and dose calculations. When four-dimensional computed tomography is available, it may not be necessary to acquire a separate CT(HEL) scan for OAR contouring or dosimetric purposes.

Motion management with phase-adapted 4D-optimization

Cancer treatment with ionizing radiation is often compromised by organ motion, in particular for lung cases. Motion uncertainties can significantly degrade an otherwise optimized treatment plan. We present a spatiotemporal optimization method, which takes into account all phases of breathing via the corresponding 4D-CTs and provides a 4D-optimal plan that can be delivered throughout all breathing phases. Monte Carlo dose calculations are employed to warrant for highest dosimetric accuracy, as pertinent to study motion effects in lung. We demonstrate the performance of this optimization method with clinical lung cancer cases and compare the outcomes to conventional gating techniques. We report significant improvements in target coverage and in healthy tissue sparing at a comparable computational expense. Furthermore, we show that the phase-adapted 4D-optimized plans are robust against irregular breathing, as opposed to gating. This technique has the potential to yield a higher delivery efficiency and a decisively shorter delivery time.

Extension of the NCAT phantom for the investigation of intra-fraction respiratory motion in IMRT using 4D Monte Carlo

The purpose of this work was to create a computational platform for studying motion in intensity modulated radiotherapy (IMRT). Specifically, the non-uniform rational B-spline (NURB) cardiac and torso (NCAT) phantom was modified for use in a four-dimensional Monte Carlo (4D-MC) simulation system to investigate the effect of respiratory-induced intra-fraction organ motion on IMRT dose distributions as a function of diaphragm motion, lesion size and lung density. Treatment plans for four clinical scenarios were designed: diaphragm peak-to-peak amplitude of 1 cm and 3 cm, and two lesion sizes--2 cm and 4 cm diameter placed in the lower lobe of the right lung. Lung density was changed for each phase using a conservation of mass calculation. Further, a new heterogeneous lung model was implemented and tested. Each lesion had an internal target volume (ITV) subsequently expanded by 15 mm isotropically to give the planning target volume (PTV). The PTV was prescribed to receive 72 Gy in 40 fractions. The MLC leaf sequence defined by the planning system for each patient was exported and used as input into the MC system. MC simulations using the dose planning method (DPM) code together with deformable image registration based on the NCAT deformation field were used to find a composite dose distribution for each phantom. These composite distributions were subsequently analyzed using information from the dose volume histograms (DVH). Lesion motion amplitude has the largest effect on the dose distribution. Tumor size was found to have a smaller effect and can be mitigated by ensuring the planning constraints are optimized for the tumor size. The use of a dynamic or heterogeneous lung density model over a respiratory cycle does not appear to be an important factor with a <or=0.6% change in the mean dose received by the ITV, PTV and right lung. The heterogeneous model increases the realism of the NCAT phantom and may provide more accurate simulations in radiation therapy investigations that use the phantom. This work further evaluates the NCAT phantom for use as a tool in radiation therapy research in addition to its extensive use in diagnostic imaging and nuclear medicine research. Our results indicate that the NCAT phantom, combined with 4D-MC simulations, is a useful tool in radiation therapy investigations and may allow the study of relative effects in many clinically relevant situations.

Monte Carlo dose mapping on deforming anatomy

This paper proposes a Monte Carlo-based energy and mass congruent mapping (EMCM) method to calculate the dose on deforming anatomy. Different from dose interpolation methods, EMCM separately maps each voxel's deposited energy and mass from a source image to a reference image with a displacement vector field (DVF) generated by deformable image registration (DIR). EMCM was compared with other dose mapping methods: energy-based dose interpolation (EBDI) and trilinear dose interpolation (TDI). These methods were implemented in EGSnrc/DOSXYZnrc, validated using a numerical deformable phantom and compared for clinical CT images. On the numerical phantom with an analytically invertible deformation map, EMCM mapped the dose exactly the same as its analytic solution, while EBDI and TDI had average dose errors of 2.5% and 6.0%. For a lung patient's IMRT treatment plan, EBDI and TDI differed from EMCM by 1.96% and 7.3% in the lung patient's entire dose region, respectively. As a 4D Monte Carlo dose calculation technique, EMCM is accurate and its speed is comparable to 3D Monte Carlo simulation. This method may serve as a valuable tool for accurate dose accumulation as well as for 4D dosimetry QA.

A radiobiological analysis of the effect of 3D versus 4D image-based planning in lung cancer radiotherapy

Dose distributions generated on a static anatomy may differ significantly from those delivered to temporally varying anatomy such as for abdominal and thoracic tumors, due largely in part to the unavoidable organ motion and deformation effects stemming from respiration. In this work, the degree of such variation for three treatment techniques, namely static conventional, gating and target tracking radiotherapy, was investigated. The actual delivered dose was approximated by planning all the phases of a 4DCT image set. Data from six (n = 6) previously treated lung cancer patients were used for this study with tumor motion ranging from 2 to 10 mm. Complete radiobiological analyses were performed to assess the clinical significance of the observed discrepancies between the 3D and 4DCT image-based dose distributions. Using the complication-free tumor control probability (P+) objective, we observed small differences in P+ between the 3D and 4DCT image-based plans (<2.0% difference on average) for the gating and static conventional regimens and higher differences in P+ (4.0% on average) for the tracking regimen. Furthermore, we observed, as a general trend, that the 3D plan underestimated the P+ values. While it is not possible to draw any general conclusions from a small patient cohort, our results suggest that there exists a patient population in which 4D planning does not provide any additional benefits beyond that afforded by 3D planning for static conventional or gated radiotherapy. This statement is consistent with previous studies based on physical dosimetric evaluations only. The higher differences observed with the tracking technique suggest that individual patient plans should be evaluated on a case-by-case basis to assess if 3D or 4D imaging is appropriate for the tracking technique.

Cumulative lung dose for several motion management strategies as a function of pretreatment patient parameters

PURPOSE

To evaluate patient parameters that may predict for relative differences in cumulative four-dimensional (4D) lung dose among several motion management strategies.

METHODS AND MATERIALS

Deformable image registration and dose accumulation were used to generate 4D treatment plans for 18 patients with 4D computed tomography scans. Three plans were generated to simulate breath hold at normal inspiration, target tracking with the beam aperture, and mid-ventilation aperture (control of the target at the mean daily position and application of an iteratively computed margin to compensate for respiration). The relative reduction in mean lung dose (MLD) between breath hold and mid-ventilation aperture (DeltaMLD(BH)) and between target tracking and mid-ventilation aperture (DeltaMLD(TT)) was calculated. Associations between these two variables and parameters of the lesion (excursion, size, location, and deformation) and dose distribution (local dose gradient near the target) were also calculated.

RESULTS

The largest absolute and percentage differences in MLD were 1.0 Gy and 21.5% between breath hold and mid-ventilation aperture. DeltaMLD(BH) was significantly associated (p < 0.05) with tumor excursion. The DeltaMLD(TT) was significantly associated with excursion, deformation, and local dose gradient. A linear model was constructed to represent DeltaMLD vs. excursion. For each 5 mm of excursion, target tracking reduced the MLD by 4% compared with the results of a mid-ventilation aperture plan. For breath hold, the reduction was 5% per 5 mm of excursion.

CONCLUSIONS

The relative difference in MLD among different motion management strategies varied with patient and tumor characteristics for a given dosimetric target coverage. Tumor excursion is useful to aid in stratifying patients according to appropriate motion management strategies.

Design and testing of a simulation framework for dosimetric motion studies integrating an anthropomorphic computational phantom into four-dimensional Monte Carlo

We have designed a simulation framework for motion studies in radiation therapy by integrating the anthropomorphic NCAT phantom into a 4D Monte Carlo dose calculation engine based on DPM. Representing an artifact-free environment, the system can be used to identify class solutions as a function of geometric and dosimetric parameters. A pilot dynamic conformal study for three lesions ( approximately 2.0 cm) in the right lung was performed (70 Gy prescription dose). Tumor motion changed as a function of tumor location, according to the anthropomorphic deformable motion model. Conformal plans were simulated with 0 to 2 cm margin for the aperture, with additional 0.5 cm for beam penumbra. The dosimetric effects of intensity modulated radiotherapy (IMRT) vs. conformal treatments were compared in a static case. Results show that the Monte Carlo simulation framework can model tumor tracking in deformable anatomy with high accuracy, providing absolute doses for IMRT and conformal radiation therapy. A target underdosage of up to 3.67 Gy (lower lung) was highlighted in the composite dose distribution mapped at exhale. Such effects depend on tumor location and treatment margin and are affected by lung deformation and ribcage motion. In summary, the complexity in the irradiation of moving targets has been reduced to a controlled simulation environment, where several treatment options can be accurately modeled and quantified The implemented tools will be utilized for extensive motion study in lung/liver irradiation.

Experimental measurements and Monte Carlo simulations for dosimetric evaluations of intrafraction motion for gated and ungated intensity modulated arc therapy deliveries

Respiratory gated radiation therapy allows for a smaller margin expansion for the planning target volume (PTV) to account for respiratory induced motion and is emerging as a common method to treat lung and liver tumors. We investigated the dosimetric effect of free motion and gated delivery for intensity modulated arc therapy (IMAT) with experimental measurements and Monte Carlo simulations. The impact of PTV margin and duty cycle for gated delivery is studied with Monte Carlo simulations. A motion phantom is used for this study. Two sets of contours were drawn on the mid-inspiration CT scan of this motion phantom. For each set of contours, an IMAT plan to be delivered with constant dose rate was created. The plans were generated on a CT scan of the phantom in the static condition with 3 mm PTV margin and applied to the motion phantom under four conditions: static, full superior-inferior (SI) motion (A = 1 cm, T = 4 s) and gating conditions (25% and 50% duty cycles) with full SI motion. A 6 by 15 cm piece of radiographic film was placed in the sagittal plane of the phantom and then irradiated under all measurement conditions. Film calibration was performed with a step-wedge method to convert optical density to dose. Gated IMAT delivery was first validated in 2D by comparing static film with that from gating and full motion. A previously verified simulation tool for IMRT that takes the log files from the multileaf collimator (MLC) controller and the gating system were adapted to simulate the delivered IMAT treatment for full 3D dosimetric analysis. The IMAT simulations were validated against the 2D film measurements. The resultant IMAT simulations were evaluated with dose criteria, dose-volume histograms and 3D gamma analysis. We validated gated IMAT deliveries when we compared the static film with the one from gating using 25% duty cycle using 2D gamma analysis. Within experimental and setup uncertainties, film measurements agreed with their corresponding simulated plans using 2D gamma analysis. Finally, when planning with margins designed for gating with 25% duty cycle and applying 50% or no gating during treatment, the dose differences in D(min,) D(99%) and D(95%) of the clinical target volume can be up to 27 cGy, 20 cGy and 18 cGy, respectively, for a plan with 200 cGy prescription dose. We have experimentally delivered gated IMAT with constant dose rate to a motion phantom and assessed their accuracies with film dosimetry and Monte Carlo simulations. Film dosimetry demonstrated that 25% gating and static plans are within 2%, 2 mm. The Monte Carlo simulation method was employed to generate dose delivered in 3D to a motion phantom, and the dosimetric results were reported. Since our film measurements agreed well with Monte Carlo simulations, we can reliably use this simulation tool to further study the dosimetric effects of target motion and effectiveness of gating for IMAT deliveries.

An overview of Monte Carlo treatment planning for radiotherapy

The implementation of Monte Carlo dose calculation algorithms in clinical radiotherapy treatment planning systems has been anticipated for many years. Despite a continuous increase of interest in Monte Carlo Treatment Planning (MCTP), its introduction into clinical practice has been delayed by the extent of calculation time required. The development of newer and faster MC codes is behind the commercialisation of the first MC-based treatment planning systems. The intended scope of this article is to provide the reader with a compact 'primer' on different approaches to MCTP with particular attention to the latest developments in the field.

Monte Carlo dose calculation of segmental IMRT delivery to a moving phantom using dynamic MLC and gating log files

Respiratory gating is emerging as a tool to limit the effect of motion for liver and lung tumors. In order to study the impact of target motion and gated intensity modulated radiation therapy (IMRT) delivery, a computer program was developed to simulate segmental IMRT delivery to a moving phantom. Two distinct plans were delivered to a rigid-motion phantom with a film insert in place under four conditions: static, sinusoidal motion, gated sinusoidal motion with a duty cycle of 25% and gated sinusoidal motion with duty cycle of 50% under motion conditions of a typical patient (A = 1 cm, T = 4 s). The MLC controller log files and gating log files were retained to perform a retrospective Monte Carlo dose calculation of the plans. Comparison of the 2D planar dose distributions between simulation and measurement demonstrated that our technique had at least 94% of the points passing gamma criteria of 3% for dose difference and 3 mm as the distance to agreement. This note demonstrates that the use of dynamic multi-leaf collimator and respiratory monitoring system log files together with a fast Monte Carlo dose calculation algorithm is an accurate and efficient way to study the dosimetric effect of motion for gated or non-gated IMRT delivery on a rigidly-moving body.