Translator Disclaimer
1 June 2011 Investigation of formation and interrelations of polarization singular structure and Mueller-matrix images of biological tissues and diagnostics of their cancer changes
Author Affiliations +
We investigate the processes of laser radiation transformation by biological crystal networks using the singular optics techniques. The obtained results show a distinct correlation between the points of "characteristic" values of coordinate distributions of Mueller matrix (Mik = 0, ±1) elements and polarization singularities (L- and C-points) in laser images of biological crystal networks with the following possibility of Mueller-matrix selection of polarization singularity. The technique of Mueller-matrix singular diagnostics of pathological changes in a woman's reproductive sphere tissue (myometrium) is proposed.



Laser polarimetry (polarization sensitive optical coherence tomography, Mueller matrix decomposition formalism, coherency matrix analysis, etc.) 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29 is able to obtain information about optical anisotropy (linear retardance, optical activity, dichroism, etc.) of biological tissues (BT). It has the potential to be an important technology for noninvasive diagnostics of organic phase-inhomogeneous layers. One model approach underlies these researches, which generalizes the optical properties of BT:5, 6, 14, 21, 30

  1. All varieties of human BT can be represented by four main types—connective, muscular, epithelial, and neural tissues.

  2. The morphological structure of any BT type is regarded as a two-component amorphous-crystalline structure (optical isotropic and optical anisotropic).

  3. The crystalline component is characterized by Mueller-matrix {M} operators of an optical anisotropy,

    Eq. 1

    [TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} \left\{ M \right\} = \left\| {\begin{array}{c@{\quad}c@{\quad}c@{\quad}c} 1 & 0 & 0 & 0 \\[3pt] 0 & {M_{22} } & {M_{23} } & {M_{24} } \\[3pt] 0 & {M_{32} } & {M_{33} } & {M_{34} } \\[3pt] 0 & {M_{42} } & {M_{43} } & {M_{44} } \\ \end{array}} \right\|. \end{equation}\end{document} M=10000M22M23M240M32M33M340M42M43M44.

Specifically, the above-mentioned model was used for finding and substantiating the interrelations between the ensemble of statistic moments of the first to fourth orders that characterize the orientation-phase structure (distribution of optical axes and phase shifts for directions of protein fibril networks) of birefringent BT architectonics and of 2-D distributions for azimuths and ellipticities in their laser images.5 It was determined 6, 7, 14, 15, 20, 23, 27, 28 that the third and fourth statistic moments for coordinate distributions of ellipticities are the most sensitive to the change (caused by dystrophic and oncological processes) of optical anisotropy inherent to protein crystals. On this basis, the criteria for early diagnostics of muscle dystrophy, precancer states of connective tissue, collagenoses, etc., were determined.

In parallel with traditional statistical investigations, formed in the most recent 10–15 years is the new optical approach to describe a structure of polarization inhomogeneous fields, in the case of scattered coherent radiation. The main feature of this approach is the analysis of definite (in contrast to continuous 2-D distributions) polarization states to determine the whole structure of coordinate distributions for azimuths and ellipticities of polarization. The so-called polarization singularities are commonly used as these states: 31, 32, 33, 34, 35, 36, 37 states with linear polarization, when the direction of rotation for the electric field vector is indefinite, the so-called L-points, and circularly polarized states, when the azimuth of polarization for the electric field vector is indefinite, the so-called C-points.

Investigations of polarization inhomogeneous object fields for BT with different morphology allowed ascertaining that they possess a developed network of L- and C-points.36 For example in Ref. 37, the authors found interrelations between conditions providing formation of polarization singular points and particularity of the orientation-phase structure of biological crystals present in the territorial matrix of the human tissue architectonic network. These interrelations served as a base to make statistical and fractal analyses of distribution densities for the number of singular points in BT images. As a result, the authors confirmed the efficiency of this method for investigation of object fields to differentiate optical properties of BT with a different morphological structure and physiological state.

It is worth noting that the singular approach is mainly used out of the analysis of the mechanisms providing formation of polarization inhomogeneous laser images of BT by an extracellular matrix. Thus, development of laser-polarimetric techniques based on determination of singular interrelations “object-field” in order to find new methods of diagnostics of transformation of the BT extracellular matrix orientation-phase structure related with precancer changes in their physiological state is very important.

To solve this problem, we should revert to the analysis of optical properties of biological crystal networks because these properties are comprehensively described by the Mueller matrix within the framework of the singular approach.


Brief Theory of the Singular Approach in the Analysis of Biological Tissue Birefringent Networks

