Investigation of a potential upstream harmonization based on image appearance matching to improve radiomics features robustness: a phantom study

Objective. Radiomics is a promising valuable analysis tool consisting in extracting quantitative information from medical images. However, the extracted radiomics features are too sensitive to variations in used image acquisition and reconstruction parameters. This limited robustness hinders the generalizable validity of radiomics-assisted models. Our aim is to investigate a possible harmonization strategy based on matching image quality to improve feature robustness.Approach.We acquired CT scans of a phantom with two scanners across different dose levels and percentages of Iterative Reconstruction algorithms. The detectability index was used as a comprehensive task-based image quality metric. A statistical analysis based on the Intraclass Correlation Coefficient was performed to determine if matching image quality/appearance could enhance the robustness of radiomics features extracted from the phantom images. Additionally, an Artificial Neural Network was trained on these features to automatically classify the scanner used for image acquisition.Main results.We found that the ICC of the features across protocols providing a similar detectability index improves with respect to the ICC of the features across protocols providing a different detectability index. This improvement was particularly noticeable in features relevant for distinguishing between scanners.Significance.This preliminary study demonstrates that a harmonization based on image quality/appearance matching could improve radiomics features robustness and heterogeneous protocols can be used to obtain a similar image appearance in terms of the detectability index. Thus protocols with a lower dose level could be selected to reduce the amount of radiation dose delivered to the patient and simultaneously obtain a more robust quantitative analysis.

Technical note: Characterization, validation, and spectral optimization of a dedicated breast CT system for contrast-enhanced imaging

BACKGROUND

The development of a new imaging modality, such as 4D dynamic contrast-enhanced dedicated breast CT (4D DCE-bCT), requires optimization of the acquisition technique, particularly within the 2D contrast-enhanced imaging modality. Given the extensive parameter space, cascade-systems analysis is commonly used for such optimization.

PURPOSE

To implement and validate a parallel-cascaded model for bCT, focusing on optimizing and characterizing system performance in the projection domain to enhance the quality of input data for image reconstruction.

METHODS

A parallel-cascaded system model of a state-of-the-art bCT system was developed and model predictions of the presampled modulation transfer function (MTF) and the normalized noise power spectrum (NNPS) were compared with empirical data collected in the projection domain. Validation was performed using the default settings of 49 kV with 1.5 mm aluminum filter and at 65 kV and 0.257 mm copper filter. A 10 mm aluminum plate was added to replicate the breast attenuation. Air kerma at the isocenter was measured at different tube current levels. Discrepancies between the measured projection domain metrics and model-predicted values were quantified using percentage error and coefficient of variation (CoV) for MTF and NNPS, respectively. The optimal filtration was for a 5 mm iodine disk detection task at 49, 55, 60, and 65 kV. The detectability index was calculated for the default aluminum filtration and for copper thicknesses ranging from 0.05 to 0.4 mm.

RESULTS

At 49 kV, MTF errors were +5.1% and -5.1% at 1 and 2 cycles/mm, respectively; NNPS CoV was 5.3% (min = 3.7%; max = 8.5%). At 65 kV, MTF errors were -0.8% and -3.2%; NNPS CoV was 13.1% (min = 11.4%; max = 16.9%). Air kerma output was linear, with 11.67 µGy/mA (R2 = 0.993) and 19.14 µGy/mA (R2 = 0.996) at 49 and 65 kV, respectively. For iodine detection, a 0.25 mm-thick copper filter at 65 kV was found optimal, outperforming the default technique by 90%.

CONCLUSION

The model accurately predicts bCT system performance, specifically in the projection domain, under varied imaging conditions, potentially contributing to the enhancement of 2D contrast-enhanced imaging in 4D DCE-bCT.

The impact of automatic tube current modulation related settings of a modern GE CT scanner on image quality and patient dose; details do matter

PURPOSE

To investigate the operation principles of the automatic tube current modulation (ATCM) of a modern GE healthcare CT scanner, and the impact of related settings on image quality and patient dose.

MATERIAL & METHODS

