A serial 4DCT study to quantify range variations in charged particle radiotherapy of thoracic cancers.

Weekly serial 4DCT scans were acquired under free breathing conditions to assess water-equivalent path length (WEL) variations due to both intrafractional and interfractional changes in tissue thickness and density and to calculate proton dose distributions resulting from anatomical variations observed in serial 4DCT. A template of region of interests (ROIs) was defined on the anterior-posterior (AP) beam's eye view, and WEL measurements were made over these ROIs to quantify chest wall thickness variations. Interfractional proton dose distributions were calculated to assess changes in the expected dose distributions caused by range variations. Mean intrafractional chest wall WEL changes during respiration varied by: -4.1 mm (<-10.2 mm), -3.6 mm (<-7.1 mm), -3.2 mm (<-5.6 mm) and -2.5 mm (<-5.1 mm) during respiration in the ITV, upper, middle and lower lung regions, respectively. The mean interfractional chest wall WEL variation at Week 6 decreased by -4.0 mm (<-8.6 mm), -9.1 mm (<-17.9 mm), -9.4 mm (<-25.3 mm) and -4.5 mm (<-15.6 mm) in the ITV, upper, middle and lower lung regions, respectively. The variations were decomposed into anterior and posterior chest wall thickness changes. Dose overshoot beyond the target was observed when the initial boli was applied throughout the treatment course. This overshoot is due to chest wall thickness variations and target positional variations. The radiological path length can vary significantly during respiration as well as over the course of several weeks of charged particle therapy. Intrafractional/interfractional chest wall thickness changes can be a significant source of range variation in treatment of lung tumors with charged particle beams, resulting in dose distribution perturbations from the initial plan. Consideration of these range variations should be made in choosing the therapeutic charged particle beam range.

A four-dimensional computed tomography analysis of multiorgan abdominal motion.

PURPOSE: To characterize and quantify multiorgan respiration-induced motion in the abdomen in liver and pancreatic cancer patients. METHODS AND MATERIALS: Four-dimensional computed tomography scans were acquired for 18 patients treated for abdominal tumors. Contours of multiple abdominal organs were drawn by the radiation oncologist at one respiratory phase; these contours were propagated to other respiratory phases by deformable registration. Three-dimensional organ models were generated from the resulting contours at each phase. Motions of the bounding box and center of mass were extracted and analyzed for the clinical target volume and organs at risk. RESULTS: On average, the center of mass motion for liver clinical target volumes was 9.7 mm (SD 5 mm) in the superior-inferior direction, with a range of 3 to 18 mm; for pancreatic tumors, the average was 5 mm (SD 1 mm) m with a range of 3 to 7 mm. Abdominal organs move in unison, but with varying amplitudes. Gating near exhale (T40-T60) reduces the range of motion by a factor of ∼10. CONCLUSION: We have used deformable registration to calculate the trajectories of abdominal organs in four dimensions, based on center of mass and bounding box motion metrics. Our results are compared with previously reported studies. Possible reasons for differences are discussed.

Proton radiography and fluoroscopy of lung tumors: a Monte Carlo study using patient-specific 4DCT phantoms.

PURPOSE: Monte Carlo methods are used to simulate and optimize a time-resolved proton range telescope (TRRT) in localization of intrafractional and interfractional motions of lung tumor and in quantification of proton range variations. METHODS: The Monte Carlo N-Particle eXtended (MCNPX) code with a particle tracking feature was employed to evaluate the TRRT performance, especially in visualizing and quantifying proton range variations during respiration. Protons of 230 MeV were tracked one by one as they pass through position detectors, patient 4DCT phantom, and finally scintillator detectors that measured residual ranges. The energy response of the scintillator telescope was investigated. Mass density and elemental composition of tissues were defined for 4DCT data. RESULTS: Proton water equivalent length (WEL) was deduced by a reconstruction algorithm that incorporates linear proton track and lateral spatial discrimination to improve the image quality. 4DCT data for three patients were used to visualize and measure tumor motion and WEL variations. The tumor trajectories extracted from the WEL map were found to be within 1 mm agreement with direct 4DCT measurement. Quantitative WEL variation studies showed that the proton radiograph is a good representation of WEL changes from entrance to distal of the target. CONCLUSIONS: MCNPX simulation results showed that TRRT can accurately track the motion of the tumor and detect the WEL variations. Image quality was optimized by choosing proton energy, testing parameters of image reconstruction algorithm, and comparing to ground truth 4DCT. The future study will demonstrate the feasibility of using the time resolved proton radiography as an imaging tool for proton treatments of lung tumors.

