PubMed Health⌕ Search

Biomedical subjects

H H Barrett

Publications and source records attributed to H H Barrett.

At least 19 recordsLinked to original sources

Ideal-observer performance under signal and background uncertainty.

We use the performance of the Bayesian ideal observer as a figure of merit for hardware optimization because this observer makes optimal use of signal-detection information. Due to the high dimensionality of certain integrals that need to be evaluated, it is difficult to compute the ideal observer test statistic, the likelihood ratio, when background variability is taken into account. Methods have been developed in our laboratory for performing this computation for fixed signals in random backgrounds. In this work, we extend these computational methods to compute the likelihood ratio in the case where both the backgrounds and the signals are random with known statistical properties. We are able to write the likelihood ratio as an integral over possible backgrounds and signals, and we have developed Markov-chain Monte Carlo (MCMC) techniques to estimate these high-dimensional integrals. We can use these results to quantify the degradation of the ideal-observer performance when signal uncertainties are present in addition to the randomness of the backgrounds. For background uncertainty, we use lumpy backgrounds. We present the performance of the ideal observer under various signal-uncertainty paradigms with different parameters of simulated parallel-hole collimator imaging systems. We are interested in any change in the rankings between different imaging systems under signal and background uncertainty compared to the background-uncertainty case. We also compare psychophysical studies to the performance of the ideal observer.

Algorithms↗

Chemical-shift imaging utilizing the positional shifts along the readout gradient direction.

In this work, we describe a method that uses the linear phase acquired during the readout period due to chemical shift to generate individual magnetic resonance (MR) images of chemically shifted species. The method utilizes sets of Fourier (or k-space) data acquired with different directions of the readout gradient and a postprocessing algorithm to generate chemical shift images. The methodology is developed for both Cartesian data acquisition and for radial data acquisition. The method is presented here for two chemically shifted species but it can be extended to more species. In this work, we present the theory, show the results in phantoms and in human images, and discuss the artifacts and signal-to-noise ratio of the images obtained with the technique.

Adipose Tissue↗

Human- and model-observer performance in ramp-spectrum noise: effects of regularization and object variability.

We consider detection of a nodule signal profile in noisy images meant to roughly simulate the statistical properties of tomographic image reconstructions in nuclear medicine. The images have two sources of variability arising from quantum noise from the imaging process and anatomical variability in the ensemble of objects being imaged. Both of these sources of variability are simulated by a stationary Gaussian random process. Sample images from this process are generated by filtering white-noise images. Human-observer performance in several signal-known-exactly detection tasks is evaluated through psychophysical studies by using the two-alternative forced-choice method. The tasks considered investigate parameters of the images that influence both the signal profile and pixel-to-pixel correlations in the images. The effect of low-pass filtering is investigated as an approximation to regularization implemented by image-reconstruction algorithms. The relative magnitudes of the quantum and the anatomical variability are investigated as an approximation to the effects of exposure time. Finally, we study the effect of the anatomical correlations in the form of an anatomical slope as an approximation to the effects of different tissue types. Human-observer performance is compared with the performance of a number of model observers computed directly from the ensemble statistics of the images used in the experiments for the purpose of finding predictive models. The model observers investigated include a number of nonprewhitening observers, the Hotelling observer (which is equivalent to the ideal observer for these studies), and six implementations of channelized-Hotelling observers. The human observers demonstrate large effects across the experimental parameters investigated. In the regularization study, performance exhibits a mild peak at intermediate levels of regularization before degrading at higher levels. The exposure-time study shows that human observers are able to detect ever more subtle lesions at increased exposure times. The anatomical slope study shows that human-observer performance degrades as anatomical variability extends into higher spatial frequencies. Of the observers tested, the channelized-Hotelling observers best capture the features of the human data.

Artifacts↗

Reconstruction of two- and three-dimensional images from synthetic-collimator data.

