Research Article | Open Access

A. G. Ushenko, I. Z. Misevich, V. Istratiy, I. Bachyns'ka, A. P. Peresunko, Omar Kamal Numan, 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", *Advances in Optical Technologies*, vol. 2010, Article ID 423145, 8 pages, 2010. https://doi.org/10.1155/2010/423145

# Evolution of Statistic Moments of 2D-Distributions of Biological Liquid Crystal Net Mueller Matrix Elements in the Process of Their Birefringent Structure Changes

**Academic Editor:**Oleg V. Angelsky

#### Abstract

This research is aimed to investigate the reliability of Mueller-matrix differentiation of birefringence change of optically thick layers of biological liquid crystals at the early stages of the change in their physiological state. This is performed by measuring the set of skewness and kurtosis values of Mueller matrix image of the phase element in various points of the object under investigation.

#### 1. Introduction

Among many methods of optical diagnostics of organic phase-inhomogeneous object a new technique—laser polarimetry [1]—has been formulated within recent 10 years. It enables to obtain information about optical anisotropy [2–5] of biological tissues (BT) in the form of coordinate distributions of BT Mueller matrix elements, azimuths, and ellipticities of their object field polarization.

To analyze this polarimetric information the following model approach was elaborated [1, 2, 6–13]:(i)all the variety of human BT can be represented as four main types—connective, muscle, epithelial, and nerve tissues;(ii)morphological structure of any BT type is regarded as a 2-component amorphous-crystalline structure; (iii)the crystalline component or extracellular matrix is an architectonic net consisting of coaxial cylindrical protein (collagen, myosin, elastin, etc.) fibrils;(iv)optically, the protein fibrils possess the properties of uniaxial birefringent crystals;(v)interaction of laser radiation with the BT layer is considered in the single scattering approximation, when the attenuation factor corresponds to .

Specifically, the above mentioned model was used for finding and substantiating the interconnections between the ensemble of statistic moments of the 1st–4th orders that characterize the orientation-phase structure (distribution of optical axes and phase shifts of protein fibrils networks directions) of birefringent BT architectonics and that of 2D distributions of the elements of the corresponding Mueller matrix [6, 10]. It was determined [8, 11] that the 3rd and the 4th statistic moments of coordinate distributions of “phase” matrix elements (, , ) are the most sensitive to the change (dystrophic and oncological processes) of optical anisotropy of protein crystals. These statistic moments characterize the BT extracellular matrix birefringence. On this basis the criteria of early diagnostics of muscle dystrophy, precancer states of connective tissue, collagenosis, and so forth were determined.

However, such techniques do not take into account the coordinate heterogeneity of orientation-phase structure of protein crystals nets of the BT layer as well as the order of scattering in its depth. Thus it is important to investigate the distribution of statistic moments of the 1st–4th orders that characterize the 2D elements of Mueller matrix not only in the section of the probing laser beam but also within the whole BT layer of various optical thickness and physiological state.

For this we shall consider the potentiality of matrix modeling of polarization properties of optically thick BT layer.

#### 2. Mueller Matrix Modeling of the Properties of Optically Thick Biological Tissue

Let us represent the layer of such a biological object as the set of successively located optical thin partial layers (Figure 1).

Polarization properties of birefringent nets of every other BT layer are described by the Mueller matrix representing a superposition of matrix operators of optically separate coaxial protein fibrils [1]: where Here is the direction of optical axis determined by the direction of packing of the birefringent fibril, is the phase shift introduced between the orthogonal components of the amplitude of laser wave with the length passing through the fibril with linear size of its geometrical section and birefringence index .

Mueller matrix elements of the net of protein fibrils of partial BT layer () are determined by the following algorithms: (i)for finite number () of fibrils (ii)for “infinite” () number of fibrils

Here and are the distribution function of orientation () and phase () parameters of biological crystals network.

Mueller matrix of optically thick BT is determined by multiplication of partial matrix operators:

To make it simpler (without decreasing the analysis depth) further we shall consider only a 2-layered BT:

In the expanded form the matrix operator elements (6) are written as follows:

