1
|
Baschnagel AM, Flakus MJ, Wallat EM, Wuschner AE, Chappell RJ, Bayliss RA, Kimple RJ, Christensen GE, Reinhardt JM, Bassetti MF, Bayouth JE. A Phase 2 Randomized Clinical Trial Evaluating 4-Dimensional Computed Tomography Ventilation-Based Functional Lung Avoidance Radiation Therapy for Non-Small Cell Lung Cancer. Int J Radiat Oncol Biol Phys 2024; 119:1393-1402. [PMID: 38387810 DOI: 10.1016/j.ijrobp.2024.02.019] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2023] [Revised: 01/10/2024] [Accepted: 02/08/2024] [Indexed: 02/24/2024]
Abstract
PURPOSE To determine whether 4-dimensional computed tomography (4DCT) ventilation-based functional lung avoidance radiation therapy preserves pulmonary function compared with standard radiation therapy for non-small cell lung cancer (NSCLC). METHODS AND MATERIALS This single center, randomized, phase 2 trial enrolled patients with NSCLC receiving curative intent radiation therapy with either stereotactic body radiation therapy or conventionally fractionated radiation therapy between 2016 and 2022. Patients were randomized 1:1 to standard of care radiation therapy or functional lung avoidance radiation therapy. The primary endpoint was the change in Jacobian-based ventilation as measured on 4DCT from baseline to 3 months postradiation. Secondary endpoints included changes in volume of high- and low-ventilating lung, pulmonary toxicity, and changes in pulmonary function tests (PFTs). RESULTS A total of 122 patients were randomized and 116 were available for analysis. Median follow up was 29.9 months. Functional avoidance plans significantly (P < .05) reduced dose to high-functioning lung without compromising target coverage or organs at risk constraints. When analyzing all patients, there was no difference in the amount of lung showing a reduction in ventilation from baseline to 3 months between the 2 arms (1.91% vs 1.87%; P = .90). Overall grade ≥2 and grade ≥3 pulmonary toxicities for all patients were 24.1% and 8.6%, respectively. There was no significant difference in pulmonary toxicity or changes in PFTs between the 2 study arms. In the conventionally fractionated cohort, there was a lower rate of grade ≥2 pneumonitis (8.2% vs 32.3%; P = .049) and less of a decline in change in forced expiratory volume in 1 second (-3 vs -5; P = .042) and forced vital capacity (1.5 vs -6; P = .005) at 3 months, favoring the functional avoidance arm. CONCLUSIONS There was no difference in posttreatment ventilation as measured by 4DCT between the arms. In the cohort of patients treated with conventionally fractionated radiation therapy with functional lung avoidance, there was reduced pulmonary toxicity, and less decline in PFTs suggesting a clinical benefit in patients with locally advanced NSCLC.
Collapse
Affiliation(s)
- Andrew M Baschnagel
- Department of Human Oncology, University of Wisconsin Hospital and Clinics, Madison, Wisconsin.
| | - Mattison J Flakus
- Department of Medical Physics, University of Wisconsin-Madison, Madison, Wisconsin
| | - Eric M Wallat
- Department of Medical Physics, University of Wisconsin-Madison, Madison, Wisconsin
| | - Antonia E Wuschner
- Department of Medical Physics, University of Wisconsin-Madison, Madison, Wisconsin
| | - Richard J Chappell
- Department of Biostatistics and Medical Informatics, University of Wisconsin-Madison, Madison, Wisconsin
| | - R Adam Bayliss
- Department of Human Oncology, University of Wisconsin Hospital and Clinics, Madison, Wisconsin
| | - Randall J Kimple
- Department of Human Oncology, University of Wisconsin Hospital and Clinics, Madison, Wisconsin
| | - Gary E Christensen
- Department of Electrical and Computer Engineering, University of Iowa, Iowa City, Iowa; Department of Radiation Oncology, University of Iowa, Iowa City, Iowa
| | - Joseph M Reinhardt
- Roy J. Carver Department of Biomedical Engineering, University of Iowa, Iowa City, Iowa
| | - Michael F Bassetti
- Department of Human Oncology, University of Wisconsin Hospital and Clinics, Madison, Wisconsin
| | - John E Bayouth
- Department of Radiation Medicine, Oregon Health & Science University, Portland, Oregon.
| |
Collapse
|
2
|
Galloy AE, Reinhardt JM, Raghavan ML. Role of lung lobar sliding on parenchymal distortion during breathing. J Appl Physiol (1985) 2023; 135:534-541. [PMID: 37439240 PMCID: PMC10538991 DOI: 10.1152/japplphysiol.00631.2022] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2022] [Revised: 06/14/2023] [Accepted: 07/10/2023] [Indexed: 07/14/2023] Open
Abstract
Sliding between lung lobes along lobar fissures is a poorly understood aspect of lung mechanics. The objective of this study was to test the hypothesis that lobar sliding helps reduce distortion in the lung parenchyma during breathing. Finite element models of left lungs with geometries and boundary conditions derived from medical images of human subjects were developed. Effect of lobar sliding was studied by comparing nonlinear finite elastic contact mechanics simulations that allowed and disallowed lobar sliding. Lung parenchymal distortion during simulated breath-holds and tidal breathing was quantified with the model's spatial mean anisotropic deformation index (ADI), a measure of directional preference in volume change that varies spatially in the lung. Models that allowed lobar sliding had significantly lower mean ADI (i.e., lesser parenchymal distortion) than models that disallowed lobar sliding under simulations of both tidal breathing (5.3% median difference, P = 0.008, n = 8) and lung deformation between breath-holds at total lung capacity and functional residual capacity (3.2% median difference, P = 0.03, n = 6). This effect was most pronounced in the lower lobe where lobar sliding reduced parenchymal distortion with statistical significance, but not in the upper lobe. In addition, more lobar sliding was correlated with greater reduction in distortion between sliding and nonsliding models in our study cohorts (Pearson's correlation coefficient of 0.95 for tidal breathing, 0.87 for breath-holds, and 0.91 for the combined dataset). These findings are consistent with the hypothesis that lung lobar sliding reduces parenchymal distortion during breathing.NEW & NOTEWORTHY The role of lobar sliding in lung mechanics is poorly understood. Delineating this role could help explain how breathing is affected by anatomical differences between subjects such as incomplete and missing lobar fissures. We used computational contact mechanics models of lungs from human subjects to delineate the effect of lobar sliding by comparing simulations that allowed and disallowed sliding. We found evidence consistent with the hypothesis that lung lobar sliding reduces parenchymal distortion during breathing.
Collapse
Affiliation(s)
- Adam E Galloy
- Roy J. Carver Department of Biomedical Engineering, University of Iowa, Iowa City, Iowa, United States
| | - Joseph M Reinhardt
- Roy J. Carver Department of Biomedical Engineering, University of Iowa, Iowa City, Iowa, United States
| | - Madhavan L Raghavan
- Roy J. Carver Department of Biomedical Engineering, University of Iowa, Iowa City, Iowa, United States
| |
Collapse
|
3
|
Huang YH, Teng X, Zhang J, Chen Z, Ma Z, Ren G, Kong FMS, Ge H, Cai J. Respiratory Invariant Textures From Static Computed Tomography Scans for Explainable Lung Function Characterization. J Thorac Imaging 2023; 38:286-296. [PMID: 37265243 DOI: 10.1097/rti.0000000000000717] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/03/2023]
Abstract
PURPOSE The inherent characteristics of lung tissue independent of breathing maneuvers may provide fundamental information for function assessment. This paper attempted to correlate textural signatures from computed tomography (CT) with pulmonary function measurements. MATERIALS AND METHODS Twenty-one lung cancer patients with thoracic 4-dimensional CT, DTPA-single-photon emission CT ventilation ( VNM ) scans, and available spirometry measurements (forced expiratory volume in 1 s, FEV 1 ; forced vital capacity, FVC; and FEV 1 /FVC) were collected. In subregional feature discovery, function-correlated candidates were identified from 79 radiomic features based on the statistical strength to differentiate defected/nondefected lung regions. Feature maps (FMs) of selected candidates were generated on 4-dimensional CT phases for a voxel-wise feature distribution study. Quantitative metrics were applied for validations, including the Spearman correlation coefficient (SCC) and the Dice similarity coefficient for FM- VNM spatial agreement assessments, intraclass correlation coefficient for FM interphase robustness evaluations, and FM-spirometry comparisons. RESULTS At the subregion level, 8 function-correlated features were identified (effect size>0.330). The FMs of candidates yielded moderate-to-strong voxel-wise correlations with the reference VNM . The FMs of gray level dependence matrix dependence nonuniformity showed the highest robust (intraclass correlation coefficient=0.96 and P <0.0001) spatial correlation, with median SCCs ranging from 0.54 to 0.59 throughout the 10 breathing phases. Its phase-averaged FM achieved a median SCC of 0.60, a median Dice similarity coefficient of 0.60 (0.65) for high (low) functional lung volumes, and a correlation of 0.565 (0.646) between the spatially averaged feature values and FEV 1 (FEV 1 /FVC). CONCLUSIONS The results provide further insight into the underlying association of specific pulmonary textures with both local ( VNM ) and global (FEV 1 /FVC, FEV 1 ) functions. Further validations of the FM generalizability and the standardization of implementation protocols are warranted before clinically relevant investigations.
Collapse
Affiliation(s)
- Yu-Hua Huang
- Department of Health Technology and Informatics, The Hong Kong Polytechnic University
| | - Xinzhi Teng
- Department of Health Technology and Informatics, The Hong Kong Polytechnic University
| | - Jiang Zhang
- Department of Health Technology and Informatics, The Hong Kong Polytechnic University
| | - Zhi Chen
- Department of Health Technology and Informatics, The Hong Kong Polytechnic University
| | - Zongrui Ma
- Department of Health Technology and Informatics, The Hong Kong Polytechnic University
| | - Ge Ren
- Department of Health Technology and Informatics, The Hong Kong Polytechnic University
| | - Feng-Ming Spring Kong
- Department of Clinical Oncology, Li Ka Shing Faculty of Medicine, The University of Hong Kong, Hong Kong SAR
- Department of Clinical Oncology, The University of Hong Kong-Shenzhen Hospital, Shenzhen
| | - Hong Ge
- Department of Radiation Oncology, The Affiliated Cancer Hospital of Zhengzhou University & Henan Cancer Hospital, Zhengzhou, China
| | - Jing Cai
- Department of Health Technology and Informatics, The Hong Kong Polytechnic University
| |
Collapse
|
4
|
Huang YH, Ren G, Xiao H, Yang D, Kong FMS, Ho WY, Cai J. Volumetric multiphase ventilation imaging based on four-dimensional computed tomography for functional lung avoidance radiotherapy. Med Phys 2022; 49:7237-7246. [PMID: 35841346 DOI: 10.1002/mp.15847] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2021] [Revised: 04/20/2022] [Accepted: 06/23/2022] [Indexed: 12/13/2022] Open
Abstract
PURPOSE Current computed tomography (CT)-based lung ventilation imaging (CTVI) techniques derive a static ventilation image without temporal information. This research aims to develop a four-dimensional CT (4DCT)-based multiphase dynamic ventilation imaging framework capable of recovering the entire ventilation process throughout the breathing cycle for functional lung avoidance radiotherapy (FLART). METHODS A total of 15 free-breathing thoracic 4DCT scans of lung or esophageal cancer patients were collected from the public datasets. The lung region of each phase image was first delineated, and then the mask-free isotropic total variation image registration algorithm was used to derive the deformation vector fields between the end-expiration (EE) phase and other phases. As a surrogate of ventilation, the voxel-wise local expansion ratio of each phase relative to the EE phase was estimated using the parameterized Integrated Jacobian Formulation method in the EE phase coordinate. Lastly, the dynamic ventilation images were generated by warping these phase-specific local expansion distributions with a same geometry into their respective breathing phases. Quantitative analysis, including interphase Spearman correlation coefficients, voxel-wise, and regional-wise expansion/contraction tracking, were performed to indirectly validate the proposed method. RESULTS The proposed method maintains the physiological meaning of ventilation on each phase and enables to recover the dynamic lung ventilation process. The mean interphase Spearman correlations ranged between 0.23 ± 0.20 and 0.93 ± 0.04 and decreased near the EE phase. Only 26.2% (2.59E + 6 out of 9.89E + 6) of lung voxels exhibited the same expansion/contraction pattern as the global lung. Qualitative and quantitative evaluations of the interphase ventilation distribution difference show that ventilation spatiotemporal heterogeneities generally exist during respiration. CONCLUSIONS In contrast to conventional CTVI metrics, our method enables to extract additional phase-resolved respiration-correlated information and reflects the generally existed ventilation spatiotemporal heterogeneities. Subsequent studies with quantitative phase-by-phase cross-modality evaluations will further explore its potential to deepen our understanding of lung function and respiration mechanics and also to facilitate more accurate implementation of FLART.
Collapse
Affiliation(s)
- Yu-Hua Huang
- Department of Health Technology and Informatics, The Hong Kong Polytechnic University, Kowloon, Hong Kong SAR
| | - Ge Ren
- Department of Health Technology and Informatics, The Hong Kong Polytechnic University, Kowloon, Hong Kong SAR
| | - Haonan Xiao
- Department of Health Technology and Informatics, The Hong Kong Polytechnic University, Kowloon, Hong Kong SAR
| | - Dongrong Yang
- Department of Health Technology and Informatics, The Hong Kong Polytechnic University, Kowloon, Hong Kong SAR
| | - Feng-Ming Spring Kong
- Department of Clinical Oncology, The University of Hong Kong, Hong Kong, Hong Kong SAR
| | - Wai Yin Ho
- Department of Nuclear Medicine, Queen Mary Hospital, Hong Kong, Hong Kong SAR
| | - Jing Cai
- Department of Health Technology and Informatics, The Hong Kong Polytechnic University, Kowloon, Hong Kong SAR
| |
Collapse
|
5
|
Inhaled Gas Magnetic Resonance Imaging: Advances, Applications, Limitations, and New Frontiers. Mol Imaging 2021. [DOI: 10.1016/b978-0-12-816386-3.00013-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022] Open
|
6
|
Liang X, Yin FF, Wang C, Cai J. A robust deformable image registration enhancement method based on radial basis function. Quant Imaging Med Surg 2019; 9:1315-1325. [PMID: 31448216 DOI: 10.21037/qims.2019.07.05] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
Background To develop and evaluate a robust deformable image registration (DIR) enhancement method based on radial basis function (RBF) expansion. Methods To improve DIR accuracy using sparsely available measured displacements, it is crucial to estimate the motion correlation between the voxels. In the proposed method, we chose to derive this correlation from the initial displacement vector fields (DVFs), and represent it in the form of RBF expansion coefficients of the voxels. The method consists of three steps: (I) convert an initial DVF to a coefficient matrix comprising expansion coefficients of the Wendland's RBF; (II) modify the coefficient matrix under the guidance of sparely distributed landmarks to generate the post-enhancement coefficient matrix; and (III) convert the post-enhancement coefficient matrix to the post-enhancement DVF. The method was tested on five DIR algorithms using a digital phantom. 3D registration errors were calculated for comparisons between the pre-/post-enhancement DVFs and the ground-truth DVFs. Effects of the number and locations of landmarks on DIR enhancement were evaluated. Results After applying the DIR enhancement method, the 3D registration errors per voxel (unit: mm) were reduced from pre-enhancement to post-enhancement by 1.3 (2.4 to 1.1, 54.2%), 0.0 (0.9 to 0.9, 0.0%), 6.1 (8.2 to 2.1, 74.4%), 3.2 (4.7 to 1.5, 68.1%), and 1.7 (2.9 to 1.2, 58.6%) for the five tested DIR algorithms respectively. The average DIR error reduction was 2.5±2.3 mm (percentage error reduction: 51.1%±29.1%). 3D registration errors decreased inverse-exponentially as the number of landmarks increased, and were insensitive to the landmarks' locations in relation to the down-sampling DVF grids. Conclusions We demonstrated the feasibility of a robust RBF-based method for enhancing DIR accuracy using sparsely distributed landmarks. This method has been shown robust and effective in reducing DVF errors using different numbers and distributions of landmarks for various DIR algorithms.
Collapse
Affiliation(s)
- Xiao Liang
- Medical Physics Graduate Program, Duke University, Durham, NC, USA
| | - Fang-Fang Yin
- Medical Physics Graduate Program, Duke University, Durham, NC, USA.,Department of Radiation Oncology, Duke University Medical Center, Durham, NC, USA
| | - Chunhao Wang
- Department of Radiation Oncology, Duke University Medical Center, Durham, NC, USA
| | - Jing Cai
- Department of Radiation Oncology, Duke University Medical Center, Durham, NC, USA.,Department of Health Technology and Informatics, The Hong Kong Polytechnic University, Hong Kong, China
| |
Collapse
|
7
|
Hu L, Huang Q, Cui T, Duarte I, Miller GW, Mugler JP, Cates GD, Mata JF, de Lange EE, Altes TA, Yin FF, Cai J. A hybrid proton and hyperpolarized gas tagging MRI technique for lung respiratory motion imaging: a feasibility study. Phys Med Biol 2019; 64:105019. [PMID: 30947154 DOI: 10.1088/1361-6560/ab160c] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
The aim of this work was to develop a novel hybrid 3D hyperpolarized (HP) gas tagging MRI (t-MRI) technique and to evaluate it for lung respiratory motion measurement with comparison to deformable image registrations (DIR) methods. Three healthy subjects underwent a hybrid MRI which combines 3D HP gas t-MRI with a low resolution (Low-R, 4.5 mm isotropic voxels) 3D proton MRI (p-MRI), plus a high resolution (High-R, 2.5 mm isotropic voxels) 3D p-MRI, during breath-holds at the end-of-inhalation (EOI) and the end-of-exhalation (EOE). Displacement vector field (DVF) of the lung motion was determined from the t-MRI images by tracking tagging grids and from the High-R p-MRI using three DIR methods (B-spline based method implemented by Velocity, Free Form Deformation by MIM, and B-spline by an open source software Elastix: denoted as A, B, and C, respectively), labeled as tDVF and dDVF, respectively. The tDVF from the HP gas t-MRI was used as ground-truth reference to evaluate performance of the three DIR methods. Differences in both magnitude and angle between the tDVF and dDVFs were analyzed. The mean lung motion of the three subjects was 37.3 mm, 8.9 mm and 12.9 mm, respectively. Relatively large discrepancies were observed between the tDVF and the dDVFs as compared to previously reported DIR errors. The mean ± standard deviation (SD) DVF magnitude difference was 8.3 ± 5.6 mm, 9.2 ± 4.5 mm, and 9.3 ± 6.1 mm, and the mean ± SD DVF angular difference was 29.1 ± 12.1°, 50.1 ± 28.6°, and 39.0 ± 6.3°, for the DIR Methods A, B, and C, respectively. These preliminary results showed that the hybrid HP gas t-MRI technique revealed different lung motion patterns as compared to the DIR methods. It may provide unique perspectives in developing and evaluating DIR of the lungs. Novelty and Significance We designed a MRI protocol that includes a novel hybrid MRI technique (3D HP gas t-MRI with a low resolution 3D p-MRI) plus a high resolution 3D p-MRI. We tested the novel hybrid MRI technique on three healthy subjects for measuring regional lung respiratory motion with comparison to deformable image registrations (DIR) methods, and observed relatively large discrepancies in lung motion between HP gas t-MRI and DIR methods.
Collapse
Affiliation(s)
- Lei Hu
- Department of Radiation Oncology, NYU Langone Health, New York, NY 10016, United States of America
| | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
8
|
Cui T, Miller GW, Mugler JP, Cates GD, Mata JF, de Lange EE, Huang Q, Altes TA, Yin FF, Cai J. An initial investigation of hyperpolarized gas tagging magnetic resonance imaging in evaluating deformable image registration-based lung ventilation. Med Phys 2018; 45:5535-5542. [PMID: 30276819 DOI: 10.1002/mp.13223] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2018] [Revised: 08/21/2018] [Accepted: 09/19/2018] [Indexed: 11/11/2022] Open
Abstract
BACKGROUND Deformable image registration (DIR)-based lung ventilation mapping is attractive due to its simplicity, and also challenging due to its susceptibility to errors and uncertainties. In this study, we explored the use of 3D Hyperpolarized (HP) gas tagging MRI to evaluate DIR-based lung ventilation. METHOD AND MATERIAL Three healthy volunteers included in this study underwent both 3D HP gas tagging MRI (t-MRI) and 3D proton MRI (p-MRI) using balanced steady-state free precession pulse sequence at end of inhalation and end of exhalation. We first obtained the reference displacement vector fields (DVFs) from the t-MRIs by tracking the motion of each tagging grid between the exhalation and the inhalation phases. Then, we determined DIR-based DVFs from the p-MRIs by registering the images at the two phases with two commercial DIR algorithms. Lung ventilations were calculated from both the reference DVFs and the DIR-based DVFs using the Jacobian method and then compared using cross correlation and mutual information. RESULTS The DIR-based lung ventilations calculated using p-MRI varied considerably from the reference lung ventilations based on t-MRI among all three subjects. The lung ventilations generated using Velocity AI were preferable for the better spatial homogeneity and accuracy compared to the ones using MIM, with higher average cross correlation (0.328 vs 0.262) and larger average mutual information (0.528 vs 0.323). CONCLUSION We demonstrated that different DIR algorithms resulted in different lung ventilation maps due to underlining differences in the DVFs. HP gas tagging MRI provides a unique platform for evaluating DIR-based lung ventilation.
Collapse
Affiliation(s)
- Taoran Cui
- Department of Radiation Oncology, Rutgers Cancer Institute of New Jersey, New Brunswick, NJ, 08901, USA
| | - G Wilson Miller
- Department of Radiology and Medical Imaging, University of Virginia, Charlottesville, VA, 22908, USA
| | - John P Mugler
- Department of Radiology and Medical Imaging, University of Virginia, Charlottesville, VA, 22908, USA
| | - Gordon D Cates
- Department of Radiology and Medical Imaging, University of Virginia, Charlottesville, VA, 22908, USA
| | - Jaime F Mata
- Department of Radiology and Medical Imaging, University of Virginia, Charlottesville, VA, 22908, USA
| | - Eduard E de Lange
- Department of Radiology and Medical Imaging, University of Virginia, Charlottesville, VA, 22908, USA
| | - Qijie Huang
- Department of Department of Medical Physics, Memorial Sloan-Kettering Cancer Center, New York, NY, 10065, USA
| | - Talissa A Altes
- Department of Radiology, University of Missouri School of Medicine, Columbia, Missouri, 65212, USA
| | - Fang-Fang Yin
- Department of Radiation Oncology, Duke University Medical Center, Durham, NC, 27710, USA
| | - Jing Cai
- Department of Radiation Oncology, Duke University Medical Center, Durham, NC, 27710, USA.,Department of Health Technology and Informatics, Hong Kong Polytechnic University, Hong Kong, China
| |
Collapse
|
9
|
Nasehi Tehrani J, McEwan A, Wang J. Lung surface deformation prediction from spirometry measurement and chest wall surface motion. Med Phys 2016; 43:5493. [PMID: 27782714 PMCID: PMC5035308 DOI: 10.1118/1.4962479] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2016] [Revised: 08/26/2016] [Accepted: 08/29/2016] [Indexed: 11/07/2022] Open
Abstract
PURPOSE The authors have developed and evaluated a method to predict lung surface motion based on spirometry measurements, and chest and abdomen motion at selected locations. METHODS A patient-specific 3D triangular surface mesh of the lung region was obtained at the end expiratory phase by the threshold-based segmentation method. Lung flow volume changes were recorded with a spirometer for each patient. A total of 192 selected points at a regular spacing of 2 × 2 cm matrix points were used to detect chest wall motion over a total area of 32 × 24 cm covering the chest and abdomen surfaces. QR factorization with column pivoting was employed to remove redundant observations of the chest and abdominal areas. To create a statistical model between the lung surface and the corresponding surrogate signals, the authors developed a predictive model based on canonical ridge regression. Two unique weighting vectors were selected for each vertex on the lung surface; they were optimized during the training process using all other 4D-CT phases except for the test inspiration phase. These parameters were employed to predict the vertex locations of a testing data set. RESULTS The position of each lung surface mesh vertex was estimated from the motion at selected positions within the chest wall surface and from spirometry measurements in ten lung cancer patients. The average estimation of the 98th error percentile for the end inspiration phase was less than 1 mm (AP = 0.9 mm, RL = 0.6 mm, and SI = 0.8 mm). The vertices located at the lower region of the lung had a larger estimation error as compared with those within the upper region of the lung. The average landmark motion errors, derived from the biomechanical modeling using real surface deformation vector fields (SDVFs), and the predicted SDVFs were 3.0 and 3.1 mm, respectively. CONCLUSIONS Our newly developed predictive model provides a noninvasive approach to derive lung boundary conditions. The proposed system can be used with personalized biomechanical respiration modeling to derive lung tumor motion during radiation therapy from noninvasive measurements.
Collapse
Affiliation(s)
- Joubin Nasehi Tehrani
- Department of Radiation Oncology, UT Southwestern Medical Center, Dallas, Texas 75235-8808
| | - Alistair McEwan
- School of Electrical and Information Engineering, University of Sydney, New South Wales 2006, Australia
| | - Jing Wang
- Department of Radiation Oncology, UT Southwestern Medical Center, Dallas, Texas 75235-8808
| |
Collapse
|
10
|
Ruppert K. Biomedical imaging with hyperpolarized noble gases. REPORTS ON PROGRESS IN PHYSICS. PHYSICAL SOCIETY (GREAT BRITAIN) 2014; 77:116701. [PMID: 25360484 DOI: 10.1088/0034-4885/77/11/116701] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
Hyperpolarized noble gases (HNGs), polarized to approximately 50% or higher, have led to major advances in magnetic resonance (MR) imaging of porous structures and air-filled cavities in human subjects, particularly the lung. By boosting the available signal to a level about 100 000 times higher than that at thermal equilibrium, air spaces that would otherwise appear as signal voids in an MR image can be revealed for structural and functional assessments. This review discusses how HNG MR imaging differs from conventional proton MR imaging, how MR pulse sequence design is affected and how the properties of gas imaging can be exploited to obtain hitherto inaccessible information in humans and animals. Current and possible future imaging techniques, and their application in the assessment of normal lung function as well as certain lung diseases, are described.
Collapse
|
11
|
Liu Z, Araki T, Okajima Y, Albert M, Hatabu H. Pulmonary hyperpolarized noble gas MRI: Recent advances and perspectives in clinical application. Eur J Radiol 2014; 83:1282-1291. [DOI: 10.1016/j.ejrad.2014.04.014] [Citation(s) in RCA: 37] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2013] [Revised: 02/21/2014] [Accepted: 04/19/2014] [Indexed: 12/01/2022]
|
12
|
Amelon RE, Cao K, Reinhardt JM, Christensen GE, Raghavan ML. A measure for characterizing sliding on lung boundaries. Ann Biomed Eng 2013; 42:642-50. [PMID: 24114112 DOI: 10.1007/s10439-013-0920-5] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2013] [Accepted: 09/30/2013] [Indexed: 11/28/2022]
Abstract
The lobes of the lung slide relative to each other during breathing. Quantifying lobar sliding can aid in better understanding lung function, better modeling of lung dynamics, and for studying phenomenon such as pleural adhesion. We propose a novel measure to characterize lobe sliding in the lung based on the displacement field obtained from image registration of CT scans. When two sliding lobes are modeled as a continuum, the discontinuity in the displacement field at the fissure will manifest as elevated maximum shear--the proposed measure--which is capable of capturing both the level and orientation of sliding. Six human lungs were analyzed using scans spanning functional residual capacity to total lung capacity. The lung lobes were segmented and registered on a lobe-by-lobe basis to obtain the displacement field from which the proposed sliding measure was calculated. The sliding measure was found to be insignificant in the parenchyma, as relatively little tissue shear occurs here. On the other hand, it was elevated along the fissures. Thus, a map of the proposed sliding measure of the entire lung clearly delineates and quantifies sliding between lung lobes. Sliding is a key aspect of lung deformation during breathing. The proposed measure may help resolve artifacts introduced by sliding in deformation analysis techniques used for radiotherapy.
Collapse
Affiliation(s)
- Ryan E Amelon
- Department of Biomedical Engineering, University of Iowa, 1136 Seamans Center, Iowa City, IA, 52242, USA
| | | | | | | | | |
Collapse
|
13
|
Chan DD, Toribio D, Neu CP. Displacement smoothing for the precise MRI-based measurement of strain in soft biological tissues. Comput Methods Biomech Biomed Engin 2012; 16:852-60. [PMID: 22292493 DOI: 10.1080/10255842.2011.641178] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/14/2022]
Abstract
Displacement and strain are fundamental quantities that describe the normal and pathological mechanical function of soft biological materials. Non-invasive imaging techniques, including displacement-encoded magnetic resonance imaging (MRI), enable the direct calculation of biomaterial displacements during the application of extrinsic mechanical forces. However, because strain is derived from measured MRI-based displacements, data processing must be accomplished to minimise the propagation and amplification of errors. Here, we evaluate smoothing methods (including averaging filters, splines, finite impulse response filters and wavelets) that enable the calculation of strain in biomaterials from MRI-based displacements for minimal error, defined in terms of bias and precision. Displacement and strain precisions were improved using all smoothing methods studied. Precision generally increased with the number of smoothing iterations (i.e. repeated applications) of a chosen smoothing method. The bias depended on the smoothing method and tended to increase with the number of smoothing iterations. A Gaussian filter characterised complex and heterogeneous strain fields with maximum precision and minimum bias. The results suggest that the optimal choice of smoothing method to compute strain for a given biomaterial or tissue application depends on a careful consideration of trade-offs between the improved precision (with increased data smoothing) and the trending increase in bias.
Collapse
Affiliation(s)
- Deva D Chan
- a Weldon School of Biomedical Engineering, Purdue University , West Lafayette , IN 47907 , USA
| | | | | |
Collapse
|
14
|
Tustison NJ, Cook TS, Song G, Gee JC. Pulmonary kinematics from image data: a review. Acad Radiol 2011; 18:402-17. [PMID: 21377592 DOI: 10.1016/j.acra.2010.10.019] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2010] [Revised: 09/02/2010] [Accepted: 10/25/2010] [Indexed: 10/18/2022]
Abstract
The effects of certain lung pathologies include alterations in lung physiology negatively affecting pulmonary compliance. Current approaches to diagnosis and treatment assessment of lung disease commonly rely on pulmonary function testing. Such testing is limited to global measures of lung function, neglecting regional measurements, which are critical for early diagnosis and localization of disease. Increased accessibility to medical image acquisition strategies with high spatiotemporal resolution coupled with the development of sophisticated intensity-based and geometric registration techniques has resulted in the recent exploration of modeling pulmonary motion for calculating local measures of deformation. In this review, the authors provide a broad overview of such research efforts for the estimation of pulmonary deformation. This includes discussion of various techniques, current trends in validation approaches, and the public availability of software and data resources.
Collapse
|
15
|
Fain S, Schiebler ML, McCormack DG, Parraga G. Imaging of lung function using hyperpolarized helium-3 magnetic resonance imaging: Review of current and emerging translational methods and applications. J Magn Reson Imaging 2011; 32:1398-408. [PMID: 21105144 DOI: 10.1002/jmri.22375] [Citation(s) in RCA: 163] [Impact Index Per Article: 12.5] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022] Open
Abstract
During the past several years there has been extensive development and application of hyperpolarized helium-3 (HP (3)He) magnetic resonance imaging (MRI) in clinical respiratory indications such as asthma, chronic obstructive pulmonary disease, cystic fibrosis, radiation-induced lung injury, and transplantation. This review focuses on the state-of-the-art of HP (3)He MRI and its application to clinical pulmonary research. This is not an overview of the physics of the method, as this topic has been covered previously. We focus here on the potential of this imaging method and its challenges in demonstrating new types of information that has the potential to influence clinical research and decision making in pulmonary medicine. Particular attention is given to functional imaging approaches related to ventilation and diffusion-weighted imaging with applications in chronic obstructive pulmonary disease, cystic fibrosis, asthma, and radiation-induced lung injury. The strengths and challenges of the application of (3)He MRI in these indications are discussed along with a comparison to established and emerging imaging techniques.
Collapse
Affiliation(s)
- Sean Fain
- Department of Medical Physics, University of Wisconsin-Madison, Madison, Wisconsin, USA
| | | | | | | |
Collapse
|
16
|
Cai J, McLawhorn R, Altes TA, de Lange E, Read PW, Larner JM, Benedict SH, Sheng K. Helical tomotherapy planning for lung cancer based on ventilation magnetic resonance imaging. Med Dosim 2011; 36:389-96. [PMID: 21377866 DOI: 10.1016/j.meddos.2010.09.008] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2010] [Revised: 08/30/2010] [Accepted: 09/21/2010] [Indexed: 12/31/2022]
Abstract
To investigate the feasibility of lung ventilation-based treatment planning, computed tomography and hyperpolarized (HP) helium-3 (He-3) magnetic resonance imaging (MRI) ventilation images of 6 subjects were coregistered for intensity-modulated radiation therapy planning in Tomotherapy. Highly-functional lungs (HFL) and less-functional lungs (LFL) were contoured based on their ventilation image intensities, and a cylindrical planning-target-volume was simulated at locations adjacent to both HFL and LFL. Annals of an anatomy-based plan (Plan 1) and a ventilation-based plan (Plan 2) were generated. The following dosimetric parameters were determined and compared between the 2 plans: percentage of total/HFL volume receiving ≥20 Gy, 15 Gy, 10 Gy, and 5 Gy (TLV(20), HFLV(20), TLV(15), HFLV(15), TLV(10), HFLV(10), TLV(5), HFLV(5)), mean total/HFL dose (MTLD/HFLD), maximum doses to all organs at risk (OARs), and target dose conformality. Compared with Plan 1, Plan 2 reduced mean HFLD (mean reduction, 0.8 Gy), MTLD (mean reduction, 0.6 Gy), HFLV(20) (mean reduction, 1.9%), TLV(20) (mean reduction, 1.5%), TLV(15) (mean reduction, 1.7%), and TLV(10) (mean reduction, 2.1%). P-values of the above comparisons are less than 0.05 using the Wilcoxon signed rank test. For HFLV(15), HFLV(10), TLV(5), and HTLV(5), Plan 2 resulted in lower values than plan 1 but the differences are not significant (P-value range, 0.063-0.219). Plan 2 did not significantly change maximum doses to OARs (P-value range, 0.063-0.563) and target conformality (P = 1.000). HP He-3 MRI of patients with lung disease shows a highly heterogeneous ventilation capacity that can be utilized for functional treatment planning. Moderate but statistically significant improvements in sparing functional lungs were achieved using helical tomotherapy plans.
Collapse
Affiliation(s)
- Jing Cai
- Department of Radiation Oncology, University of Virginia, Charlottesville, VA 22908, USA
| | | | | | | | | | | | | | | |
Collapse
|
17
|
Tustison NJ, Awate SP, Cai J, Altes TA, Miller GW, de Lange EE, Mugler JP, Gee JC. Pulmonary kinematics from tagged hyperpolarized helium-3 MRI. J Magn Reson Imaging 2010; 31:1236-41. [PMID: 20432362 DOI: 10.1002/jmri.22137] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
PURPOSE To propose and test the feasibility of a novel method for quantifying 3D regional pulmonary kinematics from hyperpolarized helium-3 tagged MRI in human subjects using a tailored image processing pipeline and a recently developed nonrigid registration framework. MATERIALS AND METHODS Following image acquisition, inspiratory and expiratory tagged (3)He magnetic resonance (MR) images were preprocessed using various image filtering techniques to enhance the tag surfaces. Segmentation of the three orthogonal sets of tag planes in each lung produced distinct point-set representations of the tag surfaces. Using these labeled point-sets, deformation fields and corresponding strain maps were obtained via nonrigid point-set registration. Kinematic analysis was performed on three volunteers. RESULTS Tag lines in inspiratory and expiratory images were coregistered producing a continuous 3D correspondence mapping. Average displacement and directional strains were calculated in three subjects in the inferior, mid, and superior portions of the right and left lungs. As expected, the predominant direction of displacements with expiration is from inferior to superior. CONCLUSION Kinematic quantitation of pulmonary motion using tagged (3)He MRI is feasible using the applied image preprocessing filtering techniques and nonrigid point-set registration. Potential benefits from regional pulmonary kinematic quantitation include the facilitation of diagnosis and local assessment of disease progression.
Collapse
Affiliation(s)
- Nicholas J Tustison
- Penn Image Computing and Science Laboratory, University of Pennsylvania, Philadelphia, PA 19104, USA.
| | | | | | | | | | | | | | | |
Collapse
|
18
|
Cai J, Sheng K, Benedict SH, Read PW, Larner JM, Mugler JP, de Lange EE, Cates GD, Miller GW. Dynamic MRI of grid-tagged hyperpolarized helium-3 for the assessment of lung motion during breathing. Int J Radiat Oncol Biol Phys 2009; 75:276-84. [PMID: 19540059 DOI: 10.1016/j.ijrobp.2009.03.051] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/26/2008] [Revised: 03/10/2009] [Accepted: 03/19/2009] [Indexed: 12/25/2022]
Abstract
PURPOSE To develop a dynamic magnetic resonance imaging (MRI) tagging technique using hyperpolarized helium-3 (HP He-3) to track lung motion. METHODS AND MATERIALS An accelerated non-Cartesian k-space trajectory was used to gain acquisition speed, at the cost of introducing image artifacts, providing a viable strategy for obtaining whole-lung coverage with adequate temporal resolution. Multiple-slice two-dimensional dynamic images of the lung were obtained in three healthy subjects after inhaling He-3 gas polarized to 35%-40%. Displacement, strain, and ventilation maps were computed from the observed motion of the grid peaks. RESULTS Both temporal and spatial variations of pulmonary mechanics were observed in normal subjects, including shear motion between different lobes of the same lung. CONCLUSION These initial results suggest that dynamic imaging of grid-tagged hyperpolarized magnetization may potentially be a powerful tool for observing and quantifying pulmonary biomechanics on a regional basis and for assessing, validating, and improving lung deformable image registration algorithms.
Collapse
Affiliation(s)
- Jing Cai
- Department of Radiation Oncology, University of Virginia, Charlottesville, VA 22908, USA
| | | | | | | | | | | | | | | | | |
Collapse
|
19
|
Cai J, Read PW, Larner JM, Jones DR, Benedict SH, Sheng K. Reproducibility of interfraction lung motion probability distribution function using dynamic MRI: statistical analysis. Int J Radiat Oncol Biol Phys 2008; 72:1228-35. [PMID: 18954717 DOI: 10.1016/j.ijrobp.2008.07.028] [Citation(s) in RCA: 27] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2007] [Revised: 06/09/2008] [Accepted: 07/10/2008] [Indexed: 10/21/2022]
Abstract
PURPOSE To investigate the statistical reproducibility of craniocaudal probability distribution function (PDF) of interfraction lung motion using dynamic magnetic resonance imaging. METHODS AND MATERIALS A total of 17 subjects, 9 healthy volunteers and 8 lung tumor patients, underwent two to three continuous 300-s magnetic resonance imaging scans in the sagittal plane, repeated 2 weeks apart. Three pulmonary vessels from different lung regions (upper, middle, and lower) in the healthy subjects and lung tumor patients were selected for tracking, and the displacement PDF reproducibility was evaluated as a function of scan time and frame rate. RESULTS For both healthy subjects and patients, the PDF reproducibility improved with increased scan time and converged to an equilibrium state during the 300-s scan. The PDF reproducibility at 300 s (mean, 0.86; range, 0.70-0.96) were significantly (p < 0.001) increased compared with those at 5 s (mean, 0.65; range, 0.25-0.79). PDF reproducibility showed less sensitivity to imaging frame rates that were >2 frames/s. CONCLUSION A statistically significant improvement in PDF reproducibility was observed with a prolonged scan time among the 17 participants. The confirmation of PDF reproducibility over times much shorter than stereotactic body radiotherapy delivery duration is a vital part of the initial validation process of probability-based treatment planning for stereotactic body radiotherapy for lung cancer.
Collapse
Affiliation(s)
- Jing Cai
- Department of Radiation Oncology, University of Virginia, Charlottesville, VA 22908, USA
| | | | | | | | | | | |
Collapse
|