Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Microscopic OCT imaging with focus extension by ultrahigh-speed acousto-optic tunable lens and stroboscopic illumination

Open Access Open Access

Abstract

We develop high-resolution optical coherence tomography (OCT) system with high-speed acousto-optic tunable lens. Stroboscopic pulsed illumination is used for the first time to perform time-resolved OCT imaging with acousto-optic tunable focusing. The operation of ultrahigh-speed tunable acousto-optic lens is demonstrated theoretically and experimentally. Focal position tuning at MHz frequency range is experimentally shown in the imaging system leading to OCT images with extended depth of focus. Imaging with active optical elements is helpful for improvement of photon collection efficiency, depth of focus and enhancement of the image quality.

© 2014 Optical Society of America

1. Introduction

Optical coherence tomography (OCT) and its microscopic version (optical coherence microscopy, OCM) form the modalities that enable acquisition of cross-sectional and volumetric data sets and reveal microstructural information on the objects [1, 2]. Different generations of OCT instruments were developed to study objects with depth ranges spanning from millimeters to few centimeters [3–6]. Unique features of the OCT technique enabled its numerous applications in biology, medicine and industry [7, 8].

Resolution is the essential factor influencing the image quality in any imaging technique. Axial resolution in OCT is governed by the wavenumber range (light source bandwidth) recorded by the detection unit. Transverse resolution is determined by the beam spot size incident on the object, and it depends on the optical design of the sample (object) arm. Another important factor influencing the image quality in OCT / OCM is the depth of focus (DOF). In a standard configuration, the beam waist should not change considerably over the entire sampled axial range. In other words, the scanning beam Rayleigh range 2zR (depth of focus) should be at least of the order of the depth range zmax defined by the detection scheme. As a consequence, high-resolution OCT / OCM systems can effectively image relatively shallow structures since tight light focusing makes collection the light from the depths beyond DOF inefficient.

Novel applications of OCT are also characterized by long depth ranges and consequently have special requirements to perform optimized imaging [5]. Unlike OCM, long range OCT utilizes weakly focused beams where the long Rayleigh range is achieved at the expense of transverse resolution. Therefore, both high-resolution OCM and long-range OCT will benefit from novel approaches of depth of focus extension.

Several methods have been implemented to optimize the performance of axial imaging. The easiest solution involves mechanical movement of the sample with respect to the objective lens (and vice versa) [9, 10]. More sophisticated approaches are based on light beam engineering techniques and include utilization of: multi-beam scanning [11, 12], microelectromechanical system technology [13], non-diffractive beams [14–17], adaptive focus [18, 19] etc. Improved transverse resolution throughout the imaging depth can be also achieved with dynamic focusing using tunable lenses. The technology providing variable focus is based on the electro-optical, electromechanical, thermo-optical or acousto-mechanical effects [20–28]. Recently, computational imaging techniques have been developed for depth of focus extension [29–32].

Laser beam control can be also obtained by the interaction of light with ultrasound [33–37]. The acoustic wave travelling in the medium causes modulation of its optical properties. Consequently, when the light propagates through such a medium, one can observe spatial, temporal and spectral modulations of light. Therefore, the acousto-optic phenomenon has been applied in many photonic technologies involving modulation, beam steering (deflection), filtering, frequency shifting and spectrometry [38].

The aim of this study is to apply high-resolution OCT / OCM system with high-speed acousto-optic tunable lens in biomedical imaging. For the first time we integrated stroboscopic (pulsed) illumination and active optical element in the OCT system to perform microscopic OCT imaging with dynamic focusing. This approach enabled obtaining microscopic OCT images with extended focus. We also characterized the operation of the designed acousto-optic cell theoretically and experimentally, and we assessed the impact of acousto-optic dynamic focusing on the OCT performance.

2. Theoretical description of tunable acousto-optic cell operation

The acousto-optic device is always composed of the following elements: (i) acoustic (ultrasonic) wave propagating in (ii) the acoustic medium, and (iii) light (electromagnetic) wave. The ultrasound is launched into the medium by a piezoelectric transducer; in our case the transducer has the shape of a hollow cylinder of the length L; Fig. 1(a). The cylinder is closed at both faces with glass windows and filled with water (acoustic medium), thus forming the acousto-optic cell (AOC). Both inner and outer surfaces of the cell are covered with electrode material (e.g. silver). Therefore, wall thickness vibrations occur when the voltage is applied to both surfaces, which results in the spatio-temporal modulation of the density of acoustic medium; Fig. 1(b).

 figure: Fig. 1

Fig. 1 Geometry of acousto-optic interaction. (a) Scheme of the acousto-optic cell. The cell is composed of cylindrical transducer of the length L. Light beam passes along the axis of the transducer. (b) The numerical approach assumes dividing the acoustic medium into M layers, each being regarded as pure phase grating. (c) Vibration of the transducer generates standing ultrasonic wave in which concentration of ultrasound is observed at the axis (r = 0). Gaussian distribution of electric field amplitude and the distribution of acoustic pressure. (d) Light beam propagation in a single layer.

Download Full Size | PDF

Accordingly, specific distribution of the acoustic pressure due to compressions and rarefactions of the medium particles is generated inside the cell. In the vibrating cylinder, the axisymmetric ultrasonic field p(r, t) is described by the Bessel function of the first kind of zero order:

p(r,t)=δpJ0(Kr)cos(Ωt),
where δp is the amplitude of the acoustic pressure, r=(x2+y2)1/2, K and Ω are wavenumber and angular frequency of the ultrasound, respectively. Equation (1) indicates that the highest amplitude can be found at the axis of the transducer (focusing of ultrasound), as also shown in Fig. 1(c). Linear case of the acousto-optic effect assumes that the pressure wave inside the medium causes proportional modulation of the optical properties of the medium. This phenomenon can be expressed in terms of refractive index of the medium n(r, t):
n(r,t)=n0+Δn(r,t)=n0+δnJ0(Kr)cos(Ωt),
where n0 is the refractive index of the non-disturbed medium (no ultrasound), and δn is amplitude of refractive index modulation.