The analysis of relations (7) shows that optical properties of anisotropic component of 2-layered BT () are described by complex superposition of orientation and phase parameters of its partial layers' biological crystals . The concrete form of such dependences for “phase” matrix elements [10] is illustrated by the following relations:

The analysis of relations (8) shows that a strict solution of the inverse problem , revealing the changes in the structure of biological crystals network of one of the layers on the basis of data about matrix elements and , is both mathematically incorrect and physically ambiguous.

Thus, it is important to elaborate approximated statistical methods of experimental solution of such diagnostic task.

Generally speaking, the increase of optical thickness of biological layer results in multiple scattering regime in its volume. Such process appeared in formation of integral depolarization of laser field and results in averaging within angular aperture of photodetector of azimuth and ellipticity random values of fully polarized isolated speckles [8]. In this case the direct interconnection between the scattered field polarization parameters and the orientation and phase structure of biological crystals will be lost. From the other hand, for such integrally depolarized fields the Mueller matrix approach remains adequate [1]. Therefore, the search of interconnections “object-laser field” on the level of indirect, statistic parameters, which characterizes their structure is topical.

#### 3. Scheme of Experimental Changes of Coordinate Distributions of Biological Tissue Mueller Matrix Elements

The histological sections of different geometrical (optical) thickness were used as the objects of investigation. The technique of obtaining of such objects is convenient: biological tissue is freezing to nitrogen temperature with the following obtaining, by means of medical microtome, the histological sections (from 10 to 100 ). Such a technique allows preserving the optical properties of biological tissues over a period of 24 hours.

Conventional optical scheme of polarimeter for measuring 2D distributions of the BT Mueller matrix elements is presented in Figure 2 [1].

Histological sections of BT were illuminated by a parallel beam of He-Ne laser (*λ* = 0.6328 , = 5.0 ) with the radius mm. Polarization illuminator consists of quarter-wave plates 3, 5 and polarizer 4, providing the formation of laser beam with random azimuth or ellipticity of polarization. Polarization images of BT by means of microobjective 7 were projected into the plane of sensitized plate ( pixels) if with CCD-camera 10. The analysis of BT images was carried out by means of polarizer 9 and quarter-wave plate 8. As a result, the Stokes vector parameters for every pixel of the BT image were determined and the set of elements of Mueller matrix was calculated according to algorithm [10]:
Indices 1–4 correspond to the following polarization states of the beam probing the BT layer: 1–; 2–; 3–; 4–*⊗* (right-hand circulation).

The technique of measuring the ensemble (—the number of probing areas) of 2D matrix elements of BT consisted in the following sequence of actions: (i)the plane of BT histological sections (10 mm 20 mm) BT was scanned line by line by the laser beam with the radius with the step of linear displacement ;(ii)within each of 50 areas ( = ) of illumination of the BT layer plane according to the algorithm (9) the local array () of values of Mueller matrix elements was determined: (iii)for every local array statistic moments of the 3rd-4th orders were calculated; (iv)the histograms of the values of statistical moments of higher orders within the whole plane of the investigated BT histological sections were determined.

#### 4. Investigation of Skewness and Kurtosis Distributions of 2D Phase Element of Mueller Matrix of the Samples of Skeletal Muscle Tissue

The following histological sections of different types of rat’s tissues were investigated—connection tissue, muscular tissue, epithelial tissue, and nervous tissue. All these tissues have common optical peculiarity—the presence of birefringent network of protein fibrils [1]. Therefore, in this work, without decrease of analysis completeness, we have concentrated on investigation of optical anisotropy of rat's skeletal muscle tissue birefringent fibrils.

Distributions of the values of skewness and kurtosis of 2D phase () element of Mueller matrix of rat's skeletal muscle tissue were investigated. Histological sections of various optical thickness ( and ) and physiological state were the object of investigation. Such choice of optical thickness of histological sections () allows us to investigate the influence of a single () and multiple () scattering on coordinate structure of Mueller matrix elements.

