International Journal of Optics

Volume 2013 (2013), Article ID 683174, 7 pages

http://dx.doi.org/10.1155/2013/683174

Research Article

Spatial-Frequency Azimuthally Stable Cartography of Biological Polycrystalline Networks

V. A. Ushenko, N. D. Pavlyukovich, and L. Trifonyuk

Correlation Optics Department, Chernivtsi National University, 2 Kotsyubinsky Street, Chernivtsi 58012, Ukraine

Received 1 May 2013; Revised 15 July 2013; Accepted 1 August 2013

Academic Editor: Nicusor Iftimia

Copyright © 2013 V. A. 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.

Abstract.

A new azimuthally stable polarimetric technique processing microscopic images of optically anisotropic structures of biological tissues histological sections is proposed. It has been used as a generalized model of phase anisotropy definition of biological tissues by using superposition of Mueller matrices of linear birefringence and optical activity. The matrix element has been chosen as the main information parameter, whose value is independent of the rotation angle of both sample and probing beam polarization plane. For the first time, the technique of concerted spatial-frequency filtration has been used in order to separate the manifestation of linear birefringence and optical activity. Thereupon, the method of azimuthally stable spatial-frequency cartography of biological tissues histological sections has been elaborated. As the analyzing tool, complex statistic, correlation, and fractal analysis of coordinate distributions of element has been performed. The possibility of using the biopsy of the uterine wall tissue in order to differentiate benign (fibromyoma) and malignant (adenocarcinoma) conditions has been estimated.

1. Introduction

The state of the art in medical science demands from the experts not only certain practical skills but also knowledge of general pathologic processes, their nomenclature and definitions, causes, pathogenesis, and results as well as their importance for human. The pathologic processes on a cellular and tissue level can be detected by means of light microscope.

Samples preparation for microscopic investigations is a complex process that involves degradation of image quality of histological samples. In order to overcome the problems connected with image deformation (artifacts), the high professionalism and complex chore are necessary for results correction. The usage of computing for the mentioned tasks allows facilitating the work with medical images.

The methods and algorithms of objects isolation in the microscopic images—segmentation of histological objects in grayscale and color images—are mainly developing at present time. The algorithms of segmentation of extended (vessels, fibers) and areal (cells and their structures) objects are being elaborated [1–8].

A new medical-physical method, namely, laser polarimetry of biological tissues histological sections, becomes an extension and an addition to the techniques of microscopic images structure processing [9–11]. The background of such method consists in measurements of coordinate distributions of polarization states (polarization maps) in the plane of microscopic images of biological tissues histological sections. Such method allows obtaining a new, inaccessible measurement for histological and mathematical methods of analysis and information about optical anisotropy (linear and circular birefringence [11]) of nonuniformly scaled structural elements of different biological objects. At the same time, a complex analysis of the polarization maps of a tissue specimen is azimuthally dependent on the probing beam polarization plane and sample rotation angle. Such a condition complicates using such a method during the comparative investigations of groups of histological sections with different pathologies.

Thus, further progress of laser polarimetry may be connected with the development of azimuthally stable methods of direct measurement of the parameters of linear and circular birefringences.

The solution of this task is the use of Mueller-matrix cartography with the so-called rotational invariants. The Mueller matrix elements appear to be azimuthally stable and sample rotation angle () independent [12–15].

It has been determined that, in the case of optically thin (attenuation coefficient ) layers, the value is preferably defined by phase anisotropy mechanisms—linear and circular birefringences [16].

The separation of such mechanisms is possible with the help of spatial-frequency Fourier selection [17].

Our research is aimed at designing the experimental method of azimuthally stable Fourier polarimetry and spatial-frequency selection of parameter distributions of linear and circular birefringences of the biopsy of the uterine wall tissue in order to differentiate benign (fibromyoma) and malignant (adenocarcinoma) conditions.

2. The Theory of Method

The following model concepts have been assumed [9–11] by us as the basis for an analysis of the processes of modulation of laser radiation by the polycrystalline network of biological crystals with linear and circular birefringences:

(i)the uterine wall tissue consists of two optical anisotropic layers of the muscular (myometrium) and connective (endometrium) layers;(ii)the optical anisotropic myometrium network is formed by large-scale (the range of the transverse dimensions ) filiform myosin fibrils with a prevalence by linear birefringence;(iii)the optical anisotropy of such structures is characterized by coordinate distribution of the matrix element in the plane of histological section, where is a direction of fibrils packing and is a phase shift between linearly polarized orthogonal components of the light beam amplitude;(iv)the optical anisotropic component of the endometrium is formed by collagen fibers regulated according to the direction of optic axes (; ) with predominantly circular birefringence , where is polarization plane rotation angle.

