Guenther Paltauf

Iterative reconstruction algorithm for optoacoustic imaging

Optoacoustic imaging is based on the generation of thermoelastic stress waves by heating an object in an optically heterogeneous medium with a short laser pulse. The stress waves contain information about the distribution of structures with preferential optical absorption. Detection of the waves with an array of broadband ultrasound detectors at the surface of the medium and applying a backprojection algorithm is used to create a map of absorbed energy inside the medium. With conventional reconstruction methods a large number of detector elements and filtering of the signals are necessary to reduce backprojection artifacts. As an alternative this study proposes an iterative procedure. The algorithm is designed to minimize the error between measured signals and signals calculated from the reconstructed image. In experiments using broadband optical ultrasound detectors and in simulations the algorithm was used to obtain three-dimensional images of multiple optoacoustic sources. With signals from a planar array of 3x3 detector elements a significant improvement was observed after about 10 iterations compared to the simple radial backprojection. Compared to conventional methods using filtered backprojection, the iterative method is computationally more intensive but requires less time and instrumentation for signal acquisition.

Photoacoustic tomography using a Mach-Zehnder interferometer as an acoustic line detector

A three-dimensional photoacoustic imaging method is presented that uses a Mach-Zehnder interferometer for measurement of acoustic waves generated in an object by irradiation with short laser pulses. The signals acquired with the interferometer correspond to line integrals over the acoustic wave field. An algorithm for reconstruction of a three-dimensional image from such signals measured at multiple positions around the object is shown that is a combination of a frequency-domain technique and the inverse Radon transform. From images of a small source scanning across the interferometer beam it is estimated that the spatial resolution of the imaging system is in the range of 100 to about 300 mum, depending on the interferometer beam width and the size of the aperture formed by the scan length divided by the source-detector distance. By taking an image of a phantom it could be shown that the imaging system in its present configuration is capable of producing three-dimensional images of objects with an overall size in the range of several millimeters to centimeters. Strategies are proposed how the technique can be scaled for imaging of smaller objects with higher resolution.

THERMOACOUSTIC TOMOGRAPHY AND THE CIRCULAR RADON TRANSFORM: EXACT INVERSION FORMULA

Thermoacoustic computed tomography (TACT) is an emerging hybrid imaging method for non-invasive medical diagnosis and fully three-dimensional visualization of biological probes. Within this modality electromagnetic illumination is used to induce acoustic waves inside an object of interest. In this paper, we assume that a cylindrical array of line detectors is used to record the acoustical data. This leads to the mathematical problem of inverting the circular Radon transform. The circular Radon transform arises in several other up-to-date imaging modalities, such as RADAR imaging or ultrasound tomography. In this paper we prove a novel stable formula for recovering a planar function from its circular Radon transform. We apply this formula to obtain an exact three-dimensional imaging algorithm for TACT. Numerical reconstructions from real and synthetic data demonstrate the potential and robustness of our algorithm.

Stealth ryanodine‐sensitive Ca2+ release contributes to activity of capacitative Ca2+ entry and nitric oxide synthase in bovine endothelial cells

1 The involvement of ryanodine‐sensitive Ca2+ release (RsCR) in bradykinin (Bk)‐induced Ca2+ release, capacitative Ca2+ entry (CCE) and nitric oxide synthase (NOS) activation was assessed in freshly isolated bovine coronary artery endothelial cells. 2 Using deconvolution microscopy fura‐2 was found throughout the whole cytosol, while the cell membrane impermeable dye FFP‐18 was exclusively in the cell membrane. Thus, perinuclear ([Ca2+]pn) and subplasmalemmal Ca2+ concentration ([Ca2+]sp) were monitored using fura‐2 and FFP‐18. 3 Inhibition of Na+−Ca2+ exchange by lowering extracellular Na+ concentration augmented the Bk‐induced [Ca2+]pn signal in Ca2+‐free solution. This effect was abolished when RsCR was prevented with 25 μmmu;mol l−1 ryanodine, while inhibition of RsCR had no effect on Bk‐induced increase in [Ca2+]pn without inhibition of Na+−Ca2+ exchange. 4 Initiating RsCR by 200 nmol l−1 ryanodine increased [Ca2+]sp, while [Ca2+]pn remained constant. However, when Na+−Ca2+ exchange was prevented, ryanodine was also able to elevate [Ca2+]pn. 5 Blockage of RsCR diminished Ca2+ extrusion in response to stimulation with Bk in normal Na+‐containing solution. 6 Inhibition of RsCR blunted Bk‐activated CCE, while inhibition of Na+−Ca2+ exchange during stimulation enhanced CCE. 7 Although direct activation of RsCR failed to activate NOS, inhibition of RsCR diminished the effect of ATP and Bk on NOS, while the effect of thapsigargin remained unchanged. 8 These data suggest that during stimulation subplasmalemmal RsCR occurs, which contributes to the activities of CCE and NOS. Thus, the function of the subplasmalemmal Ca2+ control unit must be extended as a regulator for CCE and NOS.

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.

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.

Study of different ablation models by use of high-speed-sampling photography

In our study we investigated the ablation characteristics of an aqueous dye solution with a defined absorption coefficient, irradiated by short (8 ns) and long (100 microsecond(s) ) pulses from a Nd:YAG laser (wavelength: 1064 nm). The experimental technique was schlieren photography with a second Nd:YAG laser at 532 nm as a light source and with a variable delay between the two laser pulses. With a special arrangement of the laser beams and the sample effects below and above the surface of the liquid could be simultaneously observed. We could distinguish three ablation mechanisms, depending on the pulse duration and the incident fluence. With short pulses and a fluence below the vaporization threshold the tensile pulse from the bipolar thermoelastic wave, propagating from the liquid-air interface into the sample, caused rupture and spallation of the liquid. At fluences generating a surface temperature in excess of 100 degree(s)C the short pulses caused explosive vaporization, characterized by shock wave emission both in air and in liquid. At the same fluence the long pulses caused slow vaporization, meaning that vapor and liquid ejection started during the laser pulse and was less violent than with the 8 ns pulses.

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.

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.

Design and testing of an endoscopic photoacoustic probe for determination of treatment depth after photodynamic therapy

An endoscopic photoacoustic probe is designed and tested for use in PDT treatment of esophageal cancer. The probe, measuring less than 2.5 mm in diameter, was designed to fit within the lumen of an endoscope that will be inserted into an esophagus after PDT. PDT treatment results in a blanched, necrotic layer of cancerous tissue over a healthy, deeper layer of perfused tissue. The photoacoustic probe was designed to use acoustic propagation time to determine the thickness of the blanched surface of the esophagus, which corresponds to treatment depth. A side-firing 600 micrometers fiber delivered 532 nm laser light to induce acoustic waves in the perfused layer of the esophagus beneath the blanched (treated) layer. A PVDF transducer detected the induced acoustic waves and transmitted the signal to an oscilloscope. The probe was tested on clear and turbid tissue phantom layers over an optically absorbing dye solution.