Extending the Effective Ranging Depth of Spectral Domain Optical Coherence Tomography by Spatial Frequency Domain Multiplexing

We present a spatial frequency domain multiplexing method for extending the imaging depth range of a spectral domain optical coherence tomography (SDOCT) system without any expensive device. This method uses two galvo scanners with different pivot-offset distances in two independent reference arms for spatial frequency modulation and multiplexing. The spatial frequency contents corresponding to different depth regions of the sample can be shifted to different frequency bands. The spatial frequency domain multiplexing SDOCT system provides an approximately 1.9-fold increase in the effective ranging depth compared with that of a conventional full-range SDOCT system. The reconstructed images of phantom and biological tissue demonstrate the expected increase in ranging depth. The parameters choice criterion for this method is discussed.


Introduction
Fourier domain optical coherence tomography (FDOCT) is a high-speed, high-resolution and non-invasive biomedical optical imaging technique [1].It has been applied in the field of biomedical research, clinical diagnostics, and non-destructive material inspection.FDOCT can be categorized into two embodiments.One is based on a wideband laser source and a spectrometer which is called spectral domain OCT (SDOCT), and the other one is based on a frequency sweeping laser source and a balanced photon detector which is called swept source OCT (SSOCT, also called optical frequency domain imaging, OFDI).For both types of the FDOCT, the achievable imaging depth range is an important performance indicator.In some situations, a long ranging depth is highly desirable.One such example is for the endoscopic application where the distance between the distal end of the endoscopic probe and the tissue surface cannot be accurately controlled.Another example is for imaging the whole eye from anterior to retina.
For FDOCT, the ranging depth can be calculated using the formula ∆z = λ 0 2 /4nδλ and is fundamentally limited by the center wavelength λ 0 of the light source, the spectral resolution δλ, and the refractive index n.The spectral resolution is δλ = ∆λ/N when the spectrometer is pixel-limited.In this equation, ∆λ represents the spectral bandwidth of the light source and N is the data number of the digitized axial spectral interference signal that covers the bandwidth.In fact, there is a sensitivity roll-off issue for FDOCT, which poses a further limitation to the practical imaging depth range.For the case of SDOCT, the sensitivity roll-off effect is due to the finite pixel number of the line scanning camera.

Spatial Frequency Domain Multiplexing Spectral Domain Optical Coherence Tomography (SDOCT) System
Figure 1 depicts the schematic of the spatial frequency domain multiplexing SDOCT system.It comprises a super luminescent diode (SLD), an interferometer with two pivot-offset galvo scanners implemented in the two separate reference arms, a custom-built spectrometer, and a personal computer for signal processing and image display.Reference Arms I and II have different round-trip delays, which correspond to the shallow and deep regions of the sample, shown as Z A and Z B in Figure 1.Thus, two sets of the spectral interferogram corresponding to the depths Z A and Z B can be acquired.In the sample arm, the light beam just incidents upon the pivot of the GS0 for sample scanning.The light beams in the reference arms incident upon the GS1 and GS2 with the different offset distances s 1 and s 2 away from the pivots, respectively, as shown in Figure 1.Appl. Sci.2016, 6, 360 3 of 11