The light beam propagates through the optically transparent medium of the AOC. Assuming that the collimated Gaussian beam is incident on the ultrasound, the optical field can be therefore described by:

E(r)=E0e(r/w)2,
where 2w is the diameter of the beam and E0 is the amplitude of the electric field.

Determination of both amplitude and phase of the light field at the exit plane of the AOC (z = L) is a key issue in the theoretical description of the acousto-optic effect. From the theoretical point of view, the problem is actually related to the light propagation in the inhomogeneous medium. Various analytical and numerical solutions have been developed to describe light diffraction on ultrasound [39]. However, we use successive diffraction model (SDM), which is a beam propagation method (BPM) modified to the case of inhomogeneous medium [40, 41].

In SDM, specific approach must be taken to account for optical inhomogeneity of the medium, i.e. the entire interaction length L is divided into several layers, each having the thickness Δz = L/M, where M is the number of layers; Fig. 1(b). The number of layers should be high enough to treat each layer as a pure phase grating. It has been proved that this requirement is fulfilled when for each layer [42]:

Qv0.2,
where QK2Δz/(kn0)is the Klein-Cook parameter, and vkδnΔz is the Raman-Nath parameter. The prime symbol indicates that those parameters are related to the interaction in a separate layer. The light is propagated from the entrance to exit plane of the layer of in Fourier domain (spatial frequency space):
E˜(fx,fy,z+Δz)=E˜in(fx,fy,z)exp[iΔz(kn0)2(2π)2(fx2+fy2)],
where fx and fy are the components of spatial frequency in x and y direction, E˜(fx,fy,z)=FT{E(x,y,z)} denotes the 2-D Fourier transform of signal distribution E. The evanescent waves have been neglected in Eq. (5). When the light propagation is calculated, the signal is Fourier-transformed back to Cartesian (x, y) space. Then we take into account that the layer acts as the phase object:
Eout(x,y,z+Δz)=E(x,y,z+Δz)exp[ikΔn(r,t)].
The multiplier on the right hand side of Eq. (6) refers to phase factor due to the presence of the acoustic field. Later on, the optical field in the exit plane of the layer becomes the input field of the next layer. The above-mentioned steps are repeated until the exit plane (z = L) of the ultrasound field is achieved.

The same technique can be used to simulate beam propagation through the optical elements behind the AOC.

3. Experimental set-up

The experimental configurations used in this study are shown in Fig. 2. The tunable acousto-optic lens (AOL) is composed of two elements: active AOC (vibrating shell) and passive offset lens OL of the focal length fOL. The set-up to characterize the AOL is demonstrated in Fig. 2(a). The light was emitted by a femtosecond Ti:sapphire laser (Fusion, Femtolasers GmbH, Austria). The light intensity was later modulated using electro-optic modulator (LM0202; LINOS, Qioptiq Photonics GmbH, Germany). The light was then directed to the AOC through fiber coupler. Collimated Gaussian light beam of the width 2w = 800 μm was incident on the AOC. The piezoelectric transducer (Physik Instrumente GmbH) of the length L = 50 mm and wall thickness 2 mm (outer diameter 30 mm, inner diameter 26 mm) that was closed at both faces with optical windows and filled with water to provide effective ultrasound coupling. The cell was excited by a cosine signal of frequency of 1.042 MHz from dual-channel function generator (TGA 12102; Thurlby Thandar Instruments Ltd, UK).

 figure: Fig. 2

Fig. 2 Experimental set-up. (a) Configuration for experimental verification of the operation of acousto-optic tunable lens. (b) Spectral / Fourier-domain OCT instrument with tunable lens. EOM – electro-optic intensity modulator, FC – fiber coupler, AOC – acousto-optic cylindrical cell, OL – offset lens, BP – beam profiler, FG – function generator, DDG – digital delay line, PD – high-speed photodetector, OSC – oscilloscope, DC – dispersion compensation, LSC – line-scan camera, COMP – computer. (c) Explanation of continuous and stroboscopic illumination schemes. Application of EOM and digital delay line allows for control of phase difference between light pulses and ultrasound as well as to adjust the width of light pulses.

Download Full Size | PDF

It is important to note that the mutual relation between light beam width and ultrasonic wave frequency should be properly chosen in order to obtain the desired effect of acousto-optic interaction. The width of the light beam (2w = 800 μm) incident on acoustic field should not exceed the first nodal diameter of standing cylindrical ultrasonic wave (1.12 mm). Otherwise, light diffraction would appear that would manifest in separated diffraction orders (rings). This relation is presented in Fig. 1(c).

The pulsed signal of the same frequency and variable duty cycle from the second channel of the generator was used to drive the electro-optic modulator (EOM). We excited the EOM by the signal with either 100% or 5.5% (pulses of the width of 30 ns) duty cycle. Additionally, we utilized digital delay generator (DG 645, Stanford Research Systems Inc., USA) to control the time delay τ between signal driving the electro-optic modulator and the transducer; Fig. 2(a). This experimental solution enabled stroboscopic illumination synchronized with ultrasound at variable phase of the ultrasonic wave, and made the studies of highly dynamical systems such as the AOL easier (time-resolved measurements). The illumination schemes used in this study are presented in Fig. 2(c).

Other port of the fiber coupler was used to monitor the light intensity illuminating the AOL. Signal from the photodetector (PDB-110C; Thorlabs Inc., USA) and both signals from the function generator were visualized by the oscilloscope. We placed an objective (offset) lens OL (either + 40 mm lens or 20x microscope objective) behind the AOC to transform the light to the far field. Both AOC and OL form the acousto-optic tunable lens (AOL). Light intensity distribution behind the AOL was measured using a beam profiler (WinCamD; DataRay Inc., USA) on a motorized translation stage.

