Analysis of movement of intrathoracic neoplasms using ultrafast computerized tomography

Respiration-controlled radiotherapy in lung cancer: Systematic evaluation of the optimal application practice

Background and purpose

Definitive radiochemotherapy (RCT) for non-small cell lung cancer (NSCLC) in UICC/TNM I-IVA (singular, oligometastatic) is one of the treatment methods with a potentially curative concept. However, tumour respiratory motion during RT requires exact pre-planning. There are various techniques of motion management like creating internal target volume (ITV), gating, inspiration breath-hold and tracking. The primary goal is to cover the PTV with the prescribed dose while at the same time maximizing dose reduction of surrounding normal tissues (organs at risk, OAR). In this study, two standardized online breath-controlled application techniques used alternately in our department are compared with respect to lung and heart dose.

Materials and methods

Twenty-four patients who were indicated for thoracic RT received planning CTs in voluntary deep inspiration breath-hold (DIBH) and in free shallow breathing, prospectively gated in expiration (FB-EH). A respiratory gating system by Varian (Real-time Position Management, RPM) was used for monitoring. OAR, GTV, CTV and PTV were contoured on both planning CTs. The PTV margin to the CTV was 5 mm in the axial and 6-8 mm in the cranio-caudal direction. The consistency of the contours was checked by elastic deformation (Varian Eclipse Version 15.5). RT plans were generated and compared in both breathing positions using the same technique, IMRT over fixed irradiation directions or VMAT. The patients were treated in a prospective registry study with the approval of the local ethics committee.

Results

The PTV in expiration (FB-EH) was on average significantly smaller than the PTV in inspiration (DIBH): for tumours in the lower lobe (LL) 431.5 vs. 477.6 ml (Wilcoxon test for connected samples; p = 0.004), in the upper lobe (UL) 659.5 vs. 686.8 ml (p = 0.005). The intra-patient comparison of plans in DIBH and FB-EH showed superiority of DIBH for UL-tumours and equality of DIBH and FB-EH for LL-tumours. The dose for OAR in UL-tumours was lower in DIBH than in FB-EH (mean lung dose p = 0.011; lungV20, p = 0.002; mean heart dose p = 0.016). The plans for LL-tumours in FB-EH showed no difference in OAR compared to DIBH (mean lung dose p = 0.683; V20Gy p = 0.33; mean heart dose p = 0.929). The RT setting was controlled online for each fraction and was robustly reproducible in FB-EH.

Conclusion

RT plans for treating lung tumours implemented depend on the reproducibility of the DIBH and advantages of the respiratory situation with respect to OAR. The primary tumour localization in UL correlates with advantages of RT in DIBH, compared to FB-EH. For LL-tumours there is no difference between RT in FB-EH and RT in DIBH with respect to heart or lung exposure and therefore, reproducibility is the dominant criterion. FB-EH is recommended as a very robust and efficient technique for LL-tumours.

Machine-learning-based prediction of the effectiveness of the delivered dose by exhale-gated radiotherapy for locally advanced lung cancer: The additional value of geometric over dosimetric parameters alone

Purpose

This study aimed to assess interfraction stability of the delivered dose distribution by exhale-gated volumetric modulated arc therapy (VMAT) or intensity-modulated arc therapy (IMAT) for lung cancer and to determine dominant prognostic dosimetric and geometric factors.

Methods

Clinical target volume (CTV_Plan) from the planning CT was deformed to the exhale-gated daily CBCT scans to determine CTV_i, treated by the respective dose fraction. The equivalent uniform dose of the CTV_i was determined by the power law (gEUD_i) and cell survival model (EUD_iSF) as effectiveness measure for the delivered dose distribution. The following prognostic factors were analyzed: (I) minimum dose within the CTV_i (D_{min_i}), (II) Hausdorff distance (HDD_i) between CTV_i and CTV_Plan, (III) doses and deformations at the point in CTV_Plan at which the global minimum dose over all fractions per patient occurs (PD_{min_global_i}), and (IV) deformations at the point over all CTV_i margins per patient with the largest Hausdorff distance (HDPw_orst). Prognostic value and generalizability of the prognostic factors were examined using cross-validated random forest or multilayer perceptron neural network (MLP) classifiers. Dose accumulation was performed using back deformation of the dose distribution from CTV_i to CTV_Plan.

Results

Altogether, 218 dose fractions (10 patients) were evaluated. There was a significant interpatient heterogeneity between the distributions of the normalized gEUD_i values (p<0.0001, Kruskal-Wallis tests). Accumulated gEUD over all fractions per patient was 1.004-1.023 times of the prescribed dose. Accumulation led to tolerance of ~20% of fractions with gEUD_i <93% of the prescribed dose. Normalized D_min >60% was associated with predicted gEUD values above 95%. D_min had the highest importance for predicting the gEUD over all analyzed prognostic parameters by out-of-bag loss reduction using the random forest procedure. Cross-validated random forest classifier based on D_min as the sole input had the largest Pearson correlation coefficient (R=0.897) in comparison to classifiers using additional input variables. The neural network performed better than the random forest classifier, and the gEUD values predicted by the MLP classifier with D_min as the sole input were correlated with the gEUD values characterized by R=0.933 (95% CI, 0.913-0.948). The performance of the full MLP model with all geometric input parameters was slightly better (R=0.952) than that based on D_min (p=0.0034, Z-test).

Conclusion

Accumulated dose distributions over the treatment series were robust against interfraction CTV deformations using exhale gating and online image guidance. D_min was the most important parameter for gEUD prediction for a single fraction. All other parameters did not lead to a markedly improved generalizable prediction. Dosimetric information, especially location and value of D_min within the CTV _i , are vital information for image-guided radiation treatment.

Adaptive radiation therapy: When, how and what are the benefits that literature provides?

PURPOSE

To identify from the current literature when is the right time to replan and to assign thresholds for the optimum process of replanning. Nowadays, adaptive radiotherapy (ART) for head and neck cancer plays an exceptional role consisting of an evaluation procedure of the prominent anatomical and dosimetric variations. By performing complex radiotherapy methods, the credibility of the therapeutic result is crucial. Image guided radiotherapy (IGRT) was developed to ensure locoregional control and thus changes that might occur during radiotherapy be dealt with.

MATERIALS AND METHODS

An electronic research of articles published in PubMed/MEDLINE and Science Direct databases from January 2004 to October 2020 was performed. Among a total of 127 studies assessed for eligibility, 85 articles were ultimately retained for the review.

