Fluoroscopic evaluation of diaphragmatic motion reduction with a respiratory gated radiotherapy system

Accuracy and potential improvements of surface-guided breast cancer radiotherapy in deep inspiration breath-hold with daily image-guidance

Objective. Radiotherapy of left-sided breast cancer in deep inspiration breath-hold (DIBH) reduces the heart dose. Surface guided radiotherapy (SGRT) can guide the DIBH, but the accuracy is subject to variations in the chest wall position relative to the patient surface. Approach. In this study, ten left-sided breast cancer patients received DIBH radiotherapy with tangential fields in 15–18 fractions. After initial SGRT setup in free breathing an orthogonal MV/kV image pair was acquired during SGRT-guided breath-hold. The couch was corrected to align the chest wall during another breath-hold, and a new SGRT reference surface was acquired for the gating. The chest wall position error during treatment was determined from continuous cine MV images in the imager direction perpendicular to the cranio-caudal direction. A treatment error budget was made with individual contributions from the online registration of the setup MV image, the difference in breath-hold level between setup imaging and SGRT reference surface acquisition, the SGRT level during treatment, and intra-fraction shifts of the chest wall relative to the SGRT reference surface. In addition to the original setup protocol (Scenario A), SGRT was also simulated with better integration of image-guidance by capturing either the new reference surface (Scenario B) or the SGRT positional signal (Scenario C) simultaneously with the setup MV image, and accounting for the image-guided couch correction by shifting the SGRT reference surface digitally. Main results. In general, the external SGRT signal correlated well with the internal chest wall position error (correlation coefficient >0.7 for 75% of field deliveries), but external-to-internal target position offsets above 2 mm occasionally occurred (13% of fractions). The PTV margin required to account for the treatment error was 3.5 mm (Scenario A), 3.4 mm (B), and 3.1 mm (C). Significance. Further integration of SGRT with image-guidance may improve treatment accuracy and workflow although the current study did not show large accuracy improvements of scenario B and C compared to scenario A.

Potential of a probabilistic framework for target prediction from surrogate respiratory motion during lung radiotherapy

Purpose.Respiration-induced motion introduces significant positioning uncertainties in radiotherapy treatments for thoracic sites. Accounting for this motion is a non-trivial task commonly addressed with surrogate-based strategies and latency compensating techniques. This study investigates the potential of a new unified probabilistic framework to predict both future target motion in real-time from a surrogate signal and associated uncertainty.Method.A Bayesian approach is developed, based on a Kalman filter theory adapted specifically for surrogate measurements. Breathing motions are collected simultaneously from a lung target, two external surrogates (abdominal and thoracic markers) and an internal surrogate (liver structure) for 9 volunteers during 4 min, in which severe breathing changes occur to assess the robustness of the method. A comparison with an artificial non-linear neural network (NN) is performed, although no confidence interval prediction is provided. A static worst-case scenario and a simple static design are investigated.Results.Although the NN can reduce the prediction errors from thoracic surrogate in some cases, the Bayesian framework outperforms in most cases the NN when using the other surrogates: bias on predictions is reduced by 38% and 16% on average when using respectively the liver and the abdomen for the simple scenario, and by respectively 40% and 31% for the worst-case scenario. The standard deviation of residuals is reduced on average by up to 42%. The Bayesian method is also found to be more robust to increasing latencies. The thoracic marker appears to be less reliable to predict the target position, while the liver shows to be a better surrogate. A statistical test confirms the significance of both observations.Conclusion.The proposed framework predicts both the future target position and the associated uncertainty, which can be valuably used to further assist motion management decisions. Further investigation is required to improve the predictions by using an adaptive version of the proposed framework.

Gated Radiotherapy Development and its Expansion

One of the most important challenges in treatment of patients with cancerous tumors of chest and abdominal areas is organ movement. The delivery of treatment radiation doses to tumor tissue is a challenging matter while protecting healthy and radio sensitive tissues. Since the movement of organs due to respiration causes a discrepancy in the middle of planned and delivered dose distributions. The moderation in the fatalistic effect of intra-fractional target travel on the radiation therapy correctness is necessary for cutting-edge methods of motion remote monitoring and cancerous growth irradiancy. Tracking respiratory milling and implementation of breath-hold techniques by respiratory gating systems have been used for compensation of respiratory motion negative effects. Therefore, these systems help us to deliver precise treatments and also protect healthy and critical organs. It seems aspiration should be kept under observation all over treatment period employing tracking seed markers (e.g. fiducials), skin surface scanners (e.g. camera and laser monitoring systems) and aspiration detectors (e.g. spirometers). However, these systems are not readily available for most radiotherapy centers around the word. It is believed that providing and expanding the required equipment, gated radiotherapy will be a routine technique for treatment of chest and abdominal tumors in all clinical radiotherapy centers in the world by considering benefits of respiratory gating techniques in increasing efficiency of patient treatment in the near future. This review explains the different technologies and systems as well as some strategies available for motion management in radiotherapy centers.

Technical Note: Investigating internal-external motion correlation using fast helical CT

PURPOSE

To quantify the use of anterior torso skin surface position measurement as a breathing surrogate.

METHODS

Fourteen patients were scanned 25 times in alternating directions using a free-breathing low-mA fast helical CT protocol. Simultaneously, an abdominal pneumatic bellows was used as a real-time breathing surrogate. The imaged diaphragm dome position was used as a gold standard surrogate, characterized by localizing the most superior points of the diaphragm dome in each lung. These positions were correlated against the bellows signal acquired at the corresponding scan times. The bellows system has been shown to have a slow linear drift, and the bellows-to-CT synchronization process had a small uncertainty, so the drift and time offset were determined by maximizing the correlation coefficient between the craniocaudal diaphragm position and the drift-corrected bellows signal. The corresponding fit was used to model the real-time diaphragm position. To estimate the effectiveness of skin surface positions as surrogates, the anterior torso surface position was measured from the CT scans and correlated against the diaphragm position model. The residual error was defined as the root-mean-square correlation residual with the breathing amplitude normalized to the 5^th to 95^th breathing amplitude percentiles. The fit residual errors were analyzed over the surface for the fourteen studied patients and reported as percentages of the 5^th to 95^th percentile ranges.

RESULTS

A strong correlation was measured between the diaphragm motion and the abdominal bellows signal with an average residual error of 9.21% and standard deviation of 3.77%. In contrast, the correlations between the diaphragm position model and patient surface positions varied throughout the torso and from patient to patient. However, a consistently high correlation was found near the abdomen for each patient, and the average minimum residual error relating the skin surface to the diaphragm was 11.8% with a standard deviation of 4.61%.

CONCLUSIONS

The thoracic patient surface was found to be an accurate surrogate, but the accuracy varied across the surface sufficiently that care would need to be taken to use the surface as an accurate and reliable surrogate. Future studies will use surface imaging to determine surface patch algorithms that utilize the entire chest as well as thoracic and abdominal breathing relationships.

Couch and multileaf collimator tracking: A clinical feasibility study for pancreas and liver treatment

PURPOSE

Real-time tumor tracking through active correction by the multileaf collimator or treatment couch offers a promising strategy to mitigate delivery uncertainty due to intrafractional tumor motion. This study evaluated the performance of MLC and couch tracking using the prototype iTools Tracking system in TrueBeam Developer Mode and the application for abdominal cancer treatments.

METHODS

Experiments were carried out using a phantom with embedded Calypso transponders and a motion simulation platform. Geometric evaluations were performed using a circular conformal field with sinusoidal traces and pancreatic tumor motion traces. Geometric tracking accuracy was retrospectively calculated by comparing the compensational MLC or couch motion extracted from machine log files to the target motion reconstructed from real-time MV and kV images. Dosimetric tracking accuracy was measured with radiochromic films using clinical abdominal VMAT plans and pancreatic tumor traces.

RESULTS

Geometrically, the root-mean-square errors for MLC tracking were 0.5 and 1.8 mm parallel and perpendicular to leaf travel direction, respectively. Couch tracking, in contrast, showed an average of 0.8 mm or less geometric error in all directions. Dosimetrically, both MLC and couch tracking reduced motion-induced local dose errors compared to no tracking. Evaluated with five pancreatic tumor motion traces, the average 2%/2 mm global gamma pass rate of eight clinical abdominal VMAT plans was 67.4% (range: 26.4%-92.7%) without tracking, which was improved to 86.0% (range: 67.9%-95.6%) with MLC tracking, and 98.1% (range: 94.9%-100.0%) with couch tracking. In 16 out of 40 deliveries with different plans and motion traces, MLC tracking did not achieve clinically acceptable dosimetric accuracy with 3%/3mm gamma pass rate below 95%.