The fiber-optic OCT system equipped with high-speed tunable lens is demonstrated in Fig. 2(b). The light delivery subsystem similar to that in the previous case was used. The light was split into object and reference arm. The object arm employing the AOL with 20x objective was designed. The optical elements (i.e. objective) in the reference arm were used to avoid large dispersion mismatch. The cuvette of the length L = 50 mm filled with water was included in the reference arm to compensate for the water present in the AOC. The custom-built spectrometer was attached to the detection arm of the interferometer. The line-scan CCD camera (AViiVA SM2, 2048 pixels; e2v, UK) was operating at 10 kHz (exposure time 100 μs). The spectral resolution of the spectrometer provided depth range of OCT images of 1.3 mm. The axial and transverse resolutions of the system were 2.7 μm and 20 μm, respectively. The sensitivity of the system was 91 dB when the sample was illuminated with the beam of 0.8 mW power and the system operated at the speed of 10 kHz.

4. Results

The experiments performed in this study were divided into two main parts: (i) demonstration of high-speed acousto-optic tunable lens operation and its performance, (ii) application of tunable AOL in OCT imaging.

4.1 Acousto-optic system calibration

First of all, we calibrated the acousto-optic system, i.e. we found the relation between voltage applied to the electrodes of the transducer and the Raman-Nath parameter v:

v=kδnL.
Raman-Nath parameter describes the maximum phase retardation of the light wave when passing through the acoustic field.

The light beam was not modulated prior to entering the AOC (continuous illumination; AOM off). We placed the beam profiler exactly in the focal plane of the objective (offset) lens OL; Fig. 2(a). The light intensity was captured for increasing voltage values. Although voltage was an easy variable in the experiment, theoretical approach required more objective input value such as the Raman-Nath parameter. Additionally, numerical simulations of the experiment were performed using light, ultrasound and AOC parameters. The calibration itself was based on scaling the abscissas of the experimental points to the theoretical curve. As a result, the scaling factor between the voltage and the Raman-Nath parameter was obtained.

The plot in Fig. 3 shows that light intensity at the axis of symmetry of the system in the far field decreases with increasing ultrasound amplitude. The intensity of light is normalized to the intensity when the AOC is turned off (v = 0 or U = 0 V). We also observe that the experimental results agree very well with the theoretical calculation based on the model presented in Section 2. The obtained results do not depend on the power of the offset lens inserted to the system. It means that comparison of light intensity for particular voltage to that for U = 0 V (AOC off) enables proper system calibration that results in finding corresponding the Raman-Nath parameter value.

 figure: Fig. 3

Fig. 3 Calibration of the acousto-optic cell (F = 1.042 MHz, λ = 810 nm). Normalized light intensity vs. voltage applied to the transducer.

Download Full Size | PDF

4.2 Light beam profiling

In the next experiment, we aimed at three-dimensional reconstruction of the light beam. Continuous illumination scheme was used here. The beam profiler was attached to the translation stage, and the movement of the stage was synchronized with the beam profile acquisition. Three beam profiles were acquired for each axial position to minimize intensity fluctuations.

Figure 4 demonstrates cross-sectional images of the beam behind the tunable AOL when the tunable lens is off and on. However, better visualization of depth of focus (DOF) extension can be achieved if axial light intensity profile is extracted (right panel in Fig. 4). Two offset lenses were used here for comparison: fOL = 40 mm; Fig. 4(a); and 20x microscope objective; Fig. 4(b). Although the maximum axial intensity decreases when the voltage is applied to the AOL, the total energy is conserved.

 figure: Fig. 4

Fig. 4 Profiling of the light beam behind the acousto-optic tunable lens. Time-averaged horizontal sections of the light beam (left column – theory, central column – experiment). Time-averaged central depth profile of the light beam (right column). (a) Offset lens fOL = 40 mm; (b) offset lens 20x objective fOL = 10 mm.

Download Full Size | PDF

The ability to tune the focal position offered by the AOL depends on the amplitude of ultrasound (Raman-Nath parameter) and the focal distance of the offset lens. This effect was studied theoretically by numerical calculations of time-averaged axial intensity profiles using BPM described in Section 2. In order to calculate the depth of focus (DOF), we used the definition similar to that in Gaussian optics, i.e. DOF was defined as the distance between two axial points at which the light intensity has decreased to 50% of its peak intensity (i.e. full width at half maximum of the axial intensity profile). Figure 5 shows the DOF as a function of the Raman-Nath parameter and the focal distance of the offset lens.

 figure: Fig. 5

Fig. 5 Analysis of the depth of focus (DOF) extension introduced by the dynamic focusing. The calculations were performed using BPM. (a) Surface plot of the DOF behind the AOL with respect of Raman-Nath parameter v and focal length fOL of the offset lens. (b) Time-averaged axial profiles of the normalized light beam intensity for different values of Raman-Nath parameter and for fOL = 40 mm. DOF is defined as the FWHM of the axial intensity profile. (c) Dependence of the extracted DOF on Raman-Nath parameter for fOL = 40 mm. (d) Dependence of the extracted DOF on focal length of the offset lens for v = 10.

Download Full Size | PDF

The relations from numerical simulations presented in Figs. 5(c) and 5(d) can be also confirmed qualitatively when the AOL is considered as a compound lens consisting of two phase elements (lenses): a lens of variable focus (AOC) and an offset lens. Assuming that the AOC operates in the Raman-Nath regime (pure phase object), total phase retardation of light wave passing through the AOC can be calculated using Eq. (2). Next, taking into account approximation of the Bessel function for small arguments: J0(x)1x2/4, the corresponding power (reciprocal of the focal length) of the AOC can be derived in the following form [37]:

1fAOC(t)=K2v2kcos(Ωt)=Qn0v2Lcos(Ωt),
where QK2L/(kn0)is the Klein-Cook parameter for our configuration of acousto-optic interaction. Therefore, the effective focal length of the AOL is given by the formula for compound lens:
1fAOL(t)=1fAOC+1fOL=1fOL+K2v2kcos(Ωt),
where fOL is the focal length of the offset lens. The focal point oscillates in time around the position set by the offset lens.

