Mueller Matrix-Based Approach for the Ex Vivo Detection of Riboﬂavin-Treated Transparent Biotissue

: Corneal collagen cross-linking is an established procedure for the treatment of certain eye diseases which is applied to enhance the mechanical stability of such biotissue without deteriorating its functionality. However, being transparent, the optical analysis of the outcome of such treatments is cumbersome and relies on relatively expensive experimental equipment. We aim to apply the Mueller matrix polarimetry for the detection of photo-induced collagen cross-linking in transparent biotissue after treatment with riboﬂavin and UV irradiation. A simple Mueller matrix polarimetry setup could provide a fast and non-invasive analysis of transparent media to sensitively detect small photo-induced cross-linking effects in biotissue. We demonstrated the current capabilities of the approach on non-planar porcine cornea samples ex vivo. We reported the distinction between untreated and riboﬂavin-treated samples. The differences observed were correlated with the variation of certain Mueller matrix elements and parameters derived from the decomposition. The measurement data show variation in the cross-linked and non-cross-linked samples, although the effect of the UV treatment on the riboﬂavin-treated samples was not at the same level of signiﬁcance yet and needs further investigation. The Mueller matrix measurement represents a promising approach for the detection of the effects of corneal collagen cross-linking. Further studies with a larger sample number are required to validate this approach. In the future, this could enable the reliable and non-invasive detection of photo-induced effects in biotissue and open the possibility for in vivo application, e.g., in eye disease treatment or the detection of scar collagen development.


