|
1.IntroductionLight-based techniques are becoming an increasingly popular method of probing heavily scattering media such as body tissue (for a review see Ref. 1, also see Ref. 2). For surface tissue or thin tissue samples it is possible to use the coherence or polarization properties of the light.3 4 5 6 7 8 9 At the same time, the level of depolarization of initially polarized light, as well as the mode of changing the polarization from one type to another10 11 12 13 14 15 remain the informative parameters characterizing biological tissue (BT) structure. It was determined that the length of the pathway of light, on which its complete depolarization occurs, is a characteristics of BTs of different morphological structures and pathological states. It was determined that light can propagate in skin with a length of approx. 1.2 mm, retaining its linear polarization. This enables obtaining informational images of macroinhomogeneities11 12 13 and simultaneous recording of the birefringence and optical activity14 16 17 18 19 20 21 22 of BT. A reasonable choice of the laser beam polarization state remains an urgent task. According to the data obtained13 different BTs or the same BTs in various pathological or functional states present different reactions to probing by light with linear and circular polarization. At the same time, the problems of polarization selection of images of different layers of BT with a simultaneous increase of the relation to “signal noise” in the image of optically anisotropic structures, which are screened by the other BT layer, remain inadequately explored. This can be used in novel elaborations in the field of optical coherence tomography (OCT), connected with direct recording of Mueller matrix images of successively located BT layers. Therefore, it would be important to extend the theoretical and experimental background of studies of optical properties of biotissues using polarized light at least in two ways. First, it might be important to expand our understanding of the transformation of polarized radiation by biotissues using the formalism of the Stokes vector and the Muller matrix. The second urgent task is to analyze the vector structure of random laser fields aimed at subsequent development of methods for the polarization selection of such fields. 1.1.Objective of the PaperThis paper studies the possibilities of using polarization methods to improve the SNR in images of optically anisotropic structures (collagen and myosin networks, hematomas, neoplasms, etc.) in the case when the biotissue being probed is screened by another biotissue. We consider two situations: singly and multiply scattering object scenes. 2.Theoretical Modeling2.1.Monolayers—Single Scattering2.1.1.General remarksOn the basis of information about morphological biotissue structure,23 24 25 the following model scheme is suggested. Biotissue generally consists of two phases: amorphous and optically anisotropic ones. The optically anisotropic component of the biotissue possesses two levels of organization—crystalline and architectonics ones. Optically coaxial organic fibrillae forming collagen, elastin, myosin fibers, mineralized (hydroxylapatite crystals) and organic fibers, etc. can be included in the crystallite level. Architectonic nets of such biotissues as skin derma (SD), muscle (MT), and bone tissues (BnT) are morphologically the most typical. The architectonics of SD is formed by statistically oriented bundles of collagen fibrillae (Fig. 1). The diameter of the fibrilla is between 0.5 and 2 μm. Parallel bundles of fibrillae form a fiber, the diameter of which is 5 to 7 μm (papilary layer23), and in the net layer, it reaches 30 μm. A collagen bundle, the average diameter of which changes from 100 to 200 μm, is considered to be the highest structural unit of optically active collagen (birefringence value Δn≈10−3) . BnT represents a system, consisting of the trabeculae [Fig. 2(a.1)] layer and osteons [Fig. 2(a.2)]. The optically active matrix consists of hydroxylapatite crystals [Δn≈10−1 (Ref. 24)], the long (optical) axes of which are oriented along the longitudinal axis of collagen fibers [Fig. 2(b)]. They are located between microfibrillae, fibrillae, and collagen fibers, forming a separate continuous mineral phase. Collagen fibers are spatially armored elements in the mineral matrix. The bone trabeculae fiber orientation is ordered and parallel to their plane [Fig. 2(c)]. For BnT osteons, the spatially spiral orientation of the armored collagen fibers is realized. MT is a structuralized, spatially ordered system of protein bundles, consisting of optically isotropic actin and anisotropic myosin25 (Fig. 3). 2.1.2.Optical characteristicsOptical properties of the primary (crystallite) level of BT organization are characterized by a Mueller operator of the following kind23: Here ρ is the orientation of optical axis resulting from the direction of the packing of anisotropic fibrillae (collagen myosin, elastin, hydroxyapatite, etc.); and δ is the value of phase shift, performed by their substance between ordinary and nonordinary waves.2.2.Two-Layer Scattering MediumRecently applications of the Monte Carlo technique in optics of BT have been successfully developed.26 27 28 29 The known algorithms enable us to take into account several BT layers with various optical properties. With its high accuracy and universality, the time factor remains the main drawback of Monte Carlo techniques. The elaboration of new Mueller matrix OCT (MMOCT) techniques30 31 32 33 34 35 36 37 38 requires the creation of efficient, comparatively practical methods of probing the polarization of BT, based on the mathematical formalisms of Stokes vectors and Mueller matrices. Let a light field, described by Stokes vector S (0), propagate within a two-layer optically heterogeneous medium (Fig. 4). Due to light-scattering on the heterogeneity of the first layer, the biotissue under study is irradiated by a beam with another Stokes vector S (1)=Y S (0), where Y is a matrix of radiation light scattering resulting from the screening biotissue layer. During the interaction of radiation with the subsurface layer an object field is formed, the polarization characteristic features of which can be described by the vector S (2)=XY S (0). Here X is a matrix of the light scattering of the investigated biotissue. The radiation detector measures the light field, the Stokes vector of which takes the following form: S (X)=YXY S (0) and S (Y)=Y S (0). If one takes into account the multiple acts of object reirradiation [S (3)] and the further outlet of the component of luminous flux into external medium [S ( * )], it is possible to write it as follows: S ( * )=YXYXY S (0). Proceeding from the assumption of the independence of the statistical structure of biotissue layers the resulting light field can be presented by the partial vector sum: which takes the following matrix form: where Z is the resulting matrix of light-scattering by a two-layer object scene.First let us consider the first two constituent parts, i.e., let us assume From an analysis of Eqs. (3) and (4), one can conclude that under determined conditions—i.e., having the information concerning matrices Z and Y—it is possible to define the matrix X structure statistically: Here Y −1 is the inverse matrix of Y.2.2.1.Investigation of transformation of types and forms of laser field polarization in heterogeneous object scenesLet us consider the possibilities of determining polarization structure of laser radiation scattered by biotissue layers. It is successfully possible to demonstrate the determination of polarization characteristic features of partial scattered fluxes (αxβx;αy,βy) from the following equations:23 where α0 and β0 are polarization azimuths and ellipticity of the probing laser beam, and fik are the elements of the partial light-scattering matrices.Thus, if we have the information concerning matrices Z and Y, we can “construct” matrix X, i.e., one is able to determine polarization structure of a signal (αx,βx), reflected by the object being probed. 2.2.2.Technique of polarization—phase information selection in random fields of coherent radiationThe idea of the technique of polarization-phase extraction of an object image against a background of impediments is based on determining polarization structure (αx,βx;αy,βy) of the resulting light field S (2) by using the measured matrices Y and Z with the further polarization compensation of a background signal component S (Y). One can essentially decrease the intensity value of the scattered background polarization component while the intensity of object signal remains considerable. The latter is possible by placing a polarization filter before the photodetector. The filter consists of an analyzer and quarter-wave plate. Figure 5 demonstrates the mechanism of increasing the relation of “signal to noise” (RSN) being determined as follows: where αy * is the azimuth of background signal polarization after passing through phase plate λ/4; and αx * and βx * are informative signals for the azimuth and ellipticity after passing through the quarter-wave plate. The axis of the latter is oriented in a special way, thus, the ellipticity of background signal polarization βy *=0; the constant 103 characterizes the linearity of polarizer-analyzer.By orienting the greatest speed axis at the angle Ω=αy * one can obtain the following polarization characteristic features of the informative and background signals: Thus, the RSN within the optical field under registration has the value where IS and IN are the intensities of informative and background signals, respectively (Fig. 5); and PS, ΔS, PN, and ΔN are the polarization and depolarization levels of these components of the scattered radiation field.By orienting the greatest speed axis of the phase element at the angle π/2+αx *, it is possible to increase the RSN level like this: where k=10−3 is the level of radiation intensity of a completely polarized signal component passed through the analyzer; and ζ is the angle between the fluctuation area of electrical intensity vector of completely polarized constituent part of informative signal and analyzer axis direction A (Fig. 5). In addition 50% of the intensity of the completely depolarized constituent parts is passed through the analyzer. Thus, the increase of the RSN ratio resulting from using the polarization compensation technique has the following value: where An analysis of Eqs. (15) and (16) shows that maximal level of the RNS is achieved if the following conditions are fulfilled:3.Computer Modeling, Analysis, and Discussion3.1.Single ScatteringFigure 7 presents the results of computer modeling the distributions of the coefficient η(α0,Θ) for an object scene with the following parameters of birefringence: Δnx=0.1, Δny=0.03 [Fig. 7(a)], and 0.08 [Fig. 7(b)], and the preferable directions of orientations of biotissues architectonics fibers: Δρ=ρx−ρy=20 deg [see Eq. (1)]. The results of these calculations enable us to make the following conclusions: l1.
Let us analyze these results. As follows from boundary conditions [Eq. (17)], the extrema increase in the SNR is achieved when the background signal is linearly polarized. In this case, the intensity of the background signal can be reduced to the level of μ−1IN. Such a background suppression is achieved when the azimuth α0 coincides with the anisotropy direction ρy of screening biotissue (this direction is measured from ρx) , i.e., α0=Δρ. Thus, an extreme line should exist in the distribution η(α0,Θ) for all the scattering angles Θ. On the other hand, Eqs. (9) to (12) show that the level of the background and informative signals considerably depends on the relations between the partial matrix elements {xik} and {yik}. The values of these elements fluctuate in different ways, depending on Δnx, Δny, and Θ (Ref. 23). Therefore, a linearly polarized background and the relevant extremum lines η(α0,Θ)/Θi=const may be observed within a broad range of α0 [Figs. 7(b)]. Figure 8 displays the dependences characterizing the transformation of the polarization structure of the background (IN) and informative (IS) signals in response to variations in the azimuth α0 of the probing beam. These results were obtained for the following scattering angles: Θ=0 deg [Figs. 8(a) and 8(d)], Θ=35 deg [Figs. 8(b) and 8(e)], and Θ=55 deg [Figs. 8(c) and 8(f)]. The left-hand column corresponds to anisotropy parameters Δnx=0.1 and Δny=0.03, while the right-hand column corresponds to anisotropy parameters Δnx=0, and Δny=0.08. The analysis of these data reveals a correlation between zero values of the ellipticity of the background signal βy *(α0=Δρ) (curves 3) and the extremal level of the improvement factor of the SNR ratio η(α0,Θ)→μ, which is observed for arbitrary values of polarization parameters of the informative signal αx *−αy *=Δα and βx * (curves 1 and 2 in Fig. 8). For other values of the azimuthal angle α0 of polarization of the probing beam, two field components of scattered laser radiation are generally elliptically polarized. Therefore, the level of the coefficient η(α0,Θ) becomes much lower, vanishing around βx *≈βx *. The choice of the optimal scattering angle Θ (with α0=Δρ) and the optimal polarization plane for scattering angles corresponding to the extremal situation, [η(α0,Θ)]/Θi=const , is important for the photometry of the informative signal IS. This is due to the fact that, when the extremal level η(α0,Θ)=μ is achieved, the polarization structure of the informative signal differs for different values of α0 and Θ. Figure 8 displays a series of dependences of IS calculated for the directions α0=20 deg [Fig. 8(a)], Θ=35 deg [Fig. 8(b)], and Θ=55 deg [Fig. 8(c)] of extrema of the coefficient η(α0,Θ). Analysis of these results shows that, when the minimum background level (μ−1IN) is achieved, the informative signal may vary within a broad dynamic range of three orders of magnitude (from 10−3 to 1). An optimal photometric situation is achieved when the boundary conditions for the polarization state of the object signal produced by a biotissue [Eqs. (17)] are satisfied. 3.2.Multiple ScatteringLet us reconsider the analysis of Eqs. (15) and (16). The bright dependency of η curves on the polarization condition of the probing beam (α0,β0) angle characteristics makes the necessity to solve the tasks of polarization optimization, i.e., selection of optimal characteristic features of the probe beam when T is a maximum obtained value. Let us consider first the situation of linearly polarized scattering fields. Let us study function η(α0) for its extremum: One can conclude from Eq. (18) that for polarizational probing of two-layer scattered medium the polarizationally optimal ones will be Eqs. (17). Equations (17) are possible when the polarization azimuth of the probing laser beam satisfies the following conditions: Analogous results can be achieved by using polarization ellipticities of the probing beam at the constant value of polarization azimuth (α0=const).The optimal value will be that of polarization ellipticity (Fig. 9), which is defined by the following equation solution: As we can see from Eq. (15), the depolarization parameter Δ noticeably lowers the improvement factor of the SNR in multiply scattering object scenes. This tendency is shown in Fig. 6, which illustrates the lowering in the improvement factor of the SNR with a decrease in the polarization degree on a logarithmic scale. As we can see from this figure, for Π=(1−Δ)=0.9, the coefficient η is reduced by two orders of magnitude, while for Π=0.1, the factor η decreases further, virtually by another order of magnitude.On the other hand, the analysis of Eq. (15) shows that the improvement factor of the SNR reaches its extremal level when both conditions of Eq. (17) are satisfied. 3.2.1.Computer modelingFigure 10 presents the results of computer modeling for the “spatial” and “topological” distributions of the improvement factor of the signal-to-noise ratio η obtained in the cases of single [Figs. 10(a) and 10(b), Π=1.0] and multiple [Figs. 10(c) and 10(d), Π=0.12) scattering and in the regime of polarization correction of the probing beam [Figs. 10(e) and 10(f)]. The parameters of the object scene model a sequence of layers of cancellous bony and fibrillar muscular tissues with Δnx=0.1, Δny=0.03, and Δρ=ρx−ρy=00. The results of our modeling enable us to make the following conclusions: l2.
4.Automatic Polarimeter DesignThe overall layout of the instrument is shown in Fig. 12. The main optical elements of the polarimeter are the He-Ne laser (λ=0.6328 μm, 10 mW), the polaroid P, two quartz quarter-wave (λ/4) plates (fourth order), and the polaroid A as the analyzer. By setting the polarization plane of P parallel to the optical axis of the first quarter-wave plate the object under test O (a tissue sample) is illuminated by a linearly polarized light beam, with the analyzer similarly coupled with the second quarter-wave plate. The He-Ne laser beam is expanded by the afocal system (R) up to the diameter of 3 mm. The pinhole located in the expansion optics removes the scattered light originating from secondary reflections in the resonator and optics contaminations. The role of the half-wave plate λ/2 is to match the laser beam polarization plane to the polarizer transmittance direction and, at the same time, to control the light level at the CCD matrix to avoid its saturation. The tissue sample is placed on the microscope slide. The slide with tissue is fixed in a mount, providing transverse translations in two mutually perpendicular directions in the range of 5 mm to enable us to choose various fragments of the sample. The sample image is formed on the CCD matrix by the optical system composed of two microscope objectives OL. The overall dimensions of the quarter-wave plate and analyzer driving systems do not enable us to place the imaging optics closer than 200 mm from the sample. The images registered by the CCD matrix and the frame grabber are transferred in a digital form to the computer memory. The half-wave plate, both quarter-wave plates, the polarizer, and the analyzer are rotated by means of step motors. One motor step corresponds to 22.5 arcmin of the rotation of an optical element. Before making measurements, the automatic system calibration is performed in the following order: l3.
5.Experimental ResultsFigure 13 presents the results of experimental testing of the method of polarization laser probing of biotissues. Histological slices of biostructures placed in a sequence were employed as objects of our studies. The optical and morphological parameters of these objects correlated with the data of computer modeling. As the objects under investigation, the following types of morphological sections of biotissues were used:39 40 41 42 l4.
Experiments were performed with the following object scenes, skin derma and cancellous bony tissue for the following parameters Δρ≈20 deg, α0=20 deg [Fig. 13(a)], and Θ=0 deg [Fig. 13(b)], and fibrillar muscular tissue and derma tissue with Δρ≈20 deg, α0=20 deg [Fig. 13(c)], Θ=0 deg [Fig. 13(d)], Θ=35 deg [Fig. 12(e)], and Θ=55 deg [Fig. 13(f)]. The results of measurements of the coefficient η(α0,Θ) demonstrate the efficiency of polarization probing of biotissues within a sufficiently broad range of scattering angles, Θ<20 deg. The discrepancy between the results of computer modeling (solid lines) and experimental data (points) does not exceed 10 to 20% within this range. As the observation angle Θ increases, the coefficient η(α0,Θ) noticeably decreases, while the discrepancy between theoretical predictions and experimental data may reach 50 to 80%. This effect is apparently due to the increase in the multiplicity of light scattering and the growth in the depolarization of the background signal, which is not included in theoretical modeling. However, even in this regime, the improvement factor of the SNR remains sufficiently high: η(α0,Θ)∼10 to 100. The results of experimental testing of the method of laser probing with a fixed polarization state and correction of the azimuth and ellipticity of polarization are presented in Fig. 14. In the figure, curves 1 to 3 were obtained for scattering angles Θ=0, Θ=20, and Θ=55 deg, respectively, illustrating the efficiency of polarization probing of biotissues with a high level of the depolarized background signal. The thickness of shading biotissue (skin derma) was z=250 μm. The results of laser probing with polarization correction are presented by curve 4. A comparative analysis of the experimental data reveals the growth in the improvement factor of the SNR by virtually an order of magnitude for any angle Θ. A comparison with the results of analytical modeling shows a satisfactory agreement between theoretical predictions and experimental data. The level of discrepancies varies from 10 to 30%. Lower SNRs obtained in experiments may be due to some variance of orientations of optically active elements in bony and muscular tissues, which may increase the ellipticity of the object field19 and a variance of the optical anisotropy of biotissues. REFERENCES
J. C. Hebden
,
S. R. Arridge
, and
D. T. Delpy
,
“Optical imaging in medicine: I. Experimental techniques,”
Phys. Med. Biol. , 42 825
–840
(1997). Google Scholar
Y. T. Pan
and
D. L. Farkas
,
“High-resolution imaging of living human skin with optical coherence tomography,”,”
Scanning , 21 134
–135
(1999). Google Scholar
J. Welzel
,
E. Lankenau
,
R. Birngruber
, and
R. Engelhardt
,
“Optical coherence tomography of the human skin,”,”
J. Am. Acad. Dermatol., 37 958
–963
(1997). Google Scholar
S. P. Morgan
,
M. P. Khong
, and
M. G. Somekh
,
“Effects of polarization state and scatterer concentration on optical imaging through scattering media, ,”
Appl. Opt. , 36 1560
–1565
(1997). Google Scholar
S. L. Jacques
,
J. R. Roman
, and
K. Lee
,
“Imaging superficial tissues with polarized light,”
Lasers Surg. Med. , 26 119
–129
(2000). Google Scholar
S. G. Demos
,
H. B. Radousky
, and
R. R. Alfano
,
“Deep subsurface imaging in tissues using spectral and polarization filtering,”
Opt. Express , 7 23
–28
(2000). Google Scholar
G. Yao
and
L. Wang
,
“Propagation of polarized light in turbid media: simulated animation sequences,”
Opt. Express , 7 198
–203
(2000). Google Scholar
G. D. Lewis
,
D. L. Jordan
, and
P. J. Roberts
,
“Backscattering target detection in a turbid medium by polarization discrimination,”
Appl. Opt. , 38 3937
–3944
(1999). Google Scholar
H. Dehghani
,
D. T. Delpy
, and
S. R. Arridge
,
“Photon migration in non-scattering tissue and the effects on image reconstruction,”
Phys. Med. Biol. , 44 2897
–2906
(1999). Google Scholar
H. Horinaka
,
K. Hashimoto
,
K. Wada
, and
Y. Cho
,
“Extraction of quasi-straightforward-propagating photons from diffused light transmitting through a scattering medium by polarization modulation,”
Opt. Lett. , 20 1501
–1503
(1995). Google Scholar
M. R. Ostermeyer
,
D. V. Stephens
,
L. Wang
, and
S. L. Jacques
,
“Near field polarization effects on light propagation in random media,”
OSA TOPS Biomed. Opt. Spectrosc. Diagn., 3 20
–25
(1996). Google Scholar
A. H. Hielscher
,
J. R. Mourant
, and
I. J. Bigio
, OSA TOPS Biomed. Opt. Spectros. Diagn., 3 26
–31
(1996). Google Scholar
D. Bicout
,
C. Brosseau
,
A. S. Martinez
, and
J. M. Schmitt
,
“Depolarization of multiply scattered waves by spherical diffusers: influence of the size parameter,”
Phys. Rev. E , 49 1767
–1770
(1994). Google Scholar
J. R. de Boer
,
T. E. Milner
,
M. J. C. van Gemert
, and
J. S. Nelson
,
“Determination of the depth-resolved Stokes parameters of light backscattered from turbid media by use of polarization-sensitive optical coherence tomography,”
Opt. Lett. , 22 934
–936
(1997). Google Scholar
H. Hielscher
,
A. A. Elick
,
J. R. Mourant
, and
I. J. Bigio
,
“Biomedical diagnostic with diffusely backscattered linearly and circularly polarized light,” in Biomedical Sensing, Imaging, and Tracking Technologies II, R. A. Lieberman, T. Vo-Dinh, and G. G. Vurek, Eds.,”
Proc. SPIE , 2976 298
–305
(1997). Google Scholar
A. F. Fercher
,
“Optical coherence tomography,”
J. Biomed. Opt. , 1 157
–173
(1996). Google Scholar
J. F. de Boer
,
S. M. Srinivas
,
B. H. Park
,
T. H. Pham
,
Z. Chen
,
T. E. Milner
, and
J. S. Nelson
,
“Polarization effects in optical coherence tomography of various biological tissues,”
IEEE J. Sel. Top. Quantum Electron. , 5 1200
–1203
(1999). Google Scholar
M. J. Everett
,
K. Schoenenberger
,
B. W. Colston Jr.
, and
L. B. Da Silva
,
“Birefringence characterization of biological tissue by use of optical coherence tomography,”
Opt. Lett. , 23 228
–230
(1998). Google Scholar
K. Schoenenberger
,
B. W. Colston Jr.
,
D. J. Maitland
,
L. B. Da Silva
, and
M. J. Everett
,
“Mapping of birefringence and thermal damage in tissue by use of polarization-sensitive optical coherence tomography,”
Appl. Opt. , 37 6026
–6036
(1998). Google Scholar
M. R. Hee
,
D. Huang
,
E. A. Swanson
, and
J. G. Fujimoto
,
“Polarization-sensitive low-coherence reflectometer for birefringence characterization and ranging,”
J. Opt. Soc. Am. B , 9 903
–908
(1992). Google Scholar
A. G. Ushenko
,
“Laser biospeckles fields vector structure and polarization diagnostics of skin collagen structure,”
Laser Phys. , 10
(5), 1143
–1149
(2000). Google Scholar
O. V. Angelsky
,
D. M. Burkovets
,
V. P. Pishak
,
Y. A. Ushenko
, and
O. V. Pishak
,
“Polarization-correlation investigations of biotissue multifractal structures and their pathological changes diagnostics,”
Laser Phys. , 10
(5), 1136
–1142
(2000). Google Scholar
O. V. Angelsky
,
A. G. Ushenko
, and
D. M. Burkovets
,
“Polarization-correlation analysis of anisotropic structures in bone tissue for the diagnostics of pathological changes,”
Opt. Spektrosk. , 90
(3), 458
–464
(2001). Google Scholar
S. Bartel
and
A. H. Hielscher
,
“Monte Carlo simulations of the diffuse backscattering Mueller matrix for highly scattering media,”
Appl. Opt. , 39 1580
–1588
(2000). Google Scholar
G. Yao
and
L.-H. Wang
,
“Propagation of polarized light in turbid media: simulated animation sequences,”
Opt. Express , 7 198
–203
(2000). Google Scholar
A. Ambirajan
and
D. C. Look
,
“A backward Monte Carlo study of the multiple scattering of a polarized light beam,”
J. Quantum Spectrosc. Radiat. Trans., 58 171
–192
(1997). Google Scholar
X. Wang
,
L.-H. Wang
,
C.-W. Sun
, and
C. C. Yang
,
“Polarized light propagation through scattering media: time-resolved Monte Carlo simulations and experiments,”
J. Biomed. Opt. , 8
(4), 608
–617
(2003). Google Scholar
X. Wang
,
G. Yao
, and
L.-H. Wang
,
“Monte Carlo model and single-scattering approximation of the propagation of polarized light in turbid media containing glucose,”
Appl. Opt. , 41
(4), 792
–801
(2002). Google Scholar
X. Wang
and
L.-H. Wang
,
“Propagation of polarized light in birefringent turbid media: a Monte Carlo study,”
J. Biomed. Opt. , 7
(3), 279
–290
(2002). Google Scholar
G. Yao
and
L.-H. Wang
,
“Two dimensional depth resolved Mueller matrix measurement in biological tissue with optical coherence tomography,”,”
Opt. Lett. , 24 537
–539
(1999). Google Scholar
D. Cameron
,
M. J. Rakovic
,
M. Mehrubeoglu
,
G. Kattawar
,
S. Rastegar
,
L. V. Wang
, and
G. L. Cote
,
“Measurement and calculation of the two-dimensional backscattering Mueller matrix of a turbid medium,”
Opt. Lett. , 23 485
–487
(1998). Google Scholar
J. F. de Boer
,
T. E. Milner
, and
J. S. Nelson
,
“Determination of the depth-resolved Stokes parameters of light backscattered from turbid media by use of polarization-sensitive optical coherence tomography,”
Opt. Lett. , 24 300
–302
(1999). Google Scholar
S. Jiao
,
G. Yao
, and
L. V. Wang
,
“Depth-resolved two-dimensional Stokes vectors of backscattered light and Mueller matrices of biological tissue measured with optical coherence tomography,”
Appl. Opt. , 39 6318
–6324
(2000). Google Scholar
S. Jiao
and
L.-H. Wang
,
“Two-dimensional depth-resolved Mueller matrix of biological tissue measured with double-beam polarization-sensitive optical coherence tomography,”
Opt. Lett. , 27
(2), 101
–103
(2002). Google Scholar
S. Jiao
,
W. Yu
,
G. Stoica
, and
L.-H. Wang
,
“Contrast mechanisms in polarization-sensitive Mueller-matrix optical coherence tomography and application in burn imaging,”
Appl. Opt. , 42
(25), 5191
–5197
(2003). Google Scholar
S. Jiao
,
W. Yu
,
G. Stoica
, and
L.-H. Wang
,
“Optical-fiber-based Mueller optical coherence tomography,”
Opt. Lett. , 28
(14), 1206
–1208
(2003). Google Scholar
E. Okada
,
M. Firbank
,
M. Schweiger
, et al.,
“Theoretical and experimental investigation of near-infrared light propagation in a model of the adult head,”
Appl. Opt. , 36
(1), 21
–31
(1997). Google Scholar
A. Kienle
and
R. Hibst
,
“New optical wavelength for treatment of portwine stains,”
Phys. Med. Biol. , 40 1559
–1576
(1995). Google Scholar
S. L. Jacques
,
C. A. Alter
, and
S. A. Prahl
,
“Angular dependence of HeNe laser light scattering by human dermis,”
Lasers Life Sci. , 1 309
–333
(1987). Google Scholar
R. Marchesini
,
A. Bertoni
,
S. Andreola
, et al.,
“Extinction and absorption coefficients and scattering phase functions of human tissues in vitro,”
Appl. Opt. , 28 2318
–2324
(1989). Google Scholar
|