Dataset Information

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

ABSTRACT:

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, τ , 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 0.01,0.25 . 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 τ=0.07 , the IJF and SPECT-V correlations across all 15 test cases ranged between 0.12 and 0.90. For the optimal MCVC tolerance τ=0.03 , 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 τ∈[0.03,0.15] , indicating that the clinical utility of both methods is robust to both uncertainty tolerance and DIR solution.

SUBMITTER: Castillo E

PROVIDER: S-EPMC7727923 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description: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.PurposeWe 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.MethodsNumerically, 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.ResultsThe 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.ConclusionThe 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.

Project description:A novel, non-terminal surgical procedure to remove a single placentome from the pregnant ewe for gene expression and histological analyses was recently developed in our laboratory. This technique allows for evaluation of nutritional insults on placental development at more than one stage of gestation using a single animal. Early attempts to develop a similar technique in cattle were met with complications due to inaccessibility of the gravid uterine horn because of its location and mass. One alternative is to collect a placentome from the contralateral uterine horn; however, the question remains as to whether gene expression varies among placentomes based on location relative to the fetus. Pregnant heifers were maintained on forage during early gestation and later moved into pens with a Calan gate system (American Calan, Northwood, NH). On gestational day (GD) 158, five heifers were assigned to receive a hay-based diet formulated to meet 100% of NRC requirements, and five heifers were fed 70% of NRC requirements until necropsy on GD244. At necropsy, a single representative placentome was selected for analysis from the antimesometrial side: (1) of the gravid uterine horn central to the amnion, (2) over the allantois immediately adjacent to the amnion, (3) in the tip of the gravid uterine horn, and (4) in the tip of the contralateral uterine horn. Mean placentome weight was greater (P < 0.05) for locations central to the amnion and allantois compared to locations within the tips of the ipsilateral and contralateral horns, respectively. Gene expression for angiogenic factors (FGF2, ODC1, VEGFA, and FLT1), nutrient transporters (SLC7A1 and SLC2A1), and factors associated with hormone action (ESR1, IGF1, IGFBP3, CSH1, and PAG1) were unaffected (P > 0.05) by dietary treatment or location of the placentome. Results indicate that location of the placentome in relation to the fetus does not impact gene expression, enhancing the efficacy of nonterminal methodologies for sampling gene expression in placentomes.

Project description:BackgroundQuantitative SPECT imaging in targeted radionuclide therapy with lutetium-177 holds great potential for individualized treatment based on dose assessment. The establishment of dose-effect relations requires a standardized method for SPECT quantification. The purpose of this multi-center study is to evaluate quantitative accuracy and inter-system variations of different SPECT/CT systems with corresponding commercially available quantitative reconstruction algorithms. This is an important step towards a vendor-independent standard for quantitative lutetium-177 SPECT.MethodsFour state-of-the-art SPECT/CT systems were included: Discovery™ NM/CT 670Pro (GE Healthcare), Symbia Intevo™, and two Symbia™ T16 (Siemens Healthineers). Quantitative accuracy and inter-system variations were evaluated by repeatedly scanning a cylindrical phantom with 6 spherical inserts (0.5 - 113 ml). A sphere-to-background activity concentration ratio of 10:1 was used. Acquisition settings were standardized: medium energy collimator, body contour trajectory, photon energy window of 208 keV (± 10%), adjacent 20% lower scatter window, 2 × 64 projections, 128 × 128 matrix size, and 40 s projection time. Reconstructions were performed using GE Evolution with Q.Metrix™, Siemens xSPECT Quant™, Siemens Broad Quantification™ or Siemens Flash3D™ algorithms using vendor recommended settings. In addition, projection data were reconstructed using Hermes SUV SPECT™ with standardized reconstruction settings to obtain a vendor-neutral quantitative reconstruction for all systems. Volumes of interest (VOI) for the spheres were obtained by applying a 50% threshold of the sphere maximum voxel value corrected for background activity. For each sphere, the mean and maximum recovery coefficient (RCmean and RCmax) of three repeated measurements was calculated, defined as the imaged activity concentration divided by the actual activity concentration. Inter-system variations were defined as the range of RC over all systems.ResultsRC decreased with decreasing sphere volume. Inter-system variations with vendor-specific reconstructions were between 0.06 and 0.41 for RCmean depending on sphere size (maximum 118% quantification difference), and improved to 0.02-0.19 with vendor-neutral reconstructions (maximum 38% quantification difference).ConclusionThis study shows that eliminating sources of possible variation drastically reduces inter-system variation in quantification. This means that absolute SPECT quantification for 177Lu is feasible in a multi-center and multi-vendor setting; however, close agreement between vendors and sites is key for multi-center dosimetry and quantitative biomarker studies.

