A decay-modeled compressed sensing reconstruction approach for non-Cartesian hyperpolarized 129Xe MRI.

PURPOSE: Hyperpolarized 129Xe MRI benefits from non-Cartesian acquisitions that sample k-space efficiently and rapidly. However, their reconstructions are complex and burdened by decay processes unique to hyperpolarized gas. Currently used gridded reconstructions are prone to artifacts caused by magnetization decay and are ill-suited for undersampling. We present a compressed sensing (CS) reconstruction approach that incorporates magnetization decay in the forward model, thereby producing images with increased sharpness and contrast, even in undersampled data. METHODS: Radio-frequency, T1, and T 2 * $$ {\mathrm{T}}_2^{\ast } $$ decay processes were incorporated into the forward model and solved using iterative methods including CS. The decay-modeled reconstruction was validated in simulations and then tested in 2D/3D-spiral ventilation and 3D-radial gas-exchange MRI. Quantitative metrics including apparent-SNR and sharpness were compared between gridded, CS, and twofold undersampled CS reconstructions. Observations were validated in gas-exchange data collected from 15 healthy and 25 post-hematopoietic-stem-cell-transplant participants. RESULTS: CS reconstructions in simulations yielded images with threefold increases in accuracy. CS increased sharpness and contrast for ventilation in vivo imaging and showed greater accuracy for undersampled acquisitions. CS improved gas-exchange imaging, particularly in the dissolved-phase where apparent-SNR improved, and structure was made discernable. Finally, CS showed repeatability in important global gas-exchange metrics including median dissolved-gas signal ratio and median angle between real/imaginary components. CONCLUSION: A non-Cartesian CS reconstruction approach that incorporates hyperpolarized 129Xe decay processes is presented. This approach enables improved image sharpness, contrast, and overall image quality in addition to up-to threefold undersampling. This contribution benefits all hyperpolarized gas MRI through improved accuracy and decreased scan durations.

Finite element simulations of hyperpolarized gas DWI in micro-CT meshes of acinar airways: validating the cylinder and stretched exponential models of lung microstructural length scales.

PURPOSE: This work assesses the accuracy of the stretched exponential (SEM) and cylinder models of lung microstructural length scales that can be derived from hyperpolarized gas DWI. This was achieved by simulating 3 He and 129 Xe DWI signals within two micro-CT-derived realistic acinar airspace meshes that represent healthy and idiopathic pulmonary fibrosis lungs. METHODS: The healthy and idiopathic pulmonary fibrosis acinar airway meshes were derived from segmentations of 3D micro-CT images of excised human lungs and meshed for finite element simulations of the Bloch-Torrey equations. 3 He and 129 Xe multiple b value DWI experiments across a range of diffusion times (3 He Δ = 1.6 ms; 129 Xe Δ = 5 to 20 ms) were simulated in each mesh. Global SEM mean diffusive length scale and cylinder model mean chord length value was derived from each finite element simulation and compared against each mesh's mean linear intercept length, calculated from intercept length measurements within micro-CT segmentation masks. RESULTS: The SEM-derived mean diffusive length scale was within ±10% of the mean linear intercept length for simulations with both 3 He (Δ = 1.6 ms) and 129 Xe (Δ = 7 to 13 ms) in the healthy mesh, and with 129 Xe (Δ = 13 to 20 ms) for the idiopathic pulmonary fibrosis mesh, whereas for the cylinder model-derived mean chord length the closest agreement with mean linear intercept length (11.7% and 22.6% difference) was at 129 Xe Δ = 20 ms for both healthy and IPF meshes, respectively. CONCLUSION: This work validates the use of the SEM for accurate estimation of acinar dimensions and indicates that the SEM is relatively robust across a range of experimental conditions and acinar length scales.

Incorporation of Prior Knowledge of Signal Behavior Into the Reconstruction to Accelerate the Acquisition of Diffusion MRI Data.

