Spectral domain detection in low-coherence spectroscopy.

Low-coherence spectroscopy (LCS) offers the valuable possibility to measure quantitative and wavelength resolved optical property spectra within a tissue volume of choice that is controllable both in size and in depth. Until now, only time domain detection was investigated for LCS (tdLCS), but spectral domain detection offers a theoretical speed/sensitivity advantage over tdLCS. In this article, we introduce a method for spectral domain detection in LCS (sdLCS), with optimal sensitivity as a function of measurement depth. We validate our method computationally in a simulation and experimentally on a phantom with known optical properties. The attenuation, absorption and scattering coefficient spectra from the phantom that were measured by sdLCS agree well with the expected optical properties and the measured optical properties by tdLCS.

1. Introduction

Low-coherence spectroscopy (LCS) is a method for simultaneous measurements of local attenuation, absorption, scattering and backscattering coefficient spectra in turbid media within the visible wavelength range (480 to 700 nm) [1-3]. LCS distinguishes itself from other quantitative spectroscopic techniques by controlling the photon path length of the detected light, which offers the unique possibility to control the measurement volume both in size and in depth—e.g. for the in vivo measurement of hemoglobin concentrations within distinct regions of the human dermal microcirculation [3]. Our previous studies on LCS involved a time-domain detection setup (tdLCS). However, from optical coherence tomography (OCT) studies, we know that spectral domain detection allows for more sensitive measurements, or higher acquisition speed at unaltered sensitivity [4,5]. Since clinical applications of LCS require both fast and sensitive measurements, we investigate the possibility for spectral domain detection in LCS (sdLCS), focusing primarily on validating a new method of analysis that was developed for this purpose.

Signal acquisition in sdLCS is comparable to that in spectrometer-based, spectral domain OCT (sdOCT), with the major difference that in sdLCS optical property spectra are derived from the spectral content of the signal, and in sdOCT an image is reconstructed from the backscattered intensity. To obtain optical property spectra within confined tissue volumes, sdLCS requires both high spatial and high spectral resolution. The spatial resolution for the acquisition of one backscattered spectrum in sdLCS and sdOCT is limited by the axial measurement range of the spectrograph, which depends on the spectral resolution of the spectrograph. Since the standard methods for signal acquisition in sdOCT commonly result in only one backscattered spectrum within the entire axial measurement range of the spectrograph (typically 1 to 2 mm) [4,5], the spatial resolution of these methods is not sufficient for sdLCS (typically 22 µm) [1-3]. Spectroscopic sdOCT methods solve this problem typically by windowed Fourier transforms on the measured spectrum, which result in multiple depth resolved spectra [6-8]. Due to the sensitivity decrease with depth inherent to these methods and the absence of focus tracking, these methods may not achieve optimal sensitivity in depth. In addition, precise refractive index-dependent corrections for the unwanted signal attenuation due to these effects are needed, in order to achieve quantitative measurements of the tissue attenuation. Therefore, we developed a method for sdLCS that applies focus tracking and measures with maximum spectrograph sensitivity (around zero delay) at every depth.

Our approach involves acquisition of local back scattered spectra within an axial measurement range of only 22 µm. Translation of this axial measurement range through the sample results in depth-resolved backscattered spectra acquired around zero delay with high spatial and sufficiently high spectral resolution (6 nm) while focus tracking is applied. The main challenge in our approach to sdLCS is to isolate the sample’s backscattered spectrum from the detected spectrum, in order to retrieve the spectral absorption and backscattering features from the sample only. For this purpose, the unwanted DC and modulation components in the detected spectrum can be efficiently removed by introducing a phase modulation in the reference arm [9,10], for which we use an oscillating reference mirror.