RESULTS

The most noticeable changes have been reported in the middle fraction of the treatment. Therefore, the suggested optimal time to replan is between the third and the fourth week. Anatomical deviations>1cm in the external contour, average weight loss>10%, violation in the dose coverage of the targets>5%, and violation in the dose of the peripherals were some of the thresholds that are currently used, and which lead to replanning.

CONCLUSION

ART may decrease toxicity and improve local-control. Whether it is beneficial or not, depends ultimately on each patient. However, more investigation of the changes should be performed in future prospective studies to obtain more accurate results.

3D dosimetric validation of ultrasound-guided radiotherapy with a dynamically deformable abdominal phantom

OBJECTIVES

The purpose of this study was to dosimetrically benchmark gel dosimetry measurements in a dynamically deformable abdominal phantom for intrafraction image guidance through a multi-dosimeter comparison. Once benchmarked, the study aimed to perform a proof-of-principle study for validation measurements of an ultrasound image-guided radiotherapy delivery system.

METHODS

The phantom was dosimetrically benchmarked by delivering a liver VMAT plan and measuring the 3D dose distribution with DEFGEL dosimeters. Measured doses were compared to the treatment planning system and measurements acquired with radiochromic film and an ion chamber. The ultrasound image guidance validation was performed for a hands-free ultrasound transducer for the tracking of liver motion during treatment.

RESULTS

Gel dosimeters were compared to the TPS and film measurements, showing good qualitative dose distribution matches, low γ values through most of the high dose region, and average 3%/5 mm γ-analysis pass rates of 99.2%(0.8%) and 90.1%(0.8%), respectively. Gel dosimeter measurements matched ion chamber measurements within 3%. The image guidance validation study showed the measurement of the treatment delivery improvements due to the inclusion of the ultrasound image guidance system. Good qualitative matching of dose distributions and improvements of the γ-analysis results were observed for the ultrasound-gated dosimeter compared to the ungated dosimeter.

CONCLUSIONS

DEFGEL dosimeters in phantom showed good agreement with the planned dose and other dosimeters for dosimetric benchmarking. Ultrasound image guidance validation measurements showed good proof-of-principle of the utility of the phantom system as a method of validating ultrasound-based image guidance systems and potentially other image guidance methods.

Real-time measurement of ICD lead motion during stereotactic body radiotherapy of ventricular tachycardia

Background

Here we aimed to evaluate the respiratory and cardiac-induced motion of a ICD lead used as surrogate in the heart during stereotactic body radiotherapy (SBRT) of ventricular tachycardia (VT). Data provides insight regarding motion and motion variations during treatment.

Materials and methods

We analyzed the log files of surrogate motion during SBRT of ventricular tachycardia performed in 20 patients. Evaluated parameters included the ICD lead motion amplitudes; intrafraction amplitude variability; correlation error between the ICD lead and external markers; and margin expansion in the superior-inferior (SI), latero-lateral (LL), and anterior-posterior (AP) directions to cover 90% or 95% of all amplitudes.

Results

In the SI, LL, and AP directions, respectively, the mean motion amplitudes were 5.0 ± 2.6, 3.4. ± 1.9, and 3.1 ± 1.6 mm. The mean intrafraction amplitude variability was 2.6 ± 0.9, 1.9 ± 1.3, and 1.6 ± 0.8 mm in the SI, LL, and AP directions, respectively. The margins required to cover 95% of ICD lead motion amplitudes were 9.5, 6.7, and 5.5 mm in the SI, LL, and AP directions, respectively. The mean correlation error was 2.2 ± 0.9 mm.

Conclusions

Data from online tracking indicated motion irregularities and correlation errors, necessitating an increased CTV-PTV margin of 3 mm. In 35% of cases, the motion variability exceeded 3 mm in one or more directions. We recommend verifying the correlation between CTV and surrogate individually for every patient, especially for targets with posterobasal localization where we observed the highest difference between the lead and CTV motion.

Determination of reproducibility of end-exhaled breath-holding in stereotactic body radiation therapy

Methods to evaluate the positional reproducibility of breath-hold irradiation mostly require manual operation. The purpose of this study is to propose a method to determine the reproducibility of breath-hold irradiation of lung tumors between fractions using non-artificial methods. This study included 13 patients who underwent terminal exhaled breath-hold irradiation for primary and metastatic lung cancer. All subjects received a prescribed dose of 60 Gy/8 fractions. The contours of the gross tumor volume (GTV) were extracted by threshold processing using treatment-planning computed tomography (CT) and cone-beam CT (CBCT), which was done just before the beginning of the treatment. The method proposed in this study evaluates the dice similarity coefficient (DSC) and Hausdorff distance (HD) by comparing two volumes, the GTVCTS (GTV obtained from treatment-planning CT) and GTVCBCT (GTV obtained from CBCT). The reference contours for DSC and HD are represented by GTVCTS. The results demonstrated good visual agreement for cases with a DSC of ~0.7. However, apparent misalignment occurred when the DSC was <0.5. HD was >2 mm in 3 out of 13 cases, and when the DSC was ~0.7, the HD was ~1 mm. In addition, cases with greater HD also demonstrated more significant variability. It was found that the DSC and HD evaluation methods for the positional reproducibility of breath-hold irradiation proposed in this study are straightforward and can be performed without the involvement of humans. Our study is of extreme significance in the field of radiation studies.

New method for measurement of chest surface motion in lung cancer patients: Quantification using a technique of deformable image registration

The purpose of this study was to measure the motion of the chest surface during breath-holding treatment for lung cancer using deformable image registration (DIR). Forty non-small-cell lung cancer patients treated with breath-holding stereotactic body radiation therapy were retrospectively examined. First, intensity-based DIR between 2 breath-holding computed tomography (CT) images was performed. Subsequently, deformation vector field (DVF) for all dimensions (left-right, anterior-posterior, and superior-inferior) was calculated from the result. For the analysis of chest surface, the DVF value of the only chest surface area was extracted after the chest surface was divided into 12 regions of interest (ROI) based on anatomy. Additionally, for the analysis of the correlation with the internal tumor motion, the median value of DVF for each surface ROI and the motion of the center of gravity of the tumor volume were used. It was possible to calculate the motion of chest surface without any outliers for all patients. For the average of 12 surface ROIs, the motion of 3D chest surface was within 2 mm (30 cases), 3 mm (8 cases), and 4 mm (2 cases). There was no correlation between the motion of the chest surface and that of the tumor for all 12 surface ROIs. We proposed a technique to evaluate the surface motion using DIR between multiple CT images. It could be a useful tool to calculate the motion of chest surface.

