Home | About | Journals | Submit | Contact Us | Français |

**|**HHS Author Manuscripts**|**PMC2749477

Formats

Article sections

- Abstract
- 1. Introduction
- 2. Background and Theory
- 3. Methods
- 4. Results
- 5. Discussion
- 6. Conclusion
- References and links

Authors

Related links

Opt Express. Author manuscript; available in PMC 2010 August 3.

Published in final edited form as:

PMCID: PMC2749477

NIHMSID: NIHMS143207

The publisher's final edited version of this article is available at Opt Express

See other articles in PMC that cite the published article.

Complex polarization ratio (CPR) in materials with birefringence and biattenuance is shown as a logarithmic spiral in the complex plane. A multi-state Levenberg-Marquardt nonlinear fitting algorithm using the CPR trajectory collected by polarization sensitive optical coherence tomography (PS-OCT) was developed to determine polarization properties of an anisotropic scattering medium. The Levenberg-Marquardt nonlinear fitting algorithm using the CPR trajectory is verified using simulated PS-OCT data with speckle noise. Birefringence and biattenuance of a birefringent film, *ex-vivo* rodent tail tendon and *in-vivo* primate retinal nerve fiber layer were determined using measured CPR trajectories and the Levenberg-Marquardt nonlinear fitting algorithm.

Polarization-sensitive optical coherence tomography (PS-OCT) is a non-invasive imaging modality that provides depth-resolved polarimetric information with high-resolution in anisotropic tissues [1,2]. Measurement of depth-resolved polarimetric information by PS-OCT is widely investigated for diagnosis of various pathological conditions or trauma by observing the variation of polarimetric properties between normal and abnormal tissues [3–8].

In general PS-OCT, horizontal $(\phantom{\rule{thinmathspace}{0ex}}\phantom{\rule{thinmathspace}{0ex}}{\stackrel{~}{\Gamma}}_{h}\left(\mathrm{z}\right)\phantom{\rule{thinmathspace}{0ex}}\phantom{\rule{thinmathspace}{0ex}})$ and vertical $(\phantom{\rule{thinmathspace}{0ex}}\phantom{\rule{thinmathspace}{0ex}}{\stackrel{~}{\Gamma}}_{v}\left(\mathrm{z}\right)\phantom{\rule{thinmathspace}{0ex}}\phantom{\rule{thinmathspace}{0ex}})$ interference fringe magnitudes and relative phase $(\phantom{\rule{thinmathspace}{0ex}}\mathrm{{\stackrel{~}{\Gamma}}_{vh}\left(\mathrm{z}\right)\phantom{\rule{thinmathspace}{0ex}})}$ as a function of depth (*z*) are recorded by controlling polarization elements in source, sample and reference paths of time- or frequency-domain OCT instrumentation [9,10]. Diverse analytic techniques have been used to extract polarimetric information such as depth-resolved phase retardation (*δ(z)*), birefringence (*Δn*), amplitude attenuation (*ε(z)*) (or diattenuation (*d(z)*)), biattenuance (*Δχ*), and optic axis orientation (*α*) from $\phantom{\rule{thinmathspace}{0ex}}{\stackrel{~}{\Gamma}}_{h}\left(\mathrm{z}\right)\phantom{\rule{thinmathspace}{0ex}}$, $\phantom{\rule{thinmathspace}{0ex}}{\stackrel{~}{\Gamma}}_{v}\left(\mathrm{z}\right)\phantom{\rule{thinmathspace}{0ex}}$ and $\mathrm{{\stackrel{~}{\Gamma}}_{vh}\left(\mathrm{z}\right)}$ [11–27]. In early PS-OCT studies, the phase retardation was computed from the arc tangent of the ratio of the horizontal and vertical fringe magnitudes $(\phantom{\rule{thinmathspace}{0ex}}\delta \left(\mathrm{z}\right)=\text{arctan}(\phantom{\rule{thinmathspace}{0ex}}{\stackrel{~}{\Gamma}}_{v}\left(\mathrm{z}\right)\phantom{\rule{thinmathspace}{0ex}}\phantom{\rule{thinmathspace}{0ex}}/\phantom{\rule{thinmathspace}{0ex}}\phantom{\rule{thinmathspace}{0ex}}{\stackrel{~}{\Gamma}}_{h}\left(\mathrm{z}\right)\phantom{\rule{thinmathspace}{0ex}})\phantom{\rule{thinmathspace}{0ex}})$ [1,11]. Although this technique was simple and appropriate for a coarse two-dimensional image, polarimetric values in a target region (e.g. retinal nerve fiber layer in retina) were not accurately computed due to background noise including speckle [9].