A dedicated phantom (Mercury 4.0) was scanned using two of the most frequently used clinical scanning protocols (chest and abdomen-pelvis). The preset protocol settings were used as starting points (reference conditions). Scan direction, scan mode (helical vs. axial), total beam width, tube potential (kVp), and ATCM settings were then modified individually to understand their impact on radiation dose and image quality. Regarding the ATCM settings, the SmartmA minimum and maximum mA limits, and the noise index (NI) values were varied. As surrogates of patient dose, the CTDIvol and DLP values of each scan were used. As surrogates of image quality were used the image noise and the detectability index (d') of five different materials (air, solid water, polystyrene, iodine, and bone) embedded in the Mercury phantom calculated with the ImQuest software.

RESULTS

The scanning direction did not have any effect on ATCM curves, unlike what has been observed in CT scanners from other manufacturers. Total beam width does matter, however, the SmartmA limit settings and kVp selection had the greatest impact on image quality and dose. It was seen that improper minimum mA limit settings practically invalidated the ATCM operation. In contrast, when full modulation was allowed without restrictions, noise standard deviation, and detectability index became much more consistent across the wide range of phantom diameters. For lower kVp settings an impressive dose reduction was observed that requires further investigation.

CONCLUSION

SmartmA is a tool that if not properly used may increase the patient doses considerably. Therefore, its settings should be carefully adjusted for each preset different clinical protocol.

A systematic task-based image quality assessment of photon-counting and energy integrating CT as a function of reconstruction kernel and phantom size

BACKGROUND

Image quality of photon-counting and energy integrating CT scanners changes with object size, dose to the object, and kernel selection.

PURPOSE

To comprehensively compare task-generic image quality of photon-counting CT (PCCT) and energy integrating CT (EICT) systems as a function of phantom size, dose, and reconstruction kernel.

METHODS

A size-variant phantom (Mercury Phantom 3.0) was used to characterize the image quality of PCCT and EICT systems as a function of object size. The phantom contained five cylinders attached by slanted tapered sections. Each cylinder contained two sections: one uniform for noise, and the other with inserts for spatial resolution and contrast measurements. The phantom was scanned on Siemens' SOMATOM Force and NAEOTOM Alpha at 1.18 and 7.51 mGy without tube current modulation. CTDIvol was matched across two systems by setting the required tube currents, else, all other acquisition settings were fixed. CT sinograms were reconstructed using FBP and iterative (ADMIRE2 - Force; QIR2 - Alpha) algorithms with Body regular (Br) kernels. Noise Power Spectrum (NPS), Task Transfer Function (TTF), contrast-to-noise ratio (CNR), and detectability index (d') for a task of identifying 2-mm disk were evaluated based on AAPM TG-233 metrology using imQuest, an open-source software package. Averaged noise frequency (fav ) and 50% cut-off frequency for TTF (f50 ) were used as scalar metrics to quantify noise texture and spatial resolution, respectively. The difference between image quality metrics' measurements was calculated as IQPCCT - IQEICT .

RESULTS

From Br40 (soft) to Br64 (sharp), f50 for air insert increased from 0.35 mm-1  ± 0.04 (standard deviation) to 0.76 mm-1  ± 0.17, 0.34 mm-1  ± 0.04 to 0.77 mm-1  ± 0.17, 0.37 mm-1  ± 0.02 to 0.95 mm-1  ± 0.17 for PCCT-T3D-QIR2, PCCT-70keV-QIR2, and EICT-ADMIRE2, respectively, when averaged over all sizes and dose levels. Similarly, from Br40 to Br64, noise magnitude increased from 10.86 HU ± 7.12 to 38.61 HU ± 18.84, 10.94 HU ± 7.08 to 38.82 HU ± 18.70, 13.74 HU ± 11.02 to 52.11 HU ± 26.22 for PCCT-T3D-QIR2, PCCT-70keV-QIR2, and EICT-ADMIRE2, respectively. The difference in fav was consistent across all sizes and dose levels. PCCT-70keV-VMI showed better consistency in contrast measurements for iodine and bone inserts than PCCT-T3D and EICT; however, PCCT-T3D had higher contrast for both inserts. From Br40 to Br64, considering all sizes and dose levels, CNR for iodine insert decreased from 52.30 ± 46.44 to 12.18 ± 10.07, 40.42 ± 33.42 to 9.48 ± 7.16, 39.94 ± 37.60 to 7.84 ± 6.67 for PCCT-T3D-QIR2, PCCT-70keV-QIR2, and EICT-ADMIRE2, respectively.

CONCLUSIONS

