The paper deals with investigation of the processes of laser radiation transformation by biological crystals networks using the singular optics techniques. The results obtained showed 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) of laser transformation of biological crystals networks with the following possibility of Mueller-matrix selection of polarization singularity. The technique of Mueller-matrix diagnostics of pathological changes of skin derma is proposed.
1. Introduction
Laser polarimetry [1] enabling to obtain information about optical anisotropy [2–5] of biological tissues (BT) is an important direction of noninvasive diagnostics of organic phase-heterogeneous layers. For statistic analysis of such polarimetric information a model approach has been worked out based on the following conditions [1, 2, 6–12]:
all the variety of human BT can be represented by four main types—connective, muscular epithelial, and neural tissues;
structure of any BT type is regarded as a two-component amorphous-crystalline one;
the crystalline component or extracellular matrix is formed by the network of optically uniaxial birefringent protein (collagen, myosin, elastine, etc.) fibrils or biological crystals;
the process of transformation of laser radiation polarization state by biological crystal is characterized by Mueller {M} matrix operators of an optically uniaxial crystal {M}=∥10000M22M23M240M32M33M340M42M43M44∥=∥10000(cos22ρ+sin22ρcosδ)cos2ρsin2ρ(1-cosδ)(sin2ρsinδ)0(cos2ρsin2ρ(1-cosδ))(sin22ρ+cos22ρcosδ)(cos2ρsinδ)0(-sin2ρsinδ)(-cos2ρsinδ)(cosδ)∥.
Here ρ-direction of optical axis of biological crystal with birefringence index Δn, δ=(2π/λ)Δnl-phase shift between orthogonal components of the amplitude of a probing laser beam with wave length λ.
A new approach to description of the BT laser images based on the analysis of coordinate distributions of polarization singularities became developed the above-mentioned statistical [13–22]. Linearly (L-points) and circularly (C-points) polarized states of light oscillations belong to them. For L-points the direction of the electric-intensity vector’s rotation is indefinite (singular). For a C-point, the polarization azimuth of the electric intensity vector is indefinite.
Investigation of laser images of the connective tissue layers revealed a developed network of polarization singularities [23–26], which was quantitatively estimated in the form of distribution of the amount of L- and C-points. By means of the analysis of the given distribution’s statistical moments of the 1st–4th orders (the technique of polarization mapping) the criteria of diagnostics of oncological changes of uterus neck tissue were found.
It should be pointed out that singular approach is predominantly realized out of the analysis of the mechanisms of forming polarizationally heterogeneous laser images of BT by an extracellular matrix. Thus, development of laser polarimetry techniques based on determination of singular interconnections “object-field” in order to find new methods of diagnostics of transformation of the BT extracellular matrix orientation-phase structure connected with precancer changes of their physiological state is very important.
To solve such a problem, we should revert to the analysis of optical properties of biological crystals’ nets, comprehensively described by the Mueller matrix though within a singular approach.
2. A Brief Theory of the Mueller Matrixes Approach in the Analysis of the Biological Tissue Birefringent Nets Polarization Properties
The use of the fourth parameter of the Stokes vector appears to be a suitable and widely applied means of such singularities representation
V4={0⟷L±1⟷C.
According to analysis of (1) and (2), one can see the interconnection between the polarization singular states and certain (characteristic) values of orientation ρ* and phase δ* parameters of the BT crystals’ nets of the extracellular matrix
ρ*=00;±450900;δ*=0∘,90∘,180∘.
As it can be seen, relations (3) are the necessary terms for forming polarization singular states of the laser beam (L- (δ=0∘,180∘) and C- (δ=±90∘) points) by optically coaxial birefringent crystal.
Considering expressions (1)–(3) the characteristic values Mik* were defined that determine the L- and C-points in laser image of the extracellular matrix of the BT layer:
the values M44=0 and V4=±1 determine the complete set of ± C-points (δ=±90∘);
the complete set of L-points (δ=0∘) of the laser image is caused by the terms M22=M33=M44=1 and V4=0.
Mueller-matrix analysis enables to perform the sampling of polarization singularities of the laser image, formed by biological crystals with orthogonally oriented (δ=0∘, 90∘ and δ=45∘, 135∘) optical axes to
“orthogonal” L0;90- and L45;135- points
M24,42=0-L0∘,90∘-(ρ=0∘,90∘),M34,43=0-L45∘,135∘-(ρ=45∘,135∘).
Thus, measuring the coordinate distributions of the characteristic values (Mik*=0,±1) of the BT Mueller matrix elements enables not only to foresee the scenario (Mik*→V4*) of forming the ensemble of polarization singularities (V4=0,±1) of its image, but also to additionally realize their differentiation, conditioned by the specificity of orientation structure of biological crystals.
3. The Scheme and Methods of Experimental Investigations
Figure 1 shows traditional optical scheme of polarimeter for measuring due to Gerrard technique [27] of Stokes parameters and elements of Mueller matrix of the BT histological sections.
Optical scheme of polarimeter, where 1: He-Ne laser; 2: collimator; 3: stationary quarter-wave plates; 5, 8: mechanically movable quarter-wave plates; 4, 9: polarizer and analyzer correspondingly; 6: object of investigation; 7: microobjective; 10: CCD camera; 11: personal computer.
The parallel (⌀=104μm) beam of He-Ne laser (λ=0.6328μm, W=5.0μW) was used as an illuminator. Polarization illuminator consists of quarter-wave plates 3, 5 and polarizer 4, and it sequentially forms a series of linearly polarized (I0, I45, I90, I135) 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 microobjective (4×) 7 were projected into the plane of a light-sensitive plate (m×n=800×600 pixels) of CCD-camera 10. Polarization analysis of the BT images was performed by means of polarizer 9 and quarter-wave plate 8.
The optical thin (the absorption coefficient τ<0,1) BT histological sections were used as the objects of investigation. In this situation, one has a single scattering regime of laser radiation scattered by BT network and the narrow-band scattering indicatrix is formed (95% of energy is concentrated within the angle cone ΔΩ≤150). Therefore, the speckle background formation in the BT histological section image due to scattering on optical elements is insignificant.
At the first stage the interconnections (Mik*→V4*) of matrix and polarization singularities were investigated on the example of histological section of healthy skin derma layer.
Figure 2 represents coordinate distributions of matrix elements M44,24,34(m×n) of histological section of skin derma and the fourth parameter V4(m×n) of its image’s Stokes vector with the characteristic values 0, ±1 plotted on them (within the marked 100pix×100pix sampling plot).
Networks of characteristic values M44,24,34*(m×n) of matrix elements M44,24,34 (a, b, c) and singularities of polarization image of the skin derma layer histological section V4 (“d”): “±C”-points (□) (M44 = 0); “+C”-points () (M24,34 = +1, V4 = +1); “-C”-points () (M24,34 = -1, V4 = -1); “+L”-points (Δ) (M44 = +1); “-L”-points (∇) (M44 = -1); “±L”-points (♢) (M24,34= 0, V4 = 0).
It can be seen from the data obtained that there is direct correlation between the coordinate (k,g1≤k≤m,1≤g≤n) positions of characteristic values of the matrix element M44* of skin derma and the network of L- and C-points of its laser image {M44*(k,g)={0±1}⇔V4*(k,g)={±1-±C0-L}} (Figures 2(a) and 2(d)).
Coordinate distributions of characteristic values of matrix elements M24,42*(m,n), M34,43*(m,n) and corresponding networks of “orthogonal” L0,90-, L45,135- andC0,90-, C45,135-points (relations (4) and (5)) possess individual structure. Such peculiarities of singular networks of laser image of skin derma are obviously conditioned by the asymmetry of various directions (ρ=0∘, 90∘ and ρ=45∘, 135∘) of orientation of optical axes of biological crystals in the plane of the investigated sample (Figures 2(b) and 2(c)).
Analytically substantiated and experimentally proven interconnections between the matrix and polarization singularities were used as the basis for Mueller-matrix singular diagnostics of oncological changes of the tissues of women’s reproductive sphere.
4. Mueller-Matrix Diagnostics and Differentiation of Pathological Changes of the Skin Derma
Three groups of histological sections of the main tissue of skin derma—were used as the objects of investigation:
biopsy of the sound tissue of skin derma (type “A”—Figure 3(a));
biopsy of the skin derma in precancer state (type “B”—Figure 3(b));
biopsy of the skin derma in cancer state (type “C”—Figure 3(c)).
Polarization images of the skin derma “A” (a), “B” (b), and “C” (c)-types in coaxial polarizer and analyzer.
To determine the criteria of Mueller-matrix diagnostics of skin derma oncological state and differentiation of its severity degree the following technique was used:
coordinate networks of characteristic values of matrix elements M44,24,34*(m×n)=0,±1 were scanned in the direction x≡1,…,m with the step Δx=1 pixel;
within the obtained sampling (1pix×npix)(k=1,2,…,m) for coordinate distribution of the element M44(m×n) the total amount (N(k)) of characteristic points (0, ±1), which set the complete ensemble of singular points was calculated and the dependencies N(x)≡(N(1),N(2),…,N(m)) were determined;
distributions of the number of “orthogonal” singular L- and ±C-points were determined according to the terms (4) and (5);
ρ=0∘,90∘⇔N0,90(x)=NC(M34,43=±1)+NL(M24,42=0),ρ=45∘,135∘⇔N45,135(x)=NL(M34,43=0)+NC(M24,42=±1);
statistical moments of the 1st–4th orders of the obtained distributions of N(x) amount of singularities were calculated according to the algorithms
Z1=1m×n∑i=1m×n|N(x)|,Z2=1m×n∑i=1m×n[N(x)]2,Z3=1Z231m×n∑i=1m×n[N(x)]3,Z4=1Z241m×n∑i=1m×n[N(x)]4.
Figures 4, 5, and 6 show the networks of characteristic values M44,24,34*(m×n) of coordinate distributions of matrix elements M44,24,34(m×n) of histological sections of skin derma of “A”, “B”, “C”-types.
Networks of characteristic values M44,24,34*(m×n) of matrix elements M44,24,34 (a, b, c) and singularities of polarization image of skin derma histological section of “A”-type V4 (d): “±C”-points (□) (M44 = 0); “+C”-points () (M24,34 = +1, V4 = +1); “-C”-points () (M24,34 = -1, V4 = -1); “+L”-points (Δ) (M44 = +1); “-L”-points (∇) (M44 = -1); “±L”-points (♢) (M24,34 = 0, V4 = 0).
Networks of characteristic values M44,24,34*(m×n) of matrix elements M44,24,34 (a, b, c) and singularities of polarization image of skin derma histological section of “B”-type V4 (d): “±C”-points (□) (M44 = 0); “+C”-points () (M24,34 = +1, V4 = +1); “-C”-points () (M24,34 = -1, V4 = -1); “+L”-points (Δ) (M44 = +1); “-L”-points (∇) (M44 = -1); “±L”-points (♢) (M24,34 = 0, V4 = 0).
Networks of characteristic values M44,24,34*(m×n) of matrix elements M44,24,34 (a, b, c) and singularities of polarization image of skin derma histological section of “C”-type V4 (d): “±C”-points (□) (M44 = 0); “+C”-points () (M24,34 = +1, V4 = +1); “-C”-points () (M24,34 = -1, V4 = -1); “+L”- points (Δ) (M44 = +1); “-L”- points (∇) (M44 = -1); “±L”-points (♢) (M24,34 = 0, V4 = 0).
Figure 7 illustrates the distributions of the number of characteristic values N(x), N0,90(x), N45,135(x) of skin derma tissues of “A” (left column), “B” (central column), “C” (right column) types.
Distributions of the amount of characteristic values N(x), N0,90(x), N45,135(x) of skin derma tissues of “A”-(left column), “B”-(central column), “C” (right column)-types.
The comparative analysis of the data obtained shows that
coordinate distributions of the elements M44,24,34(m×n) of Mueller matrix of skin derma tissue of all types is characterized by individual (according to quantitative and topological structure) networks of characteristic points (Figures 4–6);
total amount of ±C-points (M44*(m×n)=0) sequentially increases for the samples of skin derma of “A”, “B”, “C” types (Figures 4(a)–6(a));
dependencies N0,90(x) of the number of characteristic values of matrix elements (8) and (9) for the samples of skin derma tissue of all types are similar in their structure (Figures 7(d), 7(e) and 7(f));
distributions N45,135(x) for the samples of skin derma tissue of “B”-type are characterized by sufficient increase (by 2-3 times) of the number of characteristic values in comparison with similar dependencies N0,90(x) (Figures 7(e) and 7(h)).
The obtained results can be connected with the increase of birefringence (Δn≈1.5×10-2) of collagen fibrils of pathologically changes skin derma of “B”- and “C”-types. Besides, at early stages (precancer) the directions of the growth of newly formed fibrils are being formed. At cancer states such pathologically changed fibrils form specifically oriented network of biological crystals.
In terms of physics, such morphological processes are manifested in the increase of probability of forming the ±C-points (skin derma samples of “B” and “C”-types), as well as in appearance of asymmetry between ranges of dependences values N0,90(x) and N45,135(x), which characterize the number of orthogonal L- and C-points.
In the end, the comparative investigations of diagnostic efficiency of the potential of famous techniques of laser polarimetry (Z1,2,3,4(M44,34,24(m×n))) [27]; polarization-correlation mapping (Z1,2,3,4(V4(m×n)={0,±1)) [24, 25] and the technique of Mueller-matrix singular diagnostics Z1,2,3,4(N(x),N0,90(x),N45,135(x)) were suggested.
Table 1presents statistical averaged values within the three groups of samples of myometrium tissue (Z1,2,3,4(M44,34,24(m×n))); (Z1,2,3,4(V4(m×n)={0,±1)) and Z1,2,3,4(N(x),N0,90(x),N45;135(x)).
Values (Z1,2,3,4(M44,34,24(m×n))), (Z1,2,3,4(V4(m×n)=0,or±1)) and Z1,2,3,4(N(x),N0,90(x),N45;135(x)) statistically averaged within the three groups of skin derma samples.
Zj=1,2,3,4
skin derma “A”-type (21 samples)
skin derma “B”-type (21 samples)
skin derma “C”-type (21 samples)
M44
M34
M24
M44
M34
M24
M44
M34
M24
Z1
0,67 ± 0.059
0,32 ± 0.031
0,27 ± 0.02
0,59 ± 0.048
0,27 ± 0.02
0,24 ± 0.05
0,37 ±0.034
0,19 ± 0.015
0,18 ±0.019
Z2
0,51 ±0,046
0,29 ± 0,019
0,26 ± 0,021
0,57 ±0,05
0,23 ± 0,013
0,21 ± 0,01
0,28± 0,023
0,21 ± 0,018
0,17 ± 0,015
Z3
1,13 ± 0,12
0,88 ± 0,09
0,74 ± 0,08
0,98 ± 0,1
0,79± 0,081
0,66 ± 0,069
0,66 ± 0,071
0,49 ± 0,042
0,41 ± 0,04
Z4
3,15 ± 0,32
2,11 ± 0,29
2,27 ± 0,31
2,84 ± 0,24
1,79 ± 0,18
1,87 ± 0,17
1,57 ± 0,14
1,07 ± 0,1
1,12 ± 0,11
Zj=1,2,3,4
V4=0
V4=±1
V4=0
V4=±1
V4=0
V4=±1
Z1
0.12±0,079
0.24±0,038
0.15±0,067
0.28±0,068
0.18±0,071
0.29±0,076
Z2
0.16±0,074
0.31±0,042
0.19±0,031
0.38±0,052
0.23±0,019
0.4±0,048
Z3
0.70±0,052
0.92±0,086
0.93±0,094
1.12±0,101
1.27±0,112
1.72±0,123
Z4
1.71±0,13
2.19±0,18
2.01±0,19
3.13±0,27
3.41±0,31
4.01±0,31
Zj=1,2,3,4
N(x)
N0,90(x)
N45,135(x)
N(x)
N0,90(x)
N45,135(x)
N(x)
N0,90(x)
N45,135(x)
Z1
0,61 ± 0,052
0,43 ± 0,038
0,12 ± 0,034
0,39 ± 0,05
0,37 ± 0,042
0,23 ± 0,048
0,29 ± 0,039
0,28 ± 0,036
0,19 ± 0,038
Z2
0,75 ± 0,068
0,82 ± 0,076
0,15 ± 0,021
0,36 ± 0,042
0,69 ± 0,056
0,46 ± 0,02
0,31 ± 0,024
0,55 ± 0,049
0,23 ± 0,026
Z3
1,19 ± 0,15
0,92 ± 0,01
1,86± 0,19
0,63 ± 0,051
0,87 ± 0,07
1,86 ± 0,19
0,53 ± 0,041
0,76 ± 0,062
2,06± 0,21
Z4
1,99 ± 0,17
2,31 ± 0,19
2,32 ± 0,21
1,01 ± 0,1
2,07 ± 0,17
8,45 ± 0,73
0,8 ± 0,07
1,87 ± 0,19
2,91 ± 0,32
It follows from the data presented that:
efficiency of laser polarimetry for diagnostics and differentiation of early oncological changes of skin derma tissue is insufficient—the difference between the values of statistical moments (Z1,2,3,4(M44,34,24(m×n))) of samples “A”, “B” and “C”-types is insufficient and does not exceed 20%–45%;
the technique of polarization-correlation mapping is efficient for differentiation of optical properties of sound and oncologically changed skin derma tissue—skewness (Z3) and kurtosis (Z4) of distribution of the number of singular points of “A”- and “B”-types of laser images differ by 1.53 and 2.15 times;
the technique of Mueller-matrix singular diagnostics is efficient for differentiation of optical properties of all types of samples—statistical moments of the 3rd and 4th orders of distributions N(x) for the samples “A”, “B” and “C”-types differ by 1.7 and 2.5 times respectively;
for distributions N45,135(x) of the amount of orthogonal singular L45,135- and C45,135-points of skin derma tissue of “A” and “B” types the maximal difference (from 2.2 to 4.1 times) is observed between all statistical Zj=1,2,3,4.
5. Conclusions
Correlation between the coordinate locations of characteristic points of 2D elements of Mueller matrix of optically thin layer of biological tissue and the network of L- and C-points in its laser image is defined. The potentiality of Mueller-matrix sampling of polarization singularities formed by biological crystals with orthogonally oriented (ρ=0∘, 90∘ and ρ=45∘, 135∘) optical axes is shown. The efficiency of Mueller-matrix diagnostics not only for oncological changes of skin derma tissue but also for differentiating their severity degree is demonstrated.
CowinS. C.How is a tissue built?200012265535692-s2.0-003452652310.1115/1.1324665AngelskyO. V.UshenkoA. G.UshenkoYu. A.UshenkoYe. G.TomkaYu. Ya.PishakV. P.Polarization-correlation mapping of biological tissue coherent images20051060640252-s2.0-3364522982610.1117/1.2148251de BoerJ. F.MilnerT. E.Review of polarization sensitive optical coherence tomography and Stokes vector determination2002733593712-s2.0-003662948810.1117/1.1483879de BoerJ. F.MilnerT. E.NelsonJ. S.Trends in optics and photonics (TOPS)1998Washington, DC, USAOSAEverettM. J.SchoenenbergerK.ColstonB. W.Jr.Da SilvaL. B.Birefringence characterization of biological tissue by use of optical coherence tomography19982332282302-s2.0-0000862635AngelskyO. V.TomkaYu. Ya.UshenkoA. G.UshenkoYe. G.YermolenkoS. B.UshenkoYu. A.2-D tomography of biotissues images in pre-clinic diagnostics of their pre-cancer states5972Advanced Topics in Optoelectronics, Microelectronics, and Nanotechnologies II2004158162Proceedings of the SPIE2-s2.0-3324447461410.1117/12.639726AngelskyO. V.UshenkoA. G.BurkovetsD. N.UshenkoYu. A.Polarization visualization and selection of biotissue image two-layer scattering medium20051010140102-s2.0-1814439927410.1117/1.1854674AngelskyO. V.TomkaYu. Ya.UshenkoA. G.UshenkoYe. G.UshenkoYu. A.Investigation of 2D Mueller matrix structure of biological tissues for pre-clinical diagnostics of their pathological states20053823422742352-s2.0-2784445942010.1088/0022-3727/38/23/014JiaoS.YuW.StoicaG.WangL. V.Optical-fiber-based Mueller optical coherence tomography20032814120612082-s2.0-0042703902UshenkoYu. A.Statistical structure of polarization-inhomogeneous images of biotissues with different morphological structures200566370Angel'skiiO. V.UshenkoA. G.ArkhelyukA. D.ErmolenkoS. B.BurkovetsD. N.UshenkoYu. A.Laser polarimetry of pathological changes in biotissues20008969739782-s2.0-014167331010.1134/1.1335052AngelskyO. V.UshenkoA. G.UshenkoYu. A.Polarization reconstruction of orientation structure of biological tissues birefringent architectonic nets by using their Mueller-matrix speckle-images200527279NyeJ. F.1999Bristol, UKInstitute of PhysicsNyeJ. F.Polarization effects in the diffraction of electromagnetic waves: the role of disclinations198338717921051322-s2.0-0021098489AngelskyO. V.MaksimyakA. P.MaksimyakP. P.HansonS. G.Interference diagnostics of white-light vortices20051320817981832-s2.0-3014444014110.1364/OPEX.13.008179NyeJ. F.HajnalJ. V.The wave structure of monochromatic electromagnetic radiation1987409183621362-s2.0-0023647750BerryM. V.DennisM. R.Polarization singularities in isotropic random vector waves200145720051411552-s2.0-000280109210.1098/rspa.2000.0660DolgovA. D.DoroshkevichA. G.NovikovD. I.NovikovI. D.Classification of singular points in the polarization field of the cosmic microwave background and eigenvectors of the Stokes matrix19996964274332-s2.0-0038957582KonukhovA. I.MelnikovL. A.Optical vortices in a vector field: the general definition based on the analogy with topological solitons in a 2D ferromagnet, and examples from the polarization transverse patterns in a laser200132S139S1442-s2.0-003530568810.1088/1464-4266/3/2/358AngelskyO. V.MokhunI. I.MokhunA. I.SoskinM. S.Interferometric methods in diagnostics of polarization singularities200265350366022-s2.0-4134911352410.1103/PhysRevE.65.036602NyeJ. F.Lines of circular polarization in electromagnetic wave fields198338917972792902-s2.0-0021097922AngelskyO. V.BurkovetsD. N.MaksimyakP. P.HansonS. G.Applicability of the singular-optics concept for diagnostics of random and fractal rough surfaces20034222452945402-s2.0-0042381931AngelskyO. V.UshenkoA. G.UshenkoYe. G.2-D Stokes polarimetry of biospeckle tissues images in pre-clinic diagnostics of their pre-cancer states200522633AngelskyO. V.UshenkoA. G.UshenkoYu. A.UshenkoYe. G.Polarization singularities of the object field of skin surface20063916354735582-s2.0-3374692166210.1088/0022-3727/39/16/005AngelskyO. V.UshenkoA. G.UshenkoYe. G.TomkaY. Y.Polarization singularities of biological tissues images20061150540302-s2.0-3384535253210.1117/1.2360527AngelskyO. V.HansonS. G.MaksimyakA. P.MaksimyakP. P.On the feasibility for determining the amplitude zeroes in polychromatic fields20051312439644052-s2.0-2274445444310.1364/OPEX.13.004396UshenkoA. G.PishakV. P.2004Dordrecht, The NetherlandsKluwer Academic Publishers