A novel SPECT collimation method, termed the synthetic collimator, is proposed. The synthetic collimator employs a multiple-pinhole aperture and a high-resolution detector. The problem of multiplexing, normally associated with multiple pinholes, is reduced by obtaining projections at a number of pinhole-detector distances. Projections with little multiplexing are collected at small pinhole-detector distances and high-resolution projections are collected at greater pinhole-detector distances. These projections are then reconstructed using the ML-EM algorithm. It is demonstrated through computer simulations that the synthetic collimator has superior resolution properties to a high-resolution parallel-beam (HRPB) collimator and a specially built ultra-high-resolution parallel-beam (UHRPB) collimator designed for our 0.38-mm pixel CdZnTe detectors. It is also shown that reconstructing images in three dimensions is superior to reconstructing them in two dimensions. The advantages of a high-resolution synthetic collimator over the parallel-hole collimators are apparently reduced in the presence of statistical noise. However, a high-sensitivity synthetic collimator was designed which again shows superior properties to the parallel-hole collimators. Finally, it is demonstrated that, for the cases studied, high-resolution detectors are necessary for the proper functionality of the synthetic collimator.

Algorithms↗

Reconstruction of MR images from data acquired on a general nonregular grid by pseudoinverse calculation.

A minimum-norm least-squares image-reconstruction method for the reconstruction of magnetic resonance images from non-Cartesian sampled data is proposed. The method is based on a general formalism for continuous-to-discrete mapping and pseudoinverse calculation. It does not involve any regridding or interpolation of the data and therefore the methodology differs fundamentally from existing regridding-based methods. Moreover, the method uses a continuous representation of objects in the image domain instead of a discretized representation. Simulations and experiments show the possibilities of the method in both radial and spiral imaging. Simulations revealed that minimum-norm least-squares image reconstruction can result in a drastic decrease of artifacts compared with regridding-based reconstruction. Besides, both in vivo and phantom experiments showed that minimum-norm least-squares image reconstruction leads to contrast improvement and increased signal-to-noise ratio compared with image reconstruction based on regridding. As an appendix, an analytical calculation of the raw data corresponding to the well-known Shepp and Logan software head phantom is presented.

Humans↗

A method for approximating the density of maximum-likelihood and maximum a posteriori estimates under a Gaussian noise model.

The performance of maximum-likelihood (ML) and maximum a posteriori (MAP) estimates in non-linear problems at low data SNR is not well predicted by the Cramér-Rao or other lower bounds on variance. In order to better characterize the distribution of ML and MAP estimates under these conditions, we derive a point approximation to density values of the conditional distribution of such estimates. In an example problem, this approximate distribution captures the essential features of the distribution of ML estimates in the presence of Gaussian-distributed noise.

Algorithms↗

List-mode likelihood: EM algorithm and image quality estimation demonstrated on 2-D PET.

Using a theory of list-mode maximum-likelihood (ML) source reconstruction presented recently by Barrett et al., this paper formulates a corresponding expectation-maximization (EM) algorithm, as well as a method for estimating noise properties at the ML estimate. List-mode ML is of interest in cases where the dimensionality of the measurement space impedes a binning of the measurement data. It can be advantageous in cases where a better forward model can be obtained by including more measurement coordinates provided by a given detector. Different figures of merit for the detector performance can be computed from the Fisher information matrix (FIM). This paper uses the observed FIM, which requires a single data set, thus, avoiding costly ensemble statistics. The proposed techniques are demonstrated for an idealized two-dimensional (2-D) positron emission tomography (PET) [2-D PET] detector. We compute from simulation data the improved image quality obtained by including the time of flight of the coincident quanta.

Algorithms↗

Phase-shifting interferometer/ellipsometer capable of measuring the complex index of refraction and the surface profile of a test surface.

A novel interferometer based on a conventional phase-shifting design is presented. This interferometer is capable of measuring both the real and the imaginary parts of the complex index of refraction and the surface profile of a test surface, n, k, and h, respectively. Maximum-likelihood-estimation theory is shown to be a viable means of extracting the three parameters of interest from the measured data. A Monte Carlo simulation showed limited success in estimating the complex index parameters. The results exhibited bias or deviation from the true values for the system configuration examined. The estimate on the surface profile showed excellent agreement with the true value, although the error in the estimate was an order of magnitude worse than in the case in which only the surface profile is to be estimated.

