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

Optimization of galvanometer scanning for optical coherence tomography

Open Access Open Access

Abstract

We study experimentally the effective duty cycle of galvanometer-based scanners (GSs) with regard to three main parameters of the scanning process: theoretical/imposed duty cycle (of the input signal), scan frequency, and scan amplitude. Sawtooth and triangular input signals for the device are considered. The effects of the mechanical inertia of the oscillatory element of the GS are analyzed and their consequences are discussed in the context of optical coherence tomography (OCT) imaging. When the theoretical duty cycle and the scan amplitude are increased to the limit, the saturation of the device is demonstrated for a useful range of scan frequencies by direct measurement of the position of the galvomirror. Investigations of OCT imaging of large samples also validate this saturation, as examplified by the gaps/blurred portions obtained between neighboring images when using both triangular and sawtooth scanning at high scan frequencies. For this latter aspect, the necessary overlap between neighboring B-scans, and therefore between the corresponding volumetric reconstructions of the sample, are evaluated and implemented with regard to the same parameters of the scanning process. OCT images that are free of these artifacts are thus obtained.

© 2015 Optical Society of America

1. Introduction

Galvanometer-based scanners (GSs) are the most utilized devices for lateral scanning in various biomedical imaging techniques, including optical coherence tomography (OCT) [1,2], confocal microscopy (CM), and combinations of different imaging modalities [3]. Two-dimensional (2D) scanning of the sample, referred to as raster scanning, is usually achieved in such systems with dual-axis double GSs [4,5] or, more rarely with a combination of a polygon mirror (PM) scanner for the fast axis, and a one-axis GS to position precisely the line scan on the slow axis [6,7]. Two-dimensional GS systems have also been employed in association with Risley prisms [8,9] for volumetric (3D) scanning [10,11] or for time-included (4D) scanning in the context of real time in vivo imaging [12].

The use of unidimensional (1D) GSs has been also reported in various applications including generating the modulation necessary in en face time domain (TD) OCT with both on-axis and off-axis GSs [13], scanning transmissive delay line for Fourier domain (FD) OCT [14], phase or frequency modulation [15], and handheld scanning probes for spectral domain (SD) OCT [16,17].

The extensive use of GSs in high-end applications like biomedical imaging is related to their numerous advantages [18]: precision positioning, a large range of scan frequencies and velocities, compact constructive solutions, and accurate positioning of their mirror on the oscillatory axis. The latter two aspects are considered advantageous compared to PM scanners [4,5,1921]. However, GSs provide lower scan frequencies and velocities compared to PMs and resonant scanners; the latter are usually built as microelectromechanical systems (MEMS). Nevertheless, the adjustable frequency of the GSs with regard to the fixed frequency of resonant devices is an advantage. Thus, different input signals can be employed for GSs whereas MEMS are usually driven with sinusoidal functions. As we have shown in previous work [22], triangular scanning provides artifact-free OCT images (with appropriate scan parameters), whereas sinusoidal scanning always introduces significant distortion of images. Such aspects have to be compensated in postprocessing in order to obtain distortion-free images in the case of sinusoidal scanning with MEMS [23].

GS technology is mature in terms of constructive and functional aspects [21]: motors, mobile elements, bearings, position sensors, heatsinkers, control structures, and testing. Some issues still remain unsolved, for instance the optimal functions necessary to obtain a maximum duty cycle. We have demonstrated theoretically in this respect that, for a maximum time efficiency of the scanning process (i.e., a maximum duty cycle of the GS), one has to use—for symmetrical/triangular scanning—linear plus parabolic scanning functions [24], while in the literature the linear plus sinusoidal functions have been previously considered best from this point of view [21].

However, the biomedical imaging community does not usually rely on custom-made input signals for the GSs but on common input signals like sawtooth, purely triangular (i.e., without nonlinear portions), and sinusoidal. In previous work, we have investigated experimentally these different scanning signals and their effects on OCT imaging [22]. Rules of thumb have been proposed to allow for an optimal use of these scanning devices in OCT setups, and some of these aspects will be briefly discussed at the beginning of this study. We extended in the present work these investigations from 1D GSs to 2D dual-axis GSs in OCT applications. All OCT systems that employ GSs for lateral scanning are targeted by these studies. Thus, while in the previous study an FD-OCT system was used to explore the characteristics of 1D galvoscanning, the effect of the mechanical inertia of the GS oscillatory element on imaging has also been taken into account (based on the work in [22]) in Doppler OCT to obtain optimized images of the retina with backstitched B-scans [25].

The remainder of this paper is structured as follows: Section 2 discusses the response of the device to two different input signals (triangular and sawtooth) for the 1D GSs. The effect of the different input parameters on the scanning functions is analyzed. Particularly, the relationship between the theoretical duty cycle and the effective duty cycle determined experimentally is defined and studied. In Section 3, a mathematical model is developed for both the scanning functions and the effective duty cycle. Section 4 presents the custom Gabor domain optical coherence microscopy (GD-OCM) system that was utilized for experimental validation, while Section 5 points out the effect of triangular and sawtooth functions on the images. The necessary overlap of the individual scans when imaging larger samples is evaluated and applied to both scanning regimes in order to obtain artifact-free images. The effective duty cycles obtained with galvoscanning in the context of OCT imaging are finally compared.

2. Effective Duty Cycle of a 1D GS

A schematic illustration of the operating principle of a GS is shown in Fig. 1. The device is electrically driven with a nonresonant oscillatory element (magnet or coil [21]). It usually operates in a closed loop with different control structures designed to minimize the response time of the device and maximize its precision [26]. Thus, error corrections, including in relation to temperature variations, are addressed through appropriate feedback loops.

 figure: Fig. 1.

Fig. 1. Schematic illustrating the operating principle of a galvanometer-based scanner (GS), with constructive parameters: J, moment of inertia of the mobile element (with galvomirror); c, damping coefficient; k, elastic coefficient of the torsion springs that support the mobile element.

Download Full Size | PDF

Regardless of the complexity of the control module of the GS, the mechanical inertia does play a role in the functioning of the device. Thus, the galvomirror (together with the mobile element of the motor) cannot stop and turn instantaneously. As shown in [22], this is the reason why the scanning function of the GS (i.e., its output signal represented by the current angular position of the galvomirror) will actually differ from the input signal that usually has an ideal triangular or sawtooth profile [Fig. 2(a)]. This effect becomes significant when increasing the scan frequency (fs) above 100 Hz, especially at scan amplitudes (θm) above 1 V, yielding increased nonlinearity in the scanning function. Thus, as fs and θm are increased, the device will be less and less capable of following the input signal provided, as we further describe in this section for sawtooth input signals. A detailed study on the phenomena involved for different scanning regimes, as a consequence of the mechanical inertia, has already been reported in [22]. We then investigated three most common input signals of GSs (i.e., sawtooth, triangular, and sinusoidal) and demonstrated that sinusoidal scanning, although convenient from a mechanical point of view as it allows for a smoother functioning of the device, produces strong distortions in OCT images even at lower frequencies.

 figure: Fig. 2.

Fig. 2. Study of sawtooth scanning functions/output signals (i.e., angular positions of the galvomirror) of a GS for input signals with different theoretical duty cycles ηt equal to (a1) 50% (triangular function), (a2) 75%, and (a3) 87.5%; (b1)–(b3) output signals (positions of the galvomirror) for a scan frequency fs=50Hz; (c1)–(c3) output signals for fs=500Hz; (d1)–(d3) effective duty cycle of the GS with regard to fs for three different driving duty cycles and different scan amplitudes θm of 0.2, 0.4, 0.8, 1.6, and 3.2 V (where 1 V corresponds to about 7.6 deg of optical scan angle).

Download Full Size | PDF

Triangular scanning is the best choice, while sawtooth scanning has the flyback portion that needs to be discarded when imaging larger samples [25]. However, as previously pointed out, even triangular and sawtooth scanning regimes are affected by an increasing nonlinearity of the output signal at high frequencies and amplitudes. This nonlinearity produces distortions at the margins of the scanned field in OCT imaging. Note that these distortions can be compensated in postprocessing, as usually done with sinusoidal scans, for which equally spaced pixels of the final OCT images are thus provided [23]. However, postprocessing of the images is a process one may want to avoid, especially when large amounts of data are involved and when real-time imaging is required.

