Non-linear Recovery of Sparse Signal Representations with Applications to Temporal and Spatial Localization

Foundations of signal processing are heavily based on Shannon's sampling theorem for acquisition, representation and reconstruction. This theorem states that signals should not contain frequency components higher than the Nyquist rate, which is half of the sampling rate. Then, the signal can be perfectly reconstructed from its samples. Increasing evidence shows that the requirements imposed by Shannon's sampling theorem are too conservative for many naturally-occurring signals, which can be accurately characterized by sparse representations that require lower sampling rates closer to the signal's intrinsic information rates. Finite rate of innovation (FRI) is a new theory that allows to extract underlying sparse signal representations while operating at a reduced sampling rate. The goal of this PhD work is to advance reconstruction techniques for sparse signal representations from both theoretical and practical points of view. Specifically, the FRI framework is extended to deal with applications that involve temporal and spatial localization of events, including inverse source problems from radiating fields. We propose a novel reconstruction method using a model-fitting approach that is based on minimizing the fitting error subject to an underlying annihilation system given by the Prony's method. First, we showed that this is related to the problem known as structured low-rank matrix approximation as in structured total least squares problem. Then, we proposed to solve our problem under three different constraints using the iterative quadratic maximum likelihood algorithm. Our analysis and simulation results indicate that the proposed algorithms improve the robustness of the results with respect to common FRI reconstruction schemes. We have further developed the model-fitting approach to analyze spontaneous brain activity as measured by functional magnetic resonance imaging (fMRI). For this, we considered the noisy fMRI time course for every voxel as a convolution between an underlying activity inducing signal (i.e., a stream of Diracs) and the hemodynamic response function (HRF). We then validated this method using experimental fMRI data acquired during an event-related study. The results showed for the first time evidence for the practical usage of FRI for fMRI data analysis. We also addressed the problem of retrieving a sparse source distribution from the boundary measurements of a radiating field. First, based on Green's theorem, we proposed a sensing principle that allows to relate the boundary measurements to the source distribution. We focused on characterizing these sensing functions with particular attention for those that can be derived from holomorphic functions as they allow to control spatial decay of the sensing functions. With this selection, we developed an FRI-inspired non-iterative reconstruction algorithm. Finally, we developed an extension to the sensing principle (termed eigensensing) where we choose the spatial eigenfunctions of the Laplace operator as the sensing functions. With this extension, we showed that eigensensing principle allows to extract partial Fourier measurements of the source functions from boundary measurements. We considered photoacoustic tomography as a potential application of these theoretical developments

Compressive phase retrieval