Spatial Frequency Domain Multiplexing Spectral Domain Optical Coherence Tomography (SDOCT) System
Figure 1 depicts the schematic of the spatial frequency domain multiplexing SDOCT system.It comprises a super luminescent diode (SLD), an interferometer with two pivot-offset galvo scanners implemented in the two separate reference arms, a custom-built spectrometer, and a personal computer for signal processing and image display.Reference Arms I and II have different round-trip delays, which correspond to the shallow and deep regions of the sample, shown as ZA and ZB in Figure 1.Thus, two sets of the spectral interferogram corresponding to the depths ZA and ZB can be acquired.In the sample arm, the light beam just incidents upon the pivot of the GS0 for sample scanning.The light beams in the reference arms incident upon the GS1 and GS2 with the different offset distances s1 and s2 away from the pivots, respectively, as shown in Figure 1.When the galvo scanners in the sample arm and reference arms are driven simultaneously, the linear modulated phase will be added into the two-dimensional (2D) spectral interferograms.By ignoring the signals that do not contribute to the depth localization in the sample for simplicity, the 2D spectral interferogram detected at each wavenumber k and lateral position x can be expressed as where k represents wavenumber and  When the galvo scanners in the sample arm and reference arms are driven simultaneously, the linear modulated phase will be added into the two-dimensional (2D) spectral interferograms.By ignoring the signals that do not contribute to the depth localization in the sample for simplicity, the 2D spectral interferogram detected at each wavenumber k and lateral position x can be expressed as where k represents wavenumber and k = 2π/λ, z is the depth position coordinate of the scatters within the sample, n is the average refractive index of the sample, R (x, z) is the normalized intensity backscattered from the scatter at z, S(k) is the spectral density of the light source, Φ 1 (x) and Φ 2 (x) represent the linear modulated phase induced by GS1 and GS2, and φ A (x, z) and φ B (x, z) are the phases that relate to the optical heterogeneity of the sample in the focus coherence volume.If the linear modulation phase Φ 1 (x) and Φ 2 (x) have different slopes about the lateral position, the spatial frequency contents of the 2D OCT spectral interferogram corresponding to the depths Z A and Z B will be shifted to the different center spatial frequencies away from the zero spatial frequency.

Demodulation Algorithm for Spatial Frequency Domain Multiplexing SDOCT
As shown in Figure 2, after taking the Fourier transform of the 2D spectral interferogram I(x, k) along the x-direction (step (i) in Figure 2), we obtain the spatial spectrum of the 2D spectral interferogram, which can be expressed as where ν is the Fourier transform pair of x.C(ν, k) and Ĉ(ν, k) are the pair of Hermit conjugate term due to the Fourier transform of the real valued interference signal and contain the same information on the OCT signal in the spatial frequency domain.The subscripts A and B in Equation ( 2) indicate the spatial frequency contents corresponding to depths Z A and Z B , respectively.ν c1 and ν c2 are the center spatial carrier frequencies, which are proportional to the pivot-offset distances.For reconstruction of the complex valued spectral interferogram, two band-pass filters with center frequencies of ν c1 and ν c2 are applied to the spatial frequency spectrum (step (ii) in Figure 2).We then obtain the separated spatial frequency content C A (ν + ν c1 , k) and C B (ν + ν c2 , k). spatial frequency contents of the 2D OCT spectral interferogram corresponding to the depths ZA and ZB will be shifted to the different center spatial frequencies away from the zero spatial frequency.

Demodulation Algorithm for Spatial Frequency Domain Multiplexing SDOCT
As shown in Figure 2, after taking the Fourier transform of the 2D spectral interferogram ( , ) I x k along the x-direction (step (i) in Figure 2), we obtain the spatial spectrum of the 2D spectral interferogram, which can be expressed as ( , ) where ν is the Fourier transform pair of x.  2) indicate the spatial frequency contents corresponding to depths ZA and ZB, respectively.νc1 and νc2 are the center spatial carrier frequencies, which are proportional to the pivot-offset distances.For reconstruction of the complex valued spectral interferogram, two band-pass filters with center frequencies of νc1 and νc2 are applied to the spatial frequency spectrum (step (ii) in Figure 2).We then obtain the separated spatial frequency content  2).For reasons of energy conservation, the data have to be scaled by a factor of 2.Then, after the k space re-sampling the two full-range OCT images corresponding to the depths of ZA and ZB can be calculated by inversely Two complex-valued spectral interferogram I A (x, k) and I B (x, k) can now be reconstructed by performing the inverse Fourier transform from the spatial frequency domain back to the spatial domain for the terms C A (ν + ν c1 , k) and C B (ν + ν c2 , k) (step (iii) in Figure 2).For reasons of energy conservation, the data have to be scaled by a factor of 2.Then, after the k space re-sampling the two full-range OCT images corresponding to the depths of Z A and Z B can be calculated by inversely Fourier transforming the complex-valued linear-in-k spectral interferogram I A (x, k) and I B (x, k) along the wavenumber direction.Finally, an extended range OCT image is formed by concatenating the image regions that correspond to the effective ranging depths of the two full-range OCT images.