Uncertainties in lung motion prediction relying on external surrogate: a 4DCT study in regular vs. irregular breathers.

This paper examines the uncertainty in estimating lung motion from external surrogates for lung cancer patients with regular and irregular breathing. 4DCT data sets were analyzed using a template matching algorithm to track the spatial movement of vessel bifurcations in 12 patients. The detected internal movement of features in 3D was retrospectively synchronized with the RPM surrogate signal, and the correlation index R(2) and the prediction error were computed. Patients were classified into two groups depending on the presence or not of irregularities in their breathing pattern. Peak-to-peak values of feature motion in the SI direction ranged from 0.8 mm (upper lung) to 25.3 mm (lower lung). Some patients exhibited large motion also in the latero-lateral (10.6 mm) and anterior-posterior (12.2 mm) directions. The median +/- quartile of R(2) in SI direction was 0.89 +/- 0.09. Prediction error values were up to 4.2 mm (95th percentile) with a maximum value of 4.9 mm. Statistical differences between regular and irregular breathers were found for R(2), while prediction error depended only on the range of motion. This study is relevant for image guided radiotherapy methods that rely on external surrogates to monitor motion.

Variations in tumor size and position due to irregular breathing in 4D-CT: a simulation study.

PURPOSE: To estimate the position and volume errors in 4D-CT caused by irregular breathing. METHODS: A virtual 4D-CT scanner was designed to reproduce axial mode scans with retrospective resorting. This virtual scanner creates an artificial spherical tumor based on the specifications of the user, and recreates images that might be produced by a 4D-CT scanner using a patient breathing waveform. 155 respiratory waveforms of patients were used to test the variability of 4D-CT scans. Each breathing waveform was normalized and scaled to 1, 2, and 3 cm peak-to-peak motion, and artificial tumors with 2 and 4 cm radius were simulated for each scaled waveform. The center of mass and volume of resorted 4D-CT images were calculated and compared to the expected values of center of mass and volume for the artificial tumor. Intrasubject variability was investigated by running the virtual scanner over different subintervals of each patient's breathing waveform. RESULTS: The average error in the center of mass location of an artificial tumor was less than 2 mm standard deviation for 2 cm motion. The corresponding average error in volume was less than 4%. In the worst-case scenarios, a center of mass error of 1.0 cm standard deviation and volume errors of 30%-60% at inhale were found. Systematic errors were observed in a subset of patients due to irregular breathing, and these errors were more pronounced when the tumor volume is smaller. CONCLUSIONS: Irregular breathing during 4D-CT simulation causes systematic errors in volume and center of mass measurements. These errors are small but depend on the tumor size, motion amplitude, and degree of breathing irregularity.

Evaluation and commissioning of a surface based system for respiratory sensing in 4D CT.

The purpose of this study is to assess the temporal and reconstruction accuracy of a surface imaging system, the GateCT under ideal conditions, and compare the device with a commonly used respiratory surrogate: the Varian RPM. A clinical CT scanner, run in cine mode, was used with two optical devices, GateCT and RPM, to detect respiratory motion. A radiation detector, GM-10, triggers the X-ray on/off to GateCT system, while the RPM is directly synchronized with the CT scanner through an electronic connection. Two phantoms were imaged: the first phantom translated on a rigid plate along the anterior-posterior (AP) direction, and was used to assess the temporal synchronization of each optical system with the CT scanner. The second phantom, consisting of five spheres translating 3 cm peak-to-peak in the superior-inferior direction, was used to assess the quality of rebinned images created by GateCT and RPM. Calibration assessment showed a nearly perfect synchronization with the scanner for both the RPM and GateCT systems, thus demonstrating the good performance of the radiation detector. Results for the volume rebinning test showed discrepancies in volumes for the 3D reconstruction (compared to ground truth) of up to 36% for GateCT and up to 40% for RPM. No statistical difference was proven between the two systems in volume sorting. Errors are mainly due to phase detection inaccuracies and to the large motion of the phantom. This feasibility study assessed the consistency of two optical systems in synchronizing the respiratory signal with the image acquisition. A new patient protocol based on both RPM and GateCT will be soon started.

