Frequency filtering based analysis on the cardiac induced lung tumor motion and its impact on the radiotherapy management

Enhancing Image-Guided Radiation Therapy for Pancreatic Cancer: Utilizing Aligned Peak Response Beamforming in Flexible Array Transducers

To develop ultrasound-guided radiotherapy, we proposed an assistant structure with embedded markers along with a novel alternative method, the Aligned Peak Response (APR) method, to alter the conventional delay-and-sum (DAS) beamformer for reconstructing ultrasound images obtained from a flexible array. We simulated imaging targets in Field-II using point target phantoms with point targets at different locations. In the experimental phantom ultrasound images, image RF data were acquired with a flexible transducer with in-house assistant structures embedded with needle targets for testing the accuracy of the APR method. The lateral full width at half maximum (FWHM) values of the objective point target (OPT) in ground truth ultrasound images, APR-delayed ultrasound images with a flat shape, and images acquired with curved transducer radii of 500 mm and 700 mm were 3.96 mm, 4.95 mm, 4.96 mm, and 4.95 mm. The corresponding axial FWHM values were 1.52 mm, 4.08 mm, 5.84 mm, and 5.92 mm, respectively. These results demonstrate that the proposed assistant structure and the APR method have the potential to construct accurate delay curves without external shape sensing, thereby enabling a flexible ultrasound array for tracking pancreatic tumor targets in real time for radiotherapy.

Adaptive hypofractionted and stereotactic body radiotherapy for lung tumors with real-time MRI guidance

The treatment of central and ultracentral lung tumors with radiotherapy remains an ongoing clinical challenge. The risk of Grade 5 toxicity with ablative radiotherapy doses to these high-risk regions is significant as shown in recent prospective studies. Magnetic resonance (MR) image-guided adaptive radiotherapy (MRgART) is a new technology and may allow the delivery of ablative radiotherapy to these high-risk regions safely. MRgART is able to achieve this by utilizing small treatment margins, real-time gating/tracking and on-table plan adaptation to maintain dose to the tumor but limit dose to critical structures. The process of MRgART is complex and has nuances and challenges for the treatment of lung tumors. We outline the critical steps needed for appropriate delivery of MRgART for lung tumors safely and effectively.

CArdiac and REspiratory adaptive Computed Tomography (CARE-CT): a proof-of-concept digital phantom study

Current respiratory 4DCT imaging for high-dose rate thoracic radiotherapy treatments are negatively affected by the complex interaction of cardiac and respiratory motion. We propose an imaging method to reduce artifacts caused by thoracic motion, CArdiac and REspiratory adaptive CT (CARE-CT), that monitors respiratory motion and ECG signals in real-time, triggering CT acquisition during combined cardiac and respiratory bins. Using a digital phantom, conventional 4DCT and CARE-CT acquisitions for nineteen patient-measured physiological traces were simulated. Ten respiratory bins were acquired for conventional 4DCT scans and ten respiratory bins during cardiac diastole were acquired for CARE-CT scans. Image artifacts were quantified for 10 common thoracic organs at risk (OAR) substructures using the differential normalized cross correlation between axial slices (ΔNCC), mean squared error (MSE) and sensitivity. For all images, on average, CARE-CT improved the ΔNCC for 18/19 and the MSE and sensitivity for all patient traces. The ΔNCC was reduced for all cardiac OARs (mean reduction 21%). The MSE was reduced for all OARs (mean reduction 36%). In the digital phantom study, the average scan time was increased from 1.8 ± 0.4 min to 7.5 ± 2.2 min with a reduction in average beam on time from 98 ± 28 s to 45 s using CARE-CT compared to conventional 4DCT. The proof-of-concept study indicates the potential for CARE-CT to image the thorax in real-time during the cardiac and respiratory cycle simultaneously, to reduce image artifacts for common thoracic OARs.

Intrafraction target shift comparison using two breath-hold systems in lung stereotactic body radiotherapy

Background and purpose

In lung Stereotactic Body Radiotherapy (SBRT) respiratory management is used to reduce target motion due to respiration. This study aimed (1) to estimate intrafraction shifts through a Cone Beam Computed Tomography (CBCT) acquired during the first treatment arc when deep inspiration breath-hold (DIBH) was performed using spirometry-based (SB) or surface-guidance (SG) systems and (2) to analyze the obtained results depending on lesion localization.

Material and methods