A study of the interplay effect for VMAT SBRT using a four-axes motion phantom

Purpose

To assess the accuracy of volumetric modulated arc therapy (VMAT) stereotactic body radiation therapy (SBRT) when treating moving targets (such as lung or liver lesions), focusing on the impact of the interplay effect in the event of complex breathing motion and when a gating window is used.

Methods

A dedicated programmable motion platform was implemented. This platform can carry large quality assurance (QA) phantoms and achieve complex three‐dimensional (3D) motion. Volumetric modulated arc therapy SBRT plans were delivered with TrueBeam linac to this moving setup and the measured dose was compared to the computed one. Several parameters were assessed such as breathing period, dose rate, dose prescription, shape of the breathing pattern, the use of a planning target volume (PTV) margin, and the use of a gating window.

Results

Loss of dose coverage (D95%) was acceptable in most situations. The doses received by 95% of the CTV, D95% ( CTVm) ranged from 94 to 101% (mean 98%) and the doses received by 2% of the CTV D2% ( CTVm) ranged from 94% to 110% of the prescribed dose. A visible interplay effect was observed when no margin was used or when the number of breathing cycles during the treatment delivery was lower than 20.

Conclusions

In our clinical context, treating lung and liver lesions using VMAT SBRT is reasonable. The interplay effect was moderated and acceptable in all simulated situations.

Correlating lung tumour location and motion with respiration using 4D CT scans

Background:

Lung tumours, especially those in the lower lobes, can move a lot during respiration; this motion needs to be accounted for during radiotherapy. In cases where 4D CT simulation scans are not performed, the current protocol at our centre is to apply a generic (internal motion + setup) margin of 0·70 cm in the axial plane and 1·20 cm in the longitudinal plane to all lung tumours, regardless of location. We analyse the tumour motions of a cohort of our local patients and categorise them into different locations in the lung. We seek to assess the adequacy of the current margins and to derive a more accurate set of standard margins which are specific for lung tumour locations.

Methods:

All cases of lung tumours treated with stereotactic ablative radiotherapy between 2012 and 2016 were identified retrospectively and 4D CT scan data analysed. These tumours were grouped into the following locations: upper zone (UZ), middle zone (MZ) and lower zone (LZ). The treatment planning system was used to generate the displacements of the centre of mass of the tumours in the right–left, anterior–posterior and superior–inferior axes; these were compared with the current generic margins. Median displacements were calculated for each axis in each location. New planning target volume (PTV) margins were derived by summing the median displacement, median absolute deviation (MAD) and 0·5 cm (for setup error).

Results:

Sixty-three cases were eligible for analyses. Motion in the superior–inferior direction was the greatest for all tumour locations, ranging from a median of 0·17 cm (MAD 0·12 cm) in UZ to 0·77 cm (MAD 0·27 cm) in LZ. Median tumour displacements in the anterior–posterior and right–left axes were similar for all locations, <0·30 and 0·20 cm, respectively. The current generic margins were adequate for only one-third of the cases in this study. A new PTV margin of 2·10 cm in the superior–inferior axis may be required for LZ tumours, while an additional 1–2 mm should be added to the current radial margins.

Conclusion:

The current generic margins are inadequate for the majority of cases. Tumour motion is the greatest in LZ in the superior–inferior axis. Motion mitigation strategies are essential for large LZ tumours.

Performance Validation of In-House Developed Four-dimensional Dynamic Phantom

Objective:

The objective of this study was to validate the performance characteristics of in-house developed four-dimensional (4D) dynamic phantom (FDDP).

Materials and Methods:

There are three target inserts of 1.0, 1.5 and 2.0 cm diameter. The targets were driven in sinusoidal pattern in the longitudinal direction, using the combinations of amplitudes of 0.5, 1.0, and 1.5 cm with frequencies of 0.2 and 0.25 Hz. The amplitude and frequency of motion were measured manually, and by using Real-Time Position Management (RPM) system also. The static, free-breathing, and 4D computed tomography (CT) scans of the phantom were acquired with 1.0 mm slice thickness. The 4DCT scans were sorted into 0%–90% phase, and the maximum intensity projection (MIP) images were also generated. The static, free-breathing, and 4DCT data sets and MIP images were contoured to get V_Static, V_FB, V₀₀......V₉₀, and internal target volume ITV _MIP, respectively. The individual phase volumes were summed to obtain V_4D. The length of the target in the motion was measured using MIP image and compared with theoretical length (TL). The variation of 3D displacement vector of individual phase volume with respect to V₀₀ with the phase of motion was studied at amplitude and frequency of 1.0 cm and 0.25 Hz, respectively. The degree of similarity between V_FB and V_4D and V_FB and ITV_MIP was also studied for all the target sizes at amplitude and frequency of 1.0 cm and 0.2 Hz and 1.0 cm and 0.25 Hz, respectively.

Results:

The amplitude and frequency of motion agreed within the limits of uncertainty with the manually and RPM measured values. The length of target in the motion matched within 1.0 mm with TL. The 3D displacement of individual phase volume showed no target size dependence, and the degree of similarity between V_FB and V_4D and V_FB and ITV_MIP decreases with increase in the displacement between the two volumes.

Conclusions:

The mechanical and imaging performances of FDDP were found within the acceptable limits. Therefore, this phantom can be used for quality assurance of 4D imaging process in radiotherapy.

Optimization based trajectory planning for real-time 6DoF robotic patient motion compensation systems

Purpose

Robotic stabilization of a therapeutic radiation beam with respect to a dynamically moving tumor target can be accomplished either by moving the radiation source, the patient, or both. As the treatment beam is on during this process, the primary goal is to minimize exposure of normal tissue to radiation as much as possible when moving the target back to the desired position. Due to the complex mechanical structure of 6 degree-of-freedom (6DoF) robots, it is not intuitive as to what 6 dimensional (6D) correction trajectory is optimal in achieving such a goal. With proportional-integrative-derivative (PID) and other controls, the potential exists that the controller may generate a trajectory that is highly curved, slow, or suboptimal in that it leads to unnecessary exposure of healthy tissue to radiation. This work investigates a novel feedback planning method that takes into account a robot’s mechanical joint structure, patient safety tolerances, and other system constraints, and performs real-time optimization to search the entire 6D trajectory space in each time cycle so it can respond with an optimal 6D correction trajectory.

