Signal processing for microwave imaging systems with very sparse array

This dissertation investigates image reconstruction algorithms for near-field, two dimensional (2D) synthetic aperture radar (SAR) using compressed sensing (CS) based methods. In conventional SAR imaging systems, acquiring higher-quality images requires longer measuring time and/or more elements in an antenna array. Millimeter wave imaging systems using evenly-spaced antenna arrays also have spatial resolution constraints due to the large size of the antennas. This dissertation applies the CS principle to a bistatic antenna array that consists of separate transmitter and receiver subarrays very sparsely and non-uniformly distributed on a 2D plane. One pair of transmitter and receiver elements is turned on at a time, and different pairs are turned on in series to achieve synthetic aperture and controlled random measurements. This dissertation contributes to CS-hardware co-design by proposing several signal-processing methods, including monostatic approximation, re-gridding, adaptive interpolation, CS-based reconstruction, and image denoising. The proposed algorithms enable the successful implementation of CS-SAR hardware cameras, improve the resolution and image quality, and reduce hardware cost and experiment time. This dissertation also describes and analyzes the results for each independent method. The algorithms proposed in this dissertation break the limitations of hardware configuration. By using 16 x 16 transmit and receive elements with an average space of 16 mm, the sparse-array camera achieves the image resolution of 2 mm. This is equivalent to six percent of the λ/4 evenly-spaced array. The reconstructed images achieve similar quality as the fully-sampled array with the structure similarity (SSIM) larger than 0.8 and peak signal-to-noise ratio (PSNR) greater than 25 --Abstract, page iv

Calibrationless Multi-coil Magnetic Resonance Imaging with Compressed Sensing