CONCLUSIONS

This study demonstrated the capability of MLC and couch tracking to reduce motion-induced dose errors in abdominal cases using a prototype tracking system. Clinically significant dose errors were observed with MLC tracking for certain plans which could be attributed to the inferior MLC tracking accuracy in the direction perpendicular to leaf travel, as well as the interplay between motion tracking and plan delivery for highly modulated plans. Couch tracking outperformed MLC tracking with consistently high dosimetric accuracy in all plans evaluated, indicating its clinical potential in the treatment of abdominal cancers.

The effect of respiratory motion on electronic portal imaging device dosimetry

There is an increasing need to develop methods for in vivo verification of the delivery of radiotherapy treatments. Electronic portal imaging devices (EPID's) have been demonstrated to be of use for this application. The basic principle is relatively straightforward, the EPID is used to measure a two-dimensional (2D) planar exit or portal dose map behind the patient during the treatment delivery that can provide information on any errors in linear accelerator output or changes in the patient anatomy. In this paper we focused on the effect of intra-fraction motion, particularly respiratory motion, on the measured 2D EPID dose-response. Measurements were made with a breast phantom undergoing one-dimensional (1D) sinusoidal motion with a range of amplitudes (0.5, 1.0, and 1.5 cm) and frequencies (12, 15, and 20  cycles/min). Further measurements were made with the phantom undergoing breathing sequences measured during patient planning computed tomography simulation. We made use of the quadratic calibration method that converts the EPID images to a surrogate for dose, equivalent thickness of Plastic Water^® . Comparisons were made of the 2D thickness maps derived for the different motions compared to the static phantom case and the resulting dose difference analyzed over the "breast" region of interest. A 2D gamma analysis within the same region of interest was performed of the motion images compared to static reference image. Comparisons were made of 1D thickness profiles for the moving and static phantom. The 1D and 2D analyses show the method to be sensitive to the smallest motion amplitude of 0.5 cm tested in the phantom measurements. The results using the phantom demonstrate the method to be a potentially useful tool for monitoring intra-fraction motion during the delivery of patient radiotherapy treatments as well as more generally providing information on the effects of motion on EPID based in vivo dosimetric verification.

Variations of MRI-assessed peristaltic motions during radiation therapy

PURPOSE

Understanding complex abdominal organ motion is essential for motion management in radiation therapy (RT) of abdominal tumors. This study investigates abdominal motion induced by respiration and peristalsis, during various time durations relevant to RT, using various CT and MRI techniques acquired under free breathing (FB) and breath hold (BH).

METHODS

A series of CT and MRI images acquired with various techniques under free breathing and/or breath hold from 37 randomly-selected pancreatic or liver cancer patients were analyzed to assess the motion in various time frames. These data include FB 4DCT from 15 patients (for motion in time duration of 5 sec), FB 2D cine-MRI from 4 patients (time duration of 1.7 min, 1 second acquisition time per slice), FB cine-MRI acquired using MR-Linac from 6 patients in various fractions (acquisition time is less than 0.6 seconds per slice), FB 4DMRI from 2 patients (time duration of 2 min), respiration-gated T2 with gating at the end expiration (time duration of 3-5 min), and BH T1 with multiphase dynamic contrast in acquisition times of 17 seconds for each of five phases (pre-contrast, arterial, venous, portal venous and delayed post-contrast) from 10 patients. Motions of various organs including gallbladder (GB) and liver were measured based on these MRI data. The GB motion includes both respiration and peristalsis, while liver motion is primarily respiration. By subtracting liver motion (respiration) from GB motion (respiration and peristalsis), the peristaltic motion, along with small residual motion, was obtained.

RESULTS

From cine-MRI, the residual motion beyond the respiratory motion was found to be up to 0.6 cm in superior-inferior (SI) and 0.55 cm in anterior-posterior (AP) directions. From 2D cine-MRI acquired by the MR-Linac, different peristaltic motions were found from different fractions for each patient. The peristaltic motion was found to vary between 0.3-1 cm. From BH T1 phase images, the average motions that were primarily due to peristalsis movements were found to be 1.2 cm in SI, 0.7 cm in AP, and 0.9 cm in left-right (LR) directions. The average motions assessed from 4DCT were 1.0 cm in SI and 0.3 cm in AP directions, which were generally smaller than the motions assessed from cine-MRI, i.e., 1.8 cm in SI and 0.6 cm in AP directions, for the same patients. However, average motions from 4DMRI, which are coming from respiratory were measured to be 1.5, 0.5, and 0.4 cm in SI, AP, and LR directions, respectively.

CONCLUSION

The abdominal motion due to peristalsis can be similar in magnitude to respiratory motion as assessed. These motions can be irregular and persistent throughout the imaging and RT delivery procedures, and should be considered together with respiratory motion during RT for abdominal tumors.

Respiratory-gated PET/CT for pulmonary lesion characterisation-promises and problems

2-deoxy-2-(¹⁸Fluorine)-fluoro-D-glucose (FDG) PET/CT is an integral part of lung carcinoma staging and frequently used in the assessment of solitary pulmonary nodules. However, a limitation of conventional three-dimensional PET/CT when imaging the thorax is its susceptibility to motion artefact, which blurs the signal from the lesion resulting in inaccurate representation of size and metabolic activity. Respiratory gated (four-dimensional) PET/CT aims to negate the effects of motion artefact and provide a more accurate interpretation of pulmonary nodules and lymphadenopathy. There have been recent advances in technology and a shift from traditional hardware to more streamlined software methods for respiratory gating which should allow more widespread use of respiratory-gating in the future. The purpose of this article is to review the evidence surrounding four-dimensional PET/CT in pulmonary lesion characterisation.

Motion management strategies and technical issues associated with stereotactic body radiotherapy of thoracic and upper abdominal tumors: A review from NRG oncology

The efficacy of stereotactic body radiotherapy (SBRT) has been well demonstrated. However, it presents unique challenges for accurate planning and delivery especially in the lungs and upper abdomen where respiratory motion can be significantly confounding accurate targeting and avoidance of normal tissues. In this paper, we review the current literature on SBRT for lung and upper abdominal tumors with particular emphasis on addressing respiratory motion and its affects. We provide recommendations on strategies to manage motion for different, patient-specific situations. Some of the recommendations will potentially be adopted to guide clinical trial protocols.

Treatment planning and evaluation of gated radiotherapy in left-sided breast cancer patients using the CatalystTM/SentinelTM system for deep inspiration breath-hold (DIBH)

BACKGROUND

There is a potential for adverse cardiovascular effects in long-term breast cancer survivors following adjuvant radiotherapy (RT). For this purpose, the deep inspiration breath-hold technique (DIBH) has been introduced into clinical practice, to maximally reduce the radiation dose to the heart. However, there are a variety of DIBH delivery techniques, patient positioning and visual patient feedback mechanisms. The aim of the present study was to evaluate the application of radiotherapy in DIBH using the Catalyst^TM/Sentinel^TM system, with a special emphasis on treatment planning and dosimetric plan comparison in free breathing (FB) and DIBH.

PATIENTS AND METHODS

A total of 13 patients with left-sided breast cancer following breast conserving surgery were included in this prospective clinical trial. For treatment application the Catalyst^TM/Sentinel^TM system (C-RAD AB, Uppsala, Sweden) was used and gating control was performed by an audio-visual patient feedback system. CT and surface data were acquired in FB and DIBH and dual treatment plans were created using Pencil Beam and Collapsed Cone Convolution. Dosimetric output parameters of organs at risk were compared using Wilcoxon signed-rank test. Central lung distance (CLD) was retrieved from iView^TM portal images during treatment delivery.

RESULTS

The system contains a laser surface scanner (Sentinel^TM) and an optical surface scanner (Catalyst^TM) interconnected to the LINAC systems via a gating interface and allows for a continuous and touchless surface scanning. Overall, 225 treatment fractions with audio-visual guidance were completed without any substantial difficulties. Following initial patient training and treatment setup, radiotherapy in DIBH with the Catalyst^TM/Sentinel^TM system was time-efficient and reliable. Following dual treatment planning for all patients, nine of 13 patients were treated in DIBH. In these patients, the reduction of the mean heart dose for DIBH compared to FB was 52 % (2.73 to 1.31 Gy; p = 0.011). The maximum doses to the heart and LAD were reduced by 59 % (47.90 to 19.74 Gy; p = 0.008) and 75 % (38.55 to 9.66 Gy; p = 0.008), respectively. In six of the nine patients the heart completely moved out of the treatment field by DIBH. The standard deviation of the CLD varied between 0.12 and 0.29 cm (mean: 0.16 cm).

CONCLUSION

The Catalyst^TM/Sentinel^TM system enabled a fast and reliable application and surveillance of DIBH in daily clinical routine. Furthermore, the present data show that using the DIBH technique during RT could significantly reduce high dose areas and mean doses to the heart.

TRIAL REGISTRATION

DRKS: DRKS00010929 registered on 5. August 2016.

Lung surface deformation prediction from spirometry measurement and chest wall surface motion

PURPOSE

The authors have developed and evaluated a method to predict lung surface motion based on spirometry measurements, and chest and abdomen motion at selected locations.

METHODS

A patient-specific 3D triangular surface mesh of the lung region was obtained at the end expiratory phase by the threshold-based segmentation method. Lung flow volume changes were recorded with a spirometer for each patient. A total of 192 selected points at a regular spacing of 2 × 2 cm matrix points were used to detect chest wall motion over a total area of 32 × 24 cm covering the chest and abdomen surfaces. QR factorization with column pivoting was employed to remove redundant observations of the chest and abdominal areas. To create a statistical model between the lung surface and the corresponding surrogate signals, the authors developed a predictive model based on canonical ridge regression. Two unique weighting vectors were selected for each vertex on the lung surface; they were optimized during the training process using all other 4D-CT phases except for the test inspiration phase. These parameters were employed to predict the vertex locations of a testing data set.

RESULTS

The position of each lung surface mesh vertex was estimated from the motion at selected positions within the chest wall surface and from spirometry measurements in ten lung cancer patients. The average estimation of the 98th error percentile for the end inspiration phase was less than 1 mm (AP = 0.9 mm, RL = 0.6 mm, and SI = 0.8 mm). The vertices located at the lower region of the lung had a larger estimation error as compared with those within the upper region of the lung. The average landmark motion errors, derived from the biomechanical modeling using real surface deformation vector fields (SDVFs), and the predicted SDVFs were 3.0 and 3.1 mm, respectively.

CONCLUSIONS

Our newly developed predictive model provides a noninvasive approach to derive lung boundary conditions. The proposed system can be used with personalized biomechanical respiration modeling to derive lung tumor motion during radiation therapy from noninvasive measurements.

Optimum location of external markers using feature selection algorithms for real-time tumor tracking in external-beam radiotherapy: a virtual phantom study

In external-beam radiotherapy, using external markers is one of the most reliable tools to predict tumor position, in clinical applications. The main challenge in this approach is tumor motion tracking with highest accuracy that depends heavily on external markers location, and this issue is the objective of this study. Four commercially available feature selection algorithms entitled 1) Correlation-based Feature Selection, 2) Classifier, 3) Principal Components, and 4) Relief were proposed to find optimum location of external markers in combination with two "Genetic" and "Ranker" searching procedures. The performance of these algorithms has been evaluated using four-dimensional extended cardiac-torso anthropomorphic phantom. Six tumors in lung, three tumors in liver, and 49 points on the thorax surface were taken into account to simulate internal and external motions, respectively. The root mean square error of an adaptive neuro-fuzzy inference system (ANFIS) as prediction model was considered as metric for quantitatively evaluating the performance of proposed feature selection algorithms. To do this, the thorax surface region was divided into nine smaller segments and predefined tumors motion was predicted by ANFIS using external motion data of given markers at each small segment, separately. Our comparative results showed that all feature selection algorithms can reasonably select specific external markers from those segments where the root mean square error of the ANFIS model is minimum. Moreover, the performance accuracy of proposed feature selection algorithms was compared, separately. For this, each tumor motion was predicted using motion data of those external markers selected by each feature selection algorithm. Duncan statistical test, followed by F-test, on final results reflected that all proposed feature selection algorithms have the same performance accuracy for lung tumors. But for liver tumors, a correlation-based feature selection algorithm, in combination with a genetic search algorithm, proved to yield best performance accuracy for selecting optimum markers.