Considering extreme positions of the focal point in time, we can estimate:

DOF4kK2fOL2v4k2(KfOLv)2K2kfOL2v
since the approximation 4k2>>(KfOLv)2 is valid in our conditions. Equation (10) confirms the calculations presented in Fig. 5(c) (for v ≥ 2) and Fig. 5(d). The fact that DOF is proportional to the square of the focal length of the offset lens corresponds also to the DOF formula derived in Gaussian optics for a single lens.

4.3 Time-resolved tunable lens operation

In order to determine the temporal performance of the designed AOL, special illumination scheme was introduced. Application of EOM enabled obtaining pulsed operation at the frequency of the ultrasound that was synchronized with the ultrasound phase. In a pulsed operation, we took advantage of the stroboscopic effect by illumination in a particular, well controlled ultrasound phase; Fig. 2(c) [43]. As a result, it enabled application of the detectors that do not have high temporal resolution. In particular, illumination synchronous with the ultrasound phase allowed for precise determination of the dynamics of the light field due to acousto-optic effect.

We acquired beam profiles when the time delay (phase) between both signals was changed for different axial positions of the detector. Central light intensity from each image was extracted and time-resolved light intensity profiles are shown in Fig. 6.

 figure: Fig. 6

Fig. 6 Dynamic focusing with acousto-optic tunable lens at the frequency F = 1.042 MHz. Maps of the evolution of the axial light intensity profile within the period T of ultrasound (top) and temporal changes of light intensity for particular planes indicated in the maps (bottom; lines – theory, dots – experiment). (a) Offset lens fOL = 40 mm; (b) offset lens 20x objective fOL = 10 mm.

Download Full Size | PDF

The focus crosses a given plane twice within a period of ultrasound because the ultrasonic wave generated inside the transducer is actually a standing wave. Moreover, the results demonstrate the dynamic focusing effect due to the acousto-optic interaction of light beam with the ultrasound of cylindrical symmetry.

We also investigated theoretically how the transverse resolution changes when the AOL performs temporal wavefront modulation. Comparison of beam profiles at the extreme positions of the focus (i.e. for t = 0 and t = ½T) indicated that the beam spot size increased by ~30% (for fOL = 40 mm) or ~10% (for 20x objective) with respect to the beam width when no ultrasound is present (v = 0). This change depended on the ultrasound amplitude. Moreover, the set-up with higher numerical aperture was less prone to those spot size variations.

4.4 OCT system performance

The performance of the OCT system equipped with the high-speed AOL in one of the arms is altered by rapidly varying phase of light wave travelling through the AOL. Therefore, the impact of the light interaction with ultrasonic field on the interference fringes acquired by the OCT system is an essential part of the study. This enables demonstrating the strengths and limitations of the proposed approach with respect to OCT method. Different factors influencing OCT signal were investigated in this part of study.

We measured the fringe amplitude and total intensity of light from the object arm captured by the spectrometer for increasing values of the Raman-Nath parameter (ultrasound amplitude) when continuous (EOM off) or stroboscopic (EOM on) illumination was used. The fringe amplitude was calculated from the peak amplitude after Fourier transformation, whereas the light intensity from the object arm was obtained by integrating the intensity from all spectral channels of the spectrometer in the case of blocked reference arm. The results are presented in Fig. 7(a). The interference signal drop is higher when the ultrasound is illuminated continuously (black dots in Fig. 7(a)), and it cannot be fully explained by the lower signal coming from the object arm. This is the additional impact of the fringe averaging since all time points are illuminated and contribute to the observed spectral fringes. On the other hand, decrease in fringe amplitude coincides with the signal drop when a particular phase of ultrasound is illuminated by light pulses (red dots in Fig. 7(a)). Comparison of illumination schemes demonstrates also the advantage of stroboscopic (pulsed) illumination scheme in terms of interference signal drop with the amplitude of ultrasound (open black vs. open red dots in Fig. 7(a)).

 figure: Fig. 7

Fig. 7 Performance of the OCT system with acousto-optic tunable lens. Impact of acousto-optic lens operation on interference signal in OCT. (a) Dependence of the OCT signal (closed dots) and signal from the sample arm (open dots) on the Raman-Nath parameter for continuous and pulsed illumination. (b) Signal-to-noise ratio roll-offs for different ultrasound amplitudes (reference arm was moved). (c) OCT signal and sample arm signal vs. object mirror position for continuous illumination. (d) OCT signal and sample arm signal vs. object mirror position for different phases of pulsed (stroboscopic illumination). (e) Impact of pulse duration on the fringe wash-out.

Download Full Size | PDF

Later on, we wanted to check if the AOL influences the signal drop with depth. Spectral fringes were acquired whereas the reference mirror was moved. The plots in Fig. 7(b) show that the signal-to-noise ratio (SNR) roll-off remains the same but the SNR itself decreases with higher Raman-Nath parameters (as shown in Fig. 7(a)).

In the next experiment, the object mirror was translated from the offset lens. The results in Figs. 7(c) and 7(d) have been obtained for continuous and pulsed illumination, respectively. Both the amplitude of interference fringes and light intensity collected from the sample arm reconstruct temporal behavior of the focal point in the case of operating AOL; Fig. 7(d). Similar to the observations in Fig. 7(a), time-averaged profile of the amplitude of interference fringes does not correspond to the profile of collected light intensity, which indicates that additional fringe averaging (wash-out) is observed; Fig. 7(c).

Finally, we checked how the fringe amplitude decreases when we increase the pulse width of light illuminating the AOC. It is demonstrated in Fig. 7(e) that short illumination pulses synchronous with the ultrasound frequency make it possible to avoid OCT signal drop due to the phase averaging. When the pulse becomes longer, the effect is closer to the case of continuous illumination scheme.

4.5 Optical coherence microscopy with extended focus

