Biomedical subjects

E Somersalo

Publications and source records attributed to E Somersalo.

4 recordsLinked to original sources

PubMed2003

Statistical inversion for medical x-ray tomography with few radiographs: I. General theory.

In x-ray tomography, the structure of a three-dimensional body is reconstructed from a collection of projection images of the body. Medical CT imaging does this using an extensive set of projections from all around the body. However, in many practical imaging situations only a small number of truncated projections are available from a limited angle of view. Three-dimensional imaging using such data is complicated for two reasons: (i) typically, sparse projection data do not contain sufficient information to completely describe the 3D body, and (ii) traditional CT reconstruction algorithms, such as filtered backprojection, do not work well when applied to few irregularly spaced projections. Concerning (i), existing results about the information content of sparse projection data are reviewed and discussed. Concerning (ii), it is shown how Bayesian inversion methods can be used to incorporate a priori information into the reconstruction method, leading to improved image quality over traditional methods. Based on the discussion, a low-dose three-dimensional x-ray imaging modality is described.

Algorithms↗

PubMed2003

Statistical inversion for medical x-ray tomography with few radiographs: II. Application to dental radiology.

Diagnostic and operational tasks in dental radiology often require three-dimensional information that is difficult or impossible to see in a projection image. A CT-scan provides the dentist with comprehensive three-dimensional data. However, often CT-scan is impractical and, instead, only a few projection radiographs with sparsely distributed projection directions are available. Statistical (Bayesian) inversion is well-suited approach for reconstruction from such incomplete data. In statistical inversion, a priori information is used to compensate for the incomplete information of the data. The inverse problem is recast in the form of statistical inference from the posterior probability distribution that is based on statistical models of the projection data and the a priori information of the tissue. In this paper, a statistical model for three-dimensional imaging of dentomaxillofacial structures is proposed. Optimization and MCMC algorithms are implemented for the computation of posterior statistics. Results are given with in vitro projection data that were taken with a commercial intraoral x-ray sensor. Examples include limited-angle tomography and full-angle tomography with sparse projection data. Reconstructions with traditional tomographic reconstruction methods are given as reference for the assessment of the estimates that are based on the statistical model.

Algorithms↗

PubMed1999

Visualization of magnetoencephalographic data using minimum current estimates.

The locations of active brain areas can be estimated from the magnetic field the neural current sources produce. In this work we study a visualization method of magnetoencephalographic data that is based on minimum[symbol: see text] (1)-norm estimates. The method can represent several local or distributed sources and does not need explicit a priori information. We evaluated the performance of the method using simulation studies. In a situation resembling typical magnetoencephalographic measurement, the mean estimated source strength exceeded baseline level up to 2 cm from the simulated point-like source. The method can also visualize several sources, activated simultaneously or in a sequence, which we demonstrated by analyzing magnetic responses associated with sensory stimulation and a picture naming task.

Algorithms↗

PubMed1998

Tikhonov regularization and prior information in electrical impedance tomography.

The solution of impedance distribution in electrical impedance tomography is a nonlinear inverse problem that requires the use of a regularization method. The generalized Tikhonov regularization methods have been popular in the solution of many inverse problems. The regularization matrices that are usually used with the Tikhonov method are more or less ad hoc and the implicit prior assumptions are, thus, in many cases inappropriate. In this paper, we propose an approach to the construction of the regularization matrix that conforms to the prior assumptions on the impedance distribution. The approach is based on the construction of an approximating subspace for the expected impedance distributions. It is shown by simulations that the reconstructions obtained with the proposed method are better than with two other schemes of the same type when the prior is compatible with the true object. On the other hand, when the prior is incompatible with the true object, the method will still give reasonable estimates.

Algorithms↗