Diffusion MRI data are generally acquired using hyperpolarized gases during patient breath-hold, which yields a compromise between achievable image resolution, lung coverage, and number of -values. In this paper, we propose a novel method that accelerates the acquisition of diffusion MRI data by undersampling in both the spatial and -value dimensions and incorporating knowledge about signal decay into the reconstruction (SIDER). SIDER is compared with total variation (TV) reconstruction by assessing its effect on both the recovery of ventilation images and the estimated mean alveolar dimensions (MADs). Both methods are assessed by retrospectively undersampling diffusion data sets ( =8) of healthy volunteers and patients with Chronic Obstructive Pulmonary Disease (COPD) for acceleration factors between x2 and x10. TV led to large errors and artifacts for acceleration factors equal to or larger than x5. SIDER improved TV, with a lower solution error and MAD histograms closer to those obtained from fully sampled data for acceleration factors up to x10. SIDER preserved image quality at all acceleration factors, although images were slightly smoothed and some details were lost at x10. In conclusion, we developed and validated a novel compressed sensing method for lung MRI imaging and achieved high acceleration factors, which can be used to increase the amount of data acquired during breath-hold. This methodology is expected to improve the accuracy of estimated lung microstructure dimensions and provide more options in the study of lung diseases with MRI.

Regional Ventilation Changes in the Lung: Treatment Response Mapping by Using Hyperpolarized Gas MR Imaging as a Quantitative Biomarker.

Purpose To assess the magnitude of regional response to respiratory therapeutic agents in the lungs by using treatment response mapping (TRM) with hyperpolarized gas magnetic resonance (MR) imaging. TRM was used to quantify regional physiologic response in adults with asthma who underwent a bronchodilator challenge. Materials and Methods This study was approved by the national research ethics committee and was performed with informed consent. Imaging was performed in 20 adult patients with asthma by using hyperpolarized helium 3 (3He) ventilation MR imaging. Two sets of baseline images were acquired before inhalation of a bronchodilating agent (salbutamol 400 µg), and one set was acquired after. All images were registered for voxelwise comparison. Regional treatment response, ΔR(r), was calculated as the difference in regional gas distribution (R[r] = ratio of inhaled gas to total volume of a voxel when normalized for lung inflation volume) before and after intervention. A voxelwise activation threshold from the variability of the baseline images was applied to ΔR(r) maps. The summed global treatment response map (ΔRnet) was then used as a global lung index for comparison with metrics of bronchodilator response measured by using spirometry and the global imaging metric percentage ventilated volume (%VV). Results ΔRnet showed significant correlation (P < .01) with changes in forced expiratory volume in 1 second (r = 0.70), forced vital capacity (r = 0.84), and %VV (r = 0.56). A significant (P < .01) positive treatment effect was detected with all metrics; however, ΔRnet showed a lower intersubject coefficient of variation (64%) than all of the other tests (coefficient of variation, ≥99%). Conclusion TRM provides regional quantitative information on changes in inhaled gas ventilation in response to therapy. This method could be used as a sensitive regional outcome metric for novel respiratory interventions. © RSNA, 2017 Online supplemental material is available for this article.

Detection of early subclinical lung disease in children with cystic fibrosis by lung ventilation imaging with hyperpolarised gas MRI.

Hyperpolarised 3He ventilation-MRI, anatomical lung MRI, lung clearance index (LCI), low-dose CT and spirometry were performed on 19 children (6-16 years) with clinically stable mild cystic fibrosis (CF) (FEV1>-1.96), and 10 controls. All controls had normal spirometry, MRI and LCI. Ventilation-MRI was the most sensitive method of detecting abnormalities, present in 89% of patients with CF, compared with CT abnormalities in 68%, LCI 47% and conventional MRI 22%. Ventilation defects were present in the absence of CT abnormalities and in patients with normal physiology, including LCI. Ventilation-MRI is thus feasible in young children, highly sensitive and provides additional information about lung structure-function relationships.

Whole lung morphometry with 3D multiple b-value hyperpolarized gas MRI and compressed sensing.

PURPOSE: To demonstrate three-dimensional (3D) multiple b-value diffusion-weighted (DW) MRI of hyperpolarized 3 He gas for whole lung morphometry with compressed sensing (CS). METHODS: A fully-sampled, two b-value, 3D hyperpolarized 3 He DW-MRI dataset was acquired from the lungs of a healthy volunteer and retrospectively undersampled in the ky and kz phase-encoding directions for CS simulations. Optimal k-space undersampling patterns were determined by minimizing the mean absolute error between reconstructed and fully-sampled 3 He apparent diffusion coefficient (ADC) maps. Prospective three-fold, undersampled, 3D multiple b-value 3 He DW-MRI datasets were acquired from five healthy volunteers and one chronic obstructive pulmonary disease (COPD) patient, and the mean values of maps of ADC and mean alveolar dimension (LmD ) were validated against two-dimensional (2D) and 3D fully-sampled 3 He DW-MRI experiments. RESULTS: Reconstructed undersampled datasets showed no visual artifacts and good preservation of the main image features and quantitative information. A good agreement between fully-sampled and prospective undersampled datasets was found, with a mean difference of +3.4% and +5.1% observed in mean global ADC and LmD values, respectively. These differences were within the standard deviation range and consistent with values reported from healthy and COPD lungs. CONCLUSIONS: Accelerated CS acquisition has facilitated 3D multiple b-value 3 He DW-MRI scans in a single breath-hold, enabling whole lung morphometry mapping. Magn Reson Med 77:1916-1925, 2017. © 2016 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Finite element modeling of (129)Xe diffusive gas exchange NMR in the human alveoli.

