Single-pixel imaging of dynamic objects using multi-frame motion estimation.

2021 Scientific Reports

Abstract
Single-pixel imaging (SPI) enables the visualization of objects with a single detector by using a sequence of spatially modulated illumination patterns. For natural images, the number of illumination patterns may be smaller than the number of pixels when compressed-sensing algorithms are used. Nonetheless, the sequential nature of the SPI measurement requires that the object remains static until the signals from all the required patterns have been collected. In this paper, we present a new approach to SPI that enables imaging scenarios in which the imaged object, or parts thereof, moves within the imaging plane during data acquisition. Our algorithms estimate the motion direction from inter-frame cross-correlations and incorporate it in the reconstruction model. Moreover, when the illumination pattern is cyclic, the motion may be estimated directly from the raw data, further increasing the numerical efficiency of the algorithm. A demonstration of our approach is presented for both numerically simulated and measured data.

[Read more…]

Schematic of the two proposed algorithms (a) Flow of the algorithm in case of global motion. (b) Flow of algorithm in case of local motion..

Sagi Monin, Evgeny Hahamovich & Amir Rosenthal

2021 Scientific Reports https://doi.org/10.1038/s41598-021-83810-z

Grüneisen-relaxation photoacoustic microscopy at 1.7 µm and its application in lipid imaging.

Schematic of the GR-PAM experimental system

© 2020 Optical Society of America

Abstract
We report the first, to the best of our knowledge, demonstration of Grüneisen relaxation photoacoustic microscopy (GR-PAM) of lipid-rich tissue imaging at the 1.7 µm band, implemented with a high-energy thulium-doped fiber laser and a fiber-based delay line. GR-PAM enhances the image contrast by intensifying the region of strong absorbers and suppressing out-of-focus signals. Using GR-PAM to image swine-adipose tissue at 1725 nm, an 8.26-fold contrast enhancement is achieved in comparison to conventional PAM. GR-PAM at the 1.7 µm band is expected to be a useful tool for label-free high-resolution imaging of lipid-rich tissue, such as atherosclerotic plaque and nerves.

[Read more…]

Swine muscular tissue images

Fig. Swine muscular tissue images under (a)-(c) OR-PAM and (d)-(f) GR-PAM at (a), (d) 1700; (b), (e) 1725; and (c), (f) 1750 nm. (g) Reference confocal microscopy image at the same region.

Jiawei Shi, Can Li, Huade Mao, Yuxuan Ren, Zhi-Chao Luo, Amir Rosenthal, and Kenneth K. Y. Wong

© 2020 Optical Society of America   Vol. 45, Issue 12, pp. 3268-3271 (2020)

Increased SNR in acousto-optic imaging via coded ultrasound transmission.

© 2020 Optical Society of America

Abstract

Acousto-optic imaging (AOI) is a non-invasive method that uses acoustic modulation to map the light fluence inside biological tissue. In many AOI implementations, ultrasound pulses are used in a time-gated measurement to perform depth-resolved imaging without the need for mechanical scanning. However, to achieve high axial resolution, it is required that ultrasound pulses with few cycles are used, limiting the modulation strength. In this Letter, we develop a new approach to pulse-based AOI in which coded ultrasound transmission is used. In coded-transmission AOI (CT-AOI), one may achieve an axial resolution that corresponds to a single cycle, but with a signal-to-noise ratio (SNR) that scales as the square root of the number of cycles. Using CT-AOI with 79 cycles, we experimentally demonstrate over four-fold increase in SNR in comparison to a single-cycle AOI scheme.

One of the fundamental limitations of optical imaging of biological tissue is light scattering due to optical heterogeneity. At depths exceeding several transport lengths, scattering leads to the diffusion of light, which severely limits the imaging resolution that may be achieved [1]. Additionally, optical imaging with diffused light often requires solving nonlinear optimization problems in order to map tissue parameters.

Acousto-optic imaging (AOI) is a hybrid approach that overcomes the limitations of light diffusion by using acoustic modulation [2]. Conventionally, AOI is performed by illuminating the tissue with a highly coherent continuous-wave (CW) laser and using ultrasound to locally modulate the phase of the laser light inside the tissue. In AOI, the ultrasound-induced phase modulation is a result of two mechanisms [3]: pressure-induced modulation of the refractive index and periodic movement of the optical scatterers. When the coherence length of the laser is sufficiently long, the local ultrasound-induced phase modulation inside the tissue is translated into an intensity modulation of the speckle pattern on the tissue boundary. Thus, by measuring the modulation depth of the speckle on the tissue boundary, it is possible to quantify the light fluence within the tissue at the positions in which the acoustic modulation was performed [4].

