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

A model that predicts a real-time tumour surface using intra-treatment skin surface and end-of-expiration and end-of-inhalation planning CT images

OBJECTIVES

To develop a mapping model between skin surface motion and internal tumour motion and deformation using end-of-exhalation (EOE) and end-of-inhalation (EOI) 3D CT images for tracking lung tumours during respiration.

METHODS

Before treatment, skin and tumour surfaces were segmented and reconstructed from the EOE and the EOI 3D CT images. A non-rigid registration algorithm was used to register the EOE skin and tumour surfaces to the EOI, resulting in a displacement vector field that was then used to construct a mapping model. During treatment, the EOE skin surface was registered to the real-time, yielding a real-time skin surface displacement vector field. Using the mapping model generated, the input of a real-time skin surface can be used to calculate the real-time tumour surface. The proposed method was validated with and without simulated noise on 4D CT images from 15 patients at Léon Bérard Cancer Center and the 4D-lung dataset.

RESULTS

The average centre position error, dice similarity coefficient (DSC), 95%-Hausdorff distance and mean distance to agreement of the tumour surfaces were 1.29 mm, 0.924, 2.76 mm, and 1.13 mm without simulated noise, respectively. With simulated noise, these values were 1.33 mm, 0.920, 2.79 mm, and 1.15 mm, respectively.

CONCLUSIONS

A patient-specific model was proposed and validated that was constructed using only EOE and EOI 3D CT images and real-time skin surface images to predict internal tumour motion and deformation during respiratory motion.

ADVANCES IN KNOWLEDGE

The proposed method achieves comparable accuracy to state-of-the-art methods with fewer pre-treatment planning CT images, which holds potential for application in precise image-guided radiation therapy.

Current Research Status of Respiratory Motion for Thorax and Abdominal Treatment: A Systematic Review

Malignant tumors have become one of the serious public health problems in human safety and health, among which the chest and abdomen diseases account for the largest proportion. Early diagnosis and treatment can effectively improve the survival rate of patients. However, respiratory motion in the chest and abdomen can lead to uncertainty in the shape, volume, and location of the tumor, making treatment of the chest and abdomen difficult. Therefore, compensation for respiratory motion is very important in clinical treatment. The purpose of this review was to discuss the research and development of respiratory movement monitoring and prediction in thoracic and abdominal surgery, as well as introduce the current research status. The integration of modern respiratory motion compensation technology with advanced sensor detection technology, medical-image-guided therapy, and artificial intelligence technology is discussed and analyzed. The future research direction of intraoperative thoracic and abdominal respiratory motion compensation should be non-invasive, non-contact, use a low dose, and involve intelligent development. The complexity of the surgical environment, the constraints on the accuracy of existing image guidance devices, and the latency of data transmission are all present technical challenges.

A novel computer modeling and simulation technique for bronchi motion tracking in human lungs under respiration

In this work, we proposed a novel computer modeling and simulation technique for motion tracking of lung bronchi (or tumors) under respiration using 9 cases of computed tomography (CT)-based patient-specific finite element (FE) models and Ogden's hyperelastic model. In the fabrication of patient-specific FE models for the respiratory system, various organs such as the mediastinum, diaphragm, and thorax that could affect the lung motions during breathing were considered. To describe the nonlinear material behavior of lung parenchyma, the comparative simulation for biaxial tension-compression of lung parenchyma was carried out using several hyperelastic models in ABAQUS, and then, Ogden's model was adopted as an optimal model. Based on the aforementioned FE models and Ogden's material model, the 9 cases of respiration simulation were carried out from exhalation to inhalation, and the motion of lung bronchi (or tumors) was tracked. In addition, the changes in lung volume, lung cross-sectional area on the axial plane during breathing were calculated. Finally, the simulation results were quantitatively compared to the inhalation/exhalation CT images of 9 subjects to validate the proposed technique. Through the simulation, it was confirmed that the average relative errors of simulation to clinical data regarding to the displacement of 258 landmarks in the lung bronchi branches of total subjects were 1.10%~2.67%. In addition, the average relative errors of those with respect to the lung cross-sectional area changes and the volume changes in the superior-inferior direction were 0.20%~5.00% and 1.29 ~ 9.23%, respectively. Hence, it was considered that the simulation results were coincided well with the clinical data. The novelty of the present study is as follows: (1) The framework from fabrication of the human respiratory system to validation of the bronchi motion tracking is provided step by step. (2) The comparative simulation study for nonlinear material behavior of lung parenchyma was carried out to describe the realistic lung motion. (3) Various organs surrounding the lung parenchyma and restricting its motion were considered in respiration simulation. (4) The simulation results such as landmark displacement, lung cross-sectional area/volume changes were quantitatively compared to the clinical data of 9 subjects.