The Center Spatial Frequency and Parameters Design Criterion
where ω is the angular velocity of the galvo scanners, T is the exposure time of the camera in the spectrometer, and λ 0 is the center wavelength of the light source.along the wavenumber direction.Finally, an extended range OCT image is formed by concatenating the image regions that correspond to the effective ranging depths of the two full-range OCT images.

The Center Spatial Frequency and Parameters Design Criterion
Where ω is the angular velocity of the galvo scanners, T is the exposure time of the camera in the spectrometer, and 0 λ is the center wavelength of the light source.Considering a tomogram of N depth profiles covering a lateral range of Δx, the two center spatial carrier frequencies can be calculated by The lateral scanning range Δx is calculated by where f is the focal length of the sample objective.
The positive spatial frequency range can be calculated by For high-quality image reconstruction, another important parameter is the characteristic bandwidth (BW) of the obtained spatial frequency contents, which is determined by the inverse speckle size.The speckle size can be modeled by the transverse resolution, and BW can be expressed as [9]  Considering a tomogram of N depth profiles covering a lateral range of ∆x, the two center spatial carrier frequencies can be calculated by The lateral scanning range ∆x is calculated by where f is the focal length of the sample objective.
The positive spatial frequency range can be calculated by For high-quality image reconstruction, another important parameter is the characteristic bandwidth (BW) of the obtained spatial frequency contents, which is determined by the inverse speckle size.The speckle size can be modeled by the transverse resolution, and BW can be expressed as [9] where d is the collimated beam diameter entering the sample objective.
In theory, the proposed method could quadruple the effective imaging depth range compared with a conventional SDOCT system.However, practically some factors can affect the improvement of the range depth, which includes the crosstalk between the two spatial frequency contents of the two depths signal, the crosstalk between the spatial frequency content and its complex conjugate content, and the crosstalk between the spatial frequency contents and the interference term between the two reference arms' light fields.Crosstalk between the two spatial frequency contents corresponding to the two depths can occur when the tails of the spatial frequency contents superpose with each other as schematically represented in Figure 4a.In view of the Bedrosian theorem, spatial frequency content should be well separated from the introduced modulation frequency and the other content for successful reconstruction, as shown in Figure 4b.The crosstalk due to the reference arms' interference can be eliminated via background subtraction.According to Equations ( 3)-( 7), the parameters including the angular velocity ω, the complementary metal oxide semiconductor (CMOS) exposure time T, and the offset distance s i are all responsible to the image reconstruction quality.Thus, to avoid crosstalk and fully utilize the proposed method for extended range imaging, the center spatial carrier frequencies and the system parameters must be chosen carefully.
where d is the collimated beam diameter entering the sample objective.
In theory, the proposed method could quadruple the effective imaging depth range compared with a conventional SDOCT system.However, practically some factors can affect the improvement of the range depth, which includes the crosstalk between the two spatial frequency contents of the two depths signal, the crosstalk between the spatial frequency content and its complex conjugate content, and the crosstalk between the spatial frequency contents and the interference term between the two reference arms' light fields.Crosstalk between the two spatial frequency contents corresponding to the two depths can occur when the tails of the spatial frequency contents superpose with each other as schematically represented in Figure 4a.In view of the Bedrosian theorem, spatial frequency content should be well separated from the introduced modulation frequency and the other content for successful reconstruction, as shown in Figure 4b.The crosstalk due to the reference arms' interference can be eliminated via background subtraction.According to Equations ( 3)-( 7), the parameters including the angular velocity ω, the complementary metal oxide semiconductor (CMOS) exposure time T, and the offset distance si are all responsible to the image reconstruction quality.Thus, to avoid crosstalk and fully utilize the proposed method for extended range imaging, the center spatial carrier frequencies and the system parameters must be chosen carefully.To quantify the crosstalk effect, we introduce the SCR defined as a function of the frequencyspacing/bandwidth ratio ( BW   ) for two reference arms, where Furthermore, to avoid aliasing terms, the bandwidth of the spatial frequency content should be narrow relative to the full spatial frequency range; in this case, the small phase fluctuations will not cause a reduction of the SCR.This condition indicates that [9] The criteria shown in Equations ( 8) and ( 9) will be used to guide the choice of the values of the center spatial carrier frequencies and the other parameters in the experimental SDOCT system.

