Efficient testing of weld seam models with radii of curvature in the millimeter range using laser ultrasound.

Laser ultrasound is a widely used tool for industrial quality assurance when a contactless and fast method is required. In this work, we used a laboratory setup based on a confocal Fabry-Perot interferometer to examine weld seam models. The focus was placed on small samples with curved surfaces (small in the sense that the radius of curvature is comparable to the largest ultrasonic wavelength) and on efficient ways to detect the presence and volume of process pores, with the goal to transfer this method to industrial applications. In addition to this experimental method for investigating welds, a numerical method that models the experimental setup was implemented in MATLAB. For this purpose, first the thermal effects of the excitation process were taken into account by solving the thermal diffusion equation with an explicit scheme. Then, the elastodynamic equations were solved using the Elastodynamic Finite Integration Technique, taking into account the stresses induced by the excitation process. The B-Scans generated with this numerical model were compared with experimental B-Scans for simple test cases and good agreement was found. In a next step, the additional structures in the B-Scans resulting from air inclusions were identified and investigated with both methods using flat test specimens at first. Besides the direct echoes, structures from skimming surface waves and multiple reflections were visible. These additional structures are unwanted in defect reconstruction methods like the Synthetic Aperture Focusing Technique (SAFT) as they would lead to artifacts. In samples much larger than the largest ultrasound wavelength, however, these unwanted structures are still negligible in amplitude or can be well separated temporally, but for small samples this is no longer the case. As a result, reconstruction methods based on direct echoes like SAFT are difficult to apply. For many industrial applications, the reconstruction is not decisive at all, but only the knowledge of the total volume of process pores (TVPP). It is shown with both experimental and numerical methods, that this TVPP can be estimated from the variation in the B-Scans from various small weld seam models.

Photoacoustic computational ghost imaging.

Photoacoustic imaging with optical resolution usually requires a single-pixel raster scan. An alternative approach based on illumination with patterns obtained from a Hadamard matrix, measurement of the generated ultrasound wave with a single detector, followed by a reconstruction known from computational ghost imaging is demonstrated here. Since many pixels on the object are illuminated at the same time, thereby contributing to the recorded signal, this approach gives a better contrast-to-noise ratio compared to the raster scan, as demonstrated in a phantom experiment. Furthermore, exploiting the temporal information for depth-resolved imaging is possible. The proposed method will be beneficial in situations where the radiant exposure of a sample is limited due to either safety precautions or the properties of the available light source.

Physical Considerations for In Vitro ESWT Research Design.

In vitro investigations, which comprise the bulk of research efforts geared at identifying an underlying biomechanical mechanism for extracorporeal shock wave therapy (ESWT), are commonly hampered by inadequate descriptions of the underlying therapeutic acoustical pressure waves. We demonstrate the necessity of in-situ sound pressure measurements inside the treated samples considering the significant differences associated with available applicator technologies and cell containment. A statistical analysis of pulse-to-pulse variability in an electrohydraulic applicator yields a recommendation for a minimal pulse number of n = 300 for cell pallets and suspensions to achieve reproducible treatments. Non-linear absorption behavior of sample holders and boundary effects are shown for transient peak pressures and applied energies and may serve as a guide when in-situ measurements are not available or can be used as a controllable experimental design factor. For the use in microbiological investigations of ESWT we provide actionable identification of common problems in describing physical shockwave parameters and improving experimental setups by; (1) promoting in-situ sound field measurements, (2) statistical evaluation of applicator variability, and (3) extrapolation of treatment parameters based on focal and treatment volumes.

Conical ring array detector for large depth of field photoacoustic macroscopy.

