Sparsity-regularized photon-limited imaging

  • Authors:
  • Zachary T. Harmany;Roummel F. Marcia;Rebecca M. Willett

  • Affiliations:
  • Department of Electrical and Computer Engineering, Duke University, Durham, NC;School of Natural Sciences, University of California, Merced, Merced, CA;Department of Electrical and Computer Engineering, Duke University, Durham, NC

  • Venue:
  • ISBI'10 Proceedings of the 2010 IEEE international conference on Biomedical imaging: from nano to Macro
  • Year:
  • 2010

Quantified Score

Hi-index 0.00

Visualization

Abstract

In many medical imaging applications (e.g., SPECT, PET), the data are a count of the number of photons incident on a detector array. When the number of photons is small, the measurement process is best modeled with a Poisson distribution. The problem addressed in this paper is the estimation of an underlying intensity from photon-limited projections where the intensity admits a sparse or low-complexity representation. This approach is based on recent inroads in sparse reconstruction methods inspired by compressed sensing. However, unlike most recent advances in this area, the optimization formulation we explore uses a penalized negative Poisson log-likelihood objective function with nonnegativity constraints (since Poisson intensities are naturally nonnegative). This paper describes computational methods for solving the nonnegatively constrained sparse Poisson inverse problem. In particular, the proposed approach incorporates sequential separable quadratic approximations to the log-likelihood and computationally efficient partition-based multiscale estimation methods.