Methods

Computer simulations were created for two 6DoF robotic patient support systems: a Stewart-Gough platform for moving a patient’s head in frameless maskless stereotactic radiosurgery, and a linear accelerator treatment table for moving a patient in prostate cancer radiation therapy. Motion planning was formulated as an optimization problem and solved at real-time speeds using the L-BFGS algorithm. Three planning methods were investigated, moving the platform as fast as possible (platform-D), moving the target along a straight-line (target-S), and moving the target based on the fastest descent of position error (target-D). Both synthetic motion and prior recorded human motion were used as input data and output results were analyzed.

Results

For randomly generated 6D step-like and sinusoidal synthetic input motion, target-D planning demonstrated the smallest net trajectory error in all cases. On average, optimal planning was found to have a 45% smaller target trajectory error than platform-D control, and a 44% smaller target trajectory error than target-S planning. For patient head motion compensation, only target-D planning was able to maintain a ≤0.5mm and ≤0.5deg clinical tolerance objective for 100% of the treatment time. For prostate motion, both target-S planning and target-D planning outperformed platform-D control.

Conclusions

A general 6D target trajectory optimization framework for robotic patient motion compensation systems was investigated. The method was found to be flexible as it allows control over various performance requirements such as mechanical limits, velocities, acceleration, or other system control objectives.

Delineation of a Cardiac Planning Organ-At-Risk Volume Using Real-Time Magnetic Resonance Imaging for Cardiac Protection in Thoracic and Breast Radiation Therapy

PURPOSE

Cardiac radiation is associated with cardiotoxicity in patients with thoracic and breast malignancies. We conducted a prospective study using cine magnetic resonance imaging (MRI) scans to evaluate heart motion. We hypothesized that cine MRI could be used to define population-based cardiac planning organ-at-risk volumes (PRV).

METHODS AND MATERIALS

A total of 16 real-time acquisitions were obtained per subject on a 1.5 Tesla MRI (Philips Ingenia). Planar cine MRI was performed in 4 sequential sagittal and coronal planes at free-breathing (FB) and deep-inspiratory breath hold (DIBH). In-plane cardiac motion was assessed using a scale-invariant feature transformation-based algorithm. Subject-specific pixel motion ranges were defined in anteroposterior (AP), left-right (LR), and superoinferior (SI) planes. Averages of the 98% and 67% of the maximum ranges of pixel displacement were defined by subject, then averaged across the cohort to calculate PRV expansions at FB and DIBH.

RESULTS

Data from 20 subjects with a total of 3120 image frames collected per subject in coronal and sagittal planes at DIBH and FB, and 62,400 total frames were analyzed. Cohort averages of 98% of the maximum cardiac motion ranges comprised margin expansions of 12.5 ± 1.1 mm SI, 5.8 ± 1.2 mm AP, and 6.6 ± 1.0 mm LR at FB and 6.7 ± 1.5 mm SI, 4.7 ± 1.3 mm AP, and 5.3 ± 1.3 mm LR at DIBH. Margins for 67% of the maximum range comprised 7.7 ± 0.7 mm SI, 3.2 ± 0.6 mm AP, and 3.7 ± 0.6 mm LR at FB and 4.1 ± 0.9 mm SI, 2.7 ± 0.8 mm AP, and 3.2 ± 0.8 mm LR at DIBH. Subsequently, these margins were simplified to form PRVs for treatment planning.

CONCLUSIONS

We implemented scale-invariant feature transformation-based motion tracking for analysis of the cardiac cine MRI scans to quantify motion and create cohort-based cardiac PRVs to improve cardioprotection in breast and thoracic radiation.

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.

Real-Time External Respiratory Motion Measuring Technique Using an RGB-D Camera and Principal Component Analysis

Accurate tracking and modeling of internal and external respiratory motion in the thoracic and abdominal regions of a human body is a highly discussed topic in external beam radiotherapy treatment. Errors in target/normal tissue delineation and dose calculation and the increment of the healthy tissues being exposed to high radiation doses are some of the unsolicited problems caused due to inaccurate tracking of the respiratory motion. Many related works have been introduced for respiratory motion modeling, but a majority of them highly depend on radiography/fluoroscopy imaging, wearable markers or surgical node implanting techniques. We, in this article, propose a new respiratory motion tracking approach by exploiting the advantages of an RGB-D camera. First, we create a patient-specific respiratory motion model using principal component analysis (PCA) removing the spatial and temporal noise of the input depth data. Then, this model is utilized for real-time external respiratory motion measurement with high accuracy. Additionally, we introduce a marker-based depth frame registration technique to limit the measuring area into an anatomically consistent region that helps to handle the patient movements during the treatment. We achieved a 0.97 correlation comparing to a spirometer and 0.53 mm average error considering a laser line scanning result as the ground truth. As future work, we will use this accurate measurement of external respiratory motion to generate a correlated motion model that describes the movements of internal tumors.

Development of patient-controlled respiratory gating system based on visual guidance for magnetic-resonance image-guided radiation therapy

PURPOSE

The aim of this study is to develop a visual guidance patient-controlled (VG-PC) respiratory gating system for respiratory-gated magnetic-resonance image-guided radiation therapy (MR-IGRT) and to evaluate the performance of the developed system.

METHODS

The near-real-time cine planar MR image of a patient acquired during treatment was transmitted to a beam projector in the treatment room through an optical fiber cable. The beam projector projected the cine MR images inside the bore of the ViewRay system in order to be visible to a patient during treatment. With this visual information, patients voluntarily controlled their respiration to put the target volume into the gating boundary (gating window). The effect of the presence of the beam projector in the treatment room on the image quality of the MRI was investigated by evaluating the signal-to-noise ratio (SNR), uniformity, low-contrast detectability, high-contrast spatial resolution, and spatial integrity with the VG-PC gating system. To evaluate the performance of the developed system, we applied the VG-PC gating system to a total of seven patients; six patients received stereotactic ablative radiotherapy (SABR) and one patient received conventional fractionated radiation therapy.

RESULTS

