Sparsity-Based Recovery of Three-Dimensional Photoacoustic Images from Compressed Single-Shot Optical Detection.

Photoacoustic (PA) imaging combines optical excitation with ultrasonic detection to achieve high-resolution imaging of biological samples. A high-energy pulsed laser is often used for imaging at multi-centimeter depths in tissue. These lasers typically have a low pulse repetition rate, so to acquire images in real-time, only one pulse of the laser can be used per image. This single pulse necessitates the use of many individual detectors and receive electronics to adequately record the resulting acoustic waves and form an image. Such requirements make many PA imaging systems both costly and complex. This investigation proposes and models a method of volumetric PA imaging using a state-of-the-art compressed sensing approach to achieve real-time acquisition of the initial pressure distribution (IPD) at a reduced level of cost and complexity. In particular, a single exposure of an optical image sensor is used to capture an entire Fabry-Pérot interferometric acoustic sensor. Time resolved encoding as achieved through spatial sweeping with a galvanometer. This optical system further makes use of a random binary mask to set a predetermined subset of pixels to zero, thus enabling recovery of the time-resolved signals. The Two-Step Iterative Shrinking and Thresholding algorithm is used to reconstruct the IPD, harnessing the sparsity naturally occurring in the IPD as well as the additional structure provided by the binary mask. We conduct experiments on simulated data and analyze the performance of our new approach.

Sparsity-based photoacoustic image reconstruction with a linear array transducer and direct measurement of the forward model.

Photoacoustic (PA) imaging is an emerging imaging technique for many clinical applications. One of the challenges posed by clinical translation is that imaging systems often rely on a finite-aperture transducer rather than a full tomography system. This results in imaging artifacts arising from an underdetermined reconstruction of the initial pressure distribution (IPD). Furthermore, clinical applications often require deep imaging, resulting in a low-signal-to-noise ratio for the acquired signal because of strong light attenuation in tissue. Conventional approaches to reconstruct the IPD, such as back projection and time-reversal, do not adequately suppress the artifacts and noise. We propose a sparsity-based optimization approach that improves the reconstruction of IPD in PA imaging with a linear array ultrasound transducer. In simulation studies, the forward model matrix was measured from k-Wave simulations, and the approach was applied to reconstruct simulated point objects and the Shepp-Logan phantom. The results were compared with the conventional back projection, time-reversal approach, frequency-domain reconstruction, and the iterative least-squares approaches. In experimental studies, the forward model of our imaging system is directly measured by scanning a graphite point source through the imaging field of view. Experimental images of graphite inclusions in tissue-mimicking phantoms are reconstructed and compared with the back projection and iterative least-squares approaches. Overall these results show that our proposed optimization approach can leverage the sparsity of the PA images to improve the reconstruction of the IPD and outperform the existing popular reconstruction approaches.

Recovering fine details from under-resolved electron tomography data using higher order total variation ℓ1 regularization.

Over the last decade or so, reconstruction methods using ℓ1 regularization, often categorized as compressed sensing (CS) algorithms, have significantly improved the capabilities of high fidelity imaging in electron tomography. The most popular ℓ1 regularization approach within electron tomography has been total variation (TV) regularization. In addition to reducing unwanted noise, TV regularization encourages a piecewise constant solution with sparse boundary regions. In this paper we propose an alternative ℓ1 regularization approach for electron tomography based on higher order total variation (HOTV). Like TV, the HOTV approach promotes solutions with sparse boundary regions. In smooth regions however, the solution is not limited to piecewise constant behavior. We demonstrate that this allows for more accurate reconstruction of a broader class of images - even those for which TV was designed for - particularly when dealing with pragmatic tomographic sampling patterns and very fine image features. We develop results for an electron tomography data set as well as a phantom example, and we also make comparisons with discrete tomography approaches.

Roughness models for particle adhesion.

The effects of different surface roughness models on a previously developed van der Waals adhesion model were examined. The van der Waals adhesion model represented surface roughness with a distribution of hemispherical asperities. It was found that the constraints used to define the asperity distribution on the surface, which were determined from AFM scans, varied with scan size and thus were not constant for all surfaces examined. The greatest variation in these parameters occurred with materials that had large asperities or with materials where a large fraction of the surface was covered by asperities. These rough surfaces were modeled with fractals and also with a fast Fourier transform algorithm. When the model surfaces generated using the Fourier transforms are used in the adhesion model, the model accurately predicts the experimentally observed adhesion forces measured with the AFM.

Improving the accuracy of volumetric segmentation using pre-processing boundary detection and image reconstruction.

The concentration edge -detection and Gegenbauer image-reconstruction methods were previously shown to improve the quality of segmentation in magnetic resonance imaging. In this study, these methods are utilized as a pre-processing step to the Weibull E-SD field segmentation. It is demonstrated that the combination of the concentration edge detection and Gegenbauer reconstruction method improves the accuracy of segmentation for the simulated test data and real magnetic resonance images used in this study.

Improving tissue segmentation of human brain MRI through preprocessing by the Gegenbauer reconstruction method.

The Gegenbauer image reconstruction method, previously shown to improve the quality of magnetic resonance images, is utilized in this study as a segmentation preprocessing step. It is demonstrated that, for all simulated and real magnetic resonance images used in this study, the Gegenbauer reconstruction method improves the accuracy of segmentation. Although it is more desirable to use the k-space data for the Gegenbauer reconstruction method, only information acquired from MR images is necessary for the reconstruction, making the procedure completely self-contained and viable for all human brain segmentation algorithms.

A method to reduce the Gibbs ringing artifact in MRI scans while keeping tissue boundary integrity.

Gibbs ringing is a well known artifact that effects reconstruction of images having discontinuities. This is a problem in the reconstruction of magnetic resonance imaging (MRI) data due to the many different tissues normally present in each scan. The Gibbs ringing artifact manifests itself at the boundaries of the tissues, making it difficult to determine the structure of the brain tissue. The Gegenbauer reconstruction method has been shown to effectively eliminate the effects of Gibbs ringing in other applications. This paper presents the application of the Gegenbauer reconstruction method to neuro-imaging.