Both PCCT image types, T3D and 70-keV-VMI exhibited similar or better noise, contrast, CNR than EICT when comparing kernels with similar names. For 512 × 512 matrix, PCCT's sharp kernels had lower resolution than EICT's sharp kernels. For all image quality metrics, except extreme low, every dose condition had a similar set of IQ-matching kernels. It suggests that considering patient size and dose level to determine IQ-matching kernel pairs across PCCT and EICT systems may not be imperative while translating protocols, except when the signal to the detector is extremely low.

Impact of a Deep Learning-based Super-resolution Image Reconstruction Technique on High-contrast Computed Tomography: A Phantom Study

RATIONALE AND OBJECTIVES

Deep-learning-based super-resolution image reconstruction (DLSRR) is a novel image reconstruction technique that is expected to contribute to improvement in spatial resolution as well as noise reduction through learning from high-resolution computed tomography (CT). This study aims to evaluate image quality obtained with DLSRR and assess its clinical potential.

MATERIALS AND METHODS

CT images of a Mercury CT 4.0 phantom were obtained using a 320-row multi-detector scanner at tube currents of 100, 200, and 300 mA. Image data were reconstructed by filtered back projection (FBP), hybrid iterative reconstruction (HIR), model-based iterative reconstruction (MBIR), deep-learning-based image reconstruction (DLR), and DLSRR at image reconstruction strength levels of mild, standard, and strong. Noise power spectrum (NPS), task transfer function (TTF), and detectability index were calculated.

RESULTS

The magnitude of the noise-reducing effect in comparison with FBP was in the order MBIR

CONCLUSION

The present results suggest that DLSRR can achieve greater noise reduction and improved spatial resolution in the high-contrast region compared with conventional DLR and iterative reconstruction techniques.

Collapse

Key Words

• deep-learning-based super-resolution image reconstruction
• detectability index
• iterative reconstruction
• noise power spectrum
• task transfer function

Collapse

MESH Headings Collapse

Grants Collapse

The IAEA remote and automated quality control methodology for radiography and mammography

Radiography remains the most widely used imaging modality throughout the world. Additionally, while it has been demonstrated that a quality control (QC) program, especially in mammography, improves image quality, weekly technologist QC testing might be lacking even where there is clinical qualified medical physicist (CQMP) support. Therefore, the International Atomic Energy Agency (IAEA) developed simple QC phantoms that can easily be used on a regular basis (daily/weekly) for radiography and mammography. These are simple in design and use materials that are easily accessible in most parts of the world. A software application is also developed that automatically analyzes images and Digital Imaging and Communications in Medicine (DICOM) header information. It exports data to a comma-separated values (CSV) file that is read by a Microsoft Excel® spreadsheet for documentation and graphical analysis. The phantom and the software were tested in four institutions (in Costa Rica and the United States of America) both on computed radiography and direct digital mammography and radiography systems. Data were collected over a 3-year period. No corrective actions were taken on the data, but service was performed on two of the units. Results demonstrated noise that could be attributed to suboptimal placement of the phantom and incorrect data being put into the DICOM header. Preliminary evaluation of the IAEA methodology has demonstrated that it can provide meaningful QC data that are sensitive to changes in the imaging systems. Care must be taken at implementation to properly train personnel and ensure that the image data, including the DICOM header, are being correctly transmitted. The methodology gives the opportunity for a single CQMP to provide QC services even to remote sites where travel is prohibitive, and it is feasible and easy to implement.

Task-driven source-detector trajectories in cone-beam computed tomography: II. Application to neuroradiology

We apply the methodology detailed in "Task-driven source-detector trajectories in cone-beam computed tomography: I. Theory and methods" by Stayman et al. for task-driven optimization of source-detector orbits in cone-beam computed tomography (CBCT) to scenarios emulating imaging tasks in interventional neuroradiology. The task-driven imaging framework is used to optimize the CBCT source-detector trajectory by maximizing the detectability index, d ' . The approach was applied to simulated cases of endovascular embolization of an aneurysm and arteriovenous malformation and was translated to real data first using a CBCT test bench followed by implementation on an interventional robotic C-arm. Task-driven trajectories were found to generally favor higher fidelity (i.e., less noisy) views, with an average increase in d ' ranging from 7% to 28%. Visually, this resulted in improved conspicuity of particular stimuli by reducing the noise and altering the noise correlation to a form distinct from the spatial frequencies associated with the imaging task. The improvements in detectability and the demonstration of the task-driven workflow using a real interventional imaging system show the potential of the task-driven imaging framework to improve imaging performance on motorized, multiaxis C-arms in neuroradiology.