The choice of matrix element as an analytical parameter is explained by the fact that it is the most sensitive to the changes of birefringence of protein fibrils nets, connected with their pathological changes [1, 11]. The choice of a skeletal muscle tissue of a rat as the object of investigation is connected with the possibility of direct experimental formation of septic inflammation and monitoring the control of their optical manifestations under the condition of scattering of laser beam of various multiplication factor.

The series of Figures 3 and 4 presents the distributions of the element of optically thin histological sections of healthy (Figure 3) and inflamed (Figure 4) muscle tissue.

**(a)**

**(b)**

**(a)**

**(b)**

The following can be seen from the data obtained: (i)2D distributions of matrix elements of both types of tissues are coordinately heterogeneous;(ii)for the sample of healthy muscle tissue (Figure 3) the lesser range of values change (from −0.7 to 0.7) of the matrix element is typical in comparison with the changes (from −0.85 to 0.85) of the analogous element of Mueller matrix of the inflamed tissue (Figure 4).

The structure of 2D distributions of elements can be connected with the following peculiarities of birefringent architectonics of the skeletal tissue samples. Firstly, while a laser wave (with the wave-length ) propagates through the network of optically anisotropic (birefringence index [7]) myosin fibrils, a wide range of values of phase shifts , proportional to their geometrical sizes , is formed

Thus, phase elements distributions are coordinately heterogeneous and dependent on the peculiarities of morphological structure of extracellular matrix of muscle tissue samples.

Secondly, the inflammation process is accompanied by the increase of birefringence due to formation of myosin fibrils edema [7]. Optically, it causes the increase of the range of phase shifts changes and the growth of fluctuations of the element connected with it (Figure 4).

The histograms of distribution of statistic moments of higher-orders of coordinate distributions of investigated samples of optically thin histological sections are presented in Figure 5.

**(a)**

**(b)**

**(c)**

**(d)**

The analysis of experimentally measured histograms proved the following:(i)the change ranges of the values of skewness and kurtosis of coordinate distributions , measured for the healthy and inflamed muscle tissues, do not actually coincide; (ii)extreme values of skewness in the distribution for phase element of the healthy tissue are by 4–6 times higher than analogous values of the given statistic moment for the inflamed muscle tissue; (iii)extreme values of kurtosis in the distribution for the inflamed tissue are by 2-3 times higher than the values of the statistic moment for the healthy skeletal tissue.

Thus, it can be stated that for single differentiation of physiological state of optically thin layers of muscle tissue it is enough to measure the 2D phase matrix element in one domain () of irradiation by the laser beam and to calculate the skewness and kurtosis of its values distribution.

The series of Figures 6 and 7 presents the results of experimental measurements of phase elements of Mueller matrix of optically thick () of the rat's muscle tissue layers.

**(a)**

**(b)**

**(a)**

**(b)**

It can be seen from the data obtained the following: (i)for the samples of the healthy and inflamed muscle tissue the same range of values change (from −0,55 to 0,55) is typical in the distribution of matrix elements ;(ii)comparative visual analysis of coordinate distributions of matrix elements of histological sections of skeletal tissue of both types did not show any sufficient difference between them.

Similarity of the distribution structure of both types of rat's skeletal tissue samples can be explained by multiple light-scattering. As a result of every local (th) act of interaction between laser radiation and separate fibril random value of phase shift is formed, which is multiplied () in the process of propagation in the BT depth reaching equiprobable random values from 0 to , .

Figure 8 shows histograms of the set of values of statistic moments of skewness and kurtosis of distributions of the matrix element of optically thick samples of muscle tissue.

**(a)**

**(b)**

**(c)**

**(d)**

The analysis of experimental data showed the following: (i)the skewness values distribution of 2D phase , and changes within the whole plane of samples of the layer of the healthy and inflamed muscle tissue practically coincide and cannot be used as an objective criterion for their optical properties differentiation;(ii)histograms of the values of the statistic moment of elements distribution possess individual structure that depends on physiological state of muscle tissue;(iii)extreme values of kurtosis in the histogram of the healthy tissue (c) are by 4-5 times less than analogous values of the given statistic moment of the inflamed rat's skeletal tissue (d).

