## A model for full-field optical coherence tomography in scattering media

Ugo Tricoli 1, 2 and R´ emi Carminati 2, ∗

1 ONERA, The French Aerospace Lab, Base A´ erienne 701,13661 Salon Cedex AIR, France 2 Institut Langevin, ESPCI Paris, CNRS, PSL University, 1 rue Jussieu, 75005 Paris, France

We develop a model of full-field optical coherence tomography (FF-OCT) that includes a description of partial temporal and spatial coherence, together with a mean-field scattering theory going beyond the Born approximation. Based on explicit expressions of the FF-OCT signal, we discuss essential features of FF-OCT imaging, such as the influence of partial coherence on the optical transfer function, and on the decay of the signal with depth. We derive the conditions under which the spatially averaged signal exhibits a pure exponential decay, providing a clear frame for the use of the Beer-Lambert law for quantitative measurements of the extinction length in scattering media.

## I. INTRODUCTION

Since its initial development [1], Optical Coherence Tomography (OCT) has proven its ability to image inside scattering materials with micrometer resolution in three-dimensions [2]. The possibility to access cellular structures in tissues at millimeter depth has been a breakthrough in biomedical optics, which stimulated a fast and broad dissemination of the technique.

The OCT setup is essentially a low coherence Michelson interferometer, in which one arm collects the light backscattered from the sample, while the other arm produces a reference beam reflected on a mirror. The main feature of an OCT setup is the ability to decouple the depth (longitudinal) and transverse resolutions [3]. Depth resolution is produced by temporal coherence gating, and is controlled by the spectral width of the incident light. Transverse resolution is controlled by the numerical aperture (NA) of the microscope objective in the sample arm. In scanning OCT (S-OCT) [1, 4-7], a point-by-point image is formed by three-dimensional scanning of a focal spot. OCT systems recording en face images in planes perpendicular to the optical axis have also been developed, using spatially coherent illumination as in wide-field OCT (WD-OCT) [8-10], or spatially incoherent illumination as in full-field OCT (FF-OCT) [11, 12].

Since OCT is expected to collect the singly backscattered photons, the signal is substantially affected by multiple scattering, whose contribution overcomes the signal at depths larger than the scattering mean free path [13]. Several approaches have been followed to decrease the mutiple scattering contribution by reducing the weight of long light paths, including spatial filtering through confocal detection [14], time gating [15], or polarization gating [16, 17]. Other strategies address an inverse problem to correct a posteriori for multiple scattering, and increase resolution and penetration depth. Approaches based on interferometric synthetic aperture microscopy [18], or computational adaptive optics [19], have proven to be successful. Interestingly, it has been shown that the effect of multiple scattering depends on the degree of spatial coherence of the illuminating beam [18]. It was also demonstrated that aberrations in the sample arm do not influence the transverse resolution in FF-OCT using spatially incoherent light [20]. Recently, a method based on the measurement of the reflection matrix has demonstrated an efficient discrimination between singly and multiply scattered light, with an unprecedented increase in the OCT working depth [21], thus pushing the limits of optical microscopy in highly scattering media [22].

Despite the success of OCT, from basic to clinical studies, a comprehensive theoretical model of the OCT signal, that handles a description of partial coherence together with a realistic scattering model (beyond the Born approximation), is still missing. In this paper we develop such a model, and use it to discuss different aspects of FF-OCT imaging. We study the influence of partial coherence on the optical transfer function. We also address the question of the decay of the signal with depth, that is captured by a mean-field scattering approach. Depending on the degree of spatial coherence and on the numerical aperture of the illumination/collection optics, we derive the conditions under which the spatially averaged FF-OCT signal exhibits a pure exponential decay with depth. The provides a clear frame for the use of the Beer-Lambert law for quantitative measurements of the extinction length in scattering media.

∗ remi.carminati@espci.fr

## II. SCATTERING MODEL

We consider a sample made of large-scale inhomogeneities immersed in a scattering medium made of randomly distributed scattering centers, as represented in Fig. 1(a).The sample is characterized by a dielectric function ε ( r ) that we write as ε ( r ) = ε b + ε d ( r ) + δε ( r ), where ε b denotes a uniform background, ε d ( r ) is a large-scale deterministic dielectric function varying on a scale L ε glyph[greaterorsimilar] λ , with λ the wavelength in vacuum, and δε ( r ) is a real random variable describing a disordered distribution of small-scale scattering centers. We assume that δε ( r ) satisfies 〈 δε ( r ) 〉 = 0, the brackets denoting a statistical average over an ensemble of realizations of the disordered scattering background, and that the correlation function 〈 δε ( r ) δε ( r ′ ) 〉 is of the form f ( | r -r ′ | /glyph[lscript] ε ), where f is a positive decaying function with range close to unity. This defines glyph[lscript] ε as the microscopic length scale of the disordered scattering medium, and we assume glyph[lscript] ε glyph[lessorsimilar] λ . Note that this description is very general. For example by setting ε d ( r ) = 0 we would describe a purely scattering medium, as that represented in Fig. 1(b) (OCT is sometimes used to measure the extinction length glyph[lscript] e in such materials). By setting δε ( r ) = 0 we would describe large-scale objects in a uniform background.

Figure 1: (a) Heterogeneous medium with large scale inhomogeneities immersed in a scattering background. (b) Purely scattering medium with small scale fluctuations of the dielectric function around a uniform value ε b . (c) Effective medium characterizing the propagation of the average field.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01049v2-figure-57fc2b69-d732-46a9-9533-ae7c51272eba.jpg)

In the statistical approach, in which the scattering background is treated as a random medium, the scattered field is a random variable. In a given realization of the sample (deterministic large-scale objects immersed in one statistical realization of the scattering background), the total field can be written E = 〈 E 〉 + δE . In OCT, one measures the backscattered field E s = E -E 0 , E 0 being the incident field, and we can also write E s = 〈 E s 〉 + δE s . In the regime of weak-scattering characterized by the condition | δE s | glyph[lessmuch] |〈 E 〉| , multiple scattering in the background scattering medium can be described in the mean-field approximation. This approximation is based on well-established result in multiple scattering theory [23], that are summarized in Appendix A. In the mean-field approach, the fluctuating part δE s of the scattered field at point r D in the detector plane, and at frequency ω , reads

$$\delta E _ { s } ( { \mathbf r } _ { D } ) = k _ { 0 } ^ { 2 } \int \langle G ( { \mathbf r } _ { D }, { \mathbf r } ^ { \prime } ) \rangle \Delta \varepsilon ( { \mathbf r } ^ { \prime } ) \langle E ( { \mathbf r } ^ { \prime } ) \rangle d ^ { 3 } r ^ { \prime } \,,$$

