3D Mueller Matrix Reconstruction of the Optical Anisotropy Parameters of Myocardial Histopathology Tissue Samples

Diseases affecting myocardial tissues are currently a leading cause of death in developed nations. Fast and reliable techniques for analysing and understanding how tissues are affected by disease and respond to treatment are fundamental to combating the effects of heart disease. A 3D Mueller matrix method that reconstructs the linear and circular birefringence and dichroism parameters has been developed to image the biological structures in myocardial tissues. The required optical data is gathered using a Stokes polarimeter and then processed mathematically to recover the individual optical anisotropy parameters, expanding on existing 2D Mueller matrix implementations by combining with a digital holography approach. Changes in the different optical anisotropy parameters are rationalised with reference to the general tissue structure, such that the structures can be identified from the anisotropy distributions. The first to fourth order statistical moments characterising the distribution of the parameters of the optical anisotropy of the polycrystalline structure of the partially depolarising layer of tissues in different phase sections of their volumes are investigated and analysed. The third and fourth order statistical moments are found to be the most sensitive to changes in the phase and amplitude anisotropy. The possibility of forensic medical differentiation of death in cases of acute coronary insufficiency (ACI) and coronary heart disease (CHD) is considered as a diagnostic application. The optimal phase plane ( θ ∗ = 0.7 r a d ) has been found, in which excellent differentiation accuracy is achieved ACI and CHD - A c ( Δ Z 4 ( θ ∗ , Φ L , Δ L ) ) = 93.05 % ÷ 95.8 % . A comparative analysis of the accuracy of the Mueller-matrix reconstruction of the parameters of the optical anisotropy of the myocardium in different phase planes ( θ = 0.9 r a d and θ = 1.2 r a d ), as well as the 2D Mueller-matrix reconstruction method was carried out. This work demonstrates that a 3D Mueller matrix method can be used to effectively analyse the optical anisotropy parameters of myocardial tissues with potential for definitive diagnostics in forensic medicine.

Diseases affecting myocardial tissues are currently a leading cause of death in developed nations. Fast and reliable techniques for analysing and understanding how tissues are affected by disease and respond to treatment are fundamental to combating the effects of heart disease. A 3D Mueller matrix method that reconstructs the linear and circular birefringence and dichroism parameters has been developed to image the biological structures in myocardial tissues. The required optical data is gathered using a Stokes polarimeter and then processed mathematically to recover the individual optical anisotropy parameters, expanding on existing 2D Mueller matrix implementations by combining with a digital holography approach. Changes in the different optical anisotropy parameters are rationalised with reference to the general tissue structure, such that the structures can be identified from the anisotropy distributions. The first to fourth order statistical moments characterising the distribution of the parameters of the optical anisotropy of the polycrystalline structure of the partially depolarising layer of tissues in different phase sections of their volumes are investigated and analysed. The third and fourth order statistical moments are found to be the most sensitive to changes in the phase and amplitude anisotropy. The possibility of forensic medical differentiation of death in cases of acute coronary insufficiency (ACI) and coronary heart disease (CHD) is considered as a diagnostic application. The optimal phase plane (θ p 0.7rad) has been found, in which excellent differentiation accuracy is achieved ACI and CHD -Ac(ΔZ 4 (θ p , Φ L , Δ L )) 93.05%÷95.8%. A comparative analysis of the accuracy of the Mueller-matrix reconstruction of the parameters of the optical anisotropy of the myocardium in different phase planes (θ 0.9rad and θ 1.2rad), as well as the 2D Mueller-matrix reconstruction method was carried out. This work demonstrates that a 3D Mueller matrix method can be used to effectively analyse the optical anisotropy parameters of myocardial tissues with potential for definitive diagnostics in forensic medicine.