Technical Note: Intrafractional changes in time lag relationship between anterior-posterior external and superior-inferior internal motion signals in abdominal tumor sites

PURPOSE

To investigate constancy, within a treatment session, of the time lag relationship between implanted markers in abdominal tumors and an external motion surrogate.

METHODS

Six gastroesophageal junction and three pancreatic cancer patients (IRB-approved protocol) received two cone-beam CTs (CBCT), one before and one after treatment. Time between scans was less than 30 min. Each patient had at least one implanted fiducial marker near the tumor. In all scans, abdominal displacement (Varian RPM) was recorded as the external motion signal. Purpose-built software tracked fiducials, representing internal signal, in CBCT projection images. Time lag between superior-inferior (SI) internal and anterior-posterior external signals was found by maximizing the correlation coefficient in each breathing cycle and averaging over all cycles. Time-lag-induced discrepancy between internal SI position and that predicted from the external signal (external prediction error) was also calculated.

RESULTS

Mean ± standard deviation time lag, over all scans and patients, was 0.10 ± 0.07 s (range 0.01-0.36 s). External signal lagged the internal in 17/18 scans. Change in time lag between pre- and post-treatment CBCT was 0.06 ± 0.07 s (range 0.01-0.22 s), corresponding to 3.1% ± 3.7% (range 0.6%-10.8%) of gate width (range 1.6-3.1 s). In only one patient, change in time lag exceeded 10% of the gate width. External prediction error over all scans of all patients varied from 0.1 ± 0.1 to 1.6 ± 0.4 mm.

CONCLUSIONS

Time lag between internal motion along SI and external signals is small compared to the treatment gate width of abdominal patients examined in this study. Change in time lag within a treatment session, inferred from pre- to post-treatment measurements is also small, suggesting that a single measurement of time lag at the session start is adequate. These findings require confirmation in a larger number of patients.

Motion management in gastrointestinal cancers

The presence of tumor and organ motions complicates the planning and delivery of radiotherapy for gastrointestinal cancers. Without proper accounting of the movements, target volume could be under-dosed and the nearby normal critical organs could be over-dosed. This situation is further exacerbated by the close proximity of abdominal tumors to many normal organs at risk (OARs). A number of strategies have been developed to deal with tumor and organ motions in radiotherapy. This article presents a review of the techniques used in the evaluation, quantification, and management of tumor and organ motions for radiotherapy of gastrointestinal cancers.

Diaphragmatic motion studied by M-mode ultrasonography in combined pulmonary fibrosis and emphysema

BACKGROUND

The coexistence of emphysema and pulmonary fibrosis is known as combined pulmonary fibrosis and emphysema (CPFE). The aim of this study was to compare diaphragmatic motion measured by M-mode ultrasonography of patients with CPFE, idiopathic pulmonary fibrosis (IPF), and chronic obstructive pulmonary disease (COPD).

METHODS

Pulmonary function, high-resolution computed tomography (HRCT), and diaphragmatic motion were examined in patients with CPFE (n = 25), IPF (n = 18), and COPD (n = 60), and in healthy controls (n = 21). Diaphragmatic motions were measured on M-mode ultrasonographic images during quiet breathing and deep breathing.

RESULTS

There were no significant differences in right or left diaphragmatic motion during quiet breathing among the four groups, whereas differences were significant in right and left motion during deep breathing. Diaphragmatic motion in CPFE patients was the lowest among the four groups. COPD patients, especially those with severe COPD, showed significantly lower diaphragmatic motion than IPF patients or healthy controls. There were no differences in diaphragmatic motion between IPF patients and healthy controls. Right diaphragmatic motions during deep breathing were negatively correlated with emphysema scores (r = -0.606, p < 0.001), but were not correlated with fibrosis scores on HRCT.

CONCLUSIONS

Diaphragmatic weakness was found in CPFE patients. Emphysema but not fibrosis may be one cause of limited diaphragmatic motion in patients with CPFE. M-mode ultrasonographic evaluation of diaphragmatic motion during deep breathing may be a useful tool in diagnosing CPFE and in discriminating CPFE patients from IPF or COPD patients.

Improving radiotherapy planning, delivery accuracy, and normal tissue sparing using cutting edge technologies