A sample of 157 patients with 243 lesions was analyzed. A total of 860 and 410 fractions were treated using SB and SG. Averaged intrafraction shifts were estimated by the offsets obtained when registering a CBCT acquired during the first treatment arc with the planning CT. Offsets were recorded in superior-inferior (SI), left-right (LR) and anterior-posterior (AP). Significance tests were applied to account for differences in average offsets and variances between DIBH systems. Systematic and random errors were computed.

Results

Average offset moduli were 2.4 ± 2.2 mm and 3.5 ± 2.6 mm for SB and SG treatments (p < 0.001). When comparing SB and SG offset distributions in each direction no differences were found in average values (p > 0.3). However, variances were statistically smaller for SB than for SG (p < 0.001). The number of vector moduli offsets greater than 5 mm was 2.1 times higher for SG. Compared to other locations, lower lobe lesions moduli were at least 2.3 times higher.

Conclusions

Both systems were accuracy-equivalent but not precision-equivalent systems. Furthermore, the SB system was more precise than the SG one. Despite DIBH, patients with lower lobe lesions had larger offsets than superior lobe ones, mainly in SI.

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.

Motion correction for routine X-ray lung CT imaging

A novel motion correction algorithm for X-ray lung CT imaging has been developed recently. It was designed to perform for routine chest or thorax CT scans without gating, namely axial or helical scans with pitch around 1.0. The algorithm makes use of two conjugate partial angle reconstruction images for motion estimation via non-rigid registration which is followed by a motion compensated reconstruction. Differently from other conventional approaches, no segmentation is adopted in motion estimation. This makes motion estimation of various fine lung structures possible. The aim of this study is to explore the performance of the proposed method in correcting the lung motion artifacts which arise even under routine CT scans with breath-hold. The artifacts are known to mimic various lung diseases, so it is of great interest to address the problem. For that purpose, a moving phantom experiment and clinical study (seven cases) were conducted. We selected the entropy and positivity as figure of merits to compare the reconstructed images before and after the motion correction. Results of both phantom and clinical studies showed a statistically significant improvement by the proposed method, namely up to 53.6% (p < 0.05) and up to 35.5% (p < 0.05) improvement by means of the positivity measure, respectively. Images of the proposed method show significantly reduced motion artifacts of various lung structures such as lung parenchyma, pulmonary vessels, and airways which are prominent in FBP images. Results of two exemplary cases also showed great potential of the proposed method in correcting motion artifacts of the aorta which is known to mimic aortic dissection. Compared to other approaches, the proposed method provides an excellent performance and a fully automatic workflow. In addition, it has a great potential to handle motions in wide range of organs such as lung structures and the aorta. We expect that this would pave a way toward innovations in chest and thorax CT imaging.

[Uncertainties in the current concept of radiotherapy planning target volume]

The planning target volume is an essential notion in radiotherapy, that requires a new conceptualization. Indeed, the variability and diversity of the uncertainties involved or improved with the development of the new modern technologies and devices in radiotherapy suggest that random and systematic errors cannot be currently generalized. This article attempts to discuss these various uncertainties and tries to demonstrate that a redefinition of the concept of planning target volume toward its personalization for each patient and the robustness notion are likely an improvement basis to take into account the radiotherapy uncertainties.

Comparison of different methods for lung immobilization in an animal model

BACKGROUND AND PURPOSE

Respiratory-induced motion introduces uncertainties in the delivery of dose in radiotherapy treatments. Various methods are used clinically, e.g. breath-holding, while there is limited experience with other methods such as apneic oxygenation and high frequency jet ventilation (HFJV). This study aims to compare the latter approaches for lung immobilization and their clinical impact on gas exchange in an animal model.

MATERIALS AND METHODS

Two radiopaque tumor surrogate markers (TSM) were placed in the central (cTSM) and peripheral (dTSM) regions of the lungs in 9 anesthetized and muscle relaxed pigs undergoing 3 ventilatory interventions (1) HFJV at rates of 200 (JV200), 300 (JV300) and 400 (JV400) min^-1; (2) apnea at continuous positive airway pressure (CPAP) levels of 0, 8 and 16 cmH₂O; (3) conventional mechanical ventilation (CMV) as reference mode. cTSM and dTSM were visualized using fluoroscopy and their coordinates were computed. The ventilatory pattern was registered, and oxygen and carbon dioxide (pCO₂) partial pressures were measured.

RESULTS