Existing models of (129)Xe diffusive exchange for lung microstructural modeling with time-resolved MR spectroscopy data have considered analytical solutions to one-dimensional, homogeneous models of the lungs with specific assumptions about the alveolar geometry. In order to establish a model system for simulating the effects of physiologically-realistic changes in physical and microstructural parameters on (129)Xe exchange NMR, we have developed a 3D alveolar capillary model for finite element analysis. To account for the heterogeneity of the alveolar geometry across the lungs, we have derived realistic geometries for finite element analysis based on 2D histological samples and 3D micro-CT image volumes obtained from ex vivo biopsies of lung tissue from normal subjects and patients with interstitial lung disease. The 3D alveolar capillary model permits investigation of the impact of alveolar geometrical parameters and diffusion and perfusion coefficients on the in vivo measured (129)Xe CSSR signal response. The heterogeneity of alveolar microstructure that is accounted for in image-based models resulted in considerable alterations to the shape of the (129)Xe diffusive uptake curve when compared to 1D models. Our findings have important implications for the future design and optimization of (129)Xe MR experiments and in the interpretation of lung microstructural changes from this data.

Comparison of CT-based Lobar Ventilation with 3He MR Imaging Ventilation Measurements.

PURPOSE: To compare lobar lung ventilation computed from expiratory and inspiratory computed tomographic (CT) data with direct measurements of ventilation at hyperpolarized helium 3 ((3)He) magnetic resonance (MR) imaging by using same-breath hydrogen 1 ((1)H) MR imaging examinations to coregister the multimodality images. MATERIALS AND METHODS: The study was approved by the national research ethics committee, and written patient consent was obtained. Thirty patients with asthma underwent breath-hold CT at total lung capacity and functional residual capacity. (3)He and (1)H MR images were acquired during the same breath hold at a lung volume of functional residual capacity plus 1 L. Lobar segmentations delineated by major fissures on both CT scans were used to calculate the percentage of ventilation per lobe from the change in inspiratory and expiratory lobar volumes. CT-based ventilation was compared with (3)He MR imaging ventilation by using diffeomorphic image registration of (1)H MR imaging to CT, which enabled indirect registration of (3)He MR imaging to CT. Statistical analysis was performed by using the Wilcoxon signed-rank test, Pearson correlation coefficient, and Bland-Altman analysis. RESULTS: The mean ± standard deviation absolute difference between the CT and (3)He MR imaging percentage of ventilation volume in all lobes was 4.0% (right upper and right middle lobes, 5.4% ± 3.3; right lower lobe, 3.7% ± 3.9; left upper lobe, 2.8% ± 2.7; left lower lobe, 3.9% ± 2.6; Wilcoxon signed-rank test, P < .05). The Pearson correlation coefficient between the two techniques in all lobes was 0.65 (P < .001). Greater percentage of ventilation was seen in the upper lobes with (3)He MR imaging and in the lower lobes with CT. This was confirmed with Bland-Altman analysis, with 95% limits of agreement for right upper and middle lobes, -2.4, 12.7; right lower lobe, -11.7, 4.6; left upper lobe, -4.9, 8.7; and left lower lobe, -9.8, 2.8. CONCLUSION: The percentage of regional ventilation per lobe calculated at CT was comparable to a direct measurement of lung ventilation at hyperpolarized (3)He MR imaging. This work provides evidence for the validity of the CT model, and same-breath (1)H MR imaging enables regional interpretation of (3)He ventilation MR imaging on the underlying lung anatomy at thin-section CT.

Pulmonary ventilation and micro-structural findings in congenital diaphragmatic hernia.