The projected cine MR images were visible even when the room light was on. No image data loss or additional time delay during delivery of image data were observed. Every indicator representing MRI quality, including SNR, uniformity, low-contrast detectability, high-contrast spatial resolution, and spatial integrity exhibited values higher than the tolerance levels of the manufacturer with the VG-PC gating system; therefore, the presence of the VG-PC gating system in the treatment room did not degrade the MR image quality. The average beam-off times due to respiratory gating with and without the VG-PC gating system were 830.3 ± 278.2 s and 1264.2 ± 302.1 s respectively (P = 0.005). Consequently, the total treatment times excluding the time for patient setup with and without the VG-PC gating system were 1453.3 ± 297.3 s and 1887.2 ± 469.6 s, respectively, on average (P = 0.005). The average number of beam-off events during whole treatment session was reduced from 457 ± 154 times to 195 ± 90 times by using the VG-PC gating system (P < 0.001).

CONCLUSIONS

The developed system could improve treatment efficiency when performing respiratory-gated MR-IGRT. The VG-PC gating system could be applied to any kind of bore-type radiotherapy machine.

Robust optimization of VMAT for lung cancer: Dosimetric implications of motion compensation techniques

In inverse planning of lung radiotherapy, techniques are required to ensure dose coverage of target disease in the presence of tumor motion as a result of respiration. A range of published techniques for mitigating motion effects were compared for dose stability across 5 breath cycles of ±2 cm. Techniques included planning target volume (PTV) expansions, internal target volumes with (OITV) and without tissue override (ITV), average dataset scans (ADS), and mini-max robust optimization. Volumetric arc therapy plans were created on a thorax phantom and verified with chamber and film measurements. Dose stability was compared by DVH analysis in calculations across all geometries. The lung override technique resulted in a substantial lack of dose coverage (-10%) to the tumor in the presence of large motion. PTV, ITV and ADS techniques resulted in substantial (up to 25%) maximum dose increases where solid tissue travelled into low density optimized regions. The results highlight the need for care in optimization of highly heterogeneous where density variations may occur with motion. Robust optimization was shown to provide greater stability in both maximum (<3%) and minimum dose variations (<2%) over all other techniques.

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

Purpose

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

Method and materials

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

Results

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

Conclusions

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

Systematic evaluation of lung tumor motion using four-dimensional computed tomography

BACKGROUND

Respiratory-induced lung tumor motion may decrease robustness and outcome of radiation therapy (RT) if not accounted for. This study provides detailed information on the motion distribution of lung tumors for a group of 126 patients treated with stereotactic body RT.

MATERIAL AND METHODS

Four-dimensional computed tomography scans were reviewed to assess lung tumor motion. The tumor motion was determined by the center of mass shift based on a rigid registration of the breathing phases containing the largest positional differences in the inferior-superior (IS), left-right (LR), and anterior-posterior (AP) directions. The patients were divided into subgroups depending on tumor diameter (φ < 2.0 cm, 2.0 ≤ φ ≤ 5.0 cm, φ > 5.0 cm) and tumor location within the lung (upper, middle, or lower lobe). The observed motion distributions were evaluated for each group separately to assess the dependence on tumor size and location. For each tumor size, the motion pattern in each direction (IS, LR, and AP) was analyzed for every tumor moving >5 mm. Sinusoidal trigonometric functions were fitted to the measured data using the least mean square method to determine which type of function best describes the motion pattern. Tumor volumes between 1.6 and 52.3 cm³ were evaluated. Mann-Whitney statistical tests were used for statistical analyses.

RESULTS

The mean amplitude for the tumors in this study was 1.5 mm (LR), 2.5 mm (AP), and 6.9 mm (IS) while the maximum amplitude was 11.0 mm (LR), 9.0 mm (AP), and 53.0 mm (IS). In total, 95% of the tumors moved ≤20 mm in the IS direction, ≤3 mm in the LR direction, and ≤6 mm in the AP direction. The observed motion distributions showed no statistically significant correlation with tumor size or location within the lung except for motion in the IS direction, where the mean and maximum amplitudes significantly increased for tumors located in the middle and lower parts of the lung. The motion pattern of a tumor in any direction was best described using a squared trigonometric function of the type [Formula: see text], where A is the maximum amplitude of the motion in the current direction, t is the time of measurement, T is the total time of the breathing cycle and B is a constant used to synchronize the starting point of the breathing cycle.

CONCLUSION

Lung tumor movements were generally larger in the IS direction and the motion amplitude in this direction increased for tumors located in the middle and lower parts of the lungs. Motions in LR or AP showed no such relation. Tumor size was not found to have any correlation with the motion amplitude in any direction. The motion pattern of a lung tumor in any direction is best described with a squared sinusoidal function independently of the tumor size or tumor location.

Clinical respiratory motion correction software (reconstruct, register and averaged-RRA), for 18F-FDG-PET-CT: phantom validation, practical implications and patient evaluation

OBJECTIVE

On fluorine-18 fludeoxyglucose (¹⁸F-FDG) positron emission tomography (PET) CT of pulmonary or hepatic lesions, standard uptake value (SUV) is often underestimated due to patient breathing. The aim of this study is to validate, on phantom and patient data, a motion correction algorithm [reconstruct, register and averaged (RRA)] implemented on a PET-CT system.

METHODS

Three phantoms containing five spheres filled with ¹⁸F-FDG and suspended in a water or Styrofoam^®18F-FDG-filled tank to create different contrasts and attenuation environment were acquired on a Discovery GE710. The spheres were animated with a 2-cm longitudinal respiratory-based movement. Respiratory-gated (RRA) and ungated PET images were compared with static reference images (without movement). The optimal acquisition time, number of phases and the best phase within the respiratory cycle were investigated. The impact of irregular motion was also investigated. Quantification impact was computed on each sphere. Quantification improvement on 28 lung lesions was also investigated.

RESULTS

Phantoms: 4 min was required to obtain a stable quantification with the RRA method. The reference phase and the number of phases used for RRA did not affect the quantification which was similar on static acquisitions but different on ungated images. The results showed that the maximum standard uptake value (SUV_max) restoration is majored for the smallest spheres (≤2.1 ml).

PATIENTS

SUV_max on RRA and ungated acquisitions were statistically different to the SUV_max on whole-body images (p = 0.05) but not different from each other (mean SUV_max: 7.0 ± 7.8 vs 6.9 ± 7.8, p = 0.23 on RRA and ungated images, respectively). We observed a statistically significant correlation between SUV restoration and lesion displacement, with a real SUV quantitation improvement for lesion with movement >1.2 mm.

CONCLUSION