The highest range of TSM motion, and ventilation was found during CMV, the lowest during apnea. During HFJV the amount of motion varied inversely with increasing frequency. The reduction of TSM motion at JV300, JV400 and all CPAP levels came at the cost of increased pCO₂, however the relatively low frequency of 200 min^-1 for HFJV was the only ventilatory setting that enabled adequate CO₂ removal.

CONCLUSION

In this model, HFJV at 200 min^-1 was the best compromise between immobilization and gas exchange for sessions of 10-min duration.

Real-time intrafraction motion monitoring in external beam radiotherapy

Radiotherapy (RT) aims to deliver a spatially conformal dose of radiation to tumours while maximizing the dose sparing to healthy tissues. However, the internal patient anatomy is constantly moving due to respiratory, cardiac, gastrointestinal and urinary activity. The long term goal of the RT community to 'see what we treat, as we treat' and to act on this information instantaneously has resulted in rapid technological innovation. Specialized treatment machines, such as robotic or gimbal-steered linear accelerators (linac) with in-room imaging suites, have been developed specifically for real-time treatment adaptation. Additional equipment, such as stereoscopic kilovoltage (kV) imaging, ultrasound transducers and electromagnetic transponders, has been developed for intrafraction motion monitoring on conventional linacs. Magnetic resonance imaging (MRI) has been integrated with cobalt treatment units and more recently with linacs. In addition to hardware innovation, software development has played a substantial role in the development of motion monitoring methods based on respiratory motion surrogates and planar kV or Megavoltage (MV) imaging that is available on standard equipped linacs. In this paper, we review and compare the different intrafraction motion monitoring methods proposed in the literature and demonstrated in real-time on clinical data as well as their possible future developments. We then discuss general considerations on validation and quality assurance for clinical implementation. Besides photon RT, particle therapy is increasingly used to treat moving targets. However, transferring motion monitoring technologies from linacs to particle beam lines presents substantial challenges. Lessons learned from the implementation of real-time intrafraction monitoring for photon RT will be used as a basis to discuss the implementation of these methods for particle RT.

Characterizing spatiotemporal information loss in sparse-sampling-based dynamic MRI for monitoring respiration-induced tumor motion in radiotherapy

PURPOSE

Sparse-sampling and reconstruction techniques represent an attractive strategy to achieve faster image acquisition speeds, while maintaining adequate spatial resolution and signal-to-noise ratio in rapid magnetic resonance imaging (MRI). The authors investigate the use of one such sequence, broad-use linear acquisition speed-up technique (k-t BLAST) in monitoring tumor motion for thoracic and abdominal radiotherapy and examine the potential trade-off between increased sparsification (to increase imaging speed) and the potential loss of "true" information due to greater reliance on a priori information.

METHODS

Lung tumor motion trajectories in the superior-inferior direction, previously recorded from ten lung cancer patients, were replayed using a motion phantom module driven by an MRI-compatible motion platform. Eppendorf test tubes filled with water which serve as fiducial markers were placed in the phantom. The modeled rigid and deformable motions were collected in a coronal image slice using balanced fast field echo in conjunction with k-t BLAST. Root mean square (RMS) error was used as a metric of spatial accuracy as measured trajectories were compared to input data. The loss of spatial information was characterized for progressively increasing acceleration factor from 1 to 16; the resultant sampling frequency was increased approximately from 2.5 to 19 Hz when the principal direction of the motion was set along frequency encoding direction. In addition to the phantom study, respiration-induced tumor motions were captured from two patients (kidney tumor and lung tumor) at 13 Hz over 49 s to demonstrate the impact of high speed motion monitoring over multiple breathing cycles. For each subject, the authors compared the tumor centroid trajectory as well as the deformable motion during free breathing.

RESULTS

In the rigid and deformable phantom studies, the RMS error of target tracking at the acquisition speed of 19 Hz was approximately 0.3-0.4 mm, which was smaller than the reconstructed pixel resolution of 0.67 mm. In the patient study, the dynamic 2D MRI enabled the monitoring of cycle-to-cycle respiratory variability present in the tumor position. It was seen that the range of centroid motion as well as the area covered due to target motion during each individual respiratory cycle was underestimated compared to the entire motion range observed over multiple breathing cycles.

CONCLUSIONS

The authors' initial results demonstrate that sparse-sampling- and reconstruction-based dynamic MRI can be used to achieve adequate image acquisition speeds without significant information loss for the task of radiotherapy guidance. Such monitoring can yield spatial and temporal information superior to conventional offline and online motion capture methods used in thoracic and abdominal radiotherapy.