Algorithms↗

Objective assessment of image quality. III. ROC metrics, ideal observers, and likelihood-generating functions.

We continue the theme of previous papers [J. Opt. Soc. Am. A 7, 1266 (1990); 12, 834 (1995)] on objective (task-based) assessment of image quality. We concentrate on signal-detection tasks and figures of merit related to the ROC (receiver operating characteristic) curve. Many different expressions for the area under an ROC curve (AUC) are derived for an arbitrary discriminant function, with different assumptions on what information about the discriminant function is available. In particular, it is shown that AUC can be expressed by a principal-value integral that involves the characteristic functions of the discriminant. Then the discussion is specialized to the ideal observer, defined as one who uses the likelihood ratio (or some monotonic transformation of it, such as its logarithm) as the discriminant function. The properties of the ideal observer are examined from first principles. Several strong constraints on the moments of the likelihood ratio or the log likelihood are derived, and it is shown that the probability density functions for these test statistics are intimately related. In particular, some surprising results are presented for the case in which the log likelihood is normally distributed under one hypothesis. To unify these considerations, a new quantity called the likelihood-generating function is defined. It is shown that all moments of both the likelihood and the log likelihood under both hypotheses can be derived from this one function. Moreover, the AUC can be expressed, to an excellent approximation, in terms of the likelihood-generating function evaluated at the origin. This expression is the leading term in an asymptotic expansion of the AUC; it is exact whenever the likelihood-generating function behaves linearly near the origin. It is also shown that the likelihood-generating function at the origin sets a lower bound on the AUC in all cases.

Area Under Curve↗

Assessing and optimizing the performance of a phase-shifting interferometer capable of measuring the complex index of refraction and the surface profile of a test surface.

A novel interferometer based upon a conventional phase-shifting design is further investigated. This interferometer is capable of measuring both the real and imaginary parts of the complex index of refraction and the surface profile of a test surface. Maximum-likelihood estimation theory is shown to be an effective means of extracting the three parameters of interest from the measured data. Cramér-Rao lower bounds are introduced as a means of quantitatively assessing the performance of the system. Furthermore, it is shown that as the design parameters are optimized, the results approach the theoretical performance limit. We conclude by developing the underlying theory behind the relationship of the complex-index-of-refraction estimates to the surface-profile estimate.

Equipment Design↗

List-mode likelihood.

As photon-counting imaging systems become more complex, there is a trend toward measuring more attributes of each individual event. In various imaging systems the attributes can include several position variables, time variables, and energies. If more than about four attributes are measured for each event, it is not practical to record the data in an image matrix. Instead it is more efficient to use a simple list where every attribute is stored for every event. It is the purpose of this paper to discuss the concept of likelihood for such list-mode data. We present expressions for list-mode likelihood with an arbitrary number of attributes per photon and for both preset counts and preset time. Maximization of this likelihood can lead to a practical reconstruction algorithm with list-mode data, but that aspect is covered in a separate paper [IEEE Trans. Med. Imaging (to be published)]. An expression for lesion detectability for list-mode data is also derived and compared with the corresponding expression for conventional binned data.

Image Processing, Computer-Assisted↗

Objective assessment of image quality. II. Fisher information, Fourier crosstalk, and figures of merit for task performance.

Figures of merit for image quality are derived on the basis of the performance of mathematical observers on specific detection and estimation tasks. The tasks include detection of a known signal superimposed on a known background, detection of a known signal on a random background, estimation of Fourier coefficients of the object, and estimation of the integral of the object over a specified region of interest. The chosen observer for the detection tasks is the ideal linear discriminant, which we call the Hotelling observer. The figures of merit are based on the Fisher information matrix relevant to estimation of the Fourier coefficients and the closely related Fourier crosstalk matrix introduced earlier by Barrett and Gifford [Phys. Med. Biol. 39, 451 (1994)]. A finite submatrix of the infinite Fisher information matrix is used to set Cramer-Rao lower bounds on the variances of the estimates of the first N Fourier coefficients. The figures of merit for detection tasks are shown to be closely related to the concepts of noise-equivalent quanta (NEQ) and generalized NEQ, originally derived for linear, shift-invariant imaging systems and stationary noise. Application of these results to the design of imaging systems is discussed.

