Low-rank sparse matrix decomposition for sparsity-driven SAR image reconstruction

We consider the development of a synthetic aperture radar (SAR) image reconstruction method that decomposes the imaged field into a sparse and a low-rank component. Such a decomposition is of interest in image analysis tasks such as segmentation and background subtraction. Conventionally, such operations are performed after SAR image formation. However image formation methods may produce images that are not well suited for such interpretation tasks since they do not incorporate interpretation objectives to the SAR imaging problem. We exploit recent work on sparse and low-rank decomposition of matrices and incorporate such a decomposition into the process of SAR image formation. The outcome is a method that jointly reconstructs a SAR image and decomposes the formed image into its low-rank background and spatially sparse components. We demonstrate the effectiveness of the proposed method on both synthetic and real SAR images

Sparse and Redundant Representations for Inverse Problems and Recognition

Sparse and redundant representation of data enables the description of signals as linear combinations of a few atoms from a dictionary. In this dissertation, we study applications of sparse and redundant representations in inverse problems and object recognition. Furthermore, we propose two novel imaging modalities based on the recently introduced theory of Compressed Sensing (CS). This dissertation consists of four major parts. In the first part of the dissertation, we study a new type of deconvolution algorithm that is based on estimating the image from a shearlet decomposition. Shearlets provide a multi-directional and multi-scale decomposition that has been mathematically shown to represent distributed discontinuities such as edges better than traditional wavelets. We develop a deconvolution algorithm that allows for the approximation inversion operator to be controlled on a multi-scale and multi-directional basis. Furthermore, we develop a method for the automatic determination of the threshold values for the noise shrinkage for each scale and direction without explicit knowledge of the noise variance using a generalized cross validation method. In the second part of the dissertation, we study a reconstruction method that recovers highly undersampled images assumed to have a sparse representation in a gradient domain by using partial measurement samples that are collected in the Fourier domain. Our method makes use of a robust generalized Poisson solver that greatly aids in achieving a significantly improved performance over similar proposed methods. We will demonstrate by experiments that this new technique is more flexible to work with either random or restricted sampling scenarios better than its competitors. In the third part of the dissertation, we introduce a novel Synthetic Aperture Radar (SAR) imaging modality which can provide a high resolution map of the spatial distribution of targets and terrain using a significantly reduced number of needed transmitted and/or received electromagnetic waveforms. We demonstrate that this new imaging scheme, requires no new hardware components and allows the aperture to be compressed. Also, it presents many new applications and advantages which include strong resistance to countermesasures and interception, imaging much wider swaths and reduced on-board storage requirements. The last part of the dissertation deals with object recognition based on learning dictionaries for simultaneous sparse signal approximations and feature extraction. A dictionary is learned for each object class based on given training examples which minimize the representation error with a sparseness constraint. A novel test image is then projected onto the span of the atoms in each learned dictionary. The residual vectors along with the coefficients are then used for recognition. Applications to illumination robust face recognition and automatic target recognition are presented

Compressed sensing for enhanced through-the-wall radar imaging

Through-the-wall radar imaging (TWRI) is an emerging technology that aims to capture scenes behind walls and other visually opaque materials. The abilities to sense through walls are highly desirable for both military and civil applications, such as search and rescue missions, surveillance, and reconnaissance. TWRI systems, however, face with several challenges including prolonged data acquisition, large objects, strong wall clutter, and shadowing effects, which limit the radar imaging performances and hinder target detection and localization

SAR imaging of moving targets by subaperture based low-rank and sparse decomposition

Synthetic aperture radar (SAR) has gained significance as an indispensable instrument of remote sensing and airborne surveillance. Its applications extend to 3D terrain mapping, oil spill detection, crop yield estimation and disaster evaluation. SAR utilizes platform motion to synthesize a large antenna thus rendering a very fine spatial resolution. Nevertheless, imaging of moving targets with SAR is a challenging problem. In this thesis, we propose a moving target imaging approach for SAR which exploits the low-rank and sparse decomposition (LRSD) of the subaperture data. As a first step, multiple subapertures are constructed from the raw data using frequency domain filtering. In contrast to the stationary points, moving targets in the SAR scene shift their position in the various subapertures. This enables a successful low-rank and sparse decomposition of the subaperture data where the sparse component captures the moving targets’ phase histories and reflectivity profiles. On the other hand, the low-rank component consists of the static background due to fewer spatial variations in multiple subapertures. This framework allows the reconstruction of full-resolution sparse and low-rank images by combining the spectral information of the decomposed subapertures. Furthermore, it enhances the applicability of sparsity-driven moving target imaging frameworks to very low signal to clutter ratio (SCR) scenarios by offering a considerable SCR performance improvement. We manifest the effectiveness of our approach through experiments with synthetic as well as real SAR data. Our real SAR experiments were based on MiniSAR and EMISAR data

Dictionary learning and low-rank sparse matrix decomposition for sparsity driven SAR image reconstruction