where k 0 = ω/c is the wavenumber in free space, 〈 G ( r D , r ′ ) 〉 is the average Green's function connecting a point r ′ in the sample to the detection point r D , and the integration is over the volume of the sample. Here the average Green's function 〈 G 〉 and the average field 〈 E 〉 are defined as solutions of a propagation equation in an effective medium, described by an effective dielectric function ε eff , representing the average contribution of the random scattering background (the detailed derivation of Eq. (1) is given in Appendix A). It is important to note that this effective medium description is a rigorous result of multiple scattering theory [23, 24]. Noting that ∆ ( ε r ′ ) = ε ( r ′ ) -ε eff is the local dielectric contrast between a heterogeneity in the sample and the effective medium, a clear physical meaning can be given to Eq. (1): It expresses the fluctuating part δE s of the scattered field as the result of single scattering in the effective medium, a result valid in the weak-scattering regime | δE s | glyph[lessmuch] |〈 E 〉| . The mean-field approach differs from the Born approximation by the renormalization of the background medium into an effective medium, that accounts in particular for the decay of the average Green's function due to scattering, as we shall see. The average field is connected to the incident field E 0 in the source plane by the relation

$$\langle E ( { \mathbf r } ) \rangle = A \int \langle G ( { \mathbf r }, { \mathbf r } _ { S } ) \rangle E _ { 0 } ( { \mathbf r } _ { S } ) \, d ^ { 2 } \rho _ { S }$$

where 〈 G ( r r , S ) 〉 is the average Green's function connecting a point r S = ( ρ S , z S ) in the source plane to an arbitrary point r , A is a constant that we do not need to specify, and the integration is along the source plane. Inserting Eq. (2) into Eq. (1) leads to

$$\delta E _ { s } ( { \mathbf r } _ { D } ) = A \, k _ { 0 } ^ { 2 } \int d ^ { 3 } r ^ { \prime } \int d ^ { 2 } \rho _ { S } \, \langle G _ { T } ( { \mathbf r } _ { D }, { \mathbf r } ^ { \prime } ) \rangle \Delta \varepsilon ( { \mathbf r } ^ { \prime } ) \langle G _ { T } ( { \mathbf r } ^ { \prime }, { \mathbf r } _ { S } ) \rangle E _ { 0 } ( { \mathbf r } _ { S } ) \,.$$

In this expression, for the sake of clarity, we have denoted by 〈 G T 〉 the Green's function that accounts for propagation through the optics in the sample arm, transmission at the effective medium surface, and propagation inside the effective medium. The average scattered field at the detector is readily deduced from Eq. (2), and reads

$$\langle E _ { s } ( { \mathbf r } _ { D } ) \rangle = A \int \langle G _ { R } ( { \mathbf r } _ { D }, { \mathbf r } _ { S } ) \rangle E _ { 0 } ( { \mathbf r } _ { S } ) \, d ^ { 2 } \rho _ { S } \,.$$

Here, the average Green's function 〈 G R 〉 accounts for propagation through the optics in the sample arm, and reflection at the effective medium surface. For practical calculations, the average Green's function in transmission or reflection can be approximated using a simple model, as we shall see below.

## III. FF-OCT SIGNAL

An OCT setup is based on a Michelson interferometer, as represented schematically in Fig. 2. Starting from the source, the beam is divided by a beam-splitter to travel along two arms. The sample arm collects the field backscattered from the sample. In the second arm, a reference beam is produced by reflection on a mirror. The detector collects the intensity resulting from the interference between the sample and reference beams. In FF-OCT, a full-field illumination is used in combination with an array of detectors (in practice a CCD camera) to record the signal at multiple transverse locations in parallel, and produce an en face image. This means that both arms contain a microscope objective, not represented in Fig. 2 for simplicity. The state of coherence of the light source plays a crucial role in FF-OCT. Beyond the longitudinal sectioning given by the temporal coherence length (or equivalently the spectral bandwidth), spatial coherence influences the lateral resolution, as well as the sensitivity to aberrations in the sample arm [20]. In order to understand precisely the role of temporal and spatial coherence, we need to build a model of the FF-OCT signal that includes the scattering model introduced in the previous section, and the main features of the interferometric and broadband detection used in OCT.

Figure 2: Schematic view of a FF-OCT setup using a source with partial temporal and spatial coherence. The detector measures the interferogram between the light backscattered from the sample, and the light reflected on the mirror. Illumination and collection optics (such as microscope objectives) in the different arms of the interferometer are not represented, since the model does not rely on a particular design. The z -axis of the reference frame is chosen so that the plane z = 0 coincides with the input surface of the sample (referred to as sample plane), with z &lt; 0 inside the sample. The z -axis follows the optical axis in each arm, with the source, reference mirror and detector planes corresponding to z = z S , z = z 0 and z = z D , respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01049v2-figure-8124680c-eae5-4731-93a3-f03240a70580.jpg)

Under illumination by statistically stationary light, the intensity recorded by a Michelson interferometer is characterized by three terms: The (time) average intensity of the reference field, the average intensity of the scattered field, and the cross-correlation between them (interference term). Assuming | E s | glyph[lessmuch] | E 0 | , the signal carrying information on the sample is the interference term (the intensity of the scattered field is negligible, and the intensity of the reference field only contributes as a background signal) [25]. In the frequency domain, the OCT signal measured at a given point r D in the detector plane is

$$S ( { \mathbf r } _ { D }, \omega ) = \overline { E _ { \text{ref} } ^ { * } ( { \mathbf r } _ { D } ) E _ { s } ( { \mathbf r } _ { D } ) } \,,$$

where E ref and E s are, respectively, the complex amplitude of the reference and backscattered field at frequency ω , the superscript ∗ denotes the complex conjugate, and the overline means a time averaging over the fluctuations of the partially coherent source.

The incident field E 0 at a point r = ( ρ , z ) can be written in the form of a plane-wave expansion

$$E _ { 0 } ( { \mathbf r } ) = \int d ^ { 2 } q \, e _ { 0 } ( { \mathbf q } ) \exp [ i { \mathbf q } \cdot { \mathbf \rho } - i { k } _ { z } ( { \mathbf q } ) z ] \,,$$