According to the results obtained using phantoms, RRA method is promising, showing a real impact on the lesion quantification on phantom data. With regard to the patient study, our results showed a trend towards an increase in the SUVs and a decrease in the volume between the ungated and RRA data. We also noticed a statistically significant correlation between the quantitative restoration obtained with RRA compared with ungated data and lesion displacement, indicating that the RRA approach should be reserved to patients with small lesions or nodes moving with a displacement larger than 1.2 cm. Advances in knowledge: This article investigates the performances of motion correction software recently introduced in PET. The conclusion revealed that such respiratory motion correction approach shows a real impact on the lesion quantification but must be reserved to the patient for whom lesion displacement was confirmed and high enough to clearly impact lesion evaluation.

Modern radiotherapy using image guidance for unresectable non-small cell lung cancer can improve outcomes in patients treated with chemoradiation therapy

BACKGROUND

To investigate the impact of advances in image-guided radiotherapy (IGRT) on the outcomes of patients with non-small cell lung cancer (NSCLC) treated with chemoradiation therapy (CRT).

METHODS

We retrospectively reviewed 91 patients with NSCLC treated with definitive CRT using image guidance with daily orthogonal kilovoltage (kV) imaging compared to standard weekly megavoltage (MV) portal verifications. Kaplan-Meier curves for overall survival and locoregional failure were computed and stratified by image guidance techniques. Log-rank tests were used to compare strata. Cox Proportional Hazards models were used to identify risk factors for worse mortality and locoregional control.

RESULTS

Fifty-four percent (n=49) of patients received weekly MV portal imaging, while 46% (n=42) underwent IGRT using daily orthogonal kV imaging. kV IGRT was associated with longer median survival (36.4 months) compared to MV imaging (14.9 months; P=0.01). kV imaging was also marginally associated with lower risk of locoregional failure. Median time to local progression in patients imaged with kV was 21.4 months compared to 10.9 months (P=0.065) for those treated with MV portal imaging.

CONCLUSIONS

Daily kV imaging appears to be marginally associated with better survival and disease control when compared to MV imaging. Given the small study size and the numerable factors tested, these finding require additional confirmation.

Impact of microscopic disease extension, extra-CTV tumour islets, incidental dose and dose conformity on tumour control probability

The impact of microscopic disease extension (MDE), extra-CTV tumour islets (TIs), incidental dose and dose conformity on tumour control probability (TCP) is analyzed using insilico simulations in this study. MDE in the region in between GTV and CTV is simulated inclusive of geometric uncertainties (GE) using spherical targets and spherical dose distribution. To study the effect of incidental dose on TIs and the effect of dose-response curve (DRC) on tumour control, islets were randomly distributed and TCP was calculated for various dose levels by rescaling the dose. Further, the impact of dose conformity on required PTV margins is also studied. The required PTV margins are ~2 mm lesser than assuming a uniform clonogen density if an exponential clonogen density fall off in the GTV-CTV is assumed. However, margins are almost equal if GE is higher in both cases. This shows that GE has a profound impact on margins. The effect of TIs showed a bi-phasic relation with increasing dose, indicating that patients with islets not in the beam paths do not benefit from dose escalation. Increasing dose conformity is also found to have considerable effect on TCP loss especially for larger GE. Further, smaller margins in IGRT should be used with caution where uncertainty in CTV definition is of concern.

Analysis of precision in tumor tracking based on optical positioning system during radiotherapy

Tumor tracking is performed during patient set-up and monitoring of respiratory motion in radiotherapy. In the clinical setting, there are several types of equipment for this set-up such as the Electronic Portal imaging Device (EPID) and Cone Beam CT (CBCT). Technically, an optical positioning system tracks the difference between the infra ball reflected from body and machine isocenter. Our objective is to compare the clinical positioning error of patient setup between Cone Beam CT (CBCT) with the Optical Positioning System (OPS), and to evaluate the traditional positioning systems and OPS based on our proposed approach of patient positioning. In our experiments, a phantom was used, and we measured its setup errors in three directions. Specifically, the deviations in the left-to-right (LR), anterior-to-posterior (AP) and inferior-to-superior (IS) directions were measured by vernier caliper on a graph paper using the Varian Linear accelerator. Then, we verified the accuracy of OPS based on this experimental study. In order to verify the accuracy of phantom experiment, 40 patients were selected in our radiotherapy experiment. To illustrate the precise of optical positioning system, we designed clinical trials using EPID. From our radiotherapy procedure, we can conclude that OPS has higher precise than conventional positioning methods, and is a comparatively fast and efficient positioning method with respect to the CBCT guidance system.

Geographic miss of lung tumours due to respiratory motion: a comparison of 3D vs 4D PET/CT defined target volumes

Background

PET/CT scans acquired in the radiotherapy treatment position are typically performed without compensating for respiratory motion. The purpose of this study was to investigate geographic miss of lung tumours due to respiratory motion for target volumes defined on a standard 3D-PET/CT.

Methods

29 patients staged for pulmonary malignancy who completed both a 3D-PET/CT and 4D-PET/CT were included. A 3D-Gross Tumour Volume (GTV) was defined on the standard whole body PET/CT scan. Subsequently a 4D-GTV was defined on a 4D-PET/CT MIP. A 5 mm, 10 mm, 15 mm symmetrical and 15×10 mm asymmetrical Planning Target Volume (PTV) was created by expanding the 3D-GTV and 4D-GTV’s. A 3D conformal plan was generated and calculated to cover the 3D-PTV. The 3D plan was transferred to the 4D-PTV and analysed for geographic miss. Three types of miss were measured. Type 1: any part of the 4D-GTV outside the 3D-PTV. Type 2: any part of the 4D-PTV outside the 3D-PTV. Type 3: any part of the 4D-PTV receiving less than 95% of the prescribed dose. The lesion motion was measured to look at the association between lesion motion and geographic miss.

Results

When a standard 15 mm or asymmetrical PTV margin was used there were 1/29 (3%) Type 1 misses. This increased 7/29 (24%) for the 10 mm margin and 23/29 (79%) for a 5 mm margin. All patients for all margins had a Type 2 geographic miss. There was a Type 3 miss in 25 out of 29 cases in the 5, 10, and 15 mm PTV margin groups. The asymmetrical margin had one additional Type 3 miss. Pearson analysis showed a correlation (p < 0.01) between lesion motion and the severity of the different types of geographic miss.

Conclusion

