Feasibility of thin-slice, low noise images created using multi-kernel synthesis to replace multiple image series in head CT.

BACKGROUND: SynthesiZed Improved Resolution and Concurrent nOise reductioN (ZIRCON) is a multi-kernel synthesis method that creates a single series of thin-slice computed tomography (CT) images displaying low noise and high spatial resolution, increasing reader efficiency and minimizing partial volume averaging. PURPOSE: To compare the diagnostic performance of a single set of ZIRCON images to two routine clinical image series using conventional CT head and bone reconstruction kernels for diagnosing intracranial findings and fractures in patients with trauma or suspected acute neurologic deficit. MATERIAL AND METHODS: In total, 50 patients underwent clinically indicated head CT in the ER (15 normal, 35 abnormal cases). A non-reader neuroradiologist established the reference standard. Three neuroradiologists reviewed two routine clinical series (head and bone kernels) and a single ZIRCON series, detecting intracranial findings or fractures and rating confidence (0-100). Sensitivity, specificity, and jackknife free-response receiver operating characteristic (JAFROC) figure of merit (FOM) were compared (limit of non-inferiority: -0.10). RESULTS: ZIRCON and conventional images demonstrated comparable performance for fractures (sensitivity: 51.5% vs. 54.5%; specificity: 40.2% vs. 34.2%) and intracranial findings (sensitivity: 88.2% vs. 91.4%; specificity: 77.2% vs. 73.7%).The estimated difference of JAFROC FOM demonstrated ZIRCON non-inferiority for acute pathologies overall (0.003 [95% CI=-0.051-0.057]) and fractures (0.048 [95% CI=-0.050-0.145]) but not for intracranial findings alone (-0.024 [95% CI=-0.100-0.052]). CONCLUSION: Thin-slice, low noise, and high spatial resolution images can be created to display intracranial findings and fractures replacing multiple images series in head CT with similar performance. Future studies in more patients and further algorithmic development are warranted.

Small (< 4 cm) Renal Mass: Differentiation of Oncocytoma From Renal Cell Carcinoma on Biphasic Contrast-Enhanced CT.

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.

Small (< 4 cm) Renal Masses: Differentiation of Angiomyolipoma Without Visible Fat From Renal Cell Carcinoma Using Unenhanced and Contrast-Enhanced CT.

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.

Multiple Kernel Synthesis of Head CT Using a Task-Based Loss Function.

In CT imaging of the head, multiple image series are routinely reconstructed with different kernels and slice thicknesses. Reviewing the redundant information is an inefficient process for radiologists. We address this issue with a convolutional neural network (CNN)-based technique, synthesiZed Improved Resolution and Concurrent nOise reductioN (ZIRCON), that creates a single, thin, low-noise series that combines the favorable features from smooth and sharp head kernels. ZIRCON uses a CNN model with an autoencoder U-Net architecture that accepts two input channels (smooth- and sharp-kernel CT images) and combines their salient features to produce a single CT image. Image quality requirements are built into a task-based loss function with a smooth and sharp loss terms specific to anatomical regions. The model is trained using supervised learning with paired routine-dose clinical non-contrast head CT images as training targets and simulated low-dose (25%) images as training inputs. One hundred unique de-identified clinical exams were used for training, ten for validation, and ten for testing. Visual comparisons and contrast measurements of ZIRCON revealed that thinner slices and the smooth-kernel loss function improved gray-white matter contrast. Combined with lower noise, this increased visibility of small soft-tissue features that would be otherwise impaired by partial volume averaging or noise. Line profile analysis showed that ZIRCON images largely retained sharpness compared to the sharp-kernel input images. ZIRCON combined desirable image quality properties of both smooth and sharp input kernels into a single, thin, low-noise series suitable for both brain and skull imaging.

Adaptation of a channelized Hotelling observer model to accommodate anthropomorphic backgrounds and moving test objects in X-ray angiography.

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.