Project description:PURPOSE:To derive lobar ventilation in patients with chronic obstructive pulmonary disease (COPD) using a rapid time-series hyperpolarized xenon-129 (HPX) magnetic resonance imaging (MRI) technique and compare this to ventilation/perfusion single-photon emission computed tomography (V/Q-SPECT), correlating the results with high-resolution computed tomography (CT) and pulmonary function tests (PFTs). MATERIALS AND METHODS:Twelve COPD subjects (GOLD stages I-IV) participated in this study and underwent HPX-MRI, V/Q-SPECT/CT, high-resolution CT, and PFTs. HPX-MRI was performed using a novel time-series spiral k-space sampling approach. Relative percentage ventilations were calculated for individual lobe for comparison to the relative SPECT lobar ventilation and perfusion. The absolute HPX-MRI percentage ventilation in each lobe was compared to the absolute CT percentage emphysema score calculated using a signal threshold method. Pearson's correlation and linear regression tests were performed to compare each imaging modality. RESULTS:Strong correlations were found between the relative lobar percentage ventilation with HPX-MRI and percentage ventilation SPECT (r?=?0.644; p?<?0.001) and percentage perfusion SPECT (r?=?0.767; p?<?0.001). The absolute CT percentage emphysema and HPX percentage ventilation correlation was also statistically significant (r?=?0.695, p?<?0.001). The whole lung HPX percentage ventilation correlated with the PFT measurements (FEV1 with r?=?-?0.886, p?<?0.001*, and FEV1/FVC with r?=?-?0.861, p?<?0.001*) better than the whole lung CT percentage emphysema score (FEV1 with r?=?-?0.635, p?=?0.027; and FEV1/FVC with r?=?-?0.652, p?=?0.021). CONCLUSION:Lobar ventilation with HPX-MRI showed a strong correlation with lobar ventilation and perfusion measurements derived from SPECT/CT, and is better than the emphysema score obtained with high-resolution CT. KEY POINTS:• The ventilation hyperpolarized xenon-129 MRI correlates well with ventilation and perfusion with SPECT/CT with the advantage of higher temporal and spatial resolution. • The hyperpolarized xenon-129 MRI correlates with the PFT measurements better than the high-resolution CT with the advantage of avoiding the use of ionizing radiation.

Project description:PURPOSE:To apply a previously designed framework for predicting radiation pneumonitis by using pretreatment lung function heterogeneity metrics, anatomic dosimetry, and functional lung dosimetry derived from 2 imaging modalities within the same cohort. METHODS AND MATERIALS:Treatment planning computed tomography (CT) scans were co-registered with pretreatment [99mTc] macro-aggregated albumin perfusion single-photon positron emission tomography (SPECT)/CT scans and [18F]-fluorodeoxyglucose (FDG) positron emission tomography (PET)/CT scans of 28 patients who underwent definitive thoracic radiation. Clinical radiation pneumonitis was defined as grade ?2 (Common Terminology Criteria for Adverse Events, v. 4). Anatomic dosimetric parameters (mean lung dose [MLD], volume receiving ?20 Gy [V20]) were collected from treatment planning scans. Baseline functional lung heterogeneity parameters and functional lung dose-volume parameters were calculated from pretreatment SPECT/CT and FDG PET/CT scans. Functional heterogeneity parameters calculated over the tumor-subtracted lung included skewness, kurtosis, and coefficient of variation from perfusion SPECT and FDG PET and the global lung parenchymal glycolysis and mean standardized uptake value from FDG PET. Functional dose-volume parameters calculated in regions of highly functional lung, defined on perfusion (p) or SUV (s) images, included mean lung dose (pMLD, sMLD) and V20 (pV20, sV20). Fraction of integral lung function receiving ?20 Gy (pF20, sF20) was also calculated. Equivalent doses in 2 Gy per fraction (EQD2) were calculated to account for differences in treatment regimens and dose fractionation (EQD2Lung). RESULTS:Two anatomic dosimetric parameters (MLD, V20) and 4 functional dosimetric parameters (pMLD, pV20, pF20, sF20) were significant predictors of grade ?2 pneumonitis (area under the curve >0.84; P < .05). Dose-independent functional lung heterogeneity metrics were not associated with pneumonitis incidence. At thresholds of 100% sensitivity and 65% to 91% specificity, corresponding to maximum prediction accuracy for pneumonitis, these parameters had the following cutoff values: MLD = 13.6 Gy EQD2Lung, V20 = 25%, pMLD = 13.2 Gy EQD2Lung, pV20 = 15%, pF20 = 17%, and sF20 = 25%. Significant parameters MLD, V20, pF20, and sF20 were not cross-correlated to significant parameters pMLD and pV20, indicating that they may offer independently predictive information (Spearman ? < 0.7). CONCLUSIONS:We reported differences in anatomic and functional lung dosimetry between patients with and without pneumonitis in this limited patient cohort. Adding selected independent functional lung parameters may risk stratify patients for pneumonitis. Validation studies are ongoing in a prospective functional lung avoidance trial at our institution.