with k z ( q ) = ( k 2 0 -q 2 ) 1 / 2 . The z -axis is chosen such that the plane z = 0 coincides with the sample plane (the sample filling the half-space z &lt; 0), and the axis z &gt; 0 follows the optical axis in each arm (see Fig.1(d)). Assuming a weakly focused beam (the influence the numerical aperture is discussed in section V), we can perform a zeroth-order paraxial approximation k z ( q ) glyph[similarequal] k 0 , resulting in E 0 ( r ) = E 0 ( ρ , z = 0) exp( -ik 0 z ). The field in the reference arm is assumed to coincide with the incident field longitudinally shifted by the mirror displacement. Under the same paraxial approximation, we can write E ∗ ref ( r D ) = E ∗ 0 ( ρ D , z = 0)exp[ -ik 0 ( z D -2 z 0 )] where z 0 is the position of the mirror along the z -axis in the reference arm and z D is the position of the detection plane along the z -axis in the detection arm. From Eqs. (3), (4), (5) and the expressions of E 0 and E ∗ ref above , we can write the FF-OCT signal as

$$S ( { \mathbf r } _ { D }, \omega ) = \langle S ( { \mathbf r } _ { D }, \omega ) \rangle + \delta S ( { \mathbf r } _ { D }, \omega )$$

with

$$\langle S ( { \mathbf r } _ { D }, \omega ) \rangle \ = \ A \exp [ i k _ { 0 } ( 2 z _ { 0 } - z _ { S } - z _ { D } ) ] \int d ^ { 2 } \rho _ { S } \, \langle G _ { R } ( { \mathbf r } _ { D }, { \mathbf r } _ { S } ) \rangle W _ { 0 } ( \rho _ { S } - \rho _ { D } ) \,, \quad \ \.$$

$$\cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cd.$$

$$\delta S ( \mathbf r _ { D }, \omega ) \, = \, & \ A k _ { 0 } ^ { 2 } \exp [ i k _ { 0 } ( 2 z _ { 0 } - z _ { S } - z _ { D } ) ] \int d ^ { 3 } r ^ { \prime } \int d ^ { 2 } \rho _ { S } \, \langle G _ { T } ( \mathbf r _ { D }, \mathbf r ^ { \prime } ) \rangle \Delta \varepsilon ( \mathbf r ^ { \prime } ) \\ & \times \langle G _ { T } ( \mathbf r ^ { \prime }, \mathbf r _ { S } ) \rangle W _ { 0 } ( \rho _ { S } - \mathbf r _ { D } ) \,. \\ \intertext {throwsors in the form in the form } \, \text{throwsors in the form in the form in the form}$$

In these expressions we have introduced the cross-spectral density

$$W _ { 0 } ( \rho _ { S } - \rho _ { D } ) = \overline { E _ { 0 } ^ { * } ( \rho _ { D }, z = 0 ) E _ { 0 } ( \rho _ { S }, z = 0 ) } \,, & & ( 1 0 )$$

that characterizes the state of coherence of the incident field in the sample plane z = 0 [26]. Here we assume that the incoming intensity is uniform over the sample surface, and use a homogeneous Schell model such that W 0 is a function of ρ S -ρ D only [18].

It is clear that 〈 S ( r D , ω ) 〉 characterizes the effective medium, and does not carry information on ∆ ε ( r ′ ). In order to build the simplest model, we assume weak scattering, meaning that k glyph[lscript] 0 e glyph[greatermuch] 1 with glyph[lscript] e the extinction mean free path characterizing the decay of the average field (the exinction coefficient µ e = 1 /glyph[lscript] e can also be used equivalently). The effective medium, as seen by the average field, is described by a dielectric function ε eff = ε b + i/ k ( 0 glyph[lscript] e ) (this expression being valid to first order in the small parameter 1 / k glyph[lscript] ( 0 e )) [23, 24]. In this limit, the reflected average Green's function 〈 G R 〉 simply describes reflection at the surface of the effective medium, and contributes as a uniform background. Assuming a low index mismatch with ε b glyph[similarequal] 1, the contribution of 〈 S ( r D , ω ) 〉 becomes negligible. In the following we focus on the contribution δS ( r D , ω ) that carry the relevant information on the image formation process. It is useful to introduce the plane-wave expansion of the average Green's function:

$$\langle G _ { T } ( { \mathbf r }, { \mathbf r } ^ { \prime } ) \rangle = \frac { i } { 2 \pi } \int d ^ { 2 } q \frac { g ( { \mathbf q } ) } { k _ { z } ( q ) } \exp \left [ i { \mathbf q } \cdot ( \rho - \rho ^ { \prime } ) + i k _ { z } ( q ) | z - z ^ { \prime } | \right ].$$

We take k z ( q ) = ( k 2 0 -q 2 ) 1 / 2 for propagation outside the sample (for z, z ′ &gt; 0), and k z ( q ) = k eff z ( q ) = ( ε eff k 2 0 -q 2 ) 1 / 2 for propagation inside the sample (for z, z ′ &lt; 0). The filter g ( q ) plays the role of a pupil function that limits the transverse wavevector q within a region bounded by the numerial aperture NA of the objectives (the simplest model for g ( q ) is a disk with radius k 0 NA ) [27]. For g ( q ) = 1 (infinite pupil), we recover the Weyl expansion of the free-space Green's function. To get an expression of the FF-OCT signal relevant for an analysis in terms of optical transfer function, we insert Eq. (11) into Eq. (9), and perform again a zeroth-order paraxial approximation, using k z ( q ) glyph[similarequal] k 0 outside the sample, and k eff z ( q ) glyph[similarequal] k 0 + i/ (2 glyph[lscript] e ) inside the sample. The result is easily written in terms of the 2D Fourier transform of the signal, ˜ δS ( q , ω ) = ∫ δS ( r D , ω ) exp( -i q · ρ D ) d ρ 2 D / 4 π 2 , and reads

$$\widetilde { \delta S } ( { \mathbf q }, \omega ) = - 4 \pi ^ { 2 } A \int d z ^ { \prime } \exp [ 2 i k _ { 0 } & ( z _ { 0 } - z ^ { \prime } ) ] \, \widetilde { \Delta \varepsilon } ( { \mathbf q }, z ^ { \prime } ) \exp ( - | z ^ { \prime } | / \ell _ { e } ) \\ & \times \int d ^ { 2 } q ^ { \prime } g ( { \mathbf q } ^ { \prime } ) g ( { \mathbf q } + { \mathbf q } ^ { \prime } ) \widetilde { W } _ { 0 } ( { \mathbf q } ^ { \prime } ) \,,$$

$$( 1 2 )$$

where ∆ ( ˜ ε q , z ′ ) is the 2D Fourier transform of ∆ ε ( ρ ′ , z ′ ), and ˜ W 0 ( q ′ ) is the 2D Fourier transform of W 0 ( ρ S -ρ D ). Note that ˜ W 0 ( q ′ ) ∝ I 0 ( q ′ ), with I 0 ( q ′ ) the angular distribution of the incident intensity [26, 28]. Also note that this