BACKGROUND: With increasing survival of patients with more severe forms of congenital diaphragmatic hernia (CDH) and risk of long-term respiratory morbidity, studies on lung morphology are needed. We used hyperpolarised (3) He MRI and anatomical (1) H MRI in a cohort of young adult CDH patients to image regional lung ventilation and microstructure, focusing on morphological and micro-structural (alveolar) abnormalities. METHODS: Nine patients with left-sided CDH, born 1975-1993, were studied. Regional ventilation was imaged with hyperpolarised (3) He MRI, and the (3) He apparent diffusion coefficient (ADC) was computed separately for the ipsilateral and contralateral lungs. (1) H MRI was used to image lung anatomy, total lung volume and motion during free-breathing. RESULTS: (3) He MRI showed ventilation abnormalities in six patients, ranging from a single ipsilateral ventilation defect (3 patients) to multiple ventilation defects in both lungs (one patient treated with extra corporeal membrane oxygenation). In eight patients, (3) He ADC values for the ipsilateral lung were significantly higher than those for the contralateral lung. CONCLUSIONS: Functional and micro-structural changes persist into adulthood in most CDH patients. Ipsilateral elevated (3) He ADC values are consistent with enlargement of mean dimensions of the confining lung micro-structure at the alveolar level.

Experimental validation of the hyperpolarized 129 Xe chemical shift saturation recovery technique in healthy volunteers and subjects with interstitial lung disease.

PURPOSE: To assess the sensitivity of the hyperpolarized 129 Xe chemical shift saturation recovery (CSSR) technique for noninvasive quantification of changes to lung microstructure and function in idiopathic pulmonary fibrosis (IPF) and systemic sclerosis (SSc). METHODS: Ten healthy volunteers, four subjects with SSc and four with IPF were scanned at 1.5 T. A CSSR pulse sequence was implemented using binomial-composite radiofrequency pulses to monitor 129 Xe magnetization in tissues and blood plasma (T/P) and red blood cells (RBCs). The dynamics of 129 Xe uptake into these compartments were fitted with three existing analytical models of gas diffusion to extract parameters of lung physiology. These parameters were quantitatively compared between models. RESULTS: Uptake of xenon into the pulmonary capillaries was impaired in subjects with IPF and SSc. Statistically significant septal thickening was measured by 129 Xe CSSR in IPF patients. Preliminary data suggests age-dependent alterations to septal thickness in healthy volunteers. These findings were reproduced using each of the literature models. CSSR-derived parameters were compared with gold-standard indicators of pulmonary function; diffusing capacity of carbon monoxide and pulmonary transit-time. CONCLUSIONS: CSSR with hyperpolarized 129 Xe is sensitive to pathology-induced degradation of lung structure/function and shows promise for quantification of disease severity and monitoring treatment response. Magn Reson Med 74:196-207, 2015. © 2014 Wiley Periodicals, Inc.

A method for quantitative analysis of regional lung ventilation using deformable image registration of CT and hybrid hyperpolarized gas/1H MRI.

Hyperpolarized gas magnetic resonance imaging (MRI) generates highly detailed maps of lung ventilation and physiological function while CT provides corresponding anatomical and structural information. Fusion of such complementary images enables quantitative analysis of pulmonary structure-function. However, direct image registration of hyperpolarized gas MRI to CT is problematic, particularly in lungs whose boundaries are difficult to delineate due to ventilation heterogeneity. This study presents a novel indirect method of registering hyperpolarized gas MRI to CT utilizing (1)H-structural MR images that are acquired in the same breath-hold as the gas MRI. The feasibility of using this technique for regional quantification of ventilation of specific pulmonary structures is demonstrated for the lobes.The direct and indirect methods of hyperpolarized gas MRI to CT image registration were compared using lung images from 15 asthma patients. Both affine and diffeomorphic image transformations were implemented. Registration accuracy was evaluated using the target registration error (TRE) of anatomical landmarks identified on (1)H MRI and CT. The Wilcoxon signed-rank test was used to test statistical significance.For the affine transformation, the indirect method of image registration was significantly more accurate than the direct method (TRE = 14.7 ± 3.2 versus 19.6 ± 12.7 mm, p = 0.036). Using a deformable transformation, the indirect method was also more accurate than the direct method (TRE = 13.5 ± 3.3 versus 20.4 ± 12.8 mm, p = 0.006).Accurate image registration is critical for quantification of regional lung ventilation with hyperpolarized gas MRI within the anatomy delineated by CT. Automatic deformable image registration of hyperpolarized gas MRI to CT via same breath-hold (1)H MRI is more accurate than direct registration. Potential applications include improved multi-modality image fusion, functionally weighted radiotherapy planning, and quantification of lobar ventilation in obstructive airways disease.

