Comparative Study of Modiﬁed Mueller Matrix Transformation and Polar Decomposition Parameters for Transmission and Backscattering Tissue Polarimetries

: Mueller matrix polarimetry is widely used in biomedical studies and applications, for it can provide abundant microstructural information about tissues. Recently, several methods have been proposed to decompose the Mueller matrix into groups of parameters related to speciﬁc optical properties which can be used to reveal the microstructural information of tissue samples more clearly and quantitatively. In this study, we thoroughly compare the differences among the parameters derived from the Mueller matrix polar decomposition (MMPD) and Mueller matrix transformation (MMT), which are two popular methods in tissue polarimetry studies and applications, while applying them on different tissue samples for both backscattering and transmission imaging. Based on the Mueller matrix data obtained using the setups, we carry out a comparative analysis of the parameters derived from both methods representing the same polarization properties, namely depolarization, linear retardance, fast axis orientation and diattenuation. IN particular, we propose several modiﬁed MMT parameters, whose abilities are also analyzed for revealing the information about the speciﬁc type of tissue samples. The results presented in this study evaluate the applicability of the original and modiﬁed MMT parameters, then give the suggestions for appropriate parameter selection in tissue polarimetry, which can be helpful for future biomedical and clinical applications.


Introduction
As a non-invasive, non-contact and label-free tool, Mueller matrix polarimetry is widely used in biomedical studies and clinical applications, for it can provide abundant microstructural information about tissues and cells [1][2][3][4]. Specifically, the Mueller matrix measurement system can be used for abnormal tissue detection for both backward imaging of bulk tissue samples and forward imaging of thin tissue slices [5][6][7]. However, although a Mueller matrix contains a comprehensive description of the polarization-related structural and optical properties of the medium, it is often difficult to obtain intrinsic polarimetric characteristics such as diattenuation, retardance and depolarization encoded in the Mueller matrix elements directly. To solve this problem, in the past decades several methods have been proposed to decompose the Mueller matrix into groups of parameters related to specific optical properties which can be used to reveal the microstructural information of tissue samples more clearly and quantitatively [8][9][10][11][12][13][14].
Until now, the Mueller matrix polar decomposition (MMPD) and Mueller matrix transformation (MMT) are two popular methods in tissue polarimetry studies and applications. Nevertheless, unlike the MMPD method, which decomposes the Mueller matrix into three sub-matrices with clear physical meanings of depolarization, retardance, and diattenuation [12,13], the MMT method directly combines groups of Mueller matrix elements to derive parameters related to certain characteristics of tissue-like media with an advantage of fast calculation speed [14]. Both the MMPD and MMT methods have demonstrated the diagnostic potential in various cancerous tissue detection, including breast cancer, liver cancer, inflammatory bowel diseases and cervical cancer [15][16][17][18][19][20]. However, it is worth noting that the parameters derived from the MMPD and MMT methods are based on specific deduction models and calculation assumptions. Therefore, for accurate diagnostic purposes, their applicability and interpretability regarding different types of samples requires further discussion and analysis. Recently, some efforts have been made in this area, and some useful conclusions have been drawn by Ahmad and us [21][22][23][24].
In this study, to compare the difference among the parameters derived from the MMPD and MMT methods more comprehensively while applying them on different types of tissue samples measured by both backscattering and transmission imaging setups, we prepare seven kinds of thin tissue samples for forward measurement and five kinds of bulk tissue samples for backward measurement. Based on the Mueller matrix data obtained using the setups, we carry out thoroughly comparative analysis of the parameters derived from both methods representing the same polarization properties, namely depolarization, retardance, fast axis orientation and diattenuation. Specifically, according to tissue structures, we further proposed several modified MMT parameters, whose abilities were also analyzed for revealing information with regard to specific types of tissue samples. We discuss the applicability of the original and modified MMT parameters based on the comparisons, and then give suggestions for appropriate parameter selection in tissue polarimetry.