expression does not rely on any assumption regarding the optics in the illumination, reference and detection arms, and does not depend on the definition of a focal plane for either illumination of detection (these features drive the precise form of the Green's function g ( q )). In practice, in order to control the degree of spatial coherence in the sample plane z = 0, one could choose to control the intensity distribution in the source plane z = z S and conjugate this plane with the plane z = 0, but this practical choice does not influence the general form of Eq. (12).

A feature of OCT is to integrate the signal over a broad spectral range. For a statistically stationary source, the broadband signal is obtained by integrating the different frequency components over the source bandwidth ∆ ω . Assuming that W 0 and ∆ ε have a weak dependence on ω , the spectral integration gives

$$\int _ { \omega _ { 0 } - \Delta \omega / 2 } ^ { \omega _ { 0 } + \Delta \omega / 2 } \exp [ 2 i k _ { 0 } ( z _ { 0 } - z ^ { \prime } ) ] d \omega = \, \exp [ 2 i \bar { k } _ { 0 } ( z _ { 0 } - z ^ { \prime } ) ] \sin \left [ ( z _ { 0 } - z ^ { \prime } ) \Delta \omega / c \right ] \Delta \omega \,,$$

where ¯ k 0 = ω /c 0 , with ω 0 the central frequency of the source. The sinc function, considered as a function of z ′ , is centered at z 0 with a width glyph[lscript] ω = c/ ∆ ω corresponding to the temporal coherence length of the incident light, and is responsible for the longitudinal sectioning in OCT. When glyph[lscript] ω glyph[lessmuch] glyph[lscript] e , we can use the approximation sinc[( z 0 -z ′ ) /glyph[lscript] ω ] glyph[similarequal] πglyph[lscript] ω δ z ( 0 -z ′ ). We endup with a closed-form expression of the 2D Fourier transform of the broadband OCT signal ˜ δS ( q ) = ∫ ∆ ω ˜ δS ( q , ω ) dω , that depends on the position z 0 of the mirror and on the source bandwith ∆ ω :

$$\widetilde { \delta S } ( { \mathbf q } ) = - 4 \pi ^ { 3 } A c \, \widetilde { \Delta \varepsilon } ( { \mathbf q }, z _ { 0 } ) \exp ( - z _ { 0 } / \ell _ { e } ) \int d ^ { 2 } q ^ { \prime } g ( { \mathbf q } ^ { \prime } ) g ( { \mathbf q } + { \mathbf q } ^ { \prime } ) \widetilde { W } _ { 0 } ( { \mathbf q } ^ { \prime } ) \,.$$

Expression (14) of the the FF-OCT signal implicitly involves different length scales, whose interplay is critical in the analysis of the signal: The length scale glyph[lscript] ε characterizing the microscopic random inhomogeneities in the sample, the length scale L ε characterizing large-scale deterministic fluctuations of the dielectric function, the temporal coherence length glyph[lscript] ω and the spatial coherence length glyph[lscript] c of the incident light. Expression (14) is similar to that previously derived in Ref. [18]. The main difference is that our derivation includes a scattering model based on a mean-field approach, that only relies on the weak-scattering assumption | δE s | glyph[lessmuch] |〈 E 〉| , always satisfied when the scattered field remains much smaller than the incident field. This approach goes beyond a first-order or second-order Born approximation, by accounting rigorously for the propagation of the average field in a renormalized effective medium. An interesting consequence is that the extinction of the signal with depth, due to scattering and absorption, emerges explicitly.

## IV. SPATIAL COHERENCE AND RESPONSE FUNCTION

Expression (14) is an interesting starting point for the analysis of the role of spatial coherence on the image formation in FF-OCT. To proceed, let us rewrite it in the compact form

$$\tilde { \delta S } ( \mathbf q ) = \widetilde { H } ( \mathbf q ) \, \widetilde { \Delta \varepsilon } _ { z _ { 0 } } ( \mathbf q ) \,,$$

where ˜ H ( q ) ∼ ∫ d q g 2 ′ ( q ′ ) g ( q + q ′ ) ˜ W 0 ( q ′ ) is the FF-OCT transfer function (we omit constant prefactors for simplicity), and ∆ ˜ ε z 0 ( q ) = ∆ ( ˜ ε q , z 0 ) exp( -z /glyph[lscript] 0 e ) is the weighted 2D Fourier transform of the dielectric contrast. The width of the cross-spectral density W 0 ( ρ ) defines the spatial coherence length glyph[lscript] c of the incident field in the sample plane z = 0. When glyph[lscript] c glyph[greatermuch] L ε glyph[greatermuch] glyph[lscript] ε , the incident light can be considered fully coherent, which corresponds to ˜ W 0 ( q ′ ) ∼ I 0 δ ( q ′ ), I 0 being proportionnal to the incident intensity. The transfer function for spatially coherent illumination becomes

$$\widetilde { H } _ { c } ( { \mathbf q } ) \sim I _ { 0 } \, g ( 0 ) g ( { \mathbf q } ) \,.$$

The coherent transfer function ˜ H c ( q ) posseses the same spatial frequency content as g ( q ), and covers the spatial frequency range q ≤ ¯ k NA 0 , leading to a transverse spatial resolution limit λ/ (2 NA ), where λ = 2 π/k ¯ 0 is the central wavelength of the incident field. The regime of spatially incoherent illumination corresponds formally to W 0 ( ρ ) ∼ δ ( ρ ), or equivalently ˜ W 0 ( q ′ ) = I ′ 0 . The transfer function for spatially incoherent illumination reads as

$$\widetilde { H } _ { i } ( \mathbf q ) \sim I _ { 0 } ^ { \prime } \, \int d ^ { 2 } q ^ { \prime } g ( \mathbf q ^ { \prime } ) g ( \mathbf q + \mathbf q ^ { \prime } ) \,.$$

In practice, since glyph[lscript] c cannot be made smaller than λ/ 2 under far-field illumination, the condition glyph[lscript] c glyph[lessmuch] glyph[lscript] ε glyph[lessmuch] L ε that rigorously corresponds to W 0 ( ρ ) ∼ δ ( ρ ) is out of reach. Nevertheless, provided that FF-OCT is used to image the large-scale inhomonegeities at the scale L ε , and not to resolve the small-scale random scattering centers, the condition

of incoherent illumination can be relaxed to be glyph[lscript] c glyph[similarequal] glyph[lscript] ε glyph[lessmuch] L ε . The integral term in Eq. (17) shows that the spatial frequency range encompassed by ˜ H i ( q ) is broader by a factor of two compared to coherent illumination, which is a usual result in optical microscopy, leading to a theoretical resolution limit λ/ (4 NA ). The results above, for coherent and incoherent illumination, are identical to those previsouly established in Ref. [18], and in agreement with the qualitative analysis of coherent and incoherent FF-OCT presented in Ref. [27].

Expressions (16) and (17) may also provide a theoretical frame to study the influence of aberrations in the sample arm on the image quality, recently discussed in Ref. [20]. A precise study is beyond the scope of the present work. Qualitatively, we can simply note that weak aberrations that would tend to squeeze the Green's function g ( q ) in Fourier space (without changing the cutoff frequency) would narrow the shape of ˜ H c ( q ), thus degrading the image quality. For incoherent illumination, the convolution product in Eq. (17) reduces the sensitivity of the transfer function to change in shapes due to weak aberrations, resulting in a better protection of the image quality.

## V. DEPTH DEPENDENCE OF THE SPATIALLY AVERAGED SIGNAL

The decay of the signal with depth is a feature of OCT. For paraxial illumination and detection, and in a medium with ∆ ( ˜ ε q , z 0 ) independent of z 0 , the signal follows an exponential decay exp( -z /glyph[lscript] 0 e ), as described by Eq. (14). This exponential decay can be used for the measurement of the extinction length in weakly scattering materials [29]. At higher numerical aperture, a correction to a pure exponential decay is expected, that may also depend on the degree of spatial coherence of the incident field. To address this question, we make use of a second-order paraxial approximation

$$k _ { z } ( q ) \, \simeq \, k _ { 0 } - \frac { 1 } { 2 } \frac { q ^ { 2 } } { k _ { 0 } } \,, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \$$

$$k _ { z } ^ { \text{eff} } ( q ) \ \simeq \ \sqrt { \varepsilon _ { \text{eff} } } k _ { 0 } - \frac { 1 } { 2 } \frac { q ^ { 2 } } { \sqrt { \varepsilon _ { \text{eff} } } k _ { 0 } } \simeq k _ { 0 } + \frac { i } { 2 \ell _ { e } } - \frac { 1 } { 2 } \frac { q ^ { 2 } } { k _ { 0 } } + \frac { i q ^ { 2 } } { 4 k _ { 0 } ^ { 2 } \ell _ { e } } \,,$$

in the plane-wave expansions of the incident field [Eq. (6)] and of the average Green's function [Eq. (11)]. Note that, as a result of the mean-field approach, the paraxial approximation in the scattering medium amounts to replacing k 0 by √ ε eff k 0 . This strongly influences the dependence of the FF-OCT on the numerical aperture, as we shall see. Following again the steps leading to Eq. (14), we end up with an expression of the FF-OCT signal that explicitely accounts for the angular aperture of the illumination and detection beams. In order to discuss the depth dependence of the spatially integrated signal, we consider ˜ δS ( q = 0), which can be cast in the following form (details of the derivation are given in Appendix B):

$$\widetilde { \delta S } ( { \mathbf q } = 0 ) = - 4 \pi ^ { 3 } A c \, \widetilde { \Delta \varepsilon } ( { \mathbf q } = 0, z _ { 0 } ) \exp ( - z _ { 0 } / \ell _ { e } ) \int d ^ { 2 } q ^ { \prime 2 } g ^ { 2 } ( { \mathbf q } ^ { \prime } ) \widetilde { W } _ { 0 } ( { \mathbf q } ^ { \prime } ) \exp [ - q ^ { \prime 2 } z _ { 0 } / ( 2 \bar { k } _ { 0 } ^ { 2 } \ell _ { e } ) ] \,,$$

where ¯ k 0 = ω /c 0 is the central wavenumber of the polychromatic incident field. This expression is well suited to a discussion of the depth dependence of the integrated FF-OCT signal. Note that in practice, an exponential fit to the signal is often used to estimate the extinction mean free path glyph[lscript] e in statistically homogeneous and isotropic scattering media. From Eq. (20), it is clear that at finite NA a deviation from a pure exponential decay exp( -z /glyph[lscript] 0 e ) may be observed due to the dependence on z 0 of the integral over q ′ . This integral contains several cutoffs. First, the last exponential term gives a depth-dependent cutoff q ′ z 0 glyph[similarequal] ¯ k 0 ( glyph[lscript] e /z 0 ) 1 / 2 . For z 0 glyph[greatermuch] glyph[lscript] e , the integral is restricted to vanishingly small q ′ , and a pure exponential decay exp( -z /glyph[lscript] 0 e ) is always expected in the tail of the signal versus depth. Second, the angular spectrum of the Green's function g ( q ′ ) produces a cutoff q ′ NA glyph[similarequal] ¯ k NA 0 , showing that for NA → 0 a pure exponential decay is observed, in agreement with Eq. (14). Third, the Fourier transform of the cross spectral density ˜ W 0 ( q ′ ) introduces a cutoff q ′ c glyph[similarequal] 2 π/glyph[lscript] c that depends on the degree of spatial coherence of the incident light in the sample plane. For glyph[lscript] c → ∞ (spatially coherent illumination), a pure exponential decay is also observed. More precisely, to prevent the exp[ -q ′ 2 z / 0 (2 ¯ k glyph[lscript] 2 0 e )] term in the integral to play a role, we need either NA glyph[lessmuch] ( glyph[lscript] e /z 0 ) 1 / 2 (low numerical aperture regime), or glyph[lscript] c glyph[greatermuch] λ z /glyph[lscript] ( 0 e ) 1 / 2 (coherent illumination).

This qualitative analysis is confirmed by numerical calculations of the spatially integrated signal using Eq. (20), as shown in Fig. 3. A non-exponential decay is observed in the regime NA glyph[similarequal] ( glyph[lscript] e /z 0 ) 1 / 2 glyph[similarequal] λ/glyph[lscript] c . The deviation from a pure exponential decay is enhanced at high NA and low spatial coherence. Although the curves are displayed in the regime z 0 glyph[similarequal] glyph[lscript] e , we have verified that an exponential decay exp( -z /glyph[lscript] 0 e ) is recovered in any case when z 0 glyph[greatermuch] glyph[lscript] e . Numerical calculations based on ab initio simulations could extend the analysis beyond the weak-scattering regime and second-order paraxial approximations used in this model. This is left for future work.

Figure 3: Spatially integrated FF-OCT signal ˜ δS ( q = 0) versus the normalized depth z /λ 0 for different values of the spatial coherence length glyph[lscript] c and numerical apertures NA . The curves are normalized by their value at z 0 = 0. The central wavelength of the incident light is λ = 800 nm, the temporal coherence length is glyph[lscript] ω = 1 µ m, and the extinction length is glyph[lscript] e = 8 λ . The pupil function g ( q ) is modeled by a Gaussian profile g ( q ) ∝ exp[ -q 2 / k NA ( ¯ 0 ) 2 ]. The cross-spectral density ˜ W 0 ( q ) in the Gaussian Shell-model is ˜ W 0 ( q ) ∝ exp[ -q 2 glyph[lscript] 2 c / (4 π 2 )]. (a): Signal for different degrees of spatial coherence with NA = 1. (b): Signal for different values of glyph[lscript] c and NA . The vertical dashed line corresponds to z 0 = glyph[lscript] e .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01049v2-figure-1dc915c8-fd30-4a6b-be82-6ad8fa1ccb82.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01049v2-figure-1bf3f2b3-79c2-4605-b960-fca817a82690.jpg)