MRI Radiomics for Assessment of Molecular Subtype, Pathological Complete Response, and Residual Cancer Burden in Breast Cancer Patients Treated With Neoadjuvant Chemotherapy.

RATIONALE AND OBJECTIVES: There are limited data on pretreatment imaging features that can predict response to neoadjuvant chemotherapy (NAC). To extract volumetric pretreatment MRI radiomics features and assess corresponding associations with breast cancer molecular subtypes, pathological complete response (pCR), and residual cancer burden (RCB) in patients treated with NAC. MATERIALS AND METHODS: In this IRB-approved study, clinical and pretreatment MRI data from patients with biopsy-proven breast cancer who received NAC between September 2009 and July 2016 were retrospectively analyzed. Tumors were manually identified and semi-automatically segmented on first postcontrast images. Morphological and three-dimensional textural features were computed, including unfiltered and filtered image data, with spatial scaling factors (SSF) of 2, 4, and 6 mm. Wilcoxon rank-sum tests and area under the receiver operating characteristic curve were used for statistical analysis. RESULTS: Two hundred and fifty nine patients with unilateral breast cancer, including 73 (28.2%) HER2+, 112 (43.2%) luminal, and 74 (28.6%) triple negative breast cancers (TNBC), were included. There was a significant difference in the median volume (p = 0.008), median longest axial tumor diameter (p = 0.009), and median longest volumetric diameter (p = 0.01) among tumor subtypes. There was also a significant difference in minimum signal intensity and entropy among the tumor subtypes with SSF = 4 mm (p = 0.009 and p = 0.02 respectively) and SSF = 6 mm (p = 0.007 and p < 0.001 respectively). Additionally, sphericity (p = 0.04) in HER2+ tumors and entropy with SSF = 2, 4, 6 mm (p = 0.004, 0.02, 0.047 respectively) in luminal tumors were significantly associated with pCR. Multiple features demonstrated significant association (p < 0.05) with pCR in TNBC and with RCB in luminal tumors and TNBC, with standard deviation of intensity with SSF = 6 mm achieving the highest AUC (AUC = 0.734) for pCR in TNBC. CONCLUSION: MRI radiomics features are associated with different molecular subtypes of breast cancer, pCR, and RCB. These features may be noninvasive imaging biomarkers to identify cancer subtype and predict response to NAC.

Task-specific efficient channel selection and bias management for Gabor function channelized Hotelling observer model for the assessment of x-ray angiography system performance.

