# Agenda

## PhD Thesis Defence

- Thursday, 10 January 2019
- 15:00-17:00
- Aula Senaatszaal

**Aleksandar Jovic**

## PhD Thesis Defence

- Friday, 21 December 2018
- 12:00-14:30
- Aula Senaatszaal

### Multi-Microphone Noise Reduction for Hearing Assistive Devices

**Andreas Koutrouvelis**

The paramount importance of good hearing in everyday life has driven an exploration into the improvement of hearing capabilities of (hearing impaired) people in acoustic challenging situations using hearing assistive devices (HADs). HADs are small portable devices, which primarily aim at improving the intelligibility of an acoustic source that has drawn the attention of the HAD user. One of the most important steps to achieve this is via filtering the sound recorded using the HAD microphones, such that ideally all unwanted acoustic sources in the acoustic scene are suppressed, while the target source is maintained undistorted. Modern HAD systems often consist of two collaborative (typically wirelessly connected) HADs, each placed on a different ear. These HAD systems are commonly referred to as binaural HAD systems. The noise reduction filters designed for binaural HAD systems are referred to as binaural beamformers.

Binaural beamformers typically change the magnitude and phase relations of the microphone signals by forming a beam towards the target's direction while ideally suppressing all other directions. This may alter the spatial impression of the acoustic scene, as the filtered sources now reach both ears with possibly different relative phase and magnitude differences compared to before processing. This will appear unnatural to the HAD user. Therefore, there is an increasing interest in the preservation of the spatial information (also referred to as binaural cues) of the acoustic scene after processing. The present dissertation is mainly concerned with this particular problem and proposes several alternative binaural beamformers which try to exploit the available degrees of freedom to achieve optimal performance in both noise reduction and binaural-cue preservation.

Additional information ...

## PhD Thesis Defence

- Friday, 7 December 2018
- 12:30-13:30
- Aula Senaatszaal

### Image formation for future radio telescopes

**Shahrzad Naghibzadeh**

Fundamental scientific questions such as how the first stars were formed or how the universe came into existence and evolved to its present state drive us to observe weak radio signals impinging on the earth from the early days of the universe. During the last century, radio astronomy has been vastly advancing. Important discoveries on the formation of various celestial objects such as pulsars, neutron stars, black holes, radio galaxies and quasars are the result of radio astronomical observations. To study celestial objects and the astrophysical processes that are responsible for their radio emissions, images must be formed. This is done with the help of large radio telescope arrays.

Next generation radio telescopes such as the Low Frequency Array Radio Telescope (LOFAR) and the Square Kilometer Array (SKA), bring about increasingly more observational evidence for the study of the radio sky by generating very high resolution and high fidelity images. In this dissertation, we study radio astronomical imaging as the problem of estimating the sky spatial intensity distribution over the field of view of the radio telescope array from the incomplete and noisy array data. The increased sensitivity, resolution and sky coverage of the new instruments pose additional challenges to the current radio astronomical imaging pipeline. Namely, the large amount of data captured by the radio telescopes cannot be stored and needs to be processed quasi-realtime.

Many pixel-based imaging algorithms, such as the widely-used CLEAN [3] algorithm, are not scalable to the size of the required images and perform very slow in high resolution scenarios. Therefore, there is an urgent need for new efficient imaging algorithms. Moreover, regardless of the amount of collected data, there is an inherent loss of information in themeasurement process due to physical limitations. Therefore, to recover physically meaningful images additional information in the form of constraints and regularizing assumptions are necessary. The central objective of the current dissertation is to introduce advanced algebraic techniques together with custom-made regularization schemes to speed up the image formation pipeline of the next generation radio telescopes.

Signal processing provides powerful tools to address these issues. In the current work, following a signal processing model of the radio astronomical observation process, we first analyze the imaging system based on tools from numerical linear algebra, sampling, interpolation and filtering theory to investigate the inherent loss of information in the measurement process. Based on these results, we show that the imaging problem in radio astronomy is highly ill-posed and regularization is necessary to find a stable and physically meaningful image. We continue by deriving an adequate model for the imagingproblem in radio interferometry in the context of statistical estimation theory. Moreover, we introduce a framework to incorporate regularization assumptions into the measurement model by borrowing the concept of preconditioning from numerical linear algebra.

Radio emissions observed by radio telescopes appear either as distributed radiation from diffuse media or as compact emission from isolated point-like sources. Based on this observation, different source models need to be applied in the imaging problem formulation to obtain the best reconstruction performance. Due to the ill-posedness of the imaging problem in radio astronomy, to guarantee a reliable image reconstruction, propermodeling of the source emissions and regularizing assumptions are of utmost importance. We integrate these assumptions implementing a multi-basis dictionary based on the proposed preconditioning formalism.

In traditional radio astronomical imaging methods, the constraints and priormodels, such as positivity and sparsity, are employed for the complete image. However, large radio sky images usually manifest individual source occupancy regions in a large empty background. Based on this observation, we propose to split the field of view into multiple regions of source occupancy. Leveraging a stochastic primal dual algorithm we apply adequate regularization on each facet. We demonstrate the merits of applying facet-based regularization in terms of memory savings and computation time by realistic simulations.