Extracting lung contour deformation features with deep learning for internal target motion tracking: a preliminary study

Objective. To propose lung contour deformation features (LCDFs) as a surrogate to estimate the thoracic internal target motion, and to report their performance by correlating with the changing body using a cascade ensemble model (CEM). LCDFs, correlated to the respiration driver, are employed without patient-specific motion data sampling and additional training before treatment.Approach. LCDFs are extracted by matching lung contours via an encoder-decoder deep learning model. CEM estimates LCDFs from the currently captured body, and then uses the estimated LCDFs to track internal target motion. The accuracy of the proposed LCDFs and CEM were evaluated using 48 targets' motion data, and compared with other published methods.Main results. LCDFs estimated the internal targets with a localization error of 2.6 ± 1.0 mm (average ± standard deviation). CEM reached a localization error of 4.7 ± 0.9 mm and a real-time performance of 256.9 ± 6.0 ms. With no internal anatomy knowledge, they achieved a small accuracy difference (of 0.34∼1.10 mm for LCDFs and of 0.43∼1.75 mm for CEM at the 95% confidence level) with a patient-specific lung biomechanical model and the deformable image registration models.Significance. The results demonstrated the effectiveness of LCDFs and CEM on tracking target motion. LCDFs and CEM are non-invasive, and require no patient-specific training before treatment. They show potential for broad applications.

Developing a Lung Model in the Age of COVID-19: A Digital Image Correlation and Inverse Finite Element Analysis Framework

Pulmonary diseases, driven by pollution, industrial farming, vaping, and the infamous COVID-19 pandemic, lead morbidity and mortality rates worldwide. Computational biomechanical models can enhance predictive capabilities to understand fundamental lung physiology; however, such investigations are hindered by the lung’s complex and hierarchical structure, and the lack of mechanical experiments linking the load-bearing organ-level response to local behaviors. In this study we address these impedances by introducing a novel reduced-order surface model of the lung, combining the response of the intricate bronchial network, parenchymal tissue, and visceral pleura. The inverse finite element analysis (IFEA) framework is developed using 3-D digital image correlation (DIC) from experimentally measured non-contact strains and displacements from an ex-vivo porcine lung specimen for the first time. A custom-designed inflation device is employed to uniquely correlate the multiscale classical pressure-volume bulk breathing measures to local-level deformation topologies and principal expansion directions. Optimal material parameters are found by minimizing the error between experimental and simulation-based lung surface displacement values, using both classes of gradient-based and gradient-free optimization algorithms and by developing an adjoint formulation for efficiency. The heterogeneous and anisotropic characteristics of pulmonary breathing are represented using various hyperelastic continuum formulations to divulge compound material parameters and evaluate the best performing model. While accounting for tissue anisotropy with fibers assumed along medial-lateral direction did not benefit model calibration, allowing for regional material heterogeneity enabled accurate reconstruction of lung deformations when compared to the homogeneous model. The proof-of-concept framework established here can be readily applied to investigate the impact of assorted organ-level ventilation strategies on local pulmonary force and strain distributions, and to further explore how diseased states may alter the load-bearing material behavior of the lung. In the age of a respiratory pandemic, advancing our understanding of lung biomechanics is more pressing than ever before.