(3)He pO2 mapping is limited by delayed-ventilation and diffusion in chronic obstructive pulmonary disease.

PURPOSE: Lung pO2 mapping with (3)He MRI assumes that the sources of signal decay with time during a breath-hold are radiofrequency depolarization and oxygen-dependent T1 relaxation, but the method is sensitive to other sources of spatio-temporal signal change such as diffusion. The purpose of this work was to assess the use of (3)He pO2 mapping in patients with chronic obstructive pulmonary disease. METHODS: Ten patients with moderate to severe chronic obstructive pulmonary disease were scanned with a 3D single breath-hold pO2 mapping sequence. RESULTS: Images showed signal increasing over time in some lung regions due to delayed ventilation during breath-hold. Regions of physically unrealistic negative pO2 values were seen in all patients, and regional mean pO2 values of -0.3 bar were measured in the two patients most affected by delayed ventilation (where mean time to signal onset was 3-4 s). CONCLUSIONS: Movement of gas within the lungs during breath-hold causes regional changes in signal over time that are not related to oxygen concentration, leading to erroneous pO2 measurements using the linear oxygen-dependent signal decay model. These spatio-temporal sources of signal change cannot be reliably separated at present, making pO2 mapping using this methodology unreliable in chronic obstructive pulmonary disease patients with significant bullous emphysema or delayed ventilation.

Quantification of regional fractional ventilation in human subjects by measurement of hyperpolarized 3He washout with 2D and 3D MRI.

Multiple-breath washout hyperpolarized (3)He MRI was used to calculate regional parametric images of fractional ventilation (r) as the ratio of fresh gas entering a volume unit to the total end inspiratory volume of the unit. Using a single dose of inhaled hyperpolarized gas and a total acquisition time of under 1 min, gas washout was measured by dynamic acquisitions during successive breaths with a fixed delay. A two-dimensional (2D) imaging protocol was investigated in four healthy subjects in the supine position, and in a second protocol the capability of extending the washout imaging to a three-dimensional (3D) acquisition covering the whole lungs was tested. During both protocols, subjects were breathing comfortably, only restricted by synchronization of breathing to the sequence timings. The 3D protocol was also successfully tested on one patient with cystic fibrosis. Mean r values from each volunteer were compared with global gas volume turnover, as calculated from flow measurement at the mouth divided by total lung volume (from MRI images), and a significant correlation (r = 0.74, P < 0.05) was found. The effects of gravity on R were investigated, and an average decrease in r of 5.5%/cm (Δr = 0.016 ± 0.006 cm(-1)) from posterior to anterior was found in the right lung. Intersubject reproducibility of r imaging with the 2D and 3D protocol was tested, and a significant correlation between repeated experiments was found in a pixel-by-pixel comparison. The proposed methods can be used to measure r on a regional basis.

Response to Commentary on "The influence of lung airways branching structure and diffusion time on measurements and models of short-range 3He gas MR diffusion".

Our extensive investigation of the cylinder model theory through numerical modelling and purpose-designed experiments has demonstrated that it does produce inaccurate estimates of airway dimensions at all diffusion times currently used. This is due to a variety of effects: incomplete treatment of non-Gaussian effects, finite airway size, branching geometry, background susceptibility gradients and diffusion time dependence of the (3)He MR diffusion behaviour in acinar airways. The cylinder model is a good starting point for the development of a lung morphometry technique from (3)He diffusion MR but its limitations need to be understood and documented in the interest of reliable clinical interpretation.

Simultaneous imaging of lung structure and function with triple-nuclear hybrid MR imaging.