Project description:To evaluate and compare the theoretically achievable accuracy of two families of two-parameter photon cross-section models: basis vector model (BVM) and modified parametric fit model (mPFM).The modified PFM assumes that photoelectric absorption and scattering cross-sections can be accurately represented by power functions in effective atomic number and/or energy plus the Klein-Nishina cross-section, along with empirical corrections that enforce exact prediction of elemental cross-sections. Two mPFM variants were investigated: the widely used Torikoshi model (tPFM) and a more complex "VCU" variant (vPFM). For 43 standard soft and bony tissues and phantom materials, all consisting of elements with atomic number less than 20 (except iodine), we evaluated the theoretically achievable accuracy of tPFM and vPFM for predicting linear attenuation, photoelectric absorption, and energy-absorption coefficients, and we compared it to a previously investigated separable, linear two-parameter model, BVM.For an idealized dual-energy computed tomography (DECT) imaging scenario, the cross-section mapping process demonstrates that BVM more accurately predicts photon cross-sections of biological mixtures than either tPFM or vPFM. Maximum linear attenuation coefficient prediction errors were 15% and 5% for tPFM and BVM, respectively. The root-mean-square (RMS) prediction errors of total linear attenuation over the 20 keV to 1000 keV energy range of tPFM and BVM were 0.93% (tPFM) and 0.1% (BVM) for adipose tissue, 0.8% (tPFM) and 0.2% (BVM) for muscle tissue, and 1.6% (tPFM) and 0.2% (BVM) for cortical bone tissue. With exception of the thyroid and Teflon, the RMS error for photoelectric absorption and scattering coefficient was within 4% for the tPFM and 2% for the BVM. Neither model predicts the photon cross-sections of thyroid tissue accurately, exhibiting relative errors as large as 20%. For the energy-absorption coefficients prediction error, RMS errors for the BVM were less than 1.5%, while for the tPFM, the RMS errors were as large as 16%.Compared to modified PFMs, BVM shows superior potential to support dual-energy CT cross-section mapping. In addition, the linear, separable BVM can be more efficiently deployed by iterative model-based DECT image-reconstruction algorithms.

Project description:To determine the mean and range of location-averaged breast skin thickness using high-resolution dedicated breast CT for use in Monte Carlo-based estimation of normalized glandular dose coefficients.This study retrospectively analyzed image data from a clinical study investigating dedicated breast CT. An algorithm similar to that described by Huang et al. ["The effect of skin thickness determined using breast CT on mammographic dosimetry," Med. Phys. 35(4), 1199-1206 (2008)] was used to determine the skin thickness in 137 dedicated breast CT volumes from 136 women. The location-averaged mean breast skin thickness for each breast was estimated and the study population mean and range were determined. Pathology results were available for 132 women, and were used to investigate if the distribution of location-averaged mean breast skin thickness varied with pathology. The effect of surface fitting to account for breast curvature was also studied.The study mean (± interbreast SD) for breast skin thickness was 1.44 ± 0.25 mm (range: 0.87-2.34 mm), which was in excellent agreement with Huang et al. Based on pathology, pair-wise statistical analysis (Mann-Whitney test) indicated that at the 0.05 significance level, there were no significant difference in the location-averaged mean breast skin thickness distributions between the groups: benign vs malignant (p = 0.223), benign vs hyperplasia (p = 0.651), hyperplasia vs malignant (p = 0.229), and malignant vs nonmalignant (p = 0.172).Considering this study used a different clinical prototype system, and the study participants were from a different geographical location, the observed agreement between the two studies suggests that the choice of 1.45 mm thick skin layer comprising the epidermis and the dermis for breast dosimetry is appropriate. While some benign and malignant conditions could cause skin thickening, in this study cohort the location-averaged mean breast skin thickness distributions did not differ significantly with pathology. The study also underscored the importance of considering breast curvature in estimating breast skin thickness.