The data extracted from individual OCT scans included both quantitative and qualitative measures, and at the current level of resolution, these single factors appear insufficient for diagnosis. Our approach suggests that it may be possible to extract diagnostic data from the overall architecture of the OCT images with a reasonable diagnostic accuracy when used in combination.

A Monte Carlo (MC) method for modeling optical coherence tomography (OCT) measurements of a diffusely reflecting discontinuity embedded in a scattering medium is presented. For the first time to the authors' knowledge it is shown analytically that the applicability of an MC approach to this optical geometry is firmly justified, because, as we show, in the conjugate image plane the field reflected from the sample is delta-correlated from which it follows that the heterodyne signal is calculated from the intensity distribution only. This is not a trivial result because, in general, the light from the sample will have a finite spatial coherence that cannot be accounted for by MC simulation. To estimate this intensity distribution adequately we have developed a novel method for modeling a focused Gaussian beam in MC simulation. This approach is valid for a softly as well as for a strongly focused beam, and it is shown that in free space the full three-dimensional intensity distribution of a Gaussian beam is obtained. The OCT signal and the intensity distribution in a scattering medium have been obtained for several geometries with the suggested MC method; when this model and a recently published analytical model based on the extended Huygens-Fresnel principle are compared, excellent agreement is found.

A recently developed analytical optical coherence tomography (OCT) model [Thrane et al., J. Opt. Soc. Am. A 17, 484 (2000)] allows the extraction of optical scattering parameters from OCT images, thereby permitting attenuation compensation in those images. By expanding this theoretical model, we have developed a new method for extracting optical scattering parameters from multilayered tissue structures in vivo. To verify this, we used a Monte Carlo (MC) OCT model as a numerical phantom to simulate the OCT signal for heterogeneous multilayered tissue. Excellent agreement between the extracted values of the optical scattering properties of the different layers and the corresponding input reference values of the MC simulation was obtained, which demonstrates the feasibility of the method for in vivo applications. This is to our knowledge the first time such verification has been obtained, and the results hold promise for expanding the functional imaging capabilities of OCT.

Analytical and numerical models for describing and understanding the light propagation in samples imaged by optical coherence tomography (OCT) systems are presented. An analytical model for calculating the OCT signal based on the extended Huygens-Fresnel principle valid both for the single and multiple scattering regimes is reviewed. An advanced Monte Carlo model for calculating the OCT signal is also reviewed, and the validity of this model is shown through a mathematical proof based on the extended Huygens-Fresnel principle. Moreover, for the first time the model is verified experimentally. From the analytical model, an algorithm for enhancing OCT images is developed: the so-called true-reflection algorithm in which the OCT signal may be corrected for the attenuation caused by scattering. For the first time, the algorithm is demonstrated by using the Monte Carlo model as a numerical tissue phantom. Such algorithm holds promise for improving OCT imagery and to extend the possibility for functional imaging.

No abstract