Another parameter that has to be taken into account is the duty cycle, defined as the temporal scan efficiency, i.e., the ratio between the “active” time ta (the time interval used to scan the sample with constant velocity) and the time period T. For the scanning regimes, the theoretical or ideal duty cycle ηt that is imposed and characterizes the input signals of the GS, is defined as [4]

ηt=ta/T.

The output signals of the GS (referred to as the scanning functions) will be characterized by an effective duty cycle (Fig. 3) whose value is lower than theoretically predicted because of nonlinearities. Indeed, the active time ta that characterizes the linear portion of the function is smaller than ta (Fig. 4), and the effective duty cycle η is given as [21]

η=ta/T.

 figure: Fig. 3.

Fig. 3. Study of the effective duty cycle of a GS with regard to the theoretical/programmed duty cycle, for different scan amplitudes θm of 0.2, 0.4, 0.8,1.6, and 3.2 V (1 V corresponds to 7.6 deg optically in terms of angular scan amplitude)—each study being made for a certain scan frequency fs: (a) 10 Hz; (b) 50 Hz; (c) 100 Hz; (d) 200 Hz; (e) 300 Hz; (f) 500 Hz.

Download Full Size | PDF

 figure: Fig. 4.

Fig. 4. (a) Triangular and (b1), (b2) sawtooth scanning regimes of the GS: ideal/input signals versus scanning functions/output signals of the GS. For the latter, the current positions of the galvomirror are those determined experimentally and shown in Fig. 2.

Download Full Size | PDF

In this study, three theoretical duty cycles ηt are considered for the sawtooth function: 50%, 75%, and 87.5%. The latter value is approximately the maximum duty cycle a signal generator may provide. The 50% value of ηt corresponds to a triangular input signal in which only a unidirectional scan is used. The case of both forward and backward scanning within one cycle with a triangular input signal corresponds to a triangular scan with ηt equals 100%. In this case, the effective duty cycle η also varies from 100% to 70% with regard to the scan frequency and amplitude (Fig. 10 in [22]). In OCT, two B-scans per scan cycle may be utilized for image averaging or imaging at high speed.

For the systematic study of the scanning functions in Figs. 2 and 3, the 50% value is the necessary lower limit of ηt for sawtooth input signals.

The scan frequency fs was varied (Figs. 2 and 3) from 10 Hz, which represents point-by-point scanning, up to 500 Hz, while intermediate values of interest for GSs and for OCT applications are considered. However, while most OCT systems operate at 2–200 Hz frame rate (assuming a number of A-scans in each frame of 500 and a typical A-scan rate of 1–100 kHz), scan frequencies of around 300 Hz have also been reported. The scan amplitude θm ranges from 0.2 to 3.2 V, which correspond to about 0.76–12.16 deg (mechanical angle of the oscillation of the galvomirror [5]) or 1.52–24.32 deg (optical angle of the reflected laser beam). This range of voltage is common for OCT investigations and is considered for the GSs, but the entire range was utilized only for lower scan frequencies as it can be remarked from Fig. 2(d), by example. This is because, as previously documented (Fig. 7 in [22]), there is a characteristic scan amplitude limit of a GS for each scan frequency and for each type of input signal; over this scan amplitude limit the GS system loses stability. Thus, different figures are required for the various GS voltages that are set by the signal generator of the GS.

The input signals are presented in Fig. 2(a), while output signals are shown in Figs. 2(b) and 2(c); the latter were obtained by directly reading the GS output voltage.

The output signals, which are proportional to the scanning angles, represent accurately the scanning locations on a sample surface when the equivalent output scan angles are multiplied by the focal length of the optics for a F-Theta scan lens used in its designed specification. If the lens is not an F-Theta scan lens or is not used according to its designed specification or suffers from optical distortion, the optical mapping must be established to relate the equivalent output scan angles to the scanning locations on the sample.

Several conclusions can be drawn from the multiparametric analysis performed in Figs. 2 and 3:

  • (i) At low scan frequencies (roughly up to 50 Hz), the position of the galvomirror can still follow the input signal provided to the device, except at the highest ηt where the nonlinearity begins to manifest after the flyback portion, as shown in Fig. 2(b3).
  • (ii) When increasing fs—for instance by a factor of 10 from Fig. 2(b) to Fig. 2(c)—a strong nonlinearity appears for all the signals considered, yielding a decrease in the effective duty cycle.
  • (iii) From Fig. 2(c) a comparison can be made between the nonlinearity of the stop-and-turn portion of the maximum amplitude (after the active scan) and that of the minimum amplitude (after the flyback portion). These two portions are equal in time (i.e., the curves are symmetrical) for the triangular scan, whereas they are longer after the flyback for the sawtooth function. This observation is consistent with the fact that the flyback is performed at a higher speed, therefore the device needs more time to stop-and-turn. This phenomenon is even more significant for a higher ηt.
  • (iv) As a consequence of the latter aspect, the effective active time ta and the corresponding effective duty cycle η decrease accordingly. An evaluation of the effective duty cycle η with regard to the scan frequency fs and scan amplitudes θm was performed for each value of ηt and reported in Fig. 2(d). Results show that the effective duty cycle decreases with the increasing of frequency and amplitude.
  • (v) A “migration” of the peaks of the output signals of the GS occurs for high values of fs—Figs. 2(c2) and 2(c3). Thus, ηt decreases even more because of this effect of switching from sawtooth to triangular functions at high fs and θm, an effect that was already demonstrated in [22].
  • (vi) While the η(fs) graphs in Fig. 2(d) start at the same value of η equal to ηt (for small fs and θm), a decrease in η is produced, as pointed out previously, with the increase in fs. This decrease in η is more significant for larger samples to be imaged (i.e., for higher θm). As expected, when θm gets too large, there is a limit value of fs that can still be employed. The systematic study of this aspect provided the limit possible scan amplitudes for each fs [22]; such limitations are illustrated in Figs. 3(e) and 3(f) by the lower reachable values of θm at the highest fs considered.
  • (vii) The spread in the values of η [corresponding graphs in Fig. 2(d)] increases with ηt. One may conclude that the GS is following the input signals with increasing difficulty for a higher ηt, even if the same fs and θm parameters are set.
  • (viii) The previous remark has suggested a study of η with regard to ηt, while maintaining fs and θm constant. This study, shown in Fig. 3, reveals a most peculiar aspect that should be taken into account by GSs users: while pushing the scanning device to the limit in terms of ηt, the effect may be opposite to the expectation when high fs and θm are employed. Thus, for higher fs (i.e., from around 200 Hz), a bending of the η(ηt) graphs is produced at the higher values of θm. While only three measurements were acquired per graph, which is not sufficient to precisely locate the peak of each curve, Figs. 3(d)3(f) show that there is a turning point to the effective duty cycle as a function of the theoretical duty cycle. Specifically, there will be a different position of the peak of the η(ηt) graphs with regard to the specific GS type, but there will always be a turning/bending point, for which a maximum in η is reached, with the remark that this happens at sufficiently high values of the scan amplitude (i.e., towards the amplitude limit where, for a specific scan frequency, the GS oscillation loses stability).

    This phenomenon is more significant at intermediate values of fs (i.e., 200 Hz). At even higher values (e.g., towards 500 Hz), all the values of η drop and, although the phenomenon shown is still present, it becomes less significant. Thus, for fs>300Hz, one may say that, while increasing ηt, η remains roughly constant; a saturation of the scanning regime is reached. However, the same practical conclusion remains valid for GS users: one must not push the device to the limit, as an increase in ηt over 75% does not provide the results expected, while it only has a negative impact on the device from a mechanical point of view.

  • (ix) As an overall tendency extracted from the study in Fig. 3, one may say that the η(ηt) function starts with an identity, regardless of θm for point-by-point scan [i.e., for low fs, Fig. 3(a)]. The graphs of these functions spread progresivelly for higher fs, with a maximum gradient of the spread around 200 Hz, then they get closer again at higher fs, because of the drop of all the values of η, for the now smaller interval of possible scan amplitudes.