From the medical point of view, “optical separation” of the polarization manifestations of linear and circular birefringences of optical anisotropic layers of the uterine wall is the topical task.

The fact is that, at the initial, early stages, the pathological changes in the human organism are accompanied by a concentration increase of optically active proteins in case of the generalized preservation of the morphological macrostructures of polycrystalline protein network of the myometrium and endometrium [16, 17].

Fibrillar structural changes—the formation of the growth direction of tumors shoots and so on—correspond to a more developed pathological condition.

With the purpose of a selection of such conditions, we have applied the method of spatial-frequency filtration of the polarization inhomogeneous field of laser radiation in Fourier plane [18, 19].

The main idea of such approach is that spatial-frequency structure of Fourier form of laser image of the tissue of the uterine wall is different for its large-scale (myosin) and small-scale (collagen) protein structures. Proceeding from this, using spatial-frequency filtration, it is mainly possible to isolate either low-frequency (with linear birefringence) or high-frequency (with circular birefringence) components, which by means of Fourier conversion transform into, respectively, “separated” laser images.

3. Optical Realization of Spatial Frequency of Fourier Selection of the Uterine Tissue Optical Anisotropy of Protein Networks

Figure 1 presents a diagram of modified laser Fourier polarimeter with spatial-frequency filtration [20].

Figure 1: Optical scheme of Fourier polarimeter: 1—He-Ne laser, 2—a collimator, 3—a stationary quarter-wave plate, 5, 9—mechanically movable quarter-wave plates, 4, 10—a polarizer and analyzer, respectively, 6—object of study 7, 8—polarization microobjectives, 11—CCD camera, 12—a personal computer. —object plane, —image plane, and —Fourier plane with the diaphragms.

Illumination of a sample under study was performed by the parallel () laser beam of He-Ne ( μm,  mW). The polarization light source consisted of quarter-wave plates 3 and 5 and polarizer 4; it formed a right circularly polarized beam. The histological sections of the biopsy of the uterine wall tissue were placed in the focal plane of polarization microobjective 7 (focal distance—, aperture—0.1, magnification—). Behind the Fourier focal plane, the vignetting diaphragm was located; its size changed within the range of 2  to . Polarization microobjective 8 (focal distance— mm, aperture—0.1, and magnification—) was located at the focal length forming the frequency plane of lens 7 and, thus, performed inverse Fourier transform of a filtered image outside the polarization field of laser radiation. The coordinate distribution of intensity of such fields, polarizationally filtered by quarter-wave plate 9 and polarizer 10, was registered in the plane of CCD-camera 11 (the Imaging Source DMK 41AU02.AS, monochrome CCD, Sony ICX205AL (progressive scan); resolution—; light sensitive area size— μm; sensitivity—0.05 lx; dynamic range—8 bit; SNR—9 bit, deviation of photosensitive characteristics from linear no more than 15%). It provided the range of measuring the structural elements of polycrystalline network with the resolution of 2–2000 μm.

The matrix element was calculated for each probing beam within every pixel () as follows:

Here, and are the intensities of spatial-frequency filtered image for each polarization states of probing beams (; 0°; 90°).

In order to get a quantitative estimation of coordinate distributions , we have used, as in the case of traditional laser polarimetry [9–11, 16, 17], three mutually complementary techniques of mathematical analysis, namely, statistic, autocorrelation, and fractal.

We have calculated a set of statistical moments of the 1st–4th orders using the following algorithms:

The obtained data characterizes the histogram of the element of optically anisotropic structures in the plane of histological section.

The method of autocorrelation with the use of function was assumed as the basis of an analysis of the coordinate structure of distributions

Here, is the step of changing the coordinates () of a distribution of a set of polarization distributions for a specific horizontal line of the pixel digital camera. The resultant expression for the autocorrelation function is obtained by averaging expression (3) along all lines

The obtained dependency characterizes the coordinate interconnection between the optically anisotropic structures in the plane of histological section.

We have selected correlation moments and that determine a distribution and the excess of the autocorrelation function as quantitative parameters that characterize autocorrelation dependences .

A fractal analysis of distributions was performed using a calculation of logarithmic dependences for the power spectra with following evaluation of dispersion.

Dependences are approximated using the least-squares method into curves and are classified in such a way:

(i) are fractal on condition that the value of the slope is constant for 2-3 decades of a change of sizes ;(ii) are multifractal, on condition of the presence of several slopes in ;(iii) are random, in case of the absence of stable slopes over the whole range of a change sizes .

The obtained results characterize the scale self-similarity of the distribution of optically anisotropic structures in the plane of histological section.

4. Experimental Results and Discussion

As objects of investigation, we chose two groups of opticallythin (attenuation coefficient ), single scattered histological sections of the biopsy of benign (16 samples) and malignant (14 samples) tumors of the uterus.

The samples were prepared on a freezing microtome, using a standard technique. From the optical point of view, such samples characterized by transformation of polarization without depolarization.

The conventional microscopic images of the samples of both groups are presented in Figure 2.

(a)

(b)

Figure 2: Microscopic images of histological sections of the biopsy of benign (a) and malignant (b) tumors of the uterus.

As it can be seen from Figure 2, the coordinate and scale structures of such conventional microscopic images are similar. This fact considerably complicates the histological differentiation of benign and malignant states of uterus tissue.

4.1. Spatial-Frequency Fourier Polarimetry of Linear Birefringence of the Fibrillar Matrices of the Uterine Tissue

For the purpose of choosing the optimal conditions of spatial-frequency filtration, the following range of possible sizes of the vignetting diaphragm was used.

The criterion for diaphragm size choosing is a simultaneous change of the set of statistic moments (4) of . In our case, the optimal size was .

Such geometric size was chosen for comparative investigations of optical anisotropy of the fibrillar networks of the uterine tissue with benign and malignant changes, which characterize coordinate distributions of rotational invariant —Figure 3.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

Figure 3: “Low-frequency” coordinate maps (a), (e) histograms (b), (f) autocorrelation functions (c), (g) logarithmic dependence of the power spectra (d), (h) of the linear birefringence distribution of a histological section of benign uterine tumor (fibromyoma) biopsy ((a)–(d)), and malignant uterine tumor (adenocarcinoma) biopsy ((e)–(h)).

A comparative analysis of aggregate parameters that characterize matrix maps of linear birefringence of large-scale optically anisotropic fibrillar networks of the uterine tissue samples revealed some discrepancies between them. Namely, the main extrema of histograms of the distributions of random matrix maps, values for histological sections of both types are localized in different areas. Thus, the most probable value is for a benign tumor (Figure 3(b)), for a malignant one it is (Figure 3(f)). The revealed peculiarity, in our opinion, is connected with a more developed structure of newly formed fibers of malignant tumor than with benign fibromyoma of the similar optical anisotropic network.

Dependences of autocorrelation functions (4) of distribution of both uterine tissue samples smoothly and monotonously, (Figure 3, fragments (c) and (g)) indicate a coordinately homogeneous structure of the corresponding matrix maps formed by the mechanisms of linear birefringence of myosin fibrils. Besid such maps, a dimensionally self-similar construction has been deleted—logarithmic dependences of power distributions spectra are characterized by a similar slope (Figure 3, fragments (d), (h)).

Thus, matrix maps (Figure 3, fragments (a) and (e)) of linear birefringence of the polycrystalline network of myosin fibrils are characterized by a sufficiently similar correlation and homogeneous (fractal) structure. The main indication of oncologic changes is the formation of a developed optically anisotropic network with an increased level of linear birefringence index and a transformation of a histogram of the distribution of .

Quantitative differences between matrix maps of the protein matrices of the uterine tissue of both types illustrate average values and standard deviations of the set of statistic , correlative (), and fractal () parameters, illustrated in Table 1. Also, let us introduce the specificity parameter (—the number of the correct and incorrect diagnoses within a group) in differentiation of benign and malignant uterine tumors.

Table 1: Parameters of the statistic, correlation, and self-similar structure of of linear birefringence of polycrystalline uterine tissue networks.

Parameters

Fibromyoma
Adenocarcinoma

0.33 ± 0.059
0.72 ± 0.11
72%

0.14 ± 0.022
0.18 ± 0.031
61%

2.16 ± 0.37
0.83 ± 0.14
75%

1.95 ± 0.21
0.37 ± 0.051
77%

0.31 ± 0.049
0.33 ± 0.054
56%

0.035 ± 0.005
0.04 ± 0.006
54%

()
Fractal
Fractal
51%

0.33 ± 0.055
0.28 ± 0.051
52%

A comparative analysis of the data of laser-frequency Fourier-polarimetry of linear birefringence of protein polycrystalline networks of histological section of benign and malignant uterine tumors biopsy revealed sensibility to a differentiation of such states of all statistic moments of the 1st through the 4th orders with specificity level Sp ~ 60%–77% (printed in bold in Table 1).

4.2. Spatial-Frequency Fourier Polarimetry of Circular Birefringence of the Fibrillar Networks of Uterine Tissue

Diagnostic possibilities of differentiating of pathologically changed samples of the uterine tissue, using the method of high-frequency spatial-frequency Fourier-domain—polarimetry of circular birefringence employing opaque filter illustrates probable (fragments (b), (f)), autocorrelation (fragments (c), (g)), and spectral (fragments (d), (h)) dependences of matrix maps (fragments (a), (e)), adduced in a series in Figure 4.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

Figure 4: “High-frequency” coordinate maps (a), (e) histograms (b), (f) autocorrelation functions (c), (g) logarithmic dependence of the power spectra (d), (h) of the linear birefringence distribution of a histological section of benign uterine tumor (fibromyoma) biopsy ((a)–(d)), and malignant uterine tumor (adenocarcinoma) biopsy ((e)–(h)).

On comparing the findings about the structure of the distributions of circular birefringence index (Figure 4, fragments (a) and (e)) under the conditions of high frequency filtration of laser field radiation, transformed by polycrystalline collagen fibrils of a layer of the endometrial tissue are the following characteristics of the oncologic state, detected. Firstly, it is a range expansion of random values that changes for a histogram of the matrix map at the expense of a “pathologic” increase of the concentrations and the size of the collagenous fibrils (Figure 4, fragments (b), (f)). Secondly, the formation of growth orientations leads to the transformation of the multifractal (Figure 4, fragment (d)) distribution into a self-similar—inherent fixed inclination (Figure 4, fragment (h) for approximate curve ).

Quantitative differences between matrix maps of circular birefringence of polycrystalline collagenous networks of the endometrium with pathology of both types illustrate average values and standard deviations of statistic set , correlative (), and fractal () parameters, presented in Table 2.

Table 2: Parameters of the statistic, correlation, and self-similar structure of of circular birefringence of polycrystalline uterine tissue networks.

Parameters

Fibromyoma
Adenocarcinoma

0.13 ± 0.021
0.24 ± 0.037
75%

0.11 ± 0.018
0.22 ± 0.035
81%

0.74 ± 0.11
2.38 ± 0.39
87%

0.89 ± 0.13
1.58 ± 0.21
84%

0.19 ± 0.031
0.23 ± 0.034
53%

1.56 ± 0.28
0.76 ± 0.12
72%

()
Multifractal
Fractal
73%

0.31 ± 0.059
0.28 ± 0.047
51%

A comparative analysis revealed the following oncologic condition, with a high level of specificity () of uterus parameters (printed in bold in Table 2):

(i)statistic moments of the 1st and the 2nd orders ; differences between values are from 2 to 3, 5 times;(ii)correlation moments of the 4th order of the autocorrelation function of distributions ; differences between the values reach 2 times;(iii)transformation of multifractal distributions into fractal ones.

5. Conclusions

A model of generalized optic anisotropy of the fibrillar protein matrices has been suggested, and the method of Fourier polarimetry for the parameters of linear (phase shift between orthogonal components of the laser wave amplitude) and circular (the rotation angle of the polarization plane) birefringences of the tissue biopsy of the uterine wall a spatial frequency selection of such coordinate distributions for a differentiation of benign (fibromyoma) and malignant (adenocarcinoma) conditions has been substantiated.

A set of criteria for a polarization-phase differentiation of benign (fibromyoma) and malignant (adenocarcinoma) conditions of the uterine has been revealed and substantiated as follows:

(i)statistical moments of the 1st–4th order which characterize a distribution of azimuthally stable matrix element , stipulated by linear birefringence of myosin fibrils of the endometrium layer;(ii)statistical moments of the 1st–4th orders which characterize a distribution of azimuthally stable matrix element , stipulated by circular birefringence of the network of collagenous fibrils of the endometrium layer;(iii)correlation moments of the 4th order, which characterize autocorrelation functions of the distributions for , stipulated by circular birefringence of the network of collagenous fibrils; (iv)a slope change of approximate curve of log-log dependency for the power spectra of distribution of and transformed by optical anisotropic protein matrices.

References

