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

Research in methodologies for modelling the oral cavity

The paper aims to explore the current state of understanding surrounding in silico oral modelling. This involves exploring methodologies, technologies and approaches pertaining to the modelling of the whole oral cavity; both internally and externally visible structures that may be relevant or appropriate to oral actions. Such a model could be referred to as a 'complete model' which includes consideration of a full set of facial features (i.e. not only mouth) as well as synergistic stimuli such as audio and facial thermal data. 3D modelling technologies capable of accurately and efficiently capturing a complete representation of the mouth for an individual have broad applications in the study of oral actions, due to their cost-effectiveness and time efficiency. This review delves into the field of clinical phonetics to classify oral actions pertaining to both speech and non-speech movements, identifying how the various vocal organs play a role in the articulatory and masticatory process. Vitaly, it provides a summation of 12 articulatory recording methods, forming a tool to be used by researchers in identifying which method of recording is appropriate for their work. After addressing the cost and resource-intensive limitations of existing methods, a new system of modelling is proposed that leverages external to internal correlation modelling techniques to create a more efficient models of the oral cavity. The vision is that the outcomes will be applicable to a broad spectrum of oral functions related to physiology, health and wellbeing, including speech, oral processing of foods as well as dental health. The applications may span from speech correction, designing foods for the aging population, whilst in the dental field we would be able to gain information about patient's oral actions that would become part of creating a personalised dental treatment plan.

Association of Internal and External Motion Based on Cine MR Images Acquired During Real-Time Treatment on MRI-Guided Linear Accelerator for Patients With Lung Cancer

Purpose

With the recent clinical implementation of magnetic resonance imaging (MRI)-guided linear accelerators, a large number of real-time planar MR images has been acquired during lung cancer treatment as a standard of care. In this study, associations among lung tumor, diaphragm, and external skin movement were studied based on MR cine imaging during the entire duration of each treatment fraction.

Methods and Materials

This retrospective study used 181,798 planar MRI frames acquired over 55 treatment/imaging sessions of 13 patients with lung cancer treated on 2 MRI-guided linear accelerators. From each planar MR image frame, in-house software automatically extracted 9 features: the superior-interior (SI) and posterior-anterior (PA) positions of a lung tumor; the area of the lung (Lung_Area); the posterior (Dia_Post), dome/apex (Dia_Dome), and anterior (Dia_Ant) points of a diaphragmatic curve; the diaphragm curve point (Dia_Max); and the chest (Chest) and belly (Belly) skin points experienced the maximum range of motions. Correlation analyses were performed among the 9 features for every session. Lung tumor motion range and standard deviations were calculated based on positions obtained in cine images and compared with motion ranges obtained from 4-dimensional computed tomography images.

Results

In the study, 177,009 frames of images were successfully analyzed. For all patients, correlation coefficients were as follows: 0.91 ± 0.10 between any 2 features among Lung_Area, Dia_Post, Dia_Dome, and Dia_Max; 0.82 ± 0.21 between SI and any feature among Lung_Area, Dia_Post, Dia_Dome, and Dia_Max; 0.75 ± 0.24 between SI and Belly. Six of 13 patients were considered large amplitude motion (patients with lung tumor SI motion standard deviation >5 mm). Furthermore, 92,956 frames of images were analyzed for the 6 large-amplitude motion patients. For this set, correlation coefficients were 0.93 ± 0.07 between any 2 features among Lung_Area, Dia_Post, Dia_Dome, and Dia_Max; 0.94 ± 0.06 between SI and any feature among Lung_Area, Dia_Post, Dia_Dome, and Dia_Max; and 0.90 ± 0.09 between SI and Belly.

Conclusions

Both belly and diaphragmatic motions as assessed by cine MRI are highly correlated with large amplitude lung tumor motion in the longitudinal axis.

Machine learning for predicting accuracy of lung and liver tumor motion tracking using radiomic features

Background

Internal tumor motion is commonly predicted using external respiratory signals. However, the internal/external correlation is complex and patient-specific. The purpose of this study was to develop various models based on the radiomic features of computed tomography (CT) images to predict the accuracy of tumor motion tracking using external surrogates and to find accurate and reliable tracking algorithms.

Methods

Images obtained from a total of 108 and 71 patients pathologically diagnosed with lung and liver cancers, respectively, were examined. Real-time position monitoring motion was fitted to tumor motion, and samples with fitting errors greater than 2 mm were considered positive. Radiomic features were extracted from internal target volumes of average intensity projections, and cross-validation least absolute shrinkage and selection operator (LassoCV) was used to conduct feature selection. Based on the radiomic features, a total of 26 separate models (13 for the lung and 13 for the liver) were trained and tested. Area under the receiver operating characteristic curve (AUC), sensitivity, and specificity were used to assess performance. Relative standard deviation was used to assess stability.

Results