PURPOSE: Channelized Hotelling observer (CHO) models have been implemented to assess imaging performance in x-ray angiography systems. While current methods are appropriate for assessing unprocessed images of moving test objects upon uniform-exposure backgrounds, they are inadequate for assessing conditions which more appropriately mimic clinical imaging conditions including the combination of moving test objects, complex anthropomorphic backgrounds, and image processing. In support of this broad goal, the purpose of this work was to develop theory and methods to automatically select a subset of task-specific efficient Gabor channels from a task-generic Gabor channel base set. Also, previously described theory and methods to manage detectability index (d') bias due to nonrandom temporal variations in image electronic noise will be revisited herein. METHODS: Starting with a base set of 96 Gabor channels, backward elimination of channels was used to automatically identify an "efficient" channel subset which reduced the number of channels retained in the subset while maintaining the magnitude of the d' estimate. The concept of a pixelwise Hotelling observer (PHO) model was introduced and similarly implemented to assess the performance of the efficient-channel CHO model. Bias in d' estimates arising from temporally variable nonstationary noise was modeled as a bivariate probability density function for normal distributions, where one variable corresponds to the signal from the test object and the other variable corresponds to the signal from temporally variable nonstationary noise. Theory and methods were tested on uniform-exposure unprocessed angiography images with detector target dose (DTD) of 6, 18, and 120 nGy containing static disk-shaped test objects with diameter in the range of 0.5 to 4 mm. RESULTS: Considering all DTD levels and test object sizes, the proposed method reduced the number of Gabor channels in the final subset by 63-82% compared to the original 96 Gabor channel base set, while maintaining a mean relative performance ( ( d CHO ' / d PHO ' ) × 100 % ) of 95%  ±  4% with respect to the reference PHO model. Experimental results demonstrated that the bivariate approach to account for bias due to temporally variable nonstationary noise resulted in improved correlation between the CHO and PHO models as compared to a previously proposed univariate approach. CONCLUSIONS: Computationally efficient backward elimination can be used to select an efficient subset of Gabor channels from an initial channel base set without substantially compromising the magnitude of the d' estimate. Bias due to temporally variable nonstationary noise can be modeled through a bivariate approach leading to an improved unbiased estimate of d'.

Impact of noise reduction schemes on quantitative accuracy of CT numbers.

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.

Low-dose cone-beam CT via raw counts domain low-signal correction schemes: Performance assessment and task-based parameter optimization (Part II. Task-based parameter optimization).

PURPOSE: Different low-signal correction (LSC) methods have been shown to efficiently reduce noise streaks and noise level in CT to provide acceptable images at low-radiation dose levels. These methods usually result in CT images with highly shift-variant and anisotropic spatial resolution and noise, which makes the parameter optimization process highly nontrivial. The purpose of this work was to develop a local task-based parameter optimization framework for LSC methods. METHODS: Two well-known LSC methods, the adaptive trimmed mean (ATM) filter and the anisotropic diffusion (AD) filter, were used as examples to demonstrate how to use the task-based framework to optimize filter parameter selection. Two parameters, denoted by the set P, for each LSC method were included in the optimization problem. For the ATM filter, these parameters are the low- and high-signal threshold levels pl and ph ; for the AD filter, the parameters are the exponents δ and γ in the brightness gradient function. The detectability index d' under the non-prewhitening (NPW) mathematical observer model was selected as the metric for parameter optimization. The optimization problem was formulated as an unconstrained optimization problem that consisted of maximizing an objective function d'(P), where i and j correspond to the i-th imaging task and j-th spatial location, respectively. Since there is no explicit mathematical function to describe the dependence of d' on the set of parameters P for each LSC method, the optimization problem was solved via an experimentally measured d' map over a densely sampled parameter space. In this work, three high-contrast-high-frequency discrimination imaging tasks were defined to explore the parameter space of each of the LSC methods: a vertical bar pattern (task I), a horizontal bar pattern (task II), and a multidirectional feature (task III). Two spatial locations were considered for the analysis, a posterior region-of-interest (ROI) located within the noise streaks region and an anterior ROI, located further from the noise streaks region. Optimal results derived from the task-based detectability index metric were compared to other operating points in the parameter space with different noise and spatial resolution trade-offs. RESULTS: The optimal operating points determined through the d' metric depended on the interplay between the major spatial frequency components of each imaging task and the highly shift-variant and anisotropic noise and spatial resolution properties associated with each operating point in the LSC parameter space. This interplay influenced imaging performance the most when the major spatial frequency component of a given imaging task coincided with the direction of spatial resolution loss or with the dominant noise spatial frequency component; this was the case of imaging task II. The performance of imaging tasks I and III was influenced by this interplay in a smaller scale than imaging task II, since the major frequency component of task I was perpendicular to imaging task II, and because imaging task III did not have strong directional dependence. For both LSC methods, there was a strong dependence of the overall d' magnitude and shape of the contours on the spatial location within the phantom, particularly for imaging tasks II and III. The d' value obtained at the optimal operating point for each spatial location and imaging task was similar when comparing the LSC methods studied in this work. CONCLUSIONS: A local task-based detectability framework to optimize the selection of parameters for LSC methods was developed. The framework takes into account the potential shift-variant and anisotropic spatial resolution and noise properties to maximize the imaging performance of the CT system. Optimal parameters for a given LSC method depend strongly on the spatial location within the image object.

Low-dose cone-beam CT via raw counts domain low-signal correction schemes: Performance assessment and task-based parameter optimization (Part I: Assessment of spatial resolution and noise performance).

PURPOSE: Low-signal correction (LSC) in the raw counts domain has been shown to effectively reduce noise streaks in CT because the data inconsistency associated with photon-starved regions may be mitigated prior to the log transformation step. However, a systematic study of the performance of these raw data correction methods is still missing in literature. The purpose of this work was to provide such a systematic study for two well-known low-signal correction schemes using either the adaptive trimmed mean (ATM) filter or the anisotropic diffusion (AD) filter in the raw counts domain. METHODS: Image data were acquired experimentally using an anthropomorphic chest phantom and a benchtop cone-beam CT (CBCT) imaging system. Phantom scans were repeated 50 times at a reduced dose level of 0.5 mGy and a reference level of 1.9 mGy. The measured raw counts at 0.5 mGy underwent LSC using the ATM and AD filters. Two relevant parameters were identified for each filter and approximately one hundred operating points in each parameter space were analyzed. Following LSC and log transformation, FDK reconstruction was performed for each case. Noise and spatial resolution properties were assessed across the parameter spaces that define each LSC filter; the results were summarized through 2D contour maps to better understand the trade-offs between these competing image quality features. 2D noise power spectrum (NPS) and modulation transfer function (MTF) were measured locally at two spatial locations in the field-of-view (FOV): a posterior region contaminated by noise streaks and an anterior region away from noise streaks. An isotropy score metric was introduced to characterize the directional dependence of the NPS and MTF (viz., ϵNPS and ϵMTF , respectively), with a range from 0 for highly anisotropic to 1 for perfectly isotropic. The noise magnitude and coarseness were also measured. RESULTS: (a) Both the ATM and AD LSC methods were successful in reducing noise streaks, but their noise and spatial resolution properties were found to be highly anisotropic and shift-variant. (b) NPS isotropy scores in the posterior region were generally improved from ϵNPS = 0.09 for the images without LSC to the range ϵNPS = (0.11, 0.67) for ATM and ϵNPS = (0.06, 0.67) for AD, depending on the filter parameters used. (c) The noise magnitude was reduced across the parameter space of either LSC filter whenever a change along the axis of the controlling parameter led to stronger raw data filtration. Changes in noise magnitude were inversely related to changes in spatial resolution along the direction perpendicular to the streaks. No correlation was found, however, between the contour maps of noise magnitude and the NPS isotropy. (d) Both filters influenced the noise coarseness anisotropically, with coarser noise occurring along directions perpendicular to the noise streaks. The anisotropic noise coarseness was intrinsically and directly related to resolution losses in a given direction: coarseness plots mimic the topography of the 2D MTF, i.e., the coarser the noise, the lower the resolution. CONCLUSIONS: Both AD and ATM LSC schemes enable low-dose CBCT imaging. However, it was found that noise magnitude and overall spatial resolution vary considerably across the parameter space for each filter, and more importantly these image quality features are highly anisotropic and shift-variant.

Quantitative accuracy of CT numbers: Theoretical analyses and experimental studies.

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.

Modified ideal observer model (MIOM) for high-contrast and high-spatial resolution CT imaging tasks.

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.

Subjective and objective heterogeneity scores for differentiating small renal masses using contrast-enhanced CT.

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.

Can conclusions drawn from phantom-based image noise assessments be generalized to in vivo studies for the nonlinear model-based iterative reconstruction method?

PURPOSE: Phantom-based objective image quality assessment methods are widely used in the medical physics community. For a filtered backprojection (FBP) reconstruction-based linear or quasilinear imaging system, the use of this methodology is well justified. Many key image quality metrics acquired with phantom studies can be directly applied to in vivo human subject studies. Recently, a variety of image quality metrics have been investigated for model-based iterative image reconstruction (MBIR) methods and several novel characteristics have been discovered in phantom studies. However, the following question remains unanswered: can certain results obtained from phantom studies be generalized to in vivo animal studies and human subject studies? The purpose of this paper is to address this question. METHODS: One of the most striking results obtained from phantom studies is a novel power-law relationship between noise variance of MBIR (σ(2)) and tube current-rotation time product (mAs): σ(2) ∝ (mAs)(-0.4) [K. Li et al., "Statistical model based iterative reconstruction (MBIR) in clinical CT systems: Experimental assessment of noise performance," Med. Phys. 41, 041906 (15pp.) (2014)]. To examine whether the same power-law works for in vivo cases, experimental data from two types of in vivo studies were analyzed in this paper. All scans were performed with a 64-slice diagnostic CT scanner (Discovery CT750 HD, GE Healthcare) and reconstructed with both FBP and a MBIR method (Veo, GE Healthcare). An Institutional Animal Care and Use Committee-approved in vivo animal study was performed with an adult swine at six mAs levels (10-290). Additionally, human subject data (a total of 110 subjects) acquired from an IRB-approved clinical trial were analyzed. In this clinical trial, a reduced-mAs scan was performed immediately following the standard mAs scan; the specific mAs used for the two scans varied across human subjects and were determined based on patient size and clinical indications. The measurements of σ(2) were performed at different mAs by drawing regions-of-interest (ROIs) in the liver and the subcutaneous fat. By applying a linear least-squares regression, the ß values in the power-law relationship σ(2) ∝ (mAs)(-ß) were measured for the in vivo data and compared with the value found in phantom experiments. RESULTS: For the in vivo swine study, an exponent of ß = 0.43 was found for MBIR, and the coefficient of determination (R(2)) for the corresponding least-squares power-law regression was 0.971. As a reference, the ß and R(2) values for FBP were found to be 0.98 and 0.997, respectively, from the same study, which are consistent with the well-known σ(2) ∝ (mAs)(-1.0) relationship for linear CT systems. For the human subject study, the measured ß values for the MBIR images were 0.41 ± 0.12 in the liver and 0.37 ± 0.12 in subcutaneous fat. In comparison, the ß values for the FBP images were 1.04 ± 0.10 in the liver and 0.97 ± 0.12 in subcutaneous fat. The ß values of MBIR and FBP obtained from the in vivo studies were found to be statistically equivalent to the corresponding ß values from the phantom study within an equivalency interval of [ - 0.1, 0.1] (p < 0.05); across MBIR and FBP, the difference in ß was statistically significant (p < 0.05). CONCLUSIONS: Despite the nonlinear nature of the MBIR method, the power-law relationship, σ(2) ∝ (mAs)(-0.4), found from phantom studies can be applied to in vivo animal and human subject studies.

Impact of bowtie filter and object position on the two-dimensional noise power spectrum of a clinical MDCT system.

PURPOSE: Noise characteristics of clinical multidetector CT (MDCT) systems can be quantified by the noise power spectrum (NPS). Although the NPS of CT has been extensively studied in the past few decades, the joint impact of the bowtie filter and object position on the NPS has not been systematically investigated. This work studies the interplay of these two factors on the two dimensional (2D) local NPS of a clinical CT system that uses the filtered backprojection algorithm for image reconstruction. METHODS: A generalized NPS model was developed to account for the impact of the bowtie filter and image object location in the scan field-of-view (SFOV). For a given bowtie filter, image object, and its location in the SFOV, the shape and rotational symmetries of the 2D local NPS were directly computed from the NPS model without going through the image reconstruction process. The obtained NPS was then compared with the measured NPSs from the reconstructed noise-only CT images in both numerical phantom simulation studies and experimental phantom studies using a clinical MDCT scanner. The shape and the associated symmetry of the 2D NPS were classified by borrowing the well-known atomic spectral symbols s, p, and d, which correspond to circular, dumbbell, and cloverleaf symmetries, respectively, of the wave function of electrons in an atom. Finally, simulated bar patterns were embedded into experimentally acquired noise backgrounds to demonstrate the impact of different NPS symmetries on the visual perception of the object. RESULTS: (1) For a central region in a centered cylindrical object, an s-wave symmetry was always present in the NPS, no matter whether the bowtie filter was present or not. In contrast, for a peripheral region in a centered object, the symmetry of its NPS was highly dependent on the bowtie filter, and both p-wave symmetry and d-wave symmetry were observed in the NPS. (2) For a centered region-ofinterest (ROI) in an off-centered object, the symmetry of its NPS was found to be different from that of a peripheral ROI in the centered object, even when the physical positions of the two ROIs relative to the isocenter were the same. (3) The potential clinical impact of the highly anisotropic NPS, caused by the interplay of the bowtie filter and position of the image object, was highlighted in images of specific bar patterns oriented at different angles. The visual perception of the bar patterns was found to be strongly dependent on their orientation. CONCLUSIONS: The NPS of CT depends strongly on the bowtie filter and object position. Even if the location of the ROI with respect to the isocenter is fixed, there can be different symmetries in the NPS, which depend on the object position and the size of the bowtie filter. For an isolated off-centered object, the NPS of its CT images cannot be represented by the NPS measured from a centered object.

Hi-Res scan mode in clinical MDCT systems: Experimental assessment of spatial resolution performance.

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.

Statistical model based iterative reconstruction in clinical CT systems. Part III. Task-based kV/mAs optimization for radiation dose reduction.

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.

Influence of radiation dose and reconstruction algorithm in MDCT assessment of airway wall thickness: A phantom study.

PURPOSE: Wall thickness (WT) is an airway feature of great interest for the assessment of morphological changes in the lung parenchyma. Multidetector computed tomography (MDCT) has recently been used to evaluate airway WT, but the potential risk of radiation-induced carcinogenesis-particularly in younger patients-might limit a wider use of this imaging method in clinical practice. The recent commercial implementation of the statistical model-based iterative reconstruction (MBIR) algorithm, instead of the conventional filtered back projection (FBP) algorithm, has enabled considerable radiation dose reduction in many other clinical applications of MDCT. The purpose of this work was to study the impact of radiation dose and MBIR in the MDCT assessment of airway WT. METHODS: An airway phantom was scanned using a clinical MDCT system (Discovery CT750 HD, GE Healthcare) at 4 kV levels and 5 mAs levels. Both FBP and a commercial implementation of MBIR (Veo(TM), GE Healthcare) were used to reconstruct CT images of the airways. For each kV-mAs combination and each reconstruction algorithm, the contrast-to-noise ratio (CNR) of the airways was measured, and the WT of each airway was measured and compared with the nominal value; the relative bias and the angular standard deviation in the measured WT were calculated. For each airway and reconstruction algorithm, the overall performance of WT quantification across all of the 20 kV-mAs combinations was quantified by the sum of squares (SSQs) of the difference between the measured and nominal WT values. Finally, the particular kV-mAs combination and reconstruction algorithm that minimized radiation dose while still achieving a reference WT quantification accuracy level was chosen as the optimal acquisition and reconstruction settings. RESULTS: The wall thicknesses of seven airways of different sizes were analyzed in the study. Compared with FBP, MBIR improved the CNR of the airways, particularly at low radiation dose levels. For FBP, the relative bias and the angular standard deviation of the measured WT increased steeply with decreasing radiation dose. Except for the smallest airway, MBIR enabled significant reduction in both the relative bias and angular standard deviation of the WT, particularly at low radiation dose levels; the SSQ was reduced by 50%-96% by using MBIR. The optimal reconstruction algorithm was found to be MBIR for the seven airways being assessed, and the combined use of MBIR and optimal kV-mAs selection resulted in a radiation dose reduction of 37%-83% compared with a reference scan protocol with a dose level of 1 mGy. CONCLUSIONS: The quantification accuracy of airway WT is strongly influenced by radiation dose and reconstruction algorithm. The MBIR algorithm potentially allows the desired WT quantification accuracy to be achieved with reduced radiation dose, which may enable a wider clinical use of MDCT for the assessment of airway WT, particularly for younger patients who may be more sensitive to exposures with ionizing radiation.