K. Abe, F. Mizutani, and C. Wang, “Thinning of gray-scale images with combined sequential and parallel conditions for pixel removal,” IEEE Transactions on Systems, Man and Cybernetics, vol. 24, no. 2, pp. 294–299, 1994.

C. Arcelli and G. Ramella, “Sketching a grey-tone pattern from its distance transform,” Pattern Recognition, vol. 29, no. 12, pp. 2033–2045, 1996.

C. Arcelli and G. Ramella, “Finding grey-skeletons by iterated pixel removal,” Image and Vision Computing, vol. 13, no. 3, pp. 159–167, 1995.

A. Bieniek and A. Moga, “A connected component approach to the watershed segmentation,” in Mathematical Morphology and Applications to Image and Signal Processing, H. J. A. M. Heijmans and J. B. T. M. Roerdink, Eds., pp. 215–222, Kluwer Academic publishers, London, UK, 1998.

R. Boomgaard, Mathematical morphology: extension towards computer vision [Academisch proefschrift], Faculteit der Wiskunde en Informatica de Universiteit van Amsterdam, 1992.

L. E. Bracey Jr., “Image analysis for nanotechnology on the desktop,” Advanced Imaging, vol. 12, no. 11, pp. 54–55, 1997.

A. Nedzved and S. Ablameyko, “Extraction on cell structure in color histological images,” in Proceedings of the International Workshop on Computer Vision, Computer Graphics, New Media, pp. 41–49, Computer Gesellschaft, Graz, Austria, 2002.

A. M. Nedzved, S. V. Ablameyko, B. V. Nalibotsky, and Y. G. Iliych, “System for analysis of medical images of histological objects,” Pattern Recognition and Image Analysis, vol. 11, no. 4, pp. 732–742, 2001.

O. V. Angelsky, A. G. Ushenko, A. Yu. 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,” in Handbook of Photonics for Biomedical Science, V. V. Tuchin, Ed., pp. 283–322, CRC Press, Taylor & Francis Group, Boca Raton, Fla, USA, 2010.

O. V. Angelsky, V. P. Pishak, A. G. Ushenko, and A. Ushenko, “Statistical and fractal structure of biological tissue Mueller matrix images,” in Optical Correlation Techniques and Applications, O. Angelsky, Ed., vol. 6616 of Proceedings of SPIE, pp. 213–263, 2007.

Y. A. Ushenko, T. M. Boychuk, V. T. Bachynsky, and O. P. Mincer, “Diagnostics of structure and physiological state of birefringent biological tissues: statistical, correlation and topological approaches,” in Handbook of Coherent-Domain Optical Methods, pp. 107–148, Springer Science, New York, NY, USA, 2013.

J. M. Bueno and J. Jaronski, “Spatially resolved polarization properties for in vitro corneas,” Ophthalmic and Physiological Optics, vol. 21, no. 5, pp. 384–392, 2001.

M. Shribak and R. Oldenbourg, “Techniques for fast and sensitive measurements of two-dimensional birefringence distributions,” Applied Optics, vol. 42, no. 16, pp. 3009–3017, 2003.

M. H. Smith, “Interpreting Mueller matrix images of tissues,” in Laser-Tissue Interaction XII: Photochemical, Photothermal, and Photomechanical, vol. 4257 of Proceedings of SPIE, pp. 82–89, San Jose, Calif, USA, January 2001.

S.-Y. Lu and R. A. Chipman, “Interpretation of Mueller matrices based on polar decomposition,” Journal of the Optical Society of America A, vol. 13, no. 5, pp. 1106–1113, 1996.

O. V. Angelsky, Y. Y. Tomka, A. G. Ushenko, Y. G. Ushenko, and Y. 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.

A. G. Ushenko, D. N. Burkovets, and Y. A. Ushenko, “Polarization-phase mapping and reconstruction of biological tissue architectonics during diagnosis of pathological lesions,” Optics and Spectroscopy, vol. 93, no. 3, pp. 449–456, 2002.

J. W. Goodman, “Statistical properties of laser speckle patters,” in Laser Speckle and Related Phenomena, J. C. Dainty, Ed., pp. 9–75, Springer, Berlin, Germany, 1975.

M. Born and E. Wolf, Principles of Optics, Cambridge University Press, Cambridge, UK, 1999.

R. Józwicki, K. Patorski, O. V. Angelsky, A. G. Ushenko, D. N. Burkovets, and Y. A. Ushenko, “Automatic polarimetric system for early medical diagnosis by biotissue testing,” Optica Applicata, vol. 32, no. 4, pp. 603–612, 2003.