Search

Background: Prior implementations of the channelized Hotelling observer (CHO) model have succeeded in assessing the performance of X-ray angiography systems under a variety of imaging conditions. However, often times these conditions do not resemble those present in routine clinical imaging scenarios, such as having complex anthropomorphic backgrounds in conjunction with moving test objects., Purpose: This work builds up on prior established CHO methods and introduces a new approach to switch from the already established "multiple-sample" CHO implementation to a "single-sample" technique. The proposed implementation enables the inclusion of moving test objects upon nonuniform backgrounds by allowing only a single sample to represent the test object present condition that is to be used within the statistical test to estimate the detectability index., Methods: To assess the proposed method, two image data sets were acquired with a clinical X-ray angiography system. The first set consisted of a uniform background in combination with static test objects while the second consisted of an anthropomorphic chest phantom in conjunction with moving test objects. The first set was used to validate the proposed approach against the multiple-sample method while the second was used to assess the feasibility of the proposed method under a variety of imaging conditions, including seven object sizes and seven detector target dose (DTD) levels., Results: For the uniform background data set, considering all detectability indices greater or equal than 1, the ratio between the detectability indices of the proposed single-sample approach versus the multiple-sample method was 0.997 ± 0.056 (range 0.884-1.159). The average single-direction width of the 95% confidence intervals (CIs) of the detectability index estimates for the multiple-sample method was 0.38 ± 0.43 (range 0.03-2.20). For the single-sample approach, the average width was 2.52 ± 0.63 (range 1.11-5.44). For the anthropomorphic background image set, the results were consistent with classical quantum-limited signal-to-noise ratio (SNR) theory. The magnitude of the detectability indices varied predictably with changes in both object size and DTD, with the highest value associated with the highest dose and the largest object size. Additionally, the proposed method was able to capture differences in the imaging performance for a given test object across the field of view, which was associated with the attenuation levels provided by different features of the anthropomorphic background., Conclusions: A new single-sample variant of the CHO model to assess the performance of X-ray angiography imaging systems is proposed. The new approach is consistent with quantum-limited image quality theory and with a standard implementation of the CHO model. The proposed method enables the assessment of moving test objects in combination with complex, nonuniform image backgrounds, thereby opening the possibility to assess imaging conditions which more closely resemble those used in clinical care., (© 2023 American Association of Physicists in Medicine.)

Purpose: In previous works, it has been demonstrated that for filtered backprojection (FBP) reconstruction-based computed tomography (CT) images, the measured CT numbers are biased and the bias level decreases with increasing radiation dose. Low-dose scans typically include noise reduction schemes to reduce noise level. The purpose of this work was to investigate the potential impact of different noise reduction schemes on the CT number bias., Methods: Three different filtration methods: Gaussian, adaptive trimmed mean (ATM), and anisotropic diffusion (AD) were implemented to reduce noise. All filters were independently applied in three different domains: raw counts, log-processed sinogram, or reconstructed image domain. A quality assurance phantom was scanned on a benchtop CT cone beam CT system, at dose levels ranging from 0.6 to 4.0 mGy. The conventional FBP reconstructions were performed to reconstruct CT images for the study of CT number biases. The CT number bias of different material inserts in the phantom was then measured. To further study the overall impact of CT number bias together with the potential consequences of noise reduction schemes on both the spatial resolution and noise characteristics, the task-based detectability of a high-contrast and high spatial resolution imaging task was used as an example to assess the performance of each noise reduction scheme. To qualitatively assess the impact of these noise reduction schemes on image, an anthropomorphic head phantom was also scanned on the benchtop CT system and processed with the above noise reduction schemes to generate images for demonstration., Results: Our results demonstrated the following major findings: (a) CT number bias can be significantly reduced when the noise reduction schemes are implemented in the raw counts domain; CT number bias cannot be reduced when these noise reduction schemes are implemented either in the reconstructed image domain or in the log-processed sinogram domain. (b) The extent of CT number bias reduction is dependent on both the material composition and noise reduction parameters. (c) The overall impact of the noise reduction schemes can be studied using the task-based detectability analysis framework and this framework can be used to select the appropriate parameters in each noise reduction scheme to optimize the performance for a given imaging task., Conclusions: Noise reduction schemes can be used to considerably reduce CT number bias when they are implemented in the raw counts domain; however, their application cannot be arbitrarily extended to either the log-processed sinogram data domain or image domain. Trade-offs between bias reduction and overall image quality must be studied for an optimal performance of a given imaging task., (© 2019 American Association of Physicists in Medicine.)

