Abstract

We present a new method for computing optimized channels for estimation tasks that is feasible for high-dimensional image data. Maximum-likelihood (ML) parameter estimates are challenging to compute from high-dimensional likelihoods. The dimensionality reduction from $M$ measurements to $L$ channels is a critical advantage of channelized quadratic estimators (CQEs), since estimating likelihood moments from channelized data requires smaller sample sizes and inverting a smaller covariance matrix is easier. The channelized likelihood is then used to form ML estimates of the parameter(s). In this work we choose an imaging example in which the second-order statistics of the image data depend upon the parameter of interest: the correlation length. Correlation lengths are used to approximate background textures in many imaging applications, and in these cases an estimate of the correlation length is useful for pre-whitening. In a simulation study we compare the estimation performance, as measured by the root-mean-squared error (RMSE), of correlation length estimates from CQE and power spectral density (PSD) distribution fitting. To abide by the assumptions of the PSD method we simulate an ergodic, isotropic, stationary, and zero-mean random process. These assumptions are not part of the CQE formalism. The CQE method assumes a Gaussian channelized likelihood that can be a valid for non-Gaussian image data, since the channel outputs are formed from weighted sums of the image elements. We have shown that, for three or more channels, the RMSE of CQE estimates of correlation length is lower than conventional PSD estimates. We also show that computing CQE by using a standard nonlinear optimization method produces channels that yield RMSE within 2% of the analytic optimum. CQE estimates of anisotropic correlation length estimation are reported to demonstrate this technique on a two-parameter estimation problem.

© 2016 Optical Society of America

Method for optimizing channelized quadratic observers for binary classification of large-dimensional image datasets

M. K. Kupinski and E. Clarkson
J. Opt. Soc. Am. A 32(4) 549-565 (2015)

Maximum-likelihood methods in wavefront sensing: stochastic models and likelihood functions

Harrison H. Barrett, Christopher Dainty, and David Lara
J. Opt. Soc. Am. A 24(2) 391-414 (2007)

Channelized model observer for the detection and estimation of signals with unknown amplitude, orientation, and size

Lu Zhang, Bart Goossens, Christine Cavaro-Ménard, Patrick Le Callet, and Di Ge
J. Opt. Soc. Am. A 30(11) 2422-2432 (2013)

Estimating random signal parameters from noisy images with nuisance parameters: linear and scanning-linear methods

Meredith Kathryn Whitaker, Eric Clarkson, and Harrison H. Barrett
Opt. Express 16(11) 8150-8173 (2008)

Validating the use of channels to estimate the ideal linear observer

Brandon D. Gallas and Harrison H. Barrett
J. Opt. Soc. Am. A 20(9) 1725-1738 (2003)