Search PubMedSearch

Biomedical subjects

H H Barrett

Publications and source records attributed to H H Barrett.

At least 19 recordsLinked to original sources

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

Parallel simulated annealing for emission tomography.

A method for implementing simulated annealing in parallel to speed up the execution of emission tomography (ET) image reconstruction is presented. A high degree of parallelism can be attained by using a parallel-acceptance partitioning strategy, in which perturbations to subsets of the estimate are evaluated in parallel. However because the point spread function in ET imaging systems is globally dependent, processors cannot update the current estimate independently. Consequently, processors must be synchronized each time a perturbation is accepted to avoid introducing error. This can produce excessive communications overhead, especially when the acceptance rate is high. In this paper an energy function is constructed to reduce the synchronization requirements by using a reformulation of the log-likelihood function from the expectation maximization (EM) algorithm. The approach is to change the global dependence in the energy function from the current estimate to the estimate generated during the last iteration. The synchronization requirements for guaranteed convergence are then significantly reduced from once per acceptance to once per iteration. This parallel implementation on 54 Inmos T800 transputers connected in a ring topology resulted in execution times that were almost 50 times faster than on a VAX 8600.

Algorithms

Ideal versus human observer for long-tailed point spread functions: does deconvolution help?

The ideal observer represents a Bayesian approach to performing detection tasks. Since such tasks are frequently used as a prototype tasks for radiological imaging systems, the detectability measured at the output of an ideal detector can be used as a figure of merit to characterize the imaging system. For the detectability achieved by the ideal observer to be a good figure of merit, it should predict the ability of the human observer to perform the same detection task. Of great general interest, especially to the medical community, are imaging devices with long-tailed point spread functions (PSFs). Such PSFs may occur due to septal penetration in collimators, veiling glare in image intensifiers or scattered radiation in the body. We have investigated the effect that this type of PSF has on human visual signal detection and whether any improvement in performance can be gained by deconvolving the tails of the PSF. For the ideal observer, it is straightforward to show that the performance is independent of any linear, invertible deconvolution filter. Our psychophysical studies show, however, that performance of the human observer is indeed improved by deconvolution. The ideal observer is, therefore, not a good predictor of human observer performance for detection of a signal imaged through a long-tailed PSF. We offer some explanations for this discrepancy by using some characteristics of the visual process and suggest a standard of comparison for the human observer that takes into account these characteristics. A look at the performance of the non-prewhitening (npw) ideal observer, before and after deconvolution, also brings some good insight into this study.

Bayes Theorem

Probability modelling of a surgical probe for tumour detection.

Probability functions for the output counts from a radiation detector probe are needed to implement Bayesian detection strategies or to assess performance of the probe. This paper presents methods for simulating a surgical probe designed for tumour detection to obtain statistical information for modelling probability functions of output data. Statistical models of pharmaceutical uptake in normal organs and tumours were estimated from animal and human data, and these models were combined with a digitised human torso phantom to create a large set of simulated patients. With the simulated patients and with a spatial map of the probe response, computer simulations of intraoperative probe measurements provided a large set of simulated probe data. Probability models derived from these data using maximum-likelihood methods helped to formulate the detection strategy and to evaluate the performance of the surgical probe.

Computer Simulation

Objective assessment of image quality: effects of quantum noise and object variability.