Purpose: The CT number accuracy, that is, CT number bias, plays an important role in clinical diagnosis. When strategies to reduce radiation dose are discussed, it is important to make sure that the CT number bias is controlled within an acceptable range. The purpose of this paper was to investigate the dependence of CT number bias on radiation dose level and on image contrast (i.e., the difference in CT number between the ROI and the background) in Computed Tomography (CT)., Methods: A lesion-background model was introduced to theoretically study how the CT number bias changes with radiation exposure level and with CT number contrast when a simple linear reconstruction algorithm such as filtered backprojection (FBP) is used. The theoretical results were validated with experimental studies using a benchtop CT system equipped with a photon-counting detector (XC-HYDRA FX50, XCounter AB, Sweden) and a clinical diagnostic MDCT scanner (Discovery CT750 HD, GE Healthcare, Waukesha, WI, USA) equipped with an energy-integrating detector. The Catphan phantom (Catphan 600, the Phantom Laboratory, Salem, NY, USA) was scanned at different mAs levels and 50 scans were performed for each mAs. The bias of CT number was evaluated for each combination of mAs and ROIs with different contrast levels. An anthropomorphic phantom (ATOM 10-year-old phantom, Model 706, CIRS Inc. Norfolk, VA, USA) with much more heterogeneous object content was used to test the applicability of the theory to the more general image object cases., Results: Both theoretical and experimental studies showed that the CT number bias is inversely proportional to the radiation exposure level yet linearly dependent on the CT number contrast between the lesion and the background, that is, Bias ( μ ^ 1 FBP ) = α mAs ( 1 + β Δ H U ) ., Conclusions: The quantitative accuracy of CT numbers can be problematic and thus needs some extra attention when radiation dose is reduced. In this work, we showed that the bias of the FBP reconstruction increases as mAs is reduced; both positive and negative bias can be observed depending on the contrast difference between a targeted ROI and its surrounding background tissues., (© 2018 American Association of Physicists in Medicine.)

Purpose: Although a variety of mathematical observer models have been developed to predict human observer performance for low contrast lesion detection tasks, their predictive power for high contrast and high spatial resolution discrimination imaging tasks, including those in CT bone imaging, could be limited. The purpose of this work was to develop a modified observer model that has improved correlation with human observer performance for these tasks., Methods: The proposed observer model, referred to as the modified ideal observer model (MIOM), uses a weight function to penalize components in the task function that have less contribution to the actual human observer performance for high contrast and high spatial resolution discrimination tasks. To validate MIOM, both human observer and observer model studies were performed, each using exactly the same CT imaging task [discrimination of a connected component in a high contrast (1000 HU) high spatial resolution bone fracture model (0.3 mm)] and experimental CT image data. For the human observer studies, three physicist observers rated the connectivity of the fracture model using a five-point Likert scale; for the observer model studies, a total of five observer models, including both conventional models and the proposed MIOM, were used to calculate the discrimination capability of the CT images in resolving the connected component. Images used in the studies encompassed nine different reconstruction kernels. Correlation between human and observer model performance for these kernels were quantified using the Spearman rank correlation coefficient (ρ). After the validation study, an example application of MIOM was presented, in which the observer model was used to select the optimal reconstruction kernel for a High-Resolution (Hi-Res, GE Healthcare) CT scan technique., Results: The performance of the proposed MIOM correlated well with that of the human observers with a Spearman rank correlation coefficient ρ of 0.88 (P = 0.003). In comparison, the value of ρ was 0.05 (P = 0.904) for the ideal observer, 0.05 (P = 0.904) for the non-prewhitening observer, -0.18 (P = 0.634) for the non-prewhitening observer with eye filter and internal noise, and 0.30 (P = 0.427) for the prewhitening observer with eye filter and internal noise. Using the validated MIOM, the optimal reconstruction kernel for the Hi-Res mode to perform high spatial resolution and high contrast discrimination imaging tasks was determined to be the HD Ultra kernel at the center of the scan field of view (SFOV), or the Lung kernel at the peripheral region of the SFOV. This result was consistent with visual observations of nasal CT images of an in vivo canine subject., Conclusion: Compared with other observer models, the proposed modified ideal observer model provides significantly improved correlation with human observers for high contrast and high spatial resolution CT imaging tasks., (© 2017 American Association of Physicists in Medicine.)

