Ideal observer estimation for binary tasks with stochastic object models

<i>Objective.</i>We propose a new formulation for ideal observers (IOs) that incorporate stochastic object models (SOMs) for data acquisition optimization.<i>Approach</i>. A data acquisition system is considered as a (possibly nonlinear) discrete-to-discrete mapping from a finite-dimensional object space,x∈Rnd, to a finite-dimensional measurement space,y∈Rm. For binary tasks, the two underlying SOMs,H0and<i>H</i>₁, are specified by two probability density functions (PDFs)p0(x),p1(x). This leads to the notion of intrinsic likelihood ratio (LR)ΛI(x)=p1(x)/p0(x)and intrinsic class separability (ICS), the latter quantifies the population separability that is independent of data acquisition. With respect to ICS, the IO employs the 'extrinsic' LRΛ(y)=pr(y|H1)/pr(y|H0)of the data and quantifies the extrinsic class separability (ECS). The difference between ICS and ECS measures the efficiency of data acquisition. We show that the extrinsic LRΛ(y)is the expectation of the intrinsic LRΛI(x), where the expectation is with respect to the posterior PDFpr(x|y,H0)underH0.<i>Main results</i>. We use two examples, one to clarify the new IO and the second to demonstrate its potential for real world applications. Specifically, we apply the new IO to spectral optimization in dual-energy CT projection domain material decomposition (pMD), for which SOMs are used to describe variability of basis material line integrals. The performance rank orders obtained by IO agree with physics predictions.<i>Significance.</i>The main computation in the new IO involves sampling from the posterior PDFpr(x|y,H0), which are similar to (fully) Bayesian reconstruction. Thus our IO computation is amenable to standard techniques already familiar to CT researchers. The example of dual-energy pMD serves as a prototype for other spectral optimization problems, e.g., for photon counting CT or multi-energy CT with multi-layer detectors.