3. Mathematical Model of the Effective Duty Cycle

From the experimental study performed (Fig. 2), the output signals/scanning functions (i.e., the current positions of the galvomirror) for both the triangular and for the sawtooth regime have been derived. Figure 4 presents the mathematical model of the phenomenon described above, where output functions are shown in contrast with the input signals.

  • (a) For the triangular input signals [Fig. 4(a)], the nonlinearity that occurs for the stop-and-turn portions produces a decrease in the scan amplitude from the theoretical value θm to the effective value θa that characterizes the linear portion of θ(t). The angular velocity increases from its theoretical value
    ωt=4θm/T
    to the effective value
    ω=θa/(T/4τ),
    where 2τ is the time corresponding to the nonlinear portion of the scanning function.

    From the five portions considered in Fig. 4(a), the scanning function for the triangular scanning regime is solved in Appendix A and described as follows:

    θ(t)={ωt,t[0,T/4τ)ω2τt2+ωT4τt+θmωT232τ,t[T4τ,T4+τ)ω(t+T/2),t[T/4+τ,3T/4τ)ω2τt23ωT4τtθm9ωT232τ,t[3T4τ,3T4+τ)ω(tT),t[3T/4+τ,T)

    From Fig. 4(a), the effective duty cycle of the triangular scanning function can be described as

    η=122τT.

    This effective duty cycle would double if a bidirectional scan was considered for the GS. By considering Eqs. (4) and (A9) derived in Appendix A, Eq. (6) becomes

    η=r/22r,
    where we used a notation similar to the one introduced in [24] with
    r=θa/θm.

  • (b) For the sawtooth input signals [Figs. 4(b1) and 4(b2)], the nonlinear stop-and-turn portions are characterized by two different effective amplitudes θa and θa that characterize the portions of the scanning function after the active scan and after the flyback, respectively. An assymetry of the time intervals of the different parts of the nonlinear portions also appears in comparison with the triangular scan in Fig. 4(a).

Two cases can be distinguished for sawtooth scanning:

  • (b1) For low fs and θm, the peaks of the θ(t) output/scanning function have the same positions as those of the input function [Fig. 4(b1)] and the output function is just beginning to curve due to inertia as fs increases, as shown in Fig. 2, but also in the detailed study of the functions in [22]. For low θm (i.e., 200–400 mV), the shape of the function θ(t) is preserved even at high fs [e.g., at 500 Hz, Figs. 2(c2) and 2(c3)]; the consequence is only the slight decrease in the effective duty cycle η as fs increases, as shown in Fig. 3.
  • (b2) For higher θm, two migrations Δt1 and Δt2 of the two peaks per period T of θ(t) occur [Fig. 4(b2)], as well as an increasing nonlinearity of the respective portions. This shift of the output signal of the GS is shown in Figs. 2(c2) and 2(c3) and it has also been documented in the media files provided in [22]. At a certain scan frequency fs, for each ηt the previous phenomenon leads to the complete shift of the sawtooth profile to the triangular one. The final profile is triangular with different nonlinear portions [Fig. 4(b2)], as determined experimentally in Fig. 2(c2), for example for θm equals 0.8 V at fs equals 500 Hz.

The mathematical discussion and the resulting model will be therefore different for these two cases:

  • (b1) For the graph in Fig. 4(b1), the theoretical angular velocity of the active scan of time interval ta=ηt·T is given as
    ωt=2θm/ηtT,
    while the effective angular velocity, for the effective active scan time ta=η·T is given as
    ω=θaηtT/2τ1=θaηtT/2τ2=θa+θaηtT/2τ1τ2.
    The theoretical angular velocity of the flyback is given as
    ωft=2θm/(1ηt)T,
    while the effective angular velocity of the flyback is given as
    ωf=θa(1ηt)T2τ1=θa(1ηt)T2τ2=(θa+θa)(1ηt)T2τ1τ2.

    The scanning function for the sawtooth input signal is solved in Appendix B for the five portions considered in Fig. 4(b1) and is described as follows:

    θ(t)={ωt,t[0,ηtT/2τ)ω2τ1(tηtT2)2+θm,t[ηtT2τ1,ηtT2)ωf2τ1(tηtT2)2+θm,t[ηtT2,ηtT2+τ1)ωf(tηtT/2τ1)+θa,t[ηtT/2+τ1,(1ηt/2)Tτ2)ωf2τ2[t(1ηt2)T]2θm,t[(1ηt2)Tτ2,(1ηt2)T)ω2τ2[t(1ηt2)T]2θm,t[(1ηt2)T,(1ηt2)T+τ2)ω(tT),t[(1ηt/2)T+τ2).

    From Eqs. (B4) and (B28), the relationship between the time intervals of the nonlinear portions is obtained as

    τ1τ1=τ2τ2=1ηtηt.

    The effective duty cycle of the scanning function in this case is

    η=taT=ηtτ1+τ2T,
    which, taking into account Eqs. (10), (B14), and (B28), becomes
    η=4θmωTηt.

  • (b2) For the graph in Fig. 4(b2), as the sawtooth switches to triangular, the effective angular velocity becomes
    ω=2θaΔ=2θaΔ=2(θa+θa)Δ+Δ,
    therefore
    Δ/Δ=θa/θa,
    while from Fig. 4(b2),
    Δ+Δ=T6τ.

The scanning function for the sawtooth input signal is solved in Appendix C for the five portions considered in Fig. 4(b2) and is described as

θ(t)={ωt,t[0,Δ)ω2τ(tΔ+τ)2+θm,t[Δ,Δ+2τ)ω(t+Δ+2τ)+θa,t[Δ+2τ,T/2+Δτ)ω4τ[t(T2+Δ+τ)]2θm,t[T2+Δτ,T2+Δ+3τ)ω(tT),t[T/2+Δτ,T).
From Eq. (2), by taking into account Eqs. (17)–(19), the effective duty cycle in this case is given as
η=Δ+ΔT=θa+θaωT
and, with Eqs. (A9) and (C12), it also becomes
η=2θm3τωT=123τT.

In conclusion to this analysis, by comparing the effective duty cycle η for (a) triangular scan, Eqs. (6), (b1) sawtooth scan, Eq. (15), and for (b2) sawtooth turning to a pseudo-triangular scan, Eq. (22), one can see the confirmation of the study in Figs. 3(d)3(f): with regard to ηt, η starts from immediately below 50% (from a lower value for a higher scan frequency), it increases for the “purely sawtooth scan” (b1), and then it decreases again—to a value closer to the triangular scan—for the case (b2), where the switch from sawtooth to triangular is completed. The latter situation occurs for a sufficiently high scan amplitude to have the above “switch” completed [i.e., the highest one in Figs. 3(d)3(f)].

4. OCT System and Scan Protocols

The layout of the OCT system that was employed in this part of the experimental study is shown in Fig. 5. This system is a custom-designed GD-OCM setup that provides a high volumetric resolution of 2 μm (in depth as well as lateral resolution) across a field of view of 2mm×2mm, and an imaging depth of 1.6 mm [27]. The source is a super luminescent diode (SLD) with the center wavelength of 840 nm and a full width at half-maximum of 100 nm (BroadLighter D-840-HP-I, Superlum, Ireland). A custom spectrometer with a high-speed CMOS line camera (spl4096-70 km, Basler Inc.) was implemented to acquire the spectral information [28]. The system allows for a real-time visualization of the sample after the acquisition, using a multiple graphic processing units (GPU) recently implemented in the system [29]. The scanning probe is equipped with a dual-axis XY galvanometer (6215HSM40B, Cambridge Technology).

 figure: Fig. 5.

Fig. 5. GD-OCM setup [29].

Download Full Size | PDF

To investigate the phenomena described in Section 2 on OCT imaging, GD-OCM images were acquired at several scan frequencies fs, from 50 to 600 Hz, including at around 300 Hz, where the differences between the scanning regimes that employ different theoretical duty cycles were easier to see. For contrast, two scan amplitudes θm were employed, a small one (200 mV) and a large one (3.2 V), with three theoretical duty cycles: 50%, 75%, and 90%. Note that the 50% duty cycle corresponds to a triangular scan.

In all the experiments, the line period of the camera was kept constant with Tc equal 15 μs. For each scan frequency, the number of samples (Ns) for one period of the scanner is given as

Ns=1/(Tc·fs).

The camera is configured to acquire data only during the rising time of the scanner or equivalently the increase in voltage. Thus, the number of A-scans to acquire is set by the product of Ns and the theoretical duty cycle. The fast axis (X) is following the protocol described above, while the slow axis (Y) moves step-by-step from one frame to another. The camera and the scanner are synchronized through an external TTL trigger.

5. OCT Image Distortion Due to the Inertia of the GS

In Sections 2 and 3, investigations show that the GS scanner may not exactly follow the input signal due to inertia. To further illustrate the inertia effect in OCT images, a sample with a regular grid structure is imaged using different scanning schemes. During the imaging process, different parts of the sample are imaged and then stitched together in order to obtain the entire image of the sample. Stitching of subimages is commonly done in many domains, including holography [30,31].

To the best of our knowledge, this aspect has been applied in OCT for transversal imaging of the retina with Doppler OCT [25], taking into account the results in [22]. The issue of the GS inertia discussed in Sections 2 and 3 needs to be further addressed in terms of numerical assesment, as imaging large samples has a multitude of applications, both medical and nonmedical—the latter, for example, in materials studies [32] or for art restauration [33].

For high scan frequencies and amplitudes (and especially when sawtooth scanning is used), the nonlinear portions of the scanning function shown in Fig. 2(c) are going to produce artifacts in the image of a regular stucture image—lines and even gaps for higher theoretical duty cycles, as shown in Fig. 6, in which results are shown for both scanning regimes: triangular and sawtooth. Several scan frequencies fs have been considered in the study, from 50 to 600 Hz. In Fig. 6, the images obtained for three values of fs are shown: a low value (100 Hz) in column 1, an intermediate value (300 Hz) in column 2, and the highest value considered in the study in Fig. 2 (i.e., 500 Hz) in column 3. The numbers of A-scans considered for half of each image are given in Table 1, for each fs and ηt. A number of 200 B-scans, with a step size of 3.65 μm, were taken for each 3D image.

 figure: Fig. 6.

Fig. 6. OCT imaging of a sample with a regular structure [frontal view of 3D OCT images shown as examples in Fig. 7(b)] with two scanning regimes: (a) triangular, (b) sawtooth with a theoretical duty cycle ηt equals 75%, and (c) sawtooth with ηt equals 90%. All images considered the same scan amplitude (θm equals 1.6 V, where 1 V stands for about 7.6 deg optically) and three different scan frequencies (fs): 100 Hz (a1), (b1), (c1); 300 Hz (a2), (b2), (c2); and 500 Hz (a3), (b3), (c3). For the images in column 3, the necessary overlaps are applied (Fig. 7) and the corrected images, without the previous artifacts, are shown in column 4. The dimensions of the images are 0.73mm×0.73mm.

Download Full Size | PDF

Tables Icon

Table 1. Effective Duty Cycle Calculated from the GD-OCM Experiments

For low fs, the images show no artifacts. For the two higher values of fs, the nonlinear portions of the scan functions produced blurred and distorted parts of the images related to the loss in duty cycle from ηt to η, as previously investigated. The number of distorted A-scans (Table 1) is proportional to the time intervals τ [Figs. 4(a) and 4(b2)] and τ1 [Fig. 4(b1)] for the respective cases. This is due to the fact that each A-scan is obtained in the OCT system in an identical time interval, set from the aquisition time of the camera, as pointed out in the previous section. In the same way, the total number of A-scans in an individual OCT image is proportional to half of the time period T of the GS—as also set from the scanning algorithm of the system.

Taking into account the number of distorted A-scans, which is proportional to the nonlinear time portions of the time intervals in Fig. 4, the equations for the effective duty cycle η deduced in Section 3 were used to obtain the effective duty cycles η that result from the OCT measurements. These values are given as the ηOCT [%] row of Table 1, which are then compared with the effective duty cycles that have been obtained in Section 2 from direct voltage measurements of the position of the galvomirror for the same values of fs, θm, and ηt as above. These latter values of ηGS were given in Figs. 3(c), 3(e), and 3(f).

The values of ηOCT are affected by measurement errors in τ and T because both these time intervals are determined with a possible error equal at most to Tc, which is the line period of the camera. The equations for these errors, denoted as εOCT [%], were deduced for each of the equations of ηOCT corresponding to various scanning regimes in Appendix D, as the following: Eqs. (D5), (D7), (D9) provide εOCT for Eqs. (6), (15), and (22), respectively. These εOCT errors are provided in Table 1. As it can be remarked, they are below 5% for each of the ηOCT values determined.

As the two sets of η values (i.e., ηOCT and ηGS) are compared in the last row of Table 1, the relative errors ε=(ηOCTηGS)/ηGS between these values are calculated for each value of fs and ηt.

One may note from this comparison that the mathematical model obtained for the triangular scan is always valid and is the most precise. However, from the values in Table 1, the mathematical models developed for the sawtooth regimes are also applicable, as discussed in Section 3, for the case in Fig. 4(b1) for a low fs (i.e., 50 Hz) and for the case in Fig. 4(b2) for the highest fs considered (i.e., 500 Hz).

The limits of the mathematical modeling can be noticed from the intermediate case of fs equal 300 Hz and ηt equal 75%: the two models developed cannot be applied, as the time intervals Δt1 and Δt2 pointed out in Fig. 4(b2) for the migrations of the peaks of the output function of the GS are not null (as for a low fs) and they are not maximum (as for the highest fs).

Another remark can be extracted from the experimental values of η obtained with the OCT experiments: the saturation of the effective duty cycle can be clearly seen for the 300 and 500 Hz scan frequencies (as studied and demonstrated from direct measurement of the scanning function of the GS in Section 2), as seen from the saturated values of ηGS in Table 1.

6. Distortion Correction

From the results shown in Fig. 6, columns 1–3, the artifacts in the images were evaluated. The effective duty cycle ηOCT was then calculated by using the appropriate mathematical model from those developed in Section 3 for each scan parameter ηt, fs and a scan amplitude θm equal to 1.6 V. The equation used for each case is indicated in Table 1. A comparison is then made with the effective duty cycle η, also determined experimentally, but from the galvoscanning in Section 2. One may see that a good match is obtained for all the cases for which a model could be applied.

To correct the distortion in the images, the effective duty cycles need to be taken into account when designing the scanning protocol. In Fig. 7(a1), the nonlinear part of the GS scanning is given as

Δl=2Lτ/T,
where L is the scanning range determined by the scanner amplitude. The nonlinear scanning part is where distortion occurs. To correct for the distorted part in the imaging acquisition, the camera was configured to acquire data only in the linear scanning part of the scanner. The scanning voltage for the second subimage was then shifted by a value corresponding to Δl to compensate the field of view loss due to the removal of distortion, which is shown in Fig. 7(a2).

 figure: Fig. 7.

Fig. 7. (a) Principle of the overlapping of the adjacent images for distortion correction. (b) The same study as in Fig. 6, but only for the sawtooth scanning with the maximum reachable (i.e., 90%) theoretical duty cycle, with GD-OCM volumetric images taken with fast scan, with the scan frequency fs equal to 500 Hz: (b1) original and (b2) corrected image—the latter obtained by overlapping two individual adjacent images. The dimensions of the images are 0.73mm×0.73mm×0.35mm.

Download Full Size | PDF

OCT images were obtained to show the distortion correction. In Fig. 6, column 3, the frontal images of the structured sample are shown, for three scan frequencies: 100 Hz (c1), 300 Hz (c2), and 500 Hz (c3). Only the worst case scenario, with sawtooth scan and ηt equals 90% is considered. The final images without the artifacts are shown, for each value of the fs, in Fig. 6, column 4. The portions where the correction was applied are marked on each figure.

In Fig. 7(b), the 3D reconstructions of the same large sample are shown. A number of 200 B-scans, with a step size of 3.65 μm were taken for each 3D image. The image obtained in Fig. 7(b1) for fs equals 500 Hz was corrected to remove the artifacts. The final image without the artifacts is shown in Fig. 7(b2).

7. Conclusions

We have investigated the output signals/scanning functions of galvanometer-based scanners (GSs)—when triangular and sawtooth input signals are considered. The profiles of the scanning functions (i.e., the current angular position of the galvomirror) were studied with regard to the main functional parameters of the process: scan frequency (fs), scan amplitude (θm), and theoretical/ideal duty cycle (ηt) of the input signal—the latter ranging from 50% (for triangular scanning) to a maximum value of 90% (for sawtooth scanning).

One of the main issues of the scanning process is precisely related to the effective duty cycle not equalling the theoretical duty cycle ηt, with an alteration increasing as fs is increased—especially at higher θm. The increased nonlinearity of the stop-and-turn portions of the scanning function thus decreased the duty cycle to an effective value η, which we studied with regard to fs and θm. Thus, for the fast GS of the imaging system the triangular scanning was demonstrated to be better than sawtooth scanning, as it produced images that were distortion-free up to higher values of fs and θm. However, this is valid only with OCT systems that utilize GSs that are programmed for bidirectional scan, as with the FD-OCT system we have used in [22].

It is often advantageous to use a unidirectional sawtooth scan, as the lateral resolution can thus be increased—more A-scans can be placed in the active time ta of the fast-axis GS. A tendency that GS users may have in this respect is to set the device to the limit in terms of the theoretical duty cycle, to take advantage of this feature. In this respect, a main finding of this study was that for sawtooth input signals a saturation of the scanning regime was reached. Although ηt is set towards its maximum, the effective duty cycle (i.e., the time efficiency of the scanning process) did not increase further; in fact, it was shown to even decrease in some cases (Fig. 3). This situation occurred when fs was increased, especially at higher θm, which is desired in order to have a field of view (FOV) that is as large as possible. Approximate limits of the parameters for this general phenomenon were indicated for fs roughly higher than 200 Hz (with differences that may appear with regard to the specific GS considered) and a θm close to the limit scan amplitude specific to the respective value of fs.

Mathematical models for both the output signal/scanning function and the effective duty cycle were deduced, based on the experimental findings on the galvoscanning. Both (a) the triangular and (b) the sawtooth scanning regimes were considered. For the latter type of scan, two subcases had to be considered: (b1) for low fs and θm, the peaks of the scanning function had the same positions as those of the input signal and the duty cycle was altered slightly due to inertia as fs increased; (b2) for higher θm, migrations of the two peaks per period T of the signals occurred, and at a certain fs (for each ηt) the sawtooth profile shifted completely to a triangular one with different nonlinear portions. This shift already manifested for fs higher then 300 Hz and it was final for fs higher then 500 Hz—for each θm close to the stability limit of the GS.

An OCT experiment was set to apply these findings for the optimization of the galvoscanning algorithms for Gabor domain optical coherence microscopy (GD-OCM). An in-house GD-OCM setup was used to image a regular multilayered structure for a larger FOV than the one provided by a single 2D scan—and 3D reconstruction—of the sample. Two individual scans were thus considered for a sufficiently high value of θm (i.e., 1.6 V, which means about 12 deg optically) and for a row of fs (from 50 to 600 Hz) in order to study the artifacts produced by scanning because of the nonlinear portions of the scanning functions and of the shifts of the sawtooth ones. The equations for the effective duty cycle η were deduced for each case and they were used to obtain η from the OCT measurements. These values were then compared in Table 1 with those obtained directly from the galvoscanning discussed in Section 2. A good agreement was obtained between the two sets of η values.

The artifacts produced in the OCT imaging (i.e., deformations of the regular structure imaged) were then used to perform the necessary stitching of the two individually scanned portions. The results was an artifact-free collated image, even at high scan frequencies fs (i.e., up to 600 Hz).

Future work is envisioned with the optimized scanning algorithms on other OCT and OCM applications, both biomedical and industrial [32], including with handheld scanning probes equipped with galvoscanners [17].

Appendix A: Effective Duty Cycle of a Triangular Scanning Function [Fig. 4(a)]

A parabolic (second-order polynomial) equation is considered for the nonlinear portions of the triangular scanning function:

θ(t)=at2+bt+c,wherea,b,c=cst.,
θ˙(t)=2at+b.
The coefficients a, b, and c have to be evaluated for each time interval to fulfill the continuity and smoothness conditions of the scanning function θ(t). From Fig. 4(a), the scanning function for each portion can be derived as follows:
  • (i) For t[0,T/4τ), the equation of this portion of the scanning function is given as
    θ(t)=ωt,
    where ω represents the effective angular velocity of the galvomirror—Eq. (4).
  • (ii) For t[T/4τ,T/4τ), the value of the scanning function can be evaluated at particular points to set the boundaries conditions as
    θ(T/4τ)=θ(T/4+τ)=θa,
    θ(T/4)=θm,
    θ˙(T/4τ)=θ˙(T/4τ)=0,
    θ˙(T/4)=0.
    From Eqs. (A1), (A2), and (A4)–(A7), the coefficients are given as
    a=ω/2τ;b=ωT/4τ;c=θmωT2/32τ,
    and a supplemental designing equation is also obtained as
    ωτ=2(θmθa).
  • (iii) For t[T/4+τ,3T/4τ), the equation of this portion is given as
    θ(t)=ω(t+T/2).
  • (iv) For t[3T/4τ,3T/4+τ), the boundaries conditions of θ(t) are given as
    θ(3T/4τ)=θ(3T/4+τ)=θa,
    θ(3T/4)=θm,
    θ˙(3T/4τ)=θ˙(3T/4+τ)=ω,
    θ˙(3T/4)=0.
    From Eqs. (A1), (A2), and (A11) to (A14), the coefficients are given as
    a=ω/2τ,b=3ωT/4τ,c=θm9ωT2/32τ,
    and Eq. (A9) is also obtained as a supplemental equation.
  • (v) For t[3T/4+τ,T), this portion of the function is given as
    θ(t)=ω(tT).

Appendix B: Effective Duty Cycle of a Sawtooth Scanning Function at Low Values of the Scan Frequency [Fig. 4(b1)]

For the nonlinear portions of the scanning function, the same parabolic functions in Eq. (A1) will be considered.

  • (i) For t[0,ηtT/2τ1), we may write
    θ(t)=ωt.
  • (ii) For t[ηtT/2τ1,ηtT/2), the a, b, c constants are obtained from the conditions given as [Fig. 4(b1)]
    θ(ηtT/2τ1)=θa,
    θ(ηtT/2)=θm,
    θ˙(ηtT/2τ1)=0,
    θ˙(ηtT/2)=0,
    which were written for the continuity and smoothness of the scanning function θ(t). From Eq. (A1), using Eqs. (B2)–(B5), the three coefficients are given as
    a=ω/2τ1,b=ωηtT/2τ1,c=θmωηt2T2/8τ1,
    and a designing equation is also obtained as
    ωτ1=2(θmθa).
  • (iii) For t[ηtT/2,ηtT/2+τ1), the continuity and smoothness conditions are [Fig. 4(b1)] expressed as
    θ(ηtT/2)=θm,
    θ(ηtT/2+τ1)=θa,
    θ˙(ηtT/2)=0,
    θ˙(ηtT/2+τ1)=ωf,
    which were written for the continuity and smoothness of the scanning function θ(t). From Eq. (A1), using Eqs. (B8)–(B11), the three coefficients are then given as
    a=ωf/2τ1,b=ωfηtT/2τ1,c=θm+ωfηt2T2/8τ1,
    and a designing equation is also obtained as
    ωfτ1=2(θmθa).
    From Eqs. (B7) and (B13), taking into account Eqs. (10) and (12), a relationship between the time intervals of the previous portion was obtained as
    τ1τ1=ωfω=ηt1ηt.
  • (iv) For t[ηtT/2+τ1,(1ηt/2)Tτ2) the equation of this portion is given as
    θ(t)=ωf(tηtT/2τ1)+θa.
  • (v) For t[(1ηt/2)Tτ2),(1ηt/2)T) the continuity and smoothness conditions are [Fig. 4(b1)] given as
    θ((1ηt/2)Tτ2)=θa,
    θ((1ηt/2)T)=θm,
    θ˙((1ηt/2)Tτ2)=ωf,
    θ˙((1ηt/2)T)=0,
    and the coefficients of the functions are expressed as
    a=ωf2τ2,b=ωf(2ηt)T2τ2,c=θmωf(1ηt)T22τ2,
    where a designing equation is also obtained as
    ωfτ2=2(θmθa).
  • (vi) For t[(1ηt/2)T,(1ηt/2)T+τ2) the continuity and smoothness conditions are [Fig. 4(b1)] given as
    θ((1ηt/2)T+τ2)=θa,
    θ((1ηt/2)T)=θm,
    θ˙(((1ηt/2)T+τ2))=ω,
    θ˙((1ηt/2)T)=0,
    and the coefficients of the functions are expressed as
    a=ω2τ2,b=ω(2ηt)T2τ2,c=θm+ω(1ηt/2)T22τ2,
    where a designing equation is also obtained as
    ωτ2=2(θmθa).

    From Eqs. (B21) and (B27), taking into account Eqs. (10) and (12), a relationship between the time intervals of the previous portions is obtained as

    τ2τ2=ωωf=1ηtηt,

  • (vii) For t[1ηt/2)T+τ2,T) Eq. (A16) is valid in this case as well.

Appendix C: Effective Duty Cycle of a Sawtooth Scanning Function at High Values of the Scan Frequency [Fig. 4(b2)]

For the nonlinear portions of the scanning function, the same parabolic functions in Eq. (A1) will be considered.

  • (i) For t[0,Δ), Eq. (B1) is valid in this case as well.
  • (ii) For t[Δ,Δ+2τ), the continuity and smoothness conditions of θ(t) are [Fig. 4(b2)] given as
    θ(Δ)=θ(Δ+2τ)=θa,
    θ(Δ+τ)=θm,
    θ˙(Δ)=θ˙(Δ+2τ)=ω,
    θ˙(Δ+τ)=0.
    From Eqs. (A1), (A2), and (C1) to (C4), the coefficients are expressed as
    a=ω/2τ,b=ω(Δ+τ)/τ,c=θmω(Δ+τ)2/2τ,
    and Eq. (A9) is also obtained as a supplemental equation.
  • (iii) For t[Δ+2τ,Δ+T/2τ) the equation of this portion is given as
    θ(t)=ω(t+Δ+2τ)+θa.
  • (iv) For t[Δ+T/2τ,T/2+Δ+3τ), the continuity and smoothness conditions of θ(t) are [Fig. 4(b2)] given as
    θ(T/2+Δτ)=θ(T/2+Δ+3τ)=θa,
    θ(T/2+Δ+τ)=θm,
    θ˙(T/2+Δτ)=θ˙(T/2+Δ+3τ)=ω,
    θ˙(T/2+Δ+τ)=0.
    From Eqs. (A1), (A2), and (C7)–(C10), the coefficients are expressed as
    a=ω4τ,b=ω2τ(T2+Δ+τ),c=ω4τ(T2+Δ+τ)2θm,
    and a designing equation is also obtained as
    ωτ=θmθa.
  • (v) For t[T/2+Δ+3τ,T) Eq. (A16) is valid in this case as well.

Appendix D: Measuring Errors of the Effective Duty Cycle of Triangular and Sawtooth Scanning Functions

The effective duty cycle η for triangular scanning functions is provided by Eq. (6), while for sawtooth scanning functions by Eqs. (15) or (22), as discussed in Section 3. The maximum measuring error of each of the τ and T time intervals equals the line period of the camera Tc (equals 15 μs, see Section 4). The absolute measuring error of η is

Δη=ητΔτ+ηTΔT,
and with Δτ=ΔT=TC
Δη=(ητ+ηT)TC,
from which the relative error in measuring ηOCT is
ε=|Δηη|.

From Eqs. (D2) and (D3), taking into account that 1/T=fs, the absolute and relative errors are, respectively,

  • – for the triangular scan [Eq. (6)],
    Δη=TCT(32+η),
    and
    ε=(1+1.5η)TCfs.
  • – For the sawtooth scan [Eq. (15)], consider for simplicity that τ1τ2τ, then
    Δη=TCT(2ηt+η),
    and
    ε=(1+2ηtη)TCfs.
  • – For the sawtooth scan [Eq. (15)], if we consider that τ1τ2τ, then
    Δη=TCT(52+η),
    and
    ε=(1+2.5η)TCfs.

FUNDING INFORMATION

II-VI Foundation; National Science Foundation (NSF) (ENG 1346453); Romanian National Authority for Scientific Research (PN-II-PT-PCCA-2011-3.2-1682).

References

1. W. Drexler, M. Liu, A. Kumar, T. Kamali, A. Unterhuber, and R. A. Leitgeb, “Optical coherence tomography today: speed, contrast, and multimodality,” J. Biomed. Opt. 19, 071412 (2014). [CrossRef]  

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

3. A. Gh. Podoleanu and R. B. Rosen, “Combinations of techniques in imaging the retina with high resolution,” Prog. Retinal Eye Res. 27, 464–499 (2008). [CrossRef]  

4. G. F. Marshall, Handbook of Optical and Laser Scanning (CRC Press, 2011).

5. M. Bass, Handbook of Optics, 3rd ed. (McGraw-Hill, 2009), pp. 30.1–30.68.

6. Y. Li, “Beam deflection and scanning by two-mirror and two-axis systems of different architectures: a unified approach,” Appl. Opt. 47, 5976–5985 (2008). [CrossRef]  

7. K. H. Kim, C. Buehler, and P. T. C. So, “High-speed, two-photon scanning microscope,” Appl. Opt. 38, 6004–6009 (1999). [CrossRef]  

8. Y. Li, “Third-order theory of the Risley-prism-based beam steering system,” Appl. Opt. 50, 679–686 (2011). [CrossRef]  

9. A. Schitea, M. Tuef, and V. F. Duma, “Modeling of Risley prisms devices for exact scan patterns,” Proc. SPIE 8789, 878912 (2013).

10. D. Miyazaki, K. Shiba, K. Sotsuka, and K. Matsushita, “Volumetric display system based on three dimensional scanning of inclined optical image,” Opt. Express 14, 12760–12769 (2006). [CrossRef]  

11. X. Tao, H. Cho, and F. Janabi-Sharifi, “Optical design of a variable view imaging system with the combination of a telecentric scanner and double wedge prisms,” Appl. Opt. 49, 239–246 (2010). [CrossRef]  

12. M. W. Jenkins, O. Q. Chughtai, A. N. Basavanhally, M. Watanabe, and A. M. Rollins, “In vivo gated 4D imaging of the embryonic heart using optical coherence tomography,” J. Biomed. Opt. 12, 030505 (2007). [CrossRef]  

13. A. Gh. Podoleanu, G. M. Dobre, and D. A. Jackson, “En-face coherence imaging using galvanometer scanner modulation,” Opt. Lett. 23, 147–149 (1998). [CrossRef]  

14. C. C. Rosa, J. Rogers, and A. Gh. Podoleanu, “Fast scanning transmissive delay line for optical coherence tomography,” Opt. Lett. 30, 3263–3265 (2005). [CrossRef]  

15. B. Baumann, M. Pircher, E. Götzinger, and Ch. K. Hitzenberger, “Full range complex spectral domain optical coherence tomography without additional phase shifters,” Opt. Express 15, 13375–13387 (2007). [CrossRef]  

16. R. Cernat, T. S. Tatla, J. Pang, P. J. Tadrous, A. Bradu, G. Dobre, G. Gelikonov, V. Gelikonov, and A. Gh. Podoleanu, “Dual instrument for in vivo and ex vivo OCT imaging in an ENT department,” Biomed. Opt. Express 3, 3346–3356 (2012). [CrossRef]  

17. D. Demian, V. F. Duma, C. Sinescu, M. L. Negrutiu, R. Cernat, F. I. Topala, Gh. Hutiu, A. Bradu, and A. Gh. Podoleanu, “Design and testing of prototype handheld scanning probes for optical coherence tomography,” J. Eng. Med. 228, 743–753 (2014).

18. V. F. Duma, J. P. Rolland, and A. Gh. Podoleanu, “Perspectives of optical scanning in OCT,” Proc. SPIE 7556, 75560B (2010).

19. R. P. Aylward, “Advances and technologies of galvanometer-based optical scanners,” Proc. SPIE 3787, 158–164 (1999).

20. J. S. Gadhok, “Achieving high-duty cycle sawtooth scanning with galvanometric scanners,” Proc. SPIE 3787, 173–180 (1999).

21. J. Montagu, “Scanners—galvanometric and resonant,” in Encyclopedia of Optical Engineering, R. G. Driggers, C. Hoffman, and R. Driggers, eds. (Taylor & Francis, 2003), pp. 2465–2487. [CrossRef]  

22. V. F. Duma, K.-S. Lee, P. Meemon, and J. P. Rolland, “Experimental investigations of the scanning functions of galvanometer-based scanners with applications in OCT,” Appl. Opt. 50, 5735–5749 (2011). [CrossRef]  

23. C. D. Lu, M. F. Kraus, B. Potsaid, J. J. Liu, W. Choi, V. Jayaraman, A. E. Cable, J. Hornegger, J. S. Duker, and J. G. Fujimoto, “Handheld ultrahigh speed swept source optical coherence tomography instrument using a MEMS scanning mirror,” Biomed. Opt. Express 5, 293–311 (2014). [CrossRef]  

24. V. F. Duma, “Optimal scanning function of a galvanometer scanner for an increased duty cycle,” Opt. Eng. 49, 103001 (2010). [CrossRef]  

25. B. Braaf, K. A. Vermeer, K. V. Vienola, and J. F. de Boer, “Angiography of the retina and the choroid with phase-resolved OCT using interval-optimized backstitched B-scans,” Opt. Express 20, 20516–20534 (2012). [CrossRef]  

26. C. Mnerie and V. F. Duma, “Mathematical model of a galvanometer-based scanner: simulations and experiments,” Proc. SPIE 8789, 878915 (2013).

27. K. Lee, K. P. Thompson, P. Meemon, and J. P. Rolland, “Cellular resolution optical coherence microscopy with high acquisition speed for in-vivo human skin volumetric imaging,” Opt. Lett. 36, 2221–2223 (2011). [CrossRef]  

28. K. Lee, K. P. Thompson, and J. P. Rolland, “Broadband astigmatism-corrected Czerny–Turner spectrometer,” Opt. Express 18, 23378–23384 (2010). [CrossRef]  

29. P. Tankam, A. P. Santhanan, K. Lee, J. Won, C. Canavesi, and J. P. Rolland, “Parallelized multi-graphics processing unit framework for high-speed Gabor-domain optical coherence microscopy,” J. Biomed. Opt. 19, 071410 (2014). [CrossRef]  

30. A. Pelagotti, M. Paturzo, M. Locatelli, A. Geltrude, R. Meucci, A. Finizio, and P. Ferraro, “An automatic method for assembling a large synthetic aperture digital hologram,” Opt. Express 20, 4830–4839 (2012). [CrossRef]  

31. P. Picart and P. Tankam, “Analysis and adaptation of convolution algorithms to reconstruct extended objects in digital holography,” Appl. Opt. 52, A240–A253 (2013). [CrossRef]  

32. Gh. Hutiu, V. F. Duma, D. Demian, A. Bradu, and A. Gh. Podoleanu, “Surface imaging of metallic material fractures using optical coherence tomography,” Appl. Opt. 53, 5912–5916 (2014). [CrossRef]  

33. H. Liang, M. G. Cid, R. G. Cucu, G. M. Dobre, A. Gh. Podoleanu, J. Pedro, and D. Saunders, “En-face optical coherence tomography—a novel application of non-invasive imaging to art conservation,” Opt. Express 13, 6133–6144 (2005). [CrossRef]  

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 (7)

Fig. 1.
Fig. 1. Schematic illustrating the operating principle of a galvanometer-based scanner (GS), with constructive parameters: J , moment of inertia of the mobile element (with galvomirror); c , damping coefficient; k , elastic coefficient of the torsion springs that support the mobile element.
Fig. 2.
Fig. 2. Study of sawtooth scanning functions/output signals (i.e., angular positions of the galvomirror) of a GS for input signals with different theoretical duty cycles η t equal to (a1) 50% (triangular function), (a2) 75%, and (a3) 87.5%; (b1)–(b3) output signals (positions of the galvomirror) for a scan frequency f s = 50 Hz ; (c1)–(c3) output signals for f s = 500 Hz ; (d1)–(d3) effective duty cycle of the GS with regard to f s for three different driving duty cycles and different scan amplitudes θ m of 0.2, 0.4, 0.8, 1.6, and 3.2 V (where 1 V corresponds to about 7.6 deg of optical scan angle).
Fig. 3.
Fig. 3. Study of the effective duty cycle of a GS with regard to the theoretical/programmed duty cycle, for different scan amplitudes θ m of 0.2, 0.4, 0.8,1.6, and 3.2 V (1 V corresponds to 7.6 deg optically in terms of angular scan amplitude)—each study being made for a certain scan frequency f s : (a) 10 Hz; (b) 50 Hz; (c) 100 Hz; (d) 200 Hz; (e) 300 Hz; (f) 500 Hz.
Fig. 4.
Fig. 4. (a) Triangular and (b1), (b2) sawtooth scanning regimes of the GS: ideal/input signals versus scanning functions/output signals of the GS. For the latter, the current positions of the galvomirror are those determined experimentally and shown in Fig. 2.
Fig. 5.
Fig. 5. GD-OCM setup [29].
Fig. 6.
Fig. 6. OCT imaging of a sample with a regular structure [frontal view of 3D OCT images shown as examples in Fig. 7(b)] with two scanning regimes: (a) triangular, (b) sawtooth with a theoretical duty cycle η t equals 75%, and (c) sawtooth with η t equals 90%. All images considered the same scan amplitude ( θ m equals 1.6 V, where 1 V stands for about 7.6 deg optically) and three different scan frequencies ( f s ): 100 Hz (a1), (b1), (c1); 300 Hz (a2), (b2), (c2); and 500 Hz (a3), (b3), (c3). For the images in column 3, the necessary overlaps are applied (Fig. 7) and the corrected images, without the previous artifacts, are shown in column 4. The dimensions of the images are 0.73 mm × 0.73 mm .
Fig. 7.
Fig. 7. (a) Principle of the overlapping of the adjacent images for distortion correction. (b) The same study as in Fig. 6, but only for the sawtooth scanning with the maximum reachable (i.e., 90%) theoretical duty cycle, with GD-OCM volumetric images taken with fast scan, with the scan frequency f s equal to 500 Hz: (b1) original and (b2) corrected image—the latter obtained by overlapping two individual adjacent images. The dimensions of the images are 0.73 mm × 0.73 mm × 0.35 mm .

Tables (1)

Tables Icon

Table 1. Effective Duty Cycle Calculated from the GD-OCM Experiments

Equations (89)

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

η t = t a / T .
η = t a / T .
ω t = 4 θ m / T
ω = θ a / ( T / 4 τ ) ,
θ ( t ) = { ω t , t [ 0 , T / 4 τ ) ω 2 τ t 2 + ω T 4 τ t + θ m ω T 2 32 τ , t [ T 4 τ , T 4 + τ ) ω ( t + T / 2 ) , t [ T / 4 + τ , 3 T / 4 τ ) ω 2 τ t 2 3 ω T 4 τ t θ m 9 ω T 2 32 τ , t [ 3 T 4 τ , 3 T 4 + τ ) ω ( t T ) , t [ 3 T / 4 + τ , T )
η = 1 2 2 τ T .
η = r / 2 2 r ,
r = θ a / θ m .
ω t = 2 θ m / η t T ,
ω = θ a η t T / 2 τ 1 = θ a η t T / 2 τ 2 = θ a + θ a η t T / 2 τ 1 τ 2 .
ω f t = 2 θ m / ( 1 η t ) T ,
ω f = θ a ( 1 η t ) T 2 τ 1 = θ a ( 1 η t ) T 2 τ 2 = ( θ a + θ a ) ( 1 η t ) T 2 τ 1 τ 2 .
θ ( t ) = { ω t , t [ 0 , η t T / 2 τ ) ω 2 τ 1 ( t η t T 2 ) 2 + θ m , t [ η t T 2 τ 1 , η t T 2 ) ω f 2 τ 1 ( t η t T 2 ) 2 + θ m , t [ η t T 2 , η t T 2 + τ 1 ) ω f ( t η t T / 2 τ 1 ) + θ a , t [ η t T / 2 + τ 1 , ( 1 η t / 2 ) T τ 2 ) ω f 2 τ 2 [ t ( 1 η t 2 ) T ] 2 θ m , t [ ( 1 η t 2 ) T τ 2 , ( 1 η t 2 ) T ) ω 2 τ 2 [ t ( 1 η t 2 ) T ] 2 θ m , t [ ( 1 η t 2 ) T , ( 1 η t 2 ) T + τ 2 ) ω ( t T ) , t [ ( 1 η t / 2 ) T + τ 2 ) .
τ 1 τ 1 = τ 2 τ 2 = 1 η t η t .
η = t a T = η t τ 1 + τ 2 T ,
η = 4 θ m ω T η t .
ω = 2 θ a Δ = 2 θ a Δ = 2 ( θ a + θ a ) Δ + Δ ,
Δ / Δ = θ a / θ a ,
Δ + Δ = T 6 τ .
θ ( t ) = { ω t , t [ 0 , Δ ) ω 2 τ ( t Δ + τ ) 2 + θ m , t [ Δ , Δ + 2 τ ) ω ( t + Δ + 2 τ ) + θ a , t [ Δ + 2 τ , T / 2 + Δ τ ) ω 4 τ [ t ( T 2 + Δ + τ ) ] 2 θ m , t [ T 2 + Δ τ , T 2 + Δ + 3 τ ) ω ( t T ) , t [ T / 2 + Δ τ , T ) .
η = Δ + Δ T = θ a + θ a ω T
η = 2 θ m 3 τ ω T = 1 2 3 τ T .
N s = 1 / ( T c · f s ) .
Δ l = 2 L τ / T ,
θ ( t ) = a t 2 + b t + c , where a , b , c = cst. ,
θ ˙ ( t ) = 2 a t + b .
θ ( t ) = ω t ,
θ ( T / 4 τ ) = θ ( T / 4 + τ ) = θ a ,
θ ( T / 4 ) = θ m ,
θ ˙ ( T / 4 τ ) = θ ˙ ( T / 4 τ ) = 0 ,
θ ˙ ( T / 4 ) = 0 .
a = ω / 2 τ ; b = ω T / 4 τ ; c = θ m ω T 2 / 32 τ ,
ω τ = 2 ( θ m θ a ) .
θ ( t ) = ω ( t + T / 2 ) .
θ ( 3 T / 4 τ ) = θ ( 3 T / 4 + τ ) = θ a ,
θ ( 3 T / 4 ) = θ m ,
θ ˙ ( 3 T / 4 τ ) = θ ˙ ( 3 T / 4 + τ ) = ω ,
θ ˙ ( 3 T / 4 ) = 0 .
a = ω / 2 τ , b = 3 ω T / 4 τ , c = θ m 9 ω T 2 / 32 τ ,
θ ( t ) = ω ( t T ) .
θ ( t ) = ω t .
θ ( η t T / 2 τ 1 ) = θ a ,
θ ( η t T / 2 ) = θ m ,
θ ˙ ( η t T / 2 τ 1 ) = 0 ,
θ ˙ ( η t T / 2 ) = 0 ,
a = ω / 2 τ 1 , b = ω η t T / 2 τ 1 , c = θ m ω η t 2 T 2 / 8 τ 1 ,
ω τ 1 = 2 ( θ m θ a ) .
θ ( η t T / 2 ) = θ m ,
θ ( η t T / 2 + τ 1 ) = θ a ,
θ ˙ ( η t T / 2 ) = 0 ,
θ ˙ ( η t T / 2 + τ 1 ) = ω f ,
a = ω f / 2 τ 1 , b = ω f η t T / 2 τ 1 , c = θ m + ω f η t 2 T 2 / 8 τ 1 ,
ω f τ 1 = 2 ( θ m θ a ) .
τ 1 τ 1 = ω f ω = η t 1 η t .
θ ( t ) = ω f ( t η t T / 2 τ 1 ) + θ a .
θ ( ( 1 η t / 2 ) T τ 2 ) = θ a ,
θ ( ( 1 η t / 2 ) T ) = θ m ,
θ ˙ ( ( 1 η t / 2 ) T τ 2 ) = ω f ,
θ ˙ ( ( 1 η t / 2 ) T ) = 0 ,
a = ω f 2 τ 2 , b = ω f ( 2 η t ) T 2 τ 2 , c = θ m ω f ( 1 η t ) T 2 2 τ 2 ,
ω f τ 2 = 2 ( θ m θ a ) .
θ ( ( 1 η t / 2 ) T + τ 2 ) = θ a ,
θ ( ( 1 η t / 2 ) T ) = θ m ,
θ ˙ ( ( ( 1 η t / 2 ) T + τ 2 ) ) = ω ,
θ ˙ ( ( 1 η t / 2 ) T ) = 0 ,
a = ω 2 τ 2 , b = ω ( 2 η t ) T 2 τ 2 , c = θ m + ω ( 1 η t / 2 ) T 2 2 τ 2 ,
ω τ 2 = 2 ( θ m θ a ) .
τ 2 τ 2 = ω ω f = 1 η t η t ,
θ ( Δ ) = θ ( Δ + 2 τ ) = θ a ,
θ ( Δ + τ ) = θ m ,
θ ˙ ( Δ ) = θ ˙ ( Δ + 2 τ ) = ω ,
θ ˙ ( Δ + τ ) = 0 .
a = ω / 2 τ , b = ω ( Δ + τ ) / τ , c = θ m ω ( Δ + τ ) 2 / 2 τ ,
θ ( t ) = ω ( t + Δ + 2 τ ) + θ a .
θ ( T / 2 + Δ τ ) = θ ( T / 2 + Δ + 3 τ ) = θ a ,
θ ( T / 2 + Δ + τ ) = θ m ,
θ ˙ ( T / 2 + Δ τ ) = θ ˙ ( T / 2 + Δ + 3 τ ) = ω ,
θ ˙ ( T / 2 + Δ + τ ) = 0 .
a = ω 4 τ , b = ω 2 τ ( T 2 + Δ + τ ) , c = ω 4 τ ( T 2 + Δ + τ ) 2 θ m ,
ω τ = θ m θ a .
Δ η = η τ Δ τ + η T Δ T ,
Δ η = ( η τ + η T ) T C ,
ε = | Δ η η | .
Δ η = T C T ( 3 2 + η ) ,
ε = ( 1 + 1.5 η ) T C f s .
Δ η = T C T ( 2 η t + η ) ,
ε = ( 1 + 2 η t η ) T C f s .
Δ η = T C T ( 5 2 + η ) ,
ε = ( 1 + 2.5 η ) T C f s .
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.