AOI is capable of identifying both highly absorbing and highly scattering structures through their effect on the light fluence [5], facilitating applications such as early assessment of osteoporosis [6]. Additionally, AOI can provide information on blood flow in the acoustically modulated regions through analysis of the spectral broadening of the speckle modulation [7]. While in most applications AOI is used as an independent technique for assessing tissue parameters, it may also be used as a complimentary technique to optoacoustic tomography (OAT). In previous works [8,9], it has been shown that the information provided by AOI can remove the bias in OAT images due to light attenuation, thus enabling OAT-image quantification.

[Read more…]

1D profile of the modulated light along the ultrasound propagation path for single-pulse.

Fig. 1D profile of the modulated light along the ultrasound propagation path for single-pulse AOI (blue) and CT-AOI (red) with corresponding FWMH values of 4.32 mm and 4.07 mm, respectively.

 

Ahiad Levi, Sagi Monin, Evgeny Hahamovich, Aner Lev, Bruno G. Sfez, and Amir Rosenthal

© 2020 Optical Society of America Vol. 45, Issue 10, pp. 2858-2861 (2020)

Ultrasound Detection Arrays Via Coded Hadamard Apertures.

Experimental setup.

© 2020 IEEE

Abstract
In the medical fields, ultrasound detection is often performed with piezoelectric arrays that enable one to simultaneously map the acoustic fields at several positions. In this work, we develop a novel method for transforming a single-element ultrasound detector into an effective detectionarray by spatially filtering the incoming acoustic fields using a binary acoustic mask coded with cyclic Hadamard patterns. By scanning the mask in front of the detector, we obtain a multiplexed measurement dataset from which a map of the acoustic field is analyticallyconstructed. We experimentally demonstrate our method by transforming a single-element ultrasound detector into 1D arrays with up to 59 elements.

[Read more…]

Fig. Normalized waveform shape comprasion of the strongest signal for a single vs. multiplexed aperture detection. The plots are for (a) an aperture with 1.5 mm diameter and 31 elements and (b) an aperture with 1 mm diameter and 59 elements. 

E. Hahamovich and A. Rosenthal, “Ultrasound Detection Arrays Via Coded Hadamard Apertures,” in IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, doi: 10.1109/TUFFC.2020.2993583.

© 2020 IEEE

Burst-mode pulse interferometry for enabling low-noise multi-channel optical detection of ultrasound.

A schematic drawing of BM-PI.

© 2022 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement© 2019 Optics Express, OSA publishing.

Abstract
Ultrasound detection via optical resonators can achieve high levels of miniaturization and sensitivity as compared to piezoelectric detectors, but its scale-up from a single detector to an array is highly challenging. While the use of wideband sources may enable parallel interrogation of multiple resonators, it comes at the cost of reduction in the optical power, and ultimately in sensitivity, per channel. In this work we have developed a new interferometric approach to overcome this signal loss by using high-power bursts that are synchronized with the time window in which ultrasound detection is performed. Each burst is composed of a train of low-noise optical pulses which are sufficiently wideband to interrogate an array of resonators with non-overlapping spectra. We demonstrate our method, termed burst-mode pulse interferometry, for interrogating a single resonator in which the optical power was reduced to emulate the power loss per channel that occurs in parallel interrogation of 20 to 200 resonators. The use of bursts has led to up 25-fold improvement in sensitivity without affecting the shape of the acoustic signals, potentially enabling parallel low-noise interrogation of resonator arrays with a single source.

[Read more…]

The optical signal at the CM (blue) and BM (red)

Fig. The optical signal at the CM (blue) and BM (red) of the CRF output in the transition from the unlocked stated, in which the CRF blocks the optical pulses, to the locked state in which the CRF transmits the optical pulses and blocks only the ASE. Each of the spikes in the locked BM channel represents a burst with a width of 10 µs, shown in detail in Fig. 2(d), where the burst repletion rate was 8 kHz.