The formulation of the radio astronomical imaging problem has a direct consequence on the radio sky estimation performance. We define the astronomical imaging problem in a Bayesian-inspired regularized maximum likelihood formulation. Based on this formalism, we develop a general algorithmic framework that can handle diffuse as well as compact source models. Leveraging the linearity of radio astronomical imaging problem, we propose to directly embody the regularization operator into the system by right preconditioning. We employ an iterative method based on projections onto Krylov subspaces to solve the subsequent system. The proposed algorithmis named PRIor-conditioned Fast Iterative Radio Astronomy (PRIFIRA). We motivate the use of a beamformed image as an efficient regularizing prior-conditioner for diffuse emission recovery. Different sparsity-based regularization priors are incorporated in the algorithmic framework by generalizing the core algorithm with iterative re-weighting schemes.

We evaluate the performance of PRIFIRA based on simulated measurements as well as astronomical data and compare to the state-of-the-art imaging methods. We conclude that the proposed method maintains competitive reconstruction quality with the current techniques while remaining flexible in terms of different regularization schemes. Moreover, we show that the imaging efficiency can be greatly improved by exploiting Krylov subspace methods together with an appropriate noise-based stopping criteria.

Based on the results from this thesis we can conclude that with the help of advanced techniques from signal processing and numerical linear algebra, customized algorithms can be designed to tackle some of the challenges in the next generation radio telescope imaging. We note that since radio interferometric imaging can be considered as an instance of the broad area of inverse imaging problems, the numerical techniques as well as regularization methods developed in this dissertation have a direct impact on many different imaging application areas, such as biomedical and geophysics/seismic imaging.

Additional information ...

## Signal Processing Seminar

- Friday, 7 December 2018
- 10:00-10:45
- HB 17.150

### Image Reconstruction Using Training Images

**Per Christian Hansen**

*Technical University of Denmark*

Priors are essential for computing stable solutions to ill-posed problems, and they take many different forms. Here we consider priors in the form of cross-section images of the object, and this information must be used in a fast, reliable, and computationally efficient manner. We describe an algorithmic framework for this: From a set of training images we use techniques from machine learning to form a dictionary that captures the desired features, and we then compute a reconstruction with a sparse representation in this dictionary. We describe how to stably compute the dictionary through a regularized non-negative matrix factorization, and we study how this dictionary affects the reconstruction quality. Simulations show that for textural images our approach is superior to other methods used for limited-data problems.

### About the speaker

Professor Per Christian Hansen has worked with numerical regularization algorithms for 30 years, and he has published 4 books and 100+ papers in leading journals. He has developed a number of software packages, of which Regularization Tools (now in its 4th version) is a popular toolbox for analysis and solution of discrete inverse problems. His current research projects involve algorithms for tomographic reconstruction and iterative image deblurring algorithms. He is a SIAM fellow in recognition of his work on inverse problems and regularization.Additional information ...

## PhD Thesis Defence

- Monday, 26 November 2018
- 12:30-14:00
- Aula Senaatszaal

### Surface Acoustic Mode Aluminium Nitride Transducer for micro-size liquid sensing applications

**Thu Hang Bui**

## MSc BME Thesis Presentation

- Friday, 16 November 2018
- 14:30-15:15
- 3ME room J

### The effect of dopamine release on electrical neural activity in the prefrontal cortex

**Jack Tchimino**

How can certain oscillations be detected from the measured brain signals?

Additional information ...

## PhD Thesis Defence

- Wednesday, 24 October 2018
- 12:30-14:00
- Aula Senaatszaal

### Free standing interconnects for stretchable electronics

**Shivani Joshi**

## MSc ME Thesis Presentation

- Tuesday, 23 October 2018
- 15:00-16:00
- EKL Colloquiumroom

### Fabrication and reliability study of parylene-ceramic based flexible interconnects for implantable devices

**Diane Wu**

## Signal Processing Seminar

- Tuesday, 16 October 2018
- 13:30-16:30
- EWI LB 1.010 Snijderszaal

### Signal Processing Mini-Symposium

**Hagit Messer, KVS Hari, Andrea Simonetto**

### Talk 1: Capitalizing on the Cellular Technology Opportunities and Challenges for Near Ground Weather Monitoring

**Prof. Hagit Messer**

School of Electrical Engineering, Tel Aviv University, Israel

### Talk 2: Spatial Modulation Techniques in Wireless Systems

**Prof. K.V.S. Hari**

Dept. of ECE, Indian Institute of Science, Bangalore

### Talk 3: Time-varying optimization: algorithms and engineering applications

**Dr. Andrea Simonetto**

IBM Research Ireland, Dublin, IrelandAdditional information ...

## Signal Processing Seminar

- Thursday, 4 October 2018
- 13:30-14:30
- HB 17.150

### Machine learning in physical sciences

**Peter Gerstoft**

*UC San Diego*

Machine learning (ML) is booming thanks to efforts promoted by Google. However, ML also has use in physical sciences. I start with a general overview of ML for supervised/unsupervised learning. Then I will focus on my applications of ML in array processing in seismology and ocean acoustics. This will include source localization using neural networks or graph processing. Final example is using ML-based tomography to obtain high-resolution subsurface geophysical structure in Long Beach, CA, from seismic noise recorded on a 5200-element array. This method exploits the dense sampling obtained by ambient noise processing on large arrays by learning a dictionary of local, or small-scale, geophysical features directly from the data.

Additional information ...