PURPOSE: To re-engineer a standard clinical magnetic resonance (MR) imaging system to enable the acquisition, in the same breath hold, of lung images from two hyperpolarized gases (helium 3 [(3)He] and xenon 129 [(129)Xe]) with simultaneous registered anatomic proton (hydrogen 1 [(1)H]) MR images of lung structure. MATERIALS AND METHODS: Studies with (3)He and (129)Xe were performed with National Research Ethics Committee approval, with informed consent from the volunteer. (1)H-(3)He-(129)Xe MR imaging was achieved in the same breath by using mutually decoupled nested radiofrequency coil hardware capable of transmit and receive on each respective nucleus without power cross talk. MR pulse sequences were also developed for rapid switching between each nucleus. The system is demonstrated with triple-nuclear lung images in a healthy individual following inhalation of a mixture of (3)He and (129)Xe gases. RESULTS: Spatially and temporally registered images of all three nuclei were obtained with high signal to noise ratio and high spatial resolution in the same breath. CONCLUSION: The multinuclear technique is capable of providing registered lung images with mutually complementary functional and structural spatial information.

The influence of lung airways branching structure and diffusion time on measurements and models of short-range 3He gas MR diffusion.

Hyperpolarized (3)He diffusion experiments have been shown to be sensitive to changes in acinar structure due to emphysematous lung disease. Extracting quantitative information about lung microstructure from the diffusion signal is complicated due its dependence on a number of factors including diffusion time and the complex branching acinar geometry. A theoretical model (cylinder model) has been proposed as a means of estimating acinar airway dimensions from measured diffusivities. This model assumes that the effects of acinar branching geometry and finite airway length upon (3)He diffusion behaviour are negligible. In this work, we use finite element simulations of diffusion in a model of branching alveolar ducts to investigate in detail the effects of acinar branching structure and finite airway length on short-range (3)He diffusion measurements. The results show that branching effects have a significant influence upon (3)He diffusivity, even at short diffusion times. The expressions of the cylinder model theory do not account for significant dependences upon diffusion time, branching geometry and airway length, as a consequence of the oversimplified geometrical model used. The effect of diffusion time on (3)He ADC was also investigated through experiments with healthy human volunteers. The results demonstrate that the cylinder model can produce inaccurate estimates of the airway dimensions as a consequence of incompletely accounting for the diffusion-time dependence in the model equations and confirmed the predicted limitations of the cylinder model for reliable lung morphometry measurements. The results and models presented in this work may help in the development of a more realistic theoretical framework for 'in vivo lung morphometry' using (3)He diffusion MR.

Hyperpolarized 129Xe gas lung MRI-SNR and T2* comparisons at 1.5 T and 3 T.

In this study, the signal-to-noise ratio of hyperpolarized (129)Xe human lung magnetic resonance imaging was compared at 1.5 T and 3 T. Experiments were performed at both B(0) fields with quadrature double Helmholtz transmit-receive chest coils of the same geometry with the same subject loads. Differences in sensitivity between the two field strengths were assessed from the signal-to-noise ratio of multi-slice 2D (129)Xe ventilation lung images obtained at the two field strengths with a spatial resolution of 15 mm × 4 mm × 4 mm. There was a systematically higher signal-to-noise ratio observed at 3 T than at 1.5 T by a factor of 1.25. Mean image signal-to-noise ratio was in the range 27-44 at 1.5 T and 36-51 at 3 T. T 2* of (129)Xe gas in the partially inflated lungs was measured to be 25 ms and 18 ms at 1.5 T and 3 T, respectively. T 2* of (129)Xe gas in fully inflated lungs was measured to be 52 ms and 24 ms at 1.5 T and 3 T, respectively.

Direct visualisation of collateral ventilation in COPD with hyperpolarised gas MRI.

BACKGROUND: Collateral ventilation has been proposed as a mechanism of compensation of respiratory function in obstructive lung diseases but observations of it in vivo are limited. The assessment of collateral ventilation with an imaging technique might help to gain insight into lung physiology and assist the planning of new bronchoscopic techniques for treating emphysema. OBJECTIVE: To obtain images of delayed ventilation that might be related to collateral ventilation over the period of a single breath-hold in patients with chronic obstructive pulmonary disease (COPD). METHODS: Time-resolved breath-hold hyperpolarised (3)He MRI was used to obtain images of the progressive influx of polarised gas into initially non-ventilated defects. RESULTS: A time-series of images showed that (3)He moves into lung regions which were initially non-ventilated. Ventilation defects with delayed filling were observed in 8 of the 10 patients scanned. CONCLUSIONS: A method for direct imaging of delayed ventilation within a single breath-hold has been demonstrated in patients with COPD. Images of what is believed to be collateral ventilation and slow filling of peripheral airspaces due to increased flow resistance are presented. The technique provides 3D whole-lung coverage with sensitivity to regional information, and is non-invasive and non-ionising.