Task-driven source-detector trajectories in cone-beam computed tomography: I. Theory and methods

We develop a mathematical framework for the design of orbital trajectories that are optimal to a particular imaging task (or tasks) in advanced cone-beam computed tomography systems that have the capability of general source-detector positioning. The framework allows various parameterizations of the orbit as well as constraints based on imaging system capabilities. To accommodate nonstandard system geometries, a model-based iterative reconstruction method is applied. Such algorithms generally complicate the assessment and prediction of reconstructed image properties; however, we leverage efficient implementations of analytical predictors of local noise and spatial resolution that incorporate dependencies of the reconstruction algorithm on patient anatomy, x-ray technique, and geometry. These image property predictors serve as inputs to a task-based performance metric defined by detectability index, which is optimized with respect to the orbital parameters of data acquisition. We investigate the framework of the task-driven trajectory design in several examples to examine the dependence of optimal source-detector trajectories on the imaging task (or tasks), including location and spatial-frequency dependence. A variety of multitask objectives are also investigated, and the advantages to imaging performance are quantified in simulation studies.

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 p_l and p_h ; 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.

Estimating detectability index in vivo: development and validation of an automated methodology

This study's purpose was to develop and validate a method to estimate patient-specific detectability indices directly from patients' CT images (i.e., in vivo). The method extracts noise power spectrum (NPS) and modulation transfer function (MTF) resolution properties from each patient's CT series based on previously validated techniques. These are combined with a reference task function (10-mm disk lesion with [Formula: see text] HU contrast) to estimate detectability indices for a nonprewhitening matched filter observer model. This method was applied to CT data from a previous study in which diagnostic performance of 16 readers was measured for the task of detecting subtle, hypoattenuating liver lesions ([Formula: see text]), using a two-alternative-forced-choice (2AFC) method, over six dose levels and two reconstruction algorithms. In vivo detectability indices were estimated and compared to the human readers' binary 2AFC outcomes using a generalized linear mixed-effects statistical model. The results of this modeling showed that the in vivo detectability indices were strongly related to 2AFC outcomes ([Formula: see text]). Linear comparison between human-detection accuracy and model-predicted detection accuracy (for like conditions) resulted in Pearson and Spearman correlation coefficients exceeding 0.84. These results suggest the potential utility of using in vivo estimates of a detectability index for an automated image quality tracking system that could be implemented clinically.

Task-Driven Optimization of Fluence Field and Regularization for Model-Based Iterative Reconstruction in Computed Tomography

This paper presents a joint optimization of dynamic fluence field modulation (FFM) and regularization in quadratic penalized-likelihood reconstruction that maximizes a task-based imaging performance metric. We adopted a task-driven imaging framework for prospective designs of the imaging parameters. A maxi-min objective function was adopted to maximize the minimum detectability index ( ) throughout the image. The optimization algorithm alternates between FFM (represented by low-dimensional basis functions) and local regularization (including the regularization strength and directional penalty weights). The task-driven approach was compared with three FFM strategies commonly proposed for FBP reconstruction (as well as a task-driven TCM strategy) for a discrimination task in an abdomen phantom. The task-driven FFM assigned more fluence to less attenuating anteroposterior views and yielded approximately constant fluence behind the object. The optimal regularization was almost uniform throughout image. Furthermore, the task-driven FFM strategy redistribute fluence across detector elements in order to prescribe more fluence to the more attenuating central region of the phantom. Compared with all strategies, the task-driven FFM strategy not only improved minimum by at least 17.8%, but yielded higher over a large area inside the object. The optimal FFM was highly dependent on the amount of regularization, indicating the importance of a joint optimization. Sample reconstructions of simulated data generally support the performance estimates based on computed . The improvements in detectability show the potential of the task-driven imaging framework to improve imaging performance at a fixed dose, or, equivalently, to provide a similar level of performance at reduced dose.

Task-Based Design of Fluence Field Modulation in CT for Model-Based Iterative Reconstruction