Synthetic aperture radar (SAR) is one of the most widely used remote sensing modalities, providing images for a variety of applications including those in defense, environmental science, and weather forecasting. However, conventionally formed SAR imagery from undersampled observed data, arising in several emerging applications and sensing scenarios, suffers from artifacts that might limit effective use of such imagery in remote sensing applications. Recently, sparsity-driven SAR imaging has emerged as an e ective framework to alleviate such problems. Sparsity-based methods for SAR imaging have employed overcomplete dictionaries to represent the magnitude of the complex-valued eld sparsely. Selection of an appropriate dictionary with respect to the features of the particular type of underlying scene plays an important role in these methods. In this thesis, we develop two new sparsity-driven SAR imaging methods that significantly expand the domain of applicability of sparsity-based methods in SAR imaging. Our first contribution involves the development of a new reconstruction method that is based on learning sparsifying dictionaries and using such learned dictionaries in the reconstruction process. Adaptive dictionaries learned from data can represent the magnitude of complex-valued field more effectively and hence have the potential to widen the applicability of sparsity-based radar imaging. Our framework allows the use of both adaptive dictionaries learned offline from a training set and those learned online from the undersampled data used in image formation. We demonstrate the effectiveness of the proposed dictionary learning-based SAR imaging approach as well as the improvements it provides, on both synthetic and real data. The second contribution of this thesis involves the development of a reconstruction method that decomposes the imaged field into a sparse and a low-rank component. Such a decomposition is of interest in image analysis tasks such as segmentation and background subtraction. Conventionally, such operations are performed after SAR image formation. We exploit recent work on sparse and low-rank decomposition of matrices and incorporate such a decomposition into the process of SAR image formation. The outcome is a method that jointly reconstructs a SAR image and decomposes the formed image into its low-rank background and spatially sparse components. We demonstrate the effectiveness of the proposed method on both synthetic and real SAR images

Efficient algorithms and data structures for compressive sensing

