Browsing by Author "Myziuk, Nicholas"

Now showing 1 - 2 of 2
  • Thumbnail Image
    Item
    Quantifying pulmonary perfusion from noncontrast computed tomography
    (Wiley, 2021) Castillo, Edward; Nair, Girish; Turner‐Lawrence, Danielle; Myziuk, Nicholas; Emerson, Scott; Al‐Katib, Sayf; Westergaard, Sarah; Castillo, Richard; Vinogradskiy, Yevgeniy; Quinn, Thomas; Guerrero, Thomas; Stevens, Craig
    Purpose: Computed tomography (CT)-derived ventilation methods compute respiratory induced volume changes as a surrogate for pulmonary ventilation. Currently, there are no known methods to derive perfusion information from noncontrast CT. We introduce a novel CT-Perfusion (CT-P) method for computing the magnitude mass changes apparent on dynamic noncontrast CT as a surrogate for pulmonary perfusion. Methods: CT-Perfusion is based on a mass conservation model which describes the unknown mass change as a linear combination of spatially corresponding inhale and exhale HU estimated voxel densities. CT-P requires a deformable image registration (DIR) between the inhale/exhale lung CT pair, a preprocessing lung volume segmentation, and an estimate for the Jacobian of the DIR transformation. Given this information, the CT-P image, which provides the magnitude mass change for each voxel within the lung volume, is formulated as the solution to a constrained linear least squares problem defined by a series of subregional mean magnitude mass change measurements. Similar to previous robust CT-ventilation methods, the amount of uncertainty in a subregional sample mean measurement is related to measurement resolution and can be characterized with respect to a tolerance parameter. Spatial Spearman correlation between single photon emission CT perfusion (SPECT-P) and the proposed CT-P method was assessed in two patient cohorts via a parameter sweep of . The first cohort was comprised of 15 patients diagnosed with pulmonary embolism (PE) who had SPECT-P and 4DCT imaging acquired within 24 h of PE diagnosis. The second cohort was comprised of 15 nonsmall cell lung cancer patients who had SPECT-P and 4DCT images acquired prior to radiotherapy. For each test case, CT-P images were computed for 30 different uncertainty parameter values, uniformly sampled from the range [0.01, 0.125], and the Spearman correlation between the SPECT-P and the resulting CT-P images were computed. Results: The median correlations between CT-P and SPECT-P taken over all 30 test cases ranged between 0.49 and 0.57 across the parameter sweep. For the optimal tolerance τ = 0.0385, the CT-P and SPECT-P correlations across all 30 test cases ranged between 0.02 and 0.82. A one-sample sign test was applied separately to the PE and lung cancer cohorts. A low Spearmen correlation of 15% was set as the null median value and two-sided alternative was tested. The PE patients showed a median correlation of 0.57 (IQR = 0.305). One-sample sign test was statistically significant with 96.5 % confidence interval: 0.20–0.63, P < 0.00001. Lung cancer patients had a median correlation of 0.57(IQR = 0.230). Again, a one-sample sign test for median was statistically significant with 96.5 percent confidence interval: 0.45–0.71, P < 0.00001. Conclusion: CT-Perfusion is the first mechanistic model designed to quantify magnitude blood mass changes on noncontrast dynamic CT as a surrogate for pulmonary perfusion. While the reported correlations with SPECT-P are promising, further investigation is required to determine the optimal CT acquisition protocol and numerical method implementation for CT-P imaging.
  • Thumbnail Image
    Item
    Robust CT ventilation from the integral formulation of the Jacobian
    (Wiley, 2019) Castillo, Edward; Castillo, Richard; Vinogradskiy, Yevgeniy; Dougherty, Michele; Solis, David; Myziuk, Nicholas; Thompson, Andrew; Guerra, Rudy; Nair, Girish; Guerrero, Thomas
    Computed tomography (CT) derived ventilation algorithms estimate the apparent voxel volume changes within an inhale/exhale CT image pair. Transformation-based methods compute these estimates solely from the spatial transformation acquired by applying a deformable image registration (DIR) algorithm to the image pair. However, approaches based on finite difference approximations of the transformation's Jacobian have been shown to be numerically unstable. As a result, transformation-based CT ventilation is poorly reproducible with respect to both DIR algorithm and CT acquisition method. PURPOSE: We introduce a novel Integrated Jacobian Formulation (IJF) method for estimating voxel volume changes under a DIR-recovered spatial transformation. The method is based on computing volume estimates of DIR mapped subregions using the hit-or-miss sampling algorithm for integral approximation. The novel approach allows for regional volume change estimates that (a) respect the resolution of the digital grid and (b) are based on approximations with quantitatively characterized and controllable levels of uncertainty. As such, the IJF method is designed to be robust to variations in DIR solutions and thus overall more reproducible. METHODS: Numerically, Jacobian estimates are recovered by solving a simple constrained linear least squares problem that guarantees the recovered global volume change is equal to the global volume change obtained from the inhale and exhale lung segmentation masks. Reproducibility of the IJF method with respect to DIR solution was assessed using the expert-determined landmark point pairs and inhale/exhale phases from 10 four-dimensional computed tomographies (4DCTs) available on www.dir-lab.com. Reproducibility with respect to CT acquisition was assessed on the 4DCT and 4D cone beam CT (4DCBCT) images acquired for five lung cancer patients prior to radiotherapy. RESULTS: The ten Dir-Lab 4DCT cases were registered twice with the same DIR algorithm, but with different smoothing parameter. Finite difference Jacobian (FDJ) and IFJ images were computed for both solutions. The average spatial errors (300 landmarks per case) for the two DIR solution methods were 0.98 (1.10) and 1.02 (1.11). The average Pearson correlation between the FDJ images computed from the two DIR solutions was 0.83 (0.03), while for the IJF images it was 1.00 (0.00). For intermodality assessment, the IJF and FDJ images were computed from the 4DCT and 4DCBCT of five patients. The average Pearson correlation of the spatially aligned FDJ images was 0.27 (0.11), while it was 0.77 (0.13) for the IFJ method. CONCLUSION: The mathematical theory underpinning the IJF method allows for the generation of ventilation images that are (a) computed with respect to DIR spatial accuracy on the digital voxel grid and (b) based on DIR-measured subregional volume change estimates acquired with quantifiable and controllable levels of uncertainty. Analyses of the experiments are consistent with the mathematical theory and indicate that IJF ventilation imaging has a higher reproducibility with respect to both DIR algorithm and CT acquisition method, in comparison to the standard finite difference approach.