A task-driven imaging framework for prospective fluence field modulation (FFM) is developed in this paper. The design approach uses a system model that includes a parameterized FFM acquisition and model-based iterative reconstruction (MBIR) for image formation. Using prior anatomical knowledge (e.g. from a low-dose 3D scout image), accurate predictions of spatial resolution and noise as a function of FFM are integrated into a task-based objective function. Specifically, detectability index (d'), a common metric for task-based image quality assessment, is computed for a specific formulation of the imaging task. To optimize imaging performance in across an image volume, a maximin objective function was adopted to maximize the minimum detectability index for many locations sampled throughout the volume. To reduce the dimensionality, FFM patterns were represented using wavelet bases, the coefficients of which were optimized using the covariance matrix adaptation evolutionary strategy (CMA-ES) algorithm. The optimization was performed for a mid-frequency discrimination task involving a cluster of micro-calcifications in an abdomen phantom. The task-driven design yielded FFM patterns that were significantly different from traditional strategies proposed for FBP reconstruction. In addition to a higher minimum d' consistent with the objective function, the task-driven approach also improved d' to a greater extent over a larger area of the phantom. Results from this work suggests that FFM strategies suitable for FBP reconstruction need to be reevaluated in the context of MBIR and that a task-driven imaging framework provides a promising approach for such optimization.

Task-driven image acquisition and reconstruction in cone-beam CT

This work introduces a task-driven imaging framework that incorporates a mathematical definition of the imaging task, a model of the imaging system, and a patient-specific anatomical model to prospectively design image acquisition and reconstruction techniques to optimize task performance. The framework is applied to joint optimization of tube current modulation, view-dependent reconstruction kernel, and orbital tilt in cone-beam CT. The system model considers a cone-beam CT system incorporating a flat-panel detector and 3D filtered backprojection and accurately describes the spatially varying noise and resolution over a wide range of imaging parameters in the presence of a realistic anatomical model. Task-based detectability index (d') is incorporated as the objective function in a task-driven optimization of image acquisition and reconstruction techniques. The orbital tilt was optimized through an exhaustive search across tilt angles ranging ± 30°. For each tilt angle, the view-dependent tube current and reconstruction kernel (i.e. the modulation profiles) that maximized detectability were identified via an alternating optimization. The task-driven approach was compared with conventional unmodulated and automatic exposure control (AEC) strategies for a variety of imaging tasks and anthropomorphic phantoms. The task-driven strategy outperformed the unmodulated and AEC cases for all tasks. For example, d' for a sphere detection task in a head phantom was improved by 30% compared to the unmodulated case by using smoother kernels for noisy views and distributing mAs across less noisy views (at fixed total mAs) in a manner that was beneficial to task performance. Similarly for detection of a line-pair pattern, the task-driven approach increased d' by 80% compared to no modulation by means of view-dependent mA and kernel selection that yields modulation transfer function and noise-power spectrum optimal to the task. Optimization of orbital tilt identified the tilt angle that reduced quantum noise in the region of the stimulus by avoiding highly attenuating anatomical structures. The task-driven imaging framework offers a potentially valuable paradigm for prospective definition of acquisition and reconstruction protocols that improve task performance without increase in dose.

Task-based detectability in CT image reconstruction by filtered backprojection and penalized likelihood estimation

PURPOSE

Nonstationarity is an important aspect of imaging performance in CT and cone-beam CT (CBCT), especially for systems employing iterative reconstruction. This work presents a theoretical framework for both filtered-backprojection (FBP) and penalized-likelihood (PL) reconstruction that includes explicit descriptions of nonstationary noise, spatial resolution, and task-based detectability index. Potential utility of the model was demonstrated in the optimal selection of regularization parameters in PL reconstruction.

METHODS

Analytical models for local modulation transfer function (MTF) and noise-power spectrum (NPS) were investigated for both FBP and PL reconstruction, including explicit dependence on the object and spatial location. For FBP, a cascaded systems analysis framework was adapted to account for nonstationarity by separately calculating fluence and system gains for each ray passing through any given voxel. For PL, the point-spread function and covariance were derived using the implicit function theorem and first-order Taylor expansion according to Fessler ["Mean and variance of implicitly defined biased estimators (such as penalized maximum likelihood): Applications to tomography," IEEE Trans. Image Process. 5(3), 493-506 (1996)]. Detectability index was calculated for a variety of simple tasks. The model for PL was used in selecting the regularization strength parameter to optimize task-based performance, with both a constant and a spatially varying regularization map.

RESULTS

Theoretical models of FBP and PL were validated in 2D simulated fan-beam data and found to yield accurate predictions of local MTF and NPS as a function of the object and the spatial location. The NPS for both FBP and PL exhibit similar anisotropic nature depending on the pathlength (and therefore, the object and spatial location within the object) traversed by each ray, with the PL NPS experiencing greater smoothing along directions with higher noise. The MTF of FBP is isotropic and independent of location to a first order approximation, whereas the MTF of PL is anisotropic in a manner complementary to the NPS. Task-based detectability demonstrates dependence on the task, object, spatial location, and smoothing parameters. A spatially varying regularization "map" designed from locally optimal regularization can improve overall detectability beyond that achievable with the commonly used constant regularization parameter.

CONCLUSIONS

Analytical models for task-based FBP and PL reconstruction are predictive of nonstationary noise and resolution characteristics, providing a valuable framework for understanding and optimizing system performance in CT and CBCT.

Modeling and Control of Nonstationary Noise Characteristics in Filtered-Backprojection and Penalized Likelihood Image Reconstruction

Purpose

Nonstationarity of CT noise presents a major challenge to the assessment of image quality. This work presents models for imaging performance in both filtered backprojection (FBP) and penalized likelihood (PL) reconstruction that describe not only the dependence on the imaging chain but also the dependence on the object as well as the nonstationary characteristics of the signal and noise. The work furthermore demonstrates the ability to impart control over the imaging process by adjusting reconstruction parameters to exploit nonstationarity in a manner advantageous to a particular imaging task.

Methods

A cascaded systems analysis model was used to model the local noise-power spectrum (NPS) and modulation transfer function (MTF) for FBP reconstruction, with locality achieved by separate calculation of fluence and system gain for each view as a function of detector location. The covariance and impulse response function for PL reconstruction (quadratic penalty) were computed using the implicit function theorem and Taylor expansion. Detectability index was calculated under the assumption of local stationarity to show the variation in task-dependent image quality throughout the image for simple and complex, heterogeneous objects. Control of noise magnitude and correlation was achieved by applying a spatially varying roughness penalty in PL reconstruction in a manner that improved overall detectability.

Results

The models provide a foundation for task-based imaging performance assessment in FBP and PL image reconstruction. For both FBP and PL, noise is anisotropic and varies in a manner dependent on the path length of each view traversing the object. The anisotropy in turn affects task performance, where detectability is enhanced or diminished depending on the frequency content of the task relative to that of the NPS. Spatial variation of the roughness penalty can be exploited to control noise magnitude and correlation (and hence detectability).

Conclusions

Nonstationarity of image noise is a significant effect that can be modeled in both FBP and PL image reconstruction. Prevalent spatial-frequency-dependent metrics of spatial resolution and noise can be analyzed under assumptions of local stationarity, providing a means to analyze imaging performance as a function of location throughout the image. Knowledgeable selection of a spatially-varying roughness penalty in PL can potentially improve local noise and spatial resolution in a manner tuned to a particular imaging task.

Theoretical Framework for the Dual-Energy Cone-Beam CT Noise-Power Spectrum, NEQ, and Tasked-Based Detectability Index

The optimization of dual-energy computed tomography (DE-CT) is challenged by the lack of a theoretical foundation for image quality. This work reports a cascaded systems analysis model that was used to derive signal and noise propagation in DE-CBCT in prevalent Fourier metrics such as the noise-power spectrum (NPS) and noise-equivalent quanta (NEQ). The model was validated in comparison to measurements of the 3D NPS and NEQ in DE-CBCT images acquired using an experimental imaging bench. Task-based detectability index was derived using DE-NPS and NEQ as an objective function in optimizing DE imaging parameters such as the dose allocation factor (D_A) and kVp pair. The resulting dose allocation optimization is in agreement with the practice of assigning more dose to the high-energy image (D_A < 0.5), and the model provides a quantitative basis for examining the optimal dose allocation as a function of total dose, kVp pair, the presence of electronics noise, and the imaging task. An example optimization is shown for a breast tumor detection task. Using DE decomposition to cancel fibroglandular tissue (rendering a DE-CBCT image of breast tumor against an adipose tissue background) and assuming a total dose of 15mGy, the optimal kVp pair is identified at [45, 105]kVp with D_A=0.46. The model is sufficiently general for applications beyond this example, demonstrating utility in the optimization in a broad range of imaging parameters. The model provides a new, valuable framework for understanding the theoretical limits of DE-CBCT imaging performance and maximizing image quality while minimizing radiation dose.

Task-based modeling and optimization of a cone-beam CT scanner for musculoskeletal imaging

PURPOSE

This work applies a cascaded systems model for cone-beam CT imaging performance to the design and optimization of a system for musculoskeletal extremity imaging. The model provides a quantitative guide to the selection of system geometry, source and detector components, acquisition techniques, and reconstruction parameters.

METHODS

The model is based on cascaded systems analysis of the 3D noise-power spectrum (NPS) and noise-equivalent quanta (NEQ) combined with factors of system geometry (magnification, focal spot size, and scatter-to-primary ratio) and anatomical background clutter. The model was extended to task-based analysis of detectability index (d') for tasks ranging in contrast and frequency content, and d' was computed as a function of system magnification, detector pixel size, focal spot size, kVp, dose, electronic noise, voxel size, and reconstruction filter to examine trade-offs and optima among such factors in multivariate analysis. The model was tested quantitatively versus the measured NPS and qualitatively in cadaver images as a function of kVp, dose, pixel size, and reconstruction filter under conditions corresponding to the proposed scanner.

RESULTS

The analysis quantified trade-offs among factors of spatial resolution, noise, and dose. System magnification (M) was a critical design parameter with strong effect on spatial resolution, dose, and x-ray scatter, and a fairly robust optimum was identified at M ∼ 1.3 for the imaging tasks considered. The results suggested kVp selection in the range of ∼65-90 kVp, the lower end (65 kVp) maximizing subject contrast and the upper end maximizing NEQ (90 kVp). The analysis quantified fairly intuitive results-e.g., ∼0.1-0.2 mm pixel size (and a sharp reconstruction filter) optimal for high-frequency tasks (bone detail) compared to ∼0.4 mm pixel size (and a smooth reconstruction filter) for low-frequency (soft-tissue) tasks. This result suggests a specific protocol for 1 × 1 (full-resolution) projection data acquisition followed by full-resolution reconstruction with a sharp filter for high-frequency tasks along with 2 × 2 binning reconstruction with a smooth filter for low-frequency tasks. The analysis guided selection of specific source and detector components implemented on the proposed scanner. The analysis also quantified the potential benefits and points of diminishing return in focal spot size, reduced electronic noise, finer detector pixels, and low-dose limits of detectability. Theoretical results agreed quantitatively with the measured NPS and qualitatively with evaluation of cadaver images by a musculoskeletal radiologist.

CONCLUSIONS

A fairly comprehensive model for 3D imaging performance in cone-beam CT combines factors of quantum noise, system geometry, anatomical background, and imaging task. The analysis provided a valuable, quantitative guide to design, optimization, and technique selection for a musculoskeletal extremities imaging system under development.

Anatomical background and generalized detectability in tomosynthesis and cone-beam CT

PURPOSE

Anatomical background presents a major impediment to detectability in 2D radiography as well as 3D tomosynthesis and cone-beam CT (CBCT). This article incorporates theoretical and experimental analysis of anatomical background "noise" in cascaded systems analysis of 2D and 3D imaging performance to yield "generalized" metrics of noise-equivalent quanta (NEQ) and detectability index as a function of the orbital extent of the (circular arc) source-detector orbit.

METHODS

A physical phantom was designed based on principles of fractal self-similarity to exhibit power-law spectral density (kappa/Fbeta) comparable to various anatomical sites (e.g., breast and lung). Background power spectra [S(B)(F)] were computed as a function of source-detector orbital extent, including tomosynthesis (approximately 10 degrees -180 degrees) and CBCT (180 degrees + fan to 360 degrees) under two acquisition schemes: (1) Constant angular separation between projections (variable dose) and (2) constant total number of projections (constant dose). The resulting S(B) was incorporated in the generalized NEQ, and detectability index was computed from 3D cascaded systems analysis for a variety of imaging tasks.

RESULTS

The phantom yielded power-law spectra within the expected spatial frequency range, quantifying the dependence of clutter magnitude (kappa) and correlation (beta) with increasing tomosynthesis angle. Incorporation of S(B) in the 3D NEQ provided a useful framework for analyzing the tradeoffs among anatomical, quantum, and electronic noise with dose and orbital extent. Distinct implications are posed for breast and chest tomosynthesis imaging system design-applications varying significantly in kappa and beta, and imaging task and, therefore, in optimal selection of orbital extent, number of projections, and dose. For example, low-frequency tasks (e.g., soft-tissue masses or nodules) tend to benefit from larger orbital extent and more fully 3D tomographic imaging, whereas high-frequency tasks (e.g., microcalcifications) require careful, application-specific selection of orbital extent and number of projections to minimize negative effects of quantum and electronic noise.

CONCLUSIONS

The complex tradeoffs among anatomical background, quantum noise, and electronic noise in projection imaging, tomosynthesis, and CBCT can be described by generalized cascaded systems analysis, providing a useful framework for system design and optimization.

The Generalized NEQ and Detectability Index for Tomosynthesis and Cone-Beam CT: From Cascaded Systems Analysis to Human Observers

PURPOSE

In the early development of new imaging modalities - such as tomosynthesis and cone-beam CT (CBCT) - an accurate predictive model for imaging performance is particularly valuable in identifying the physical factors that govern image quality and guiding system optimization. In this work, a task-based cascaded systems model for detectability index is proposed that describes not only the signal and noise propagation in the 2D (projection) and 3D (reconstruction) imaging chain but also the influence of background anatomical noise. The extent to which generalized detectability index provides a valid metric for imaging performance was assessed through direct comparison to human observer experiments.

METHODS

Detectability index (d') was generalized to include anatomical background noise in the same manner as the generalized noise-equivalent quanta (NEQ) proposed by Barrett et al. (Proc. SPIE Med. Imaging, Vol. 1090, 1989). Anatomical background noise was measured from a custom phantom designed to present power-law spectral density comparable to various anatomical sites (e.g., breast and lung). Theoretical calculations of d' as a function of the source-detector orbital extent (θtot) was obtained from a 3D cascaded systems analysis model for tomosynthesis and cone-beam CT (CBCT). Four model observers were considered in the calculation of d': prewhitening (PW), non-prewhitening (NPW), prewhitening with eye filter and internal noise (PWE), and non-prewhitening with eye filter and internal noise (NPWE). Human observer performance was measured from 9AFC tests for a variety of idealized imaging tasks presented within a clutter phantom. Theoretical results (d') were converted to area under the ROC curve (Az ) and compared directly to human observer performance as a function of imaging task and orbital extent.

RESULTS

Theoretical results demonstrated reasonable correspondence with human observer response for all tasks across the continuum in θtot ranging from low-angle tomosynthesis (θtot ~10°) to CBCT (θtot ~180°). Both theoretical and experimental Az were found to increase with acquisition angle, consistent with increased rejection of out-of-plane clutter for larger tomosynthesis angle. Of the four theoretical model observers considered, the prewhitening models tended to overestimate real observer performance, while the non-prewhitening models demonstrated reasonable agreement.

CONCLUSIONS

Generalized detectability index was shown to provide a meaningful metric for imaging performance, helping to bridge the gap between real observer performance and prevalent Fourier-based metrics based in first principles of spatial-frequency-dependent NEQ and imaging task.

Detectability index measures of binaural masking level difference across populations of inferior colliculus neurons

In everyday life we continually need to detect signals against a background of interfering noise (the "cocktail party effect"): a task that is much easier to accomplish using two ears. The binaural masking level difference (BMLD) measures the ability of listeners to use a difference in binaural attributes to segregate sound sources and thus improve their discriminability against interfering noises. By computing the detectability of tones from rate-versus-level functions in the presence of a suprathreshold noise, we previously demonstrated that individual low-frequency delay-sensitive neurons in the inferior colliculus are able to show BMLDs. Here we consider the responses of a population of such neurons when the noise level is held constant (as conventionally in psychophysical paradigms). We have sampled the responses of 121 units in the inferior colliculi of five guinea pigs to identical noise and 500 Hz tones at both ears (NoSo) and to identical noise but with the 500 Hz tone at one ear inverted (NoSpi). The result suggests that the neurons subserving detection of So tones in No (identical noise at the two ears) noise are those neurons with best frequencies (BFs) close to 500 Hz that respond to So tones with an increase in their discharge rate from that attributable to the noise. The detection of the inverted (Spi) signal is also attributable to neurons with BFs close to 500 Hz. However, among these neurons, the presence of the Spi tone was indicated by an increased discharge rate in some neurons and by a decreased discharge rate in others.