Experiments and Results
For experimental verification of the method, we constructed a fiber-based spatial frequency To quantify the crosstalk effect, we introduce the SCR defined as a function of the frequency-spacing/bandwidth ratio (∆ν/BW) for two reference arms, where ∆ν = ν c2 − ν c1 as shown in Figure 4.According to the theoretical value, the SCR can be higher than 100 dB when the frequency-spacing/bandwidth ratio τ satisfies the condition [12] Furthermore, to avoid aliasing terms, the bandwidth of the spatial frequency content should be narrow relative to the full spatial frequency range; in this case, the small phase fluctuations will not cause a reduction of the SCR.This condition indicates that [9] The criteria shown in Equations ( 8) and ( 9) will be used to guide the choice of the values of the center spatial carrier frequencies and the other parameters in the experimental SDOCT system.

Experiments and Results
For experimental verification of the method, we constructed a fiber-based spatial frequency domain multiplexing SDOCT system as depicted in Figure 1.The imaging system employs a super luminescent diode light source (SLD-371-HP2, Superlum Inc., Carrigtwohill, Ireland) with a center wavelength at 837 nm and a full width at half maximum (FWHM) of 54 nm.The beam diameter d in the sample arm is 1.8 mm, which will produce a 44 µm 1/e 2 diameter focal spot on the sample with a confocal parameter of 3.7 mm.A 50/50 fiber coupler splits the light beam into the sample arm and the reference arms.In the sample arm, the light beam is deflected by a 2D galvo scanner and focused via a 75 mm achromatic lens L0 to the sample.Within the reference arm, a wide bandwidth non-polarizing beam splitter (NPBS, Daheng Optics Inc., Beijing, China) is used to divide the reference optical power into the two separate reference arms.The fiber collimator on a translation stage in the reference arm is used to adjust the overall delay of the two channels.The single-pass OPD between Reference Arms I and II is set to 3 mm.The GS1 and GS2 are the X and Y mirror of a pair of large beam 2D GS system (GVS012, Thorlabs Inc., Newton, NJ, USA), respectively.The focus lengths of achromatic lens L1and L2 are both 30 mm.The plane mirror M1 and M2 are the wide bandwidth dielectric coated mirrors covering the wavelength range of 800 to 900 nm.The spectral interference signal is detected by a custom-built spectrometer, which consists of a collimator (f = 60 mm, OZ optics Inc., Ottawa, ON, Canada), a transmissive diffraction grating (1800 lines/mm), a camera lens (f =105 mm, Nikon Inc., Tokyo, Japan), and a line-scan CMOS camera (sp2048-70km, Basler AG Inc., Ahrensburg, German).The spectral resolution of the spectrometer is 0.065 nm, allowing for a theoretical effective imaging depth range of 3.5 mm.
In view of avoiding aliasing terms and according to Equation ( 9), the number of A-scans per B-scan is set to 500, and the lateral scanning range is set to 1.3 mm.The highest positive spatial frequency calculated by Equation ( 6) is 192 mm −1 .According to the values of the sample beam diameter d, the center wavelength, the focus length of the sample objective and Equation (7), the theoretical BW of the spatial frequency content is calculated to be 22.5 mm −1 .In order to suppress the crosstalk of the two spatial frequency contents of the two depths signal, the difference of the two center spatial carrier frequencies must be greater than 67.5 mm −1 according to Equation (8).Considering that the optimal conjugate suppression ratio occurs when δϕ is π/2 [9], each of the two center spatial carrier frequencies needs to be close to the central region of the positive spatial frequency range.Thus, the two center spatial carrier frequencies ν c1 and ν c2 are set to 56.5 and 123.5 mm −1 , respectively.According to Equation ( 4), the phase shift δϕ 1 and δϕ 2 between successive axial scans are calculated to be 0.3π and 0.7π.The angular velocity of the galvo scanners is set to 0.93 rad/s.The exposure time of the CMOS camera is set to 18 µs.According to Equation (3) and the above parameters, the pivot-offset distances on the GS1 and GS2 mirrors are calculated to be 1.8 and 3.9 mm, respectively.The pivot-offset distance is adjusted via the diagonal translation of the scanners to keep the light beam at the center of the focusing optics.
To experimentally characterize the extension of the imaging depth range of the spatial frequency domain multiplexing SDOCT system, the point spread function is measured by using a plane mirror at various depths in the sample arm.Figure 5 shows the experimentally measured normalized sensitivity as a function of depth for the two reference arms.The measured value of the absolute sensitivity is 59.7 dB.From Figure 5, the separation between the two sensitivity peaks can be used to calibrate the actual relative OPD between the two reference arms, and it can be shown that the actual relative OPD is about 3 mm.The benefit from the use of the spatial frequency domain multiplexing SDOCT system is also demonstrated, and the effective imaging depth range for the proposed system is increased from about 3.4 to 6.4 mm, compared with a conventional full-range SDOCT without depth range enhancement.The depth enhancement factor is 1.9.In this demonstration, the two sensitivity curves are crossed at the −5.5 dB from the sensitivity peaks.Therefore, the effective imaging depth range would be doubled if the relative optical path delay was precisely tuned to make the two sensitivity curves cross at the −6 dB point from the sensitivity peaks.To demonstrate the ranging depth enhancement, we compare image results of the OCT images acquired by the proposed two-reference-arm spatial frequency domain multiplexing approach and the conventional single-reference-arm approach.Firstly, an artificial phantom is imaged by the system.Figure 6a shows the photograph of the phantom taken with a digital camera.The phantom is made by sticking two blocks of sellotape to generate a step 3 mm high.Figure 6b-d are the complex conjugate ambiguity free OCT images of the phantom reconstructed by the data processing algorithm described previously.Figure 6b shows the image acquired with Reference Arm I open while blocking Reference Arm II.As we can see, the upper part of the step is clearly imaged, while the image of the lower part is hard to be discerned.The OPD between the upper part and Reference Arm I resides in the effective imaging range, while that of the lower part and Reference Arm I is larger than 6 mm, out of the effective imaging range.Figure 6c shows the OCT image acquired with Reference Arm II open while blocking Reference Arm I. Likewise, the lower part of the step can be clearly imaged, while the upper part is not clear.Figure 6d shows the reconstructed OCT image of the phantom acquired with both Reference Arms I and II open.The whole phantom including the upper and lower part of the step can be clearly seen.The imaging depth-enhanced OCT image of the artificial phantom demonstrates the feasibility of the proposed method for extending the OCT ranging depth.To demonstrate the ranging depth enhancement, we compare image results of the OCT images acquired by the proposed two-reference-arm spatial frequency domain multiplexing approach and the conventional single-reference-arm approach.Firstly, an artificial phantom is imaged by the system.Figure 6a shows the photograph of the phantom taken with a digital camera.The phantom is made by sticking two blocks of sellotape to generate a step 3 mm high.Figure 6b-d are the complex conjugate ambiguity free OCT images of the phantom reconstructed by the data processing algorithm described previously.Figure 6b shows the image acquired with Reference Arm I open while blocking Reference Arm II.As we can see, the upper part of the step is clearly imaged, while the image of the lower part is hard to be discerned.The OPD between the upper part and Reference Arm I resides in the effective imaging range, while that of the lower part and Reference Arm I is larger than 6 mm, out of the effective imaging range.Figure 6c shows the OCT image acquired with Reference Arm II open while blocking Reference Arm I. Likewise, the lower part of the step can be clearly imaged, while the upper part is not clear.Figure 6d shows the reconstructed OCT image of the phantom acquired with both Reference Arms I and II open.The whole phantom including the upper and lower part of the step can be clearly seen.The imaging depth-enhanced OCT image of the artificial phantom demonstrates the feasibility of the proposed method for extending the OCT ranging depth.
The biological tissue of the swine adipose is also in vitro, imaged and compared with the results acquired by the proposed two-reference-arm spatial frequency domain multiplexing approach and the conventional single-reference-arm approach.Figure 7a,b show the images recorded using the two Reference Arms I and II separately.In each case, the area of the tissue around the DC signal component has maximum sensitivity.Figure 7c shows the OCT image of the tissue acquired with both of Reference Arms I and II open.As we can see, the upper and lower part of the biological sample is clearly appeared at the same time.Figure 7d plots an A-scan profile from the same sample position indicated by the red arrows in Figure 7a,b.It can be seen that the signal is about 13 dB stronger than that of the signal 3 mm away due to the sensitivity roll-off.The extended range OCT image of the biological tissue obtained with the system demonstrates the feasibility of the proposed method for extending the OCT ranging depth.Because the galvanometer scanners in the system are used for phase modulation of the lateral interference spectral signal, the phase noise does not affect the quality of the results.In future, the proposed method for extending the OCT ranging depth can be combined with the depth of the focus extension method to improve the quality of the extended range OCT image.
the effective imaging range, while that of the lower part and Reference Arm I is larger than 6 mm, out of the effective imaging range.Figure 6c shows the OCT image acquired with Reference Arm II open while blocking Reference Arm I. Likewise, the lower part of the step can be clearly imaged, while the upper part is not clear.Figure 6d shows the reconstructed OCT image of the phantom acquired with both Reference Arms I and II open.The whole phantom including the upper and lower part of the step can be clearly seen.The imaging depth-enhanced OCT image of the artificial phantom demonstrates the feasibility of the proposed method for extending the OCT ranging depth.The biological tissue of the swine adipose is also in vitro, imaged and compared with the results acquired by the proposed two-reference-arm spatial frequency domain multiplexing approach and the conventional single-reference-arm approach.Figure 7a,b show the images recorded using the two Reference Arms I and II separately.In each case, the area of the tissue around the DC signal component has maximum sensitivity.Figure 7c shows the OCT image of the tissue acquired with both of Reference Arms I and II open.As we can see, the upper and lower part of the biological sample is clearly appeared at the same time.Figure 7d plots an A-scan profile from the same sample position indicated by the red arrows in Figure 7a,b.It can be seen that the signal is about 13 dB stronger than that of the signal 3 mm away due to the sensitivity roll-off.The extended range OCT image of the biological tissue obtained with the system demonstrates the feasibility of the proposed method for extending the OCT ranging depth.Because the galvanometer scanners in the system are used for phase modulation of the lateral interference spectral signal, the phase noise does not affect the quality of the results.In future, the proposed method for extending the OCT ranging depth can be combined with the depth of the focus extension method to improve the quality of the extended range OCT image.