Thesis (Ph. D.)--Massachusetts Institute of Technology, Dept. of Mechanical Engineering, 2013.This electronic version was submitted by the student author. The certified thesis is available in the Institute Archives and Special Collections.Cataloged from student-submitted PDF version of thesis.Includes bibliographical references (p. 129-138).Recovering a full description of a wave from limited intensity measurements remains a central problem in optics. Optical waves oscillate too fast for detectors to measure anything but time{averaged intensities. This is unfortunate since the phase can reveal important information about the object. When the light is partially coherent, a complete description of the phase requires knowledge about the statistical correlations for each pair of points in space. Recovery of the correlation function is a much more challenging problem since the number of pairs grows much more rapidly than the number of points. In this thesis, quantitative phase imaging techniques that works for partially coherent illuminations are investigated. In order to recover the phase information with few measurements, the sparsity in each underly problem and ecient inversion methods are explored under the framework of compressed sensing. In each phase retrieval technique under study, diffraction during spatial propagation is exploited as an effective and convenient mechanism to uniformly distribute the information about the unknown signal into the measurement space. Holography is useful to record the scattered field from a sparse distribution of particles; the ability of localizing each particles using compressive reconstruction method is studied. When a thin sample is illuminated with partially coherent waves, the transport of intensity phase retrieval method is shown to be eective to recover the optical path length of the sample and remove the eect of the illumination. This technique is particularly suitable for X-ray phase imaging since it does not require a coherent source or any optical components. Compressive tomographic reconstruction, which makes full use of the priors that the sample consists of piecewise constant refractive indices, are demonstrated to make up missing data. The third technique, known as the phase space tomography (PST), addresses the correlation function recovery problem. Implementing the PST involves measuring many intensity images under spatial propagation. Experimental demonstration of a compressive reconstruction method, which finds the sparse solution by decomposing the correlation function into a few mutually uncorrelated coherent modes, is presented to produce accurate reconstruction even when the measurement suers from the 'missing cone' problem in the Fourier domain.by Lei Tian.Ph.D

Listening to Distances and Hearing Shapes:Inverse Problems in Room Acoustics and Beyond

A central theme of this thesis is using echoes to achieve useful, interesting, and sometimes surprising results. One should have no doubts about the echoes' constructive potential; it is, after all, demonstrated masterfully by Nature. Just think about the bat's intriguing ability to navigate in unknown spaces and hunt for insects by listening to echoes of its calls, or about similar (albeit less well-known) abilities of toothed whales, some birds, shrews, and ultimately people. We show that, perhaps contrary to conventional wisdom, multipath propagation resulting from echoes is our friend. When we think about it the right way, it reveals essential geometric information about the sources--channel--receivers system. The key idea is to think of echoes as being more than just delayed and attenuated peaks in 1D impulse responses; they are actually additional sources with their corresponding 3D locations. This transformation allows us to forget about the abstract \emph{room}, and to replace it by more familiar \emph{point sets}. We can then engage the powerful machinery of Euclidean distance geometry. A problem that always arises is that we do not know \emph{a priori} the matching between the peaks and the points in space, and solving the inverse problem is achieved by \emph{echo sorting}---a tool we developed for learning correct labelings of echoes. This has applications beyond acoustics, whenever one deals with waves and reflections, or more generally, time-of-flight measurements. Equipped with this perspective, we first address the ``Can one hear the shape of a room?'' question, and we answer it with a qualified ``yes''. Even a single impulse response uniquely describes a convex polyhedral room, whereas a more practical algorithm to reconstruct the room's geometry uses only first-order echoes and a few microphones. Next, we show how different problems of localization benefit from echoes. The first one is multiple indoor sound source localization. Assuming the room is known, we show that discretizing the Helmholtz equation yields a system of sparse reconstruction problems linked by the common sparsity pattern. By exploiting the full bandwidth of the sources, we show that it is possible to localize multiple unknown sound sources using only a single microphone. We then look at indoor localization with known pulses from the geometric echo perspective introduced previously. Echo sorting enables localization in non-convex rooms without a line-of-sight path, and localization with a single omni-directional sensor, which is impossible without echoes. A closely related problem is microphone position calibration; we show that echoes can help even without assuming that the room is known. Using echoes, we can localize arbitrary numbers of microphones at unknown locations in an unknown room using only one source at an unknown location---for example a finger snap---and get the room's geometry as a byproduct. Our study of source localization outgrew the initial form factor when we looked at source localization with spherical microphone arrays. Spherical signals appear well beyond spherical microphone arrays; for example, any signal defined on Earth's surface lives on a sphere. This resulted in the first slight departure from the main theme: We develop the theory and algorithms for sampling sparse signals on the sphere using finite rate-of-innovation principles and apply it to various signal processing problems on the sphere

Roadmap on signal processing for next generation measurement systems

Signal processing is a fundamental component of almost any sensor-enabled system, with a wide range of applications across different scientific disciplines. Time series data, images, and video sequences comprise representative forms of signals that can be enhanced and analysed for information extraction and quantification. The recent advances in artificial intelligence and machine learning are shifting the research attention towards intelligent, data-driven, signal processing. This roadmap presents a critical overview of the state-of-the-art methods and applications aiming to highlight future challenges and research opportunities towards next generation measurement systems. It covers a broad spectrum of topics ranging from basic to industrial research, organized in concise thematic sections that reflect the trends and the impacts of current and future developments per research field. Furthermore, it offers guidance to researchers and funding agencies in identifying new prospects.AerodynamicsMicrowave Sensing, Signals & System

Microwave Sensing and Imaging

In recent years, microwave sensing and imaging have acquired an ever-growing importance in several applicative fields, such as non-destructive evaluations in industry and civil engineering, subsurface prospection, security, and biomedical imaging. Indeed, microwave techniques allow, in principle, for information to be obtained directly regarding the physical parameters of the inspected targets (dielectric properties, shape, etc.) by using safe electromagnetic radiations and cost-effective systems. Consequently, a great deal of research activity has recently been devoted to the development of efficient/reliable measurement systems, which are effective data processing algorithms that can be used to solve the underlying electromagnetic inverse scattering problem, and efficient forward solvers to model electromagnetic interactions. Within this framework, this Special Issue aims to provide some insights into recent microwave sensing and imaging systems and techniques

Online monitoring instantaneous 2D temperature distributions in a furnace using acoustic tomography based on frequency division multiplexing

The online and accurate capture of dynamic changes in furnace temperature distribution is crucial for production efficiency improvement and international environmental policy compliance in power plants. To achieve this, a measurement system with a reliable online reconstruction capability and high temporal resolution is necessary. This paper presents a novel technique that can improve the temporal resolution of the currently existing acoustic tomography (AT) system using frequency division multiplexing (FDM). This method allows for concurrent transmissions of acoustic signals in several different frequency bands instead of a sequential manner, which leads to more efficient channel utilization and allows all acoustic signals to be acquired at the same time, so that a better temporal uniformity of multipath acoustic signals can be realized. Theoretical analysis and experiments have been conducted to verify the effectiveness of this technique. The results prove that the proposed method can significantly improve the temporal resolution of the AT system while maintaining the accuracy and robustness of the reconstruction

Diffuse Optical Imaging with Ultrasound Priors and Deep Learning

Diffuse Optical Imaging (DOI) techniques are an ever growing field of research as they are noninvasive, compact, cost-effective and can furnish functional information about human tissues. Among others, they include techniques such as Tomography, which solves an inverse reconstruction problem in a tissue volume, and Mapping which only seeks to find values on a tissue surface. Limitations in reliability and resolution, due to the ill-posedness of the underlying inverse problems, have hindered the clinical uptake of this medical imaging modality. Multimodal imaging and Deep Learning present themselves as two promising solutions to further research in DOI. In relation to the first idea, we implement and assess here a set of methods for SOLUS, a combined Ultrasound (US) and Diffuse Optical Tomography (DOT) probe for breast cancer diagnosis. An ad hoc morphological prior is extracted from US B-mode images and utilised for the regularisation of the inverse problem in DOT. Combination of the latter in reconstruction with a linearised forward model for DOT is assessed on specifically designed dual phantoms. The same reconstruction approach with the incorporation of a spectral model has been assessed on meat phantoms for reconstruction of functional properties. A simulation study with realistic digital phantoms is presented for an assessment of a non-linear model in reconstruction for the quantification of optical properties of breast lesions. A set of machine learning tools is presented for diagnosis breast lesions based on the reconstructed optical properties. A preliminary clinical study with the SOLUS probe is presented. Finally, a specifically designed deep learning architecture for diffusion is applied to mapping on the brain cortex or Diffuse Optical Cortical Mapping (DOCM). An assessment of its performances is presented on simulated and experimental data