A Mueller matrix formalism was applied to time- and frequency-domain PS-OCT to obtain the polarimetric properties of biological tissue [12–14]. A two-dimensional depth-resolved 4 × 4 Mueller matrix of the tissue sample was measured by open-air and fiber-based PS-OCT systems. Birefringence, diattenuation and optic axis orientation were extracted from the measured Mueller matrix. Although the Mueller matrix can provide the complete polarimetric transformation of the tissue specimen, each element of the Mueller matrix is not easily interpreted.

A Jones matrix formalism was applied to analyze polarization state of light backscattered from tissue and recorded by PS-OCT [15–19]. Light propagation into many components (including tissue specimen) was described by products of Jones matrices in time-domain fiber-based PS-OCT. Birefringence, diattenuation and optic axis orientation of the samples were determined by calculation of the Jones matrices [15]. Similar to Mueller matrix images using PS-OCT, two-dimensional depth-resolved Jones matrix images were demonstrated and local polarization properties such as phase retardation, diattenuation, and optic axis orientation were computed in open-air [16,17] and fiber-based [18,19] Fourier-domain PS-OCT.

Stokes parameters are also utilized to analyze polarimetric properties of biological tissue using PS-OCT [20–27]. In early studies using Stokes parameters, two-dimensional images of depth-resolved Stokes parameters were generated and compared to intensity images for identifying tissue characteristics [20,21]. As the Stokes parameters are represented geometrically as a three-dimensional vector (Stokes vector), polarization properties of tissue in PS-OCT can be visually interpreted compared with other polarization analysis techniques such as Jones vectors. Depth-resolved Stokes vectors on the Poincaré sphere allow visualization of the polarization state of light backscattered by a tissue specimen. Birefringence and optic axis orientation were obtained by vector calculation using the depth-resolved Stokes vectors in a fiber-based PS-OCT instrument [22,23]. Trajectory of the depth-resolved Stokes vectors on the Poincaré sphere was theoretically and experimentally investigated corresponding to light propagation in anisotropic tissues. Numerical expressions of the trajectory and associated differential geometry were derived for materials with arbitrary birefringence, biattenuance and optic axis orientation [24]. A Levenberg-Marquardt nonlinear fitting algorithm with multi-incident polarization states of light were applied to analyze the depth-resolved Stokes parameters of backscattered light from tissue specimens recorded by an open-air PS-OCT instrument. Tissue birefringence, biattenuance, and optic axis orientation were determined by estimating Stokes parameters from speckle-noise corrupted PS-OCT data. Multiple incident polarization states of light were used to suppress noise and increase polarimetric signal to noise ratio (PSNR) of PS-OCT data [25–27].

Complex valued analytic signals are utilized to solve a variety of problems arising in science and engineering. Even though real numbers are natural for representing recorded data, complex numbers provide a useful approach to analyze many engineering problems. In control theory, systems are transformed from time-domain to frequency-domain, and vice versa using Laplace or Z-transforms. Characteristics of systems are analyzed by poles and zeros using signals represented as complex numbers in the complex plane [28]. In other fields such as fluid dynamics, quantum mechanics, relativity and applied mathematics [29], complex numbers are routinely employed to represent real phenomena.

We demonstrate a new approach which analyzes the trajectory of the complex polarization ratio (CPR) in the complex plane of PS-OCT data to determine polarimetric properties of biological tissues. The technique using CPR combines the advantages of matrix formalisms (Jones and Mueller matrices) and Stokes parameters. Similar to Jones and Mueller matrices, CPR can mathematically express the polarization state of light backscattered from tissue as a single complex number. Computations relating to the polarization state of light are simplified, however, because the CPR is a single number. Similar to Stokes vectors on the Poincaré sphere, CPR may be displayed geometrically to visualize the polarization state of light backscattered from tissue on a complex plane. Moreover, two-dimensional display of CPR on the flat complex plane is more convenient than three-dimensional display of Stokes vectors on the Poincaré sphere.