Oleg Volodarsky, Yoav Hazan, Michael Nagli, and Amir Rosenthal.
Optics Express, Volume 30, Issue 6, Pages 8959-8973
https://opg.optica.org/oe/abstract.cfm?URI=oe-30-6-8959

 

Simultaneous multi-channel ultrasound detection via phase modulated pulse interferometry.

The acoustic setup used to test the performance of PM-PI.

© 2019 Optics Express, OSA publishing.

Abstract
In optical detection of ultrasound, resonators with high Q-factors are often used to maximize sensitivity. However, in order to perform parallel interrogation, conventional interferometric techniques require an overlap between the spectra of all the resonators, which is difficult to achieve with high Q-factor resonators. In this paper, a new method is developed for parallel interrogation of optical resonators with non-overlapping spectra. The method is based on a phase-modulation scheme for pulse interferometry (PM-PI) and requires only a single photodetector and sampling channel per ultrasound detector. Using PM-PI, parallel ultrasound detection is demonstrated with four high Q-factor resonators.

[Read more…]

A schematic drawing of the PM-PI.

Fig. A schematic drawing of the PM-PI system used in this work to interrogate 4 resonators, implemented with Ï€-phase shifted fiber Bragg gratings (Ï€-FBGs). A wideband pulse laser with band-pass filters (BPFs) and an erbium-doped fiber amplifier (EDFA) create a source with a high spectral power density and sufficient bandwidth to cover the spectra of all the resonators. The modulation unit is an unbalanced Mach-Zehnder interferometer (MZI), composed of optical fiber couplers (FC) and a phase modulator (PM). The input phase signal to the PM, shown in the top-right plot, alternates between two values with a difference of ðœ‹/2. For each phase value, the pulses interfere differently at the output of each resonator depending on the phase difference in the MZI for the specific resonance wavelength of that resonator. The bottom-right plot, shows a typical voltage signal measured for one of the resonators, which alternates between two states that correspond to the two phase values. As the bottom-right plot shows, in the current implementation, the duration of each phase value delivered to the PM corresponded to 5 laser pulses. We note that the limited bandwidth of our measurement did not allow full separation between the pulses in the bottom-right plot.

Yoav Hazan and Amir Rosenthal. 
Optics Express Vol. 27, Issue 20, pp. 28844-28854 (2019) https://doi.org/10.1364/OE.27.028844

 

Enhanced Sensitivity of Silicon-Photonics-Based Ultrasound Detection via BCB Coating.

Sio and Bcb over clading

© 2019 IEEE Photonics Journal

Impact Statement:
This paper gives a solution to one of the fundamental limitations of silicon-photonics based ultrasound detectors: the low photo-elastic response of silicon and silica. By using a BCB over-cladding, 5-fold increase in the acoustic sensitivity is achieved. We additionally provide a detailed analysis of the sensing mechanism, quantifying the different effects that contribute to the enhanced sensitivity.
Abstract:
Ultrasound detection via silicon waveguides relies on the ability of acoustic waves to modulate the effective refractive index of the guided modes. However, the low photo-elastic response of silicon and silica limits the sensitivity of conventional silicon-on-insulator sensors, in which the silicon core is surrounded by a silica cladding. In this paper, we demonstrate that the sensitivity of silicon waveguides to ultrasound may be significantly enhanced by replacing the silica over-cladding with bisbenzocyclobutene (BCB)-a transparent polymer with a high photo-elastic coefficient. In our experimental study, the response to ultrasound, in terms of the induced modulation in the effective refractive index, achieved for a BCB-coated silicon waveguide with TM polarization was comparable to values previously reported for polymer waveguides and an order of magnitude higher than the response achieved by an optical fiber. In addition, in our study, the susceptibility of the sensors to surface acoustic waves and reverberations was reduced for both TE and TM modes when the BCB over-cladding was used.

[Read more…]

setup schematic

Fig. – The measurement setup. For each of the polarizations, a Mach-Zehnder interferometer was constructed, where in each of the interferometer arms a chip with a different over-cladding material (BCB or silica) was connected. An ultrasound transducer was used to generate acoustic waves that impinged on only one of the chips, which were separated by more than 10 cm.

Resmi Ravi Kumar, Evgeny Hahamovich, Shai Tsesses, Yoav Hazan, Assaf Grinberg, Amir Rosenthal. IEEE Photonics Journal ( Volume: 11 , Issue: 3 , June 2019 )