To analytically describe the L- and C-states of polarization, the most suitable is to use the extreme values of the fourth Stokes vector parameter. It is a widely applied mean and appears to be reasonable to represent these singularities as follows:

Eq. 2

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} \left\{ \begin{array}{l} V_4 = \sin 2\beta, \\[3pt] V_4 = 0 \leftrightarrow L(\beta = 0), \\[3pt] V_4 = \pm 1 \leftrightarrow \pm C\left({\beta = \pm \frac{\pi }{4}} \right)\hspace*{-1.5pt}. \\ \end{array} \right. \end{equation}\end{document} V4=sin2β,V4=0L(β=0),V4=±1±Cβ=±π4.

Here, β is the value of polarization ellipticity; the +C-point is assumed to designate a right-circulated polarization state, which can be characterized by the phase shift δ = π/2 between the orthogonal components of laser beam amplitude; the –C-point is assumed to designate a left-circular polarization state [δ = −(π/2)]. Thus, ±C-points are the orthogonal states of a circularly polarized wave. For the L-points, the phase shift reaches δ = 0; π.

Using Eqs. 1, 2, it is possible to determine the interrelations between the characteristic values [TeX:] $M_{ik}^ *$ Mik* of Mueller-matrix elements, which correspond to the “extreme” values of BT optical anisotropy, and to the extreme values of V 4. Thus, in order to characterize the BT structure, we obtain the coordinate network generated by a finite number of characteristic values of Mueller-matrix elements [Eq. 1] (see Table 1).

Fig. 7

Distributions of the amount of singularities N(x), N +(x), N (x) of myometrium tissues of (a,d,g) A, (b,e,h) B, and (c,f,i) C types.


Table 1

Interrelation between the characteristic values of Mueller-matrix M ik elements of biological tissues and polarization singularities described by the fourth parameter of the Stokes vector V 4.

MikV4Polarization state

Thus, measuring the coordinate distributions of the characteristic values ( [TeX:] $M_{ik}^ * = 0, \pm 1$ Mik*=0,±1 ) of the BT Mueller-matrix elements enables not only to foresee the scenario ( [TeX:] $M_{ik}^ * \to V_4^ *$ Mik*V4* ) of forming the ensemble of polarization singularities (V 4 = 0, ±1) of its image, but also to additionally realize their differentiation, conditioned by the specificity of orientation structure of biological crystals.

Here, +L-point corresponds to a linear polarization state with δ = 0; –L-point corresponds to a linear polarization state with δ = π. In this case, as for ±C-points, we will identify ±L-points as orthogonal.


Scheme and Methods of Experimental Investigations

Figure 1 shows the traditional optical scheme of polarimeter for measuring the elements of Mueller matrix of the BT histological sections.5, 20

Fig. 1

Optical scheme of the polarimeter: 1 is the He-Ne laser; 2, the collimator; 3, the stationary quarter-wave plates; 5 and 8, mechanically movable quarter-wave plates; 4 and 9, polarizer and analyzer, correspondingly; 6, object of the investigation; 7, micro-objective; 10, the CCD camera; and 11, a personal computer.


The parallel (ø = 104 μm) beam of He-Ne laser (λ = 0.6328 μm, W = 5.0 μW) was used as an illuminator. The polarization illuminator consists of the quarter-wave plates 3, 5, and polarizer 4, and it sequentially forms a series of linearly polarized (I 0, I 45, I 90, I 135) with azimuths 0°, 90°, 45°, 135°, and right-hand (I ) and left-hand (I ) circularly polarized probing BT laser beams. The BT images made by the micro-objective 7 (4×) were projected into the plane of a light-sensitive area (800 × 600 pixels) of CCD camera 10.

Polarization images of BT were projected by means of the micro-objective 7 (focal distance – 1.5 cm, aperture – 0.2, magnification – 4×) into the plane of light-sensitive area of the CCD camera (overall amount of pixels, 800 × 600; light-sensitive area size, 4000 × 3000 μm; deviation of photosensitive characteristics from the linear one was no more than 15%), which provided the range of measuring the structural elements of BT with the resolution 2–2000 μm. Maximal resolution verification (2 μm) was performed using the stage micrometer (linear scale), where the image was projected into the light-sensitive area of the CCD camera using the micro-objective 7. Minimal resolution (2000 μm) corresponds to the situation when the light-sensitive area of the CCD camera is entirely filled by two equal-sized structural elements (light and dark) of a stage micrometer. The experimental conditions were chosen in such a way that it enabled one to reduce the space-angular aperture filtering when forming the BT images. It was ensured by conformance of angular characteristics of the indicatrix of light scattering by the BT samples (Ω ≈ 160) and the angular aperture of the micro-objective (Δω = 200). Here, Ω is the solid angle within which 98% of all the energy of light-scattered radiation is concentrated.

Polarization analysis of the BT images was performed using the polarizer 9 and quarter-wave plate 8 according to the following technique:

  1. Within the section of the illuminating laser beam, the array (m × n = 800 × 600) of values for the Stokes vector V j = 1,2,3,4 parameters and elements of Mueller matrix M ik(m × n) were determined in accord with the following algorithms:

    Eq. 3

    [TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} \left\{ {\begin{array}{*{20}c} {V_1 = I_0 + I_{90},} & {V_2 = I_0 - I_{90},} & {V_3 = I_{45} - I_{135},} & {V_4 = I_ \otimes - I_ \oplus,} \\[5pt] {M_{i1} = 0.5\big[ {V_i^{(1)} + V_i^{(2)} } \big],} & {M_{i2} = 0.5\big[ {V_i^{(1)} - V_i^{(2)} } \big],} & {M_{i3} = V_i^{(3)} - M_{i1},} & {M_{i3} = V_i^{(4)} - M_{i1},\,\,\,\,i = 1,2,3,4.} \\ \end{array}} \right. \end{equation}\end{document} V1=I0+I90,V2=I0I90,V3=I45I135,V4=II,Mi1=0.5Vi(1)+Vi(2),Mi2=0.5Vi(1)Vi(2),Mi3=Vi(3)Mi1,Mi3=Vi(4)Mi1,i=1,2,3,4.

  2. In each array M ik(m × n) and V 4(m × n), coordinate distributions of characteristic (singular) values 0, ±1 were determined.

At the first stage, the interrelations ( [TeX:] $M_{ik}^ * \to V_4^ *$ Mik*V4* ) of matrix and polarization singularities were investigated on the sample of a histological section prepared from healthy skin derma. Figure 2 represents coordinate distributions of matrix elements M 44,24,34(m × n) corresponding to the skin derma histological section and the fourth Stokes vector parameter V 4(m × n) of its image with the characteristic values (0, ±1) plotted on them [within the marked (100 × 100)-pixel sampling plot].

Fig. 2

Networks of (a–c) characteristic values [TeX:] $M_{44,24,34}^* (m \times n)$ M44,24,34*(m×n) of matrix elements M 44,24,34 and (d) singularities of polarization image inherent to the skin derma layer histological section V 4: ± C-points (❑) (M 44 = 0) + L-points (Δ) (M 44 = + 1); + C-points (⧅) [TeX:] $ (M_{24,34} = + 1,\;V_4 = + 1); { \; - {\rm L} {\text -}{\rm points }\; (\nabla)\;(M_{44} = - 1);} $ (M24,34=+1,V4=+1);L-points()(M44=1); −C-points (⧄) (M 24,34 = −1, [TeX:] $V_4 = - 1); \; \pm {\rm L} {\text -} {\rm points }\; (\diamondsuit)\; (M_{24,34} = 0,$ V4=1);±L-points()(M24,34=0, V 4 = 0).


It can be seen from the data obtained that there is a direct correlation between the coordinate (k,g 1 ≤ km, 1 ≤ gn) positions of characteristic values of the matrix element [TeX:] $M_{44}^ *$ M44* for skin derma and the network of L- and C-points in its laser image [TeX:] $ {M_{44}^ * \left({k,g} \right) = \left\{ \begin{array}{l} 0 \\ \pm 1 \\ \end{array} \right\} \Leftrightarrow V_4^ * \left({k,g} \right) = \left\{ {\begin{array}{*{20}c} { \pm 1} & - & { \pm C} \\ 0 & - & L \\\end{array}} \right\}} $ M44*k,g=0±1V4*k,g=±1±C0L [Figs. 2a, 2d].

Coordinate distributions of characteristic values of matrix elements [TeX:] $M_{24,42}^ * \left({m,n} \right)$ M24,42*m,n , [TeX:] $M_{34,43}^ * \left({m,n} \right)$ M34,43*m,n and corresponding networks of orthogonal ±L- and ±C-points (Table 1) possess the individual structure [Figs. 2b, 2c].

Analytically substantiated and experimentally proven interrelations between the matrix and polarization singularities were used as the basis for Mueller-matrix singular diagnostics of oncological changes in the tissues of a woman's reproductive sphere.


Mueller-Matrix Singular Diagnostics and Differentiation of Pathological Changes in the Tissues of Woman Reproductive Sphere

Three groups of histological sections of the main tissue of a woman's reproductive sphere—myometrium—were used as the objects of the investigation: biopsy of the healthy tissue of a woman's reproductive sphere [type A, Fig. 3a], biopsy of the inflamed tissue (ectonia) [type B, Fig. 3b], and biopsy of the tissue in the state of dysplasia (precancer state) [type C, Fig. 3c].

Fig. 3

Polarization images of woman reproductive sphere tissue – myometrium (a) A, (b) B, and (c) C types in the scheme of coaxial polarizer and analyzer.


To determine the criteria of Mueller-matrix singular diagnostics of the myometrium oncological state and differentiate its severity degree, the following technique was used:

  1. Coordinate networks of characteristic values for matrix elements [TeX:] $M_{44,24,34}^ * (m \times n) = 0{\rm,}\; \pm 1$ M44,24,34*(m×n)=0,±1 were scanned in the direction x ≡ 1,. . .,m with the step Δx = 1 pixel.

  2. Within the obtained sampling (1pixel × n pixel)(k = 1,2,. . .,m) for coordinate distribution of the element M 44(m × n), the total amount [N (k)] of characteristic points (0, ±1) that set the complete ensemble of singular points was calculated, and the dependences N(x) ≡ [N (1), N (2),…,N (m)] were determined.

  3. Distributions of the number of singular ±L- and ±C-points were determined using the following expressions:

    Eq. 4

    [TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} \left\{ \begin{array}{l} N^{\left(+ \right)} (x) = N_C ({M_{34,43} = + 1}) + N_L ({M_{24,42} = 0})\hspace*{-1pt}, \\[5pt] N^{\left(- \right)} (x) = N_L({M_{34,43} = 0}) + N_C ({M_{24,42} = - 1})\hspace*{-1pt}. \\ \end{array} \right. \end{equation}\end{document} N+(x)=NC(M34,43=+1)+NL(M24,42=0),N(x)=NL(M34,43=0)+NC(M24,42=1).

  4. Statistical moments of the first to fourth orders for the obtained distributions of N(x) amount of singularities were calculated according to the following algorithms:

    Eq. 5

    [TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} \begin{array}{*{20}c} {Z_1 = \displaystyle\frac{1}{{m \times n}}\sum\limits_{i = 1}^{m \times n} {\left| {N\left(x \right)} \right|,} } \\[12pt] {Z_2 = \sqrt {\displaystyle\frac{1}{{m \times n}}\sum\limits_{i = 1}^{m \times n} {\left[ {N\left(x \right)} \right]^2 } },} \\[12pt] {Z_3 = \displaystyle\frac{1}{{Z_2^3 }}\displaystyle\frac{1}{{m \times n}}\sum\limits_{i = 1}^{m \times n} {\left[ {N\left(x \right)} \right]^3,} } \\[12pt] {Z_4 = \displaystyle\frac{1}{{Z_2^4 }}\displaystyle\frac{1}{{m \times n}}\sum\limits_{i = 1}^{m \times n} {\left[ {N\left(x \right)} \right]^4.} } \\ \end{array} \end{equation}\end{document} Z1=1m×ni=1m×nNx,Z2=1m×ni=1m×nNx2,Z3=1Z231m×ni=1m×nNx3,Z4=1Z241m×ni=1m×nNx4.

Figures 4, 5, 6 show the networks of characteristic values [TeX:] $M_{44,24,34}^ * (m \times n)$ M44,24,34*(m×n) inherent to coordinate distributions for matrix elements M 44, 24, 34(m × n) of histological sections of myometrium tissues of A, B, and C types.

Fig. 4

Networks of (a–c) characteristic values [TeX:] $M_{44,24,34}^ * (m \times n)$ M44,24,34*(m×n) of matrix elements M 44,24,34 and (d) singularities of polarization image of myometrium histological section of A type V 4: ±C-points (❑) (M 44 = 0); + L-points (Δ) (M 44 = + 1); + C-points (⧅) (M 24,34 = +1, V 4 = + 1); −L-points (∇) (M 44 = −1); −C-points (⧄) (M 24,34 = −1, [TeX:] $ V_4 = - 1); \; \pm {\rm L} {\text -} {\rm points }\;(\diamondsuit)\;(M_{24,34} = 0,\;V_4 = 0).$ V4=1);±L-points()(M24,34=0,V4=0).


Fig. 5

Networks of (a–c) characteristic values [TeX:] $M_{44,24,34}^ * (m \times n)$ M44,24,34*(m×n) of matrix elements M 44,24,34 and (d) singularities of polarization image of myometrium histological section of B type V 4: ±C-points (❑) (M 44 = 0); [TeX:] $+ {\rm L} {\text -} {\rm points }\;(\Delta)\;(M_{44}\, =\, + 1)$ +L-points(Δ)(M44=+1) ; +C-points (⧅) (M 24,34 = + 1, V 4 = + 1); [TeX:] $- {\rm L} {\text -} {\rm points }\;(\nabla)\;(M_{44} = - 1)$ L-points()(M44=1) ; −C-points (⧄) (M 24,34 = −1, [TeX:] $V_4\break = - 1); \pm {\rm L} {\text -} $ V4=1);±L- points (◇) (M 24, 34 = 0, V 4 = 0).


Fig. 6

Networks of (a–c) characteristic values [TeX:] $M_{44,24,34}^ * (m \times n)$ M44,24,34*(m×n) of matrix elements M 44,24,34 and (d) singularities of polarization image of myometrium histological section of C type V 4: ±C-points (❑) (M 44 = 0); [TeX:] $+ {\rm L} {\text -} {\rm points }\;(\Delta)\;(M_{44} = + 1)$ +L-points(Δ)(M44=+1) ; +C-points (⧅) (M 24, 34 = +1, V 4 = +1); [TeX:] $- {\rm L} {\text -} {\rm points }\;(\nabla)\;(M_{44} = - 1);$ L-points()(M44=1); −C-points (⧄) (M 24,34 = −1, V 4 = −1); ±L-points (◇) (M 24,34 = 0, V 4 = 0).


Figure 7 illustrates the distributions of the number of singularities N(x), N +(x), N (x) of myometrium tissues of (a,d,g) A, (b,e,h) B, and (c,f,i) C types. The comparative analysis of the data obtained shows that the following:

  1. The coordinate distributions of the elements M 44, 24, 34(m × n) of the Mueller matrix corresponding to myometrium tissue of all types is characterized by individual (according to quantitative and topological structure) networks of singular points (Figs. 4, 5, 6).

  2. The total amount of ±C-points [ [TeX:] $M_{44}^ * (m \times n) = 0$ M44*(m×n)=0 ] sequentially increases for the samples of myometrium of A, B, and C types [Figs. 4a, 4b, 4c, 6a].

  3. The dependences N +(x) of the number of singular values of matrix elements for the samples of myometrium tissue of all the types are similar in their structure [Figs. 7d, 7e, 7f].

  4. The distributions N (x) for the samples of myometrium tissue of B type are characterized by sufficient increase (two- to threefold) of the number of singular values in comparison to similar dependencies N +(x) [Figs. 7e, 7f, 7g, 7h].

The obtained results can be related with the increase of birefringence (Δn ≈ 1.5 × 10−2) of collagen and myosin fibrils in pathologically changed myometrium of B and C types. Besides, at early stages (ectonia), the directions of the growth of newly formed fibrils are being formed. In dysplasia states, these pathologically changed fibrils form a specifically oriented network of biological crystals.

In terms of physics, these morphological processes are manifested in the increase of probability of forming the ±C-points (myometrium samples of B and C types), as well as in appearance of asymmetry between ranges of dependences values N +(x) and N (x), which characterize the number of orthogonal L- and C-points (Table 1).

Finally, comparative investigations of the diagnostic efficiency typical to famous techniques of laser polarimetry {Z 1, 2, 3, 4[M 44, 34, 24(m × n)]};14 polarization-correlation mapping26, 27

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation*} {Z_{1,2,3,4} \left[ {V_4 \left({m \times n} \right) = \left\{ \begin{array}{l} 0, \\[3pt] \pm 1 \\ \end{array} \right.} \right]} \end{equation*}\end{document} Z1,2,3,4V4m×n=0,±1
and the technique of Mueller-matrix singular diagnostics Z 1,2,3,4[N(x), N +(x), N (x)] were suggested.

Summarized in Table 2 are the statistical averaged values within these three groups of samples of myometrium tissue {Z 1, 2, 3, 4[M 44, 34, 24(m × n)]},

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation*} {Z_{1,2,3,4} \left[ {V_4 \left({m \times n} \right) = \left\{ \begin{array}{l} 0, \\[3pt] \pm 1 \\ \end{array} \right.} \right]} \hspace*{-1.5pt}, \end{equation*}\end{document} Z1,2,3,4V4m×n=0,±1,
and Z 1,2,3,4[N(x), N +(x), N (x)].

Table 2

Values {Z 1, 2, 3, 4[M 44, 34, 24(m × n)]}, $ {Z_{1,2,3,4} \left[ {V_4 \left({m \times n} \right) = \left\{ \begin{array}{l} 0, \\ \pm 1 \\ \end{array} \right.} \right]} $ Z1,2,3,4V4m×n=0,±1 , and Z 1,2,3,4[N(x), N +(x), N −(x)] statistically averaged within the three groups of myometrium tissue samples.

A typeB typeC type
(21 samples)(21 samples)(21 samples)
Zj = 1,2,3,4M44M34M24M44M34M24M44M34M24
Z10.67 ± 0.0590.32 ± 0.0310.27 ± 0.020.59 ± 0.0480.27 ± 0.020.24 ± 0.050.37 ± 0.030.19 ± 0.0150.18 ± 0.019
Z20.51 ± 0.0460.29 ± 0.0190.26 ± 0.0210.57 ± 0.050,23 ± 0.010.21 ± 0.010.28 ± 0.020.21 ± 0.0180.17 ± 0.015
Z31.13 ± 0.120.88 ± 0.090.74 ± 0.080.98 ± 0.10.79 ± 0.080.66 ± 0.0690.66 ± 0.070.49 ± 0.0420.41 ± 0.04
Z43.15 ± 0.322.11 ± 0.292.27 ± 0.312.84 ± 0.241.79 ± 0.181.87 ± 0.171.57 ± 0.141.07 ± 0.11.12 ± 0.11
Zj = 1,2,3,4V4 = 0V4 = ±1V4 = 0V4 = ±1V4 = 0V4 = ±1
Z10.12 ± 0.0790.24 ± 0.0380.15 ± 0.0670.28 ± 0.0680.18 ± 0.0710.29 ± 0.076
Z20.16 ± 0.0740.31 ± 0.0420.19 ± 0.0310.38 ± 0.0520.23 ± 0.0190.4 ± 0.048
Z30.70 ± 0.0520.92 ± 0.0860.93 ± 0.0941.12 ± 0.1011.27 ± 0,1121.72 ± 0.123
Z41.71 ± 0.132.19 ± 0.182.01 ± 0.193.13 ± 0.273.41 ± 0.314.01 ± 0.31
Zj = 1,2,3,4N(x)N+(x)N(x)N(x)N+(x)N(x)N(x)N+(x)N(x)
Z10.61 ± 0.0520.43 ± 0.0380.12 ± 0.0340.39 ± 0.050.37 ± 0.0420.23 ± 0.040.29 ± 0.0390.28 ± 0.0360.19 ± 0.038
Z20.75 ± 0.0680.82 ± 0.0760.15 ± 0.020.36 ± 0.0420.69 ± 0.0560.46 ± 0.020.1 ± 0.020.55 ± 0.0490.23 ± 0.026
Z31.19 ± 0.150.92 ± 0.011.86 ± 0.190.63 ± 0.0510.87 ± 0.071.86 ± 0.190.53 ± 0,040.76 ± 0.0622,06 ± 0.21
Z41.99 ± 0.172.31 ± 0.192.32 ± 0.211.01 ± 0.12.07 ± 0.178.45 ± 0.730.8 ± 0.071.87 ± 0.192.91 ± 0.32

It follows from the data presented that:

  1. Efficiency of laser polarimetry for diagnostics and differentiation of early oncological changes of myometrium tissue is insufficient—the difference between the values of statistical moments {Z 1, 2, 3, 4[M 44, 34, 24(m × n)]} of samples A, B, and C types is insignificant and does not exceed 20–45%;

  2. The technique of polarization-correlation mapping is efficient for differentiation of optical properties of healthy and oncological changed myometrium tissue—the skewness (Z 3) and the kurtosis (Z 4) of distribution of the number of singular points of A and B types of laser images differ by 1.53 and 2.15 times.

  3. The technique of Mueller-matrix singular diagnostics is efficient for differentiation of optical properties of all types of samples—statistical moments of the third and fourth orders of distributions N(x) for samples A, B, and C types differ by 1.7 and 2.5 times, respectively.

  4. For distributions N (x) of myometrium tissue of A and B types, the maximal difference (from 2.2 to 4.1 times) is observed between all the statistical Z j = 1, 2, 3, 4 moments.



It has been ascertained that correlation between coordinate locations of characteristic points for 2-D elements of the Mueller matrix corresponding to an optically thin layer of biological tissue and the network of L- and C-points in its laser image. Shown is the potentiality of Mueller-matrix sampling for polarization singularities formed by biological crystals. The efficiency of Mueller-matrix singular diagnostics has been demonstrated, not only for oncological changes of myometrium tissue, but also for differentiating their degree of severity.



J. F. de Boer and T. E. Milner, “Review of polarization sensitive optical coherence tomography and Stokes vector determination,” J. Biomed. Opt., 7 359 –371 (2002). Google Scholar


J. F. de Boer, T. E. Milner, and J. S. Nelson, Trends in Optics and Photonics (TOPS): Advances in Optical Imaging and Photon Migration, OSA, Washington (1998). Google Scholar


M. J. Everett, K. Shoenenberger, B. W. Colston, and L. B. Da Silva, “Birefringence characterization of biological tissue by use of optical coherence tomography,” Opt. Lett., 23 228 –230 (1998). Google Scholar


Sh. Jiao, W. Yu, G. Stoica, and L. V. Wang, “Optical-fiber based Mueller optical coherence tomography,” Opt. Lett., 28 1206 –1208 (2003). Google Scholar


A. G. Ushenko and V. P. Pishak, “Laser Polarimetry of Biological Tissue: Principles and Applications,” Handbook of Coherent-Domain Optical Methods: Biomedical Diagnostics, Environmental and Material Science, 93 –138 Kluwer Academic Publishers, Boston (2004). Google Scholar


O. V. Angelsky, A. G. Ushenko, Yu. A. Ushenko, and V. P. Pishak, “Statistical and Fractal Structure of Biological Tissue Mueller Matrix Images,” Optical Correlation Techniques and Applications, 213 –266 SPIE Press, Bellingham, WA (2007). Google Scholar


O. V. Angelsky, A. G. Ushenko, Yu. A. Ushenko, V. P. Pishak, and A. P. Peresunko, “Statistical, correlation, and topological approaches in diagnostics of the structure and physiological state of birefringent biological tissues,” Handbook of Photonics for Biomedical Science, 21 –67 CRC Press, Boca Raton (2010). Google Scholar


A. G. Ushenko, “Polarization structure of laser scattering fields,” Opt. Eng., 34 1088 –1093 (1995). Google Scholar


A. G. Ushenko, “Laser diagnostics of biofractals,” Quantum Electron., 29 (12), 1078 –1084 (1999). Google Scholar


O. V. Angel'skii, A. G. Ushenko, A. D. Arkhelyuk, S. B. Ermolenko, and D. N. Burkovets, “Structure of matrices for the transformation of laser radiation by biofractals,” Quantum Electron., 29 (12), 1074 –1077 (1999). Google Scholar


O. V. Angel'skii, A. G. Ushenko, A. D. Arheluk, S. B. Ermolenko, and D. N. Burkovets, “Scattering of laser radiation by multifractal biological structures,” Opt. Spectrosc., 88 (3), 444 –448 (2000). Google Scholar


A. G. Ushenko, “Polarization structure of biospeckles and the depolarization of laser radiation,” Opt. Spectrosc., 89 (4), 597 –601 (2000). Google Scholar


A. G. Ushenko, “Stokes correlometry of biotissues,” Laser Phys., 10 (5), 1286 –1292 (2000). Google Scholar


A. G. Ushenko, “The vector structure of laser biospeckle fields and polarization diagnostics of collagen skin structures,” Laser Phys., 10 (5), 1143 –1149 (2000). Google Scholar


A. G. Ushenko, “Laser polarimetry of polarization-phase statistical moments of the object field of optically anisotropic scattering layers,” Opt. Spectrosc., 91 (2), 313 –316 (2001). Google Scholar


A. G. Ushenko, “Polarization contrast enhancement of images of biological tissues under the conditions of multiple scattering,” Opt. Spectrosc., 91 (6), 937 –940 (2001). Google Scholar


A. G. Ushenko, “Laser probing of biological tissues and the polarization selection of their images,” Opt. Spectrosc., 91 (6), 932 –936 (2001). Google Scholar


A. G. Ushenko, “Correlation processing and wavelet analysis of polarization images of biological tissues,” Opt. Spectrosc., 91 (5), 773 –778 (2002). Google Scholar


A. G. Ushenko, “Polarization correlometry of angular structure in the microrelief pattern or rough surfaces,” Opt. Spectrosc., 92 (2), 227 –229 (2002). Google Scholar


O. V. Angelsky, A. G. Ushenko, Ye. G. Ushenko, “2-D Stokes polarimetry of biospeckle tissues images in pre-clinic diagnostics of their pre-cancer states,” J. Holography Speckle, 2 (1), 26 –33 (2005). Google Scholar


O. V. Angelsky, A. G. Ushenko, and Y. G. Ushenko, “Complex degree of mutual polarization of biological tissue coherent images for the diagnostics of their physiological state,” J. Biomed. Opt., 10 (6), 060502 (2005). Google Scholar


O. V. Angelsky, A. G. Ushenko, and Ye. G. Ushenko, “Investigation of the correlation structure of biological tissue polarization images during the diagnostics of their oncological changes,” Phys. Med. Biol., 50 (20), 4811 –4822 (2005). Google Scholar


A. Ushenko, S. Yermolenko, A. Prydij, S. Guminetsky, I. Gruia, O. Toma, and K. Vladychenko, “Statistical and fractal approaches in laser polarimetry diagnostics of the cancer prostate tissues,” Proc. SPIE, 7008 70082C (2008). Google Scholar


A. G. Ushenko, A. I. Fediv, and Yu. F. Marchuk, “Correlation and fractal structure of Jones matrices of human bile secret,” Proc. SPIE, 7368 73681Q (2009). Google Scholar


O. V. Angelsky, A. G. Ushenko, Y. G. Ushenko, and Y. Y. Tomka, “Polarization singularities of biological tissues images,” J. Biomed. Opt., 11 (5), 054030 (2006). Google Scholar


S. B. Yermolenko, A. G. Ushenko, and P. Ivashko, “Spectropolarimetry of cancer change of biotissues,” Proc. SPIE, 7388 73881D (2009). Google Scholar


A. G. Ushenko, I. Z. Misevich, V. Istratiy, I. Bachyns'ka, A. P. Peresunko, O. K. Numan, and T. G. Moiysuk, “Evolution of statistic moments of 2D-distributions of biological liquid crystal net Mueller matrix elements in the process of their birefringent structure changes,” Adv. Opt. Technol., 2010 423145 (2010). Google Scholar


Yu. O. Ushenko, I. Z. Misevich, A. P. Angelsky, V. T. Bachinsky, O. Yu. Telen'ga, and O. I. Olar, “Polarization-singular structure in laser images of phase-inhomogeneous layers to diagnose and classify their optical properties,” Semicond. Phys., Quantum Electron. Optoelectron., 13 (3), 248 –258 (2010). Google Scholar


V. T. Bachinsky, Yu. O. Ushenko, Yu. Ya. Tomka, O. V. Dubolazov, V. O. Balanets'ka, and A. V. Karachevtsev, “Wavelet analysis for polarization maps of networks formed by liquid biological crystals in blood plasma: statistical and fractal approaches,” Semicond. Phys., Quantum Electron. Optoelectron., 13 (2), 189 –201 (2010). Google Scholar


S. C. Cowin, “How is a tissue built,” J. Biomed. Eng., 122 553 –568 (2000). Google Scholar


J. F. Nye, Natural Focusing and the Fine Structure of Light, Institute of Physics, Bristol (1999). Google Scholar


J. F. Nye, “Polarization effects in the diffraction of electromagnetic waves: the role of disclinations,” Proc. R. Soc. A, 387 105 –132 (1983). Google Scholar


M. V. Berry and M. R. Dennis, “Polarization singularities in isotropic random vector waves,” Proc. R. Soc. A, 457 141 –155 (2001). Google Scholar


J. F. Nye, “Lines of circular polarization in electromagnetic wave fields,” Proc. R. Soc. A, 389 279 –290 (1983). Google Scholar


Coherency matrix description of optical polarization singularities,” J. Opt. A., 6 S229 –S234 (2004). Google Scholar


O. V. Dubolazov, A. G. Ushenko, V. T. Bachynsky, A. P. Peresunko, and O. Ya. Vanchulyak, “On the feasibilities of using the wavelet analysis of Mueller matrix images of biological crystals,” Adv. Opt. Technol., 2010 162832 (2010). Google Scholar


A. G. Ushenko, A. I. Fediv, and Yu. F. Marchuk, “Singular structure of polarization images of bile secret in diagnostics of human physiological state,” Proc. SPIE, 7368 73681S (2009). Google Scholar
©(2011) Society of Photo-Optical Instrumentation Engineers (SPIE)
Yuriy A. Ushenko "Investigation of formation and interrelations of polarization singular structure and Mueller-matrix images of biological tissues and diagnostics of their cancer changes," Journal of Biomedical Optics 16(6), 066006 (1 June 2011).
Published: 1 June 2011

Back to Top