Effects of interfractional anatomical changes on water-equivalent pathlength in charged-particle radiotherapy of lung cancer.

Intrafractional motion and interfractional changes affect the accuracy of the delivered dose in radiotherapy, particularly in charged-particle radiotherapy. Most recent studies are focused on intrafractional motion (respiratory motion). Here, we report a quantitative simulation analysis of the effects of interfractional changes on water-equivalent pathlength (WEL) in charged-particle lung therapy. Serial four-dimensional (4D) CT scans were performed under free breathing conditions; the time span between the first and second 4DCT scans was five weeks. We quantified WEL changes between the first and second CT scans due to interfractional changes (tumor shrinkage and tissue density changes) and compared the particle-beam-stopping point between the serial 4DCT scans with use of the same initial bolus. Both tumor-shrinkage and lung-density changes were observed in a single patient over the course of therapy. The lung density decreased by approximately 0.1 g/cm(3) between the first and second-CT scans, resulting in a 1.5 cm WEL changes. Tumor shrinkage resulted in approximately 3 cm WEL changes. If the same initial bolus and plan were used through the treatment course, an unexpected significant beam overshoot would occur by interfractional changes due to tumor shrinkage and lung density variation.

Accuracy in breast shape alignment with 3D surface fitting algorithms.

Surface imaging is in use in radiotherapy clinical practice for patient setup optimization and monitoring. Breast alignment is accomplished by searching for a tentative spatial correspondence between the reference and daily surface shape models. In this study, the authors quantify whole breast shape alignment by relying on texture features digitized on 3D surface models. Texture feature localization was validated through repeated measurements in a silicone breast phantom, mounted on a high precision mechanical stage. Clinical investigations on breast shape alignment included 133 fractions in 18 patients treated with accelerated partial breast irradiation. The breast shape was detected with a 3D video based surface imaging system so that breathing was compensated. An in-house algorithm for breast alignment, based on surface fitting constrained by nipple matching (constrained surface fitting), was applied. Results were compared with a commercial software where no constraints are utilized (unconstrained surface fitting). Texture feature localization was validated within 2 mm in each anatomical direction. Clinical data show that unconstrained surface fitting achieves adequate accuracy in most cases, though nipple mismatch is considerably higher than residual surface distances (3.9 mm vs 0.6 mm on average). Outliers beyond 1 cm can be experienced as the result of a degenerate surface fit, where unconstrained surface fitting is not sufficient to establish spatial correspondence. In the constrained surface fitting algorithm, average surface mismatch within 1 mm was obtained when nipple position was forced to match in the [1.5; 5] mm range. In conclusion, optimal results can be obtained by trading off the desired overall surface congruence vs matching of selected landmarks (constraint). Constrained surface fitting is put forward to represent an improvement in setup accuracy for those applications where whole breast positional reproducibility is an issue.

A review of image-guided radiotherapy.

Image-guided radiotherapy (IGRT) is in the midst of a strong development and implementation cycle, stimulated by pioneering work performed in Japan. We present a review of the rationale, technology, and methodology of image guidance, as well as an overview of current work in IGRT at the Massachusetts General Hospital. The technology is rapidly evolving, and synergisms between the various acquisition approaches are converging to provide unparalleled information on target and normal tissue location and motion. With these new approaches to patient localization, we expect improved clinical results to be forthcoming.