Noise reduction in resonator-based ultrasound sensors by using a CW laser and phase detection.

© 2019 Optical Society of America

The detection of ultrasound via optical resonators is conventionally performed by tuning a continuous-wave (CW) laser to the linear slope of the resonance and monitoring the intensity modulation at the resonator output. While intensity monitoring offers the advantage of simplicity, its sensitivity is often limited by the frequency noise of the CW laser. In this work, we develop an alternative CW technique that can significantly reduce measurement noise by monitoring variations in the phase, rather than intensity, at the resonator output. In our current implementation, which is based on a balanced Mach–Zehnder interferometer for phase detection, we demonstrate a 24-fold increase in the signal-to-noise ratio of the detected ultrasound signal over the conventional, intensity-monitoring approach.[Read more…]

Signals detected using the IM and PM scheme

Fig. – Signals detected using the IM and PM interrogation methods with an optical power of 0.1 mW. In the PM measurement, the OPD was set to zero to maximize the SNR.

Lucas Riobó, Yoav Hazan, Francisco Veiras, María Garea, Patricio Sorichetti, and Amir Rosenthal.  Optics Letters Vol. 44, Issue 11, pp. 2677-2680 (2019) •https://doi.org/10.1364/OL.44.002677
© 2019 Optical Society of America

The Impulse Response of Negatively Focused Spherical Ultrasound Detectors and its Effect on Tomographic Optoacoustic Reconstruction.

The Impulse Response

source: © 2019  IEEE Transactions on Medical Imaging. 

In optoacoustic tomography, negatively focused detectors have been shown to improve the tangential image resolution without sacrificing sensitivity. Since no exact inversion formulae exist for optoacoustic image reconstruction with negatively focused detectors, image reconstruction in such cases is based on using the virtual-detector approximation, in which it is assumed that the response of the negatively focused detector is identical, up to a constant time delay, to that of a point-like detector positioned in the detector’s center of curvature. In this work, we analyze the response of negatively focused spherical ultrasound detectors in three dimensions and demonstrate how their properties affect the optoacoustic reconstruction. Our analysis sheds new light on commonly reported experimental reconstruction artifacts in optoacoustic systems that employ negatively focused detectors. Based on our analysis, we introduce a simple correction to the virtual-detector approximation that significantly enhances image contrast and reduces artifacts.  [Read more…]

The Impulse Response

Fig. (a) The geometry of the negative acoustic lens studied with full
acoustic simulations. The speed of sound of the surrounding medium and lens
material were 1500 m/s and 2757 m/s. respectively. (b) The detected acoustic
signals obtained when the lens was acoustically matched to the surrounding
medium (solid-blue curve) and when its acoustic impedance was 1.xx times that
of the surrounding medium (dashed-red curve), leading to internal reflections
in the lens structure. The reflection from the detection surface was 50% of the
pressure signal. 

Gilad Drozdov, Ahiad Levi, Amir Rosenthal.
IEEE Transactions on Medical Imaging.DOI: 10.1109/TMI.2019.2897588

 

Algebraic determination of back-projection operators for optoacoustic tomography

source: © 2018 Optical Society of America

The simplicity and computational efficiency of back-projection formulae have made them a popular choice in optoacoustic tomography. Nonetheless, exact back-projection formulae exist for only a small set of tomographic problems. This limitation is overcome by algebraic algorithms, but at the cost of higher numerical complexity. In this paper, we present a generic algebraic framework for calculating back-projection operators in optoacoustic tomography. We demonstrate our approach in a two-dimensional optoacoustic-tomography example and show that once the algebraic back-projection operator has been found, it achieves a comparable run time to that of the conventional back-projection algorithm, but with the superior image quality of algebraic methods.[Read More…]

Fig. 1 (a) The grid of the image and detector locations used for calculating the model matrix?. The image is divided into ??×?? square pixels with a pixel area of ???? and the acoustic signals are sampled at ?? positions over a line with a distance of?? between them. (b) The image grid on which the projection operator ? is calculated. Here, only a single back-projection is calculated, and the number of pixels in the x directions is increased to ??=??+??−1.

Amir Rosenthal, “Algebraic determination of back-projection operators for optoacoustic tomography,” Biomed. Opt. Express 9, 5173-5193 (2018)