Purpose: The aim of this study was to assess the effect of denoising on objective heterogeneity scores and its diagnostic capability for the diagnosis of angiomyolipoma (AML) and renal cell carcinoma (RCC)., Materials and Methods: A total of 158 resected renal masses ≤4 cm [98 clear cell (cc) RCCs, 36 papillary (pap)-RCCs, and 24 AMLs] from 139 patients were evaluated. A representative contrast-enhanced computed tomography (CT) image for each mass was selected by a genitourinary radiologist. A largest possible region of interest was drawn on each mass by the radiologist, from which three objective heterogeneity indices were calculated: standard deviation (SD), entropy (Ent), and uniformity (Uni). Objective heterogeneity indices were also calculated after images were processed with a denoising algorithm (non-local means) at three strengths: weak, medium, and strong. Two genitourinary radiologists also subjectively scored each mass independently using a three-point scale (1-3; with 1 the least and 3 the most heterogeneous), which were added to represent the final subjective heterogeneity score of each mass. Heterogeneity scores were compared among mass types, and area under the ROC curve (AUC) was calculated., Results: For all heterogeneity indices, cc-RCC was significantly more heterogeneous than pap-RCC and AML (p < 0.001), but no significant difference was found between pap-RCC and AML (p > 0.01). For cc-RCC and pap-RCC differentiation, AUCs were 0.91, 0.81, 0.78, and 0.78 for the subjective score, SD, Ent, and Uni, respectively, using original images. The corresponding AUC values were 0.84, 0.74, 0.79, and 0.80 for differentiation of AML and cc-RCC. Noise reduction at weak setting improves AUC values by 0.03, 0.05, and 0.05 for SD, entropy, and uniformity for differentiation of cc-RCC from pap-RCC. Further increase of filtering strength did not improve AUC values. For differentiation of AML vs. cc-RCC, the AUC values stayed relatively flat using the noise reduction technique at different strengths for all three indices., Conclusions: Both subjective and objective heterogeneity indices can differentiate cc-RCC from pap-RCC and AML. Noise reduction improved differentiation of cc-RCC from pap-RCC, but not differentiation of AML from cc-RCC.

Purpose: The introduction of a High-Resolution (Hi-Res) scan mode and another associated option that combines Hi-Res mode with the so-called High Definition (HD) reconstruction kernels (referred to as a Hi-Res/HD mode in this paper) in some multi-detector CT (MDCT) systems offers new opportunities to increase spatial resolution for some clinical applications that demand high spatial resolution. The purpose of this work was to quantify the in-plane spatial resolution along both the radial direction and tangential direction for the Hi-Res and Hi-Res/HD scan modes at different off-center positions., Methods: A technique was introduced and validated to address the signal saturation problem encountered in the attempt to quantify spatial resolution for the Hi-Res and Hi-Res/HD scan modes. Using the proposed method, the modulation transfer functions (MTFs) of a 64-slice MDCT system (Discovery CT750 HD, GE Healthcare) equipped with both Hi-Res and Hi-Res/HD modes were measured using a metal bead at nine different off-centered positions (0-16 cm with a step size of 2 cm); at each position, both conventional scans and Hi-Res scans were performed. For each type of scan and position, 80 repeated acquisitions were performed to reduce noise induced uncertainties in the MTF measurements. A total of 15 reconstruction kernels, including eight conventional kernels and seven HD kernels, were used to reconstruct CT images of the bead. An ex vivo animal study consisting of a bone fracture model was performed to corroborate the MTF results, as the detection of this high-contrast and high frequency task is predominantly determined by spatial resolution. Images of this animal model generated by different scan modes and reconstruction kernels were qualitatively compared with the MTF results., Results: At the centered position, the use of Hi-Res mode resulted in a slight improvement in the MTF; each HD kernel generated higher spatial resolution than its counterpart conventional kernel. However, the MTF along the tangential direction of the scan field of view (SFOV) was significantly degraded at off-centered positions, yet the combined Hi-Res/HD mode reduced this azimuthal MTF degradation. Images of the animal bone fracture model confirmed the improved spatial resolution at the off-centered positions through the use of the Hi-Res mode and HD kernels., Conclusions: The Hi-Res/HD scan improve spatial resolution of MDCT systems at both centered and off-centered positions.