Quantification and visualization of charged particle range variations.

PURPOSE: Range variations during respiration affect the penetration of charged particle beams and can result in beam overshoot or undershoot to the target. We have developed analysis tools to quantify the water equivalent pathlength (WEL) variations resulting from respiration (WEL analyzer, Aqualyzer), as well as a data explorer to view WEL variations interactively. METHODS AND MATERIALS: The metrics to characterize and quantify penetration of a charged particle beam during respiration were calculated semiautomatically. The analysis involved the generation of images that (1) encode the radiologic pathlength across a beam's eye view image during the respiratory phase, (2) display the variation of the radiologic pathlength relative to a reference respiratory phase, (3) display isopenetration as a function of breathing, and (4) show range fluctuation images for a compensating bolus when applied to four-dimensional computed tomography. Additional quantities relevant to the analysis of charged particle beams in a breathing patient that are calculated include the beam overshoot volume and beam overshoot distance. These quantities are calculated as a function of time, gantry angle, and position. RESULTS: The software was applied to test cases to illustrate its utility in the analysis of range variations of charged particle beams in the treatment of lung tumors. CONCLUSION: WEL analysis is useful in the rapid assessment of range variations in the treatment of lung tumors and in determining the optimal gantry angle and respiratory gating window. The extension of encoding range fluctuations to a beam's eye view display is helpful in designing plans that are more robust in the presence of motion.

Maximum-intensity volumes for fast contouring of lung tumors including respiratory motion in 4DCT planning.

PURPOSE: To assess the accuracy of maximum-intensity volumes (MIV) for fast contouring of lung tumors including respiratory motion. METHODS AND MATERIALS: Four-dimensional computed tomography (4DCT) data of 10 patients were acquired. Maximum-intensity volumes were constructed by assigning the maximum Hounsfield unit in all CT volumes per geometric voxel to a new, synthetic volume. Gross tumor volumes (GTVs) were contoured on all CT volumes, and their union was constructed. The GTV with all its respiratory motion was contoured on the MIV as well. Union GTVs and GTVs including motion were compared visually. Furthermore, planning target volumes (PTVs) were constructed for the union of GTVs and the GTV on MIV. These PTVs were compared by centroid position, volume, geometric extent, and surface distance. RESULTS: Visual comparison of GTVs demonstrated failure of the MIV technique for 5 of 10 patients. For adequate GTV(MIV)s, differences between PTVs were <1.0 mm in centroid position, 5% in volume, +/-5 mm in geometric extent, and +/-0.5 +/- 2.0 mm in surface distance. These values represent the uncertainties for successful MIV contouring. CONCLUSION: Maximum-intensity volumes are a good first estimate for target volume definition including respiratory motion. However, it seems mandatory to validate each individual MIV by overlaying it on a movie loop displaying the 4DCT data and editing it for possible inadequate coverage of GTVs on additional 4DCT motion states.

Comparison of target registration errors for multiple image-guided techniques in accelerated partial breast irradiation.

PURPOSE: External beam accelerated partial breast irradiation requires accurate localization of the target volume for each treatment fraction. Using the concept of target registration error (TRE), the performance of several methods of target localization was compared. METHODS AND MATERIALS: Twelve patients who underwent external beam accelerated partial breast irradiation were included in this study. TRE was quantified for four methods of image guidance: standard laser-based setup, kilovoltage imaging of the chest wall, kilovoltage imaging of surgically implanted clips, and three-dimensional surface imaging of the breast. The use of a reference surface created from a free-breathing computed tomography scan and a reference surface directly captured with three-dimensional video imaging were compared. The effects of respiratory motion were also considered, and gating was used for 8 of 12 patients. RESULTS: The median value of the TRE for the laser, chest wall, and clip alignment was 7.1 mm (n=94), 5.4 mm (n=81), and 2.4 mm (n=93), respectively. The median TRE for gated surface imaging based on the first fraction reference surface was 3.2 mm (n=49), and the TRE for gated surface imaging using the computed tomography-based reference surface was 4.9 mm (n=56). The TRE for nongated surface imaging using the first fraction reference surface was 6.2 mm (n=25). CONCLUSIONS: The TRE of surface imaging using a reference surface directly captured with three-dimensional video and the TRE for clip-based setup were within 1 mm. Gated capture is important for surface imaging to reduce the effects of respiratory motion in accelerated partial breast irradiation.

