Model‐based optoacoustic inversion with arbitrary‐shape detectors

source: © 2011 American Association of Physicists in Medicine

Purpose:
Optoacoustic imaging enables mapping the optical absorption of biological tissue using optical excitation and acoustic detection. Although most image‐reconstruction algorithms are based on the assumption of a detector with an isotropic sensitivity, the geometry of the detector often leads to a response with spatially dependent magnitude and bandwidth. This effect may lead to attenuation or distortion in the recorded signal and, consequently, in the reconstructed image.

Methods:
Herein, an accurate numerical method for simulating the spatially dependent response of an arbitrary‐shape acoustic transducer is presented. The method is based on an analytical solution obtained for a two‐dimensional line detector. The calculated response is incorporated in the forward model matrix of an optoacoustic imaging setup using temporal convolution, and image reconstruction is performed by inverting the matrix relation.  [Read more…]

Fig. 8 Experimental reconstructions of a point optoacoustic source detected by a flat detector with a width of 1.3 cm obtained using (a) the back‐projection algorithm (b) IMMI modeled with a point detector (c) IMMI modeled with a 1.3‐mm flat detector using spatial convolution and (d) IMMI modeled with a 1.3‐mm flat detector using temporal convolution. The point source was obtained by applying plane‐selective illuminating on a black hair embedded in a clear agar phantom, as shown in Fig. 6(a). Although both the spatial‐ and temporal‐convolution methods managed enhancing the reconstruction resolution, the temporal‐convolution method yielded a more accurate reconstruction with less background texture.

Results:
The method was numerically and experimentally demonstrated in two dimensions for both flat and focused transducers and compared to the spatial‐convolution method. In forward simulations, the developed method did not suffer from the numerical errors exhibited by the spatial‐convolution method. In reconstruction simulations and experiments, the use of both temporal‐convolution and spatial‐convolution methods lead to an enhancement in resolution compared to a reconstruction with a point detector model. However, because of its higher modeling accuracy, the temporal‐convolution method achieved a noise figure approximated three times lower than the spatial‐convolution method.

Conclusions:
The demonstrated performance of the spatial‐convolution method shows it is a powerful tool for reducing reconstruction artifacts originating from the detector finite size and improving the quality of optoacoustic reconstructions. Furthermore, the method may be used for assessing new system designs. Specifically, detectors with nonstandard shapes may be investigated.

Amir Rosenthal Vasilis Ntziachristos Daniel Razansky, “Model‐based optoacoustic inversion with arbitrary‐shape detectors,” Medical Physics Volume38, Issue7,July 2011,Pages 4285-4295