Wegen der kontinuierlich anwachsenden Anzahl von Sensoren, und den stetig wachsenden Datenmengen, die jene produzieren, stößt die konventielle Art Signale zu verarbeiten, beruhend auf dem Nyquist-Kriterium, auf immer mehr Hindernisse und Probleme. Die kürzlich entwickelte Theorie des Compressive Sensing (CS) formuliert das Versprechen einige dieser Hindernisse zu beseitigen, indem hier allgemeinere Signalaufnahme und -rekonstruktionsverfahren zum Einsatz kommen können. Dies erlaubt, dass hierbei einzelne Abtastwerte komplexer strukturierte Informationen über das Signal enthalten können als dies bei konventiellem Nyquistsampling der Fall ist. Gleichzeitig verändert sich die Signalrekonstruktion notwendigerweise zu einem nicht-linearen Vorgang und ebenso müssen viele Hardwarekonzepte für praktische Anwendungen neu überdacht werden. Das heißt, dass man zwischen der Menge an Information, die man über Signale gewinnen kann, und dem Aufwand für das Design und Betreiben eines Signalverarbeitungssystems abwägen kann und muss. Die hier vorgestellte Arbeit trägt dazu bei, dass bei diesem Abwägen CS mehr begünstigt werden kann, indem neue Resultate vorgestellt werden, die es erlauben, dass CS einfacher in der Praxis Anwendung finden kann, wobei die zu erwartende Leistungsfähigkeit des Systems theoretisch fundiert ist. Beispielsweise spielt das Konzept der Sparsity eine zentrale Rolle, weshalb diese Arbeit eine Methode präsentiert, womit der Grad der Sparsity eines Vektors mittels einer einzelnen Beobachtung geschätzt werden kann. Wir zeigen auf, dass dieser Ansatz für Sparsity Order Estimation zu einem niedrigeren Rekonstruktionsfehler führt, wenn man diesen mit einer Rekonstruktion vergleicht, welcher die Sparsity des Vektors unbekannt ist. Um die Modellierung von Signalen und deren Rekonstruktion effizienter zu gestalten, stellen wir das Konzept von der matrixfreien Darstellung linearer Operatoren vor. Für die einfachere Anwendung dieser Darstellung präsentieren wir eine freie Softwarearchitektur und demonstrieren deren Vorzüge, wenn sie für die Rekonstruktion in einem CS-System genutzt wird. Konkret wird der Nutzen dieser Bibliothek, einerseits für das Ermitteln von Defektpositionen in Prüfkörpern mittels Ultraschall, und andererseits für das Schätzen von Streuern in einem Funkkanal aus Ultrabreitbanddaten, demonstriert. Darüber hinaus stellen wir für die Verarbeitung der Ultraschalldaten eine Rekonstruktionspipeline vor, welche Daten verarbeitet, die im Frequenzbereich Unterabtastung erfahren haben. Wir beschreiben effiziente Algorithmen, die bei der Modellierung und der Rekonstruktion zum Einsatz kommen und wir leiten asymptotische Resultate für die benötigte Anzahl von Messwerten, sowie die zu erwartenden Lokalisierungsgenauigkeiten der Defekte her. Wir zeigen auf, dass das vorgestellte System starke Kompression zulässt, ohne die Bildgebung und Defektlokalisierung maßgeblich zu beeinträchtigen. Für die Lokalisierung von Streuern mittels Ultrabreitbandradaren stellen wir ein CS-System vor, welches auf einem Random Demodulators basiert. Im Vergleich zu existierenden Messverfahren ist die hieraus resultierende Schätzung der Kanalimpulsantwort robuster gegen die Effekte von zeitvarianten Funkkanälen. Um den inhärenten Modellfehler, den gitterbasiertes CS begehen muss, zu beseitigen, zeigen wir auf wie Atomic Norm Minimierung es erlaubt ohne die Einschränkung auf ein endliches und diskretes Gitter R-dimensionale spektrale Komponenten aus komprimierten Beobachtungen zu schätzen. Hierzu leiten wir eine R-dimensionale Variante des ADMM her, welcher dazu in der Lage ist die Signalkovarianz in diesem allgemeinen Szenario zu schätzen. Weiterhin zeigen wir, wie dieser Ansatz zur Richtungsschätzung mit realistischen Antennenarraygeometrien genutzt werden kann. In diesem Zusammenhang präsentieren wir auch eine Methode, welche mittels Stochastic gradient descent Messmatrizen ermitteln kann, die sich gut für Parameterschätzung eignen. Die hieraus resultierenden Kompressionsverfahren haben die Eigenschaft, dass die Schätzgenauigkeit über den gesamten Parameterraum ein möglichst uniformes Verhalten zeigt. Zuletzt zeigen wir auf, dass die Kombination des ADMM und des Stochastic Gradient descent das Design eines CS-Systems ermöglicht, welches in diesem gitterfreien Szenario wünschenswerte Eigenschaften hat.Along with the ever increasing number of sensors, which are also generating rapidly growing amounts of data, the traditional paradigm of sampling adhering the Nyquist criterion is facing an equally increasing number of obstacles. The rather recent theory of Compressive Sensing (CS) promises to alleviate some of these drawbacks by proposing to generalize the sampling and reconstruction schemes such that the acquired samples can contain more complex information about the signal than Nyquist samples. The proposed measurement process is more complex and the reconstruction algorithms necessarily need to be nonlinear. Additionally, the hardware design process needs to be revisited as well in order to account for this new acquisition scheme. Hence, one can identify a trade-off between information that is contained in individual samples of a signal and effort during development and operation of the sensing system. This thesis addresses the necessary steps to shift the mentioned trade-off more to the favor of CS. We do so by providing new results that make CS easier to deploy in practice while also maintaining the performance indicated by theoretical results. The sparsity order of a signal plays a central role in any CS system. Hence, we present a method to estimate this crucial quantity prior to recovery from a single snapshot. As we show, this proposed Sparsity Order Estimation method allows to improve the reconstruction error compared to an unguided reconstruction. During the development of the theory we notice that the matrix-free view on the involved linear mappings offers a lot of possibilities to render the reconstruction and modeling stage much more efficient. Hence, we present an open source software architecture to construct these matrix-free representations and showcase its ease of use and performance when used for sparse recovery to detect defects from ultrasound data as well as estimating scatterers in a radio channel using ultra-wideband impulse responses. For the former of these two applications, we present a complete reconstruction pipeline when the ultrasound data is compressed by means of sub-sampling in the frequency domain. Here, we present the algorithms for the forward model, the reconstruction stage and we give asymptotic bounds for the number of measurements and the expected reconstruction error. We show that our proposed system allows significant compression levels without substantially deteriorating the imaging quality. For the second application, we develop a sampling scheme to acquire the channel Impulse Response (IR) based on a Random Demodulator that allows to capture enough information in the recorded samples to reliably estimate the IR when exploiting sparsity. Compared to the state of the art, this in turn allows to improve the robustness to the effects of time-variant radar channels while also outperforming state of the art methods based on Nyquist sampling in terms of reconstruction error. In order to circumvent the inherent model mismatch of early grid-based compressive sensing theory, we make use of the Atomic Norm Minimization framework and show how it can be used for the estimation of the signal covariance with R-dimensional parameters from multiple compressive snapshots. To this end, we derive a variant of the ADMM that can estimate this covariance in a very general setting and we show how to use this for direction finding with realistic antenna geometries. In this context we also present a method based on a Stochastic gradient descent iteration scheme to find compression schemes that are well suited for parameter estimation, since the resulting sub-sampling has a uniform effect on the whole parameter space. Finally, we show numerically that the combination of these two approaches yields a well performing grid-free CS pipeline

PLACE: Post LANDSAT D Advanced Concept Evaluation

There are no author-identified significant results in this report