Quantitative assessment of range fluctuations in charged particle lung irradiation.

PURPOSE: Water equivalent path length (WEL) variations due to respiration can change the range of a charged particle beam and result in beam overshoot to critical organs or beam undershoot to tumor. We have studied range fluctuations by analyzing four-dimensional computed tomography data and quantitatively assessing potential beam overshoot. METHODS AND MATERIALS: The maximal intensity volume is calculated by combining the gross tumor volume contours at each respiratory phase in the four-dimensional computed tomography study. The first target volume calculates the maximal intensity volume for the entire respiratory cycle (internal target volume [ITV]-radiotherapy [RT]), and the second target volume is the maximal intensity volume corresponding to gated RT (gated-RT, approximately 30% phase window around exhalation). A compensator at each respiratory phase is calculated. Two "composite" compensators for ITV-RT and gated-RT are then designed by selecting the minimal compensator depth at the respective respiratory phase. These compensators are then applied to the four-dimensional computed tomography data to estimate beam penetration. Analysis metrics include range fluctuation and overshoot volume, both as a function of gantry angle. We compared WEL fluctuations observed in treating the ITV-RT versus gated-RT in 11 lung patients. RESULTS: The WEL fluctuations were <21.8 mm-WEL and 9.5 mm-WEL for ITV-RT and gated-RT, respectively for all patients. Gated-RT reduced the beam overshoot volume by approximately a factor of four compared with ITV-RT. Such range fluctuations can affect the efficacy of treatment and result in an excessive dose to a distal critical organ. CONCLUSION: Time varying range fluctuation analysis provides information useful for determining appropriate patient-specific treatment parameters in charged particle RT. This analysis can also be useful for optimizing planning and delivery.

Effects of intrafractional motion on water equivalent pathlength in respiratory-gated heavy charged particle beam radiotherapy.

PURPOSE: To analyze the water equivalent pathlength (WEL) fluctuations resulting from cardiac motion and display these variations on a beam's-eye-view image; the analysis provides insight into the accuracy of lung tumor irradiation with heavy charged particle beams. MATERIALS AND METHODS: Volumetric cine computed tomography (CT) images were obtained on 7 lung cancer patients under free-breathing conditions with a 256-multislice CT scanner. Cardiac phase was determined by selecting systole and diastole. A WEL difference image (DeltaWEL) was calculated by subtracting the WEL image at end-systole from that at end-diastole at respiratory exhalation phase. Two calculation regions were defined: Region 1 was limited to the volume defined by planes bounding the heart; Region 2 included the entire body thickness for a given beam's-eye-view angle. RESULTS: The DeltaWEL values observed in Region 1 showed fluctuations at the periphery of the heart that varied from 20.4 (SD, 5.2) mm WEL to -15.6 (3.2) mm WEL. The areas over which these range perturbation values were observed were 36.8 (32.4) mm(2) and 6.0 (2.8) mm(2) for positive and negative WEL, respectively. The WEL fluctuations in Region 2 increased by approximately 3-4 mm WEL, whereas negative WEL fluctuations changed by approximately -4 to -5 mm WEL, compared with WEL for Region 1; areas over 20 mm WEL changes in Region 2 increased by 9 mm(2) for positive DeltaWEL and 2 mm(2) for negative DeltaWEL. CONCLUSIONS: Cine CT with a 256-multislice CT scanner captures both volumetric cardiac and respiratory motion with a temporal resolution sufficient to estimate range fluctuations by these motions. This information can be used to assess the range perturbations that charged particle beams may experience in irradiation of lung or esophageal tumors adjacent to the heart.

Four-dimensional imaging and treatment planning of moving targets.

