Nam Huynh

Accelerated high-resolution photoacoustic tomography via compressed sensing.

Current 3D photoacoustic tomography (PAT) systems offer either high image quality or high frame rates but are not able to deliver high spatial and temporal resolution simultaneously, which limits their ability to image dynamic processes in living tissue (4D PAT). A particular example is the planar Fabry-Pérot (FP) photoacoustic scanner, which yields high-resolution 3D images but takes several minutes to sequentially map the incident photoacoustic field on the 2D sensor plane, point-by-point. However, as the spatio-temporal complexity of many absorbing tissue structures is rather low, the data recorded in such a conventional, regularly sampled fashion is often highly redundant. We demonstrate that combining model-based, variational image reconstruction methods using spatial sparsity constraints with the development of novel PAT acquisition systems capable of sub-sampling the acoustic wave field can dramatically increase the acquisition speed while maintaining a good spatial resolution: first, we describe and model two general spatial sub-sampling schemes. Then, we discuss how to implement them using the FP interferometer and demonstrate the potential of these novel compressed sensing PAT devices through simulated data from a realistic numerical phantom and through measured data from a dynamic experimental phantom as well as from in vivo experiments. Our results show that images with good spatial resolution and contrast can be obtained from highly sub-sampled PAT data if variational image reconstruction techniques that describe the tissues structures with suitable sparsity-constraints are used. In particular, we examine the use of total variation (TV) regularization enhanced by Bregman iterations. These novel reconstruction strategies offer new opportunities to dramatically increase the acquisition speed of photoacoustic scanners that employ point-by-point sequential scanning as well as reducing the channel count of parallelized schemes that use detector arrays.

Single-pixel optical camera for video rate ultrasonic imaging

A coherent-light single-pixel camera was used to interrogate a Fabry–Perot polymer film ultrasound sensor, thereby serially encoding a time-varying 2D ultrasonic field onto a single optical channel. By facilitating compressive sensing, this device enabled video rate imaging of ultrasound fields. In experimental demonstrations, this compressed sensing capability was exploited to reduce motion blur and capture dynamic features in the data. This relatively simple and inexpensive proof-of-principle device offers a route to high pixel count, high frame rate, broadband 2D ultrasound field mapping.

Photoacoustic imaging using an 8-beam Fabry-Perot scanner

The planar Fabry Perot (FP) photoacoustic scanner has been shown to provide exquisite high resolution 3D images of soft tissue structures in vivo to depths up to approximately 10mm. However a significant limitation of current embodiments of the concept is low image acquisition speed. To increase acquisition speed, a novel multi-beam scanner architecture has been developed. This enables a line of equally spaced 8 interrogation beams to be scanned simultaneously across the FP sensor and the photoacoustic signals detected in parallel. In addition, an excitation laser operating at 200Hz was used. The combination of parallelising the detection and the high pulse repetition frequency (PRF) of the excitation laser has enabled dramatic reductions in image acquisition time to be achieved. A 3D image can now be acquired in 10 seconds and 2D images at video rates are now possible.

Ultrasound modulated imaging of luminescence generated within a scattering medium

Abstract. Ultrasound modulated optical tomography modulates scattered light within tissue by deterministically altering the optical properties of the sample with the ultrasonic pressure. This allows the light to be “tagged” and the degradation in spatial resolution associated with light scattering to be reduced. To our knowledge, this is the first demonstration of ultrasound modulated imaging of light generated within a scattering medium without an external light source. The technique has the potential to improve the spatial resolution of chemi- or bioluminescence imaging of tissue. Experimental results show that ultrasound modulated luminescence imaging can resolve two chemiluminescent objects separated by 5 mm at a 7 mm depth within a tissue phantom with a scattering coefficient of 30  cm−1. The lateral resolution is estimated to be 3 mm. Monte Carlo simulations indicate that, with the current system signal to noise ratio, it is feasible to apply the approach to bioluminescence imaging when the concentration of bacteria in the animal organ is above 3.4×105/μL.

Model-Based Learning for Accelerated, Limited-View 3-D Photoacoustic Tomography