Without any form of motion suppression, the current standard of a 3D- PET/CT and 15 mm PTV margin employed for lung lesions has an increasing risk of significant geographic miss when tumour motion increases. Use of smaller asymmetric margins in the cranio-caudal direction does not comprise tumour coverage. Reducing PTV margins for volumes defined on 3D-PET/CT will greatly increase the chance and severity of a geometric miss due to respiratory motion. 4D-imaging reduces the risk of geographic miss across the population of tumour sizes and magnitude of motion investigated in the study.

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.

The road to clinical use of high-intensity focused ultrasound for liver cancer: technical and clinical consensus

Clinical use of high-intensity focused ultrasound (HIFU) under ultrasound or MR guidance as a non-invasive method for treating tumors is rapidly increasing. Tens of thousands of patients have been treated for uterine fibroid, benign prostate hyperplasia, bone metastases, or prostate cancer. Despite the methods' clinical potential, the liver is a particularly challenging organ for HIFU treatment due to the combined effect of respiratory-induced liver motion, partial blocking by the rib cage, and high perfusion/flow. Several technical and clinical solutions have been developed by various groups during the past 15 years to compensate for these problems. A review of current unmet clinical needs is given here, as well as a consensus from a panel of experts about technical and clinical requirements for upcoming pilot and pivotal studies in order to accelerate the development and adoption of focused ultrasound for the treatment of primary and secondary liver cancer.

A novel technique for markerless, self-sorted 4D-CBCT: feasibility study

PURPOSE

Four-dimensional CBCT (4D-CBCT) imaging in the treatment room can provide verification of moving targets, facilitating the potential for margin reduction and consequent dose escalation. Reconstruction of 4D-CBCT images requires correlation of respiratory phase with projection acquisition, which is often achieved with external surrogate measures of respiration. However, external measures may not be a direct representation of the motion of the internal anatomy and it is therefore the aim of this work to develop a novel technique for markerless, self-sorted 4D-CBCT reconstruction.

METHODS

A novel 4D-CBCT reconstruction technique based on the principles of Fourier transform (FT) theory was investigated for markerless extraction of respiratory phase directly from projection data. In this FT technique, both phase information (FT-phase) and magnitude information (FT-magnitude) were separately implemented in order to discern projections corresponding to peak inspiration, which then facilitated the proceeding sort and bin processes involved in retrospective 4D image reconstruction. In order to quantitatively evaluate the accuracy of the Fourier methods, peak-inspiration projections identified each by FT-phase and FT-magnitude were compared to those manually identified by visual tracking of structures. The average phase difference as assigned by each method vs the manual technique was calculated per projection dataset. The percentage of projections that were assigned within 10% phase of each other was also computed. Both Fourier methods were tested on two phantom datasets, programmed to exhibit sinusoidal respiratory cycles of 2.0 cm in amplitude with respiratory cycle lengths of 3 and 6 s, respectively. Additionally, three sets of patient projections were studied. All of the data were previously acquired at slow-gantry speeds ranging between 0.6°/s and 0.7°/s over a 200° rotation. Ten phase bins with 10% phase windows were selected for 4D-CBCT reconstruction of one phantom and one patient case for visual and quantitative comparison. Line profiles were plotted for the 0% and 50% phase images as reconstructed by the manual technique and each of the Fourier methods.

RESULTS

As compared with the manual technique, the FT-phase method resulted in average phase differences of 1.8% for the phantom with the 3 s respiratory cycle, 3.9% for the phantom with the 6 s respiratory cycle, 2.9% for patient 1, 5.0% for patient 2, and 3.8% for patient 3. For the FT-magnitude method, these numbers were 2.1%, 4.0%, 2.9%, 5.3%, and 3.5%, respectively. The percentage of projections that were assigned within 10% phase by the FT-phase method as compared to the manual technique for the five datasets were 100.0%, 100.0%, 97.6%, 93.4%, and 94.1%, respectively, whereas for the FT-magnitude method these percentages were 98.1%, 92.3%, 98.7%, 87.3%, and 95.7%. Reconstructed 4D phase images for both the phantom and patient case were visually and quantitatively equivalent between each of the Fourier methods vs the manual technique.

CONCLUSIONS

A novel technique employing the basics of Fourier transform theory was investigated and demonstrated to be feasible in achieving markerless, self-sorted 4D-CBCT reconstruction.

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

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

Potential benefits and pitfalls of respiratory-gated radiotherapy in the treatment of thoracic malignancy

AIM

Despite advances in radiotherapy delivery, the prognosis of lung cancer remains poor. Higher doses of radiation have been associated with improved outcomes but may result in higher toxicities. Respiratory gated radiotherapy (RGRT) has the potential to reduce pulmonary toxicity but there are significant limitations and pitfalls to its use. The aim of this article is to (i) describe the RGRT technique currently employed at Nepean and Westmead Hospitals; (ii) discuss the practical issues of implementing such a program; (iii) present the results of our RGRT program and (iv) review the potential uncertainties in using this technique and the methods we have used to overcome these.

METHODS

A retrospective review of all patients who had a 4D-computed tomography (4D-CT) scan was undertaken. Records from treatment planning systems were used to assess the prospective gating program.

RESULTS

Between September 2007 and June 2011, 53 patients at Nepean and 26 patients at Westmead Hospital underwent a 4D-CT. Between April and August 2011, 26 patients at Westmead Hospital underwent a prospective 4D-CT scan as treatment verification. Two of the 26 patients (7.7%) were found to have incomplete coverage of the planning target volume. Both patients underwent respiratory re-coaching, alleviating the need for replanning.

CONCLUSION

RGRT may reduce doses to organs at risk with the potential for dose escalation. However its implementation requires significant staff training, treatment time and resources. Treatment verification with image guided radiation therapy are essential for safe delivery.

Development and validation of a prototypal neural networks-based tumor tracking method

In radiotherapy, intra-fractional organ motion introduces uncertainties in target localization, leading to unacceptable inaccuracy in dose delivery. Especially in highly selective treatments, such as those delivered with particles beams instead of photons, organ motion may results in severe side effects and/or limited tumor control. Tumor tracking is a motion mitigation strategy that allows an almost continuous dose delivery while the beam is dynamically steered to match the position of the moving target in real-time. Currently, tumor tracking is applied clinically only in the CyberKnife system for photon radiotherapy, whereas neither clinical solutions nor dedicated methodologies are available for particle therapy. Consequently, the aim of the proposed study is to develop a neural networks-based prototypal tracking algorithm intended for particle therapy. We developed a method that exploits three independent neural networks to estimate the internal target position as a function of external surrogate signals. This method was tested on data relative to 20 patients treated with CyberKnife, whose performance was used as benchmark. Results show that the developed algorithm allows targeting error reduction with respect to the CyberKnife system, thus proving the potential value of artificial neural networks for the implementation of tumor tracking methodologies.

