Bradley E. Treeby

Bio: Bradley E. Treeby is an academic researcher from University College London. The author has contributed to research in topics: Time domain & Iterative reconstruction. The author has an hindex of 26, co-authored 131 publications receiving 4412 citations. Previous affiliations of Bradley E. Treeby include University College of Engineering & University of Western Australia.

k-Wave: MATLAB toolbox for the simulation and reconstruction of photoacoustic wave fields

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.

Modeling nonlinear ultrasound propagation in heterogeneous media with power law absorption using a k-space pseudospectral method

Abstract: The simulation of nonlinear ultrasound propagation through tissue realistic media has a wide range of practical applications. However, this is a computationally difficult problem due to the large size of the computational domain compared to the acoustic wavelength. Here, the k-space pseudospectral method is used to reduce the number of grid points required per wavelength for accurate simulations. The model is based on coupled first-order acoustic equations valid for nonlinear wave propagation in heterogeneous media with power law absorption. These are derived from the equations of fluid mechanics and include a pressure-density relation that incorporates the effects of nonlinearity, power law absorption, and medium heterogeneities. The additional terms accounting for convective nonlinearity and power law absorption are expressed as spatial gradients making them efficient to numerically encode. The governing equations are then discretized using a k-space pseudospectral technique in which the spatial gradients are computed using the Fourier-collocation method. This increases the accuracy of the gradient calculation and thus relaxes the requirement for dense computational grids compared to conventional finite difference methods. The accuracy and utility of the developed model is demonstrated via several numerical experiments, including the 3D simulation of the beam pattern from a clinical ultrasound probe.

Deep in vivo photoacoustic imaging of mammalian tissues using a tyrosinase-based genetic reporter

Abstract: Deep photoacoustic imaging of mammalian cells featuring genetically encoded contrast is reported.

Photoacoustic tomography in absorbing acoustic media using time reversal

Abstract: The reconstruction of photoacoustic images typically neglects the effect of acoustic absorption on the measured time domain signals. Here, a method to compensate for acoustic absorption in photoacoustic tomography is described. The approach is based on time-reversal image reconstruction and an absorbing equation of state which separately accounts for acoustic absorption and dispersion following a frequency power law. Absorption compensation in the inverse problem is achieved by reversing the absorption proportionality coefficient in sign but leaving the equivalent dispersion parameter unchanged. The reconstruction is regularized by filtering the absorption and dispersion terms in the spatial frequency domain using a Tukey window. This maintains the correct frequency dependence of these parameters within the filter pass band. The method is valid in one, two and three dimensions, and for arbitrary power law absorption parameters. The approach is verified through several numerical experiments. The reconstruction of a carbon fibre phantom and the vasculature in the abdomen of a mouse are also presented. When absorption compensation is included, a general improvement in the image magnitude and resolution is seen, particularly for deeper features.

In vivo preclinical photoacoustic imaging of tumor vasculature development and therapy.

Abstract: The use of a novel all-optical photoacoustic scanner for imaging the development of tumor vasculature and its response to a therapeutic vascular disrupting agent is described. The scanner employs a Fabry-Perot polymer film ultrasound sensor for mapping the photoacoustic waves and an image reconstruction algorithm based upon attenuation-compensated acoustic time reversal. The system was used to noninvasively image human colorectal tumor xenografts implanted subcutaneously in mice. Label-free three-dimensional in vivo images of whole tumors to depths of almost 10 mm with sub-100-micron spatial resolution were acquired in a longitudinal manner. This enabled the development of tumor-related vascular features, such as vessel tortuosity, feeding vessel recruitment, and necrosis to be visualized over time. The system was also used to study the temporal evolution of the response of the tumor vasculature following the administration of a therapeutic vascular disrupting agent (OXi4503). This revealed the well-known destruction and recovery phases associated with this agent. These studies illustrate the broader potential of this technology as an imaging tool for the preclinical and clinical study of tumors and other pathologies characterized by changes in the vasculature.

Photoacoustic Tomography: In Vivo Imaging from Organelles to Organs

Abstract: Photoacoustic tomography (PAT) can create multiscale multicontrast images of living biological structures ranging from organelles to organs. This emerging technology overcomes the high degree of scattering of optical photons in biological tissue by making use of the photoacoustic effect. Light absorption by molecules creates a thermally induced pressure jump that launches ultrasonic waves, which are received by acoustic detectors to form images. Different implementations of PAT allow the spatial resolution to be scaled with the desired imaging depth in tissue while a high depth-to-resolution ratio is maintained. As a rule of thumb, the achievable spatial resolution is on the order of 1/200 of the desired imaging depth, which can reach up to 7 centimeters. PAT provides anatomical, functional, metabolic, molecular, and genetic contrasts of vasculature, hemodynamics, oxygen metabolism, biomarkers, and gene expression. We review the state of the art of PAT for both biological and clinical studies and discuss future prospects.

Biomedical photoacoustic imaging

Abstract: Photoacoustic (PA) imaging, also called optoacoustic imaging, is a new biomedical imaging modality based on the use of laser-generated ultrasound that has emerged over the last decade. It is a hybrid modality, combining the high-contrast and spectroscopic-based specificity of optical imaging with the high spatial resolution of ultrasound imaging. In essence, a PA image can be regarded as an ultrasound image in which the contrast depends not on the mechanical and elastic properties of the tissue, but its optical properties, specifically optical absorption. As a consequence, it offers greater specificity than conventional ultrasound imaging with the ability to detect haemoglobin, lipids, water and other light-absorbing chomophores, but with greater penetration depth than purely optical imaging modalities that rely on ballistic photons. As well as visualizing anatomical structures such as the microvasculature, it can also provide functional information in the form of blood oxygenation, blood flow and temperature. All of this can be achieved over a wide range of length scales from micrometres to centimetres with scalable spatial resolution. These attributes lend PA imaging to a wide variety of applications in clinical medicine, preclinical research and basic biology for studying cancer, cardiovascular disease, abnormalities of the microcirculation and other conditions. With the emergence of a variety of truly compelling in vivo images obtained by a number of groups around the world in the last 2–3 years, the technique has come of age and the promise of PA imaging is now beginning to be realized. Recent highlights include the demonstration of whole-body small-animal imaging, the first demonstrations of molecular imaging, the introduction of new microscopy modes and the first steps towards clinical breast imaging being taken as well as a myriad of in vivo preclinical imaging studies. In this article, the underlying physical principles of the technique, its practical implementation, and a range of clinical and preclinical applications are reviewed.

k-Wave: MATLAB toolbox for the simulation and reconstruction of photoacoustic wave fields

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.

Measurement of intraocular distances by backscattering spectral interferometry

Abstract: The diffraction tomography theorem is adapted to one-dimensional length measurement. The resulting spectral interferometry technique is described and the first length measurements using this technique on a model eye and on a human eye in vivo are presented.