Estimating lung ventilation directly from 4D CT Hounsfield unit values

• Published in: Medical physics (Lancaster) 2016-01, Vol.43 (1), p.33-43
• Main Authors: Kipritidis, John, Hofman, Michael S., Siva, Shankar, Callahan, Jason, Le Roux, Pierre‐Yves, Woodruff, Henry C., Counter, William B., Keall, Paul J.
• Format: Article
• Language: English
• Subjects: 4DCT; Biological material, e.g. blood, urine; Haemocytometers; Biomedical modeling; cancer; Carcinoma, Non-Small-Cell Lung - diagnostic imaging; Carcinoma, Non-Small-Cell Lung - physiopathology; Computed tomography; Computerised tomographs; computerised tomography; CT ventilation imaging; deformable image registration; Digital computing or data processing equipment or methods, specially adapted for specific applications; Dosimetry; Four-Dimensional Computed Tomography; Humans; Image data processing or generation, in general; Image Processing, Computer-Assisted - methods; image registration; lung; Lung - diagnostic imaging; Lung - physiopathology; lung cancer; Lung Neoplasms - diagnostic imaging; Lung Neoplasms - physiopathology; Lungs; Measuring half‐life of a radioactive substance; Medical image artifacts; medical image processing; Medical image segmentation; PET ventilation; Pneumodyamics, respiration; pneumodynamics; Positron emission tomography (PET); Positron-Emission Tomography; Pulmonary Ventilation; radiation therapy; Scintigraphy; Time Factors; Tissues; tumours
• ISSN: 0094-2405; 2473-4209
• DOI: 10.1118/1.4937599

Purpose: Computed tomography ventilation imaging (CTVI) aims to visualize air‐volume changes in the lung by quantifying respiratory motion in 4DCT using deformable image registration (DIR). A problem is that DIR‐based CTVI is sensitive both to 4DCT image artifacts and DIR parameters, hindering clinical validation of the technique. To address this, the authors present a streamlined CTVI approach that estimates blood–gas exchange in terms of time‐averaged 4DCT Hounsfield unit (HU) values without relying on DIR. The purpose of this study is to quantify the accuracy of the HU‐based CTVI method using high‐resolution 68Ga positron emission tomography (“Galligas PET”) scans in lung cancer patients. Methods: The authors analyzed Galligas 4D‐PET/CT scans acquired for 25 lung cancer patients at up to three imaging timepoints during lung cancer radiation therapy. For each 4DCT scan, the authors produced three types of CTVIs: (i) the new method (CTV IHU¯), which takes the 4D time‐averaged product of regional air and tissue densities at each voxel, and compared this to DIR‐based estimates of (ii) breathing‐induced density changes (CTV IDIR‐HU), and (iii) breathing‐induced volume changes (CTV IDIR‐Jac) between the exhale/inhale phase images. The authors quantified the accuracy of CTV IHU¯, CTV IDIR‐HU and CTV IDIR‐Jac versus Galligas PET in terms of voxel‐wise Spearman correlation (r) and the separation of mean voxel values between clinically defined defect/nondefect regions. Results: Averaged over 62 scans, CTV IHU¯ showed better accuracy than CTV IDIR‐HU and CTV IDIR‐Jac in terms of Spearman correlation with Galligas PET, with (mean ± SD) r values of (0.50 ± 0.17), (0.42 ± 0.20), and (0.19 ± 0.23), respectively. A two‐sample Kolmogorov–Smirnov test indicates that CTV IHU¯ shows statistically significant separation of mean ventilation values between clinical defect/nondefect regions. Qualitatively, CTV IHU¯ appears concordant with Galligas PET for emphysema related defects, but differences arise in tumor‐obstructed regions (where aeration is overestimated due to motion blur) and for other abnormal morphology (e.g., fluid‐filled or peritumoral lung with HU ≳ − 600) where the assumptions of the HU model may break down. Conclusions: The HU‐based CTVI method can improve voxel‐wise correlations with Galligas PET compared to DIR‐based methods and may be a useful approximation for voxels with HU values in the range (−1000,   − 600). With further clinical verification, HU‐base