Thermoacoustic tomography with detectors on an open curve: an efficient reconstruction algorithm

Practical applications of thermoacoustic tomography require numerical inversion of the spherical mean Radon transform with the centers of integration spheres occupying an open surface. Solution of this problem is needed (both in 2-D and 3-D) because frequently the region of interest cannot be completely surrounded by the detectors, as it happens, for example, in breast imaging. We present an efficient numerical algorithm for solving this problem in 2-D (similar methods are applicable in the 3-D case). Our method is based on the numerical approximation of plane waves by certain single layer potentials related to the acquisition geometry. After the densities of these potentials have been precomputed, each subsequent image reconstruction has the complexity of the regular filtration backprojection algorithm for the classical Radon transform. The peformance of the method is demonstrated in several numerical examples: one can see that the algorithm produces very accurate reconstructions if the data are accurate and sufficiently well sampled, on the other hand, it is sufficiently stable with respect to noise in the data

Efficient Compressive Sampling of Spatially Sparse Fields in Wireless Sensor Networks

Wireless sensor networks (WSN), i.e. networks of autonomous, wireless sensing nodes spatially deployed over a geographical area, are often faced with acquisition of spatially sparse fields. In this paper, we present a novel bandwidth/energy efficient CS scheme for acquisition of spatially sparse fields in a WSN. The paper contribution is twofold. Firstly, we introduce a sparse, structured CS matrix and we analytically show that it allows accurate reconstruction of bidimensional spatially sparse signals, such as those occurring in several surveillance application. Secondly, we analytically evaluate the energy and bandwidth consumption of our CS scheme when it is applied to data acquisition in a WSN. Numerical results demonstrate that our CS scheme achieves significant energy and bandwidth savings wrt state-of-the-art approaches when employed for sensing a spatially sparse field by means of a WSN.Comment: Submitted to EURASIP Journal on Advances in Signal Processin

Fast Hadamard transforms for compressive sensing of joint systems: measurement of a 3.2 million-dimensional bi-photon probability distribution

We demonstrate how to efficiently implement extremely high-dimensional compressive imaging of a bi-photon probability distribution. Our method uses fast-Hadamard-transform Kronecker-based compressive sensing to acquire the joint space distribution. We list, in detail, the operations necessary to enable fast-transform-based matrix-vector operations in the joint space to reconstruct a 16.8 million-dimensional image in less than 10 minutes. Within a subspace of that image exists a 3.2 million-dimensional bi-photon probability distribution. In addition, we demonstrate how the marginal distributions can aid in the accuracy of joint space distribution reconstructions

Non-convex image reconstruction via Expectation Propagation

Tomographic image reconstruction can be mapped to a problem of finding solutions to a large system of linear equations which maximize a function that includes \textit{a priori} knowledge regarding features of typical images such as smoothness or sharpness. This maximization can be performed with standard local optimization tools when the function is concave, but it is generally intractable for realistic priors, which are non-concave. We introduce a new method to reconstruct images obtained from Radon projections by using Expectation Propagation, which allows us to reframe the problem from an Bayesian inference perspective. We show, by means of extensive simulations, that, compared to state-of-the-art algorithms for this task, Expectation Propagation paired with very simple but non log-concave priors, is often able to reconstruct images up to a smaller error while using a lower amount of information per pixel. We provide estimates for the critical rate of information per pixel above which recovery is error-free by means of simulations on ensembles of phantom and real images.Comment: 12 pages, 6 figure

Hierarchical Bayesian sparse image reconstruction with application to MRFM

This paper presents a hierarchical Bayesian model to reconstruct sparse images when the observations are obtained from linear transformations and corrupted by an additive white Gaussian noise. Our hierarchical Bayes model is well suited to such naturally sparse image applications as it seamlessly accounts for properties such as sparsity and positivity of the image via appropriate Bayes priors. We propose a prior that is based on a weighted mixture of a positive exponential distribution and a mass at zero. The prior has hyperparameters that are tuned automatically by marginalization over the hierarchical Bayesian model. To overcome the complexity of the posterior distribution, a Gibbs sampling strategy is proposed. The Gibbs samples can be used to estimate the image to be recovered, e.g. by maximizing the estimated posterior distribution. In our fully Bayesian approach the posteriors of all the parameters are available. Thus our algorithm provides more information than other previously proposed sparse reconstruction methods that only give a point estimate. The performance of our hierarchical Bayesian sparse reconstruction method is illustrated on synthetic and real data collected from a tobacco virus sample using a prototype MRFM instrument.Comment: v2: final version; IEEE Trans. Image Processing, 200

Three-dimensional magnetization structures revealed with X-ray vector nanotomography

In soft ferromagnetic materials, the smoothly varying magnetization leads to the formation of fundamental patterns such as domains, vortices and domain walls<sup>1</sup>. These have been studied extensively in thin films of thicknesses up to around 200 nanometres, in which the magnetization is accessible with current transmission imaging methods that make use of electrons or soft X-rays. In thicker samples, however, in which the magnetization structure varies throughout the thickness and is intrinsically three dimensional, determining the complex magnetic structure directly still represents a challenge<sup>1, 3</sup>. We have developed hard-X-ray vector nanotomography with which to determine the three-dimensional magnetic configuration at the nanoscale within micrometre-sized samples. We imaged the structure of the magnetization within a soft magnetic pillar of diameter 5 micrometres with a spatial resolution of 100 nanometres and, within the bulk, observed a complex magnetic configuration that consists of vortices and antivortices that form cross-tie walls and vortex walls along intersecting planes. At the intersections of these structures, magnetic singularities—Bloch points—occur. These were predicted more than fifty years ago<sup>4</sup> but have so far not been directly observed. Here we image the three-dimensional magnetic structure in the vicinity of the Bloch points, which until now has been accessible only through micromagnetic simulations, and identify two possible magnetization configurations: a circulating magnetization structure<sup>5</sup> and a twisted state that appears to correspond to an ‘anti-Bloch point’. Our imaging method enables the nanoscale study of topological magnetic structures<sup>6</sup> in systems with sizes of the order of tens of micrometres. Knowledge of internal nanomagnetic textures is critical for understanding macroscopic magnetic properties and for designing bulk magnets for technological applications<sup>7</sup>