We utilized a Levenberg-Marquardt nonlinear fitting algorithm to determine the polarization properties of a tissue specimen from CPR trajectories in the complex plane. The algorithm was verified using CPR trajectories of simulated PS-OCT data corrupted with polarimetric speckle noise. In addition, the CPR algorithm was applied to PS-OCT data recorded from a birefringent film, *ex-vivo* rodent tail tendon and *in-vivo* primate retinal nerve fiber layer (RNFL).

Ellipsometry is an optical technique for determining polarization properties of a material by observing the change of polarization state of light reflected from a sample surface. Jones vector and complex polarization ratio (CPR) are used to describe purely polarized light. However, Stokes vector and coherence matrix formalisms are the most general representations of the polarization state of incident light. Both representations describe both polarized and unpolarized components of light, whereas Jones vector and CPR are used to describe purely polarized components [30,31].

The CPR (* C_{yx}*) is defined by the ratio of Jones vector components (

$${\mathit{C}}_{yx}=\frac{{\mathit{E}}_{y}}{{\mathit{E}}_{x}}=\frac{{\mathit{E}}_{y}{e}^{j{\theta}_{y}}{\mathit{E}}_{x}{e}^{j{\theta}_{x}}=\frac{{\mathit{E}}_{y}{\mathit{E}}_{x}{e}^{j({\theta}_{y}-{\theta}_{x})}}{}}{}$$

(1)

where | * C_{yx}* |=|

An arbitrary CPR representing a purely polarized state of light can be assigned to a point on the Cartesian complex plane. All purely polarized states of light can be transformed from the CPR representation on the Cartesian complex plane (more simply the complex plane). For example, if linearly horizontal (*h*) and vertical (*v*) polarization states of light are used as an orthonormal basis set, the origin (* C_{vh}* = 0) and the point-at-infinity (

Although a horizontal and vertical linear polarization basis set (*h, v*) is commonly used in CPR polarization analysis, any two fixed elliptic states (*a, b*) can be selected as a basis set. Position of CPR in the complex plane depends on the selected basis set, and numerous displays of CPR on the complex plane are possible. An algebraic expression for transformation of the CPR due to a change in basis set is easily derived from Jones vector calculus. An arbitrary state of polarized light expressed by Jones vector (* E_{xy}*) with the basis set (

$$\left[\begin{array}{c}\hfill {E}_{x}\hfill \\ \hfill {E}_{y}\hfill \end{array}\right]=\left[\begin{array}{cc}\hfill {f}_{11}\hfill & \hfill {f}_{12}\hfill \\ \hfill {f}_{21}\hfill & \hfill {f}_{22}\hfill \end{array}\right]\left[\begin{array}{c}\hfill {E}_{a}\hfill \\ \hfill {E}_{b}\hfill \end{array}\right],$$

(2)

the transformation of CPR between basis sets (*x, y*) and (*a, b*) is computed by taking the ratios in both sides of Eq. (2), and gives

$${\mathit{C}}_{yx}=\frac{{f}_{22}{\mathit{C}}_{ba}+{f}_{21}}{{f}_{12}{\mathit{C}}_{ba}+{f}_{11}}.$$

(3)

The ratio of two linear functions (Eq. (3)) is a Möbious (or linear fractional) transformation. Lines and circles in the basis states (*x, y*) are mapped to lines and circles in the basis set (*a, b*) (and vice versa) by the characteristics of the Möbious transformation [30,31].

We consider a birefringent tissue with a polarization state corresponding to an optic axis given by CPR *a* and an orthogonal state *b*. The orthonormal basis pair is represented by (*a, b*). We consider this basis set and examine the CPR trajectory corresponding to light propagation in the birefringent tissue. After forward-scattered light propagates a distance (*z*) through tissue, the CPR with an arbitrary polarization basis set (*a, b*) is given by

$${\mathit{C}}_{ba}\left(\mathrm{z}\right)={E}_{b}\left(\mathrm{z}\right)\phantom{\rule{thinmathspace}{0ex}}/\phantom{\rule{thinmathspace}{0ex}}{E}_{a}\left(\mathrm{z}\right)=\text{exp}(-\epsilon \left(\mathrm{z}\right))\text{exp}\left(i\delta \left(\mathrm{z}\right)\right)$$

(4)

where *δ(z)* and *ε(z)* are the depth-resolved phase retardation and amplitude attenuation, respectively [24]. When the optic axis corresponds to either states *a* or *b*, a depth-resolved trajectory of the CPR in the (*a, b*) basis complex plane is generated by rotation *δ(z)* and attenuation *ε(z)* (Fig. 2). Considering only *δ(z)*, the trajectory of the CPR uniformly rotates around the origin and forms a circular arc (Fig. 2(a)). Diattenuation (or biattenuance [27]) leads to *ε(z)* exponentially collapsing the trajectory for increasing depths (*z*). Therefore, trajectory in the complex plane with both *δ(z)* and *ε(z)* becomes a logarithmic spiral converging toward the origin in the (*a, b*) basis set (Fig. 2(b)).

Multi-state, noise-free PS-OCT data were generated to verify polarization analysis using the CPR representation. The *M* incident polarization states were uniformly distributed on a meridian of the Poincaré sphere. A Jones matrix formalism was used to transform each incident polarization state due to light propagation through the material [14,15]. The transform gave noise-free PS-OCT data as a function of depth (*z*) for the *M* polarization states. For practical simulations, speckle-noise corrupted PS-OCT data were also generated by adding random speckle noise onto noise-free PS-OCT data. Random speckle noise had a uniformly distributed phase variation between 0 and 2π radians and Gaussian distributed magnitude variation with zero mean and standard deviation (*σ*).

An open-air time-domain PS-OCT instrument with a mode-locked Ti:Al_{2}O_{3} laser source (*λ*_{0} = 830 nm, *λ*_{FWHM} = 55 nm) was employed to record PS-OCT data of specimens including birefringent film, *ex-vivo* rat tail tendon and *in-vivo* primate retinal nerve fiber layer (RNFL). PS-OCT data of specimens in the *M* polarization states were acquired by controlling the polarization state of light incident on tissue specimens using a liquid-crystal variable retarder positioned in the sample arm of the interferometer. Details of the PS-OCT system were described previously [9].

PS-OCT data (horizontal $(\phantom{\rule{thinmathspace}{0ex}}\phantom{\rule{thinmathspace}{0ex}}{\stackrel{~}{\Gamma}}_{h}\left(\mathrm{z}\right)\phantom{\rule{thinmathspace}{0ex}}\phantom{\rule{thinmathspace}{0ex}})$ and vertical $(\phantom{\rule{thinmathspace}{0ex}}\phantom{\rule{thinmathspace}{0ex}}{\stackrel{~}{\Gamma}}_{v}\left(\mathrm{z}\right)\phantom{\rule{thinmathspace}{0ex}}\phantom{\rule{thinmathspace}{0ex}})$ interference fringe magnitudes and relative phase $(\phantom{\rule{thinmathspace}{0ex}}\mathrm{{\stackrel{~}{\Gamma}}_{vh}\left(\mathrm{z}\right)\phantom{\rule{thinmathspace}{0ex}})}$) were converted to CPRs (* C_{vh}* (z)) as a function of depth (

A basis transformation between (*h, v*) and (*x, y*) is given by Eq. (3). If the (*x, y*) basis set from the optic axis of specimen are determined, * C_{vh}* (z) can be transformed to the CPR (

A Levenberg-Marquardt nonlinear fitting algorithm with single incident polarization state was developed for determining an unbiased estimate of double-pass phase retardation (2*δ*), double-pass amplitude attenuation (2*ε*), CPR of optic axis (* C_{vh _ oa}*), noise-free CPRs at the surface of specimen in (

$${\mathit{R}}_{\mathit{o}}=\underset{z}{\Sigma}{{\mathit{c}}_{yx}\left(\mathrm{z}\right)-{\mathcal{T}}_{yx}\left[{\mathit{C}}_{vh}\{\mathrm{z};2\delta ,2\epsilon ,{\mathit{C}}_{vh\_oa},{\mathit{C}}_{vh}\left(0\right)\}\right]2}^{}$$

(5)

where * T_{yx}*[

A multi-state residual function (* R_{M}*) was computed as an algebraic sum of single-state residual function (

$${\mathit{R}}_{\mathit{M}}\phantom{\rule{thickmathspace}{0ex}}=\underset{\mathit{m}\phantom{\rule{thickmathspace}{0ex}}=\phantom{\rule{thickmathspace}{0ex}}1}{\overset{\mathit{M}}{\Sigma}}{\mathit{R}}_{\mathit{o}}[{\mathit{c}}_{yx\left(m\right)}\left(\mathrm{z}\right);2\delta ,2\epsilon ,{\mathit{C}}_{vh\_oa},{\mathit{C}}_{vh\left(m\right)}\left(0\right)]$$

(6)

More accurate 2*δ*, 2*ε*, * C_{vh _ oa}* , and noise-free CPRs at the surface of the specimen in (

Trajectory of CPRs using single-state simulated noise-free PS-OCT data (* C_{vh}* (z)) is displayed in (

Trajectories of CPRs in the (*h, v*) and (*x, y*) basis complex planes. (a) Trajectory with double-pass phase retardation (2*δ(z)* = 360°) and an optic axis (**c**_{vh _ oa} = 0.3exp(*j*45°), black dot) in the (*h, v*) basis complex plane (b) Trajectory **...**

Multi-state trajectories of simulated noise-free PS-OCT data (**C**_{vh(M =6)} (z)) with six incident polarization states (*M* = 6) are plotted in (*h, v*) basis complex plane (Fig. 4(a)). Each trajectory is generated with 2*δ(z)* = 60°, 2*ε(z)* = 6.0° and optic axis (* C_{vh _ oa}* = 0.3exp(

Multi-state trajectories (*M* = 6) of speckle-noise corrupted CPRs (**C**_{vh(M =6)} (z)) in (*h, v*) basis complex plane are applied to verify the ability of a multi-state Levenberg-Marquardt nonlinear algorithm to determine *δ* and *ε*. Polarization values (2*δ(z)* = 60°, 2*ε(z)* = 6.0° and * C_{vh _ oa}* = 0.3exp(

Phase retardation (*δ _{film}(z)*) of a turbid birefringent film (New Focus, #5842) are determined by the multi-state (M = 6) Levenberg-Marquardt nonlinear fitting algorithm using the CPRs. The thickness was 80μm (

*Ex-vivo* rodent tail tendon is an excellent specimen to study both phase retardation (*δ _{tendon}*

*In-vivo* primate RNFL was imaged to obtain the phase retardation (*δ _{RNFL}(z)*) in a thick region (1mm inferior to the center of the optic nerve head, thickness

We have demonstrated the use of CPR to analyze polarized light backscattered from a medium and measured with PS-OCT instrumentation. CPR provides a mathematically simple framework to easily visualize polarization states in the complex plane. Application of CPR to determine the polarimetric properties of a tissue specimen using a Levenberg-Marquardt nonlinear fitting algorithm was evaluated. Three main considerations such as complexity of the Levenberg-Marquardt nonlinear fitting algorithm, processing-time efficiency and accuracy of fitted polarization properties are investigated in greater detail by comparing Levenberg-Marquardt nonlinear fitting algorithms using CPR to Stokes vectors.

Stokes vectors represented by three real numbers are computed from horizontal $(\phantom{\rule{thinmathspace}{0ex}}\phantom{\rule{thinmathspace}{0ex}}{\stackrel{~}{\Gamma}}_{h}\left(\mathrm{z}\right)\phantom{\rule{thinmathspace}{0ex}}\phantom{\rule{thinmathspace}{0ex}})$ and vertical $(\phantom{\rule{thinmathspace}{0ex}}\phantom{\rule{thinmathspace}{0ex}}{\stackrel{~}{\Gamma}}_{v}\left(\mathrm{z}\right)\phantom{\rule{thinmathspace}{0ex}}\phantom{\rule{thinmathspace}{0ex}})$ interference fringe magnitudes and relative phase $(\phantom{\rule{thinmathspace}{0ex}}\mathrm{{\stackrel{~}{\Gamma}}_{vh}\left(\mathrm{z}\right)\phantom{\rule{thinmathspace}{0ex}})}$ as a function of depth (*z*). This calculation requires eight real multiplications and two real additions. In contrast, CPR represented by a complex number requires only three multiplications for the same computation. When *M* multiple polarization states with *N* depth samples in each state are considered to determine polarization properties of a specimen, 8 × *M* × *N* real multiplications and 2 × *M* × *N* real additions for the Stokes vectors, and 3 × *M* × *N* real multiplications for CPR are required. Therefore, CPR is a more efficient representation than Stokes vectors for the conversion from the recorded PS-OCT data ($\phantom{\rule{thinmathspace}{0ex}}{\stackrel{~}{\Gamma}}_{h}\left(\mathrm{z}\right)\phantom{\rule{thinmathspace}{0ex}},\phantom{\rule{thickmathspace}{0ex}}\phantom{\rule{thinmathspace}{0ex}}{\stackrel{~}{\Gamma}}_{v}\left(\mathrm{z}\right)\phantom{\rule{thinmathspace}{0ex}}$ and $\mathrm{{\stackrel{~}{\Gamma}}_{vh}\left(\mathrm{z}\right)}$). In the multi-state Levenberg-Marquardt nonlinear fitting algorithms using polarization data, a trajectory of noise-free polarization data is generated using model parameters to minimize a residual function. Complexity of the residual function is highly correlated to generation of a noise-free trajectory. As the Stokes vectors are applied to the nonlinear fitting algorithms, the trajectory is three-dimensional and requires use of intricate geometrical expressions on the Poincaré sphere with approximately 22 × *M* × *N* real multiplications and 10 × *M* × *N* real additions [24]. Although the conversion of polarization basis set with about 16 × *M* × *N* real (or 4 × *M* × *N* complex) multiplications and 8 × *M* × *N* real (or 4 × *M* × *N* complex) additions slightly complicates generation of CPR trajectories, total number of operations with 20 × *M* × *N* real (or 5 × *M* × *N* complex) multiplications and 8 × *M* × *N* real (or 4 × *M* × *N* complex) additions are required. For example, as six polarization states (*M* = 6) and 1000 depth samples (*N* = 1000) are used, a reduction of 24,000 computational operations are realized using CPRs that can provide to a substantial reduction in processing time as discussed in the next section.

Processing times of the Levenberg-Marquardt nonlinear fitting algorithms were evaluated using simulated CPRs and Stokes vectors. 100 independent estimates using equivalent polarization conditions are measured and averaged. The nonlinear fitting algorithms used a Microsoft Windows XP operating system and hence measured processing times included Windows OS calls. Assuming that the Windows operation time is identical for both algorithms, relative efficiency can be investigated by computing the relative processing times between two algorithms. Figure 8(a) shows the relative processing time which is computed using the measured processing time divided by a maximum processing time in both algorithms. Phase retardation (*δ(z)* = 10°, 30°, 50°, 70° and 90°) and high speckle noise with standard deviation (*σ* = 5°) were applied to generate the speckle noise-corrupted CPRs and Stokes vectors. For the Stokes vector calculation, as the retardation increases a monotonic reduction of the relative processing time is observed. However, for CPR as the retardation increases a more rapid reduction of the relative processing time is observed. Difference of the relative processing time apparently increases, and an over 20% time difference is observed when the retardation is 90° (Fig. 8(b)). Therefore, the nonlinear algorithm using CPRs is faster by approximately 25% than the algorithm using Stokes vectors. Processing times with respect to different speckle noise (*σ* = 1°, 2°, 3°, 4° and 5°) at a phase retardation (*δ(z)* = 30°) are also studied and displayed in Fig. 8(c). The processing time in the Stokes vectors is less affected than that in CPRs by magnitude of speckle noise, and a 10% time difference is observed when the speckle nose is 5° (Fig. 8(d)). As speckle noise is reduced, difference between processing times increases. In other words, determination of polarimetric properties is increasingly faster as the speckle noise is reduced in the nonlinear fitting algorithm using CPRs compared to Stokes vectors.

Accuracy of Levenberg-Marquardt nonlinear algorithms using CPRs and Stokes vectors was compared by a statistical analysis using simulated PS-OCT data. In order to know distribution characteristics of the unknown variables, we computed confidence intervals in both nonlinear algorithms. Based on the estimated polarization properties determined by the Levenberg-Marquardt nonlinear fitting algorithm, we calculated the estimate (*s ^{2}*) of the variance of the residual function with

$$C{I}_{i}=\pm t\sqrt{{s}^{2}{P}_{ii}}$$

(7)

where *t* is a value of the t-distribution, and *P _{ii}* is the

The confidence intervals of two algorithms were tested by varying depth-resolved phase retardation and speckle noise using simulated PS-OCT data. 100 independent estimate sets were obtained and averaged under the same condition. Figure 9(a) shows 95% confidence intervals computed by changing phase retardation (*δ(z)* = 10°, 30°, 50°, 70° and 90°) at a fixed optic axis (* c_{vh _ oa}* = 0.3exp(

CPR is a mathematical expression to represent polarization state of light as a single complex number. Complex number computation of CPRs and two-dimensional visualization in a complex plane are competitive advantages compared to methods to analyze PS-OCT data using Jones, Mueller matrices and Stokes vectors. A Levenberg-Marquardt nonlinear fitting algorithm using CPR was developed to determine polarization properties including phase retardation (or birefringence), relative amplitude attenuation (or biattenuance) and optic axis orientation of anisotropic tissues from PS-OCT data. After the algorithm was tested by simulated PS-OCT data, polarization properties of several anisotropic specimens including birefringent film, *ex-vivo* rat tail tendon and *in-vivo* primate RNFL are determined. The Levenberg-Marquardt nonlinear fitting algorithm using CPR has less complexity, faster processing time than that using Stokes vectors [9,25]. Although uncertainty of estimated polarization properties are slightly greater when using CPR, the benefits provides by easy visualization, reduced complexity and fast processing time are important practical advantages using CPR.

The authors gratefully acknowledge the Texas Higher Education Coordinating Board Advanced Technology Program and the National Institutes of Health for financial support of this research (R01 EY 016462-04).

1. Hee MR, Huang D, Swanson EA, Fujimoto JG. Polarization-sensitive low-coherence reflectometer for birefringence characterization and ranging. J. Opt. Soc. Am. B. 1992;9(6):903–908.

2. de Boer JF, Milner TE. Review of polarization sensitive optical coherence tomography and Stokes vector determination. J. Biomed. Opt. 2002;7(3):359–371. [PubMed]

3. Hitzenberger CK, Götzinger E, Pircher M. Birefringence properties of the human cornea measured with polarization sensitive optical coherence tomography. Bull. Soc. Belge Ophtalmol. 2006;302(302):153–168. [PubMed]

4. Cense B, Chen TC, Park BH, Pierce MC, de Boer JF. Thickness and birefringence of healthy retinal nerve fiber layer tissue measured with polarization-sensitive optical coherence tomography. Invest. Ophthalmol. Vis. Sci. 2004;45(8):2606–2612. [PubMed]

5. Pierce MC, Sheridan RL, Hyle Park B, Cense B, de Boer JF. Collagen denaturation can be quantified in burned human skin using polarization-sensitive optical coherence tomography. Burns. 2004;30(6):511–517. [PubMed]

6. Patel NA, Zoeller J, Stamper DL, Fujimoto JG, Brezinski ME. Monitoring osteoarthritis in the rat model using optical coherence tomography. IEEE Trans. Med. Imaging. 2005;24(2):155–159. [PubMed]

7. Giattina SD, Courtney BK, Herz PR, Harman M, Shortkroff S, Stamper DL, Liu B, Fujimoto JG, Brezinski ME. Assessment of coronary plaque collagen with polarization sensitive optical coherence tomography (PS-OCT) Int. J. Cardiol. 2006;107(3):400–409. [PubMed]

8. Jones RS, Darling CL, Featherstone JDB, Fried D. Remineralization of in vitro dental caries assessed with polarization-sensitive optical coherence tomography. J. Biomed. Opt. 2006;11(1):014016. [PubMed]

9. Kemp NJ, Park J, Zaatari HN, Rylander HG, Milner TE. High-sensitivity determination of birefringence in turbid media with enhanced polarization-sensitive optical coherence tomography. J. Opt. Soc. Am. A. 2005;22(3):552–560. [PubMed]

10. Götzinger E, Pircher M, Hitzenberger CK. High speed spectral domain polarization sensitive optical coherence tomography of the human retina. Opt. Express. 2005. pp. 10217–10229. http://www.opticsinfobase.org/abstract.cfm?URI=oe-13-25-10217. [PMC free article] [PubMed]

11. Hitzenberger CK, Goetzinger E, Sticker M, Pircher M, Fercher AF. Measurement and imaging of birefringence and optic axis orientation by phase resolved polarization sensitive optical coherence tomography. Opt. Express. 2001. pp. 780–790. http://www.opticsinfobase.org/abstract.cfm?URI=oe-9-13-780. [PubMed]

12. Jiao SL, Wang LV. Two-dimensional depth-resolved Mueller matrix of biological tissue measured with double-beam polarization-sensitive optical coherence tomography. Opt. Lett. 2002;27(2):101–103. [PubMed]

13. Jiao SL, Todorović M, Stoica G, Wang LV. Fiber-based polarization-sensitive Mueller matrix optical coherence tomography with continuous source polarization modulation. Appl. Opt. 2005;44(26):5463–5467. [PubMed]

14. Yasuno Y, Makita S, Sutoh Y, Itoh M, Yatagai T. Birefringence imaging of human skin by polarization-sensitive spectral interferometric optical coherence tomography. Opt. Lett. 2002;27(20):1803–1805. [PubMed]

15. Park BH, Pierce MC, Cense B, de Boer JF. Jones matrix analysis for a polarization-sensitive optical coherence tomography system using fiber-optic components. Opt. Lett. 2004;29(21):2512–2514. [PMC free article] [PubMed]

16. Yasuno Y, Makita S, Endo T, Itoh M, Yatagai T, Takahashi M, Katada C, Mutoh M. Polarization-sensitive complex Fourier domain optical coherence tomography for Jones matrix imaging of biological samples. Appl. Phys. Lett. 2004;85(15):3023–3025.

17. Makita S, Yasuno Y, Endo T, Itoh M, Yatagai T. Polarization contrast imaging of biological tissues by polarization-sensitive Fourier-domain optical coherence tomography. Appl. Opt. 2006;45(6):1142–1147. [PubMed]

18. Yamanari M, Makita S, Madjarova VD, Yatagai T, Yasuno Y. Fiber-based polarization-sensitive Fourier domain optical coherence tomography using B-scan-oriented polarization modulation method. Opt. Express. 2006. pp. 6502–6515. http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-14-14-6502. [PubMed]

19. Yamanari M, Makita S, Yasuno Y. Polarization-sensitive swept-source optical coherence tomography with continuous source polarization modulation. Opt. Express. 2008. pp. 5892–5906. http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-16-8-5892. [PubMed]

20. de Boer JF, Milner TE, Nelson JS. Determination of the depth-resolved Stokes parameters of light backscattered from turbid media by use of polarization-sensitive optical coherence tomography. Opt. Lett. 1999;24(5):300–302. [PubMed]

21. Ducros MG, Marsack JD, Rylander HG, III, Thomsen SL, Milner TE. Primate retina imaging with polarization-sensitive optical coherence tomography. J. Opt. Soc. Am. A. 2001;18(12):2945–2956. [PubMed]

22. Park BH, Saxer C, Srinivas SM, Nelson JS, de Boer JF. In vivo burn depth determination by high-speed fiber-based polarization sensitive optical coherence tomography. J. Biomed. Opt. 2001;6(4):474–479. [PubMed]

23. Park BH, Pierce MC, Cense B, de Boer JF. Optic axis determination accuracy for fiber-based polarization-sensitive optical coherence tomography. Opt. Lett. 2005;30(19):2587–2589. [PubMed]

24. Park J, Kemp NJ, Zaatari HN, Rylander HG, III, Milner TE. Differential geometry of normalized Stokes vector trajectories in anisotropic media. J. Opt. Soc. Am. A. 2006;23(3):679–690. [PubMed]

25. Rylander HG, 3rd, Kemp NJ, Park JS, Zaatari HN, Milner TE. Birefringence of the primate retinal nerve fiber layer. Exp. Eye Res. 2005;81(1):81–89. [PubMed]

26. Kemp NJ, Zaatari HN, Park J, Rylander HG, III, Milner TE. Depth-resolved optic axis orientation in multiple layered anisotropic tissues measured with enhanced polarization-sensitive optical coherence tomography (EPS-OCT) Opt. Express. 2005. pp. 4507–4518. http://www.opticsinfobase.org/abstract.cfm?URI=oe-13-12-4507. [PubMed]

27. Kemp NJ, Zaatari HN, Park J, Rylander HG, III, Milner TE. Form-biattenuance in fibrous tissues measured with polarization-sensitive optical coherence tomography (PS-OCT) Opt. Express. 2005. pp. 4611–4628. http://www.opticsinfobase.org/abstract.cfm?URI=oe-13-12-4611. [PubMed]

28. Lathi BP. Linear systems and signals. Oxford University Press; New York: 2005.

29. Paliouras JD. Complex variables for scientists and engineers. Macmillan; New York: 1975.

30. Azzam RMA, Bashara NM. Ellipsometry and polarized light. North-Holland Pub. Co., Amsterdam; New York: 1977.

31. Brosseau C. Fundamentals of polarized light: a statistical optics approach. Wiley; New York: 1998.

32. Gallant AR. Nonlinear-Regression. Am. Stat. 1975;29(2):73–81.

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |