Comparison of Piezoelectric and Optical Projection Imaging for Three-Dimensional In Vivo Photoacoustic Tomography.

Ultrasound sensor arrays for photoacoustic tomography (PAT) are investigated that create line projections of the pressure generated in an object by pulsed light illumination. Projections over a range of viewing angles enable the reconstruction of a three-dimensional image. Two line-integrating arrays are compared in this study for the in vivo imaging of vasculature, a piezoelectric array, and a camera-based setup that captures snapshots of the acoustic field emanating from the sample. An array consisting of 64 line-shaped sensors made of piezoelectric polymer film, which was arranged on a half-cylindrical area, was used to acquire spatiotemporal data from a human finger. The optical setup used phase contrast to visualize the acoustic field generated in the leg of a mouse after a selected delay time. Time-domain back projection and frequency-domain back propagation were used for image reconstruction from the piezoelectric and optical data, respectively. The comparison yielded an about threefold higher resolution for the optical setup and an about 13-fold higher sensitivity of the piezoelectric array. Due to the high density of data in the camera images, the optical technique gave images without streak artifacts, which were visible in the piezo array images due to the discrete detector positions. Overall, both detection concepts are suited for almost real-time projection imaging and three-dimensional imaging with a data acquisition time of less than a minute without averaging, which was limited by the repetition rate of the laser.

Simultaneous three-dimensional photoacoustic and laser-ultrasound tomography.

A tomographic setup that provides the co-registration of photoacoustic (PA) and ultrasound (US) images is presented. For pulse-echo US-tomography laser-induced broadband plane ultrasonic waves are produced by illuminating an optically absorbing target with a short near-infrared laser pulse. Part of the same pulse is frequency doubled and used for the generation of PA waves within the object of interest. The laser-generated plane waves are scattered at the imaging object and measured with the same interferometric detector that also acquires the photoacoustic signals. After collection and separation of the data image reconstruction is done using back-projection resulting in three-dimensional, co-registered PA and US images. The setup is characterized and the resolution in PA and US mode is estimated to be about 85 µm and 40 µm, respectively. Besides measurements on phantoms the performance is also tested on a biological sample.

On the use of frequency-domain reconstruction algorithms for photoacoustic imaging.

We investigate the use of a frequency-domain reconstruction algorithm based on the nonuniform fast Fourier transform (NUFFT) for photoacoustic imaging (PAI). Standard algorithms based on the fast Fourier transform (FFT) are computationally efficient, but compromise the image quality by artifacts. In our previous work we have developed an algorithm for PAI based on the NUFFT which is computationally efficient and can reconstruct images with the quality known from temporal backprojection algorithms. In this paper we review imaging qualities, such as resolution, signal-to-noise ratio, and the effects of artifacts in real-world situations. Reconstruction examples show that artifacts are reduced significantly. In particular, image details with a larger distance from the detectors can be resolved more accurately than with standard FFT algorithms.

Attenuation of ultrasound in severely plastically deformed nickel.

Ultrasound attenuation was measured in nickel specimens of about 30 mm diameter prepared using the high pressure torsion technique. The cold working process produced an equivalent shear strain increasing from zero at the center up to 1000% at the edge of the specimen. The fragmentation of the grains due to multiple dislocations led to an ultrafine microstructure with large angle grain boundaries. The mean value of the grain size distribution gradually decreased from ∼50 µm at the center to 0.2 µm at the edge. Laser pulses of 5 ns were employed for the excitation of broadband ultrasound pulses covering the spectral range of 0.1-150 MHz. The ultrasound pulses were measured from the opposite side of the specimen by means of an optical interferometer and a piezoelectric foil transducer in two experimental setups. The features of the detected signal forms are discussed. The absolute value of the attenuation decreases from the center to the edge of the specimen showing nearly linear frequency dependence. The variation of the phase velocity was measured in a 6 mm-thick high pressure torsion nickel sample, revealing a velocity increase from the center to the edge.

Photoacoustic tomography of ex vivo mouse hearts with myocardial infarction.

In the present study, we evaluated the applicability of ex vivo photoacoustic imaging (PAI) on small animal organs. We used photoacoustic tomography (PAT) to visualize infarcted areas within murine hearts and compared these data to other imaging techniques [magnetic resonance imaging (MRI), micro-computed tomography] and histological slices. In order to induce ischemia, an in vivo ligation of the left anterior descending artery was performed on nine wild-type mice. After varying survival periods, the hearts were excised and fixed in formaldehyde. Samples were illuminated with nanosecond laser pulses delivered by a Nd:YAG pumped optical parametric oscillator. Ultrasound detection was achieved using a Mach-Zehnder interferometer (MZI) working as an integrating line detector. The voxel data were computed using a Fourier-domain based reconstruction algorithm, followed by inverse Radon transforms. The results clearly showed the capability of PAI to visualize myocardial infarction and to produce three-dimensional images with a spatial resolution of approximately 120 µm. Regions of affected muscle tissue in PAI corresponded well with the results of MRI and histology. Photoacoustic tomography utilizing a MZI for ultrasound detection allows for imaging of small tissue samples. Due to its high spatial resolution, good soft tissue contrast and comparatively low cost, PAT offers great potentials for imaging.

Photoacoustic microtomography using optical interferometric detection.

