Spatial priors for tomographic reconstructions from limited data

Tomografie is het reconstrueren van het inwendige van een object a.d.h.v externe metingen, b.v. beelden verkregen met X-stralen of microgolven. Deze thesis bekijkt de specifieke aspecten van microgolftomografie en magnetische resonantie beeldvorming (Magnetic Resonance Imaging – MRI); beide technieken zijn onschadelijk voor de mens. Terwijl het gebruik van MRI wijdverspreid is voor veel klinische toepassingen, is microgolftomografie nog niet in klinisch gebruik ondanks zijn potentiële voordelen. Door de lage kost en draagbaarheid van de toestellen is het een waardevolle aanvulling aan het assortiment

Compressed Sensing for Open-ended Waveguide Non-Destructive Testing and Evaluation

Ph. D. ThesisNon-destructive testing and evaluation (NDT&E) systems using open-ended waveguide (OEW) suffer from critical challenges. In the sensing stage, data acquisition is time-consuming by raster scan, which is difficult for on-line detection. Sensing stage also disregards demand for the latter feature extraction process, leading to an excessive amount of data and processing overhead for feature extraction. In the feature extraction stage, efficient and robust defect region segmentation in the obtained image is challenging for a complex image background. Compressed sensing (CS) demonstrates impressive data compression ability in various applications using sparse models. How to develop CS models in OEW NDT&E that jointly consider sensing & processing for fast data acquisition, data compression, efficient and robust feature extraction is remaining challenges. This thesis develops integrated sensing-processing CS models to address the drawbacks in OEW NDT systems and carries out their case studies in low-energy impact damage detection for carbon fibre reinforced plastics (CFRP) materials. The major contributions are: (1) For the challenge of fast data acquisition, an online CS model is developed to offer faster data acquisition and reduce data amount without any hardware modification. The images obtained with OEW are usually smooth which can be sparsely represented with discrete cosine transform (DCT) basis. Based on this information, a customised 0/1 Bernoulli matrix for CS measurement is designed for downsampling. The full data is reconstructed with orthogonal matching pursuit algorithm using the downsampling data, DCT basis, and the customised 0/1 Bernoulli matrix. It is hard to determine the sampling pixel numbers for sparse reconstruction when lacking training data, to address this issue, an accumulated sampling and recovery process is developed in this CS model. The defect region can be extracted with the proposed histogram threshold edge detection (HTED) algorithm after each recovery, which forms an online process. A case study in impact damage detection on CFRP materials is carried out for validation. The results show that the data acquisition time is reduced by one order of magnitude while maintaining equivalent image quality and defect region as raster scan. (2) For the challenge of efficient data compression that considers the later feature extraction, a feature-supervised CS data acquisition method is proposed and evaluated. It reserves interested features while reducing the data amount. The frequencies which reveal the feature only occupy a small part of the frequency band, this method finds these sparse frequency range firstly to supervise the later sampling process. Subsequently, based on joint sparsity of neighbour frame and the extracted frequency band, an aligned spatial-spectrum sampling scheme is proposed. The scheme only samples interested frequency range for required features by using a customised 0/1 Bernoulli measurement matrix. The interested spectral-spatial data are reconstructed jointly, which has much faster speed than frame-by-frame methods. The proposed feature-supervised CS data acquisition is implemented and compared with raster scan and the traditional CS reconstruction in impact damage detection on CFRP materials. The results show that the data amount is reduced greatly without compromising feature quality, and the gain in reconstruction speed is improved linearly with the number of measurements. (3) Based on the above CS-based data acquisition methods, CS models are developed to directly detect defect from CS data rather than using the reconstructed full spatial data. This method is robust to texture background and more time-efficient that HTED algorithm. Firstly, based on the histogram is invariant to down-sampling using the customised 0/1 Bernoulli measurement matrix, a qualitative method which only gives binary judgement of defect is developed. High probability of detection and accuracy is achieved compared to other methods. Secondly, a new greedy algorithm of sparse orthogonal matching pursuit (spOMP)-based defect region segmentation method is developed to quantitatively extract the defect region, because the conventional sparse reconstruction algorithms cannot properly use the sparse character of correlation between the measurement matrix and CS data. The proposed algorithms are faster and more robust to interference than other algorithms.China Scholarship Counci

Compensated Row-Column Ultrasound Imaging System

