Sparsity-driven ideal observer for computed medical imaging systems

The Bayesian ideal observer (IO) has been widely advocated to guide hardware optimization. However, except for special cases, computation of the IO test statistic is computationally burdensome and requires an appropriate stochastic object model that may be difficult to determine in practice. Modern reconstruction methods, referred to as sparse reconstruction methods, exploit the fact that objects of interest typically possess sparse representations and have proven to be highly effective at reconstructing images from under-sampled measurement data. Moreover, in computed imaging approaches that employ compressive sensing concepts, imaging hardware and image reconstruction are innately coupled technologies. In this work, we propose a sparsity-driven IO (SD-IO) to guide the optimization of data acquisition parameters for modern computed imaging systems. The SD-IO employs a variational Bayesian inference method to estimate the posterior distribution and calculates an approximate likelihood ratio analytically as its test statistic. Since it assumes knowledge of low-level statistical properties of the object that are related to sparsity, the SD-IO exploits the same statistical information regarding the object that is utilized by highly effective sparse image reconstruction methods. Preliminary simulation results are presented to demonstrate the feasibility of the SD-IO calculation.

[1]  M. Lustig,et al.  Compressed Sensing MRI , 2008, IEEE Signal Processing Magazine.

[2]  Bernhard Schölkopf,et al.  Optimization of k‐space trajectories for compressed sensing by Bayesian experimental design , 2010, Magnetic resonance in medicine.

[3]  Matthias W. Seeger,et al.  Large Scale Bayesian Inference and Experimental Design for Sparse Linear Models , 2011, SIAM J. Imaging Sci..

[4]  E. Sidky,et al.  Accurate image reconstruction from few-views and limited-angle data in divergent-beam CT , 2009, 0904.4495.

[5]  H H Barrett,et al.  Effect of random background inhomogeneity on observer detection performance. , 1992, Journal of the Optical Society of America. A, Optics and image science.

[6]  D. Louis Collins,et al.  Design and construction of a realistic digital brain phantom , 1998, IEEE Transactions on Medical Imaging.

[7]  Michael Elad,et al.  Sparse and Redundant Representations - From Theory to Applications in Signal and Image Processing , 2010 .

[8]  Guido Gerig,et al.  Simulation of brain tumors in MR images for evaluation of segmentation efficacy , 2009, Medical Image Anal..

[9]  David P. Wipf,et al.  Variational Bayesian Inference Techniques , 2010, IEEE Signal Processing Magazine.

[10]  Harrison H. Barrett,et al.  Foundations of Image Science , 2003, J. Electronic Imaging.

[11]  H. Barrett,et al.  Ideal-observer computation in medical imaging with use of Markov-chain Monte Carlo techniques. , 2003, Journal of the Optical Society of America. A, Optics, image science, and vision.

[12]  David L Donoho,et al.  Compressed sensing , 2006, IEEE Transactions on Information Theory.

[13]  Harrison H. Barrett Objective assessment of image qualityand its implications for system design , 2001 .

[14]  D. Field,et al.  Natural image statistics and efficient coding. , 1996, Network.