Computer Simulation↗

Cone-beam tomography with discrete data sets.

Sufficiency conditions for cone-beam data are well known for the case of continuous data collection along a cone-vertex curve with continuous detectors. These continuous conditions are inadequate for real-world data where discrete vertex geometries and discrete detector arrays are used. In this paper we present a theoretical formulation of cone-beam tomography with arbitrary discrete arrays of detectors and vertices. The theory models the imaging system as a linear continuous-to-discrete mapping and represents the continuous object exactly as a Fourier series. The reconstruction problem is posed as the estimation of some subset of the Fourier coefficients. The main goal of the theory is to determine which Fourier coefficients can be reliably determined from the data delivered by a specific discrete design. A Fourier component will be well determined by the data if it satisfies two conditions: it makes a strong contribution to the data, and this contribution is relatively independent of the contribution of other Fourier components. To make these considerations precise, we introduce a concept called the cross-talk matrix. A diagonal element of this matrix measures the strength of a Fourier component in the data, while an off-diagonal element quantifies the dependence or aliasing of two different components. One reasonable approach to system design is to attempt to make the diagonal elements of this matrix large and the off-diagonal elements small for some set of Fourier components. If this goal can be achieved, simple linear reconstruction algorithms are available for estimating the Fourier coefficients. To illustrate the usefulness of this approach, numerical results on the cross-talk matrix are presented for different discrete geometries derived from a continuous helical vertex orbit, and simulated images reconstructed with two linear algorithms are presented.

Algorithms↗

Noise properties of the EM algorithm: I. Theory.

The expectation-maximization (EM) algorithm is an important tool for maximum-likelihood (ML) estimation and image reconstruction, especially in medical imaging. It is a non-linear iterative algorithm that attempts to find the ML estimate of the object that produced a data set. The convergence of the algorithm and other deterministic properties are well established, but relatively little is known about how noise in the data influences noise in the final reconstructed image. In this paper we present a detailed treatment of these statistical properties. The specific application we have in mind is image reconstruction in emission tomography, but the results are valid for any application of the EM algorithm in which the data set can be described by Poisson statistics. We show that the probability density function for the grey level at a pixel in the image is well approximated by a log-normal law. An expression is derived for the variance of the grey level and for pixel-to-pixel covariance. The variance increases rapidly with iteration number at first, but eventually saturates as the ML estimate is approached. Moreover, the variance at any iteration number has a factor proportional to the square of the mean image (though other factors may also depend on the mean image), so a map of the standard deviation resembles the object itself. Thus low-intensity regions of the image tend to have low noise. By contrast, linear reconstruction methods, such as filtered back-projection in tomography, show a much more global noise pattern, with high-intensity regions of the object contributing to noise at rather distant low-intensity regions. The theoretical results of this paper depend on two approximations, but in the second paper in this series we demonstrate through Monte Carlo simulation that the approximations are justified over a wide range of conditions in emission tomography. The theory can, therefore, be used as a basis for calculation of objective figures of merit for image quality.

Algorithms↗

Noise properties of the EM algorithm: II. Monte Carlo simulations.

In an earlier paper we derived a theoretical formulation for estimating the statistical properties of images reconstructed using the iterative ML-EM algorithm. To gain insight into this complex problem, two levels of approximation were considered in the theory. These techniques revealed the dependence of the variance and covariance of the reconstructed image noise on the source distribution, imaging system transfer function, and iteration number. In this paper a Monte Carlo approach was taken to study the noise properties of the ML-EM algorithm and to test the predictions of the theory. The study also served to evaluate the approximations used in the theory. Simulated data from phantoms were used in the Monte Carlo experiments. The ML-EM statistical properties were calculated from sample averages of a large number of images with different noise realizations. The agreement between the more exact form of the theoretical formulation and the Monte Carlo formulation was better than 10% in most cases examined, and for many situations the agreement was within the expected error of the Monte Carlo experiments. Results from the studies provide valuable information about the noise characteristics of ML-EM reconstructed images. Furthermore, the studies demonstrate the power of the theoretical and Monte Carlo approaches for investigating noise properties of statistical reconstruction algorithms.