INTRODUCTION
Myocardial tissues form the muscles in the heart and are hence critical to human life [1][2][3]; diseases affecting myocardial tissues are currently the leading cause of death in developed nations [4]. Myocardial tissue sections provide an excellent platform for wider understanding of cardiac function [5], and hence imaging methods that can be applied to myocardial tissues are of significant interest. As exemplified by the recent rapid spread of COVID-19, new threats with potential to affect the heart are continually emerging [6]. To keep up, fast and reliable techniques for analysing and understanding how tissues are affected by disease, and respond to treatment, are fundamental necessities. 3D imaging of structures within tissues can provide such insights. Computed tomography scanning, magnetic resonance imaging, and X-rays are widely used. However, 3D imaging techniques can be limited by relatively low sample throughput, and high cost [7][8][9]. Furthermore, as medical diagnostics move away from human analysis to more automated, AI-driven approaches, different techniques for looking at tissues become more viable. In particular, statistical analysis of tissue properties can give rapid results with high degrees of accuracy [10,11]. One possibility of emerging interest is to look at the optical anisotropy of the tissue, from which one can then infer the tissue structure and other properties [12][13][14][15].
Optical anisotropy is a result of a material interacting differently with different polarisations of incident light, such that the different polarisations are absorbed, transmitted, reflected and refracted with different intensities [16]. There are four optical anisotropy properties that can be considered: linear birefringence, linear dichroism, circular birefringence, and circular dichroism. The birefringence and dichroism of a material can be determined by measuring changes in the polarisation of light passing through the material, by so-called polarimetry measurements [16][17][18][19]. Polarimetry is a relatively easy technique to implement -at its most simple, requiring only: a light source (laser), polarising filters (and quarter wave plate for circular birefringence and dichroism measurements), and a detector [20,21]. To achieve greater sensitivity of measurements, one can employ interferometric techniques.
Looking at the optical anisotropy of myocardial tissues is promising for two reasons: firstly, the tissues contain spatially ordered protein fibrils which should demonstrate clear linear anisotropy; secondly, the tissues are partially depolarising. After polarimetric measurements of the four anisotropy properties, it is necessary to correlate the experimental data with the physical sample under investigation. Mueller Matrix Polarimetric (MMP) diagnostics are the tool of choice for this purpose, with many distinct directions being considered, including: the investigation of scattering matrices [1,[22][23][24][25]; Mueller matrix polarimetry [17,[26][27][28][29]; polar decomposition of Mueller matrices [30,31]; and two-dimensional Mueller matrix mapping [18,19,32,33]. A Mueller matrix is a 4 × 4 matrix representing the effect of a specific optical element on the polarisation of light. The multiplication of the initial Stokes vector of light passing through the element by the Mueller matrix of the element gives the Stokes vector of the output light. MMP methods and tools have evolved around two limiting approximations. The first is the search for relationships between the angular indicatrices (1D), the coordinate Mueller matrix distributions (2D), and the structure of diffuse (depolarising) layers [1,[24][25][26][27]. The second is the MMP of optically thin, nondepolarising layers [17,28,29] with subsequent reconstruction of the distributions of the phase and amplitude anisotropy parameters [10,11,[34][35][36][37][38][39][40][41].
However, while such 2D methods maybe useful for imaging and characterising surfaces or materials which are isotropic in at least one direction, they are of limited use for more complex systems. It is generally desirable to understand variations across a structure fully in three dimensions. Hence, there is a need to expand MMP diagnostic techniques to three dimensions. Additionally, most biological objects are partially depolarising. They have spatially inhomogeneous, optically anisotropic structures. Therefore, it is necessary to further develop and generalise existing MMP techniques to consider such partially depolarising structures. The theoretical basis of this direction can be established from the synthesis of methods of differential Mueller matrix [12,13,[42][43][44][45] and holographic mapping of phase-inhomogeneous layers [46,47]. Here, we develop and experimentally demonstrate a technique for the 3D Mueller matrix reconstruction of the phase and amplitude anisotropy parameters. We first outline the method for Mueller matrix analysis in two dimensions, before demonstrating the further development to three dimensions. We then consider the functional possibility of 3D reconstruction of each of the four parameters of the optical anisotropy of a myocardial tissue layer. Finally, the approbation of the approach is carried out for definitive diagnosis of myocardium tissue death as a result of acute coronary insufficiency (ACI) and coronary heart disease (CHD).

2D Mueller Matrix Reconstruction
The 2D Mueller matrix reconstruction of the distributions of linear and circular birefringence and dichroism within biological layers has previously been considered in detail [48]. The matrix operator D characterises the distribution of the mean values of the parameters of the phase (Φ) and amplitude (Δ) anisotropy. The second-order differential matrix D determines the changes in the polarisation due to fluctuations of the linear and circular birefringence ( Φ) and dichroism ( Δ). The Mueller matrix resulting from the superposition of the first order (fully polarised part D) and second order (fully depolarised part D) differential matrices [42,43] gives algorithms for the interrelation between the distributions of phase and amplitude anisotropy and Mueller matrix images. Here we confine ourselves to the consideration of the completely polarised component of Frontiers in Physics | www.frontiersin.org October 2021 | Volume 9 | Article 737866 2 the Mueller matrix -the coordinate distributions of the elements of the first-order differential matrix d ik (x, y). The matrix operator D is defined as It contains six distinct, non-zero parameters: • Δ 0,90 , Δ 45,135the linear dichroism between the orthogonal components 0°− 90°and 45°− 135°, respectively; • Φ 0,90 , Φ 45,135the linear birefringence between the orthogonal components 0°− 90°and 45°− 135°, respectively; • Δ ⊗,⊕the circular dichroism for between right-(⊗) and left-(⊕) circularly polarised components.
These phase and amplitude anisotropy parameters are in turn determined by the following identities: Here: n j and τ j are the refractive indices and absorption coefficients for the j-polarised components (where j: 0°− 90°, 45°− 135°, and ⊗ − ⊕) of the incident laser radiation; λ is the wavelength of the laser radiation; and l is the thickness of the sample under investigation. For a partially depolarising layer, the following inter-relations between the individual elements of the first-order differential matrix d ik and the individual elements of the Mueller matrix F ik exist: By combining Eqs 2-7 with Eq. 8, we obtain algorithms determining the birefringence and dichroism from the elements of the Mueller matrix: These analytical expressions are the basis of Mueller matrix reconstruction of mean values of the optical anisotropy parameters of the layer. Without reducing the completeness of the analysis [42], we can reduce the system to consider only the generalised parameters of linear birefringence (Φ) and linear dichroism (Δ): Also noting the identities Eq 4, 7, we now have a complete set of the four anisotropy parameters Δ L , Δ C , Φ L , and Φ C . We note that the 2D Mueller matrix reconstruction obtains the averaged distributions of linear and circular birefringence and dichroism over the entire thickness l of the sample, without considering variations in the z direction. However, a complete analysis requires equal consideration of all three dimensions. Therefore, it is important to obtain 3D

3D Expansion
The key principle for the determination of a threedimensional layered series of distributions of Mueller matrix elements F ik (x, y, z) is the use of a reference laser radiation wave [48,49]. The six distinct polarization states are obtained in both the illuminating (Zond) and reference (Ref) beams: {Zond-Ref} > 0°; 90°; 45°; 135°; ⊗; ⊕ utilizing polarizer. The detection of partial interference patterns at the CCD camera (14, see Figure 1) through the polarizeranalyzer is defined with the orientation of the transmission plane at angles p 0°; 90°. A two-dimensional discrete Fourier transform DF(υ, ]) is applied to each partial interference distribution/image. The DF(υ, ]) of a two-dimensional Frontiers in Physics | www.frontiersin.org October 2021 | Volume 9 | Article 737866 3 array I P 0 0 ;90 0 (x, y) is a function of two discrete variables coordinates (x, y) is defined as [23]: where (DF p ) 0 0 ;90 0 ;45 0 ;135 0 ;⊗;⊕ P 0 0 ;90 0 (a, b) ≡ E 0 0 ;90 0 ;45 0 ;135 0 ;⊗;⊕ P 0 0 ;90 0 (a, b). Therefore, a distribution of complex amplitudes: is obtained for each state of polarization {Zond-Ref} in different phase planes θ i (δ 90 0 − δ 0 0 ). In single scattering approximation the phase of the fields of complex amplitudes Eq. 15 relates to physical depth z i in the volume of an optical anisotropy of biological layer While in multiple scattering the physical or effective depth z p i becomes proportional to the thickness of biological layer z The corresponding parameters of Stokes vector and polarization parameters of the object field for each phase plane θ i are defined as: Based on relations Eq 14, 15, 18, the elements of Mueller matrix {F} are calculated using the following Stokes-based polarimetric relation: is obtained by applying Eq 9, 10 to distributions Eq. 19: Thus, such a polarization-interference-based cultivation Eqs 13-15, 18, 19 of elements of the first-order differential matrix Eqs 1-10 provides layer-by-layer maps of linear and circular birefringence and dichroism of the myocardial fibrillary networks. This approach extends significantly the functionality of the 3D Mueller-matrix imaging technique for depolarization mapping of diffuse biological layers [48,49], that presently utilized only diagonal elements of the resulting Mueller matrix Eq. 19. Figure 1 shows the optical arrangement of the 3D Mueller matrixbased polarimetry system developed in-house. The parallel beam (/ 2 × 10 3 μm) of the He-Ne laser (λ 0.6328 μm), formed by the collimator, and the beam splitter, is divided into two equally intense beams. These are denoted as the irradiating and reference beams, respectively. The irradiating beam is guided through the polarisation filters 5-7 to the sample 8. A polarisation-inhomogeneous image of the object 8 is projected into the plane of the digital camera 14 (The Imaging Source DMK 41AU02.AS, monochrome 1/2 " CCD, Sony ICX205AL (progressive scan); resolution-1280 × 960; sensor area -7600 × 6200 μm; sensitivity -0.05 lx; dynamic range -8 bit; signal-to-noise ratio -9 bit) by the lens 9 (Nikon CFI Achromat P, focal distance -30mm, numerical aperture -0.1, magnification -×4). The reference beam is guided, by the mirror 4, through the polarisation filters 10-12 into the plane of the polarisation-inhomogeneous image of the object 9. As a result, an interference pattern is formed which is recorded by the digital camera 14. The formation of the required polarisation states of the irradiating and reference beams is carried out using polarisation filters 5-7 and 10-12, each of which contains two linear polarisers (B + W Kaesemann XS-Pro Polariser MRC Nano) and a quarter-wave plate (Achromatic True Zero-Order Waveplate).

Assessment of Optical Anisotropy
The layer-by-layer assessment of optical anisotropy W of myocardium histological sections is utilized by using first (Z 1 ) second (Z 2 ), third (Z 3 ) and fourth (Z 4 ) order statistical moments [14,15]: where C A × B -number of pixels of the photosensitive area of the CCD camera.

Biological Samples
Three groups of myocardial histological sections were utilized in the study. Control group-Group 1 with the myocardial tissues that have no relation to myocardial diseases, whereas Group 2 and Group 3 with the structural malformations in myocardial tissues caused by ACI and CHD, respectively. All the groups consisted of equal number of histological sections:k 36 .The histological sections of the myocardium were prepared according to the standard technique on a microtome with quick freezing. Figure 2 shows the original images of the histological sections of the myocardium from all three groups. The coordinate intensity distribution I(a, b) was normalized according to the maximum value I max in the image plane.
The obtained images (see Figure 2) show for all groups a fibrillar morphology structure consisting of a network of protein fibers, formed by optically active molecules of myosin and by optically isotropic molecules of actin [1].
Qualitative (visual) and quantitative (statistical) analysis of given microscopic images did not reveal significant differences between the three groups of myocardium histological samples. Table 1 presents the optical, geometrical and statistical parameters of the myocardium histological samples from each of the groups.
The attenuation (extinction) coefficient (τ, cm −1 ) of the myocardium histological samples was assessed by standard photometry spectral approach [50], utilizing the integrating sphere [51]. The degree of depolarization (DEP, %) was measured with standard Mueller-matrix polarimeter [14,15]. In terms of statistical significance [52] the standard deviation (ϑ 2 ≤ 0.025) corresponds to a confidence interval p30.05, which demonstrates the statistical reliability of the 3D Mueller-matrix mapping method. The quantitative differentiation of samples of histological sections of the myocardium of different groups by the method of statistical analysis of coordinate distributions I(a,b) Imax turned out to be statistically unreliable p_0.05.
The studies conducted in accordance with the principles of the Declaration of Helsinki, and in compliance with the International Conference on Harmonization-Good Clinical Practice and local regulatory requirements. Ethical approval was obtained from the Ethics Committee of the Bureau of Forensic Medicine of the Chernivtsi National University and the Bukovinian State Medical University (Chernivtsi, Ukraine), and written informed consent was obtained from all subjects prior to study initiation. The spatial order of fibrillary network that formed polycrystalline structure of myocardium, is optically discernible as linear birefringences and dichroism, whereas the optical activity of molecular domains of myosin forms a circular birefringence and dichroism [1-6, 17, 25]. The most optically expressed such These phenomena are seeing more clearly at the low level of depolarization background that is formed due to multiple scattering of light within the sample of histological sections of myocardium.
Therefore, polarization-based layer-by-layer detection of optical anisotropy and its variations provides a new quantitative approach of the evaluation of myocardial samples morphological structure and its pathological changes.

RESULTS AND DISCUSSION
A key feature of the morphological structure of myocardial tissue is the presence of a spatially ordered network of protein fibrils. These fibrils are formed by optically active molecules of myosin and isotropic actin proteins [1]. In terms of optics such a network has two types of anisotropy, namely the structural anisotropy leading to linear birefringence (Φ L ) and dichroism (Δ L ), and correspondingly to the circular birefringence (Φ C ) and dichroism (Δ C ) [25].
3D Mueller matrix-based polarimetry approach (see Figure 1) is utilized for preconstruction of the phase and amplitude anisotropy parameters to the functional diagnostic imaging of myocardial tissues in vitro. In particular functional capabilities of 3D Mueller-matrix layer-by-layer reconstruction of optical anisotropy parameters were investigated utilizing sample of myocardium histological section from Group 1. Figures 3, 4 show structural optical anisotropy obtained, respectively, for linear and circular birefringence and linear and circular dichroism for the different physical depths Eq. 16, and, correspondingly, Figures 5, 6 demonstrate structural optical anisotropy for the effective depths Eq. 17.
The next stage is devoted to demarcation of the samples of myocardium histological sections from Group 2 and Group 3 according the most effective phase plan Eq. 15 and the corresponding depth Eq 16, 17, guided with the following parameters: Δn ≈ 1.5 × 10 −3 ,h 50μm, λ 0.63μm.
A single pass of laser light through the histological section of the myocardium corresponds to the value θ(K 1) ≈ 0.75 rad 5 z(K 1) ≈ 50 µm double θ(K 2) ≈ 1.5 rad 5 z(K 2) ≈ 50 µm etc. In other words, the phase shifts correspond to the predominantly single scattering or low order of multiplicity of scattering. Whereas for θ ≥ 1.5rad and over the multiple scattering prevails. Therefor the complex amplitudes phase scanning are considered for both single scattering: θ 1 0.4rad5z 1 ≈ 25μm; θ 2 0.75rad5z 2 ≈ 50μm(see
The analysis of the embossed topographic maps of myocardial histological section with distinctive single scattering shows a peculiar spatial distribution of linear and circular birefringence (see Figure 3) and dichroism (see Figure 4) at the variety of physical depths (Z i ), defined by Eq. 16. The large-scale domains, the structure of which correlates with the size and directions of folding of protein fibrillar fibers, are seeing clearly both at the embossed topographic maps of circularly birefringence Φ L in Figures  3A,B and dichroism Δ L in Figures 4A,B. At the same time the small-scale domains that correspond to the coordinate positions of optically active molecular domains are observed at the embossed topographic maps of circularly birefringence Φ C (see Figures 3C,D) and dichroism Δ C (see Figures 4C,D).
With the increase of physical depth (Z i ) the magnitude and spread of random values of the parameters of linear and circular birefringence are also growing progressively.
The analysis of the embossed topographic maps of myocardial histological section at multiple scattering shows the significantly lower alterations in the structure of the spatial distributions of linear and circular birefringence (see Figure 5) and dichroism (see Figure 6) at the variety of physical depths (Z i *), defined by Eq. 17. The scale of domain structure of the embossed topographic maps of linear birefringence Φ L (see Figures 5A,B) and dichroism Δ L (see Figures 6A,B) notably decreases and no correlation with the size and directions of laying of protein fibrillar fibres is observed. The scale of the domain structure of the embossed topographic maps of circular birefringence Φ C (see Figures 5C,D) and dichroism Δ C (see Figures 6C,D) is reduced as well. An increase of effective depth (Z i *) From a physical point of view, this can be related to the peculiarities of the morphological structure of the myocardial tissue. Each partial fibril has a long-range geometric order that determines the structural anisotropy Φ L and Δ L along the direction of optical axis and the phase shift between the linearly and orthogonally polarised components of the laser wave [17,28,29]. For small values of the phase plane (θ 1 0.4rad5z 1 ≈ 25μm; θ 2 0.75rad5z 2 ≈ 50μm) within the volume of the histological section of the myocardium, single scattering events are dominant. Therefore, within the corresponding ϕ L (θ i , x, y) and Δ L (θ i , x, y) distribution's cross-sections, there is a direct relationship between the morphological structure of the myocardial fibrillar network and the parameters of linear anisotropy. An increase in the phase plane level θ 3 1.5rad5z p 3 ≈ 100μm; θ 4 2.3rad5z p 4 ≈ 150μm is accompanied by an increase in the average multiplicity of light scattering. As a result, the parameters of linear birefringence and dichroism are averaged. Quantitatively, this process manifests itself in an increase in the magnitude of the mean ϕ L (θ i , x, y) and the dispersion Δ L (θ i , x, y) of the distributions. The distributions of circular birefringence ϕ C (θ i , x, y) and dichroism Δ C (θ i , x, y) of protein complexes of the histological section of the myocardium have a different structure. The layered maps of these parameters present small-scale island structures with weak coordinate fluctuations. Physically, this can be attributed to the fact that the level of anisotropy of this type is determined by the concentration of optically active molecules of myosin, which are equally distributed in different phase sections of the myocardial tissue. Figure 7 shows the series of dependences of the statistical moments orders (Z 1;2;3;4 ), characterising the layer-by-layer distributions of the parameters of the anisotropy phase and amplitude of the partially depolarising layer of the myocardium tissue. To display the complete dynamics of changes in the value of statistical moments of the 3rd and 4th orders, which characterize the distributions W(θ i ), the corresponding individual range of values was selected along the ordinate for each type of optical anisotropy. The values of the statistical moments of the 1st and 2nd orders are presented with a coefficient ×1000.
Analysis of the obtained data revealed the following trends, which characterises the changes in the distributions of the parameters of the optical anisotropy of the myocardium layer: The observed regularities can be related to the fact that as the multiplicity of light scattering in the volume of the Frontiers in Physics | www.frontiersin.org October 2021 | Volume 9 | Article 737866 myocardium increases, in the limit of small θ i the asymmetric (Z 3;4 ≥ Z 1;2 ) distributions ϕ L , ϕ C , Δ L , Δ C (x, y), tend to normal (Z 1;2 ↑; Z 3;4 → 0) in accordance with the central limit theorem [50,51]. The statistical moments of the third and fourth orders are the most sensitive to changes in the polarisation manifestations of phase and amplitude anisotropy. The range of changes for these two cases spans an order of magnitude.

DIFFERENTIATION OF ACI AND CHD MYOCARDIUM HISTOLOGICAL SAMPLES
Finally, for the definitive diagnosis of ACI and CHD by the layer-by-layer Mueller-matrix approach described above the following protocol of reconstruction of myocardium optical anisotropy parameters has been developed. 1) As soon as the position of the phase plane most sensitive to pathological changes in the parameters of the optical anisotropy structure is determined θ p with the selected step of discrete phase "macro" scanning Δθ max The sensitivity (Se q q+g 100%), specificity (Sp j j+s 100%) and balanced accuracy (Ac 0.5(Se + Sp)) of the approach is estimated [53]. Here, q and g are the number of correct and incorrect diagnoses within Group 3, and j and s are the same within Group 2.
The experimental results of 3D layer-by-layer Mueller-matrix reconstruction are presented as embossed topographic maps [48] of linear and circular birefringence ( Figure 8) and dichroism (Figure 9) of ACI and CHD myocardial histological tissue samples.
As one can see the average magnitude of linear birefringence and dichroism (see Figures 8A,B, 9A,B) is prevail up to 3 times  compare to the magnitude of circular (see Figures 8C,D, 9C,D) birefringence and dichroism for both ACI and CHD myocardial histological samples. A decrease (↓) in the level of structural anisotropy (ϕ L , Δ L ) of the CHD myocardium is associated with degenerative-dystrophic changes in the spatial-angular ordering of myosin fibrillary networks. While changes in concentrations of optically active molecules (ϕ L , Δ L ) of myosin seems to be insignificant. Thus, the values of intergroup difference can be assessed quantitatively by handling the embossed maps of optical anisotropy myocardial histological tissue samples from Group 2 and Group 3 utilizing higher order statistical moments Z k (26). These values presented in Table 2. 4th order statistical moment (ΔZ 4 ), defining kurtosis of distributions of anisotropy phase and amplitude Φ L (θ p , a, b); Δ L (θ p , a, b) for myocardial histological tissue samples, demonstrates the highest sensitivity (see Table 2). The sensitivity Se, specificity Sp and balanced accuracy Ac of practical use of 4th order statistical moment for ACI and CHD diagnosis are presented in Table 3.
In comparison Table 4 presents the operational characteristics of diagnostic performance of the 3D Mueller-matrix reconstruction method of optical anisotropy of myocardial histological sections in other phase planes θ 0.9rad and θ 1.2rad, as well as the 2D Mueller-matrix reconstruction approach.

SUMMARY AND CONCLUSIONS
A method for 3D Mueller matrix reconstruction of layer-bylayer distributions of the linear and circular birefringence and dichroism has been, for the first time to the best of our knowledge, introduced theoretically in line with the protocol of experimental measurements. The proposed approach has been tested for the case of partially depolarising layers of myocardial tissues. The dynamics of the change in the magnitude of the statistical moments of the first to fourth orders that characterise the distribution of the parameters of the optical anisotropy of the polycrystalline structure of the partially depolarising layer (τ 1.02; Λ 57%) of the myocardial tissue sample in different phase sections of its volume was successfully investigated and analysed. The third and fourth order statistical moments are the most sensitive to changes in the polarisation manifestations of the phase and amplitude anisotropy. The range of changes for these two cases spans an order of magnitude. This study paves the way for the wider use of 3D Mueller matrix approach to the analysis and morphological imaging of optically anisotropic structures. The 3rd and 4th order statistical moments (Z 3 and Z 4 ) are found to be the most sensitive to the changes in phase anisotropy of myocardium histological sections (in a range θ i 0.01-3 rad) both for linear Φ L (Z 3 ≡ 1.27÷0.31; Z 4 ≡ 2.26÷0.39) and circular Φ C (Z 3 ≡ 4.13÷1.12; Z 4 ≡ 7.86÷1.36) birefringence; linear Δ L (Z 3 ≡ 1.23÷0.28; Z 4 ≡ 2.41÷0.37) and circular Δ C (Z 3 ≡ 1.49÷0.43; Z 4 ≡ 2.99÷0.78) dichroism. The optimal phase plane (θ p 0.7rad) has been found, for which an excellent ACI-CHD differentiation accuracy is achieved Ac(ΔZ 4 (θ p , Φ L , Δ L )) 93.05%÷95.8%.
Current study demonstrates that a 3D Mueller matrix method can be used to effectively analyse the optical anisotropy parameters of myocardial tissues with a strong potential for definitive diagnosis in forensic medicine.