Mueller Matrix Polar Decomposition (MMPD)
The Mueller matrix polar decomposition method proposed by Lu and Chipman describes the process of interaction between light and medium by decomposing a Mueller matrix into three sub-matrices with clear physical meanings: diattenuation (M D ), retardance (M R ) and depolarization (M ∆ ), as shown as Equation (1) [12].
Based on Equation (1), a group of MMPD parameters D, R and ∆ can be calculated by Equations (2)-(4), representing diattenuation, retardance and depolarization properties of the samples, respectively. It should be noted that the lower case m ij (i,j = 1, 2, 3, 4) in Equation (2) represents the corresponding Mueller matrix elements before performing any decomposition, while the lower case m ∆ in Equation (4)  (2) The parameters δ and θ derived from the retardance matrix M R are usually used to describe the value and fast axis orientation of linear birefringence induced by fibrous Appl. Sci. 2021, 11, 10416 3 of 20 tissue structures [13]. Since the birefringent fibrous structures including collagen, elastics, and muscle fibers that are prevalent in tissues, these two MMPD parameters can provide important information of the distribution behaviors of fibrous structures during many pathological processes [15][16][17][18][19][20][25][26][27][28]. It should be noted that compared with the original orientation parameters θ, a modified MMPD parameter θ e proposed in our recent study has been demonstrated as a more suitable indicator for describing the orientation of complicated tissue samples containing layered fibrous structures [26]. The parameters δ and θ e can be calculated according to Equations (5)-(8) [13,26]. The M R (i,j) (i,j = 2, 3) in Equation (5) represents the corresponding elements of the 4 × 4 retardance matrix M R , while R is the total retardance, ε ijk is the Levi-Cività permutation symbol, δ ij is the Kronecker delta, and m R is the 3 × 3 retardance sub-matrix striking out the first row and the first column of M R .

Mueller Matrix Transformation (MMT) Parameters
The Mueller matrix transformation method was proposed to derive groups of quantitative parameters for describing the microstructural features of biomedical tissue samples by fitting the Mueller matrix elements to the trigonometric functions [14]. Several MMT parameters such as t 2 , t 3 , x 3 and 1 − b were introduced in our previous studies [14,29], where t 2 (also denoted as t 1213 shown in Equation (9)) reveals the linear diattenuation property, t 3 (also denoted as t 2434 and t 4243 shown in Equations (10) and (11)) is related to the linear retardance of the tissue samples, while x 3 (also denoted as φ 2434 and φ 4243 shown in Equations (12) and (13)) can characterize the orientation of the birefringence structure and 1 − b (shown in Equation (14)) represents the linear depolarization of the tissue samples.
However, as these parameters are derived based on specific models [14], when we apply these parameters on complex tissue samples, their applicability needs further discussion. Our previous studies showed that there is a deviation between the values of t 2434 and t 4243 when multiple crossing linear birefringence effects coexist, as the information of linear retardance are distributed among the m 24 , m 34 , m 42 and m 43 elements [14,26,30]. Such deviation also occurs in the birefringence orientation parameters φ 2434 and φ 4243 .
Here, according to tissue structures, we propose two modified MMT parameters t qr and φ qr (shown in Equations (15) and (16)), which take m 24 , m 34 , m 42 and m 43 elements into account, to calculate the linear retardance and birefringence orientation of the layered tissue sample more accurately. Besides, we also propose another two MMT parameters t 121314 and 1 − b m (Equations (17) and (18)) to describe the overall diattenuation and depolarization properties of the tissue samples by taking the elements m 14 and m 44 into account, which contain the information of circular diattenuation and circular depolarization.
Among the available Mueller matrix analyzing methods, the calculation speed of MMT is relatively fast, which is advantageous for clinical applications [5,24]. However, the principle of MMT calculations is based on the fitting parameters obtained by experiments and simulations, and its accuracy needs to be further studied when applied to quantitative characterization of specific tissues and cells. On the other hand, the MMPD parameters are widely used in the field of biomedical and pre-clinical detection. The applicability of the MMPD method on tissue samples has been proven in various ex vivo and in vivo measurements [5,6].