A device for three-dimensional (3-D) photoacoustic tomography with resolution in the range of tens of micrometers is presented that uses a light beam for interferometric detection of acoustic waves. Reconstruction of the 3-D initial pressure distribution from the signals representing line integrals of the acoustic field is a two-step process. It uses an inversion of 2-D wave propagation to obtain line projections of the initial pressure distribution and the inverse Radon transform. The light beam, propagating freely in a water bath, is scanned either in an arc- or box-shaped curve around the object. Simulations are performed to compare the two scanning procedures. The projection images are obtained either using the filtered back projection algorithm for the pi-arc scanning mode or the frequency domain algorithm for the box scanning mode. While the former algorithm provides slightly better image quality, the latter is about 20 times faster. The ability of the photoacoustic tomography device to create 3-D images with constant resolution throughout the reconstruction volume is demonstrated experimentally using a human hair phantom. These measurements revealed a 3-D resolution below 100 mum. In a second experiment, 3-D imaging of an isolated mouse heart is demonstrated to show the applicability for preclinical and biological research.

Three-dimensional photoacoustic imaging using fiber-based line detectors.

For photoacoustic imaging, usually point-like detectors are used. As a special sensing technology for photoacoustic imaging, integrating detectors have been investigated that integrate the acoustic pressure over an area or line that is larger than the imaged object. Different kinds of optical fiber-based detectors are compared regarding their sensitivity and resolution in three-dimensional photoacoustic tomography. In the same type of interferometer, polymer optical fibers yielded much higher sensitivity than glass fibers. Fabry-Perot glass-fiber interferometers in turn gave higher sensitivity than Mach-Zehnder-type interferometers. Regarding imaging resolution, the single-mode glass fiber showed the best performance. Last, three-dimensional images of phantoms and insects using a glass-fiber-based Fabry-Perot interferometer as integrating line detector are presented.

Full field detection in photoacoustic tomography.

Imaging the full acoustic field around an object by use of an optical phase contrast method is used to accelerate the data acquisition in photoacoustic tomography. Images obtained with a CCD-camera at a certain time show a projection of the instantaneous pressure field in a given direction. In this work a reconstruction method is presented to obtain the two-dimensional initial pressure distribution by back propagating the observed wave pattern in Radon space. Numerical simulations are used to prove the accuracy of the reconstruction algorithm and to demonstrate a method for correcting limited data artifacts. Finally, the overall performance is shown with experimentally obtained data.

Femtosecond-laser-induced nanocavitation in water: implications for optical breakdown threshold and cell surgery.

We determined the bubble radius R_(max) for femtosecond optical breakdown in water at 347, 520, and 1040 nm with an unprecedented accuracy (+/-10 nm). At threshold, R_(max) was smaller than the diffraction-limited focus radius and ranged from 190 nm to 320 nm. The increase of R_(max) with laser energy E_(L) is slowest at 347 nm, providing optimum control of cell surgery. Experimental results agree with a model of bubble formation in heated and thermoelastically stretched liquids. Theory predicts a threshold temperature T_(th) approximately equal to 168 degrees C. For T>300 degrees C, a phase explosion sets in, and R_(max) increases rapidly with E_(L).

Exact and approximative imaging methods for photoacoustic tomography using an arbitrary detection surface.

Two universal reconstruction methods for photoacoustic (also called optoacoustic or thermoacoustic) computed tomography are derived, applicable to an arbitrarily shaped detection surface. In photoacoustic tomography acoustic pressure waves are induced by illuminating a semitransparent sample with pulsed electromagnetic radiation and are measured on a detection surface outside the sample. The imaging problem consists in reconstructing the initial pressure sources from those measurements. The first solution to this problem is based on the time reversal of the acoustic pressure field with a second order embedded boundary method. The pressure on the arbitrarily shaped detection surface is set to coincide with the measured data in reversed temporal order. In the second approach the reconstruction problem is solved by calculating the far-field approximation, a concept well known in physics, where the generated acoustic wave is approximated by an outgoing spherical wave with the reconstruction point as center. Numerical simulations are used to compare the proposed universal reconstruction methods with existing algorithms.

Thermoacoustic tomography with integrating area and line detectors.

Thermoacoustic (optoacoustic, photoacoustic) tomography is based on the generation of acoustic waves by illumination of a sample with a short electromagnetic pulse. The absorption density inside the sample is reconstructed from the acoustic pressure measured outside the illuminated sample. So far measurement data have been collected with small detectors as approximations of point detectors. Here, a novel measurement setup applying integrating detectors (e.g., lines or planes made of piezoelectric films) is presented. That way, the pressure is integrated along one or two dimensions, enabling the use of numerically efficient algorithms, such as algorithms for the inverse radon transformation, for thermoacoustic tomography. To reconstruct a three-dimensional sample, either an area detector has to be moved tangential around a sphere that encloses the sample or an array of line detectors is rotated around a single axis. The line detectors can be focused on cross sections perpendicular to the rotation axis using a synthetic aperture (SAFT) or by scanning with a cylindrical lens detector. Measurements were made with piezoelectric polyvinylidene fluoride film detectors and evaluated by comparison with numerical simulations. The resolution achieved in the resulting tomography images is demonstrated on the example of the reconstructed cross section of a grape.

[Optoacoustic spectroscopy and imaging]. / Optoakustische Spektroskopie und Bildgebung.

Optoacoustic methods are based on the thermoelastic effect, according to which a short laser pulse is absorbed in a medium and generates an acoustic wave following rapid heating and thermal expansion. This study it sought to demonstrate how the optical properties and structure of the sample can be obtained from acoustic signals that are measured on the surface of the medium. This was pursued on the hand by using broadband piezo-electric and optical sensors, and on the other hand by means of theoretical simulations. For optically one-dimensional media, where the optical properties depend only on one direction in space, the measurement of a plane wave in the acoustical near field of the sample yielded directly the depth distribution of absorbed energy. In the far field, this simple relationship was not valid, due to acoustic diffraction. The distribution, however, could be reconstructed mathematically. Far field signals could also be used to image the contours of optical structures.