Acoustic Wave Field Reconstruction From Compressed Measurements With Application in Photoacoustic Tomography

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.

[1]  Ting Sun,et al.  Single-pixel imaging via compressive sampling , 2008, IEEE Signal Process. Mag..

[2]  Markus Haltmeier,et al.  Compressed sensing and sparsity in photoacoustic tomography , 2016, 1605.09249.

[3]  Lihong V Wang,et al.  Universal back-projection algorithm for photoacoustic computed tomography , 2005, SPIE BiOS.

[4]  E. Candès.  The restricted isometry property and its implications for compressed sensing , 2008 .

[5]  B T Cox,et al.  Fast calculation of pulsed photoacoustic fields in fluids using k-space methods. , 2005, The Journal of the Acoustical Society of America.

[6]  Massimo Fornasier,et al.  Compressive Sensing , 2015, Handbook of Mathematical Methods in Imaging.

[7]  Linh V. Nguyen.  On singularities and instability of reconstruction in thermoacoustic tomography , 2009 .

[8]  S. Frick,et al.  Compressed Sensing , 2014, Computer Vision, A Reference Guide.

[9]  Otmar Scherzer,et al.  A direct method for photoacoustic tomography with inhomogeneous sound speed , 2015, 1507.01741.

[10]  Markus Haltmeier,et al.  Universal Inversion Formulas for Recovering a Function from Spherical Means , 2012, SIAM J. Math. Anal..

[11]  Peter Kuchment,et al.  Mathematics of Photoacoustic and Thermoacoustic Tomography , 2009, Handbook of Mathematical Methods in Imaging.

[12]  E. Candès,et al.  Stable signal recovery from incomplete and inaccurate measurements , 2005, math/0503066.

[13]  E. Candès,et al.  New tight frames of curvelets and optimal representations of objects with piecewise C2 singularities , 2004 .

[14]  B T Cox,et al.  k-Wave: MATLAB toolbox for the simulation and reconstruction of photoacoustic wave fields. , 2010, Journal of biomedical optics.

[15]  Otmar Scherzer,et al.  Handbook of Mathematical Methods in Imaging , 2015, Handbook of Mathematical Methods in Imaging.

[16]  G. Uhlmann,et al.  Thermoacoustic tomography with variable sound speed , 2009, 0902.1973.

[17]  Eric Todd Quinto,et al.  Artifacts in Incomplete Data Tomography with Applications to Photoacoustic Tomography and Sonar , 2014, SIAM J. Appl. Math..

[18]  Mark A. Anastasio,et al.  Full-Wave Iterative Image Reconstruction in Photoacoustic Tomography With Acoustically Inhomogeneous Media , 2013, IEEE Transactions on Medical Imaging.

[19]  Lihong V. Wang.  Photoacoustic imaging and spectroscopy , 2009 .

[20]  Paul C. Beard,et al.  A real-time ultrasonic field mapping system using a Fabry Pérot single pixel camera for 3D photoacoustic imaging , 2015, Photonics West - Biomedical Optics.

[21]  P. Beard.  Biomedical photoacoustic imaging , 2011, Interface Focus.

[22]  Aswin C. Sankaranarayanan,et al.  Compressive Sensing , 2008, Computer Vision, A Reference Guide.

[23]  E. Candès,et al.  The curvelet representation of wave propagators is optimally sparse , 2004, math/0407210.

[24]  Yuan Xu,et al.  Exact frequency-domain reconstruction for thermoacoustic tomography. I. Planar geometry , 2002, IEEE Transactions on Medical Imaging.

[25]  Marta Betcke,et al.  Accelerated high-resolution photoacoustic tomography via compressed sensing , 2016, Physics in medicine and biology.

[26]  Alistair P. Rendell,et al.  Modeling nonlinear ultrasound propagation in heterogeneous media with power law absorption using a k-space pseudospectral method. , 2012, The Journal of the Acoustical Society of America.

[27]  Frédéric Lesage,et al.  The Application of Compressed Sensing for Photo-Acoustic Tomography , 2009, IEEE Transactions on Medical Imaging.

[28]  Paul C. Beard,et al.  Patterned interrogation scheme for compressed sensing photoacoustic imaging using a Fabry Perot planar sensor , 2014, Photonics West - Biomedical Optics.

[29]  Jan Laufer,et al.  Backward-mode multiwavelength photoacoustic scanner using a planar Fabry-Perot polymer film ultrasound sensor for high-resolution three-dimensional imaging of biological tissues. , 2008, Applied optics.

[30]  José M. Bioucas-Dias,et al.  Fast Image Recovery Using Variable Splitting and Constrained Optimization , 2009, IEEE Transactions on Image Processing.

[31]  Simon R. Arridge,et al.  On the adjoint operator in photoacoustic tomography , 2016, 1602.02027.

[32]  Yong Zhou,et al.  Tutorial on photoacoustic tomography , 2016, Journal of biomedical optics.

[33]  Lihong V Wang,et al.  Compressed sensing in photoacoustic tomography in vivo. , 2010, Journal of biomedical optics.

[34]  K. Valluru,et al.  Photoacoustic Imaging in Oncology: Translational Preclinical and Early Clinical Experience. , 2016, Radiology.

[35]  Rakesh,et al.  Determining a Function from Its Mean Values Over a Family of Spheres , 2004, SIAM J. Math. Anal..

[36]  B T Cox,et al.  A k-space Green's function solution for acoustic initial value problems in homogeneous media with power law absorption. , 2011, The Journal of the Acoustical Society of America.

[37]  H. Weber,et al.  Temporal backward projection of optoacoustic pressure transients using fourier transform methods. , 2001, Physics in medicine and biology.

[38]  Lihong V. Wang,et al.  Small-Animal Whole-Body Photoacoustic Tomography: A Review , 2014, IEEE Transactions on Biomedical Engineering.

[39]  Paul C. Beard,et al.  Single-pixel optical camera for video rate ultrasonic imaging , 2016 .

[40]  Laurent Demanet,et al.  Fast Discrete Curvelet Transforms , 2006, Multiscale Model. Simul..

[41]  D. Donoho.  Sparse Components of Images and Optimal Atomic Decompositions , 2001 .

[42]  T. D. Mast,et al.  A k-space method for coupled first-order acoustic propagation equations. , 2002, The Journal of the Acoustical Society of America.

[43]  Emmanuel J. Candès,et al.  Robust uncertainty principles: exact signal reconstruction from highly incomplete frequency information , 2004, IEEE Transactions on Information Theory.