Conclusions
We have proposed and demonstrated a spatial frequency domain multiplexing method to extend the effective imaging depth range for SDOCT system without any expensive acousto-optic or

Conclusions
We have proposed and demonstrated a spatial frequency domain multiplexing method to extend the effective imaging depth range for SDOCT system without any expensive acousto-optic or electro-optic modulation device.The spatial frequency domain multiplexing is realized by using two pivot-offset galvo scanners in the two reference arms.The images of shallow and deep regions of a sample can be multiplexed and acquired simultaneously with no image acquisition rate penalty.The proposed two-reference-arm spatial frequency domain multiplexing SDOCT system provides an approximately 1.9-fold increase in the effective imaging range compared with that of a conventional single-reference-arm full-range SDOCT system.We present key system parameter design criteria for ensuring a high signal/crosstalk ratio imaging.The images of an artificial phantom and the swine adipose tissue reconstructed with the proposed method demonstrate the expected increase in the effective imaging depth range.This simple and low-cost method is also applicable to the SSOCT system, and desirable in the in vivo applications requiring a relatively long imaging range.

2 λk
  , z is the depth position coordinate of the scatters within the sample, n is the average refractive index of the sample,   , R x z is the normalized intensity backscattered from the scatter at z, S(k) is the spectral density of the light source, represent the linear modulated phase induced by GS1 and GS2, and ( , ) A x z  and ( , ) B x z  are the phases that relate to the optical heterogeneity of the sample in the focus coherence volume.If the linear modulation phase about the lateral position, the