Photoacoustic microscopy and macroscopy (PAM) using focused detector scanning are emerging imaging methods for biological tissue, providing high resolution and high sensitivity for structures with optical absorption contrast. However, achieving a constant lateral resolution over a large depth of field for deeply penetrating photoacoustic macroscopy is still a challenge. In this work, a detector design for scanning photoacoustic macroscopy is presented. Based on simulation results, a sensor array geometry is developed and fabricated that consists of concentric ring elements made of polyvinylidene fluoride (PVDF) film in a geometry that combines a centered planar ring with several inclined outer ring elements. The reconstruction algorithm, which uses dynamic focusing and coherence weighting, is explained and its capability to reduce artefacts occurring for single element conical sensors is demonstrated. Several phantoms are manufactured to evaluate the performance of the array in experimental measurements. The sensor array provides a constant axial and lateral resolution of 95 µm and 285 µm, respectively, over a depth of field of 20 mm. The depth of field corresponds approximately to the maximum imaging depth in biological tissue, estimated from the sensitivity of the array. With its ability to achieve the maximum resolution even with a very small scanning range, the array is believed to have applications in the imaging of limited regions of interest buried in biological tissue.

TGV-regularized inversion of the Radon transform for photoacoustic tomography.

We propose and study a reconstruction method for photoacoustic tomography (PAT) based on total generalized variation (TGV) regularization for the inversion of the slice-wise 2D-Radon transform in 3D. The latter problem occurs for recently-developed PAT imaging techniques with parallelized integrating ultrasound detection where projection data from various directions is sequentially acquired. As the imaging speed is presently limited to 20 seconds per 3D image, the reconstruction of temporally-resolved 3D sequences of, e.g., one heartbeat or breathing cycle, is very challenging and currently, the presence of motion artifacts in the reconstructions obstructs the applicability for biomedical research. In order to push these techniques forward towards real time, it thus becomes necessary to reconstruct from less measured data such as few-projection data and consequently, to employ sophisticated reconstruction methods in order to avoid typical artifacts. The proposed TGV-regularized Radon inversion is a variational method that is shown to be capable of such artifact-free inversion. It is validated by numerical simulations, compared to filtered back projection (FBP), and performance-tested on real data from phantom as well as in-vivo mouse experiments. The results indicate that a speed-up factor of four is possible without compromising reconstruction quality.

Quantitative comparison of frequency-domain and delay-and-sum optoacoustic image reconstruction including the effect of coherence factor weighting.

Image reconstruction in optoacoustic imaging is often based on a delay-and-sum (DAS) or a frequency domain (FD) algorithm. In this study, we performed a comprehensive comparison of these two algorithms together with coherence factor (CF) weighting using phantom and in-vivo mouse data obtained with optoacoustic microscopy. For this purpose we developed an FD based definition of the CF. Our results reveal the equivalence of DAS and FD, with and without CF weighting, in terms of spatial resolution and contrast-to-noise ratio (CNR) but highlight the clear advantage of FD in terms of computational cost, making it preferable for 3D reconstruction or real-time applications. An important additional result of this research is that, contradictory to previous studies, CF weighting does not lead to any improvement in lateral resolution.

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.

Modeling photoacoustic imaging with a scanning focused detector using Monte Carlo simulation of energy deposition.

Photoacoustic imaging using a focused, scanning detector in combination with a pulsed light source is a common technique to visualize light-absorbing structures in biological tissue. In the acoustic resolution mode, where the imaging resolution is given by the properties of the transducer, there are various challenges related to the choice of sensors and the optimization of the illumination. These are addressed by linking a Monte Carlo simulation of energy deposition to a time-domain model of acoustic propagation and detection. In this model, the spatial and electrical impulse responses of the focused transducer are combined with a model of acoustic attenuation in a single response matrix, which is used to calculate detector signals from a volumetric distribution of absorbed energy density. Using the radial symmetry of the detector, the calculation yields a single signal in less than a second on a standard personal computer. Various simulation results are shown, comparing different illumination geometries and demonstrating spectral imaging. Finally, simulation results and experimental images of an optically characterized phantom are compared, validating the accuracy of the model. The proposed method will facilitate the design of photoacoustic imaging devices and will be used as an accurate forward model for iterative reconstruction techniques.

Piezoelectric line detector array for photoacoustic tomography.