In the United States, more than half of all new invasive cancers diagnosed are non-small cell lung cancer, with a significant number of these cases presenting at locally advanced stages, resulting in about one-third of all cancer deaths. While the advent of stereotactic ablative radiation therapy (SABR, also known as stereotactic body radiotherapy, or SBRT) for early-staged patients has improved local tumor control to >90%, survival results for locally advanced stage lung cancer remain grim. Significant challenges exist in lung cancer radiation therapy including tumor motion, accurate dose calculation in low density media, limiting dose to nearby organs at risk, and changing anatomy over the treatment course. However, many recent technological advancements have been introduced that can meet these challenges, including four-dimensional computed tomography (4DCT) and volumetric cone-beam computed tomography (CBCT) to enable more accurate target definition and precise tumor localization during radiation, respectively. In addition, advances in dose calculation algorithms have allowed for more accurate dosimetry in heterogeneous media, and intensity modulated and arc delivery techniques can help spare organs at risk. New delivery approaches, such as tumor tracking and gating, offer additional potential for further reducing target margins. Image-guided adaptive radiation therapy (IGART) introduces the potential for individualized plan adaptation based on imaging feedback, including bulky residual disease, tumor progression, and physiological changes that occur during the treatment course. This review provides an overview of the current state of the art technology for lung cancer volume definition, treatment planning, localization, and treatment plan adaptation.

Development of real-time motion verification system using in-room optical images for respiratory-gated radiotherapy

Phase‐based respiratory‐gated radiotherapy relies on the reproducibility of patient breathing during the treatment. To monitor the positional reproducibility of patient breathing against a 4D CT simulation, we developed a real‐time motion verification system (RMVS) using an optical tracking technology. The system in the treatment room was integrated with a real‐time position management system. To test the system, an anthropomorphic phantom that was mounted on a motion platform moved on a programmed breathing pattern and then underwent a 4D CT simulation with RPM. The phase‐resolved anterior surface lines were extracted from the 4D CT data to constitute 4D reference lines. In the treatment room, three infrared reflective markers were attached on the superior, middle, and inferior parts of the phantom along with the body midline and then RMVS could track those markers using an optical camera system. The real‐time phase information extracted from RPM was delivered to RMVS via in‐house network software. Thus, the real‐time anterior‐posterior positions of the markers were simultaneously compared with the 4D reference lines. The technical feasibility of RMVS was evaluated by repeating the above procedure under several scenarios such as ideal case (with identical motion parameters between simulation and treatment), cycle change, baseline shift, displacement change, and breathing type changes (abdominal or chest breathing). The system capability for operating under irregular breathing was also investigated using real patient data. The evaluation results showed that RMVS has a competence to detect phase‐matching errors between patient's motion during the treatment and 4D CT simulation. Thus, we concluded that RMVS could be used as an online quality assurance tool for phase‐based gating treatments.

PACS number: 87.55.Qr

Abdominal and pancreatic motion correlation using 4D CT, 4D transponders, and a gating belt

The correlation between the pancreatic and external abdominal motion due to respiration was investigated on two patients. These studies utilized four dimensional computer tomography (4D CT), a four dimensional (4D) electromagnetic transponder system, and a gating belt system. One 4D CT study was performed during simulation to quantify the pancreatic motion using computer tomography images at eight breathing phases. The motion under free breathing and breath‐hold were analyzed for the 4D electromagnetic transponder system and the gating belt system during treatment. A linear curve was fitted for all data sets and correlation factors were evaluated between the 4D electromagnetic transponder system and the gating belt system data. The 4D CT study demonstrated a modest correlation between the external marker and the pancreatic motion with R‐square values larger than 0.8 for the inferior–superior (inf‐sup). Then, the relative pressure from the belt gating system correlated well with the 4D electromagnetic transponder system's motion in the anterior–posterior (ant‐post) and the inf–post directions. These directions have a correlation value of −0.93 and 0.76, while the lateral only had a 0.03 correlation coefficient. Based on our limited study, external surrogates can be used as predictors of the pancreatic motion in the inf–sup and the ant–post directions. Although there is a low correlation on the lateral direction, its motion is significantly shorter. In conclusion, an appropriate treatment delivery can be used for pancreatic cancer when an internal tracking system, such as the 4D electromagnetic transponder system, is unavailable.

PACS number: 87.55.kh

Impact of PET - CT motion correction in minimizing the gross tumor volume in non-small cell lung cancer

OBJECTIVES

To investigate the impact of respiratory motion on localization, and quantification of lung lesions for the Gross Tumor Volume utilizing a fully automated Auto3Dreg program and dynamic NURBS-based cardiac-torso digitized phantom (NCAT).

METHODS

Respiratory motion may result in more than 30% underestimation of the SUV values of lung, liver and kidney tumor lesions. The motion correction technique adopted in this study was an image-based motion correction approach using, a voxel-intensity-based and a multi-resolution multi-optimization (MRMO) algorithm. The NCAT phantom was used to generate CT attenuation maps and activity distribution volumes for the lung regions. All the generated frames were co-registered to a reference frame using a time efficient scheme. Quantitative assessment including Region of Interest (ROI), image fidelity and image correlation techniques, as well as semi-quantitative line profile analysis and qualitatively overlaying non-motion and motion corrected image frames were performed.

RESULTS

The largest motion was observed in the Z-direction. The greatest translation was for the frame 3, end inspiration, and the smallest for the frame 5 which was closet frame to the reference frame at 67% expiration. Visual assessment of the lesion sizes, 20-60mm at 3 different locations, apex, mid and base of lung showed noticeable improvement for all the foci and their locations. The maximum improvements for the image fidelity were from 0.395 to 0.930 within the lesion volume of interest. The greatest improvement in activity concentration underestimation was 7.7% below the true activity for the 20 mm lesion in comparison to 34.4% below, prior to correction. The discrepancies in activity underestimation were reduced with increasing the lesion sizes. Overlaying activity distribution on the attenuation map showed improved localization of the PET metabolic information to the anatomical CT images.

CONCLUSION

The respiratory motion correction for the lung lesions has led to an improvement in the lesion size, localization and activity quantification with a potential application in reducing the size of the PET GTV for radiotherapy treatment planning applications and hence improving the accuracy of the regime in treatment of lung cancer.

Challenges and opportunities in patient-specific, motion-managed and PET/CT-guided radiation therapy of lung cancer: review and perspective

The increasing interest in combined positron emission tomography (PET) and computed tomography (CT) to guide lung cancer radiation therapy planning has been well documented. Motion management strategies during treatment simulation PET/CT imaging and treatment delivery have been proposed to improve the precision and accuracy of radiotherapy. In light of these research advances, why has translation of motion-managed PET/CT to clinical radiotherapy been slow and infrequent? Solutions to this problem are as complex as they are numerous, driven by large inter-patient variability in tumor motion trajectories across a highly heterogeneous population. Such variation dictates a comprehensive and patient-specific incorporation of motion management strategies into PET/CT-guided radiotherapy rather than a one-size-fits-all tactic. This review summarizes challenges and opportunities for clinical translation of advances in PET/CT-guided radiotherapy, as well as in respiratory motion-managed radiotherapy of lung cancer. These two concepts are then integrated into proposed patient-specific workflows that span classification schemes, PET/CT image formation, treatment planning, and adaptive image-guided radiotherapy delivery techniques.

Magnitude of shift of tumor position as a function of moderated deep inspiration breath-hold: An analysis of pooled data of lung patients with active breath control in image-guided radiotherapy

The purpose of this study was to evaluate the reproducibility and magnitude of shift of tumor position by using active breathing control and iView-GT for patients with lung cancer with moderate deep-inspiration breath-hold (mDIBH) technique. Eight patients with 10 lung tumors were studied. CT scans were performed in the breath-holding phase. Moderate deep-inspiration breath-hold under spirometer-based monitoring system was used. Few important bony anatomic details were delineated by the radiation oncologist. To evaluate the interbreath-hold reproducibility of the tumor position, we compared the digital reconstruction radiographs (DRRs) from planning system with the DRRs from the iView-GT in the machine room. We measured the shift in x, y, and z directions. The reproducibility was defined as the difference between the bony landmarks from the DRR of the planning system and those from the DRR of the iView-GT. The maximum shift of the tumor position was 3.2 mm, 3.0 mm, and 2.9 mm in the longitudinal, lateral, and vertical directions. In conclusion, the moderated deep-inspiration breath-hold method using a spirometer is feasible, with relatively good reproducibility of the tumor position for image-guided radiotherapy in lung cancers.

Characterization and identification of spatial artifacts during 4D-CT imaging

PURPOSE

The purpose of this work is twofold: First, to characterize the artifacts occurring in helical 4D-CT imaging; second, to propose a method that can automatically identify the artifacts in 4D-CT images. The authors have designed a process that can automatically identify the artifacts in 4D-CT images, which may be invaluable in quantifying the quality of 4D-CT images and reducing the artifacts from the reconstructed images on a large dataset.

