Accelerated High-Resolution Photoacoustic Tomography via Compressed Sensing

Current 3D photoacoustic tomography (PAT) systems offer either high image quality or high frame rates but are not able to deliver high spatial and temporal resolution simultaneously, which limits their ability to image dynamic processes in living tissue. A particular example is the planar Fabry-Perot (FP) scanner, which yields high-resolution images but takes several minutes to sequentially map the photoacoustic field on the sensor plane, point-by-point. However, as the spatio-temporal complexity of many absorbing tissue structures is rather low, the data recorded in such a conventional, regularly sampled fashion is often highly redundant. We demonstrate that combining variational image reconstruction methods using spatial sparsity constraints with the development of novel PAT acquisition systems capable of sub-sampling the acoustic wave field can dramatically increase the acquisition speed while maintaining a good spatial resolution: First, we describe and model two general spatial sub-sampling schemes. Then, we discuss how to implement them using the FP scanner and demonstrate the potential of these novel compressed sensing PAT devices through simulated data from a realistic numerical phantom and through measured data from a dynamic experimental phantom as well as from in-vivo experiments. Our results show that images with good spatial resolution and contrast can be obtained from highly sub-sampled PAT data if variational image reconstruction methods that describe the tissues structures with suitable sparsity-constraints are used. In particular, we examine the use of total variation regularization enhanced by Bregman iterations. These novel reconstruction strategies offer new opportunities to dramatically increase the acquisition speed of PAT scanners that employ point-by-point sequential scanning as well as reducing the channel count of parallelized schemes that use detector arrays.Comment: submitted to "Physics in Medicine and Biology

A sub-Nyquist co-prime sampling music spectral approach for natural frequency identification of white-noise excited structures

Motivated by practical needs to reduce data transmission payloads in wireless sensors for vibration-based monitoring of civil engineering structures, this paper proposes a novel approach for identifying resonant frequencies of white-noise excited structures using acceleration measurements acquired at rates significantly below the Nyquist rate. The approach adopts the deterministic co-prime sub-Nyquist sampling scheme, originally developed to facilitate telecommunication applications, to estimate the autocorrelation function of response acceleration time-histories of low-amplitude white-noise excited structures treated as realizations of a stationary stochastic process. This is achieved without posing any sparsity conditions to the signals. Next, the standard MUSIC algorithm is applied to the estimated autocorrelation function to derive a denoised super-resolution pseudo-spectrum in which natural frequencies are marked by prominent spikes. The accuracy and applicability of the proposed approach is numerically assessed using computer-generated noise-corrupted acceleration time-history data obtained by a simulation-based framework pertaining to a white-noise excited structural system with two closely-spaced modes of vibration carrying the same amount of energy, and a third isolated weakly excited vibrating mode. All three natural frequencies are accurately identified by sampling at as low as 78% below Nyquist rate for signal to noise ratio as low as 0dB (i.e., energy of additive white noise equal to the signal energy), suggesting that the proposed approach is robust and noise-immune while it can reduce data transmission requirements in acceleration wireless sensors for natural frequency identification of engineering structures

Accelerated Cardiac Diffusion Tensor Imaging Using Joint Low-Rank and Sparsity Constraints

Objective: The purpose of this manuscript is to accelerate cardiac diffusion tensor imaging (CDTI) by integrating low-rankness and compressed sensing. Methods: Diffusion-weighted images exhibit both transform sparsity and low-rankness. These properties can jointly be exploited to accelerate CDTI, especially when a phase map is applied to correct for the phase inconsistency across diffusion directions, thereby enhancing low-rankness. The proposed method is evaluated both ex vivo and in vivo, and is compared to methods using either a low-rank or sparsity constraint alone. Results: Compared to using a low-rank or sparsity constraint alone, the proposed method preserves more accurate helix angle features, the transmural continuum across the myocardium wall, and mean diffusivity at higher acceleration, while yielding significantly lower bias and higher intraclass correlation coefficient. Conclusion: Low-rankness and compressed sensing together facilitate acceleration for both ex vivo and in vivo CDTI, improving reconstruction accuracy compared to employing either constraint alone. Significance: Compared to previous methods for accelerating CDTI, the proposed method has the potential to reach higher acceleration while preserving myofiber architecture features which may allow more spatial coverage, higher spatial resolution and shorter temporal footprint in the future.Comment: 11 pages, 16 figures, published on IEEE Transactions on Biomedical Engineerin