Ultrasound imaging is a valuable tool in many applications ranging from material science to medical imaging. While 2-D ultrasound imaging is more commonly used, 3-D ultrasound imaging offers unique opportunities that can only be found with the help of the extra dimension. Acquiring a 3-D ultrasound image can be done in two main ways: mechanically moving a transducer over a region of interest and using a fixed 2-D transducer. Mechanical motion introduces unwanted artifacts and increases image acquisition time, so a fixed 2-D is usually preferred. However, a fully addressed 2-D array will require a significant amount of connections and data to handle. This motivated the exploration of different simplification schemes to make 2-D arrays for 3-D ultrasound imaging feasible. A method that received a lot of attention for making real-time volumetric ultrasound imaging possible is the row-column method. The row-column method simplifies the fully addressed 2-D array by utilizing a set of 1-D arrays arranged in rows and another set in columns, one set will be responsible for transmit beamforming, while the other for receive beamforming. Using this setup, only $N+N$ connections are needed instead of $N\times N$. This simplification comes at the cost of image quality. Recent advances in row-column ultrasound imaging systems were largely focused on transducer design. However, these imaging systems face a few intrinsic challenges which cannot be addressed through transducer design alone: the issues of sparsity, speckle noise inherent to ultrasound, the spatially varying point spread function, and the ghosting artifacts inherent to the row-column method must all be taken into account. As such, strategies for tackling these intrinsic challenges in row-column imaging would be highly desired to improve imaging quality. In this thesis, we propose a novel compensated row-column ultrasound imaging system where the intrinsic characteristics of the transducer and other aspects of the physical row-column imaging apparatus are leveraged to computationally produce high quality ultrasound imagery. More specifically, the proposed system incorporates a novel conditional random field-driven computational image reconstruction component consisting of two phases: i) characterization and ii) compensation. In the characterization phase, a joint statistical image formation and noise model is introduced for characterizing the intrinsic properties of the physical row-column ultrasound imaging system. In the compensation phase, the developed joint image formation and noise model is incorporated alongside a conditional random field model within an energy minimization framework to reconstruct the compensated row-column ultrasound imagery. To explore the efficacy of the proposed concept, we introduced three different realizations of the proposed compensated row-column ultrasound imaging system. First, we introduce a compensated row-column imaging system based on a novel multilayered conditional random field driven framework to better account for local spatial relationships in the captured data. Second, we incorporated more global relationships by introducing a compensated row-column imaging system based around a novel edge-guided stochastically fully connected random field framework. Third, accounting for the case where the analytical image formation model may not optimally reflect the real-world physical system, we introduce a compensated row-column imaging system based around a data-driven spatially varying point-spread-function learning framework to better characterize the true physical image formation characteristics. While these different realizations of the compensated row-column system have their advantages and disadvantages, which will be discussed throughout this thesis, they all manage to boost the performance of the row-column method to comparable and often higher levels than the fully addressed 2-D array

Advanced Computational Methods for Oncological Image Analysis