## VI. CONCLUSION

In summary, we have developed a model of FF-OCT that includes a mean-field scattering theory, in addition to a precise description of temporal and spatial coherence. The model describes the decay of the signal with depth due to scattering and absorption, an essential feature in OCT imaging. It also permits an analysis of the interplay between different length scales characterizing the scattering medium and the degree of coherence of the incident field. Based on explicit expressions of the FF-OCT signal, we have discussed several features of FF-OCT imaging. We have analyzed the influence of spatial coherence on the optical transfer function, and discussed the particular cases of fully coherent and incoherent illuminations. We have also studied the depth dependence of the signal integrated over the transverse directions, that not only limits the penetration depth in OCT, but is also used to measure the extinction length in scattering materials. For spatially incoherent illumination, and/or with high numerical aperture of the illumination/detection optics, deviations from a pure exponential decay exp( -z /glyph[lscript] 0 e ) can be observed. A pure exponential decay is always recovered in the tail of the signal versus depth when z 0 glyph[greatermuch] glyph[lscript] e . Our analysis provides a clear frame for the use of the Beer-Lambert law for quantitative measurements of the extinction length in scattering media. The model also provides a framework for a precise analysis of the role of aberrations generated by the optics in the sample arm, or by the scattering medum itself, and for the development of advanced inverse reconstruction procedures going beyond the first and second-order Born approximations.