Algorithms↗

Model observers for assessment of image quality.

Image quality can be defined objectively in terms of the performance of some "observer" (either a human or a mathematical model) for some task of practical interest. If the end user of the image will be a human, model observers are used to predict the task performance of the human, as measured by psychophysical studies, and hence to serve as the basis for optimization of image quality. In this paper, we consider the task of detection of a weak signal in a noisy image. The mathematical observers considered include the ideal Bayesian, the nonprewhitening matched filter, a model based on linear-discriminant analysis and referred to as the Hotelling observer, and the Hotelling and Bayesian observers modified to account for the spatial-frequency-selective channels in the human visual system. The theory behind these observer models is briefly reviewed, and several psychophysical studies relating to the choice among them are summarized. Only the Hotelling model with channels is mathematically tractable in all cases considered here and capable of accounting for all of these data. This model requires no adjustment of parameters to fit the data and is relatively insensitive to the details of the channel mechanism. We therefore suggest it as a useful model observer for the purpose of assessing and optimizing image quality with respect to simple detection tasks.

Discrimination, Psychological↗

A stationary hemispherical SPECT imager for three-dimensional brain imaging.

A completely stationary, hemispherical-coded aperture SPECT imaging system was designed to produce three-dimensional images of the brain. The system consisted of a hemispherical multiple-pinhole coded aperture and 20 small (100 x 100 mm crystal area) digital gamma cameras. Reconstructions and measured performance specifications from two laboratory versions of the imager are presented. The reconstructed field of view of these systems was an ellipsoidal region with semi-diameters of 100 x 100 x 50 mm. The reconstructed spatial resolution for a point source in air at the center of this field was found to be 4.8 mm FWHM and the corresponding system sensitivity was 36 cps/microCi. An analysis using an ideal-observer model indicated that the multiplexed projection data suffered a 21% degradation relative to similar, but nonmultiplexed SPECT data. Therefore, by this measure, the effective sensitivity of the brain imager was 79% of the measured value.

Brain↗

Effect of random background inhomogeneity on observer detection performance.

Many psychophysical studies of the ability of the human observer to detect a signal superimposed upon a uniform background, where both the signal and the background are known exactly, have been reported in the literature. In such cases, the ideal or the Bayesian observer is often used as a mathematical model of human performance since it can be readily calculated and is a good predictor of human performance for the task at hand. If, however, the background is spatially inhomogeneous (lumpy), the ideal observer becomes nonlinear, and its performance becomes difficult to evaluate. Since inhomogeneous backgrounds are commonly encountered in many practical applications, we have investigated the effects of background inhomogeneities on human performance. The task was detection of a two-dimensional Gaussian signal superimposed upon an inhomogeneous background and imaged through a pinhole imaging system. Poisson noise corresponding to a certain exposure time and aperture size was added to the detected image. A six-point rating scale technique was used to measure human performance as a function of the strength of the nonuniformities (lumpiness) in the background, the amount of blur of the imaging system, and the amount of Poisson noise in the image. The results of this study were compared with earlier theoretical predictions by Myers et al. [J. Opt. Soc. Am. A 7, 1279 (1990)] for two observer models: the optimum linear discriminant, also known as the Hotelling observer, and a nonprewhitening matched filter. Although the efficiency of the human observer relative to the Hotelling observer was only approximately 10%, the variation in human performance with respect to varying aperture size and exposure time was well predicted by the Hotelling model. The nonprewhitening model, on the other hand, fails to predict human performance in lumpy backgrounds in this study. In particular, this model predicts that performance will saturate with increasing exposure time and drop precipitously with increasing lumpiness; neither effect is observed with human observers.

Diagnostic Imaging↗