Thirty-three and 22 radiomic features were selected for the lung and liver, respectively. For the lung, the AUC varied from 0.848 (decision tree) to 0.941 [support vector classifier (SVC), logistic regression]; sensitivity varied from 0.723 (extreme gradient boosting) to 0.848 [linear support vector classifier (linearSVC)]; specificity varied from 0.834 (gaussian naive bayes) to 0.936 [multilayer perceptron (MLP), wide and deep (W&D)]; and MLP and W&D had better performance and stability than the median. For the liver, the AUC varied from 0.677 [light gradient boosting machine (Light)] to 0.892 (logistic regression); sensitivity varied from 0.717 (W&D) to 0.862 (MLP); specificity varied from 0.566 (Light) to 0.829 (linearSVC); and logistic regression, MLP, and SVC had better performance and stability than the median.

Conclusions

Respiratory-sensitive radiomic features extracted from CT images of lung and liver tumors were proved to contain sufficient information to establish an external/internal motion relationship. We developed a rapid and accurate method based on radiomics to classify the accuracy of monitoring a patient's external surface for lung and liver tumor tracking. Several machine learning algorithms-in particular, MLP-demonstrated excellent classification performance and stability.

Focused Ultrasound and Ultrasound Stimulated Microbubbles in Radiotherapy Enhancement for Cancer Treatment

Radiation therapy (RT) has been the standard of care for treating a multitude of cancer types. However, ionizing radiation has adverse short and long-term side effects which have resulted in treatment complications for decades. Thus, advances in enhancing the effects of RT have been the primary focus of research in radiation oncology. To avoid the usage of high radiation doses, treatment modalities such as high-intensity focused ultrasound can be implemented to reduce the radiation doses required to destroy cancer cells. In the past few years, the use of focused ultrasound (FUS) has demonstrated immense success in a number of applications as it capitalizes on spatial specificity. It allows ultrasound energy to be delivered to a targeted focal area without harming the surrounding tissue. FUS combined with RT has specifically demonstrated experimental evidence in its application resulting in enhanced cell death and tumor cure. Ultrasound-stimulated microbubbles have recently proved to be a novel way of enhancing RT as a radioenhancing agent on its own, or as a delivery vector for radiosensitizing agents such as oxygen. In this mini-review article, we discuss the bio-effects of FUS and RT in various preclinical models and highlight the applicability of this combined therapy in clinical settings.

Association Between Internal Organ/Liver Tumor and External Surface Motion From Cine MR Images on an MRI-Linac

Purposes/Objectives

Historically, motion correlation between internal tumor and external surrogates have been based on limited sets of X-ray or magnetic resonance (MR) images. With the recent clinical implementation of MR-guided linear accelerators, a vast quantity of continuous planar real-time MR imaging data is acquired. In this study, information was extracted from MR cine imaging during liver cancer treatments to establish associations between internal tumor/diaphragm and external surface/skin movement.

Methods and Materials

This retrospective study used 305,644 MR image frames acquired over 118 treatment/imaging sessions of the first 23 liver cancer patients treated on an MRI-linac. 9 features were automatically determined on each MR image frame: Lung_Area, the posterior (Dia_Post), dome (Dia_Dome), and anterior (Dia_Ant) points of a diaphragmatic curve and the diaphragm curve point (Dia_Max), the chest (Chest) and the belly (Belly) skin points experiencing the maximum motion ranges; the superior-interior (SI) and posterior-anterior (PA) positions of a target. For every session, correlation analyses were performed twice among the 9 features: 1) over a breath-hold (BH) set and 2) on a pseudo free-breathing (PFB) generated by removing breath-holding frames.

Results

303,123 frames of images were successfully analyzed. For BH set analysis, correlation coefficients were as follows: 0.94 ± 0.07 between any two features among Dia_Post, Dia_Dome, Dia_Max, and Lung_Area; 0.95 ± 0.06 between SI and any feature among Dia_Post, Dia_Dome, Dia_Max, or Lung_Area; 0.76 ± 0.29 between SI and Belly (with 50% of correlations ≥ 0.87). The PFB set had 142,862 frames of images. For this set, correlation coefficients were 0.96 ± 0.06 between any two features among Dia_Post, Dia_Dome, Dia_Max, and Lung_Area; 0.95 ± 0.06 between SI and any feature among Dia_Post, Dia_Dome, Dia_Max, or Lung_Area; 0.80 ± 0.26 between SI and Belly (with 50% of correlations ≥ 0.91).

Conclusion

Diaphragmatic motion as assessed by cine MR imaging is highly correlated with liver tumor motion. Belly vertical motion is highly correlated with liver tumor longitudinal motion in approximately half of the cases. More detailed analyses of those cases displaying weak correlations are in progress.

Quantification of the geometric uncertainty when using implanted markers as a surrogate for lung tumor motion

BACKGROUND

