University College London, Department of Medical Physics and Bioengineering, Gower Street, London, WC1E 6BT, United Kingdom.
J Biomed Opt. 2010 Mar-Apr;15(2):021314. doi: 10.1117/1.3360308.
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.
描述了一种新的、免费可用的第三方 MATLAB 工具箱，用于模拟和重建光声波场。该工具箱名为 k-Wave，旨在使逼真的光声建模变得简单快捷。正向模拟基于 k 空间伪谱时域解耦的一阶声学方程，适用于一维、二维和三维的均匀或非均匀介质。模拟函数还可以用作任意形状测量表面的灵活时间反转图像重建算法。还包括一种基于快速傅里叶变换 (FFT) 的用于平面探测器几何形状的一步图像重建算法。本文描述了该工具箱的体系结构和使用方法，并给出了几个新的建模示例。首先，通过数据插值的使用，当测量表面只有稀疏的探测器点阵列时，时间反转重建得到了显著改善。其次，通过与一步、基于 FFT 的重建进行比较，时间反转被证明足够通用，因此它也可以用于有限大小的平面测量表面。最后，通过使用图形处理单元进行并行执行来演示计算速度的优化。