Introduction
The fast and non-invasive detection of tissue properties can be relevant to investigations in the life sciences or medical diagnostics and could provide interesting alternatives to invasive approaches.Optical techniques have already proven their impressive capabilities for the detection of subtle effects in biotissues, in many cases also for in vivo measurements.However, techniques such as Brillouin microscopy [1], optical coherence tomography (OCT) [2], optical coherence elastography [3,4], laser scanning microscopy [5], hyperspectral imaging [6], multi-photon imaging [7] or diffuse reflectance methods [8,9] rely on relatively sophisticated equipment and might not be applicable to certain use cases.With regard to achievable contrast or sensitivity, all approaches have their individual strengths and weaknesses.In particular, the visualization of effects in transparent biological media with reasonable sensitivity and speed, ideally in vivo and under real conditions as Appl.Sci.2021, 11, 11515.https://doi.org/10.3390/app112311515https://www.mdpi.com/journal/applsciclinical settings represent an important challenge to be addressed.Additionally, as laser light is employed as the coherent radiation in most cases, safety aspects must be considered.
The highly ordered structure of some bio-tissues [10,11] provides measurable depolarization values in polarimetry setups [12].The human cornea is one such bio-tissue consisting of successive layers of collagen fibrils immersed in an amorphous ground substance containing water, proteins and salts [13].Such a formation is optically anisotropic and is therefore a suitable sample for polarimetry measurements [14].Corneal collagen cross-linking (CXL) is a promising treatment that increases the mechanical stiffness of the cornea.This is particularly useful for several diseases in ophthalmology such as keratoconus in which the cornea thins over time.The treatment combines the application of ultraviolet-A light and a chromophore-like riboflavin following the Dresden protocol to induce polymerization [15].The quantitative assessment of the treatment was performed via corneal stroma depth and biomechanical stiffness measurements.To date, optical methods such as OCT [16], confocal microscopy [17] and non-linear microscopy [18] were also applied for this purpose.Mueller matrix (MM) polarimetry is promising for the evaluation of CXL treatment because the degree of orientation in the samples can in principle be detected by polarization-sensitive measurements [19,20].Compared to the optical methods mentioned above, it provides a relatively simple and straightforward approach that delivers information about the polarization changing properties of a sample [9,[21][22][23].This relies on the measurement of the laser intensity reflected from or transmitted through a sample for different polarization states of the radiation before and after the interaction with the biotissue.Among many others, Mueller matrix polarimetry was utilized for the examination of skin tissue properties in dermatology [24], the assessment of the quality of tissue implants [20] or the identification of cancerous tissue sections from excised samples [25].In addition, the technique was applied for the evaluation of β-amyloid progression in Alzheimer's disease [21].Depending on the application, measurements in reflection or transmission mode are required and can be realized.As the Mueller matrix contains the full polarization information of a sample, it is advantageous over other biotissue imaging techniques such as non-contact dermoscopy which uses cross-polarization configurations for imaging [19,26,27].
In this work, we employed Mueller matrix polarimetry for the fast and non-invasive analysis of transparent biological media.As an example, we investigated porcine cornea samples ex vivo to demonstrate that small effects originating from photo-induced collagen cross-linking processes in the tissue can be detected.In the long term, this technique might represent an interesting alternative to the detection of cross-linking processes in human cornea in vivo and could lead to a relatively simple and cost-effective approach for the examination of therapy outcomes in ophthalmology.Moreover, it might be applicable for the quantification of collagen deposition in scars, which is relevant for wound healing, or pathological scar treatment and where the in vivo measurement of small changes of the tissues properties is desirable [28].In order to evaluate the capability of the technique, in a first attempt, we applied an established protocol to the porcine cornea samples consisting of the application of riboflavin followed by UV irradiation to induce the cross-linking of the collagen fibers contained therein [29].

Photo-Induced Corneal Collagen Cross-Linking
The ex vivo porcine eye cornea samples were obtained from a local slaughterhouse.The cornea was cut and stored in a MEM (minimal essential medium with Earle's salts) solution with 6% Dextran 500, 2% FBS (fetal bovine serum) and antibiotic-antimycotic solution overnight.On the measurement day, we placed the cornea in an artificial chamber where it was stored in a 15% Dextran 500 solution for at least 90 min to minimize the influence of dehydration.Stable hydration conditions were reached if the thickness of the sample differed by less than 3% in three consecutive measurements taken 5 min apart from each other.The pachymetry was controlled before every step of the measurement procedure.For the MM measurements, as mentioned above, we considered corneal collagen cross-linking (CXL) as a relevant use case for a potential future application [29].In this case, the aim was to change the rigidity of the cornea by cross-linking the collagen fiber contained therein through the generation of oxygen radicals.This method was used as a therapy for corneal degeneration [30].To do this, the cornea samples were exposed to 0.1% riboflavin (vitamin B2) in 15% Dextran 500 solution which diffuses into the tissue.The samples are then irradiated with UV irradiation according to the established Dresden protocol [29].Here, this procedure was applied to the ex vivo porcine cornea samples.For this purpose, the samples were first stored in a riboflavin solution (0.1%) for 30 min and then exposed to UV irradiation at 365 nm at a fluence of 3 mW/cm 2 for another 30 min.As in everyday clinical practice, the protocol was often adapted and the UV exposure time was reduced to 10 min when applying 9 mW/cm 2 [31].In this work, we used the modified protocol for preparation.The aim was to identify whether the cross-linking was induced inside the cornea which is not observable by eye but can be detected by the Mueller matrix measurements.

Mueller Matrix Measurement Setup
Mueller matrix (MM) polarimetry delivers information about the polarization changing properties of the samples under study.This information can be derived from the measurement of the intensity of light which is reflected or transmitted through a sample following illumination at different polarization states.The Mueller matrix relates to the Stokes formalism used to describe the polarization state of light before and after interaction with the sample [32][33][34][35].The Stokes vector S o of the light field after interaction with a sample can be calculated from the Stokes vector of the incoming light S i for known MM as follows: Here, the 4 × 4 MM contains all information about the polarization-changing properties of the samples.The Stokes vectors consist of defined intensity values of different polarization states: The indices for the intensity values I stand for the different polarization states of the light (H: horizontal polarization; V: vertical polarization; P: light polarized at an angle of 45 • ; M: light polarized at −45 • ; R: right circular polarization; L: left circular polarization).For the measurement of the MM, the change in the Stokes vector after the interaction of the incoming light field with the sample was recorded.In total, this typically requires either 16 or 36 different measurements depending on the polarization states which need to be generated, i.e., 4 states (H, V, P, R) or 6 states (H, V, P, M, R, L), respectively, as can be seen in Equation (2).Details on the formalism can be found in reference [36].In this work, we determined 6 polarization states and performed 36 different intensity measurements, in order to compensate for the calibration and measurement errors and increase the accuracy.This, however, increases the measurement times and needs to be individually considered with respect to the demands of each application, particularly when in vivo measurements are required.As the experimentally obtained values for the MM elements cannot be directly related to the physical properties of a sample, the decomposition of the MM into matrices of known contributions is usually performed.Here, we employed the polar decomposition by Lu and Chipman, which was widely used in the literature as [37]: M exp represents the experimentally measured MM, M ∆ a depolarizer, M R a retarder and M D a diattenuator.For a pure depolarizer, M ∆ can be written in the matrix form as Further to that, key figures are the parameter ∆, the depolarization power; R, the total retardance (combination of the effect of linear and circular birefringence), and P, the total polarizance.These parameters can be calculated to determine specific sample properties by Equations ( 5)- (7) [37] as follows: For the measurements reported in this work, we used our in-house polarimetry system which was described in detail in [20].A sketch of the system is shown in Figure 1.Here, we briefly summarize the main aspects of the system and methodology.The MM measurement system generates the polarization states of the incoming light using a polarization state generator (PSG) and records the outgoing states upon the interaction with the samples and passing through a polarization state analyzer (PSA).Using a CCD camera (BFS-U3-32S4M-C, FLIR Integrated Imaging Solutions Inc., Richmond, BC, Canada), 2D measurements of the MM of a sample were obtained in one shot [20].The PSG and PSA components were assembled from fixed polarizers (LPVISE100-A2", Thorlabs, Newton, NJ, USA) and liquid crystal retarders (LCRs) (LCC1223T, Thorlabs, Newton, NJ, USA).By changing the voltage applied to the LCR, all required polarization states can be generated in approximately 15 s in total.The system was designed to operate in reflection and transmission mode and uses laser light at a wavelength of 532 nm (CW532-04-1, Roithner Lasertechnik GmbH, Vienna, Austria) for the measurement.From the measured intensity values, the MM matrix elements were calculated for each pixel, as can be seen in Table 1.The setup was calibrated to ensure the optical system is adjusted and behaves similarly to an ideal system.This requires the LCRs to be oriented at the correct angle with respect to each another.For this, the retardation to be set at the 532 nm wavelength needs to be determined.The detailed calibration procedure and steps are described in refs.[38,39].Here, we summarized the main aspects in brief.The calibration relies on applying suitable values for the supply voltages of the LCRs for achieving λ/4 and λ/2 phase delays, respectively.All angles are aligned with respect to the first polarizer in the PSG which does not change its orientation during the whole process.The second polarizer is mounted in a motorized rotating holder and is then oriented perpendicular to the first one.Then, the first polarizer is manually set to 45 • and the following steps are implemented: (i) the fast axes of all LCRs are calibrated to 0 • ; (ii) a λ/2 retardation is calibrated for all LCRs by the analysis of the relation between the applied supply voltage and the intensity at the sensor; and (iii) a λ/4 retardation is calibrated accordingly for LCR2 and LCR3 by using a quarter-wave-plate.The procedure is based on intensity measurements with the imaging camera acting as the sensor.The displayed inhomogeneity of the Mueller matrix entries in Figure 2 can be partly due to changes in the required phase delay over the aperture of the LCRs [40].However, measurements of the homogeneous samples such as air or polarizers (not displayed in this work) show homogeneity over the illuminated field-of-view.However, the size of the effect is not homogeneously distributed over the surface of the cornea, which is evident from the individual entries and probably due to non-optimal sample fixation.

Polarimetry of Porcine Cornea in Transmission and Reflection Mode
In a first step, we recorded the MM of unmodified porcine cornea ex vivo by using our polarimetry system in transmission mode [20].For this purpose, the cornea was carefully placed in a 3D printed ring holder.The result of a MM measurement with illumination at a wavelength of 532 nm is shown in Figure 2.
A clear non-zero signal can be seen in some diagonal and off-diagonal entries in the MM, with M(2, 2), M(2, 3), M(3, 2), M(3, 3) and M(4, 4) having the largest mean values.These values indicate the polarization-changing properties of the sample.The signals are not homogeneously distributed over the whole surface area of the cornea.This may be because the cornea cannot be held with tension in the current holder.Therefore, the surface remains curved.The birefringent properties of the cornea are already known from the literature [41].Spatially resolved features and the anisotropy of the human cornea were shown using an imaging polarimeter, demonstrating negative biaxial material features [42,43].It was clearly observed here that our polarimetry setup was capable of resolving spatial features on such samples.In order to improve similar measurements in the future, the special cornea sample holders consisting of spherically shaped quartz glass elements could be used.
The measurements, discussed in the remainder of this work, were carried out in reflection mode as shown in Figure 3 [38].The effect of non-optimal sample fixation is smaller in this mode.Another reason is that the measurements in reflection mode are the only option for the needed in vivo measurements.
A photograph of the experimental setup used for this study is shown in Figure 3.The cornea is clamped in a specific holder that prevents it drying out.To do this, an artificial anterior chamber behind the cornea was filled with Dextran 500 soluble salt.Figure 4 shows a measurement image obtained with the polarization state HH of the illuminating source and detector, as ensured by the PSG and PSA elements.
An overexposed area can be seen in the center of the cornea.Since the cornea is not planar, there is always a part of the surface for which the angle of incidence of the exposure source is equal to the angle of reflection to the detector.From this area, a part of the light is directly reflected from the surface of the sample and causes the overexposure of the CCD camera sensor.A circular, dark area located on the left side of the cornea image stems from the opening to the rear chamber through which fluid exchange can be realized.
From the setup it is clear that the hemispheric shape of the cornea needs to be considered.The effects of polarized light incident to a curved surface can be a source of measurement error.For the following measurements, we tried to reduce the influence of the hemispheric shape by maintaining the same relative orientation for each sample.We assumed that the radius of the curvature of the cornea was approximately constant for the different samples.

Mueller Matrix Measurements Ex Vivo
Porcine samples prepared as described in Section 3.1 were used for the study presented in this work.A distinction was made between the four measurement conditions as specified in Table 2.
Table 2.The different measurement conditions specify the different treatments of the porcine cornea samples with UV irradiation and riboflavin.

Measurement Condition
Riboflavin UV Irradiation  The cornea is mounted on a specific polymer chamber to prevent it drying out which is achieved by the refilling of Dextran 500 solution.The back of the holder can also be seen as the cornea is transparent.In the central part, it becomes obvious that a region of the image is overexposed which needs to be excluded from later analysis.The black, nearly circular areas in the image show the liquid access to the chamber.The two boxes indicate from which areas of the images the used measurement data are.The blue area consists of the cornea without the reflection while the red area also contains the reflection spot.
Noticeable corneal collagen cross-linking (CXL) should only take place in the last group of samples (IV), i.e., the established procedure for photo-induced cross-linking.The remaining conditions are intended for reference and cross-checking.A total of nine cornea samples were available in the first series of measurements.Through the multiple sample use, a sample size of n = 4 can be achieved in all measurement conditions.Five cornea samples were available for the group without riboflavin exposure but with UV radiation.This is because one sample was not used several times to rule out the influence of multiple use on the measurement.The other cornea samples were measured first without UV irradiation and then with UV irradiation.Care was taken to keep the timing for the preparation of the measurements constant.This was important to minimize the influence of drying out of the cornea samples during the measurement.
The results of the measurements are displayed in Figure 5. There, the mean values obtained for each MM entry are plotted with the standard deviation.For the calculation, we used the whole images as the input.The series of measurements were grouped for samples that were not prepared with riboflavin (conditions I and II) and those that were prepared with riboflavin (conditions III and IV).
It is obvious that even for the small number of samples investigated, the MM entries for the elements M(2, 1), M(2, 2), M(2, 3), M(3, 2), M(3, 3), and M(4, 4) differed between the two groups as a consequence of CXL due to the application of riboflavin and UV irradiation, whereas they are quite similar for all other MM entries.Matching to the results of the measurement in transmission mode, these elements also show the largest mean values, with M(2, 2), M(2, 3), M(3, 2), M(3, 3) representing the entries constructed from linear polarized light states.In contrast, M(4, 4) was obtained from circular polarized light states.Especially in M(2, 2) and M (3,3), the results between the measurements with and without riboflavin differed even for the small number samples.However, the measurement data for conditions III and IV should differ from each other as well, as in the former case, UV irradiation was not applied and the effects of the cross-linking should be smaller.
As the area of the holder does not contain information about the CXL process, it can be omitted.For the calculation of the MM entries, we crop the images as shown in Figure 4.The blue area is the largest continuous rectangle that does not contain the reflection region.Thus, the maximal number of pixels can be evaluated.For the processing, it is beneficial in the current approach to evaluate rectangular areas of the same size only.If only the measurement data from the blue area is used, the entries M(2, 1), M(2, 2) and M (3,3) show a clear difference between the measurements with and without riboflavin (Figure 6a-c).If only the measurement data from the red area are taken into account, the same difference is less obvious (Figure 6d-f).
For the calculated R, P and ∆ metrics, a difference between the measurement conditions I, II and III, IV can be found for P and ∆ after a reduction in the measurement data to the certain areas specified in Figure 4.The results are displayed in Figure 7.
This indicates that the signal-to-noise-ratio can be improved by the careful choice of the region of interest segmentation for the input images and that the effect of the CXL is measurable with the MM setup.The polar decomposed MMs can show the effect of CXL after adding riboflavin and treating the samples with UV irradiation as well.This can be derived from the comparison of the matrix entries of M D and M ∆ .For the entries of M D (not displayed), the elements M(1, 2), M(1, 4), M(2, 1), M(2, 4), M(3, 3), M(3, 4), M(4, 1), M(4, 2), M(4, 3) and M (4,4) show differences between the samples with and without riboflavin treatment.Similarly to the results for M D , there is a difference between the samples with and without riboflavin treatment for the entries of M ∆ in the elements M(2, 1), M(2, 2), M(3, 3) and M(4, 4) as displayed in Figure 8.
The analysis of the decomposed MMs shows that the difference between measurement conditions I, II, and III, IV can also be derived after the polar decomposition of the MMs.As shown in [20] the degree of orientation of the samples, in this case polymer fibers, is related to the Mueller matrix and its decomposition.In the present work, the changes displayed in Figures 7 and 8 might be attributed to the changes in the orientation of the collagen fibers in the cornea.The structural changes due to the induced corneal collagen cross-linking, leading to a higher degree of orientation, might increase the polarizance and decrease the depolarization power-both of which are observable as macroscopic effects.As a practical way to quantify and to distinguish the polarization states of corneal collagen cross-linking, a Poincaré sphere is used [44].Such representation has been used for the diagnosis and classification of cancer tissues [45].To visualize the measured MMs from the four different measurement conditions, we generated 50,000 Stokes vectors which are uniformly distributed within the Poincaré sphere.Each of these Stokes vectors is multiplied with the visualized MM.The resulting Stokes vectors are plotted in the Poincaré sphere.The obtained surfaces in the Poincaré sphere are displayed in Figure 9. Every color refers to the MM of a different sample.
In Figure 9, it can be seen that the MMs occupy different spaces within the Poincaré sphere.The influence of the riboflavin is visible.The radius of the ellipsoids can be related to the degree of polarization.The larger the radii, the higher the degree of polarization is.We cannot yet separate measurement conditions III and IV, which would also be the case if the MM occupies a certain space in the Poincaré sphere solely.In further work, the standard deviation can be taken into account if the sample size can be increased.There might be room for the further optimization of the choice of the region of interest for the measurement data that is used for the calculation of MMs as well as the metrics derived from the MM.Another approach would be the calculation of another metric that enhances the contrast between measurements III and IV.The visualization within the Poincaré sphere remains a promising approach.Due to the scatter of the data and the small number of samples available, the difference between the two cases of riboflavin with and without UV irradiation is not yet obvious for the whole measurement data as well as selected measurement data.We attribute the observed behavior to the possible cross-linking from ambient UV irradiation that we cannot completely rule out.This will be studied in detail in the next steps using an improved fixation and housing setup and on a larger amount of samples.The measurement results, however, indicate that the effects of the CXL treatment, which is established in medicine for strengthening the mechanical properties of the cornea, can be investigated with the MM polarimetry, as proposed here.

Summary and Outlook
In this work, we reported MM measurements on non-planar collagen containing animal tissue, i.e., ex vivo porcine corneas.We developed a setup which allowed for the investigation of such samples in transmission and reflection mode capable of spatially resolving the inhomogeneities of the sample or the effect of photo-induced processes.The application of riboflavin followed by UV exposure along the Dresden protocol leads

Summary and Outlook
In this work, we reported MM measurements on non-planar collagen containing animal tissue, i.e., ex vivo porcine corneas.We developed a setup which allowed for the investigation of such samples in transmission and reflection mode capable of spatially resolving the inhomogeneities of the sample or the effect of photo-induced processes.The application of riboflavin followed by UV exposure along the Dresden protocol leads to noticeable changes in the MM entries indicating that these can be correlated with the photo-induced cross-linking.Furthermore, we found that the careful choice of the regionof-interest segmentation, i.e., removing overexposed areas as well as bright patches on the surface originating from impurities, leads to an improved analysis.The polarization decomposition metrics point towards changes induced by the CXL.This indicates that the results of subtle processes occurring in such media can be observed using non-invasive MM polarimetry in reflection mode.One potential application of such polarimetric measurements may address the detection of the photo-induced cross-linking of collagen fibers in the human corneas in vivo [46].This procedure, applied to increase the mechanical stability of the cornea, was employed as therapy for various diseases in ophthalmology and demands for better monitoring technologies.Furthermore, the studies of the change of collagen texture in scars to monitor wound healing appear feasible [47].However, better systematic quantification is required in the future steps to improve the insight into photo-induced processes as well as the dynamics of cross-linking.For the measurement setup and procedure, several improvements are reasonable to overcome the current limitations.For example, an improved mechanical holder can be developed including a spherically shaped quartz glass form to induce tension in the cornea sample tissue [48].In the future, the analysis of the statistical moments of the distributions of individual MM elements may reveal static and dynamic changes for the cornea samples and enhance accuracy [49].In addition, exposure bracketing can be applied to reduce overexposed areas.Also, the visualization of the MMs within the Poincaré sphere remains a promising approach for the quantitative manifestation of the physical changes of sample targets.Experimentally, the polarimetric measurements could be compared with data from Brillouin spectroscopy to better correlate the optical and mechanical properties of biotissue for further validation [50].

Figure 1 .
Figure 1.Sketch of the Mueller matrix polarimetry system.The light from the laser passes through the polarization state generator (PSG) which consists of two liquid crystal retarders (LCR) and a linear polarizer.After the interaction with the sample, the polarization states were analyzed by the polarization state analyzer (PSA) which also consists of LCRs and a linear polarizer.The LCR pairs are placed at different angles with respect to their principal axes.

Figure 2 .
Figure 2. Result of a MM measurement in the transmission mode of a porcine cornea sample at a wavelength of 532 nm.The sample has a polarization-changing effect, as obvious from the non-zero entries of the different matrix elements, essentially all, which are normalized to the M(1, 1) entry.However, the size of the effect is not homogeneously distributed over the surface of the cornea, which is evident from the individual entries and probably due to non-optimal sample fixation.

Figure 3 .
Figure 3. Photograph of the setup for MM measurements of porcine corneas in reflection mode.Since the cornea has a curved surface, the sample cannot be aligned in parallel to the illumination beam.Reproduced from[38].

Figure 4 .
Figure 4. Raw image taken with the polarization state HH for PSG (first letter) and PSA (second letter) as an example.The cornea is mounted on a specific polymer chamber to prevent it drying out which is achieved by the refilling of Dextran 500 solution.The back of the holder can also be seen as the cornea is transparent.In the central part, it becomes obvious that a region of the image is overexposed which needs to be excluded from later analysis.The black, nearly circular areas in the image show the liquid access to the chamber.The two boxes indicate from which areas of the images the used measurement data are.The blue area consists of the cornea without the reflection while the red area also contains the reflection spot.

Figure 5 .Figure 6 .Figure 7 .Figure 8 .
Figure 5. Mean values of the MM images from a series of measurements on porcine corneas.The measurements were performed on corneas neither prepared with riboflavin nor irradiated with UV (I), not prepared with riboflavin but irradiated with UV (II), prepared with riboflavin but not irradiated with UV (III) and both prepared with riboflavin and irradiated with UV (IV).The error bars represent the standard deviation for the sample size used as derived from multiple measurements.The elements M(2, 1), M(2, 2), M(2, 3), M(3, 2), M(3, 3) and M(4, 4), highlighted by bar plots show a clear difference between the samples with and without riboflavin treatment.

Figure 9 .
Figure 9. Poincaré sphere with every measurement included in the same sphere.Surfaces in the Poincaré sphere for measurement conditions I (a), II (b), III (c) and IV (d).(e) Color legend.

Figure 9 .
Figure 9. Poincaré sphere with every measurement included in the same sphere.Surfaces in the Poincaré sphere for measurement conditions I (a), II (b), III (c) and IV (d).(e) Color legend.