Objective: The purpose of this study was to evaluate if small (< 4 cm) angiomyolipoma without visible fat can be differentiated from renal cell carcinoma (RCC) using contrast-enhanced CT alone and using unenhanced and contrast-enhanced CT in combination., Materials and Methods: Twenty-three patients with 24 angiomyolipomas without visible fat and 130 patients with 148 RCCs underwent unenhanced and contrast-enhanced CT. Demographic data and size, shape, CT attenuation, and heterogeneity (entropy and subjective score) of the renal mass on unenhanced CT and contrast-enhanced CT were recorded. Multivariate logistic regression models were constructed for parameters obtained by contrast-enhanced CT alone and by both unenhanced and contrast-enhanced CT. Demographic data and size and shape of renal mass were used in each model. Sensitivity and specificity were calculated., Results: Logistic regression model from contrast-enhanced CT data included sex, percentage of exophytic growth, entropy, and CT attenuation on contrast-enhanced CT. Model from both unenhanced and contrast-enhanced CT data included age, sex, short-axis diameter, percentage of exophytic growth, lesion-to-kidney CT attenuation difference on unenhanced CT, and CT attenuation on contrast-enhanced CT. The contrast-enhanced CT-based model and combined unenhanced and contrast-enhanced CT-based model differentiated angiomyolipoma from RCC with sensitivity and specificity of 42% and 98% versus 50% and 98%, respectively., Conclusion: Combinations of various CT and demographic findings allowed differentiation of angiomyolipoma from RCC.

Objective: The purpose of this study was to evaluate whether small (< 4 cm) oncocytomas can be differentiated from renal cell carcinomas (RCCs) on biphasic contrast-enhanced CT., Materials and Methods: Forty-three patients with 53 oncocytomas and 123 patients with 128 RCCs (24 papillary subtype and 104 clear cell and other subtypes) who underwent biphasic contrast-enhanced CT were included in the study. Patient demographics and CT tumor characteristics were evaluated in each case. A multinomial logistic regression model was then constructed for differentiating oncocytoma from clear cell and other subtype RCCs, oncocytoma from papillary RCCs, and clear cell and other subtype RCCs from papillary RCCs. The probability of each group was calculated from the model. Diagnostic performance among three pairwise diagnoses and between oncocytoma and any RCC (clear cell and other subtypes and papillary) were assessed by AUC values., Results: Patient age, tumor CT attenuation values and skewness (i.e., histogram analysis of CT values) in both the corticomedullary and nephrographic phases, and subjective tumor heterogeneity were statistically significant variables in the multinomial logistic regression analysis. The logistic regression model using the variables yielded AUCs of 0.82, 0.95, 0.91, and 0.84 for differentiating oncocytomas from clear cell and other subtype RCCs, oncocytomas from papillary RCCs, clear cell and other subtype RCCs from papillary RCCs, and oncocytomas from any RCC (clear cell and other subtypes and papillary), respectively., Conclusion: A combination of imaging features on biphasic CT, including tumor CT attenuation values and tumor texture (heterogeneity and skewness), can help differentiate oncocytoma from RCC.