Since the primary aim of this study is to validate whether we can use this new approach for the determination of local absorption and scattering coefficient spectra, we will provide proof of principle with an on-the-shelf spectrograph (USB4000, Ocean Optics, USA). Only part of the feasible speed advantage will be demonstrated, since the quantum efficiency and acquisition speed of this spectrograph are not optimized for sdLCS. We will validate our algorithm with a simulation of an sdLCS absorption measurement, while simulating the USB4000 as the detecting spectrograph. Subsequently, we experimentally demonstrate the ability of sdLCS to measure the attenuation, absorption and scattering coefficient spectra from a polystyrene-dye phantom, and we compare our results to those of a tdLCS measurement.

2. Theory

2.1. Background

For sdLCS, we use a setup that is based on a Michelson interferometer and a broad band light source (Fig. 1 ). A spectrograph detects the spectrum iD, which (for a single reflector) is modulated by the optical path length difference between the sample arm and the reference arm ΔOPL = 2(xS – xR), with sample arm length xS and reference arm length xR:

Fig. 1

Schematic illustration of setup and parameters in sdLCS/tdLCS. λ0: center wavelength, λBW: spectral bandwidth, xR: reference arm length, xS: sample arm length, ΔOPL: path length difference, ΔR: reference mirror scanning range, vR: reference mirror scanning velocity, fR: reference mirror scanning frequency, ℓ: geometrical path length in sample, ΔℓR: scanning window in sample induced by ΔR, n: refractive index, µt: attenuation coefficient, iD: detector current, τ: integration time, Np: number of pixels, δλ: pixel width, Δλ: wavelength resolution.

where k is the wavenumber (k = 2π/λ) and IS and IR are the signal intensities from the sample and reference arm, respectively.

Since we use sdLCS for doing spectroscopy, we are interested in the spectrum IS(k). The sample arm contribution IS(k) to iD(k) is related to the sample’s backscattered power spectrum S, which can be obtained from iD(k) using the procedure described in Section 2.2. Note that all bold faced characters in this article represent wavelength (λ) dependent parameters. By adjusting xR, we can control the geometrical round-trip path length ℓ in the sample at which we obtain S and by step-wise alteration of ℓ, we obtain a data set S(ℓ) [1-3]. Subsequently, we use a single exponential decay model to obtain the attenuation coefficient spectrum µ of the sample:

assuming local sample homogeneity for all path lengths in the data set S(ℓ). The model is fitted to S(ℓ) with free running fit parameters (α and µ), similar to our analysis for tdLCS [1-3]. The contributions of the scattering coefficient µ and the absorption coefficient µ can be obtained by fitting their individual contributions to the measured µ, as described in Ref. [3].

The axial measurement range of the spectrograph is limited by its maximal measurable path length difference and sensitivity roll-off in depth [4]. Therefore, the spectra S(ℓ) are acquired within a path length window Δℓ, which determines the spatial resolution for the µ-determination (see Section 2.3). Since a small axial measurement range is related to low spectral resolution [11], a tradeoff exists between spectral and spatial resolution in sdLCS. However, we will show in this article that by selecting the right spectrograph, both sufficiently high spectral (6 nm) and spatial resolution (axial: 22 µm) can be achieved.

2.2. Removal of the DC and modulation terms

