PubMed Health⌕ Search

Biomedical subjects

Yoram Bresler

Publications and source records attributed to Yoram Bresler.

4 recordsLinked to original sources

Asymptotic global confidence regions for 3-D parametric shape estimation in inverse problems.

This paper derives fundamental performance bounds for statistical estimation of parametric surfaces embedded in R3. Unlike conventional pixel-based image reconstruction approaches, our problem is reconstruction of the shape of binary or homogeneous objects. The fundamental uncertainty of such estimation problems can be represented by global confidenceregions, which facilitate geometric inference and optimization ofthe imaging system. Compared to our previous work on global confidence region analysis for curves [two-dimensional (2-D) shapes], computation of the probability that the entire surface estimate lies within the confidence region is more challenging because a surface estimate is an inhomogeneous random field continuously indexed by a 2-D variable. We derive an asymptotic lower bound to this probability by relating it to the exceedence probability of a higher dimensional Gaussian random field, which can, in turn, be evaluated using the tube formula due to Sun. Simulation results demonstrate the tightness of the resulting bound and the usefulness of the three-dimensional global confidence region approach.

Algorithms↗

O(N3 log N) backprojection algorithm for the 3-D radon transform.

We present a novel backprojection algorithm for three-dimensional (3-D) radon transform data that requires O(N3 log2 N) operations for reconstruction of an N x N x N volume from O(N2) plane-integral projections. Our algorithm uses a hierarchical decomposition of the 3-D radon transform to recursively decompose the backprojection operation. Simulations are presented demonstrating reconstruction quality comparable to the standard filtered backprojection, which requires O(N5) computations under the same circumstances.

Algorithms↗

Electromagnetic impedance tomography (EMIT): a new method for impedance imaging.

We propose a new impedance imaging method, electromagnetic impedance tomography (EMIT), in which the boundary electric potential measurements in electrical impedance tomography (EIT) are augmented by measurements of the exterior magnetic field induced by the currents excited in the object by the standard EIT procedures. These magnetic measurements can be obtained reliably and inexpensively by simple pickup coils located around the imaged cross section. We derive expressions for the forward problem and for the Jacobian of the measurements, and propose an iterative reconstruction algorithm using a squared error cost function. The performance of EMIT and EIT is compared in numerical simulations using a finite-element model for the conductivity distribution of several phantoms. Evaluation of the rank and condition of the Jacobian demonstrates that the additional magnetic measurements provided by a few pickup coils in EMIT turn an underdetermined EIT problem into a well-posed one with reasonable condition, or significantly improve the conditioning of the EIT problem when it is already fully determined. Reconstructions of various phantoms reveal that EMIT provides particularly significant visual and quantitative improvement (threefold to tenfold reduction in the root-mean-squared error) in the sensitivity at the center of the object, which is the area most difficult to image using EIT.

Computer Simulation↗

Level-set algorithm for the reconstruction of functional activation in near-infrared spectroscopic imaging.

We introduce a new algorithm for the reconstruction of functional brain activations from near-infrared spectroscopic imaging (NIRSI) data. While NIRSI offers remarkable biochemical specificity, the attainable spatial resolution with this technique is rather limited, mainly due to the highly scattering nature of brain tissue and the low number of measurement channels. Our approach exploits the support-limited (spatially concentrated) nature of the activations to make the reconstruction problem well-posed. The new algorithm considers both the support and the function values of the activations as unknowns and estimates them from the data. The support of the activations is represented using a level-set scheme. We use a two-step alternating iterative scheme to solve for the activations. Since our approach uses the inherent nature of functional activations to make the problem well-posed, it provides reconstructions with better spatial resolution, fewer artifacts, and is more robust to noise than existing techniques. Numerical simulations and experimental data indicate a significant improvement in the quality (resolution and robustness to noise) over standard techniques such as truncated conjugate gradients (TCG) and simultaneous iterative reconstruction technique (SIRT) algorithms. Furthermore, results on experimental data obtained from simultaneous functional magnetic resonance imaging (fMRI) and optical measurements show much closer agreement of the optical reconstruction using the new approach with fMRI images than TCG and SIRT.

Algorithms↗