Purpose: For a given imaging task and patient size, the optimal selection of x-ray tube potential (kV) and tube current-rotation time product (mAs) is pivotal in achieving the maximal radiation dose reduction while maintaining the needed diagnostic performance. Although contrast-to-noise (CNR)-based strategies can be used to optimize kV/mAs for computed tomography (CT) imaging systems employing the linear filtered backprojection (FBP) reconstruction method, a more general framework needs to be developed for systems using the nonlinear statistical model-based iterative reconstruction (MBIR) method. The purpose of this paper is to present such a unified framework for the optimization of kV/mAs selection for both FBP- and MBIR-based CT systems., Methods: The optimal selection of kV and mAs was formulated as a constrained optimization problem to minimize the objective function, Dose(kV,mAs), under the constraint that the achievable detectability index d'(kV,mAs) is not lower than the prescribed value of d'R for a given imaging task. Since it is difficult to analytically model the dependence of d' on kV and mAs for the highly nonlinear MBIR method, this constrained optimization problem is solved with comprehensive measurements of Dose(kV,mAs) and d'(kV,mAs) at a variety of kV-mAs combinations, after which the overlay of the dose contours and d' contours is used to graphically determine the optimal kV-mAs combination to achieve the lowest dose while maintaining the needed detectability for the given imaging task. As an example, d' for a 17 mm hypoattenuating liver lesion detection task was experimentally measured with an anthropomorphic abdominal phantom at four tube potentials (80, 100, 120, and 140 kV) and fifteen mA levels (25 and 50-700) with a sampling interval of 50 mA at a fixed rotation time of 0.5 s, which corresponded to a dose (CTDIvol) range of [0.6, 70] mGy. Using the proposed method, the optimal kV and mA that minimized dose for the prescribed detectability level of d'R=16 were determined. As another example, the optimal kV and mA for an 8 mm hyperattenuating liver lesion detection task were also measured using the developed framework. Both an in vivo animal and human subject study were used as demonstrations of how the developed framework can be applied to the clinical work flow., Results: For the first task, the optimal kV and mAs were measured to be 100 and 500, respectively, for FBP, which corresponded to a dose level of 24 mGy. In comparison, the optimal kV and mAs for MBIR were 80 and 150, respectively, which corresponded to a dose level of 4 mGy. The topographies of the iso-d' map and the iso-CNR map were the same for FBP; thus, the use of d'- and CNR-based optimization methods generated the same results for FBP. However, the topographies of the iso-d' and iso-CNR map were significantly different in MBIR; the CNR-based method overestimated the performance of MBIR, predicting an overly aggressive dose reduction factor. For the second task, the developed framework generated the following optimization results: for FBP, kV = 140, mA = 350, dose = 37.5 mGy; for MBIR, kV = 120, mA = 250, dose = 18.8 mGy. Again, the CNR-based method overestimated the performance of MBIR. Results of the preliminary in vivo studies were consistent with those of the phantom experiments., Conclusions: A unified and task-driven kV/mAs optimization framework has been developed in this work. The framework is applicable to both linear and nonlinear CT systems such as those using the MBIR method. As expected, the developed framework can be reduced to the conventional CNR-based kV/mAs optimization frameworks if the system is linear. For MBIR-based nonlinear CT systems, however, the developed task-based kV/mAs optimization framework is needed to achieve the maximal dose reduction while maintaining the desired diagnostic performance.

Objective: The purpose of this article is to assess the effect of radiation dose reduction in dual-energy CT (DECT) on the performance of renal stone characterization using a patient cohort., Materials and Methods: CT data from 39 unenhanced DECT examinations performed for stone characterization were retrospectively analyzed in this study. Reduced-dose images were simulated at 75%, 50%, and 25% of the routine dose using a previously validated noise-insertion algorithm. Differentiation between uric acid (UA) and non-UA stones was performed using a fixed cutoff value for the dual-energy ratio. ROC analysis was performed to determine optimal cutoff values and the associated sensitivity and specificity., Results: Of the 206 stones found, 43 were UA and 163 were non-UA. The mean (± SD) volume CT dose index (CTDIvol) was 16.0 ± 4.0 mGy at the 100% dose level. The mean noise in 100-kV images increased from 40.9 ± 6.8 HU at 100% dose to 46.8 ± 8.8 HU, 57.7 ± 12.5 HU, and 85.4 ± 22.9 HU at 75%, 50%, and 25% dose levels, respectively. Using the default cutoff value, for stones 10 mm(3) or larger, the sensitivity/specificity were 100.0%/98.8%, 82.8%/98.8%, and 89.3%/98.7%, at 75%, 50%, and 25% dose levels, respectively. ROC analysis showed varying optimal cutoff values at different dose levels. The sensitivity and specificity improved with use of these optimal cutoff values. Differentiation capability decreased for stones smaller than 10 mm(3)., Conclusion: At 75% of the 16-mGy routine dose, the sensitivity and specificity for differentiating UA from non-UA stones were minimally affected for stones 10 mm(3) or larger. The use of optimal cutoff values for dual-energy ratio as dose decreased (and noise increased) provided improved performance.