1
|
Zhou W, Villa U, Anastasio MA. Ideal Observer Computation by Use of Markov-Chain Monte Carlo With Generative Adversarial Networks. IEEE TRANSACTIONS ON MEDICAL IMAGING 2023; 42:3715-3724. [PMID: 37578916 PMCID: PMC10769588 DOI: 10.1109/tmi.2023.3304907] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 08/16/2023]
Abstract
Medical imaging systems are often evaluated and optimized via objective, or task-specific, measures of image quality (IQ) that quantify the performance of an observer on a specific clinically-relevant task. The performance of the Bayesian Ideal Observer (IO) sets an upper limit among all observers, numerical or human, and has been advocated for use as a figure-of-merit (FOM) for evaluating and optimizing medical imaging systems. However, the IO test statistic corresponds to the likelihood ratio that is intractable to compute in the majority of cases. A sampling-based method that employs Markov-chain Monte Carlo (MCMC) techniques was previously proposed to estimate the IO performance. However, current applications of MCMC methods for IO approximation have been limited to a small number of situations where the considered distribution of to-be-imaged objects can be described by a relatively simple stochastic object model (SOM). As such, there remains an important need to extend the domain of applicability of MCMC methods to address a large variety of scenarios where IO-based assessments are needed but the associated SOMs have not been available. In this study, a novel MCMC method that employs a generative adversarial network (GAN)-based SOM, referred to as MCMC-GAN, is described and evaluated. The MCMC-GAN method was quantitatively validated by use of test-cases for which reference solutions were available. The results demonstrate that the MCMC-GAN method can extend the domain of applicability of MCMC methods for conducting IO analyses of medical imaging systems.
Collapse
|
2
|
Granstedt JL, Zhou W, Anastasio MA. Approximating the Hotelling observer with autoencoder-learned efficient channels for binary signal detection tasks. J Med Imaging (Bellingham) 2023; 10:055501. [PMID: 37767114 PMCID: PMC10520791 DOI: 10.1117/1.jmi.10.5.055501] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2023] [Revised: 08/31/2023] [Accepted: 09/05/2023] [Indexed: 09/29/2023] Open
Abstract
Purpose The objective assessment of image quality (IQ) has been advocated for the analysis and optimization of medical imaging systems. One method of computing such IQ metrics is through a numerical observer. The Hotelling observer (HO) is the optimal linear observer, but conventional methods for obtaining the HO can become intractable due to large image sizes or insufficient data. Channelized methods are sometimes employed in such circumstances to approximate the HO. The performance of such channelized methods varies, with different methods obtaining superior performance to others depending on the imaging conditions and detection task. A channelized HO method using an AE is presented and implemented across several tasks to characterize its performance. Approach The process for training an AE is demonstrated to be equivalent to developing a set of channels for approximating the HO. The efficiency of the learned AE-channels is increased by modifying the conventional AE loss function to incorporate task-relevant information. Multiple binary detection tasks involving lumpy and breast phantom backgrounds across varying dataset sizes are considered to evaluate the performance of the proposed method and compare to current state-of-the-art channelized methods. Additionally, the ability of the channelized methods to generalize to images outside of the training dataset is investigated. Results AE-learned channels are demonstrated to have comparable performance with other state-of-the-art channel methods in the detection studies and superior performance in the generalization studies. Incorporating a cleaner estimate of the signal for the detection task is also demonstrated to significantly improve the performance of the proposed method, particularly in datasets with fewer images. Conclusions AEs are demonstrated to be capable of learning efficient channels for the HO. The resulting significant increase in detection performance for small dataset sizes when incorporating a signal prior holds promising implications for future assessments of imaging technologies.
Collapse
Affiliation(s)
- Jason L. Granstedt
- University of Illinois Urbana-Champaign, Department of Computer Science, Champaign, Illinois, United States
| | - Weimin Zhou
- Shanghai Jiao Tong University, Global Institute of Future Technology, Shanghai, China
| | - Mark A. Anastasio
- University of Illinois Urbana-Champaign, Department of Computer Science, Champaign, Illinois, United States
- University of Illinois Urbana-Champaign, Department of Bioengineering, Champaign, Illinois, United States
| |
Collapse
|
3
|
Zhang X, Kelkar VA, Granstedt J, Li H, Anastasio MA. Impact of deep learning-based image super-resolution on binary signal detection. J Med Imaging (Bellingham) 2021; 8:065501. [PMID: 34796251 PMCID: PMC8594450 DOI: 10.1117/1.jmi.8.6.065501] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2021] [Accepted: 10/27/2021] [Indexed: 11/14/2022] Open
Abstract
Purpose: Deep learning-based image super-resolution (DL-SR) has shown great promise in medical imaging applications. To date, most of the proposed methods for DL-SR have only been assessed using traditional measures of image quality (IQ) that are commonly employed in the field of computer vision. However, the impact of these methods on objective measures of IQ that are relevant to medical imaging tasks remains largely unexplored. We investigate the impact of DL-SR methods on binary signal detection performance. Approach: Two popular DL-SR methods, the super-resolution convolutional neural network and the super-resolution generative adversarial network, were trained using simulated medical image data. Binary signal-known-exactly with background-known-statistically and signal-known-statistically with background-known-statistically detection tasks were formulated. Numerical observers (NOs), which included a neural network-approximated ideal observer and common linear NOs, were employed to assess the impact of DL-SR on task performance. The impact of the complexity of the DL-SR network architectures on task performance was quantified. In addition, the utility of DL-SR for improving the task performance of suboptimal observers was investigated. Results: Our numerical experiments confirmed that, as expected, DL-SR improved traditional measures of IQ. However, for many of the study designs considered, the DL-SR methods provided little or no improvement in task performance and even degraded it. It was observed that DL-SR improved the task performance of suboptimal observers under certain conditions. Conclusions: Our study highlights the urgent need for the objective assessment of DL-SR methods and suggests avenues for improving their efficacy in medical imaging applications.
Collapse
Affiliation(s)
- Xiaohui Zhang
- University of Illinois at Urbana–Champaign, Department of Bioengineering, Urbana, Illinois, United States
| | - Varun A. Kelkar
- University of Illinois at Urbana–Champaign, Department of Electrical and Computer Engineering, Urbana, Illinois, United States
| | - Jason Granstedt
- University of Illinois at Urbana–Champaign, Department of Computer Science, Urbana, Illinois, United States
| | - Hua Li
- University of Illinois at Urbana–Champaign, Department of Bioengineering, Urbana, Illinois, United States
- University of Illinois at Urbana–Champaign, Cancer Center at Illinois, Urbana, Illinois, United States
- Carle Foundation Hospital, Carle Cancer Center, Urbana, Illinois, United States
| | - Mark A. Anastasio
- University of Illinois at Urbana–Champaign, Department of Bioengineering, Urbana, Illinois, United States
- University of Illinois at Urbana–Champaign, Department of Electrical and Computer Engineering, Urbana, Illinois, United States
- University of Illinois at Urbana–Champaign, Department of Computer Science, Urbana, Illinois, United States
| |
Collapse
|
4
|
Hernandez AM, Becker AE, Hyun Lyu S, Abbey CK, Boone JM. High-resolution μ CT imaging for characterizing microcalcification detection performance in breast CT. J Med Imaging (Bellingham) 2021; 8:052107. [PMID: 34307737 PMCID: PMC8291078 DOI: 10.1117/1.jmi.8.5.052107] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2021] [Accepted: 06/28/2021] [Indexed: 01/07/2023] Open
Abstract
Purpose: To demonstrate the utility of high-resolution micro-computed tomography ( μ CT ) for determining ground-truth size and shape properties of calcium grains for evaluation of detection performance in breast CT (bCT). Approach: Calcium carbonate grains ( ∼ 200 μ m ) were suspended in 1% agar solution to emulate microcalcifications ( μ Calcs ) within a fibroglandular tissue background. Ground-truth imaging was performed on a commercial μ CT scanner and was used for assessing calcium-grain size and shape, and for generating μ Calc signal profiles. Calcium grains were placed within a realistic breast-shaped phantom and imaged on a prototype bCT system at 3- and 6-mGy mean glandular dose (MGD) levels, and the non-prewhitening detectability was assessed. Additionally, the μ CT -derived signal profiles were used in conjunction with the bCT system characterization (MTF and NPS) to obtain predictions of bCT detectability. Results: Estimated detectability of the calcium grains on the bCT system ranged from 2.5 to 10.6 for 3 mGy and from 3.8 to 15.3 for 6 mGy with large fractions of the grains meeting the Rose criterion for visibility. Segmentation of μ CT images based on morphological operations produced accurate results in terms of segmentation boundaries and segmented region size. A regression model linking bCT detectability to μ Calc parameters indicated significant effects of μ Calc size and vertical position within the breast phantom. Detectability using μ CT -derived detection templates and bCT statistical properties (MTF and NPS) were in good correspondence with those measured directly from bCT ( R 2 > 0.88 ). Conclusions: Parameters derived from μ CT ground-truth data were shown to produce useful characterizations of detectability when compared to estimates derived directly from bCT. Signal profiles derived from μ CT imaging can be used in conjunction with measured or hypothesized statistical properties to evaluate the performance of a system, or system component, that may not currently be available.
Collapse
Affiliation(s)
- Andrew M. Hernandez
- University of California Davis, Department of Radiology, Sacramento, California, United States,Address all correspondence to Andrew M. Hernandez,
| | - Amy E. Becker
- University of California Davis, Biomedical Engineering Graduate Group, Davis, California, United States
| | - Su Hyun Lyu
- University of California Davis, Biomedical Engineering Graduate Group, Davis, California, United States
| | - Craig K. Abbey
- University of California Santa Barbara, Psychological and Brain Sciences, Santa Barbara, California, United States
| | - John M. Boone
- University of California Davis, Department of Radiology, Sacramento, California, United States,University of California Davis, Biomedical Engineering, Davis, California, United States
| |
Collapse
|
5
|
Sidky EY, Phillips JP, Zhou W, Ongie G, Cruz-Bastida JP, Reiser IS, Anastasio MA, Pan X. A signal detection model for quantifying overregularization in nonlinear image reconstruction. Med Phys 2021; 48:6312-6323. [PMID: 34169538 DOI: 10.1002/mp.14703] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2020] [Revised: 12/09/2020] [Accepted: 12/21/2020] [Indexed: 11/08/2022] Open
Abstract
Many useful image quality metrics for evaluating linear image reconstruction techniques do not apply to or are difficult to interpret for nonlinear image reconstruction. The vast majority of metrics employed for evaluating nonlinear image reconstruction are based on some form of global image fidelity, such as image root mean square error (RMSE). Use of such metrics can lead to overregularization in the sense that they can favor removal of subtle details in the image. To address this shortcoming, we develop an image quality metric based on signal detection that serves as a surrogate to the qualitative loss of fine image details. The metric is demonstrated in the context of a breast CT simulation, where different equal-dose configurations are considered. The configurations differ in the number of projections acquired. Image reconstruction is performed with a nonlinear algorithm based on total variation constrained least-squares (TV-LSQ). The resulting images are studied as a function of three parameters: number of views acquired, total variation constraint value, and number of iterations. The images are evaluated visually, with image RMSE, and with the proposed signal-detection-based metric. The latter uses a small signal, and computes detectability in the sinogram and in the reconstructed image. Loss of signal detectability through the image reconstruction process is taken as a quantitative measure of loss of fine details in the image. Loss of signal detectability is seen to correlate well with the blocky or patchy appearance due to overregularization with TV-LSQ, and this trend runs counter to the image RMSE metric, which tends to favor the over-regularized images. The proposed signal detection-based metric provides an image quality assessment that is complimentary to that of image RMSE. Using the two metrics in concert may yield a useful prescription for determining CT algorithm and configuration parameters when nonlinear image reconstruction is used.
Collapse
Affiliation(s)
- Emil Y Sidky
- Department of Radiology, The University of Chicago, 5841 S. Maryland Ave., Chicago, IL, 60637, USA
| | - John Paul Phillips
- Department of Radiology, The University of Chicago, 5841 S. Maryland Ave., Chicago, IL, 60637, USA
| | - Weimin Zhou
- Department of Bioengineering, University of Illinois at Urbana-Champaign, 1406 W. Green St., Urbana, IL, 61801, USA
| | - Greg Ongie
- Department of Mathematical and Statistical Sciences, Marquette University, 1313 W. Wisconsin Ave., Milwaukee, WI, 53233, USA
| | - Juan P Cruz-Bastida
- Department of Radiology, The University of Chicago, 5841 S. Maryland Ave., Chicago, IL, 60637, USA
| | - Ingrid S Reiser
- Department of Radiology, The University of Chicago, 5841 S. Maryland Ave., Chicago, IL, 60637, USA
| | - Mark A Anastasio
- Department of Bioengineering, University of Illinois at Urbana-Champaign, 1406 W. Green St., Urbana, IL, 61801, USA
| | - Xiaochuan Pan
- Department of Radiology, The University of Chicago, 5841 S. Maryland Ave., Chicago, IL, 60637, USA
| |
Collapse
|
6
|
Chen L, Chen J. Deep Neural Network for Automatic Classification of Pathological Voice Signals. J Voice 2020; 36:288.e15-288.e24. [PMID: 32660846 DOI: 10.1016/j.jvoice.2020.05.029] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2020] [Revised: 05/17/2020] [Accepted: 05/26/2020] [Indexed: 10/23/2022]
Abstract
OBJECTIVES Computer-aided pathological voice detection is efficient for initial screening of pathological voice, and has received high academic and clinical attention. This paper proposes an automatic diagnosis method of pathological voice based on deep neural network (DNN). Other two classification models (support vector machines and random forests) were used to verify the effectiveness of DNN. METHODS In this paper, we extracted 12 Mel frequency cepstral coefficients of each voice sample as row features. The constructed DNN consists a two-layer stacked sparse autoencoders network and a softmax layer. The stacked sparse autoencoders layer can learn high-level features from raw Mel frequency cepstral coefficients features. Then, the softmax layer can diagnose pathological voice according to high-level features. The DNN and the other two comparison models used the same train set and test set for the experiment. RESULTS Experimental results reveal that the value of sensitivity, specificity, precision, accuracy, and F1 score of the DNN can reach 97.8%, 99.4%, 99.4%, 98.6%, and 98.4%, respectively. The five indexes of DNN classification results are at least 6.2%, 5%, 5.6%, 5.7%, and 6.2% higher than the comparison models (support vector machine and random forest). CONCLUSIONS The proposed DNN can learn advanced features from raw acoustic features, and distinguish pathological voice from healthy voice. To the extent of this preliminary study, future studies can further explore the application of DNN in other experiments and clinical practice.
Collapse
Affiliation(s)
- Lili Chen
- School of Mechatronics and Vehicle Engineering, Chongqing Jiaotong University, Chongqing, China; Chongqing Survey Institute, Chongqing, China.
| | - Junjiang Chen
- School of Mechatronics and Vehicle Engineering, Chongqing Jiaotong University, Chongqing, China
| |
Collapse
|
7
|
Zhou W, Li H, Anastasio MA. Approximating the Ideal Observer and Hotelling Observer for Binary Signal Detection Tasks by Use of Supervised Learning Methods. IEEE TRANSACTIONS ON MEDICAL IMAGING 2019; 38:2456-2468. [PMID: 30990425 PMCID: PMC6858982 DOI: 10.1109/tmi.2019.2911211] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/26/2023]
Abstract
It is widely accepted that the optimization of medical imaging system performance should be guided by task-based measures of image quality (IQ). Task-based measures of IQ quantify the ability of an observer to perform a specific task, such as detection or estimation of a signal (e.g., a tumor). For binary signal detection tasks, the Bayesian Ideal Observer (IO) sets an upper limit of observer performance and has been advocated for use in optimizing medical imaging systems and data-acquisition designs. Except in special cases, the determination of the IO test statistic is analytically intractable. Markov-chain Monte Carlo (MCMC) techniques can be employed to approximate the IO detection performance, but their reported applications have been limited to relatively simple object models. In cases where the IO test statistic is difficult to compute, the Hotelling Observer (HO) can be employed. To compute the HO test statistic, potentially large covariance matrices must be accurately estimated and subsequently inverted, which can present computational challenges. This paper investigates the supervised learning-based methodologies for approximating the IO and HO test statistics. Convolutional neural networks (CNNs) and single-layer neural networks (SLNNs) are employed to approximate the IO and HO test statistics, respectively. The numerical simulations were conducted for both signal-known-exactly (SKE) and signal-known-statistically (SKS) signal detection tasks. The considered background models include the lumpy object model and the clustered lumpy object model. The measurement noise models considered are Gaussian, Laplacian, and mixed Poisson-Gaussian. The performances of the supervised learning methods are assessed via receiver operating characteristic (ROC) analysis, and the results are compared to those produced by the use of traditional numerical methods or analytical calculations when feasible. The potential advantages of the proposed supervised learning approaches for approximating the IO and HO test statistics are discussed.
Collapse
|
8
|
Chen Y, Lou Y, Wang K, Kupinski MA, Anastasio MA. Reconstruction-Aware Imaging System Ranking by Use of a Sparsity-Driven Numerical Observer Enabled by Variational Bayesian Inference. IEEE TRANSACTIONS ON MEDICAL IMAGING 2019; 38:1251-1262. [PMID: 30475713 PMCID: PMC6559219 DOI: 10.1109/tmi.2018.2880870] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
It is widely accepted that optimization of imaging system performance should be guided by task-based measures of image quality. It has been advocated that imaging hardware or data-acquisition designs should be optimized by use of an ideal observer that exploits full statistical knowledge of the measurement noise and class of objects to be imaged, without consideration of the reconstruction method. In practice, accurate and tractable models of the complete object statistics are often difficult to determine. Moreover, in imaging systems that employ compressive sensing concepts, imaging hardware and sparse image reconstruction are innately coupled technologies. In this paper, a sparsity-driven observer (SDO) that can be employed to optimize hardware by use of a stochastic object model describing object sparsity is described and investigated. The SDO and sparse reconstruction method can, therefore, be "matched" in the sense that they both utilize the same statistical information regarding the class of objects to be imaged. To efficiently compute the SDO test statistic, computational tools developed recently for variational Bayesian inference with sparse linear models are adopted. The use of the SDO to rank data-acquisition designs in a stylized example as motivated by magnetic resonance imaging is demonstrated. This paper reveals that the SDO can produce rankings that are consistent with visual assessments of the reconstructed images but different from those produced by use of the traditionally employed Hotelling observer.
Collapse
|
9
|
Rose SD, Sanchez AA, Sidky EY, Pan X. Investigating simulation-based metrics for characterizing linear iterative reconstruction in digital breast tomosynthesis. Med Phys 2018; 44:e279-e296. [PMID: 28901614 DOI: 10.1002/mp.12445] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2016] [Revised: 05/29/2017] [Accepted: 06/21/2017] [Indexed: 11/08/2022] Open
Abstract
PURPOSE Simulation-based image quality metrics are adapted and investigated for characterizing the parameter dependences of linear iterative image reconstruction for DBT. METHODS Three metrics based on a 2D DBT simulation are investigated: (1) a root-mean-square-error (RMSE) between the test phantom and reconstructed image, (2) a gradient RMSE where the comparison is made after taking a spatial gradient of both image and phantom, and (3) a region-of-interest (ROI) Hotelling observer (HO) for signal-known-exactly/background-known-exactly (SKE/BKE) and signal-known-exactly/background-known-statistically (SKE/BKS) detection tasks. Two simulation studies are performed using the aforementioned metrics, varying voxel aspect ratio, and regularization strength for two types of Tikhonov-regularized least-squares optimization. The RMSE metrics are applied to a 2D test phantom with resolution bar patterns at varying angles, and the ROI-HO metric is applied to two tasks relevant to DBT: lesion detection, modeled by use of a large, low-contrast signal, and microcalcification detection, modeled by use of a small, high-contrast signal. The RMSE metric trends are compared with visual assessment of the reconstructed bar-pattern phantom. The ROI-HO metric trends are compared with 3D reconstructed images from ACR phantom data acquired with a Hologic Selenia Dimensions DBT system. RESULTS Sensitivity of the image RMSE to mean pixel value is found to limit its applicability to the assessment of DBT image reconstruction. The image gradient RMSE is insensitive to mean pixel value and appears to track better with subjective visualization of the reconstructed bar-pattern phantom. The ROI-HO metric shows an increasing trend with regularization strength for both forms of Tikhonov-regularized least-squares; however, this metric saturates at intermediate regularization strength indicating a point of diminishing returns for signal detection. Visualization with the reconstructed ACR phantom images appear to show a similar dependence with regularization strength. CONCLUSIONS From the limited studies presented it appears that image gradient RMSE trends correspond with visual assessment better than image RMSE for DBT image reconstruction. The ROI-HO metric for both detection tasks also appears to reflect visual trends in the ACR phantom reconstructions as a function of regularization strength. We point out, however, that the true utility of these metrics can only be assessed after amassing more data.
Collapse
Affiliation(s)
- Sean D Rose
- University of Chicago, Department of Radiology MC-2026, 5841 S. Maryland Avenue, Chicago, IL, 60637, USA
| | - Adrian A Sanchez
- University of Chicago, Department of Radiology MC-2026, 5841 S. Maryland Avenue, Chicago, IL, 60637, USA
| | - Emil Y Sidky
- University of Chicago, Department of Radiology MC-2026, 5841 S. Maryland Avenue, Chicago, IL, 60637, USA
| | - Xiaochuan Pan
- University of Chicago, Department of Radiology MC-2026, 5841 S. Maryland Avenue, Chicago, IL, 60637, USA
| |
Collapse
|
10
|
Gang GJ, Siewerdsen JH, Webster Stayman J. Task-driven optimization of CT tube current modulation and regularization in model-based iterative reconstruction. Phys Med Biol 2017. [PMID: 28362638 DOI: 10.1088/1361-6560/aa6a97/meta] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/30/2023]
Abstract
Tube current modulation (TCM) is routinely adopted on diagnostic CT scanners for dose reduction. Conventional TCM strategies are generally designed for filtered-backprojection (FBP) reconstruction to satisfy simple image quality requirements based on noise. This work investigates TCM designs for model-based iterative reconstruction (MBIR) to achieve optimal imaging performance as determined by a task-based image quality metric. Additionally, regularization is an important aspect of MBIR that is jointly optimized with TCM, and includes both the regularization strength that controls overall smoothness as well as directional weights that permits control of the isotropy/anisotropy of the local noise and resolution properties. Initial investigations focus on a known imaging task at a single location in the image volume. The framework adopts Fourier and analytical approximations for fast estimation of the local noise power spectrum (NPS) and modulation transfer function (MTF)-each carrying dependencies on TCM and regularization. For the single location optimization, the local detectability index (d') of the specific task was directly adopted as the objective function. A covariance matrix adaptation evolution strategy (CMA-ES) algorithm was employed to identify the optimal combination of imaging parameters. Evaluations of both conventional and task-driven approaches were performed in an abdomen phantom for a mid-frequency discrimination task in the kidney. Among the conventional strategies, the TCM pattern optimal for FBP using a minimum variance criterion yielded a worse task-based performance compared to an unmodulated strategy when applied to MBIR. Moreover, task-driven TCM designs for MBIR were found to have the opposite behavior from conventional designs for FBP, with greater fluence assigned to the less attenuating views of the abdomen and less fluence to the more attenuating lateral views. Such TCM patterns exaggerate the intrinsic anisotropy of the MTF and NPS as a result of the data weighting in MBIR. Directional penalty design was found to reinforce the same trend. The task-driven approaches outperform conventional approaches, with the maximum improvement in d' of 13% given by the joint optimization of TCM and regularization. This work demonstrates that the TCM optimal for MBIR is distinct from conventional strategies proposed for FBP reconstruction and strategies optimal for FBP are suboptimal and may even reduce performance when applied to MBIR. The task-driven imaging framework offers a promising approach for optimizing acquisition and reconstruction for MBIR that can improve imaging performance and/or dose utilization beyond conventional imaging strategies.
Collapse
Affiliation(s)
- Grace J Gang
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD 21205, United States of America
| | | | | |
Collapse
|
11
|
Gang GJ, Siewerdsen JH, Webster Stayman J. Task-driven optimization of CT tube current modulation and regularization in model-based iterative reconstruction. Phys Med Biol 2017; 62:4777-4797. [PMID: 28362638 DOI: 10.1088/1361-6560/aa6a97] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
Tube current modulation (TCM) is routinely adopted on diagnostic CT scanners for dose reduction. Conventional TCM strategies are generally designed for filtered-backprojection (FBP) reconstruction to satisfy simple image quality requirements based on noise. This work investigates TCM designs for model-based iterative reconstruction (MBIR) to achieve optimal imaging performance as determined by a task-based image quality metric. Additionally, regularization is an important aspect of MBIR that is jointly optimized with TCM, and includes both the regularization strength that controls overall smoothness as well as directional weights that permits control of the isotropy/anisotropy of the local noise and resolution properties. Initial investigations focus on a known imaging task at a single location in the image volume. The framework adopts Fourier and analytical approximations for fast estimation of the local noise power spectrum (NPS) and modulation transfer function (MTF)-each carrying dependencies on TCM and regularization. For the single location optimization, the local detectability index (d') of the specific task was directly adopted as the objective function. A covariance matrix adaptation evolution strategy (CMA-ES) algorithm was employed to identify the optimal combination of imaging parameters. Evaluations of both conventional and task-driven approaches were performed in an abdomen phantom for a mid-frequency discrimination task in the kidney. Among the conventional strategies, the TCM pattern optimal for FBP using a minimum variance criterion yielded a worse task-based performance compared to an unmodulated strategy when applied to MBIR. Moreover, task-driven TCM designs for MBIR were found to have the opposite behavior from conventional designs for FBP, with greater fluence assigned to the less attenuating views of the abdomen and less fluence to the more attenuating lateral views. Such TCM patterns exaggerate the intrinsic anisotropy of the MTF and NPS as a result of the data weighting in MBIR. Directional penalty design was found to reinforce the same trend. The task-driven approaches outperform conventional approaches, with the maximum improvement in d' of 13% given by the joint optimization of TCM and regularization. This work demonstrates that the TCM optimal for MBIR is distinct from conventional strategies proposed for FBP reconstruction and strategies optimal for FBP are suboptimal and may even reduce performance when applied to MBIR. The task-driven imaging framework offers a promising approach for optimizing acquisition and reconstruction for MBIR that can improve imaging performance and/or dose utilization beyond conventional imaging strategies.
Collapse
Affiliation(s)
- Grace J Gang
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD 21205, United States of America
| | | | | |
Collapse
|
12
|
Shi L, Vedantham S, Karellas A, Zhu L. Library based x-ray scatter correction for dedicated cone beam breast CT. Med Phys 2016; 43:4529. [PMID: 27487870 PMCID: PMC4947049 DOI: 10.1118/1.4955121] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/20/2016] [Revised: 06/15/2016] [Accepted: 06/21/2016] [Indexed: 01/02/2023] Open
Abstract
PURPOSE The image quality of dedicated cone beam breast CT (CBBCT) is limited by substantial scatter contamination, resulting in cupping artifacts and contrast-loss in reconstructed images. Such effects obscure the visibility of soft-tissue lesions and calcifications, which hinders breast cancer detection and diagnosis. In this work, we propose a library-based software approach to suppress scatter on CBBCT images with high efficiency, accuracy, and reliability. METHODS The authors precompute a scatter library on simplified breast models with different sizes using the geant4-based Monte Carlo (MC) toolkit. The breast is approximated as a semiellipsoid with homogeneous glandular/adipose tissue mixture. For scatter correction on real clinical data, the authors estimate the breast size from a first-pass breast CT reconstruction and then select the corresponding scatter distribution from the library. The selected scatter distribution from simplified breast models is spatially translated to match the projection data from the clinical scan and is subtracted from the measured projection for effective scatter correction. The method performance was evaluated using 15 sets of patient data, with a wide range of breast sizes representing about 95% of general population. Spatial nonuniformity (SNU) and contrast to signal deviation ratio (CDR) were used as metrics for evaluation. RESULTS Since the time-consuming MC simulation for library generation is precomputed, the authors' method efficiently corrects for scatter with minimal processing time. Furthermore, the authors find that a scatter library on a simple breast model with only one input parameter, i.e., the breast diameter, sufficiently guarantees improvements in SNU and CDR. For the 15 clinical datasets, the authors' method reduces the average SNU from 7.14% to 2.47% in coronal views and from 10.14% to 3.02% in sagittal views. On average, the CDR is improved by a factor of 1.49 in coronal views and 2.12 in sagittal views. CONCLUSIONS The library-based scatter correction does not require increase in radiation dose or hardware modifications, and it improves over the existing methods on implementation simplicity and computational efficiency. As demonstrated through patient studies, the authors' approach is effective and stable, and is therefore clinically attractive for CBBCT imaging.
Collapse
Affiliation(s)
- Linxi Shi
- Nuclear and Radiological Engineering and Medical Physics Programs, The George W. Woodruff School of Mechanical Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332
| | - Srinivasan Vedantham
- Department of Radiology, University of Massachusetts Medical School, Worcester, Massachusetts 01655
| | - Andrew Karellas
- Department of Radiology, University of Massachusetts Medical School, Worcester, Massachusetts 01655
| | - Lei Zhu
- Nuclear and Radiological Engineering and Medical Physics Programs, The George W. Woodruff School of Mechanical Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332
| |
Collapse
|
13
|
Lorsakul A, Fakhri GE, Worstell W, Ouyang J, Rakvongthai Y, Laine AF, Li Q. Numerical observer for atherosclerotic plaque classification in spectral computed tomography. J Med Imaging (Bellingham) 2016; 3:035501. [PMID: 27429999 PMCID: PMC4940624 DOI: 10.1117/1.jmi.3.3.035501] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2016] [Accepted: 06/20/2016] [Indexed: 11/14/2022] Open
Abstract
Spectral computed tomography (SCT) generates better image quality than conventional computed tomography (CT). It has overcome several limitations for imaging atherosclerotic plaque. However, the literature evaluating the performance of SCT based on objective image assessment is very limited for the task of discriminating plaques. We developed a numerical-observer method and used it to assess performance on discrimination vulnerable-plaque features and compared the performance among multienergy CT (MECT), dual-energy CT (DECT), and conventional CT methods. Our numerical observer was designed to incorporate all spectral information and comprised two-processing stages. First, each energy-window domain was preprocessed by a set of localized channelized Hotelling observers (CHO). In this step, the spectral image in each energy bin was decorrelated using localized prewhitening and matched filtering with a set of Laguerre-Gaussian channel functions. Second, the series of the intermediate scores computed from all the CHOs were integrated by a Hotelling observer with an additional prewhitening and matched filter. The overall signal-to-noise ratio (SNR) and the area under the receiver operating characteristic curve (AUC) were obtained, yielding an overall discrimination performance metric. The performance of our new observer was evaluated for the particular binary classification task of differentiating between alternative plaque characterizations in carotid arteries. A clinically realistic model of signal variability was also included in our simulation of the discrimination tasks. The inclusion of signal variation is a key to applying the proposed observer method to spectral CT data. Hence, the task-based approaches based on the signal-known-exactly/background-known-exactly (SKE/BKE) framework and the clinical-relevant signal-known-statistically/background-known-exactly (SKS/BKE) framework were applied for analytical computation of figures of merit (FOM). Simulated data of a carotid-atherosclerosis patient were used to validate our methods. We used an extended cardiac-torso anthropomorphic digital phantom and three simulated plaque types (i.e., calcified plaque, fatty-mixed plaque, and iodine-mixed blood). The images were reconstructed using a standard filtered backprojection (FBP) algorithm for all the acquisition methods and were applied to perform two different discrimination tasks of: (1) calcified plaque versus fatty-mixed plaque and (2) calcified plaque versus iodine-mixed blood. MECT outperformed DECT and conventional CT systems for all cases of the SKE/BKE and SKS/BKE tasks (all [Formula: see text]). On average of signal variability, MECT yielded the SNR improvements over other acquisition methods in the range of 46.8% to 65.3% (all [Formula: see text]) for FBP-Ramp images and 53.2% to 67.7% (all [Formula: see text]) for FBP-Hanning images for both identification tasks. This proposed numerical observer combined with our signal variability framework is promising for assessing material characterization obtained through the additional energy-dependent attenuation information of SCT. These methods can be further extended to other clinical tasks such as kidney or urinary stone identification applications.
Collapse
Affiliation(s)
- Auranuch Lorsakul
- Massachusetts General Hospital, Division of Nuclear Medicine and Molecular Imaging, Gordon Center for Medical Imaging, 55 Fruit Street, White 427, Boston, Massachusetts 02114, United States
- Columbia University, Department of Biomedical Engineering, 1210 Amsterdam Avenue, New York, New York 10027, United States
| | - Georges El Fakhri
- Massachusetts General Hospital, Division of Nuclear Medicine and Molecular Imaging, Gordon Center for Medical Imaging, 55 Fruit Street, White 427, Boston, Massachusetts 02114, United States
- Harvard Medical School, Department of Radiology, 25 Shattuck Street, Boston, Massachusetts 02115, United States
| | - William Worstell
- PhotoDiagnostic System Inc., 85 Swanson Road, Boxborough, Massachusetts 01719, United States
| | - Jinsong Ouyang
- Massachusetts General Hospital, Division of Nuclear Medicine and Molecular Imaging, Gordon Center for Medical Imaging, 55 Fruit Street, White 427, Boston, Massachusetts 02114, United States
- Harvard Medical School, Department of Radiology, 25 Shattuck Street, Boston, Massachusetts 02115, United States
| | - Yothin Rakvongthai
- Massachusetts General Hospital, Division of Nuclear Medicine and Molecular Imaging, Gordon Center for Medical Imaging, 55 Fruit Street, White 427, Boston, Massachusetts 02114, United States
- Chulalongkorn University, Department of Radiology, Faculty of Medicine, 1873 Rama 4 Road, Pathumwan, Bangkok 10330, Thailand
| | - Andrew F. Laine
- Columbia University, Department of Biomedical Engineering, 1210 Amsterdam Avenue, New York, New York 10027, United States
| | - Quanzheng Li
- Massachusetts General Hospital, Division of Nuclear Medicine and Molecular Imaging, Gordon Center for Medical Imaging, 55 Fruit Street, White 427, Boston, Massachusetts 02114, United States
- Harvard Medical School, Department of Radiology, 25 Shattuck Street, Boston, Massachusetts 02115, United States
| |
Collapse
|
14
|
Rigie DS, La Rivière PJ. Optimizing spectral CT parameters for material classification tasks. Phys Med Biol 2016; 61:4599-622. [PMID: 27227430 DOI: 10.1088/0031-9155/61/12/4599] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
In this work, we propose a framework for optimizing spectral CT imaging parameters and hardware design with regard to material classification tasks. Compared with conventional CT, many more parameters must be considered when designing spectral CT systems and protocols. These choices will impact material classification performance in a non-obvious, task-dependent way with direct implications for radiation dose reduction. In light of this, we adapt Hotelling Observer formalisms typically applied to signal detection tasks to the spectral CT, material-classification problem. The result is a rapidly computable metric that makes it possible to sweep out many system configurations, generating parameter optimization curves (POC's) that can be used to select optimal settings. The proposed model avoids restrictive assumptions about the basis-material decomposition (e.g. linearity) and incorporates signal uncertainty with a stochastic object model. This technique is demonstrated on dual-kVp and photon-counting systems for two different, clinically motivated material classification tasks (kidney stone classification and plaque removal). We show that the POC's predicted with the proposed analytic model agree well with those derived from computationally intensive numerical simulation studies.
Collapse
Affiliation(s)
- D S Rigie
- Department of Radiology, University of Chicago, Chicago, IL, USA
| | | |
Collapse
|
15
|
Sánchez AA, Sidky EY, Pan X. Use of the Hotelling observer to optimize image reconstruction in digital breast tomosynthesis. J Med Imaging (Bellingham) 2015; 3:011008. [PMID: 26702408 DOI: 10.1117/1.jmi.3.1.011008] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2015] [Accepted: 11/16/2015] [Indexed: 11/14/2022] Open
Abstract
We propose an implementation of the Hotelling observer that can be applied to the optimization of linear image reconstruction algorithms in digital breast tomosynthesis. The method is based on considering information within a specific region of interest, and it is applied to the optimization of algorithms for detectability of microcalcifications. Several linear algorithms are considered: simple back-projection, filtered back-projection, back-projection filtration, and [Formula: see text]-tomography. The optimized algorithms are then evaluated through the reconstruction of phantom data. The method appears robust across algorithms and parameters and leads to the generation of algorithm implementations which subjectively appear optimized for the task of interest.
Collapse
Affiliation(s)
- Adrian A Sánchez
- University of Chicago , Department of Radiology, 5841 South Maryland Avenue, Chicago 60615, United States
| | - Emil Y Sidky
- University of Chicago , Department of Radiology, 5841 South Maryland Avenue, Chicago 60615, United States
| | - Xiaochuan Pan
- University of Chicago , Department of Radiology, 5841 South Maryland Avenue, Chicago 60615, United States ; University of Chicago , Department of Radiation and Cellular Oncology, 5758 South Maryland Avenue, Chicago 60615, United States
| |
Collapse
|
16
|
Gang GJ, Stayman JW, Ehtiati T, Siewerdsen JH. Task-driven image acquisition and reconstruction in cone-beam CT. Phys Med Biol 2015; 60:3129-50. [PMID: 25803361 PMCID: PMC4539970 DOI: 10.1088/0031-9155/60/8/3129] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
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.
Collapse
Affiliation(s)
- Grace J. Gang
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD 21205, USA
| | - J. Webster Stayman
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD 21205, USA
| | - Tina Ehtiati
- Siemens Healthcare AX Division, Erlangen, Germany
| | - Jeffrey H. Siewerdsen
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD 21205, USA
| |
Collapse
|
17
|
Sanchez AA, Sidky EY, Pan X. Region of interest based Hotelling observer for computed tomography with comparison to alternative methods. J Med Imaging (Bellingham) 2014; 1:031010. [PMID: 25685825 DOI: 10.1117/1.jmi.1.3.031010] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We compare several approaches to estimation of Hotelling observer (HO) performance in x-ray computed tomography (CT). We consider the case where the signal of interest is small so that the reconstructed image can be restricted to a small region of interest (ROI) surrounding the signal. This reduces the dimensionality of the image covariance matrix so that direct computation of HO metrics within the ROI is feasible. We propose that this approach is directly applicable to systems optimization in CT; however, many alternative approaches exist, which make computation of HO performance tractable through a range of approximations, assumptions, or estimation strategies. Here, we compare several of these methods, including the use of Laguerre-Gauss channels, discrete Fourier domain computation of the HO (which assumes noise stationarity), and two approaches to HO estimation through samples of noisy images. Since our method computes HO performance exactly within an ROI, this allows us to investigate the validity of the assumptions inherent in various common approaches to HO estimation, such as the stationarity assumption in the case of the discrete Fourier transform domain method.
Collapse
Affiliation(s)
- Adrian A Sanchez
- The University of Chicago, Department of Radiology, 5841 South Maryland Avenue, Chicago, Illinois 60615, United States
| | - Emil Y Sidky
- The University of Chicago, Department of Radiology, 5841 South Maryland Avenue, Chicago, Illinois 60615, United States
| | - Xiaochuan Pan
- The University of Chicago, Department of Radiology, 5841 South Maryland Avenue, Chicago, Illinois 60615, United States ; The University of Chicago, Department of Radiation and Cellular Oncology, 5758 South Maryland Avenue, Chicago, Illinois 60615, United States
| |
Collapse
|