Multiple scattering of ultrasound in weakly inhomogeneous media: application to human soft tissues

Waves scattered by a weakly inhomogeneous random medium contain a predominant single scattering contribution as well as a multiple scattering contribution which is usually neglected, especially for imaging purposes. A method based on random matrix theory is proposed to separate the single and multiple scattering contributions. The experimental set up uses an array of sources/receivers placed in front of the medium. The impulse responses between every couple of transducers are measured and form a matrix. Single-scattering contributions are shown to exhibit a deterministic coherence along the antidiagonals of the array response matrix, whatever the distribution of inhomogeneities. This property is taken advantage of to discriminate single from multiple-scattered waves. This allows one to evaluate the absorption losses and the scattering losses separately, by comparing the multiple scattering intensity with a radiative transfer model. Moreover, the relative contribution of multiple scattering in the backscattered wave can be estimated, which serves as a validity test for the Born approximation. Experimental results are presented with ultrasonic waves in the MHz range, on a synthetic sample (agar-gelatine gel) as well as on breast tissues. Interestingly, the multiple scattering contribution is found to be far from negligible in the breast around 4.3 MHz.Comment: 35 pages, 11 figures, final version, contains the appendix of the original articl

EM Algorithm for Multiple Wideband Source Localization

A computationally efficient algorithm using the expectation-maximization (EM) algorithm for multiple wideband source localization in the near field of a sensor array/area is addressed in this thesis. Our idea is to decompose the observed sensor data, which is a superimposition of multiple sources, into the individual components in the frequency domain and then estimate the corresponding location parameters associated with each component separately. Instead of the conventional alternating projection (AP) method, we propose to adopt the EM algorithm in this work; our new method involves two steps, namely Expectation (E-step) and Maximization (M-step). In the E-step, the individual incident source waveforms are estimated. Then, in the M-step, the maximum likelihood estimates of the source location parameters are obtained. These two steps are executed iteratively and alternatively until the pre-defined convergence is reached. The computational complexity comparison between our proposed EM algorithm and the existing AP scheme is investigated. It is shown through Monte Carlo simulations that the computational complexity of the proposed EM algorithm is significantly lower than that of the existing AP algorithm

Proceedings of the second "international Traveling Workshop on Interactions between Sparse models and Technology" (iTWIST'14)

The implicit objective of the biennial "international - Traveling Workshop on Interactions between Sparse models and Technology" (iTWIST) is to foster collaboration between international scientific teams by disseminating ideas through both specific oral/poster presentations and free discussions. For its second edition, the iTWIST workshop took place in the medieval and picturesque town of Namur in Belgium, from Wednesday August 27th till Friday August 29th, 2014. The workshop was conveniently located in "The Arsenal" building within walking distance of both hotels and town center. iTWIST'14 has gathered about 70 international participants and has featured 9 invited talks, 10 oral presentations, and 14 posters on the following themes, all related to the theory, application and generalization of the "sparsity paradigm": Sparsity-driven data sensing and processing; Union of low dimensional subspaces; Beyond linear and convex inverse problem; Matrix/manifold/graph sensing/processing; Blind inverse problems and dictionary learning; Sparsity and computational neuroscience; Information theory, geometry and randomness; Complexity/accuracy tradeoffs in numerical methods; Sparsity? What's next?; Sparse machine learning and inference.Comment: 69 pages, 24 extended abstracts, iTWIST'14 website: http://sites.google.com/site/itwist1

The influence of random element displacement on DOA estimates obtained with (Khatri-Rao-)root-MUSIC

Although a wide range of direction of arrival (DOA) estimation algorithms has been described for a diverse range of array configurations, no specific stochastic analysis framework has been established to assess the probability density function of the error on DOA estimates due to random errors in the array geometry. Therefore, we propose a stochastic collocation method that relies on a generalized polynomial chaos expansion to connect the statistical distribution of random position errors to the resulting distribution of the DOA estimates. We apply this technique to the conventional root-MUSIC and the Khatri-Rao-root-MUSIC methods. According to Monte-Carlo simulations, this novel approach yields a speedup by a factor of more than 100 in terms of CPU-time for a one-dimensional case and by a factor of 56 for a two-dimensional case

Source localization and denoising: a perspective from the TDOA space

In this manuscript, we formulate the problem of denoising Time Differences of Arrival (TDOAs) in the TDOA space, i.e. the Euclidean space spanned by TDOA measurements. The method consists of pre-processing the TDOAs with the purpose of reducing the measurement noise. The complete set of TDOAs (i.e., TDOAs computed at all microphone pairs) is known to form a redundant set, which lies on a linear subspace in the TDOA space. Noise, however, prevents TDOAs from lying exactly on this subspace. We therefore show that TDOA denoising can be seen as a projection operation that suppresses the component of the noise that is orthogonal to that linear subspace. We then generalize the projection operator also to the cases where the set of TDOAs is incomplete. We analytically show that this operator improves the localization accuracy, and we further confirm that via simulation.Comment: 25 pages, 9 figure