## Funding

This work was supported by the European Union's Seventh Framework Program under HELMHOLTZ grant agreement No 610110, and by LABEX WIFI (Laboratory of Excellence ANR-10-LABX-24) within the French Program 'Investments for the Future' under reference ANR-10-IDEX-0001-02 PSL . ∗

## Acknowledgments

We thank A. Aubry, V. Barolle, C. Boccara, R. Bocheux, K. Irsch and K. Plamann for stimulating discussions.

## Appendix A: Elements of multiple scattering theory

In order to build a mean-field expression of the scattered field, we first consider a purely scattering medium with dielectric function ε b + δε ( r ), as represented in Fig. 1(b). We assume that the random variable δε ( r ), that describes a spatial distribution of small-scale scattering centers, is statistically homogeneous and isotropic. An important result of multiple scattering theory, states that the average field obeys a propagation equation in an effective homogeneous medium, as represented in Fig. 1(c). We take this result as our starting point (for a derivation, see for example Refs. [23, 24]). More precisely, the average field obeys

$$\nabla ^ { 2 } \langle E ( { \mathbf r } ) \rangle + k _ { 0 } ^ { 2 } \varepsilon _ { \text{eff} } \left \langle E ( { \mathbf r } ) \right \rangle = 0 \,,$$

where ε eff is the effective dielectric function that describes the average contributio of the random scattering medium. Equivalently, the average Green's function satisfies the Dyson equation

$$\nabla ^ { 2 } \langle G ( { \mathbf r }, { \mathbf r } ^ { \prime } ) \rangle + k _ { 0 } ^ { 2 } \varepsilon _ { \text{eff} } \left \langle G ( { \mathbf r }, { \mathbf r } ^ { \prime } ) \right \rangle = - \delta ( { \mathbf r } - { \mathbf r } ^ { \prime } ) \,,$$

with an outgoing wave condition at infinity. This effective medium approach is strictly valid provided that glyph[lscript] ε glyph[lessmuch] λ , meaning that non-locality in the effective dielectric function can be neglected, glyph[lscript] ε being the microscopic clength scale associted to δε ( r ) [23]. Interestingly, the imaginary part of the effective dielectric constant describes the attenuation of the average field by scattering (and absorption), and defines the extinction mean free path glyph[lscript] e such that Im ε eff = ( k glyph[lscript] 0 e ) -1 . In absence of absorption, the extinction mean free path coincides with the scattering mean free path glyph[lscript] s .

In the general situation represented in Fig. 1(a), in which large-scale deterministic objects are superimposed to the random scattering background, the dielectric function is ε ( r ) = ε b + ε d ( r ) + δε ( r ). The total field obeys

$$\nabla ^ { 2 } E ( { \mathbf r } ) + k _ { 0 } ^ { 2 } \varepsilon ( { \mathbf r } ) E ( { \mathbf r } ) = 0 \,.$$

In order to obtain the equation satisfied by the fluctuating scattered field δE s = E -〈 E 〉 , we subtract Eq. (A1) to Eq. (A3), which leads to

$$\nabla ^ { 2 } \delta E _ { s } ( { \mathbf r } ) + k _ { 0 } ^ { 2 } \varepsilon _ { \text{eff} } \, \delta E _ { s } ( { \mathbf r } ) = - k _ { 0 } ^ { 2 } [ \varepsilon ( { \mathbf r } ) - \varepsilon _ { \text{eff} } ] E ( { \mathbf r } ) \.$$

Using the average Green's function 〈 G ( r r , ′ ) 〉 defined in Eq. (A2), the solution to Eq. (A4) is shown to obey the following integral equation:

$$\delta E _ { s } ( { \mathbf r } ) & = k _ { 0 } ^ { 2 } \int \langle G ( { \mathbf r }, { \mathbf r } ^ { \prime } ) \rangle [ \varepsilon ( { \mathbf r } ^ { \prime } ) - \varepsilon _ { \text{eff} } ] E ( { \mathbf r } ^ { \prime } ) d ^ { 3 } r ^ { \prime } \. & ( A 5 )$$

Assuming that scattering is sufficiently weak for the condition | δE s | glyph[lessmuch] |〈 E 〉| to be valid, the above expression can be approximated by