Four-dimensional CT acquisition is commercially available, and provides important information on the shape and trajectory of the tumor and normal tissues. The primary advantage of four-dimensional imaging over light breathing helical scans is the reduction of motion artifacts during scanning that can significantly alter tumor appearance. Segmentation, image registration, visualization are new challenges associated with four-dimensional data sets because of the overwhelming increase in the number of images. Four-dimensional dose calculations, while currently laborious, provide insights into dose perturbations due to organ motion. Imaging before treatment (image guidance) improves accuracy of radiation delivery, and recording transmission images can provide a means of verifying gated delivery.

Assessing residual motion for gated proton-beam radiotherapy.

Gated radiation therapy is a promising method for improving the dose conformality of treatments to moving targets and reducing the total volume of irradiated tissue. Target motion is of particular concern in proton beam radiotherapy, due to the finite range of proton dose deposition in tissue. Gating allows one to reduce the extent of variation, due to respiration, of the radiological depth to target during treatment delivery. However, respiratory surrogates typically used for gating do not always accurately reflect the position of the internal target. For instance, a phase delay often exists between the internal motion and the motion of the surrogate. Another phenomenon, baseline drifting refers to a gradual change in the exhale position over time, which generally affects the external and internal markers differently. This study examines the influence of these two physiological phenomena on gated radiotherapy using an external surrogate.

Deformable registration of 4D computed tomography data.

Four-dimensional radiotherapy requires deformable registration to track delivered dose across varying anatomical states. Deformable registration based on B-splines was implemented to register 4D computed tomography data to a reference respiratory phase. To assess registration performance, anatomical landmarks were selected across ten respiratory phases in five patients. These point landmarks were transformed according to global registration parameters between different respiratory phases. Registration uncertainties were computed by subtraction of transformed and reference landmark positions. The selection of appropriate registration masks to separate independently moving anatomical subunits is crucial to registration performance. The average registration error for five landmarks for each of five patients was 2.1 mm. This level of accuracy is acceptable for most radiotherapy applications.

The susceptibility of IMRT dose distributions to intrafraction organ motion: an investigation into smoothing filters derived from four dimensional computed tomography data.

This study investigated the sensitivity of static planning of intensity-modulated beams (IMBs) to intrafraction deformable organ motion and assessed whether smoothing of the IMBs at the treatment-planning stage can reduce this sensitivity. The study was performed with a 4D computed tomography (CT) data set for an IMRT treatment of a patient with liver cancer. Fluence profiles obtained from inverse-planning calculations on a standard reference CT scan were redelivered on a CT scan from the 4D data set at a different part of the breathing cycle. The use of a nonrigid registration model on the 4D data set additionally enabled detailed analysis of the overall intrafraction motion effects on the IMRT delivery during free breathing. Smoothing filters were then applied to the beam profiles within the optimization process to investigate whether this could reduce the sensitivity of IMBs to intrafraction organ motion. In addition, optimal fluence profiles from calculations on each individual phase of the breathing cycle were averaged to mimic the convolution of a static dose distribution with a motion probability kernel and assess its usefulness. Results from nonrigid registrations of the CT scan data showed a maximum liver motion of 7 mm in superior-inferior direction for this patient. Dose-volume histogram (DVH) comparison indicated a systematic shift when planning treatment on a motion-frozen, standard CT scan but delivering over a full breathing cycle. The ratio of the dose to 50% of the normal liver to 50% of the planning target volume (PTV) changed up to 28% between different phases. Smoothing beam profiles with a median-window filter did not overcome the substantial shift in dose due to a difference in breathing phase between planning and delivery of treatment. Averaging of optimal beam profiles at different phases of the breathing cycle mainly resulted in an increase in dose to the organs at risk (OAR) and did not seem beneficial to compensate for organ motion compared with using a large margin. Additionally, the results emphasized the need for 4D CT scans when aiming to reduce the internal margin (IM). Using only a single planning scan introduces a systematic shift in the dose distribution during delivery. Smoothing beam profiles either based on a single scan or over the different breathing phases was not beneficial for reducing this shift.