Direct measurement of lung motion using hyperpolarized helium-3 MR tagging

Role of lung lobar sliding on parenchymal distortion during breathing

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.

Respiratory Invariant Textures From Static Computed Tomography Scans for Explainable Lung Function Characterization

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.

Volumetric multiphase ventilation imaging based on four-dimensional computed tomography for functional lung avoidance radiotherapy

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.

A robust deformable image registration enhancement method based on radial basis function

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.

A hybrid proton and hyperpolarized gas tagging MRI technique for lung respiratory motion imaging: a feasibility study

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.

An initial investigation of hyperpolarized gas tagging magnetic resonance imaging in evaluating deformable image registration-based lung ventilation

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.

Lung surface deformation prediction from spirometry measurement and chest wall surface motion

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.

Biomedical imaging with hyperpolarized noble gases

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.

A measure for characterizing sliding on lung boundaries

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.

Displacement smoothing for the precise MRI-based measurement of strain in soft biological tissues

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.

Pulmonary kinematics from image data: a review

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.

Imaging of lung function using hyperpolarized helium-3 magnetic resonance imaging: Review of current and emerging translational methods and applications

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.

Helical tomotherapy planning for lung cancer based on ventilation magnetic resonance imaging

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.

Pulmonary kinematics from tagged hyperpolarized helium-3 MRI

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.

Dynamic MRI of grid-tagged hyperpolarized helium-3 for the assessment of lung motion during breathing

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.

Reproducibility of interfraction lung motion probability distribution function using dynamic MRI: statistical analysis

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.