METHODS

Given two adjacent stacks obtained from the same respiration phase, the authors determine if there are artifacts between them. The proposed method uses a "bridge" stack strategy to connect the two stacks. Using normalized cross correlation convolution (NCCC), the two stacks are mapped to the bridge stack and the best matching positions can be located. Using this position information, the authors can then determine if there are artifacts between the two stacks. By combining the matching positions with NCCC values, the performance can be improved.

RESULTS

To validate the method, three expert observers independently labeled over 600 stacks on five patients. The results confirmed that high performance was obtained using the proposed method. The average sensitivity was about 0.87 and the average specificity was 0.82. The proposed method also outperformed the method of using respiratory signal (sensitivity increased from 0.50 to 0.87 and specificity increased from 0.70 to 0.82).

CONCLUSIONS

This study shows that the spatial artifacts during 4D-CT imaging are characterized and can be located automatically by the proposed method. The method is relatively simple but effective. It provides a way to evaluate the artifacts more objectively and accurately. The reported approach has promising potential for automatically identifying the types and frequency of artifacts on large scale 4D-CT image dataset.

To study tumor motion and planning target volume margins using four dimensional computed tomography for cancer of the thorax and abdomen regions

In this study, four dimensional computed tomography (4DCT) scanning was performed during free breathing on a 16-slice Positron emission tomography PET /computed tomography (CT) for abdomen and thoracic patients. Images were sorted into 10 phases based on the temporal correlation between surface motion and data acquisition with an Advantage Workstation. Gross tumor volume gross tumor volume (GTV) s were manually contoured on all 10 phases of the 4DCT scan. GTVs in the multiple CT phases were called GTV4D. GTV4D plus an isotropic margin of 1.0 cm was called CTV4D. Two sets of planning target volume (PTV) 4D (PTV4D) were derived from the CTV4D, i.e. PTV4D_2cm = CTV4D plus 1 cm setup margin (SM) and 1 cm internal margin (IM) and PTV4D_1.5cm = CTV4D plus 1 cm SM and 0.5cm IM. PTV3D was derived from a CTV3D of the helical CT scan plus conventional margins of 2 cm. PTV_gated was generated only selecting three CT phases, with a total margin of 1.5 cm. All four volumes were compared. To quantify the extent of the motion, we selected the two phases where the tumor exhibited the greatest range of motion. We also studied the effect of different PTV volumes on dose to the surrounding critical structures. Volume of CTV4D was greater than that of CTV3D. We found, on an average, a reduction of 14% volume of PTV4D_1.5cm as compared with PTV3D and reduction of 10% volume of PTV_gated as compared with PTV4D_1.5cm. We found that 2 cm of margin was inadequate if true motion of tumor was not known. We observed greater sparing of critical structures for PTVs drawn taking into account the tumor motion.

The effect of an inconsistent breathing amplitude on the relationship between an external marker and internal lung deformation in a porcine model

PURPOSE

Investigate the relationship between the motion of the Varian Real-time Position Management (RPM) device and the internal motion of a pig during induced inconsistencies in the amplitude of breathing.

METHODS

Twelve studies were performed on four ventilated female Landrace cross pigs using a GE Healthcare, Discovery CT 750 HD scanner. In each study, a 4.0 cm section (64 slices) of the pig's lungs was repeatedly scanned 20 times using cine mode, each time lasting more than one breathing cycle. During these cine scans, a Varian RPM device was used to collect respiratory amplitudes and the ventilator air return tube was periodically crimped to induce inconsistent breathing amplitudes. Each breathing cycle and its associated cine scan were categorized as either consistent or inconsistent, based on thresholds of the minimum expiration and maximum inspiration amplitudes. From the group of consistent amplitude cine scans in a study, a reference scan was chosen. The effect of inconsistent breathing amplitudes on the relationship between the motion of the RPM marker and the motion within three regions of interest (in each lung and the chest wall) was investigated with two methods: (1) A 4D-CT sorting algorithm based on RPM amplitude was used to sort volumes into 4D-CT phase bins. Within each phase bin, the nonlinear deformation of volumes collected during consistent and inconsistent breathing amplitude was calculated with respect to the reference volume. The magnitude of the deformations (in mm) were compared to determine if inconsistent breathing amplitude caused greater deformations. (2) Nonlinear deformations between each CT volume from a cine scan and the maximum expiration volume of the reference scan were calculated. Regression analyses between the nonlinear deformations within three regions of interest (in each lung and the chest wall) and the RPM amplitudes were performed to test the effect of inconsistent breathing amplitudes on the linearity of the relationship between the 3D motion of internal anatomy and the 1D motion of the RPM external marker.

RESULTS

(1) Inconsistent versus consistent breathing amplitudes caused a significant increase in deformation relative to the reference scan within the left lung (1.40 +/- 0.42 versus 1.29 +/- 0.36 mm, p < 0.05). (2) One-to-one correspondences between motions of internal anatomies and motion of the RPM external marker did not exist. The regression lines between the two types of motions did not yield an identity relationship (unity slope and zero intercept). Inconsistent breathing produced significantly different regression lines than consistent breathing in ten of the 12 studies within a left lung region of interest.

CONCLUSIONS

The results of these two studies indicate that inconsistency in the amplitude of breathing disrupted the correspondence between the motion of the external marker and internal anatomies. As a consequence, radiation therapy of tumors embedded in lung tissue may be prone to significant errors if inconsistent breathing amplitudes occur during treatment.

Evaluation of motion measurement using cine MRI for image guided stereotactic body radiotherapy on a new phantom platform

The objective of this study is to investigate accuracy of motion tracking of cine magnetic resonance imaging (MRI) for image-guided stereotactic body radiotherapy. A phantom platform was developed in this work to fulfill the goal. The motion phantom consisted of a platform, a solid thread, a motor and a control system that can simulate motion in various modes. To validate its reproducibility, the phantom platform was setup three times and imaged with fluoroscopy using an electronic portal imaging device (EPID) for the same motion profile. After the validation test, the phantom platform was evaluated using cine MRI at 2.5 frames/second on a 1.5T GE scanner using five different artificial profiles and five patient profiles. The above profiles were again measured with EPID fluoroscopy and used as references. Discrepancies between measured profiles from cine MRI and EPID were quantified using root-mean-square (RMS) and standard deviation (SD). Pearson's product moment correlational analysis was used to test correlation. The standard deviation for the reproducibility test was 0.28 mm. The discrepancies (RMS) between all profiles measured by cine MRI and EPID fluoroscopy ranged from 0.30 to 0.49 mm for artificial profiles and ranged from 0.75 to 0.91 mm for five patient profiles. The cine MRI sequence could precisely track phantom motion and the proposed motion phantom was feasible to evaluate cine MRI accuracy.

Stereotactic body radiation therapy: the report of AAPM Task Group 101

Task Group 101 of the AAPM has prepared this report for medical physicists, clinicians, and therapists in order to outline the best practice guidelines for the external-beam radiation therapy technique referred to as stereotactic body radiation therapy (SBRT). The task group report includes a review of the literature to identify reported clinical findings and expected outcomes for this treatment modality. Information is provided for establishing a SBRT program, including protocols, equipment, resources, and QA procedures. Additionally, suggestions for developing consistent documentation for prescribing, reporting, and recording SBRT treatment delivery is provided.

A novel analytical approach to the prediction of respiratory diaphragm motion based on external torso volume change

An analytical approach to predict respiratory diaphragm motion should have advantages over a correlation-based method, which cannot adapt to breathing pattern changes without re-calibration for a changing correlation and/or linear coefficient. To quantitatively calculate the diaphragm motion, a new expandable 'piston' respiratory (EPR) model was proposed and tested using 4DCT torso images of 14 patients. The EPR model allows two orthogonal lung motions (with a few volumetric constraints): (1) the lungs expand (DeltaV(EXP)) with the same anterior height variation as the thoracic surface, and (2) the lungs extend (DeltaV(EXT)) with the same inferior distance as the volumetrically equivalent 'piston' diaphragm. A volume conservation rule (VCR) established previously (Li et al 2009 Phys. Med. Biol. 54 1963-78) was applied to link the external torso volume change (TVC) to internal lung volume change (LVC) via lung air volume change (AVC). As the diaphragm moves inferiorly, the vacant space above the diaphragm inside the rib cage should be filled by lung tissue with a volume equal to DeltaV(EXT) (=LVC-DeltaV(EXP)), while the volume of non-lung tissues in the thoracic cavity should conserve. It was found that DeltaV(EXP) accounted for 3-24% of the LVC in these patients. The volumetric shape of the rib cage, characterized by the variation of cavity volume per slice over the piston motion range, deviated from a hollow cylinder by -1.1% to 6.0%, and correction was made iteratively if the variation is >3%. The predictions based on the LVC and TVC (with a conversion factor) were compared with measured diaphragm displacements (averaged from six pivot points), showing excellent agreements (0.2 +/- 0.7 mm and 0.2 +/- 1.2 mm, respectively), which are within clinically acceptable tolerance. Assuming motion synchronization between the piston and points of interest along the diaphragm, point motion was estimated but at higher uncertainty ( approximately 10% +/- 4%). This analytical approach provides a patient-independent technique to calculate the patient-specific diaphragm motion, using the anatomical and respiratory volumetric constraints.

