Technical Note: On the spatial correlation between robust CT-ventilation methods and SPECT ventilation

Edward Castillo, William Beaumont Hospital
Richard Castillo, Emory University
Yevgeniy Vinogradskiy, University of Colorado Denver
Girish Nair, William Beaumont Hospital
Inga Grills, William Beaumont Hospital
Thomas Guerrero, William Beaumont Hospital
Craig Stevens, William Beaumont Hospital


© 2020 The Authors. Medical Physics published by Wiley Periodicals LLC on behalf of American Association of Physicists in Medicine. Purpose: The computed tomography (CT)-derived ventilation imaging methodology employs deformable image registration (DIR) to recover respiratory motion-induced volume changes from an inhale/exhale CT image pair, as a surrogate for ventilation. The Integrated Jacobian Formulation (IJF) and Mass Conserving Volume Change (MCVC) numerical methods for volume change estimation represent two classes of ventilation methods, namely transformation based and intensity (Hounsfield Unit) based, respectively. Both the IJF and MCVC methods utilize subregional volume change measurements that satisfy a specified uncertainty tolerance. In previous publications, the ventilation images resulting from this numerical strategy demonstrated robustness to DIR variations. However, the reduced measurement uncertainty comes at the expense of measurement resolution. The purpose of this study was to examine the spatial correlation between robust CT-ventilation images and single photon emission CT-ventilation (SPECT-V). Methods: Previously described implementations of IJF and MCVC require the solution of a large scale, constrained linear least squares problem defined by a series of robust subregional volume change measurements. We introduce a simpler parameterized implementation that reduces the number of unknowns while increasing the number of data points in the resulting least squares problem. A parameter sweep of the measurement uncertainty tolerance, (Formula presented.), was conducted using the 4DCT and SPECT-V images acquired for 15 non-small cell lung cancer patients prior to radiotherapy. For each test case, MCVC and IJF CT-ventilation images were created for 30 different uncertainty parameter values, uniformly sampled from the range (Formula presented.). Voxel-wise Spearman correlation between the SPECT-V and the resulting CT-ventilation images was computed. Results: The median correlations between MCVC and SPECT-V ranged from 0.20 to 0.48 across the parameter sweep, while the median correlations for IJF and SPECT-V ranged between 0.79 and 0.82. For the optimal IJF tolerance (Formula presented.), the IJF and SPECT-V correlations across all 15 test cases ranged between 0.12 and 0.90. For the optimal MCVC tolerance (Formula presented.), the MCVC and SPECT-V correlations across all 15 test cases ranged between −0.06 and 0.84. Conclusion: The reported correlations indicate that robust methods generate ventilation images that are spatially consistent with SPECT-V, with the transformation-based IJF method yielding higher correlations than those previously reported in the literature. For both methods, overall correlations were found to marginally vary for (Formula presented.), indicating that the clinical utility of both methods is robust to both uncertainty tolerance and DIR solution.