The above-mentioned steps enabled optimization of our OCT set-up for imaging. In the final stage of our experiments, we performed imaging of different samples with the OCT system equipped with the tunable AOL; Fig. 2(b). The phantoms were selected and designed to visualize tunable focusing with the AOL. The samples / phantoms were scanned laterally using a motorized stage. The scanning was repeated for different phase values between the ultrasound and illumination pulse (delay τ in Fig. 2) so that one obtained a set of cross-sectional images of the same object with different axial position of the focus. In the post-processing, we selected the centers (depth positions) of the focus for each image and performed Gaussian windowing of the image. The width of Gaussian window depended on the number of partial image acquisitions (number of phase steps). The final (combined) image with extended illumination depth was obtained by combining the images filtered with Gaussian windows.

Figure 8(a) shows cross-sectional images of the stack of microscopic cover glasses. The focal position tuning allowed for enhancing light reflection from a particular glass interface.

 figure: Fig. 8

Fig. 8 OCT imaging with dynamic focusing. (a) Imaging of the stack of microscopic cover glasses. (b) Imaging of the elastomer phantom with titanium dioxide. Averaged depth profiles demonstrating increased light collection efficiency from selected depths during different phases of focus tuning (Media 1).

Download Full Size | PDF

Furthermore, a scattering phantom (an elastomer with titanium dioxide) was used to demonstrate the impact of moving focus on the amplitude of collected light from specific depths. High-definition B-scans at different phase delays between the pulse and ultrasound wave were acquired. The images presented in Fig. 8(b) show that the focus of the beam can be controlled and adjusted for the sake of image improvement. Stroboscopic illumination allows for enhancement of photon collection efficiency from particular depth. This effect is even better visualized when one generates an averaged depth profile (A-scan), as indicated in the plot in Fig. 8(b). The amplitude of light in deeper depths is lower due to the multiple scattering that limits light penetration into the sample.

Finally, to demonstrate the applicability of our new set-up we imaged a lemon as an example of the biological object. The scanning protocol was similar to that used in previous experiments. The images in Fig. 9 clearly demonstrate that the image quality in selected depths can be enhanced using dynamic focusing. The signal from the out-of-focus depths is blurred and reduced due to the low collection efficiency. When we control the position of the focus, it is possible to improve the image quality at particular depth. In the panel demonstrating the combined image we were able to reconstruct sharp cross-sectional image of the lime internal structure with fairly good and homogenous backscattered intensity. High contrast reconstruction was obtained for the total imaging depth two times larger than that in the original image without extended DOF.

 figure: Fig. 9

Fig. 9 OCT imaging of the lime. Selected frames of the movie showing focal position tuning by stroboscopic illumination. The left cross-section shows the image without dynamic focusing. The right image is the composite cross-sectional image with extended DOF obtained by averaging different depth positions (Media 2).

Download Full Size | PDF

5. Discussion

The tunable acousto-optic lens presented in this study operates in MHz frequency range. Therefore, it allows for ultrahigh-speed tuning of the focus position. The ultrahigh tuning rate of the focus demonstrated here is faster than commercially available tunable lenses. This feature enables development of novel approaches to scanning scenarios. Although the tuning speed is fixed by the resonance mode of shell wall vibrations (thus, it depends on the thickness of the walls and properties of the piezoelectric material), the acousto-optic lens enables full control of the axial tuning range. To our best knowledge, this is the fastest method for dynamic focusing. Another advantage of using axisymmetric acousto-optic element is that it does not require high electrical power levels to obtain the desired effect. This is due to the fact that the ultrasound energy is accumulated at the axis of the vibrating piezoelectric shell. Moreover, due to the nature of the ultrasonic field, the AOL presented in this study does not introduce any non-axisymmetric aberrations like tunable electro-optic liquid lenses. We did not observe any chromatic effects when the AOL was operating although spectrally broadband light was used in this study.

The implemented successive diffraction model of the interaction of ultrasound with light demonstrates general numerical approach to solving the problem of light propagation through non-homogenous medium. We obtained very good agreement between numerical simulations of the AOL performance and experimental results, which indicates that the numerical calculation technique based on Fourier optics is correct.

Tunable focusing generates time-averaged light intensity distribution that is characterized by the extended DOF. The degree of DOF extension depends linearly on square of the focal length of the offset lens and on the amplitude of ultrasound (Raman-Nath parameter). This is the feature enabling improvement of capabilities of imaging devices.

We benefited from the pulsed illumination strategy implemented in our study. Although similar approach was applied in reducing motion artifacts [43, 44] we applied illumination synchronous with the focus tuning frequency in order to demonstrate cross-sectional imaging where focus position is well-defined. Analogical effect can be observed when the detector (camera) is gated e.g. with the shutter. However, this is not optimum since most of light is not detected. Stroboscopic effect may be also applied when ultrafast periodic changes in the sample properties are expected to be observed. Therefore, it is especially suited for studies of rapid phenomena.

We have presented careful assessment of utilization of high-speed tunable lens in the OCT technology. Although slower rates of focus tuning (ultrasound frequency) would be better for OCT to demonstrate its applicability, advances in OCT technology might fully take advantage of the presented approach. Dynamic focusing (or focus position adjustment) in OCT was demonstrated with different technologies [13, 18, 20, 21, 25]. Most of them enable static tuning and the tuning speed is limited but this report explores impact of temporal focus behavior on OCT signal when tuning at much faster time scales is provided. Although acousto-optic tunable lenses have been already employed in other imaging techniques [37, 45–47], this is the first demonstration of interference-based imaging with fast acousto-optic lens.

The control and tuning of focus position along with stroboscopic illumination enable better collection of light from deeper depths of the scattering sample. Furthermore, dynamic focusing applied to imaging allows also for increase in effective depth of focus when the system with high numerical aperture is used. However, the proposed method does not overcome the light penetration limits. Another challenging issue is the fact that the operation of the AOL impacts the interferometric signal in coherence-based imaging, which requires careful choice of the illumination method.

7. Summary and conclusions