Determination of internal target volume from a single positron emission tomography/computed tomography scan in lung cancer

PURPOSE

The use of four-dimensional computed tomography (4D-CT) to determine the tumor internal target volume (ITV) is usually characterized by high patient radiation exposure. The objective of this study was to propose and evaluate an approach that relies on a single static positron emission tomography (PET)/CT scan to determine the ITV, thereby eliminating the need for 4D-CT and thus reduce patient radiation dose.

METHODS AND MATERIALS

The proposed approach is based on the concept that the observed PET image is the result of a joint convolution of an ideal PET image (free from motion and partial volume effect) with a motion-blurring kernel (MBK) and partial volume effect. In this regard, the MBK and tumor ITV are then estimated from the deconvolution of this joint model. To test this technique, phantom and patient studies were performed using different sphere/tumor sizes and motion trajectories. In all studies, a 4D-CT and a PET/CT image of the sphere/tumor were acquired. The ITV from the proposed technique was then compared to the maximum intensity projection (MIP) volume of the 4D-CT images. A Dice coefficient of the two volumes was calculated to represent the similarity between the two ITVs.

RESULTS

The average ITVs of the proposed technique were 97.2% ± 0.3% and 81.0% ± 16.7% similar to the MIP volume in the phantom and patient studies, respectively. The average dice coefficients were 0.87 ± 0.05 and 0.73 ± 0.16, respectively, for the two studies.

CONCLUSION

Using the proposed approach, a single static PET/CT scan has the potential to replace a 4D-CT to determine the tumor ITV. This approach has the added advantage of reducing patient radiation exposure and determining the tumor MBK compared to 4D-CT/MIP-CT.

Evaluation and treatment of high-risk patients with early-stage lung cancer

Standard therapy for early-stage non-small cell lung cancer is lobectomy for patients who are able to tolerate such surgery. However, the risk of postoperative morbidity is not trivial, with a 30% to 40% incidence of postoperative complications and a 1% to 5% incidence of operative mortality. Some patients, though technically resectable, refuse surgery or are considered medically inoperable because of insufficient respiratory reserve, cardiovascular disease, or general frailty. This group is considered either "high risk" or "medically inoperable."

The management of tumor motions in the stereotactic irradiation to lung cancer under the use of Abches to control active breathing

PURPOSE

Breathing control is crucial to ensuring the accuracy of stereotactic irradiation for lung cancer. This study monitored respiration in patients with inoperable nonsmall-cell lung cancer using a respiration-monitoring apparatus, Abches, and investigated the reproducibility of tumor position in these patients.

METHODS

Subjects comprised 32 patients with nonsmall-cell lung cancer who were administered stereotactic radiotherapy under breath-holding conditions monitored by Abches. Computed tomography (CT) was performed under breath-holding conditions using Abches (Abches scan) for treatment planning. A free-breathing scan was performed to determine the range of tumor motions in a given position. After the free-breathing scan, Abches scan was repeated and the tumor position thus defined was taken as the intrafraction tumor position. Abches scan was also performed just before treatment, and the tumor position thus defined was taken as the interfraction tumor position. To calculate the errors, tumor positions were compared based on Abches scan for the initial treatment plan. The error in tumor position was measured using the BrainSCAN treatment-planning device, then compared for each lung lobe.

RESULTS

Displacements in tumor position were calculated in three dimensions (i.e., superior-inferior (S-I), left-right (L-R), and anterior-posterior (A-P) dimensions) and recorded as absolute values. For the whole lung, average intrafraction tumor displacement was 1.1 mm (L-R), 1.9 mm (A-P), and 2.0 mm (S-I); the average interfraction tumor displacement was 1.1 mm (L-R), 2.1 mm (A-P), and 2.0 mm (S-I); and the average free-breathing tumor displacement was 2.3 mm (L-R), 3.5 mm (A-P), and 7.9 mm (S-I). The difference between using Abches and free breathing could be reduced from approximately 20 mm at the maximum to approximately 3 mm in the S-I direction for both intrafraction and interfraction positions in the lower lobe. In addition, maximum intrafraction tumor displacement with the use of Abches was 4.5 mm (S-I) in the lingular segment. These results suggest that use of the Abches system can reduce deviations in tumor position to levels below those achieved under free breathing, irrespective of the tumor location.

CONCLUSIONS

Respiratory control with high accuracy and reproducibility is required for high-precision radiotherapy of inoperable nonsmall-cell lung cancer and was achieved using Abches in this study.

Measurement of time delay for a prospectively gated CT simulator

For the management of mobile tumors, respiratory gating is the ideal option, both during imaging and during therapy. The major advantage of respiratory gating during imaging is that it is possible to create a single artifact-free CT data-set during a selected phase of the patient's breathing cycle. The purpose of the present work is to present a simple technique to measure the time delay during acquisition of a prospectively gated CT. The time delay of a Philips Brilliance BigBore™ (Philips Medical Systems, Madison, WI) scanner attached to a Varian Real-Time Position Management™ (RPM) system (Varian Medical Systems, Palo Alto, CA) was measured. Two methods were used to measure the CT time delay: using a motion phantom and using a recorded data file from the RPM system. In the first technique, a rotating wheel phantom was altered by placing two plastic balls on its axis and rim, respectively. For a desired gate, the relative positions of the balls were measured from the acquired CT data and converted into corresponding phases. Phase difference was calculated between the measured phases and the desired phases. Using period of motion, the phase difference was converted into time delay. The Varian RPM system provides an external breathing signal; it also records transistor-transistor logic (TTL) ‘X-Ray ON’ status signal from the CT scanner in a text file. The TTL ‘X-Ray ON’ indicates the start of CT image acquisition. Thus, knowledge of the start time of CT acquisition, combined with the real-time phase and amplitude data from the external respiratory signal, provides time-stamping of all images in an axial CT scan. The TTL signal with time-stamp was used to calculate when (during the breathing cycle) a slice was recorded. Using the two approaches, the time delay between the prospective gating signal and CT simulator has been determined to be 367 ± 40 ms. The delay requires corrections both at image acquisition and while setting gates for the treatment delivery; otherwise the simulation and treatment may not be correlated with the patient's breathing.