The sample’s backscattered spectrum S(ℓ) in Eq. (2) cannot be obtained directly from iD(k), because iD(k) contains unwanted DC (i.e. non-modulated components) and modulation components (Eq. (1). Since the spectral features in S—originating from the sample’s characteristic absorption and backscattering—may exhibit oscillation frequencies similar to the modulation on iD(k), S cannot be directly obtained from iD(k). Similar to the method described for OCT in Ref. [9], we will isolate S from iD(k) by phase modulation and frequency domain filtering of the signal.

We start by modulating iD(k) using an oscillating mirror in the reference arm of our LCS system. This introduces a Doppler shift fD for every wavelength: fD = 2vR/λ = 2ΔR∙fR∙k/π, with vR the velocity, fR the oscillation frequency and ΔR the scanning range in one direction of the reference mirror. Acquisition of N spectra iD(k) over a time interval N(τ + τd) with integration time τ and dead time τd results in a data set iD(k,t), as illustrated in Fig. 2(a) . Since τd (12 µs) << τ (6 ms) for the measurements described in this manuscript, we apply N(τ + τd) ≈Nτ in the remaining part of our analysis. The frequency of the modulation on iD(k,t) changes as a function of time, because the movement of the reference mirror induces a change in ΔOPL (vertical direction in Fig. 2(a)). Figure 2(b) illustrates the time-modulation on iD(k,t) by fD for one wavelength (horizontal direction in Fig. 2(a)). Since the sensitivity of the measurement is largest around ΔOPL = 0 or ‘zero delay’ (see Eq. (3), the range of ΔOPL is chosen such that it crosses zero delay at ½Nτ (Fig. 2(c)). Adequate sampling of the time-modulation on iD(k,t) is determined by the Nyquist criterion, i.e. the sampling frequency fs = 1/τ needs to be larger than 2fD,max (the largest value of fD within the investigated spectral range). For the integration time, this results in τ < π/(2vRkmax) = λmin/(4vR), with kmax the largest wave number and λmin the shortest wavelength within the investigated spectral range.

Fig. 2

Signal acquisition and processing in sdLCS. See Section 2.2 for details.

Fourier transformation on iD(k,t) with respect to time provides the frequency content of the signal iD(k,f), which contains the modulation frequency fD, its mirror image and a DC component (Fig. 2(d)). By filtering out the part of iD(k,f) containing only the positive frequencies + fD, the unwanted DC component and the negative frequencies –fD are lost. Inverse Fourier transformation on the filtered iD(k,f) with respect to f provides a smooth, non-modulated iD(k,t) that is confined in both amplitude and phase (Fig. 2(e)). Averaging of iD(k,t) over time results in the final spectrum S that we can use for the determination of µ (Fig. 2(f)). A computational validation of this method will be shown by means of a simulation in Section 3.3.

2.3. The path length window Δℓ in sdLCS

In our approach for sdLCS, the path length window Δℓ over which S(λ) is acquired results from a convolution of the reference mirror scanning window ΔℓR in the sample (ΔℓR = 2ΔR/n, with n the refractive index of the sample), with the path length window that is probed by the spectrograph ΔℓS. In practice, this path length window can be approached with Δℓ = ΔℓR + ΔℓS (inset of Fig. 3 ). The spectrograph probing window ΔℓS depends on the maximal measurable path length difference by the spectrograph ΔOPLmax (or imaging depth dmax = ½ΔOPLmax) and the sensitivity roll-off with ΔOPL of the spectrograph. Due to the Nyquist criterion, ΔOPLmax = 2π/δk with ‘pixel resolution’ δk = (kmax-kmin)/Np and Np the number of pixels of the spectrograph. The sensitivity roll-off of the spectrograph can be calculated using [11]

Fig. 3

Theoretical and measured sensitivity roll-off of the spectrograph. Inset: illustration of the path length window Δℓ in sdLCS (Section 2.3).

in which Δk is the spectral resolution of the spectrograph, assuming a Gaussian beam profile.

For our sdLCS experiments, we use a spectrograph with a sensitivity roll-off that is two orders of magnitude smaller than ΔOPLmax (Section 3.1). Therefore, we can approach ΔℓS with the full width half max (FWHM) of the roll-off function: ΔℓS = FWHMroll-off/n. Hence, the path length window in sdLCS that determines the spatial resolution of the measurements is determined by: Δℓ = ΔℓR + ΔℓS = (2ΔR + FWHMroll-off)/n.

2.4. Sensitivity/speed analysis

The SNR of any shot-noise limited LCS system is given by

Here, ε is the detector quantum efficiency; hv is the photon energy and fBW is the electronic detection bandwidth of the system [5]. For a spectral domain (sdLCS) system, fBW = 1/τ [4]. For a time domain (tdLCS) system, fBW is determined by the imposed Doppler frequency f = 2vR/λ = 2Lscan,TD∙fscan,TD /λ of the scanning reference mirror [1], with vR its velocity, Lscan,TD = 2ΔR its scanning range in air and fscan,TD = fR its scanning frequency. A commonly used comparison results in an improved SNR with a factor Np (number of spectrograph pixels) for sdLCS compared to tdLCS [4,5,12], provided that ε, Lscan and the fscan (i.e. acquisition speed) are equal for both methods.

For a more generalized comparison between sdLCS and tdLCS, we take into account the differences in quantum efficiency, scan speed and scan range of the two methods (λ0 and λBW are equal since we employ the same light source):

with ΔLmax the maximal measurable path length difference by the spectrograph after removal of the DC and modulation terms.

3. Materials and methods

In order to validate our data acquisition and analysis approach for sdLCS, we simulated an absorption measurement with the data acquisition parameters of our experimental setup. In addition, we measured attenuation spectra µ from 480 to 700 nm on a phantom with known absorption µ and scattering µ coefficient spectra and we separated the individual contributions of µ and µ from the measured µ. To compare sdLCS to tdLCS, we performed a signal to noise analysis for both systems and repeated the phantom measurement with the tdLCS system described in detail in Ref. [1].

3.1. System and acquisition

The sdLCS system described in this article is identical to the tdLCS system with a low-pass filtered supercontinuum light source (SC430-4, Fianium Ltd., UK) described in Ref. [1-3], except for the detection end of the system, which consists of a spectrograph rather than a photodiode/lock-in amplifier combination. Hence, the multimode detection fiber is connected to the detecting spectrograph (USB4000, Ocean Optics, USA), with λmin = 345 nm, λmax = 1042 nm and Np = 3648.

For both the sdLCS and tdLCS measurements, we controlled ℓ (ℓ = 0 – 2000 µm, with ℓ = 0 the sample surface) by translating the reference mirror in steps of 27 µm. By translating the sample in the axial direction, focus tracking of the spot size (r = 4.5 μm) in the sample was achieved. At every ℓ, back scattered power spectra S(ℓ) were obtained over a path length window in the sample of Δℓ = 43 μm for sdLCS (see Section 3.2) and Δℓ = 44 μm for tdLCS [1]. For both sdLCS and tdLCS, S(ℓ) was temporally averaged (N = 250 spectra per ℓ) and corrected for the background (i.e. the LCS spectrum obtained from the non-scattering glass of the inner cuvette wall of the sample). The integration time was set at the minimal stable integration time of this spectrograph, τ = 6 ms. To meet the sampling frequency requirement (Section 2.2), the velocity of the reference mirror was set at vR = 12 µm/s (frequency fR = 0.3 Hz, scanning range ΔR = 20 µm). The acquisition settings for the tdLCS measurements were identical to the settings in Ref. [1-3]. For both sdLCS and tdLCS, the most important acquisition parameters have been summarized in Table 1 . Note that the acquisition speed for one temporally averaged spectrum S(ℓ) is 3.6 times faster for sdLCS.

Table 1

Acquisition settings for sdLCS and tdLCS

Acquisition parameter	sdLCS	tdLCS
Δℓ	43 μm	44 μm
Δλ	6 nm (all λ)	4 nm @ λ = 480 nm
		9 nm @ λ = 700 nm
vR	0.006 mm/s	1.84 mm/s
fR	0.3 Hz	23 Hz
ΔR	20 μm	40 μm
N = # averages per ℓ	250	250
Acquisition time of S(ℓ)	Nτ = 1.5 s	N/(2fR) = 5.4 s

As described in Section 2, fitting the single exponential decay model S(ℓ) = α∙exp(-μ∙ℓ) (free running fit parameters α and μ) to the sdLCS and tdLCS acquired S(ℓ) vs. ℓ, results in a μ spectrum for both detection methods [1-3]. Uncertainties in α and μ are estimated by their 95% confidence intervals (c.i.). Since the spectral resolution in sdLCS (Δλ = 6 nm) is higher than the pixel width of 0.2 nm, S(ℓ) was binned into wavelength regions of 6 nm for the sdLCS measurements. The values of µ and µ were obtained by fitting their individual contributions to the measured µ, as described in Ref. [3].

3.2. System characterization

Using the specifications of the spectrograph and the theory described in Section 2.3, we calculated that ΔOPLmax = 1880 µm. After measuring the spectral resolution of Δλ = 6 nm with a 543 nm and 633 nm HeNe laser line, we used Eq. (3) to calculate that the FWHMroll-off is 18 µm (‘theoretical’ roll-off in Fig. 3). This FWHMroll-off corresponds well with the experimentally measured roll-off function, measured from the reflection of a glass slide in the sample arm as a function of ΔOPL between the sample arm and the reference arm (Fig. 3). Since the FWHMroll-off of the spectrograph is two orders of magnitude smaller than ΔOPLmax, the path length window that is probed by the spectrograph can be approximated with ΔℓS = FWHMroll-off/n = 13 µm. Consequently, the total path length window that is probed in our sdLCS measurements is Δℓ = ΔℓR + ΔℓS = (2ΔR + FWHMroll-off)/n = 43 µm, with a fixed ΔR of 20 µm.

Based on Eq. (5), we can calculate the expected difference in SNR between both systems. The sdLCS spectrum covers approximately Np = 1150 pixels on the USB4000. Based on a maximum number of 60 photons/count and a pixel well depth of 1.105 electrons given by the manufacturer, the quantum efficiency (QE = #electrons/photon) of the 16-bits USB4000 is 2.5% at 600 nm. Using the photodiode response of 0.43 A/W (2001, New Focus, USA), we can calculate that our tdLCS detector has a QE of 89% at 600 nm. With τ = 6 ms for our sdLCS system, the sensitivity advantage of the sdLCS measurements with the current equipment reduces to 0.38. Hence, when using the USB4000 as a detector for sdLCS with the settings in Table 1, we can expect a decrease in SNR of 4.2 dB at 600 nm with respect to tdLCS.

In order to compare the sensitivity of the current sdLCS system to the tdLCS system, we measured the SNR of the two detection methods, using a mirror behind an attenuating (2.0 optical density) filter in the sample arm. The optical power reference arm was optimized for (near) shot noise limited detection and the acquisition settings in Table 1 were applied. For both sdLCS and tdLCS, the spectrally resolved SNR was determined from the measured spectrum S at ΔOPL = 0 using [5]: SNR = 10·log(S/var[noise]). Here, var[noise] denotes the variance of the noise spectrum, measured with a ΔOPL large enough to observe no interference signal within the path length window Δℓ. For both detection methods, the measured SNR is around 75 dB per wavelength bin between 480 and 550 nm. After 550 nm, the SNR of tdLCS increases to a value of 85 dB at 700 nm, which can be explained by the larger quantum efficiency of the photodiode in tdLCS for longer wavelengths. The SNR of sdLCS remains spectrally flat over the full wavelength region of 480 to 700 nm, resulting in an SNR of ~3 dB lower compared to tdLCS at 600 nm, which is close to the predicted value of 4.2 dB.

3.3. Simulation

To validate the analysis algorithm described in Section 2.2, we simulated an absorption measurement on a sample containing oxyhemoglobin [13]. With this simulation, we tried to mimic the reality of the acquisition as closely as possible; therefore we incorporated summation over the integration time for every spectrum while simulating reference mirror movement, a discrete wavelength axis and a limited spectral resolution. As a consequence, dynamic effects on such as fringe washout are taken into account. The contributions of scattering and noise were neglected in the simulation.

The detecting spectrograph was simulated with the properties of the spectrograph used for our measurements, as well as the values for the integration time and reference mirror velocity. The source spectrum S was simulated as a Gaussian with center wavelength λ0 = 550 nm and bandwidth λFWHM = 140 nm. To account for the spectral resolution (Δλ) of the spectrograph, we convolved the spectrum with a Gaussian with a full-width at half maximum of 6 nm. The simulated sample arm power fraction was ηS = 0.1, resulting in a reference arm power fraction of ηR = 0.9. Within each integration time, 16 time-separated spectra at different reference mirror positions were added to simulate an acquired spectrum. A total number of N = 125 acquired spectra was simulated over a time interval of Nτ = 0.75 s. The initial ΔOPL was set at −4.5 µm to achieve zero delay at ½Nτ.

Since scattering was neglected in the simulation, the attenuation coefficient µ of the sample can be simplified to the absorption coefficient µ. The spectrum S(ℓ) was calculated at ℓ1 = 2 mm and ℓ2 = 4 mm inside the sample, corresponding to depths d1 = 1 mm and at d2 = 2 mm, respectively. We obtained the absorption coefficient of the sample using Beer’s law (Eq. (2): µ = –[½(ℓ2 – ℓ1)]−1∙ln[S(ℓ2)/ S(ℓ2)].

3.4. Phantom

We prepared a phantom, consisting of 0.096 vol% NIST-certified polystyrene spheres (⌀ 602 ± 6 nm, Thermo Scientific, USA) and 10% magenta dye (Ecoline #337, Royal Talens, The Netherlands). The µ of the polystyrene spheres was calculated using Mie theory and integrated over the size distribution of the spheres (2*SD), as described in Ref. [2]. and the µ of the dye was determined in a separate transmission measurement on the dye only, as described in Ref. [1]. The measured tdLCS spectra are broadened by Doppler frequency shifts induced by the Brownian motion of the polystyrene spheres. Therefore, both reference spectra for µ and µ were convolved with a Lorentzian with a line width of 9 nm for adequate comparison to the LCS spectra [1-3]. For the sdLCS measurements, Brownian motion does not influence the spectrum S, because the Doppler shifted frequencies remain within the bounds of the frequency domain filter (Fig. 2(d)). Since the reference spectrum of µ was obtained with the same spectral resolution as the LCS spectra, we did not broaden this spectrum. The reference spectrum of µ was convolved with a Lorentzian with a line width of 6 nm (Δλ of the spectrograph).

4. Results

4.1. Simulation

Figure 4 shows the result of the simulation described in Section 3.3. The iD(λ) at ℓ1 and ℓ2 show little influence of the sample’s absorption (Fig. 4(a)). However, after frequency domain filtering of iD(λ), the resulting spectra S at ℓ1 and ℓ2 clearly show the presence of the oxyhemoglobin absorption peaks (Fig. 4(b)). From Fig. 4(c), we can conclude that our method of data acquisition and analysis for sdLCS fully recovers the input µ of this simulation. Even though the spectral features of the oxyhemoglobin absorption exhibit oscillations with a frequency comparable to that of the modulation on iD(λ), µ is recovered without any artifacts.

Fig. 4

Simulation of a µ measurement in sdLCS. a.) input spectra iD and b.) filtered spectra (S) at path lengths ℓ1 and ℓ2 inside the sample. c.) input µ and recovered µ.

4.2. Phantom

Figure 5 shows the measured µ spectra by sdLCS and tdLCS on the polystyrene-dye phantom. The µ spectra from both detection methods agree within the estimated uncertainties (error bars) for nearly all wavelengths. Also the fits on µ and the µ-contributions to the fits agree well, and the µ-contributions are in good agreement with the Mie calculated µ. The fitted dye concentrations were 10.1 ± 0.4% for the sdLCS measurement and 9.5 ± 0.4% for the tdLCS measurement, which are both very close to the expected dye concentration of 10%. Note that part of the difference in spectral shape is caused by the presence of Brownian motion influences for tdLCS, versus the absence of it for sdLCS (Section 3.4).

Fig. 5

Comparison of an sdLCS measurement to a tdLCS measurement of µ, µ and µ on a polystyrene-dye phantom.

5. Discussion and conclusion

In this article, we demonstrated and validated a new detection method and analysis algorithm for spectral domain detection in LCS. We showed that this method can be used for measuring local µ, µ and µ spectra in turbid media. The results of the sdLCS phantom measurement agreed well with the expected optical properties of the phantom and were comparable to the measured optical properties by tdLCS (Fig. 5). Whereas we did not explicitly show the localization aspect of our method by measuring on homogeneous media, previous studies with tdLCS proved that localized measurements of optical property spectra are possible within distinct tissue volumes [3].

The reason to validate this new detection method, is its theoretical sensitivity advantage compared to time domain detection (Section 2.4). As a consequence, the acquisition speed of the measurement can be reduced with respect to time domain detection, without reducing the signal to noise ratio. The sdLCS system that was validated in this article was more than three times faster than the tdLCS system, but not more sensitive (Section 3.2). A sensitivity advantage for sdLCS could not be demonstrated, because the quantum efficiency of the spectrograph is not optimized for these experiments. However, if the quantum efficiency would be similar to that of our tdLCS detector, the sensitivity advantage would be 11 dB at equal acquisition times for tdLCS and sdLCS (Eq. (5).

An additional shortcoming of the used spectrograph is its relatively slow acquisition speed (167 spectra/sec.). Slow acquisition speeds result in a decrease of the modulation depth in iD(λ) (‘fringe wash-out’), due to the reference mirror movement induced change of ΔOPL within the integration time. As a consequence, an additional decrease in sensitivity can occur, accounting for approximately 0.5 dB in our sdLCS system. Hence, not only higher quantum efficiency, but also faster acquisition speed is required to further improve the performance of our sdLCS system. Since high quantum efficiency line scan cameras are available with line rates up to 140 kHz, enhanced sensitivity and/or acquisition speed are realizable for sdLCS.

As we showed before for tdLCS, it may be advantageous to support the selection of the volume for optical property determination with an image of the investigated tissue volume [3]. Similar to tdLCS, this image can be reconstructed from the backscattered intensity as a function of depth. Our method for removal of the DC and modulation components in sdLCS also removes the ‘mirror image’ (complex ambiguity) when ΔOPL = 0 is situated inside the sample. As a consequence, this method for sdLCS can be used to reconstruct an image without cross-talk between iD(ΔOPL) and iD(–ΔOPL) within the investigated path length window inside the sample.

The only technique that is reported to have comparable performance to our sdLCS system in terms of localized measurements of optical properties, is dual window sOCT that was developed by Robles et al. [8,14]. Whereas Robles et al. report that their method achieves higher resolution in both the spectral and the spatial domain by avoiding the inherent tradeoff between the two, their system may suffer from the effects of the sensitivity roll-off of the detecting spectrograph and the absence of focus tracking, which influence measurement depth. An advantage of our sdLCS system is that it does not rely on refractive index-dependent corrections for the unwanted signal attenuation due to these effects. This facilitates the exact determination of µ and µ contributions to the measured µ. Since the current spatial and spectral resolution of our sdLCS system are comparable to the tdLCS system (Table 1), the resolution of our system is sufficient to obtain relevant determinations of local tissue chromophore concentrations, e.g. within the dermis, or epidermis only as shown in Ref. [3].

In conclusion, we have demonstrated and validated a new approach for sdLCS, which has high potential to improve the accuracy and speed of localized optical property measurements by LCS. Undoubtedly, this will lead to improved clinical utility of the technique, e.g. for the non-invasive determination of blood composition (hemoglobin/bilirubin concentration and oxygen saturation).