Materials and Measurement Setup
In previous studies, we have developed the transmission Mueller matrix microscope for measuring thin tissue slices and backscattering Mueller matrix imaging set-ups for measuring bulk tissue samples [19,[31][32][33][34][35]. The measurement schemes are based on the dualrotating retarder method. For this method, the PSG (polarization state generator) and PSA (polarization state analyzer) consist of a pair of polarizers fixed in the horizontal direction, which is also parallel to the source-sample-camera triangle plane for the backscattering imaging, and two retarders rotating with a fixed ratio of angles. It should be noted that the source, sample and camera are colinear in transmission imaging geometry. Then the 16 Mueller matrix elements can be calculated by using the Fourier coefficients α n and β n shown as Equation (19) [36].
where I is the intensity, θ 1 is the rotation angle of the retarder 1 for the PSG. The two retarders (retarder 1 for PSG, retarder 2 for PSA) rotate with a fixed rate θ 1 = 5θ 2 . The Mueller matrix measurement setups used in this study were calibrated to ensure that the maximum errors of the measured Mueller matrix elements are less than 1%. More details on this Mueller matrix imaging method and calibration process can be found in [32,37,38].
Here, we prepared two types of tissue samples (n = 149) for the comparative analysis: thin pathological tissue slices for transmission measurement and bulk tissue samples for backscattering measurement. The unstained thin tissue samples include human colorectal, intestinal tuberculosis (ITB), liver, Crohn's disease (CD), bladder, breast, adenoma tissue slices whose thickness are 12 µm, while the bulk tissue samples include human breast, rat-skin, porcine intestine, porcine stomach, and porcine liver tissues. The H-E stained 4-µm-thick tissue slices of the samples were also prepared for pathological observations. More details of the samples used in this study are shown in Table 1. It should be noted that the signal-to-noise ratio for the camera is related to the energy obtained from the source, which depends on both the power of LED and the exposure of the camera. As the power of LED for our Mueller matrix imaging is commonly within 1-3 W [24,25,34], we have selected an appropriate exposure time for the camera to achieve a good signal-to-noise ratio.
Our previous studies showed that the LED with the power of 1 W did not generate a lot of heat to the thick sample when an appropriate exposure time was selected [24,25,34]. In this work, the human tissue samples are provided by the Department of Gastroenterology, Zhujiang Hospital, Southern Medical University, and Shenzhen Sixth People's (Nanshan) Hospital. This work was approved by the Ethics Committee of the Shenzhen International Graduate School, Tsinghua University.

Quantitative Indicators for MMT and MMPD Parameters Comparison
Here, we evaluate the difference between MMT and MMPD parameters in two aspects: how their linear relationship is and whether there exits statistical significance. For the linear relationship comparison, we take the square of the Pearson correlation coefficient (R 2 or r 2 ) as a measure between two sets of parameters, as it can provide a quantitative value for their linear relationship evaluation (range 0-1). For a correlation between variables x and y, the formula for calculating the sample Pearson's correlation coefficient (R of r) is given by Equation (20) [39,40].
where x i and y i are the values of x and y for the ith individual. It is also worth mentioning that the larger the value is, the more similar the changing trends of MMT and MMPD parameters are. Based on the correlation analysis, we conduct the statistical analysis between MMT and MMPD with p-value, which can assess whether both parameters are interchangeable without linear fitting according to their statistical differences [41].

Results and Discussion
To compare the MMT and MMPD parameters more comprehensively, in the following sections we analyze the three main polarization properties (diattenuation, retardance, and depolarization) of both the transmission imaging results of thin tissue samples, and backscattering imaging results of bulk tissue samples.

Quantitative Comparison of Diattenuation Related MMPD and MMT Parameters for Transmission Imaging of Thin Tissue Samples
We measured all seven different kinds of thin pathological tissue slices (86 in total) using the forward transmission Mueller matrix microscope. Then, for each kind of tissue sample, we calculated the diattenuation related MMT and MMPD parameters according to Equations (2), (9) and (17). Specifically, to find out whether the modified MMT parameters can be used to accurately obtain the structural properties of tissue samples, the quantitative analyses were conducted to show the correlations and differences between the MMT and the corresponding MMPD parameters.
The correlation coefficient R 2 and statistical analysis p values of each kind of thin tissue slice samples are presented in Table 2. The comparison of the MMPD and MMT diattenuation parameters was performed in two groups: (D, t 1213 ) and (D, t 121314 ). As we can see from Table 2, for the correlation analysis, compared with the group of (D, t 1213 ), the group (D, t 121314 ), for each kind of sample has a higher correlation coefficient value of 1. More thorough and intuitive comparisons are illustrated in Figure 1a, c by plotting out all the sample data, which confirm that the modified MMT diattenuation parameter t 121314 has a good accordance with the MMPD parameter D in the transmission measurement. On the other hand, as can be noticed from Table 2 and Figure 1b, the ITB tissue samples show statistically significant differences, with p < 0.01 between the parameters D and t 1213 . While for the modified parameter t 121314 , the p-values of all kinds of tissue samples are above 0.05, meaning that the modified MMT parameter t 121314 can be used as a precise indicator for the diattenuation property of thin tissue samples.

Quantitative Comparison of Diattenuation Related MMPD and MMT Parameters for Backscattering Imaging of Bulk Tissue Samples
To study the relationship between the diattenuation related MMPD and MMT parameters for backscattering polarimetric measurement, we prepared five different kinds of bulk tissue samples (63 in total) to attain their backscattering Mueller matrix data. Then, for each kind of tissue sample, we also calculated the diattenuation-related parameters in Section 3.1 for the correlation and statistical analysis whose results are shown in Table 3 and Figure 2.

Quantitative Comparison of Diattenuation Related MMPD and MMT Parameters for Backscattering Imaging of Bulk Tissue Samples
To study the relationship between the diattenuation related MMPD and MMT parameters for backscattering polarimetric measurement, we prepared five different kinds of bulk tissue samples (63 in total) to attain their backscattering Mueller matrix data. Then, for each kind of tissue sample, we also calculated the diattenuation-related parameters in Section 3.1 for the correlation and statistical analysis whose results are shown in Table 3 and Figure 2.
We can observe from Table 3 that, for the backscattering measurement results, the correlation coefficient values of parameters (D, t121314) for each kind of bulk samples also have higher values of 1 compared with that of parameters (D, t1213), meaning that the modified MMT diattenuation parameter t121314 also has a good accordance with the MMPD parameter D in backscattering measurement. Meanwhile, for the statistical analysis of each kind of sample, it can be seen from Table 3 and Figure 2b that the rat-skin, porcine-intestine, porcine-stomach and liver tissue samples show statistically significant differences with p < 0.05 between the parameters D and t1213. However, the p-values of all kinds of tissue samples are larger than 0.05 for parameters (D, t121314) shown in Table 3 and Figure 2d, which also confirms that the modified MMT parameter t121314 can be used    We can observe from Table 3 that, for the backscattering measurement results, the correlation coefficient values of parameters (D, t 121314 ) for each kind of bulk samples also have higher values of 1 compared with that of parameters (D, t 1213 ), meaning that the modified MMT diattenuation parameter t 121314 also has a good accordance with the MMPD parameter D in backscattering measurement. Meanwhile, for the statistical analysis of each kind of sample, it can be seen from Table 3 and Figure 2b that the rat-skin, porcineintestine, porcine-stomach and liver tissue samples show statistically significant differences with p < 0.05 between the parameters D and t 1213 . However, the p-values of all kinds of tissue samples are larger than 0.05 for parameters (D, t 121314 ) shown in Table 3 and Figure 2d, which also confirms that the modified MMT parameter t 121314 can be used to accurately characterize the diattenuation property of the bulk tissue samples measured by backscattering polarimetry.
In summary, we can see from the comparative analyzing results presented in Sections 3.1 and 3.2 that the modified MMT parameter t 121314 shows a better accordance with the MMPD parameter D in both the transmission and backscattering measurements than the original parameter t 1213 , which only reflects the linear diattenuation of the samples. Though the circular diattenuation effects of the tissue samples are often very limited, taking the element m 14 (encoded the circular diattenuation property) [42] into account for tissue polarimetric calculations can produce more accurate results.

Quantitative Comparison of Linear Retardance Related MMPD and MMT Parameters for Transmission Imaging of Thin Tissue Samples
Linear retardance is viewed as one of the most important polarization properties of tissue samples in terms of revealing the location and density of the birefringent fibrous structures [15][16][17][18][19][20]. Our previous studies showed that the MMT parameters t 2434 and t 4243 can be used to reflect the linear retardance of the tissues [14,29]. However, both of the parameters may slightly deviate from each other when measuring complicated tissue samples containing layered anisotropic structures, reducing the accuracy of Mueller matrix polarimetry for quantitative tissue evaluation. To deal with this problem, we proposed a modified MMT parameter t qr based on the elements from both the fourth column and fourth row shown in Equation (15) to describe the linear retardance property. To verify its effectiveness, we conducted the comparisons among the linear retardance related MMPD and MMT parameters in three groups: (δ, t 2434 ), (δ, t 4243 ), and (δ, t qr ). The correlation and statistical analysis values of MMT and MMPD parameters groups for each kind of thin tissue samples are shown in Table 4 and Figure 3. For more intuitive comparisons, all the parameters are normalized to the [0, 1] value range. For the correlation analysis, Figure 3a,c,e show that, compared with the parameter groups (δ, t 2434 ) and (δ, t 4243 ), the group (δ, t qr ) for all the samples have higher correlation coefficient values. We can also see the differences more clearly in Table 4, where the correlation coefficient values of (δ, t 2434 ) and (δ, t 4243 ) deviate from each other ranging between 0.90 to 0.99, while the correlation coefficients of (δ, t qr ) are more stable with higher values for all the samples. Moreover, for the statistical analysis, all kinds of thin tissue samples show statistically significant differences with p < 0.01 among the parameter groups according to Table 4 and Figure 3b,d,f, meaning that the MMT linear retardance parameters cannot be substituted for MMPD parameter δ directly without linear fitting. However, the modified MMT linear retardance parameter t qr shows a good accordance with the MMPD parameter δ, indicating that it can be used to well present linear retardance distributions of tissues by linear fitting with MMPD parameter δ.

Quantitative Comparison of Linear Retardance Related MMPD and MMT Parameters for Backscattering Imaging of Bulk Tissue Samples
As shown in Table 5, for the backscattering measurement results of bulk tissue samples, the linear retardance related parameters derived from MMT and MMPD methods do not show a good correlation. The correlation coefficient values of each kind of bulk tissue samples for parameter groups (δ, t 2434 ), (δ, t 4243 ) and (δ, t qr ) are between 0.16 and 0.89, and we can also see the data points of the tissue samples distribute randomly in Figure 4a,c,e. Besides, Figure 4b,d,f show that the MMPD linear retardance parameters are larger than the MMT linear retardance parameters, and they all have statistically significant differences with the p-value smaller than 0.05. Therefore, it is worth mentioning that when we try to measure linear retardance property of bulk tissues sample using backscattering polarimetry, the MMPD parameter δ will provide the information more accurately, while the MMT linear retardance parameters can be suitable for quickly evaluating the existence and location of birefringent structures.

Quantitative Comparison of Linear Birefringence Orientation Related MMPD and MMT Parameters for Transmission Imaging of Thin Tissue Samples
Besides the linear retardance value, the fast axis orientation of linear birefringence is another important indicator to reveal the structural features of tissue, such as the distribution behaviors of collagen and elastic fibers in pathological tissue samples, which is helpful for clinical diagnosis. In our previous study, considering that the layered fibrous structures can result in deviations of MMPD linear birefringence orientation parameters, we proposed a modified MMPD parameter θ e to more accurately represent the fast axis orientation of the complicated linear birefringent structures in tissue samples [26]. Meanwhile, it was also found that while applying the MMT parameters φ 2434 , φ 4243 to the fibrous tissue samples, they exist with positive or negative deviations from the ideal values. Hence, according to the experimental results, we proposed a modified MMT parameter, φ qr shown in Equation (16), for layered linear birefringent structures. To verify its effectiveness, the correlation and statistical analysis are conducted with the variance of parameters groups (θ e , φ 2434 ), (θ e , φ 4243 ) and (θ e , φ qr ), whose results are shown in Table 6 and Figure 5.     For the correlation analysis, Figure 5a,c,e show that compared with the parameter groups (θ e , φ 2434 ), (θ e , φ 4243 ) and (θ e , φ qr ), the correlation coefficient values of (θ e , φ qr ) for all the samples are larger. We can also find the difference more clearly in Table 6, where the correlation coefficient values of (θ e , φ 2434 ) and (θ e , φ 4243 ) deviate with each other ranging between 0.44 to 0.98, while for (θ e , φ qr ) the values are more stable and closer to 1 for all kinds of thin tissue samples. For the statistical analysis, the group (θ e , φ qr ) shows no statistically significant difference with p > 0.05 according to Table 6 and Figure 5b,d,f, which means that the modified MMT parameters φ qr has the similar ability as the MMPD parameter θ e for revealing the linear birefringence orientation information of complicated thin tissue samples.

Quantitative Comparison of Linear Birefringence Orientation Related MMPD and MMT Parameters for Backscattering Imaging of Bulk Tissue Samples
For the backscattering measurement results of bulk tissue samples, the modified linear birefringence orientation MMT parameter φ qr also shows a good fit with the MMPD parameter θ e . As shown in Table 7, the correlation coefficient values of each kind of bulk tissue samples for the parameters group (θ e , φ qr ) are all equal to 0.99, while those for (θ e , φ 2434 ), (θ e , φ 4243 ) range between 0.58 and 0.97. We can also see that the data points of all tissue samples are evenly distributed along a straight line in Figure 6e. Particularly, Table 7 and Figure 6f also confirm that the parameters (θ e , φ qr ) show no statistically significant differences with the p > 0.05. Therefore, when measuring linear birefringence orientation of bulk tissue samples using backscattering polarimetry, both the parameters φ qr and θ e can provide quantitatively similar information, and φ qr can be calculated faster. In summary, the results shown in Sections 3.5 and 3.6 demonstrate that among the available MMT linear birefringence orientation parameters, φ qr shows the best accordance with the MMPD parameter θ e . The parameters φ 2434 and φ 4243 only contain partial linear birefringence orientation information of tissue samples with layered fibrous structures. Thus, when taking the elements m 24 , m 34 , m 42 , and m 43 from both the fourth column and fourth row of the Mueller matrix into account, the modified parameter φ qr can well present the linear birefringence orientation information more thoroughly.

Quantitative Comparison of Depolarization MMPD and MMT Parameters for Transmission Imaging of Thin Tissue Samples
In previous studies, both the MMT parameter 1 − b shown in Equation (14) and MMPD parameter ∆ shown in Equation (4) have been used to represent the depolarization property of tissue samples [5]. However, there would be a deviation between these two parameters when the medium contains a circular depolarization property encoded in the element m 44 , which the parameter 1-b does not take into account. To more accurately represent the overall depolarization of tissues, we proposed a modified parameter 1 − b m as Equation (18) shows. The comparison results of the MMT and MMPD depolarization parameters for 7 kinds of thin tissue samples are shown in Table 8 and Figure 7.

Type Tissues
Depolarization

Quantitative Comparison of Linear Birefringence Orientation Related MMPD and MMT Parameters for Backscattering Imaging of Bulk Tissue Samples
For the bulk tissue samples, both the MMT parameters 1 − b and 1 − bm show a good accordance with the MMPD parameter Δ. As shown in Table 9, the correlation coefficient values of all kinds of bulk tissue samples for parameters group (Δ, 1 − bm) are equal to 0.99, which differ from that for the group (Δ, 1 − b) ranging between 0.96 and 0.98. For the correlation analysis, the values of parameters groups of (∆, 1 − b) and (∆, 1 − b m ) are close and much smaller than 1 for most of the tissue samples indicated in Table 8, which means both groups of parameters hardly exist in a linear relationship. More intuitive comparisons illustrated in Figure 7a,c confirm that all the sample data points distribute randomly. For the statistical analysis, most of the samples show statistically significant differences with p < 0.01 for the parameters groups (∆, 1 − b) and (∆, 1 − b m ) shown in Table 8. From Figure 7b,d, we can find that the MMPD and MMT depolarization parameters values of thin tissue samples are smaller than 0.02. This is the reason why the parameters (∆, 1 − b) and (∆, 1 − b m ) show no significant difference in correlation analysis for the thin tissue samples with a limited thickness.

Quantitative Comparison of Linear Birefringence Orientation Related MMPD and MMT Parameters for Backscattering Imaging of Bulk Tissue Samples
For the bulk tissue samples, both the MMT parameters 1 − b and 1 − b m show a good accordance with the MMPD parameter ∆. As shown in Table 9, the correlation coefficient values of all kinds of bulk tissue samples for parameters group (∆, 1 − b m ) are equal to 0.99, which differ from that for the group (∆, 1 − b) ranging between 0.96 and 0.98. Moreover, we can also see that the data points of the tissue samples evenly distributed along a straight line in Figure 8a,c. Besides, Figure 8b,d indicate that the MMPD depolarization parameter ∆ is approximately equal to the MMT parameters. Specifically, the parameters group (∆, 1 − b m ) shows no statistically significant difference with p > 0.05. Therefore, it is worth noting that when bulk tissue samples are measured using backscattering equipment, the parameter 1 − b m can be used to extract the depolarization information quickly and accurately.
Appl. Sci. 2021, 11, x FOR PEER REVIEW 19 of 22 Table 9. Correlation and statistical analysis comparison results between backward measured MMT and MMPD depolarization parameters for each kind of bulk tissue samples.

Type Tissues
Depolarization

Conclusions
In this study we compared the parameters derived from the MMT and MMPD methods thoroughly by measuring two types of tissue samples: seven kinds of thin tissue slices and five kinds of bulk tissue samples, with a transmission Mueller matrix microscope and backscattering Mueller matrix measurement setup. After grouping the pa-  Table 9. Correlation and statistical analysis comparison results between backward measured MMT and MMPD depolarization parameters for each kind of bulk tissue samples.

Type Tissues
Depolarization

Conclusions
In this study we compared the parameters derived from the MMT and MMPD methods thoroughly by measuring two types of tissue samples: seven kinds of thin tissue slices and five kinds of bulk tissue samples, with a transmission Mueller matrix microscope and backscattering Mueller matrix measurement setup. After grouping the parameters with the same polarization properties, namely diattenuation, linear retardance, linear birefringence fast axis orientation and depolarization, we performed both correlation and statistical analysis. Our preliminary experimental results showed that: (1) for obtaining the diattenuation property of complex tissues, we can calculate it through the modified MMT parameter t 121314 , with a fast speed and a high correlation with the MMPD parameter D; (2) for revealing the linear retardance value related to the density of birefringent structures like layered fibers of thin tissues, the modified MMT parameter t qr can reach a high linear correlation with the MMPD parameter δ. While for the bulk tissue samples, it is more reasonable to use the MMPD parameter δ to extract the linear retardance value accurately; (3) for characterizing the orientation distribution of the birefringent structures, the modified MMT parameter φ qr and MMPD parameter θ e can both reach similar accurate results for thin and thick tissues; (4) the modified MMT parameter 1 − b m shows a high degree of consistency with the MMPD parameter ∆ for the bulk tissue samples with strong depolarization, so it can be used to quickly evaluate the tissue depolarization property. In summary, based on the analysis and discussion regarding the applicability of the MMPD and MMT parameters, this study gave suggestions for the appropriate selection of parameters in Mueller matrix imaging for different types of tissue samples, which can be useful for biomedical and clinical polarimetry.