To summarize, ultrahigh-speed tunable lens based on acousto-optic element was designed and implemented into OCT system. The results demonstrated the ability of the device to obtain ultrahigh-speed (MHz frequency range) dynamic focusing in OCT with stroboscopic illumination. We believe that imaging with active optical elements can improve photon collection efficiency, depth of focus and enhance image quality in OCM as well as in long depth range OCT.

Acknowledgments

IG and MW acknowledge support from the National Science Center within the MAESTRO Programme (#2011/02/A/ST2/00302). KSz and MW acknowledge the TEAM project (TEAM/2011-8/8) financed by the European Union within the frame of Operational Programme Innovative Economy coordinated by the Foundation for Polish Science. Additionally, IG is supported by the scholarship of the Polish Ministry of Science and Higher Education for 2013-2016. Finally, the authors thank Szymon Tamborski, MSc for technical help.

References and links

1. D. Huang, E. A. Swanson, C. P. Lin, J. S. Schuman, W. G. Stinson, W. Chang, M. R. Hee, T. Flotte, K. Gregory, C. A. Puliafito, and J. G. Fujimoto, “Optical coherence tomography,” Science 254(5035), 1178–1181 (1991). [CrossRef]   [PubMed]  

2. J. A. Izatt, M. R. Hee, G. M. Owen, E. A. Swanson, and J. G. Fujimoto, “Optical coherence microscopy in scattering media,” Opt. Lett. 19(8), 590–592 (1994). [CrossRef]   [PubMed]  

3. C. Zhou, D. W. Cohen, Y. Wang, H.-C. Lee, A. E. Mondelblatt, T.-H. Tsai, A. D. Aguirre, J. G. Fujimoto, and J. L. Connolly, “Integrated Optical Coherence Tomography and Microscopy for Ex Vivo Multiscale Evaluation of Human Breast Tissues,” Cancer Res. 70(24), 10071–10079 (2010). [CrossRef]   [PubMed]  

4. V. J. Srinivasan, H. Radhakrishnan, J. Y. Jiang, S. Barry, and A. E. Cable, “Optical coherence microscopy for deep tissue imaging of the cerebral cortex with intrinsic contrast,” Opt. Express 20(3), 2220–2239 (2012). [PubMed]  

5. I. Grulkowski, J. J. Liu, B. Potsaid, V. Jayaraman, J. Jiang, J. G. Fujimoto, and A. E. Cable, “High-precision, high-accuracy ultralong-range swept-source optical coherence tomography using vertical cavity surface emitting laser light source,” Opt. Lett. 38(5), 673–675 (2013). [CrossRef]   [PubMed]  

6. H.-C. Lee, J. J. Liu, Y. Sheikine, A. D. Aguirre, J. L. Connolly, and J. G. Fujimoto, “Ultrahigh speed spectral-domain optical coherence microscopy,” Biomed. Opt. Express 4(8), 1236–1254 (2013). [CrossRef]   [PubMed]  

7. J. G. Fujimoto and W. Drexler, eds., Optical Coherence Tomography. Technology and Applications (Springer, 2008).

8. M. Wojtkowski, “High-speed optical coherence tomography: basics and applications,” Appl. Opt. 49(16), D30–D61 (2010). [CrossRef]   [PubMed]  

9. J. M. Schmitt, S. L. Lee, and K. M. Yung, “An optical coherence microscope with enhanced resolving power in thick tissue,” Opt. Commun. 142(4-6), 203–207 (1997). [CrossRef]  

10. F. Lexer, C. K. Hitzenberger, W. Drexler, S. Molebny, H. Sattmann, M. Sticker, and A. F. Fercher, “Dynamic coherent focus OCT with depth-independent transversal resolution,” J. Mod. Opt. 46(3), 541–553 (1999). [CrossRef]  

11. J. Holmes, “Theory and applications of multi-beam OCT,” Proc. SPIE 7139, 713908 (2008). [CrossRef]  

12. J. Holmes and S. Hattersley, “Image blending and speckle noise reduction in multi-beam OCT,” Proc. SPIE 7168, 71681N (2009). [CrossRef]  

13. B. Qi, A. Phillip Himmer, L. Maggie Gordon, X. D. Victor Yang, L. David Dickensheets, and I. Alex Vitkin, “Dynamic focus control in high-speed optical coherence tomography based on a microelectromechanical mirror,” Opt. Commun. 232(1-6), 123–128 (2004). [CrossRef]  

14. Z. Ding, H. Ren, Y. Zhao, J. S. Nelson, and Z. Chen, “High-resolution optical coherence tomography over a large depth range with an axicon lens,” Opt. Lett. 27(4), 243–245 (2002). [CrossRef]   [PubMed]  

15. R. A. Leitgeb, M. Villiger, A. H. Bachmann, L. Steinmann, and T. Lasser, “Extended focus depth for Fourier domain optical coherence microscopy,” Opt. Lett. 31(16), 2450–2452 (2006). [CrossRef]   [PubMed]  

16. C. Blatter, B. Grajciar, C. M. Eigenwillig, W. Wieser, B. R. Biedermann, R. Huber, and R. A. Leitgeb, “Extended focus high-speed swept source OCT with self-reconstructive illumination,” Opt. Express 19(13), 12141–12155 (2011). [CrossRef]   [PubMed]  

17. D. Lorenser, C. Christian Singe, A. Curatolo, and D. D. Sampson, “Energy-efficient low-Fresnel-number Bessel beams and their application in optical coherence tomography,” Opt. Lett. 39(3), 548–551 (2014). [CrossRef]   [PubMed]  

18. M. Pircher, E. Götzinger, and C. K. Hitzenberger, “Dynamic focus in optical coherence tomography for retinal imaging,” J. Biomed. Opt. 11(5), 054013 (2006). [CrossRef]   [PubMed]  

19. Y.-K. Fuh and M.-X. Lin, “Adaptive optics correction of a tunable fluidic lens for ophthalmic applications,” Opt. Commun. 308, 100–104 (2013). [CrossRef]  

20. A. Divetia, T.-H. Hsieh, J. Zhang, Z. Chen, M. Bachman, and G.-P. Li, “Dynamically focused optical coherence tomography for endoscopic applications,” Appl. Phys. Lett. 86(10), 103902 (2005). [CrossRef]  

21. P. Meemon, K.-S. Lee, S. Murali, and J. Rolland, “Optical design of a dynamic focus catheter for high-resolution endoscopic optical coherence tomography,” Appl. Opt. 47(13), 2452–2457 (2008). [CrossRef]   [PubMed]  

22. C. A. Lopez and A. H. Hirsa, “Fast focusing using a pinned-contact oscillating liquid lens,” Nat. Photonics 2(10), 610–613 (2008). [CrossRef]  

23. J. Kang, H. Yu, and H. Chen, “Liquid tunable lens integrated with a rotational symmetric surface for long depth of focus,” Appl. Opt. 49(28), 5493–5500 (2010). [CrossRef]   [PubMed]  

24. M. Duocastella, B. Sun, and C. B. Arnold, “Simultaneous imaging of multiple focal planes for three-dimensional microscopy using ultra-high-speed adaptive optics,” J. Biomed. Opt. 17(5), 050505 (2012). [CrossRef]   [PubMed]  

25. H.-S. Chen and Y.-H. Lin, “An endoscopic system adopting a liquid crystal lens with an electrically tunable depth-of-field,” Opt. Express 21(15), 18079–18088 (2013). [CrossRef]   [PubMed]  

26. F. O. Fahrbach, F. F. Voigt, B. Schmid, F. Helmchen, and J. Huisken, “Rapid 3D light-sheet microscopy with a tunable lens,” Opt. Express 21(18), 21010–21026 (2013). [CrossRef]   [PubMed]  

27. K. C. Heo, S. H. Yu, J. H. Kwon, and J. S. Gwag, “Thermally tunable-focus lenticular lens using liquid crystal,” Appl. Opt. 52(35), 8460–8464 (2013). [CrossRef]   [PubMed]  

28. J. M. Jabbour, B. H. Malik, C. Olsovsky, R. Cuenca, S. Cheng, J. A. Jo, Y.-S. L. Cheng, J. M. Wright, and K. C. Maitland, “Optical axial scanning in confocal microscopy using an electrically tunable lens,” Biomed. Opt. Express 5(2), 645–652 (2014). [CrossRef]   [PubMed]  

29. D. Hillmann, C. Lührs, T. Bonin, P. Koch, and G. Hüttmann, “Holoscopy--holographic optical coherence tomography,” Opt. Lett. 36(13), 2390–2392 (2011). [CrossRef]   [PubMed]  

30. J. Mo, M. de Groot, and J. F. de Boer, “Focus-extension by depth-encoded synthetic aperture in Optical Coherence Tomography,” Opt. Express 21(8), 10048–10061 (2013). [PubMed]  

31. A. Kumar, W. Drexler, and R. A. Leitgeb, “Subaperture correlation based digital adaptive optics for full field optical coherence tomography,” Opt. Express 21(9), 10850–10866 (2013). [CrossRef]   [PubMed]  

32. Y.-Z. Liu, N. D. Shemonski, S. G. Adie, A. Ahmad, A. J. Bower, P. S. Carney, and S. A. Boppart, “Computed optical interferometric tomography for high-speed volumetric cellular imaging,” Biomed. Opt. Express 5(9), 2988–3000 (2014). [CrossRef]   [PubMed]  

33. I. Grulkowski and P. Kwiek, “Experimental study of light diffraction by standing ultrasonic wave with cylindrical symmetry,” Opt. Commun. 267(1), 14–19 (2006). [CrossRef]  

34. E. McLeod, A. B. Hopkins, and C. B. Arnold, “Multiscale Bessel beams generated by a tunable acoustic gradient index of refraction lens,” Opt. Lett. 31(21), 3155–3157 (2006). [CrossRef]   [PubMed]  

35. E. McLeod and C. B. Arnold, “Mechanics and refractive power optimization of tunable acoustic gradient lenses,” J. Appl. Phys. 102(3), 033104 (2007). [CrossRef]  

36. I. Grulkowski, D. Jankowski, and P. Kwiek, “Acousto-optic interaction of a Gaussian laser beam with an ultrasonic wave of cylindrical symmetry,” Appl. Opt. 46(23), 5870–5876 (2007). [CrossRef]   [PubMed]  

37. A. Mermillod-Blondin, E. McLeod, and C. B. Arnold, “High-speed varifocal imaging with a tunable acoustic gradient index of refraction lens,” Opt. Lett. 33(18), 2146–2148 (2008). [CrossRef]   [PubMed]  

38. M. Bass, E. W. Van Stryland, D. R. Williams, and W. L. Wolfe, Handbook of optics (McGraw-Hill, 2001), Vol. 2.

39. A. Korpel, Acousto-optics (CRC, 1996).

40. I. Grulkowski and P. Kwiek, “Successive diffraction model based on Fourier optics as a tool for the studies of light interaction with arbitrary ultrasonic field,” Eur. Phys. J. Spec. Top. 154(1), 77–83 (2008). [CrossRef]  

41. C. Koch and R. Reibold, (personal communication, 2007).

42. R. Reibold and P. Kwiek, “Optical Near-field Investigation into the Raman-Nath and KML Regimes of Diffraction by Ultrasonic Waves,” Acta Acust. Acustica 70, 223–229 (1990).

43. S. H. Yun, G. Tearney, J. de Boer, and B. Bouma, “Pulsed-source and swept-source spectral-domain optical coherence tomography with reduced motion artifacts,” Opt. Express 12(23), 5614–5624 (2004). [CrossRef]   [PubMed]  

44. G. Moneron, A. C. Boccara, and A. Dubois, “Stroboscopic ultrahigh-resolution full-field optical coherence tomography,” Opt. Lett. 30(11), 1351–1353 (2005). [CrossRef]   [PubMed]  

45. J. Yan, A. Mermillod-Blondin, E. McLeod, and C. B. Arnold, “Rapidly tunable acoustic gradient index lenses for pulsed imaging and laser processing,” in Conference on Lasers and Electro-Optics and Quantum Electronics and Laser Science Conference, (Institute of Electrical and Electronics Engineers, 2009), 106–107. [CrossRef]  

46. N. Olivier, A. Mermillod-Blondin, C. B. Arnold, and E. Beaurepaire, “Two-photon microscopy with simultaneous standard and extended depth of field using a tunable acoustic gradient-index lens,” Opt. Lett. 34(11), 1684–1686 (2009). [CrossRef]   [PubMed]  

47. M. Duocastella and C. B. Arnold, “Enhanced depth of field laser processing using an ultra-high-speed axial scanner,” Appl. Phys. Lett. 102(6), 061113 (2013). [CrossRef]  

Supplementary Material (2)

Media 1: MP4 (96 KB)     
Media 2: MP4 (2468 KB)     

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (9)

Fig. 1
Fig. 1 Geometry of acousto-optic interaction. (a) Scheme of the acousto-optic cell. The cell is composed of cylindrical transducer of the length L. Light beam passes along the axis of the transducer. (b) The numerical approach assumes dividing the acoustic medium into M layers, each being regarded as pure phase grating. (c) Vibration of the transducer generates standing ultrasonic wave in which concentration of ultrasound is observed at the axis (r = 0). Gaussian distribution of electric field amplitude and the distribution of acoustic pressure. (d) Light beam propagation in a single layer.
Fig. 2
Fig. 2 Experimental set-up. (a) Configuration for experimental verification of the operation of acousto-optic tunable lens. (b) Spectral / Fourier-domain OCT instrument with tunable lens. EOM – electro-optic intensity modulator, FC – fiber coupler, AOC – acousto-optic cylindrical cell, OL – offset lens, BP – beam profiler, FG – function generator, DDG – digital delay line, PD – high-speed photodetector, OSC – oscilloscope, DC – dispersion compensation, LSC – line-scan camera, COMP – computer. (c) Explanation of continuous and stroboscopic illumination schemes. Application of EOM and digital delay line allows for control of phase difference between light pulses and ultrasound as well as to adjust the width of light pulses.
Fig. 3
Fig. 3 Calibration of the acousto-optic cell (F = 1.042 MHz, λ = 810 nm). Normalized light intensity vs. voltage applied to the transducer.
Fig. 4
Fig. 4 Profiling of the light beam behind the acousto-optic tunable lens. Time-averaged horizontal sections of the light beam (left column – theory, central column – experiment). Time-averaged central depth profile of the light beam (right column). (a) Offset lens fOL = 40 mm; (b) offset lens 20x objective fOL = 10 mm.
Fig. 5
Fig. 5 Analysis of the depth of focus (DOF) extension introduced by the dynamic focusing. The calculations were performed using BPM. (a) Surface plot of the DOF behind the AOL with respect of Raman-Nath parameter v and focal length fOL of the offset lens. (b) Time-averaged axial profiles of the normalized light beam intensity for different values of Raman-Nath parameter and for fOL = 40 mm. DOF is defined as the FWHM of the axial intensity profile. (c) Dependence of the extracted DOF on Raman-Nath parameter for fOL = 40 mm. (d) Dependence of the extracted DOF on focal length of the offset lens for v = 10.
Fig. 6
Fig. 6 Dynamic focusing with acousto-optic tunable lens at the frequency F = 1.042 MHz. Maps of the evolution of the axial light intensity profile within the period T of ultrasound (top) and temporal changes of light intensity for particular planes indicated in the maps (bottom; lines – theory, dots – experiment). (a) Offset lens fOL = 40 mm; (b) offset lens 20x objective fOL = 10 mm.
Fig. 7
Fig. 7 Performance of the OCT system with acousto-optic tunable lens. Impact of acousto-optic lens operation on interference signal in OCT. (a) Dependence of the OCT signal (closed dots) and signal from the sample arm (open dots) on the Raman-Nath parameter for continuous and pulsed illumination. (b) Signal-to-noise ratio roll-offs for different ultrasound amplitudes (reference arm was moved). (c) OCT signal and sample arm signal vs. object mirror position for continuous illumination. (d) OCT signal and sample arm signal vs. object mirror position for different phases of pulsed (stroboscopic illumination). (e) Impact of pulse duration on the fringe wash-out.
Fig. 8
Fig. 8 OCT imaging with dynamic focusing. (a) Imaging of the stack of microscopic cover glasses. (b) Imaging of the elastomer phantom with titanium dioxide. Averaged depth profiles demonstrating increased light collection efficiency from selected depths during different phases of focus tuning (Media 1).
Fig. 9
Fig. 9 OCT imaging of the lime. Selected frames of the movie showing focal position tuning by stroboscopic illumination. The left cross-section shows the image without dynamic focusing. The right image is the composite cross-sectional image with extended DOF obtained by averaging different depth positions (Media 2).

Equations (10)

Equations on this page are rendered with MathJax. Learn more.

p( r,t )=δp J 0 ( Kr )cos( Ωt ),
n( r,t )= n 0 +Δn( r,t )= n 0 +δn J 0 ( Kr )cos( Ωt ),
E( r )= E 0 e ( r/w ) 2 ,
Q v 0.2,
E ˜ ( f x , f y ,z+Δz )= E ˜ in ( f x , f y ,z )exp[ iΔz ( k n 0 ) 2 ( 2π ) 2 ( f x 2 + f y 2 ) ],
E out ( x,y,z+Δz )=E( x,y,z+Δz )exp[ ikΔn( r,t ) ].
v=kδnL.
1 f AOC ( t ) = K 2 v 2k cos( Ωt )= Q n 0 v 2L cos( Ωt ),
1 f AOL ( t ) = 1 f AOC + 1 f OL = 1 f OL + K 2 v 2k cos( Ωt ),
DOF 4k K 2 f OL 2 v 4 k 2 ( K f OL v ) 2 K 2 k f OL 2 v
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.