We present a method for combining the data retrieved by multiple coils of a Magnetic Resonance Imaging (MRI) system with the a priori assumption of compressed sensing to reconstruct a single image. The final image is the result of an optimization problem that only includes constraints based on fundamental physics (Maxwell's equations and the Biot-Savart law) and accepted phenomena (e.g. sparsity in the Wavelet domain). The problem is solved using an alternating minimization approach: two convex optimization problems are alternately solved, one with the Fast Iterative Shrinkage Threshold Algorithm (FISTA) and the other with the Primal-Dual Hybrid Gradient (PDHG) method. We show results on simulated data as well as data of the knee, brain, and ankle. In all cases studied, results from the new algorithm show higher quality and increased detail when compared to conventional reconstruction algorithms

Compressive sensing for 3D microwave imaging systems

Compressed sensing (CS) image reconstruction techniques are developed and experimentally implemented for wideband microwave synthetic aperture radar (SAR) imaging systems with applications to nondestructive testing and evaluation. These techniques significantly reduce the number of spatial measurement points and, consequently, the acquisition time by sampling at a level lower than the Nyquist-Shannon rate. Benefiting from a reduced number of samples, this work successfully implemented two scanning procedures: the nonuniform raster and the optimum path. Three CS reconstruction approaches are also proposed for the wideband microwave SAR-based imaging systems. The first approach reconstructs a full-set of raw data from undersampled measurements via L1-norm optimization and consequently applies 3D forward SAR on the reconstructed raw data. The second proposed approach employs forward SAR and reverse SAR (R-SAR) transforms in each L1-norm optimization iteration reconstructing images directly. This dissertation proposes a simple, elegant truncation repair method to combat the truncation error which is a critical obstacle to the convergence of the CS iterative algorithm. The third proposed CS reconstruction algorithm is the adaptive basis selection (ABS) compressed sensing. Rather than a fixed sparsifying basis, the proposed ABS method adaptively selects the best basis from a set of bases in each iteration of the L1-norm optimization according to a proposed decision metric that is derived from the sparsity of the image and the coherence between the measurement and sparsifying matrices. The results of several experiments indicate that the proposed algorithms recover 2D and 3D SAR images with only 20% of the spatial points and reduce the acquisition time by up to 66% of that of conventional methods while maintaining or improving the quality of the SAR images --Abstract, page iv

Under-Sampled Reconstruction Techniques for Accelerated Magnetic Resonance Imaging

Due to physical and biological constraints and requirements on the minimum resolution and SNR, the acquisition time is relatively long in magnetic resonance imaging (MRI). Consequently, a limited number of pulse sequences can be run in a clinical MRI session because of constraints on the total acquisition time due to patient comfort and cost considerations. Therefore, it is strongly desired to reduce the acquisition time without compromising the reconstruction quality. This thesis concerns under-sampled reconstruction techniques for acceleration of MRI acquisitions, i.e., parallel imaging and compressed sensing. While compressed sensing MRI reconstructions are commonly regularized by penalizing the decimated wavelet transform coefficients, it is shown in this thesis that the visual artifacts, associated with the lack of translation-invariance of the wavelet basis in the decimated form, can be avoided by penalizing the undecimated wavelet transform coefficients, i.e., the stationary wavelet transform (SWT). An iterative SWT thresholding algorithm for combined SWT-regularized compressed sensing and parallel imaging reconstruction is presented. Additionally, it is shown that in MRI applications involving multiple sequential acquisitions, e.g., quantitative T1/T2 mapping, the correlation between the successive acquisitions can be incorporated as an additional constraint for joint under-sampled reconstruction, resulting in improved reconstruction performance. While quantitative measures of quality, e.g., reconstruction error with respect to the fully-sampled reference, are commonly used for performance evaluation and comparison of under-sampled reconstructions, this thesis shows that such quantitative measures do not necessarily correlate with the subjective quality of reconstruction as perceived by radiologists and other expert end users. Therefore, unless accompanied by subjective evaluations, quantitative quality measurements/comparisons will be of limited clinical impact. The results of experiments aimed at subjective evaluation/comparison of different under-sampled reconstructions for specific clinical neuroimaging MRI applications are presented in this thesis. One motivation behind the current work was to reduce the acquisition time for relaxation mapping techniques DESPOT1 and DESPOT2. This work also includes a modification to the Driven Equilibrium Single Pulse Observation of T1 with high-speed incorporation of RF field inhomogeneities (DESPOT1-HIFI), resulting in more accurate estimation of T1 values at high strength (3T and higher) magnetic fields

Investigation of Sparsifying Transforms in Compressed Sensing for Magnetic Resonance Imaging with Fasttestcs

The goal of this contribution is to achieve higher reduction factors for faster Magnetic Resonance Imaging (MRI) scans with better Image Quality (IQ) by using Compressed Sensing (CS). This can be accomplished by adopting and understanding better sparsifying transforms for CS in MRI. There is a tremendous number of transforms and optional settings potentially available. Additionally, the amount of research in CS is growing, with possible duplication and difficult practical evaluation and comparison. However, no in-depth analysis of the effectiveness of different redundant sparsifying transforms on MRI images with CS has been undertaken until this work. New theoretical sparsity bounds for the dictionary restricted isometry property constants in CS are presented with mathematical proof. In order to verify the sparsifying transforms in this setting, the experiments focus on several redundant transforms contrasting them with orthogonal transforms. The transforms investigated are Wavelet (WT), Cosine (CT), contourlet, curvelet, k-means singular value decomposition, and Gabor. Several variations of these transforms with corresponding filter options are developed and tested in compression and CS simulations. Translation Invariance (TI) in transforms is found to be a key contributing factor in producing good IQ because any particular translation of the signal will not effect the transform representation. Some transforms tested here are TI and many others are made TI by transforming small overlapping image patches. These transforms are tested by comparing different under-sampling patterns and reduction ratios with varying image types including MRI data. Radial, spiral, and various random patterns are implemented and demonstrate that the TIWT is very robust across all under-sampling patterns. Results of the TIWT simulations show improvements in de-noising and artifact suppression over that of individual orthogonal wavelets and total variation ell-1 minimization in CS simulations. Some of these transforms add considerable time to the CS simulations and prohibit extensive testing of large 3D MRI datasets. Therefore, the FastTestCS software simulation framework is developed and customized for testing images, under-samping patterns and sparsifying transforms. This novel software is offered as a practical, robust, universal framework for evaluating and developing simulations in order to quickly test sparsifying transforms for CS MRI

Integrated Segmentation and Interpolation of Sparse Data

This paper addresses the two inherently related problems of segmentation and interpolation of 3D and 4D sparse data by integrating integrate these stages in a level set framework. The method supports any spatial configurations of sets of 2D slices having arbitrary positions and orientations. We introduce a new level set scheme based on the interpolation of the level set function by radial basis functions. The proposed method is validated quantitatively and/or subjectively on artificial data and MRI and CT scans and is compared against the traditional sequential approach