$$\delta E _ { s } ( { \mathbf r } ) = k _ { 0 } ^ { 2 } \int \langle G ( { \mathbf r }, { \mathbf r } ^ { \prime } ) \rangle [ \varepsilon ( { \mathbf r } ^ { \prime } ) - \varepsilon _ { \text{eff} } ] \langle E ( { \mathbf r } ^ { \prime } ) \rangle d ^ { 3 } r ^ { \prime } \.$$

Equation (A6) is a mean-field approximation to Eq. (A5). Physically, it describes the scattered field in one realization of the medium as resulting from a single scattering process in the effective medium.

## Appendix B: FF-OCT signal in the second-order paraxial approximation

In this appendix we develop the steps leading to Eq. (20). We start by performing the expansion (18) in Eq. (6), resulting in the following expressions of the incident and reference fields:

$$E _ { 0 } ( \mathbf r _ { S } ) = \exp ( \frac { - i k _ { 0 } z _ { S } ) } { \dots } \int d ^ { 2 } q \, e _ { 0 } ( \mathbf q ) \exp ( i \mathbf q \cdot \mathbf \rho _ { S } ) \exp [ i q ^ { 2 } z _ { S } / ( 2 k _ { 0 } ) ] \,,$$

$$E _ { r e f } ^ { * } ( \mathbf r _ { D } ) = \exp [ - i k _ { 0 } ( z _ { D } - 2 z _ { 0 } ) ] \int d ^ { 2 } q \, e _ { 0 } ^ { * } ( \mathbf q ) \exp ( - i \mathbf q \cdot \mathbf \rho _ { D } ) \exp [ i \mathbf q ^ { 2 } ( z _ { D } - 2 z _ { 0 } ) / ( 2 k _ { 0 } ) ] \,.$$

Next we use expansions (18) and (19) in the expressions (11) of the average Green's function. Here the expansion to second order in q is performed only in the exponential term, keeping a zeroth-order expansion k z ( q ) glyph[similarequal] k 0 in the

denominator that plays a minor role in the overall z dependence. Using the resulting expressions of E 0 , E ∗ ref and 〈 G T 〉 into Eqs. (3) and (5), we obtain

$$^ { \prime } \quad \cdot \quad \cdot \quad \cdot \quad \cdot \\ \widetilde { \delta S } ( { \mathbf q }, \omega ) = - 4 \pi ^ { 2 } A \int d z ^ { \prime } \exp [ 2 i k _ { 0 } ( z _ { 0 } - z ^ { \prime } ) ] \widetilde { \Delta e } ( { \mathbf q }, z ^ { \prime } ) \exp ( - | z ^ { \prime } | / \ell _ { e } ) \\ \times \exp [ i q ^ { 2 } ( z ^ { \prime } - z _ { D } ) / ( 2 k _ { 0 } ) ] \exp [ - q ^ { 2 } | z ^ { \prime } | / ( 4 k _ { 0 } ^ { 2 } \ell _ { e } ) ] \int d ^ { 2 } q ^ { \prime } g ( { \mathbf q } ^ { \prime } ) g ( { \mathbf q } + { \mathbf q } ^ { \prime } ) \widetilde { W } _ { 0 } ( { \mathbf q } ^ { \prime } ) \\ \times \exp [ i q q ^ { \prime } ( z ^ { \prime } - z _ { D } ) / k _ { 0 } ] \exp [ i q ^ { \prime } ( z ^ { \prime } - z _ { 0 } ) / k _ { 0 } ] \exp [ - ( q q ^ { \prime } + q ^ { \prime } ) | z ^ { \prime } | / ( 2 k _ { 0 } ^ { 2 } \ell _ { e } ) ] \,. \\ \colon \text{expression extends Eq. (12) to incident and detection beams with non-negligible numerical apertures.  We now}$$

This expression extends Eq. (12) to incident and detection beams with non-negligible numerical apertures. We now focus on the expression of the signal integrated over the transverse direction, which is obtained by taking q = 0:

$$\widetilde { \delta S } ( { \mathbf q } = 0, \omega ) & = - 4 \pi ^ { 2 } A \int d z ^ { \prime } \exp [ 2 i k _ { 0 } ( z _ { 0 } - z ^ { \prime } ) ] \, \widetilde { \Delta \varepsilon } ( { \mathbf q } = 0, z ^ { \prime } ) \exp ( - | z ^ { \prime } | / \ell _ { e } ) \\ & \int d ^ { 2 } q ^ { \prime } g ^ { 2 } ( { \mathbf q } ^ { \prime } ) \widetilde { W } _ { 0 } ( { \mathbf q } ^ { \prime } ) \exp [ i q ^ { \prime 2 } ( z ^ { \prime } - z _ { 0 } ) / k _ { 0 } ] \exp [ - q ^ { \prime 2 } | z ^ { \prime } | / ( 2 k _ { 0 } ^ { 2 } \ell _ { e } ) ] \,.$$

It can be verified by a numerical evaluation that the integral over q ′ weakly depends on frequency over a bandwidth ∆ ω corresponding to a depth resolution on the order of one micrometer, which corresponds to glyph[lscript] ω glyph[similarequal] 1 µm . Making use of Eq. (13) in the limit glyph[lscript] ω glyph[lessmuch] glyph[lscript] e , the signal integrated over frequencies becomes

$$\widetilde { \delta S } ( { \mathbf q } = 0 ) = - 4 \pi ^ { 3 } A c \widetilde { \Delta \varepsilon } ( { \mathbf q } = 0, z _ { 0 } ) \exp ( - z _ { 0 } / \ell _ { e } ) \int d ^ { 2 } q ^ { \prime } g ^ { 2 } ( { \mathbf q } ^ { \prime } ) \widetilde { W } _ { 0 } ( { \mathbf q } ^ { \prime } ) \exp [ - q ^ { \prime 2 } z _ { 0 } / ( 2 \bar { k } _ { 0 } ^ { 2 } \ell _ { e } ) ] \,,$$

where ¯ k 0 = ω /c 0 is the central wavenumber of the incident field. This expression is Eq. (20) of the main text.