Assessing four-dimensional radiotherapy planning and respiratory motion-induced dose difference based on biologically effective uniform dose

Four-dimensional (4D) radiotherapy is considered as a feasible and ideal solution to accommodate intra-fractional respiratory motion during conformal radiation therapy. With explicit inclusion of the temporal changes in anatomy during the imaging, planning, and delivery of radiotherapy, 4D treatment planning in principle provides better dose conformity. However, the clinical benefits of developing 4D treatment plans in terms of tumor control rate and normal tissue complication probability as compared to other treatment plans based on CT images of a fixed respiratory phase remains mostly unproven. The aim of our study is to comprehensively evaluate 4D treatment planning for nine lung tumor cases with both physical and biological measures using biologically effective uniform dose (D =) together with complication-free tumor control probability, P+. Based on the examined lung cancer patients and PTV margin applied, we found similar but not identical curves of DVH, and slightly different mean doses in tumor (up to 1.5%) and normal tissue in all cases when comparing 4D, P0%, and P50% plans. When it comes to biological evaluations, we did not observe definitively PTV size dependence in P+ among these nine lung cancer patients with various sizes of PTV. Moreover, it is not necessary that 4D plans would have better target coverage or higher P+ as compared to a fixed phase IMRT plan. However, on the contrary to significant deviations in P+ (up to 14.7%) observed if delivering the IMRT plan made at end-inhalation incorrectly at end-exhalation phase, we estimated the overall P+, PB, and PI for 4D composite plans that have accounted for intra-fractional respiratory motion.

Lung tumor tracking in fluoroscopic video based on optical flow

Respiratory gating and tumor tracking for dynamic multileaf collimator delivery require accurate and real-time localization of the lung tumor position during treatment. Deriving tumor position from external surrogates such as abdominal surface motion may have large uncertainties due to the intra- and interfraction variations of the correlation between the external surrogates and internal tumor motion. Implanted fiducial markers can be used to track tumors fluoroscopically in real time with sufficient accuracy. However, it may not be a practical procedure when implanting fiducials bronchoscopically. In this work, a method is presented to track the lung tumor mass or relevant anatomic features projected in fluoroscopic images without implanted fiducial markers based on an optical flow algorithm. The algorithm generates the centroid position of the tracked target and ignores shape changes of the tumor mass shadow. The tracking starts with a segmented tumor projection in an initial image frame. Then, the optical flow between this and all incoming frames acquired during treatment delivery is computed as initial estimations of tumor centroid displacements. The tumor contour in the initial frame is transferred to the incoming frames based on the average of the motion vectors, and its positions in the incoming frames are determined by fine-tuning the contour positions using a template matching algorithm with a small search range. The tracking results were validated by comparing with clinician determined contours on each frame. The position difference in 95% of the frames was found to be less than 1.4 pixels (approximately 0.7 mm) in the best case and 2.8 pixels (approximately 1.4 mm) in the worst case for the five patients studied.

Water-equivalent pathlength reproducibility due to respiratory pattern variation in charged-particle pancreatic radiotherapy

We evaluated the water-equivalent length (WEL) reproducibility due to variation in the external respiratory marker position when using a 4DCT scan in respiratory-gated charged-particle treatment. Two sets of pancreatic 4DCT data from two patients were acquired under free breathing conditions with 256-slice CT. The 4DCT data included two exhalation phases and the respiratory patterns in each patient differed, one being regular and the other irregular. The WEL calculation region is defined in the first respiratory cycle by two planes, one at the patient entrance surface and the other behind the target in the anterior-posterior (AP) and posterior-anterior (PA) directions. In the regular respiratory pattern, the WEL variation within the target region was less than 1.7 mm between the first and second exhalations in both AP and PA calculation directions. However, in the irregular breathing pattern, the respiratory amplitude at the second exhalation was 20% lower than that at the first exhalation; therefore, WEL variations from 8.1 to -9.1 mm and from 3.1 to -3.4 mm were observed within the target region in the AP and PA calculation directions, respectively. The WEL variation in the PA direction was smaller than that in the AP direction because the abdominal thickness is affected more in the AP direction. Respiratory pattern variation even affects WEL values in the respiratory-gated phase. This variation should be considered in treatment planning, and necessary improvements in respiratory reproducibility should be made.

Retrospective analysis of artifacts in four-dimensional CT images of 50 abdominal and thoracic radiotherapy patients

PURPOSE

To quantify the type, frequency, and magnitude of artifacts in four-dimensional (4D) CT images acquired using a multislice cine method.

METHODS AND MATERIALS

Fifty consecutive patients who underwent 4D-CT scanning and radiotherapy for thoracic or abdominal cancers were included in this study. All the 4D-CT scans were performed on the GE multislice PET/CT scanner with the Varian Real-time Position Management system in cine mode. The GE Advantage 4D software was used to create 4D-CT data sets. The artifacts were then visually and quantitatively analyzed. We performed statistical analyses to evaluate the relationships between patient- or breathing-pattern-related parameters and the occurrence as well as magnitude of artifacts.

RESULTS

It was found that 45 of 50 patients (90%) had at least one artifact (other than blurring) with a mean magnitude of 11.6 mm (range, 4.4-56.0 mm) in the diaphragm or heart. We also observed at least one artifact in 6 of 20 lung or mediastinal tumors (30%). Statistical analysis revealed that there were significant differences between several breathing-pattern-related parameters, including abdominal displacement (p < 0.01), for the subgroups of patients with and without artifacts. The magnitude of an artifact was found to be significantly but weakly correlated with the abdominal displacement difference between two adjacent couch positions (R = 0.34, p < 0.01).

CONCLUSIONS

This study has identified that the frequency and magnitude of artifacts in 4D-CT is alarmingly high. Significant improvement is needed in 4D-CT imaging.

Tumor trailing strategy for intensity-modulated radiation therapy of moving targets

Internal organ motion during the course of radiation therapy of cancer affects the distribution of the delivered dose and, generally, reduces its conformality to the targeted volume. Previously proposed approaches aimed at mitigating the effect of internal motion in intensity-modulated radiation therapy (IMRT) included expansion of the target margins, motion-correlated delivery (e.g., respiratory gating, tumor tracking), and adaptive treatment plan optimization employing a probabilistic description of motion. We describe and test the tumor trailing strategy, which utilizes the synergy of motion-adaptive treatment planning and delivery methods. We regard the (rigid) target motion as a superposition of a relatively fast cyclic component (e.g., respiratory) and slow aperiodic trends (e.g., the drift of exhalation baseline). In the trailing approach, these two components of motion are decoupled and dealt with separately. Real-time motion monitoring is employed to identify the "slow" shifts, which are then corrected by applying setup adjustments. The delivery does not track the target position exactly, but trails the systematic trend due to the delay between the time a shift occurs, is reliably detected, and, subsequently, corrected. The "fast" cyclic motion is accounted for with a robust motion-adaptive treatment planning, which allows for variability in motion parameters (e.g., mean and extrema of the tidal volume, variable period of respiration, and expiratory duration). Motion-surrogate data from gated IMRT treatments were used to provide probability distribution data for motion-adaptive planning and to test algorithms that identified systematic trends in the character of motion. Sample IMRT fields were delivered on a clinical linear accelerator to a programmable moving phantom. Dose measurements were performed with a commercial two-dimensional ion-chamber array. The results indicate that by reducing intrafractional motion variability, the trailing strategy enhances relevance and applicability of motion-adaptive planning methods, and improves conformality of the delivered dose to the target in the presence of irregular motion. Trailing strategy can be applied to respiratory-gated treatments, in which the correction for the slow motion can increase the duty cycle, while robust probabilistic planning can improve management of the residual motion within the gate window. Similarly, trailing may improve the dose conformality in treatment of patients who exhibit detectable target motion of low amplitude, which is considered insufficient to provide a clinical indication for the use of respiratory-gated treatment (e.g., peak-to-peak motion of less than 10 mm). The mechanical limitations of implementing tumor trailing are less rigorous than those of real-time tracking, and the same technology could be used for both.