Photoacoustic tomography relies on a dense coverage of the surface surrounding the imaged object with ultrasound sensors in order to enable an accurate reconstruction. A curved arrangement of integrating line sensors is proposed that is able to acquire data for a linear projection image of the absorbed energy density distribution in the object. Upon rotation of the object relative to the array, three-dimensional (3D) images can be obtained. The proposed design is based on the cost-effective piezoelectric polymer film technology with 64 line shaped sensors arranged on a half-cylindrical surface. It is combined with an optical parametric oscillator for the near infrared as a source for laser pulses. Image reconstruction from recorded signals consists of two-dimensional (2D) back projection followed by an inverse Radon transform. The tomograph exhibits a spatial resolution on the order of 200 to 250 µm. In a phantom experiment, the steps from acquisition of a single, 2D projection image to a full 3D image are demonstrated. Finally, in vivo projection images of a human finger are shown, revealing the near real-time imaging capability of the device in 2D.

Combined photoacoustic, pulse-echo laser ultrasound, and speed-of-sound imaging using integrating optical detection.

A purely optical setup for the coregistration of photoacoustic (PA), ultrasound (US), and speed-of-sound (SOS) section images is presented. It extends a previously developed method for simultaneous PA and laser-US (LUS) pulse-echo imaging with a LUS transmission imaging setup providing two-dimensional (2-D) SOS maps. For transmission imaging, the sound waves traversing the investigated object are generated instantaneously by illuminating optically absorbing targets that are arranged at various distances in front of the sample. All signals are recorded by an optical beam which is part of a Mach­Zehnder interferometer that integrates the acoustic field along its path. Due to the cascaded arrangement of LUS sources, a single-recorded signal yields information for a projection of the SOS distribution. After collection of data from all directions, an inverse Radon transform is applied to this set of projections to obtain a 2-D SOS image. The setup is characterized and its performance is tested on phantom experiments. In addition to providing additional contrast, it is also shown that the resolution of the coregistered PA and LUS images can be improved by implementing the knowledge of the SOS distribution in the reconstruction.

High resolution three-dimensional photoacoutic tomography with CCD-camera based ultrasound detection.

A photoacoustic tomograph based on optical ultrasound detection is demonstrated, which is capable of high resolution real-time projection imaging and fast three-dimensional (3D) imaging. Snapshots of the pressure field outside the imaged object are taken at defined delay times after photoacoustic excitation by use of a charge coupled device (CCD) camera in combination with an optical phase contrast method. From the obtained wave patterns photoacoustic projection images are reconstructed using a back propagation Fourier domain reconstruction algorithm. Applying the inverse Radon transform to a set of projections recorded over a half rotation of the sample provides 3D photoacoustic tomography images in less than one minute with a resolution below 100 µm. The sensitivity of the device was experimentally determined to be 5.1 kPa over a projection length of 1 mm. In vivo images of the vasculature of a mouse demonstrate the potential of the developed method for biomedical applications.

Deblurring algorithms accounting for the finite detector size in photoacoustic tomography.

Most reconstruction algorithms for photoacoustic tomography, like back projection or time reversal, work ideally for point-like detectors. For real detectors, which integrate the pressure over their finite size, images reconstructed by these algorithms show some blurring. Iterative reconstruction algorithms using an imaging matrix can take the finite size of real detectors directly into account, but the numerical effort is significantly higher compared to the use of direct algorithms. For spherical or cylindrical detection surfaces, the blurring caused by a finite detector size is proportional to the distance from the rotation center (spin blur) and is equal to the detector size at the detection surface. In this work, we apply deconvolution algorithms to reduce this type of blurring on simulated and on experimental data. Two particular deconvolution methods are compared, which both utilize the fact that a representation of the blurred image in polar coordinates decouples pixels at different radii from the rotation center. Experimental data have been obtained with a flat, rectangular piezoelectric detector measuring signals around a plastisol cylinder containing various small photoacoustic sources with variable distance from the center. Both simulated and experimental results demonstrate a nearly complete elimination of spin blur.

