Open Access
1 March 2010 k-Wave: MATLAB toolbox for the simulation and reconstruction of photoacoustic wave fields
Author Affiliations +
Abstract
A new, freely available third party MATLAB toolbox for the simulation and reconstruction of photoacoustic wave fields is described. The toolbox, named k-Wave, is designed to make realistic photoacoustic modeling simple and fast. The forward simulations are based on a k-space pseudo-spectral time domain solution to coupled first-order acoustic equations for homogeneous or heterogeneous media in one, two, and three dimensions. The simulation functions can additionally be used as a flexible time reversal image reconstruction algorithm for an arbitrarily shaped measurement surface. A one-step image reconstruction algorithm for a planar detector geometry based on the fast Fourier transform (FFT) is also included. The architecture and use of the toolbox are described, and several novel modeling examples are given. First, the use of data interpolation is shown to considerably improve time reversal reconstructions when the measurement surface has only a sparse array of detector points. Second, by comparison with one-step, FFT-based reconstruction, time reversal is shown to be sufficiently general that it can also be used for finite-sized planar measurement surfaces. Last, the optimization of computational speed is demonstrated through parallel execution using a graphics processing unit.
©(2010) Society of Photo-Optical Instrumentation Engineers (SPIE)
Bradley E. Treeby and Benjamin T. Cox "k-Wave: MATLAB toolbox for the simulation and reconstruction of photoacoustic wave fields," Journal of Biomedical Optics 15(2), 021314 (1 March 2010). https://doi.org/10.1117/1.3360308
Published: 1 March 2010
Lens.org Logo
CITATIONS
Cited by 1629 scholarly publications and 19 patents.
Advertisement
Advertisement
RIGHTS & PERMISSIONS
Get copyright permission  Get copyright permission on Copyright Marketplace
KEYWORDS
Sensors

Photoacoustic spectroscopy

Acoustics

Image restoration

Reconstruction algorithms

MATLAB

Computer simulations

Back to Top