Fiducial markers are used as surrogates for tumor location during radiation therapy treatment. Developments in lung fiducial marker and implantation technology have provided a means to insert markers endobronchially for tracking of lung tumors. This study quantifies the surrogacy uncertainty (SU) when using endobronchially implanted markers as a surrogate for lung tumor position.

METHODS

We evaluated SU for 17 patients treated in a prospective electromagnetic-guided MLC tracking trial. Tumor and markers were segmented on all phases of treatment planning 4DCTs and all frames of pretreatment kilovoltage fluoroscopy acquired from lateral and frontal views. The difference in tumor and marker position relative to end-exhale position was calculated as the SU for both imaging methods and the distributions of uncertainties analyzed.

RESULTS

The mean (range) tumor motion amplitude in the 4DCT scan was 5.9 mm (1.7-11.7 mm) in the superior-inferior (SI) direction, 2.2 mm (0.9-5.5 mm) in the left-right (LR) direction, and 3.9 mm (1.2-12.9 mm) in the anterior-posterior (AP) direction. Population-based analysis indicated symmetric SU centered close to 0 mm, with maximum 5th/95th percentile values over all axes of -2.0 mm/2.1 mm with 4DCT, and -2.3/1.3 mm for fluoroscopy. There was poor correlation between the SU measured with 4DCT and that measured with fluoroscopy on a per-patient basis. We observed increasing SU with increasing surrogate motion. Based on fluoroscopy analysis, the mean (95% CI) SU was 5% (2%-8%) of the motion magnitude in the SI direction, 16% (6%-26%) of the motion magnitude in the LR direction, and 33% (23%-42%) of the motion magnitude in the AP direction. There was no dependence of SU on marker distance from the tumor.

CONCLUSION

We have quantified SU due to use of implanted markers as surrogates for lung tumor motion. Population 95th percentile range are up to 2.3 mm, indicating the approximate contribution of SU to total geometric uncertainty. SU was relatively small compared with the SI motion, but substantial compared with LR and AP motion. Due to uncertainty in estimations of patient-specific SU, it is recommended that population-based margins are used to account for this component of the total geometric uncertainty.

Towards Non-Invasive Lung Tumor Tracking Based on Patient Specific Model of Respiratory System

The goal of this paper is to calculate a complex internal respiratory and tumoral movements by measuring respiratory air flows and thorax movements. In this context, we present a new lung tumor tracking approach based on a patient-specific biomechanical model of the respiratory system, which takes into account the physiology of respiratory motion to simulate the real non-reproducible motion. The behavior of the lungs, is directly driven by the simulated actions of the breathing muscles, i.e. the diaphragm and the intercostal muscles (the rib cage). In this paper, the lung model is monitored and controlled by a personalized lung pressure/volume relationship during a whole respiratory cycle. The lung pressure and rib kinematics are patient specific and obtained by surrogate measurement. The rib displacement corresponding to the transformation which was computed by finite helical axis method from the end of exhalation (EE) to the end of inhalation (EI). The lung pressure is calculated by an optimization framework based on inverse finite element analysis, by minimizing the lung volume errors, between the respiratory volume (respiratory airflow exchange) and the simulated volume (calculated by biomechanical simulation). We have evaluated the model accuracy on five public datasets. We have also evaluated the lung tumor motion identified in 4D CT scan images and compared it with the trajectory that was obtained by finite element simulation. The effects of rib kinematics on lung tumor trajectory were investigated. Over all phases of respiration, our developed model is able to predict the lung tumor motion with an average landmark error of [Formula: see text]. The results demonstrate the effectiveness of our physics-based model. We believe that this model can be potentially used in 4D dose computation, removal of breathing motion artifacts in positron emission tomography (PET) or gamma prompt image reconstruction.

Generation of a local lung respiratory motion model using a weighted sparse algorithm and motion prior-based registration

Respiration-introduced tumor location uncertainty is a challenge in lung percutaneous interventions, especially for the respiratory motion estimation of the tumor and surrounding vessel structures. In this work, a local motion modeling method is proposed based on whole-chest computed tomography (CT) and CT-fluoroscopy (CTF) scans. A weighted sparse statistical modeling (WSSM) method that can accurately capture location errors for each landmark point is proposed for lung motion prediction. By varying the sparse weight coefficients of the WSSM method, newly input motion information is approximately represented by a sparse linear combination of the respiratory motion repository and employed to serve as prior knowledge for the following registration process. We have also proposed an adaptive motion prior-based registration method to improve the motion prediction accuracy of the motion model in the region of interest (ROI). This registration method adopts a B-spline scheme to interactively weight the relative influence of the prior knowledge, model surface and image intensity information by locally controlling the deformation in the CTF image region. The proposed method has been evaluated on 15 image pairs between the end-expiratory (EE) and end-inspiratory (EI) phases and 31 four-dimensional CT (4DCT) datasets. The results reveal that the proposed WSSM method achieved a better motion prediction performance than other existing lung statistical motion modeling methods, and the motion prior-based registration method can generate more accurate local motion information in the ROI.