The investigation on the location effect of external markers in respiratory-gated radiotherapy

PURPOSE

To investigate the effect of the marker placement on the correlation relationship between the motions of external markers and the internal target under different breathing patterns for several lung cancer patients.

METHOD AND MATERIAL

To monitor and record simultaneous motions of internal target and associated surrogate markers during respiratory gated radiotherapy, an infrared camera system synchronized with a medical simulator was installed in our institute. Multiple external markers were placed on the patients' chest wall with proper geometrical arrangement in closely monitoring the motion of skin near tumor. The motion signals of three breathing sessions (free breathing, breath-holding, and free breathing after breath-holding) were recorded and the quality of correlation between them was analyzed. For a single marker motion, its correlation with the internal target was analyzed using cross-covariance function. For the multiple markers, their correlation with the internal target was analyzed based on additive model.

RESULT

Seven patients undergoing radiotherapy with right upper or middle lobe lesions were enrolled in this study. Statistic analysis based on the internal-external motion signals shows that the effect of marker location on the quality of its correlation with the internal target is varied from patient to patient. There was no specific marker location where consistently demonstrated superior quality of correlation with the internal target motion over three breathing sessions for all patients. As the composite surrogate signal which was generated from the motions of multiple external markers was used to correlate the internal target motion, significant improvement of the quality of correlation was achieved.

CONCLUSION

The correlation of external marker to the internal target could be influenced by several factors such as patient population, marker locations, and breathing patterns, considerably. The quality of correlation and predictability to the internal target furnished by a single external marker is inferior to that of the composite signal generated from multiple external markers. The use of composite signal shows great potential in improving the predictability of internal target motion and presents an effective way to track tumor more accurately.

Effect of novel amplitude/phase binning algorithm on commercial four-dimensional computed tomography quality

PURPOSE

Respiratory motion is a significant source of anatomic uncertainty in radiotherapy planning and can result in errors of portal size and the subsequent radiation dose. Although four-dimensional computed tomography allows for more accurate analysis of the respiratory cycle, breathing irregularities during data acquisition can cause considerable image distortions. The aim of this study was to examine the effect of respiratory irregularities on four-dimensional computed tomography, and to evaluate a novel image reconstruction algorithm using percentile-based tagging of the respiratory cycle.

METHODS AND MATERIALS

Respiratory-correlated helical computed tomography scans were acquired for 11 consecutive patients. The inspiration and expiration data sets were reconstructed using the default phase-based method, as well as a novel respiration percentile-based method with patient-specific metrics to define the ranges of the reconstruction. The image output was analyzed in a blinded fashion for the phase- and percentile-based reconstructions to determine the prevalence and severity of the image artifacts.

RESULTS

The percentile-based algorithm resulted in a significant reduction in artifact severity compared with the phase-based algorithm, although the overall artifact prevalence did not differ between the two algorithms. The magnitude of differences in respiratory tag placement between the phase- and percentile-based algorithms correlated with the presence of image artifacts.

CONCLUSION

The results of our study have indicated that our novel four-dimensional computed tomography reconstruction method could be useful in detecting clinically relevant image distortions that might otherwise go unnoticed and to reduce the image distortion associated with some respiratory irregularities. Additional work is necessary to assess the clinical impact on areas of possible irregular breathing.

Verifying 4D gated radiotherapy using time-integrated electronic portal imaging: a phantom and clinical study

Background

Respiration-gated radiotherapy (RGRT) can decrease treatment toxicity by allowing for smaller treatment volumes for mobile tumors. RGRT is commonly performed using external surrogates of tumor motion. We describe the use of time-integrated electronic portal imaging (TI-EPI) to verify the position of internal structures during RGRT delivery

Methods

TI-EPI portals were generated by continuously collecting exit dose data (aSi500 EPID, Portal vision, Varian Medical Systems) when a respiratory motion phantom was irradiated during expiration, inspiration and free breathing phases. RGRT was delivered using the Varian RPM system, and grey value profile plots over a fixed trajectory were used to study object positions. Time-related positional information was derived by subtracting grey values from TI-EPI portals sharing the pixel matrix. TI-EPI portals were also collected in 2 patients undergoing RPM-triggered RGRT for a lung and hepatic tumor (with fiducial markers), and corresponding planning 4-dimensional CT (4DCT) scans were analyzed for motion amplitude.

Results

Integral grey values of phantom TI-EPI portals correlated well with mean object position in all respiratory phases. Cranio-caudal motion of internal structures ranged from 17.5–20.0 mm on planning 4DCT scans. TI-EPI of bronchial images reproduced with a mean value of 5.3 mm (1 SD 3.0 mm) located cranial to planned position. Mean hepatic fiducial markers reproduced with 3.2 mm (SD 2.2 mm) caudal to planned position. After bony alignment to exclude set-up errors, mean displacement in the two structures was 2.8 mm and 1.4 mm, respectively, and corresponding reproducibility in anatomy improved to 1.6 mm (1 SD).

Conclusion

TI-EPI appears to be a promising method for verifying delivery of RGRT. The RPM system was a good indirect surrogate of internal anatomy, but use of TI-EPI allowed for a direct link between anatomy and breathing patterns.

The correlation evaluation of a tumor tracking system using multiple external markers

The purpose of this study is to evaluate the correlations between external markers and internal targets for radiation therapy of lung cancer patients. Using an infrared camera system coupled with a clinical simulator, the simultaneous motions of multiple external markers and an internal target were obtained. The correlation between external and internal signals was analyzed using a cross-covariance function. A linear regression model was employed to generate a composite signal from multiple external markers in order to predict the internal target motion. The external and internal signals, and their correlations, demonstrated a wide range of variation with respect to marker location, motion dimension, and breathing pattern. The performance of the composite signal indicates that when more external signals were taken into account, the mean correlation between the composite signal and internal signal was improved. This implies that a combination of multiple external signals might be an improved way to predict internal target motion. Also, since the characteristics of respiratory signals can vary significantly, certain methods of preprocessing and external signal combination are necessary.

The management of respiratory motion in radiation oncology report of AAPM Task Group 76

This document is the report of a task group of the AAPM and has been prepared primarily to advise medical physicists involved in the external-beam radiation therapy of patients with thoracic, abdominal, and pelvic tumors affected by respiratory motion. This report describes the magnitude of respiratory motion, discusses radiotherapy specific problems caused by respiratory motion, explains techniques that explicitly manage respiratory motion during radiotherapy and gives recommendations in the application of these techniques for patient care, including quality assurance (QA) guidelines for these devices and their use with conformal and intensity modulated radiotherapy. The technologies covered by this report are motion-encompassing methods, respiratory gated techniques, breath-hold techniques, forced shallow-breathing methods, and respiration-synchronized techniques. The main outcome of this report is a clinical process guide for managing respiratory motion. Included in this guide is the recommendation that tumor motion should be measured (when possible) for each patient for whom respiratory motion is a concern. If target motion is greater than 5 mm, a method of respiratory motion management is available, and if the patient can tolerate the procedure, respiratory motion management technology is appropriate. Respiratory motion management is also appropriate when the procedure will increase normal tissue sparing. Respiratory motion management involves further resources, education and the development of and adherence to QA procedures.

Initial evaluation of a four-dimensional computed tomography system using a programmable motor

A dynamic lung tumor phantom was used to investigate the geometric reconstruction accuracy of a commercial four-dimensional computed tomography (4D-CT) system. A ball filled with resin, embedded in a cork cube, was placed on a moving platform. Various realistic antero-posterior (AP) motions were programmed to reproduce the respiratory motion of a lung tumor. Several three-dimensional (3D) CT and 4D-CT images of this moving object were acquired and compared using different acquisition parameters. Apparent volume and diameter of the ball were measured and compared to the real values. The position of two points (the AP limits of the ball) during the motion in the coordinate system of the CT scanner were also compared with the expected values. Volume error was shown to increase with object speed. However, although the volume error was associated with intraslice artifacts, it did not reflect large interslice inconstancies in object position and should not be used as an indicator of image accuracy. The 3D-CT gave a random position of the tumor along the phantom excursion; accuracy in the assessment of position by 4D-CT ranged from 0.4 mm to 2.6 mm during extreme phases of breathing. We used average projection (AVE) and maximum intensity projection (MIP) algorithms available on the commercial software to create internal target volumes (ITVs) by merging gross tumor volume (GTV) images at various respiratory phases. The ITVs were compared to a theoretical value computed from the programmed ball excursion. The ITVs created from the MIP algorithm were closer to the theoretical value (within 12%) than were those created from the AVE algorithm (within 40%).