- [1] D. Huang, E. A. Swanson, C. P. Lin, J. S. Schuman, W. G. Stinson, W. Chang, M. R. Hee, T. Flotte, K. Gregory, C. A. Puliafito et al. , 'Optical coherence tomography,' Science 254 , 1178-1181 (1991).
- [2] J. M. Schmitt, 'Optical coherence tomography (oct): a review,' IEEE J. Select. Topics Quantum Electron. 5 , 1205-1215 (1999).
- [3] A. F. Fercher, W. Drexler, C. K. Hitzenberger, and T. Lasser, 'Optical coherence tomography-principles and applications,' Rep. Prog. Phys. 66 , 239 (2003).
- [4] N. Nassif, B. Cense, B. H. Park, S. H. Yun, T. C. Chen, B. E. Bouma, G. J. Tearney, and J. F. de Boer, 'In vivo human retinal imaging by ultrahigh-speed spectral domain optical coherence tomography,' Opt. Lett. 29 , 480-482 (2004).
- [5] R. Huber, M. Wojtkowski, and J. Fujimoto, 'Fourier domain mode locking (fdml): A new laser operating regime and applications for optical coherence tomography,' Opt. express 14 , 3225-3237 (2006).
- [6] B. M. Hoeling, A. D. Fernandez, R. C. Haskell, E. Huang, W. R. Myers, D. C. Petersen, S. E. Ungersma, R. Wang, M. E. Williams, and S. E. Fraser, 'An optical coherence microscope for 3-dimensional imaging in developmental biology,' Opt. Express 6 , 136-146 (2000).
- [7] Y. Zhang, J. Rha, R. S. Jonnal, and D. T. Miller, 'Adaptive optics parallel spectral domain optical coherence tomography for imaging the living retina,' Opt. Express 13 , 4792-4811 (2005).
- [8] S. Bourquin, P. Seitz, and R. Salath´, 'Optical coherence topography based on a two-dimensional smart detector array,' e Opt. Lett. 26 , 512-514 (2001).
- [9] E. Bordenave, E. Abraham, G. Jonusauskas, N. Tsurumachi, J. Oberle, C. Rulliere, P. Minot, M. Lassegues, and J. S. Bazeille, 'Wide-field optical coherence tomography: imaging of biological tissues,' Appl. Opt. 41 , 2059-2064 (2002).
- [10] M. Laubscher, M. Ducros, B. Karamata, T. Lasser, and R. Salathe, 'Video-rate three-dimensional optical coherence tomography,' Opt. Express 10 , 429-435 (2002).
- [11] L. Vabre, A. Dubois, and A. C. Boccara, 'Thermal-light full-field optical coherence tomography,' Opt. Lett. 27 , 530-532 (2002).
- [12] A. Fercher, C. Hitzenberger, M. Sticker, E. Moreno-Barriuso, R. Leitgeb, W. Drexler, and H. Sattmann, 'A thermal light source technique for optical coherence tomography,' Opt. Commun. 185 , 57-64 (2000).
- [13] J. M. Schmitt and A. Kn¨ttel, 'Model of optical coherence tomography of heterogeneous tissue,' u J. Opt. Soc. Am. A 14 , 1231-1242 (1997).
- [14] R. H. Webb, 'Confocal optical microscopy,' Rep. Prog. Phys. 59 , 427 (1996).
- [15] M. R. Hee, J. A. Izatt, J. M. Jacobson, J. G. Fujimoto, and E. A. Swanson, 'Femtosecond transillumination optical coherence tomography,' Opt. Lett. 18 , 950-952 (1993).
- [16] J. Schmitt, A. Gandjbakhche, and R. Bonner, 'Use of polarized light to discriminate short-path photons in a multiply scattering medium,' Appl. Opt. 31 , 6535-6546 (1992).

- [17] C. M. Macdonald, U. Tricoli, A. Da Silva, and V. A. Markel, 'Numerical investigation of polarization filtering for direct optical imaging within scattering media,' J. Opt. Soc. Am. A 34 , 1330-1338 (2017).
- [18] D. L. Marks, B. J. Davis, S. A. Boppart, and P. S. Carney, 'Partially coherent illumination in full-field interferometric synthetic aperture microscopy,' J. Opt. Soc. Am. A 26 , 376-386 (2009).
- [19] S. G. Adie, B. W. Graf, A. Ahmad, P. S. Carney, and S. A. Boppart, 'Computational adaptive optics for broadband optical interferometric tomography of biological tissue,' Proc. Natl. Acad. Sci. (USA) (2012).
- [20] P. Xiao, M. Fink, and A. C. Boccara, 'Full-field spatially incoherent illumination interferometry: a spatial resolution almost insensitive to aberrations,' Optics letters 41 , 3920-3923 (2016).
- [21] A. Badon, D. Li, G. Lerosey, A. C. Boccara, M. Fink, and A. Aubry, 'Smart optical coherence tomography for ultra-deep imaging through highly scattering media,' Science advances 2 , e1600370 (2016).
- [22] A. Badon, A. C. Boccara, G. Lerosey, M. Fink, and A. Aubry, 'Multiple scattering limit in optical microscopy,' Opt. Express 25 , 28914-28934 (2017).
- [23] P. Sheng, Introduction to Wave Scattering, Localization and Mesoscopic Phenomena (Springer, 2010).
- [24] E. Akkermans and G. Montambaux, Mesoscopic Physics of Electrons and Photons (Cambridge University Press, 2007).
- [25] B. J. Davis, T. S. Ralston, D. L. Marks, S. A. Boppart, and P. S. Carney, 'Autocorrelation artifacts in optical coherence tomography and interferometric synthetic aperture microscopy,' Opt. Lett. 32 , 1441-1443 (2007).
- [26] L. Mandel and E. Wolf, Optical Coherence and Quantum Optics (Cambridge University Press, 1995).
- [27] A. Sentenac and J. Mertz, 'Unified description of three-dimensional optical diffraction microscopy: from transmission microscopy to optical coherence tomography: tutorial,' J. Opt. Soc. Am. A 35 , 748-754 (2018).
- [28] Y. N. Barabanenkov and V. Finkelberg, 'Radiation transport equation for correlated scatterers,' Sov. Phys. JETP 26 , 587-591 (1968).
- [29] R. Bocheux, P. Pernot, V. Borderie, K. Plamann and K. Irsch, 'Quantitative measures of corneal transparency, derived from objective analysis of depth-resolved corneal images, demonstrated with full-field optical coherence tomographic microscopy,' PLoS ONE 14 , e0221707 (2019).
Ugo Tricoli Rémi Carminati
A model for full-field optical coherence tomography in scattering media
We develop a model of full-field optical coherence tomography (FF-OCT) that
includes a description of partial temporal and spatial coherence, together with
a mean-field scattering theory going beyond the Born approximation. Based on
explicit expressions of the FF-OCT signal, we discuss essential features of
FF-OCT imaging, such as the influence of partial coherence on the optical
transfer function, and on the decay of the signal with depth that is captured
by the model. We derive the conditions under which the spatially averaged
signal exhibits a pure exponential decay with depth, providing a clear frame
for the use of the Beer-Lambert law for quantitative measurements of the
extinction length in scattering media.