[Cancer is the second most common cause of death worldwide and encompasses highly variable clinical and biological scenarios. Some of the current clinical challenges are (i) early diagnosis of the disease and (ii) precision medicine, which allows for treatments targeted to specific clinical cases. The ultimate goal is to optimize the clinical workflow by combining accurate diagnosis with the most suitable therapies. Toward this, large-scale machine learning research can define associations among clinical, imaging, and multi-omics studies, making it possible to provide reliable diagnostic and prognostic biomarkers for precision oncology. Such reliable computer-assisted methods (i.e., artificial intelligence) together with clinicians’ unique knowledge can be used to properly handle typical issues in evaluation/quantification procedures (i.e., operator dependence and time-consuming tasks). These technical advances can significantly improve result repeatability in disease diagnosis and guide toward appropriate cancer care. Indeed, the need to apply machine learning and computational intelligence techniques has steadily increased to effectively perform image processing operations—such as segmentation, co-registration, classification, and dimensionality reduction—and multi-omics data integration.

Earth resources: A continuing bibliography with indexes, issue 50

This bibliography lists 523 reports, articles and other documents introduced into the NASA scientific and technical information system between April 1 and June 30, 1986. Emphasis is placed on the use of remote sensing and geophysical instrumentation in spacecraft and aircraft to survey and inventory natural resources and urban areas. Subject matter is grouped according to agriculture and forestry, environmental changes and cultural resources, geodesy and cartography, geology and mineral resources, hydrology and water management, data processing and distribution systems, instrumentation and sensors, and economic analysis

A precise bare simulation approach to the minimization of some distances. Foundations

In information theory -- as well as in the adjacent fields of statistics, machine learning, artificial intelligence, signal processing and pattern recognition -- many flexibilizations of the omnipresent Kullback-Leibler information distance (relative entropy) and of the closely related Shannon entropy have become frequently used tools. To tackle corresponding constrained minimization (respectively maximization) problems by a newly developed dimension-free bare (pure) simulation method, is the main goal of this paper. Almost no assumptions (like convexity) on the set of constraints are needed, within our discrete setup of arbitrary dimension, and our method is precise (i.e., converges in the limit). As a side effect, we also derive an innovative way of constructing new useful distances/divergences. To illustrate the core of our approach, we present numerous examples. The potential for widespread applicability is indicated, too; in particular, we deliver many recent references for uses of the involved distances/divergences and entropies in various different research fields (which may also serve as an interdisciplinary interface)

Digital Holography Microscopy at Lab-on-a-Chip scale: novel algorithms and recording strategies

Il lavoro presentato è mirato allo sviluppo di nuove tecniche di microscopia olografica digitale (Digital Holography Microscopy, DHM), e di opportuni algoritmi numerici per lo studio di biomateriali in ambiente microfluidico. Nello specifico vengono affrontate due problematiche di imaging particolarmente rilevanti nello studio di sistemi Lab-on-a-Chip (LoC). Dapprima è stato studiato il problema della microscopia quantitativa di oggetti biologici osservati attraverso mezzi complessi, come soluzioni torbide e substrati diffondenti, dove la formazione dell’immagine è ostacolata da processi di scattering. Lo studio condotto è stato mirato all’analisi di processi di diffusione da layer statico e da mezzo liquido di tipo colloidale, in regime quasi-statico e dinamico. Sono stati sviluppati a tale scopo dei metodi di registrazione e nuovi algoritmi di ricostruzione dell’immagine olografica (Multi-Look Digital Holography, MLDH) che consentono di fornire un imaging quantitativo dei campioni in esame. Di particolare interesse è il caso di volumi di liquido costituiti da globuli rossi: nel lavoro presentato viene dimostrata la possibilità di studiare, mediante MLDH, processi di adesione cellulare di materiale biologico situato in presenza di flussi di globuli rossi ad alta concentrazione. La possibilità di visualizzare e analizzare quantitativamente materiale biologico all’interno di un capillare o una vena, compensando l’effetto di diffusione del sangue, potrebbe in futuro consentire di studiare la formazione all’interno del vaso di coaguli e placche di colesterolo, sintomatici dell’insorgere di malattie cardiache. La stessa tecnica è in grado di recuperare l’informazione distorta a causa della presenza all’interno del canale di ostacoli statici o quasi-statici (dovuti alla formazione di bio-film o sospensioni batteriche, o causata da processi di fabbricazione del canale microfluidico), aumentando così notevolmente la varietà dei processi biologici analizzabili su piattaforme LoC. Nel lavoro viene anche dimostrato come la presenza di un mezzo torbido possa essere sfruttata vantaggiosamente al fine di migliorare la qualità dell’immagine in sistemi di imaging basati su luce coerente. Parallelamente è stata messa a punto una tecnica interferometrica che, sfruttando il movimento dei campioni nei canali microfluidici, consente di sostituire un sensore convenzionale 2D con un sensore lineare, più compatto e integrabile a bordo del chip, e capace di fornire prestazioni superiori in termini di velocità di acquisizione. Il lavoro presentato descrive il processo di sintesi di un nuovo tipo di ologramma (Space-Time Digital Hologram, STDH), che consente di ottenere un Field-of-View (FoV) illimitato nella direzione del flusso e, quindi, di superare il trade-off esistente tra fattore di ingrandimento e FoV, comune ad ogni tecnica di microscopia convenzionale. Viene inoltre dimostrato che un STDH mantiene le caratteristiche e i vantaggi di un ologramma digitale standard, quali la focalizzazione numerica flessibile, che permette di analizzare contemporaneamente tutti gli oggetti presenti in un volume di liquido, e la possibilità di estrarre la segnatura di fase degli stessi

MS FT-2-2 7 Orthogonal polynomials and quadrature: Theory, computation, and applications

Quadrature rules find many applications in science and engineering. Their analysis is a classical area of applied mathematics and continues to attract considerable attention. This seminar brings together speakers with expertise in a large variety of quadrature rules. It is the aim of the seminar to provide an overview of recent developments in the analysis of quadrature rules. The computation of error estimates and novel applications also are described

Generalized averaged Gaussian quadrature and applications

A simple numerical method for constructing the optimal generalized averaged Gaussian quadrature formulas will be presented. These formulas exist in many cases in which real positive GaussKronrod formulas do not exist, and can be used as an adequate alternative in order to estimate the error of a Gaussian rule. We also investigate the conditions under which the optimal averaged Gaussian quadrature formulas and their truncated variants are internal