Performance and characteristics of an IR localizing system for radiation therapy

We report the development of a new system for interactive patient posture, position and respiratory control during radiation therapy treatment. The system consists of an infrared camera, retro-reflective markers and dedicated software that makes it practical to use in the clinic. The system is designed to be used with multiple retro-reflective markers to monitor not only position, but also the posture of the patient in real time. Specific features of the system include: 1. The system reports an absolute misalignment at several points on a patients, and also provides feedbacks on any necessary adjustments in terms of site specific set-up parameters, such as focus to surface distance (PIN), superior and inferior alignment, chest-wall angle, etc. 2. The system is based on the set of predefined templates containing number and position of control points and feedback parameters developed for different treatment sites. 3. A "virtual portal vision" procedure is developed to project organ contours in the beams-eye-view (BEV) based on the marker locations obtained in real time and compare them with digitally reconstructed radiographs (DRRs) from CT simulation. Assuming good correlation between external markers and internal anatomy, the system offers the possibility of mimicking a verification procedure without taking port-films, which can potentially reduce the setup time. In this paper, we concentrate on system properties and performance, while initial applications on a number of clinical sites is ongoing. Accuracy and precision of this system are evaluated in the context of breast/chest treatments using rigid phantoms. The system has an intrinsic uncertainty of +/- 1 mm; and when two systems in different rooms (CT and treatment room) are used for correlating positional information, the uncertainty is less than 2 mm.

Dosimetric characteristics of a new linear accelerator under gated operation

Respiratory gated radiotherapy may allow reduction of the treatment margins, thus sparing healthy tissue and/or allowing dose escalation to the tumor. However, current commissioning and quality assurance of linear accelerators do not include evaluation of gated delivery. The purpose of this study is to test gated photon delivery of a Siemens ONCOR Avant‐Garde linear accelerator. Dosimetric characteristics for gated and nongated delivery of 6‐MV and 15‐MV photons were compared for the range of doses, dose rates, and for several gating regimes. Dose profiles were also compared using Kodak EDR2 and X‐Omat V films for 6‐MV and 15‐MV photons for several dose rates and gating regimes. Results showed that deviation is less than or equal to 0.6% for all dose levels evaluated with the exception of the lowest dose delivered at 25 MU at an unrealistically high gating frequency of 0.5 Hz. At 400 MU, dose profile deviations along the central axes in in‐plane and cross‐plane directions within 80% of the field size are below 0.7%. No unequivocally detectable dose profile deviation was observed for 50 MU. Based on the comparison with widely accepted standards for conventional delivery, our results indicate that this LINAC is well suited for gated delivery of nondynamic fields.

PACS numbers: 87.56‐By, 87.66‐Cd, 87.66‐Jj

Interfractional anatomic variation in patients treated with respiration-gated radiotherapy

As quality assurance for respiration‐gated treatments using the Varian RPM™ system, we monitor interfractional diaphragm variation throughout treatment using extra anterior‐posterior (AP) portal images. We measure the superior‐inferior (SI) distance between one or more bony landmarks and the ipsilateral diaphragm dome in each such radiograph and calculate its difference, D, from the corresponding distance in a planning CT scan digitally reconstructed radiograph (DRR). For each patient, the mean of D represents the systematic diaphragm displacement, and the standard deviation of D represents random diaphragm variations and is a measure of interfractional gating reproducibility. We present results for 31 sequential patients (21 lung, 10 liver tumors), each with at least 8 such portal images. For all patients, the gate included end‐exhale. The patient‐specific duty cycle ranged from 30% to 60%. All patients received customized audio prompting for simulation and treatment, and 14 patients also received visual prompting. Respiration‐synchronized fluoroscopic movies taken at a conventional simulator revealed patient‐specific diaphragm excursions from 1.0 cm to 5.0 cm and diaphragm excursion within the gate from 0.5 cm to 1.0 cm, demonstrating a significant reduction of intra‐fractional diaphragm (and by inference tumor) motion by respiratory gating. One standard deviation of the systematic displacement (the mean of D) was 0.63 cm and 0.48 cm for the lung and liver patient groups, respectively. The average ±1 SD of the random displacements (i.e., the average of the standard deviations of D) was 0.42±0.11 cm and 0.50±0.19 for the two groups, respectively. The similar magnitude of the systematic and random displacements suggests that both derive from a common distribution of interfractional variations. Combining visual with audio prompting did not significantly improve performance, as judged by D. Guided by these portal images, field changes were made during the course of treatment for 6 patients (1 lung, 5 liver).

PACS numbers: 87.53.‐j, 87.53.Oq

Comparison of helical and cine acquisitions for 4D-CT imaging with multislice CT

We proposed a data sufficiency condition (DSC) for four-dimensional-CT (4D-CT) imaging on a multislice CT scanner, designed a pitch factor for a helical 4D-CT, and compared the acquisition time, slice sensitivity profile (SSP), effective dose, ability to cope with an irregular breathing cycle, and gating technique (retrospective or prospective) of the helical 4D-CT and the cine 4D-CT on the General Electric (GE) LightSpeed RT (4-slice), Plus (4-slice), Ultra (8-slice) and 16 (16-slice) multislice CT scanners. To satisfy the DSC, a helical or cine 4D-CT acquisition has to collect data at each location for the duration of a breathing cycle plus the duration of data acquisition for an image reconstruction. The conditions for the comparison were 20 cm coverage in the cranial-caudal direction, a 4 s breathing cycle, and half-scan reconstruction. We found that the helical 4D-CT has the advantage of a shorter scan time that is 10% shorter than that of the cine 4D-CT, and the disadvantages of 1.8 times broadening of SSP and requires an additional breathing cycle of scanning to ensure an adequate sampling at the start and end locations. The cine 4D-CT has the advantages of maintaining the same SSP as slice collimation (e.g., 8 x 2.5 mm slice collimation generates 2.5 mm SSP in the cine 4D-CT as opposed to 4.5 mm in the helical 4D-CT) and a lower dose by 4% on the 8- and 16-slice systems, and 8% on the 4-slice system. The advantage of faster scanning in the helical 4D-CT will diminish if a repeat scan at the location of a breathing irregularity becomes necessary. The cine 4D-CT performs better than the helical 4D-CT in the repeat scan because it can scan faster and is more dose efficient.

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.

Application of the spirometer in respiratory gated radiotherapy

The signal from a spirometer is directly correlated with respiratory motion and is ideal for target respiratory motion tracking. However, its susceptibility to signal drift deters its application in radiotherapy. In this work, a few approaches are investigated to control spirometer signal drift for a Bernoulli-type spirometer. A method is presented for rapid daily calibration of the spirometer to obtain a flow sensitivity function. Daily calibration assures accurate airflow measurement and also reduces signal drift. Dynamic baseline adjustment further controls the signal drift. The accuracy of these techniques was studied and it was found that the spirometer is able to provide a long-term drift-free breathing signal. The tracking error is comprised of two components: calibration error and stochastic signal baseline variation error. The calibration error is very small (1% of 3 l) and therefore negligible. The stochastic baseline variation error can be as large as 20% of the normal breathing amplitude. In view of these uncertainties, the applications of spirometers in treatment techniques that rely on breathing monitoring are discussed. Spirometer-based monitoring is noted most suitable for deep inspiration breath-hold but less important for free breathing gating techniques.

New radiotherapy technologies

Conventional radiation therapy has had limited success in curing inoperable lung cancer due to poor local control. There is evidence to suggest that higher doses of radiation will improve local control. In order to safely deliver higher doses of thoracic radiation, advanced treatment techniques are required. Different biologic indices have been utilized to determine whether dose escalation can be safely accomplished, and the results have been reported from many institutions. Tumor motion control aids in treatment since it allows radiation oncologists to more accurately target tumors and therefore to spare more normal tissue from the radiation field. The imaging information from 18-FDG-PET scans also improves target delineation. Advanced treatment delivery techniques, such as three-dimensional conformal radiation therapy, intensity modulated radiation therapy, and stereotactic radiosurgery are also being used to safely escalate the radiation dose. This article explores the current literature on these issues and other advanced radiation therapy techniques.