#### 5. Conclusion

Thus for realization of Mueller-matrix differentiation of optically thick BT layers at early stages of changing their physiological state the measurement of distributions of the kurtosis of the ensemble of phase elements is objective.

The Mueller-matrix approach in diagnostics and differentiation of the biological crystals networks appears to be effective for optically thick histological sections layers of human tissues (skin derma, tissues of women reproductive sphere, parenchymatous tissues) also.

#### References

- A. G. Ushenko and V. P. Pishak,
*Coherent-Domain Optical Methods. Biomedical Diagnostics, Environmental and Material Science*, Kluwer Academic Publishers, Dordrecht, The Netherlands, 2004. - J. F. de Boer and T. E. Milner, “Review of polarization sensitive optical coherence tomography and Stokes vector determination,”
*Journal of Biomedical Optics*, vol. 7, no. 3, pp. 359–371, 2002. View at: Publisher Site | 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, DC, USA, 1998. - 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,”
*Optics Letters*, vol. 23, no. 3, pp. 228–230, 1998. View at: Google Scholar - S. Jiao, W. Yu, G. Stoica, and L. V. Wang, “Optical-fiber-based Mueller optical coherence tomography,”
*Optics Letters*, vol. 28, no. 14, pp. 1206–1208, 2003. View at: Google Scholar - O. V. Angelsky, G. V. Demianovsky, A. G. Ushenko, D. N. Burkovets, and Yu. A. Ushenko, “Wavelet analysis of two-dimensional birefringence images of architectonics in biotissues for diagnosing pathological changes,”
*Journal of Biomedical Optics*, vol. 9, no. 4, pp. 679–690, 2004. View at: Publisher Site | Google Scholar - O. V. Angelsky, Yu. Ya. Tomka, A. G. Ushenko, Ye. G. Ushenko, S. B. Yermolenko, and Yu. A. Ushenko, “2-D tomography of biotissues images in pre-clinic diagnostics of their pre-cancer states,” in
*Advanced Topics in Optoelectronics, Microelectronics, and Nanotechnologies II*, Proceedings of the SPIE, 2004. View at: Publisher Site | Google Scholar - O. V. Angelsky, A. G. Ushenko, D. N. Burkovets, and Yu. A. Ushenko, “Polarization visualization and selection of biotissue image two-layer scattering medium,”
*Journal of Biomedical Optics*, vol. 10, no. 1, Article ID 014010, 2005. View at: Publisher Site | Google Scholar - O. V. Angelsky, Yu. Ya. Tomka, A. G. Ushenko, Ye. G. Ushenko, and Yu. A. Ushenko, “Investigation of 2D Mueller matrix structure of biological tissues for pre-clinical diagnostics of their pathological states,”
*Journal of Physics D*, vol. 38, no. 23, pp. 4227–4235, 2005. View at: Publisher Site | Google Scholar - Yu. A. Ushenko, “Statistical structure of polarization-inhomogeneous images of biotissues with different morphological structures,”
*Ukrainian Journal of Physical Optics*, vol. 6, pp. 63–70, 2005. View at: Google Scholar - O. V. Angelsky, A. G. Ushenko, and Yu. A. Ushenko, “Polarization reconstruction of orientation structure of biological tissues birefringent architectonic nets by using their Mueller-matrix speckle-images,”
*Journal of Holography and Speckle*, vol. 2, pp. 72–79, 2005. View at: Google Scholar - O. V. Angelsky, P. P. Maksimyak, V. V. Ryukhtin, and S. G. Hanson, “New feasibilities for characterizing rough surfaces by optical-correlation techniques,”
*Applied Optics*, vol. 40, no. 31, pp. 5693–5707, 2001. View at: Google Scholar - O. V. Angelsky, P. P. Maksimyak, and T. O. Perun, “Optical correlation method for measuring spatial complexity in optical fields,”
*Optics Letters*, vol. 18, no. 2, pp. 90–92, 1993. View at: Google Scholar

#### Copyright

Copyright © 2010 A. G. Ushenko et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.