are the pair of Hermit conjugate term due to the Fourier transform of the real valued interference signal and contain the same information on the OCT signal in the spatial frequency domain.The subscripts A and B in Equation (

Figure 2 .
Figure 2. Data processing flowchart of the spatial frequency domain multiplexing SDOCT system.FT: Fourier Transform

Figure 2 .
Figure 2. Data processing flowchart of the spatial frequency domain multiplexing SDOCT system.FT: Fourier Transform.

Figure 3
Figure3is the sketch of the two reference arms with the pivot-offset galvo scanners.The color lines represent the light beams when the galvo scanners are at different scan angles.The pivot-offset distances are indicated as s 1 and s 2 .The phase shift between successive axial scans (A-scan) induced by the galvo scanner can be expressed as[9]

Figure 3
Figure3is the sketch of the two reference arms with the pivot-offset galvo scanners.The color lines represent the light beams when the galvo scanners are at different scan angles.The pivot-offset distances are indicated as s1 and s2.The phase shift between successive axial scans (A-scan) induced by the galvo scanner can be expressed as[9]

Figure 3 .
Figure 3. Sketch of the two reference arms in the spatial frequency domain multiplexing SDOCT system.s1 and s2: the pivot-offset distances; L1 and L2: the focal lens in the reference arms; M1 and M2: the plane mirrors; f: focal length of the focal lens.

Figure 3 .
Figure 3. Sketch of the two reference arms in the spatial frequency domain multiplexing SDOCT system.s 1 and s 2 : the pivot-offset distances; L1 and L2: the focal lens in the reference arms; M1 and M2: the plane mirrors; f: focal length of the focal lens.

Figure 4 .
Figure 4. Schematic of the choice of the center spatial carrier frequencies.Crosstalk between the spatial spectra occurs in (a), and the two spatial spectra can be well separated through parameters design as shown in (b).

4 .
According to the theoretical value, the SCR can be higher than 100 dB when the frequencyspacing/bandwidth ratio τ satisfies the condition[12]

Figure 4 .
Figure 4. Schematic of the choice of the center spatial carrier frequencies.Crosstalk between the spatial spectra occurs in (a), and the two spatial spectra can be well separated through parameters design as shown in (b).

Figure 5 .
Figure 5.The measured sensitivity curves corresponding to the Reference Arm I [1] and Reference Arm II [2].

Figure 5 .
Figure 5.The measured sensitivity curves corresponding to the Reference Arm I and Reference Arm II.

Figure 6 .
Figure 6.(a) The photograph of the artificial phantom of a step comprised of two blocks of sellotape.The red line is the lateral scan position on the phantom; (b,c) the cross-sectional optical coherence tomography (OCT) image of the phantom acquired with the conventional single-reference-arm approach; (d) the extended range OCT image acquired with the proposed two-reference-arm spatial frequency domain multiplexing approach.The blurry horizontal white line is the residual direct current (DC) term.

Figure 6 .
Figure 6.(a) The photograph of the artificial phantom of a step comprised of two blocks of sellotape.The red line is the lateral scan position on the phantom;(b, c) the cross-sectional optical coherence tomography (OCT) image of the phantom acquired with the conventional single-reference-arm approach; (d) the extended range OCT image acquired with the proposed two-reference-arm spatial frequency domain multiplexing approach.The blurry horizontal white line is the residual direct current (DC) term.

Figure 7 .
Figure 7.Cross-sectionalOCT image of the swine adipose tissue, (a, b) acquired with the conventional single-reference-arm approach.(c)The extended range OCT image is obtained with the tworeference-arm spatial frequency domain multiplexing approach.The blurry horizontal white line is the residual DC terms.(d) A-scans at the position marked by the red arrows in (a, b) demonstrate the advantage of the ranging depth extension.

Figure 7 .
Figure 7. Cross-sectionalOCT image of the swine adipose tissue, (a,b) acquired with the conventional single-reference-arm approach.(c) The extended range OCT image is obtained with the two-reference-arm spatial frequency domain multiplexing approach.The blurry horizontal white line is the residual DC terms.(d) A-scans at the position marked by the red arrows in (a,b) demonstrate the advantage of the ranging depth extension.