Recent advances in deep learning for tomographic reconstructions have shown great potential to create accurate and high quality images with a considerable speed up. In this paper, we present a deep neural network that is specifically designed to provide high resolution 3-D images from restricted photoacoustic measurements. The network is designed to represent an iterative scheme and incorporates gradient information of the data fit to compensate for limited view artifacts. Due to the high complexity of the photoacoustic forward operator, we separate training and computation of the gradient information. A suitable prior for the desired image structures is learned as part of the training. The resulting network is trained and tested on a set of segmented vessels from lung computed tomography scans and then applied to in-vivo photoacoustic measurement data.

Effect of object size and acoustic wavelength on pulsed ultrasound modulated fluorescence signals

Detection of ultrasound (US)-modulated fluorescence in turbid media is a challenge because of the low level of fluorescent light and the weak modulation of incoherent light. A very limited number of theoretical and experimental investigations have been performed, and this is, to our knowledge, the first demonstration of pulsed US-modulated fluorescence tomography. Experimental results show that the detected signal depends on the acoustic frequency and the fluorescent targets size along the ultrasonic propagation axis. The modulation depth of the detected signal is greatest when the length of the object along the acoustic axis is an odd number of half wavelengths and is weakest when the object is an integer multiple of an acoustic wavelength. Images of a fluorescent tube embedded within a 22- by 13- by 30 mm scattering gel phantom (μ(s)∼15  cm(-1), g=0.93) with 1-, 1.5-, and 2 MHz frequency US are presented. The modulation depth of the detected signal changes by a factor of 5 depending on the relative size of the object and the frequency. The approach is also verified by some simple experiments in a nonscattering gel and using a theoretical model.

A real-time ultrasonic field mapping system using a Fabry Pérot single pixel camera for 3D photoacoustic imaging

A system for dynamic mapping of broadband ultrasound fields has been designed, with high frame rate photoacoustic imaging in mind. A Fabry-Pérot interferometric ultrasound sensor was interrogated using a coherent light single-pixel camera. Scrambled Hadamard measurement patterns were used to sample the acoustic field at the sensor, and either a fast Hadamard transform or a compressed sensing reconstruction algorithm were used to recover the acoustic pressure data. Frame rates of 80 Hz were achieved for 32x32 images even though no specialist hardware was used for the on-the-fly reconstructions. The ability of the system to obtain photocacoustic images with data compressions as low as 10% was also demonstrated.

Patterned interrogation scheme for compressed sensing photoacoustic imaging using a Fabry Perot planar sensor

Photoacoustic tomography (PAT) has become a powerful tool for biomedical imaging, particularly pre-clinical small animal imaging. Several different measurement systems have been demonstrated, in particular, optically addressed Fabry-Perot interferometer (FPI) sensors have been shown to provide exquisite images when a planar geometry is suitable. However, in its current incarnation the measurements must be made at each point sequentially, so these devices therefore suffer from slow data acquisition time. An alternative to this point-by-point interrogation scheme, is to interrogate the whole sensor with a series of independent patterns, so each measurement is the spatial integral of the product of the pattern and the acoustic field (as in the single-pixel Rice camera). Such an interrogation scheme allows compressed sensing to be used. This enables the number of measurements to be reduced significantly, leading to much faster data acquisition. An experimental implementation will be described, which employs a wide NIR tunable laser beam to interrogate the FPI sensor. The reflected beam is patterned by a digital micro-mirror device, and then focused to a single photodiode. To demonstrate the idea of patterned and compressed sensing for ultrasound detection, a scrambled Hadamard operator is used in the experiments. Photoacoustic imaging experiments of phantoms shows good reconstructed results with 20% compression.

Sub-sampled Fabry-Perot photoacoustic scanner for fast 3D imaging

The planar Fabry Perot (FP) photoacoustic scanner provides exquisite high resolution 3D images of soft tissue structures for sub-cm penetration depths. However, as the FP sensor is optically addressed by sequentially scanning an interrogation laser beam over its surface, the acquisition speed is low. To address this, a novel scanner architecture employing 8 interrogation beams and an optimised sub-sampling framework have been developed that increase the data acquisition speed significantly. With a 200Hz repetition rate excitation laser, full 3D images can be obtained within 10 seconds. Further increases in imaging speed with only minor decreases in image quality can be obtained by applying sub-sampling techniques with rates as low as 12.5%. This paper shows 3D images reconstructed from sub-sampled data for an ex vivo dataset, and results from a dynamic phantom imaging experiment.

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