A number of task-specific approaches to the assessment of image quality are treated. Both estimation and classification tasks are considered, but only linear estimators or classifiers are permitted. Performance on these tasks is limited by both quantum noise and object variability, and the effects of postprocessing or image-reconstruction algorithms are explicitly included. The results are expressed as signal-to-noise ratios (SNR's). The interrelationships among these SNR's are considered, and an SNR for a classification task is expressed as the SNR for a related estimation task times four factors. These factors show the effects of signal size and contrast, conspicuity of the signal, bias in the estimation task, and noise correlation. Ways of choosing and calculating appropriate SNR's for system evaluation and optimization are also discussed.

Algorithms

Aperture optimization for emission imaging: effect of a spatially varying background.

A method for optimizing the aperture size in emission imaging is presented that takes into account limitations due to the Poisson nature of the detected radiation stream as well as the conspicuity limitation imposed by a spatially varying background. System assessment is based on the calculated performance of two model observers: the best linear observer, also called the Hotelling observer, and the nonprewhitening matched-filter observer. The tasks are the detection of a Gaussian signal and the discrimination of a single from a double Gaussian signal. When the background is specified, detection is optimized by enlarging the aperture; an inhomogeneous background results in an optimum aperture size matched naturally to the signal. The discrimination task has a finite optimum aperture for a flat background; a nonuniform background drives the optimum toward still-finer resolution.

Image Processing, Computer-Assisted

A full-field modular gamma camera.

A modular gamma ray camera is described that gives useful image information over its entire crystal face. The lack of dead area on the periphery of the camera is made possible by a unique application of digital electronics and optimal position estimation using maximum likelihood (ML) estimates. The ML estimates are calculated directly from photomultiplier tube responses and stored in a lookup table, so the restriction of calculating the position estimates in separate circuitry is removed. Each module is designed to be optically and electronically independent, so that many modules can be combined in a large system. Results from a prototypical module, which has an active crystal area of 10 cm X 10 cm, are presented.

Equipment Design

Comparison of in vivo scintillation probes and gamma cameras for detection of small, deep tumours.

An in vivo probe was compared with a gamma camera for the task of detecting radiolabelled tumour models in a water phantom. The probe, which contained a 1 cm diameter NaI(T1) detector, was designed and used by us for surgical staging studies of gynaecological patients. Tumours were spherical sources of different sizes and activities per unit volume of cobalt-57. The phantom was a tank of water (451) containing dissolved radioactivity to simulate background activity. Detector-to-source separations and tank depths were also varied. Camera images of 10 min duration were compared with probe counts of 15 s. For each configuration a large number of source and background runs were analysed using an ideal-observer ROC technique. Area under the ROC curve was used as the figure of merit. Results show that for approximately uniform background the probe should perform substantially better than the gamma camera in detecting small, deep tumours provided that the probe can be manoeuvred to within a few centimetres of the tumour. Mathematical modelling of our results indicates that this conclusion is not dependent on radiopharmaceutical dose, tumour uptake or camera type.

Cobalt Radioisotopes

Linear estimation theory applied to the evaluation of a priori information and system optimization in coded-aperture imaging.

Linear estimation theory is developed in the context of object reconstruction from data obtained by a general shift-variant imaging system. The formalism adopts nonstationary first- and second-order statistics of the object and noise classes as priori information. In addition, a metric for system optimization that depends on the a priori information is presented. The role of this a priori information as derived from several different training sets is then studied with respect to reconstruction performance for various noise levels in the data, using a tomographic coded-aperture system as the model. In a separate experiment, a simple coded-aperture system is optimized to a particular object class, and the results are compared with those from an earlier optimization experiment.

Diagnostic Imaging

Dual-detector probe for surgical tumor staging.

A hand-held, dual-detector probe has been developed for surgical tumor staging. This dual probe simultaneously monitors counts from a possible tumor site along with counts from adjacent normal tissue using two concentric, collimated scintillation detectors. A comparison of counts from the detectors can distinguish a small tumor directly in front of the probe from variations in background activity. The probe was tested in computer simulations of surgical staging of metastases to para-aortic and iliac lymph nodes using a spatial response map of the probe, a numerical torso phantom, and organ activity data for [57Co]bleomycin in rabbits. Results show that the dual probe performs better than a single-detector probe in detecting tumors and solves the problem caused by spatial variations in the background source distribution.

Computer Simulation

Hotelling trace criterion and its correlation with human-observer performance.

The Hotelling trace criterion (HTC) is used to find a set of linear features that optimally separate two classes of objects. The objects used in our study were simulated livers with and without tumors, with noise, blur, and object variability. Using the receiver-operating-characteristic parameter da as our measure, we have found that the ability of the HTC to separate these objects into their correct classes, by detecting the presence or absence of a tumor, has a correlation of 0.988 with the ability of humans to separate the same two classes of objects. This suggests, therefore, that the HTC can be used as a figure of merit for optimizing system parameters, since it calculates a single, scalar figure of merit that has a high correlation with human-observer performance.

Humans

Finite-length line-spread function.

In this paper we derive a formula for calculating the point-spread function (PSF) of a rotationally symmetric imaging system from measurements along a line through the image of an arbitrary separable input object. An important special case of this formula is when the input object is a finite-length slit. The set of measurements in this case is called the finite-length line-spread function (FLSF). The FLSF differs from the infinite-length line-spread function (LSF) only in the assumed finite length of the line that is input into the system. This difference between the FLSF and the LSF becomes important for imaging systems for which the PSF is large in extent and in which the isoplanatic patch is relatively small. The usual LSF-to-PSF conversion formulas cannot be applied accurately to such systems.

Fourier Analysis

Addition of a channel mechanism to the ideal-observer model.

Several authors have measured the detection ability of human observers for objects in correlated (nonwhite) noise. These studies have shown that the human observer has approximately constant efficiency when compared with a nonprewhitening ideal observer. In this paper we add a frequency-selective mechanism to the ideal-observer model, similar to the channel mechanism that has been demonstrated through experiments that measure a subject's ability to detect grating stimuli. For a number of detection and discrimination tasks, the nonprewhitening ideal-observer model and the channelized ideal-observer model yield similar performance predictions. Thus both models seem equally capable of explaining a considerable body of psychophysical data, and it would be difficult to devise an experiment to determine which model is more nearly correct.

Humans

Image reconstruction from coded data: I. Reconstruction algorithms and experimental results.

Two algorithms have been developed for reconstructing objects from their coded images and a priori knowledge of the object class. Reconstructions from both algorithms are presented, but the results appear to be largely independent of the algorithm used. One of the algorithms, a Monte Carlo approach, is used to investigate the quality of the reconstruction of two- and three-dimensional objects from simulated coded-image data with respect to viewing geometry and multiplexing (mixing) of the data. The cases examined include reconstructions from data with and without signal-dependent photon noise. It is found that reconstructing from multiplexed data is not so serious a problem as reconstructing from data obtained with a limited viewing angle. Also, when photon noise is included in the data, reconstructions obtained from multiplexed data are better than those obtained from unmultiplexed data because of the higher photon count made available by multiplexing. It appears that the fidelity of a reconstruction depends much more strongly on the design of the data-taking system (the coded apertures) than on the reconstruction algorithm.

Models, Structural

Image reconstruction from coded data: II. Code design.

A strategy is given for the design of coded apertures with respect to a given class of objects that are to be imaged. Previous knowledge of the first- and second-order statistics for the object class is assumed. The object class is characterized by its Karhunen-Loève eigenvectors and eigenvalues, whereas the imaging system is characterized by its singular-value decomposition. We introduce the concept of alignment in which the aperture parameters are adjusted until the system is tuned to measure the given object class well. A mean-square-error figure of merit that indicates degree of alignment is given, and alignment is performed by standard optimization techniques. We illustrate this technique with a simple proof-of-principle experiment. These concepts are general and may be applied to any linear imaging system.

Models, Structural

Effect of noise correlation on detectability of disk signals in medical imaging.

Pixel signal-to-noise ratio is one accepted measure of image quality for predicting observer performance in medical imaging. We have found, however, that images with equal pixel signal-to-noise ratio (SNRp) but different correlation properties give quite different observer-performance measures for a simple detection experiment. The SNR at the output of an ideal detector with the ability to prewhiten the noise is also a poor predictor of human performance for disk signals in high-pass noise. We have found constant observer efficiencies for humans relative to the performance of a nonprewhitening detector for this task.

Fourier Analysis