Abstract

We present a method for the recovery of compressively sensed acoustic fields using patterned, instead of point-by-point, detection. From a limited number of such compressed measurements, we propose to reconstruct the field on the sensor plane in each time step independently assuming its sparsity in a Curvelet frame. A modification of the Curvelet frame is proposed to account for the smoothing effects of data acquisition and motivated by a frequency domain model for photoacoustic tomography. An ADMM type algorithm, split augmented Lagrangian shrinkage algorithm, is used to recover the pointwise data in each individual time step from the patterned measurements. For photoacoustic applications, the photoacoustic image of the initial pressure is reconstructed using time reversal in $ {\mathbf k}$ -Wave Toolbox.

Highlights

  • C OMPRESSED Sensing (CS) is a new measurement paradigm, which allows for the reconstruction of sparse signals sampled at sub-Nyquist rates

  • We discuss the motivation for using a Curvelet tranform as the sparsifying transformation for the acoustic field at the detector and we propose its modification: a low-frequency Curvelet tranform tailored to the frequency range of the acoustic field on the detector

  • An appropriate multiscale representation of the time series photoacoustic tomography (PAT) data is considered in Section IV, where we derive the frequency model of sensor data and propose a modified version of Curvelet transform tailored to the range of frequencies of the acoustic field

Read more

Summary

INTRODUCTION

C OMPRESSED Sensing (CS) is a new measurement paradigm, which allows for the reconstruction of sparse signals sampled at sub-Nyquist rates. It is common understanding that many digital signals and images admit an adequate representation with far fewer coefficients than their actual length. This phenomena is known as compressibility and it has been a driving force in many image processing applications, most notably the image compression algorithms JPEG and its successor JPEG 2000. As the Curvelets essentially describe the wave front sets their propagation is well approximated through geometrical optics (high frequency asymptotic solution to the wave equation). Of the corresponding Hamilton-Jacobi equation resulting from the high frequency asymptotic Motivated by this result we investigate the Curvelet representation of the cross-section of the wave field by the planar ultrasound sensor. While the arguments in [5] do not directly apply to this situation, the planar cross-section through the acoustic wave front constitutes a singularity along a smooth curve for which Curvelets have been demonstrated to be a nearly optimal representation [6]

Contribution
Related work
Outline
PHOTOACOUSTIC TOMOGRAPHY
Frequency model of sensor data
COMPRESSED SENSING Let Ψ be a sparsifying transform
Low-frequency Curvelets
Single-pixel optical camera
Sparse representation of the sensor data
Sensing patterns
Recovery of the sensor data
PAT RECONSTRUCTION WITH THE RECOVERED
Simulated data: clock phantom
Synthesized pattern data: knotted tubes filled with ink
Experimental pattern data: hair knot
Findings
CONCLUSIONS AND DISCUSSION

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call

Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.