Artifact removal in photoacoustic section imaging by combining an integrating cylindrical detector with model-based reconstruction.

Photoacoustic section imaging reveals optically absorbing structures within a thin slice of an object. It requires measuring acoustic waves excited by absorption of short laser pulses with a cylindrical acoustic lens detector rotating around the object. Owing to the finite detector size and its limited depth of focus, various artifacts arise, seen as distortions within the imaging slice and cross-talk from neighboring areas of the object. The presented solution aims at avoiding these artifacts by a special design of the sensor and by use of a model-based reconstruction algorithm that improves section images by incorporating information from neighboring sections. The integrating property of the cylindrical detector, which exceeds in direction of the cylinder axis the size of the imaged object, avoids the lateral blurring that normally results from the finite width of a small detector. Applying a maximum likelihood reconstruction method for the inversion of the imaging system matrix to the temporal pressure signals yields line projections of the initial energy distribution, from which section images are obtained by applying the inverse Radon transform. By using data from few sections, a significant reduction of artifacts related to the imperfections of the sensor is demonstrated both in simulations and in phantom experiments.

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.

Hybrid photoacoustic and ultrasound section imaging with optical ultrasound detection.

A setup is proposed that provides perfectly co-registered photoacoustic (PA) and ultrasound (US) section images. Photoacoustic and ultrasound backscatter signals are generated by laser pulses coming from the same laser system, the latter by absorption of some of the laser energy on an optically absorbing target near the imaged object. By measuring both signals with the same optical detector, which is focused into the selected section by use of a cylindrical acoustic mirror, the information for both images is acquired simultaneously. Co-registered PA and US images are obtained after applying the inverse Radon transform to the data, which are gathered while rotating the object relative to the detector. Phantom experiments demonstrate a resolution of 1.1 mm between the sections of both imaging modalities and a in-plane resolution of about 60 µm and 120 µm for the US and PA modes, respectively. The complementary contrast mechanisms of the two modalities are shown by images of a zebrafish.

Photoacoustic section imaging using an elliptical acoustic mirror and optical detection.

A method is proposed that utilizes the advantages of optical ultrasound detection in two-dimensional photoacoustic section imaging, combining an optical interferometer with an acoustic mirror. The concave mirror has the shape of an elliptical cylinder and concentrates the acoustic wave generated around one focal line in the other one, where an optical beam probes the temporal evolution of acoustic pressure. This yields line projections of the acoustic sources at distances corresponding to the time of flight, which, after rotating the sample about an axis perpendicular to the optical detector, allows reconstruction of a section using the inverse Radon transform. A resolution of 120 [micro sign]m within and 1.5 mm between the sections can be obtained with the setup. Compared to a bare optical probe beam, the signal-to-noise ratio (SNR) is seven times higher with the mirror. Furthermore, the imaging system is tested on a biological sample.

Photoacoustic section imaging with an integrating cylindrical detector.

A piezoelectric detector with a cylindrical shape is investigated for photoacoustic section imaging. Images are acquired by rotating a sample in front of the cylindrical detector. With its length exceeding the size of the imaging object, it works as an integrating sensor and therefore allows reconstructing section images with the inverse Radon transform. Prior to the reconstruction the Abel transform is applied to the measured signals to improve the accuracy of the image. A resolution of about 100 µm within a section and of 500 µm between sections is obtained. Additionally, a series of images of a zebra fish is shown.

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.

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.

Downstream Fabry-Perot interferometer for acoustic wave monitoring in photoacoustic tomography.

An optical detection setup consisting of a focused laser beam fed into a downstream Fabry-Perot interferometer (FPI) for demodulation of acoustically generated optical phase variations is investigated for its applicability in photoacoustic tomography. The device measures the time derivative of acoustic signals integrated along the beam. Compared to a setup where the detection beam is part of a Mach-Zehnder interferometer, the signal-to-noise ratio of the FPI is lower, but the image quality of the two devices is similar. Using the FPI in a photoacoustic tomograph allows scanning the probe beam around the imaging object without moving the latter.