Blind Poissonian reconstruction algorithm via curvelet regularization for an FTIR spectrometer

Opt Express. 2018 Sep 3;26(18):22837-22856. doi: 10.1364/OE.26.022837.

Abstract

An FTIR spectrometer often suffers from common problems of band overlap and Poisson noises. In this paper, we show that the issue of infrared (IR) spectrum degradation can be considered as a maximum a posterior (MAP) problem and solved by minimized a cost function that includes a likelihood term and two prior terms. In the MAP framework, the likelihood probability density function (PDF) is constructed based on the observed Poisson noise model. A fitted distribution of curvelet transform coefficient is used as spectral prior PDF, and the instrument response function (IRF) prior is described based on a Gauss-Markov function. Moreover, the split Bregman iteration method is employed to solve the resulting minimization problem, which highly reduces the computational load. As a result, the Poisson noises are perfectly removed, while the spectral structure information is well preserved. The novelty of the proposed method lies in its ability to estimate the IRF and latent spectrum in a joint framework, thus eliminating the degradation effects to a large extent. The reconstructed IR spectrum is more convenient for extracting the spectral feature and interpreting the unknown chemical or biological materials.