Ahmad, I. et al. Polarimetric assessment of healthy and radiofrequency ablated porcine myocardial tissue. J. Biophotonics 9, 750–759 (2016).

Compain, E., Poirier, S. & Drevillon, B. General and self-consistent method for the calibration of polarization modulators, polarimeters, and Mueller-matrix ellipsometers. Appl. Opt. 38, 3490–3502 (1999).

Lim, N. S. J. et al. Early detection of biomolecular changes in disrupted porcine cartilage using polarized Raman spectroscopy. J. Biomed. Opt. 16, 017003 (2011).

Fiberoptic lighting for homes

Tyo, J. S. Design of optimal polarimeters: maximization of signal-to-noise ratio and minimization of systematic error. Appl. Opt. 41, 619–630 (2002).

Jiao, S. L., Yao, G. & Wang, L. V. Depth-resolved two-dimensional Stokes vectors of backscattered light and Mueller matrices of biological tissue measured with optical coherence tomography. Appl. Opt. 39, 6318–6324 (2000).

a Time-sequenced techniques: (i) Partial Stokes vector polarimetry; the PSA can be a rotating polariser. (ii) Partial MM polarimetry; the PSG and PSA can both be a rotating polariser. (iii) Full Stokes vector polarimetry; the PSA and PSG can both be a tuneable retarder (rotating quarter waveplate or SLM assembly) followed by a fixed polariser28,63,79. (iv) Full MM polarimetry; the PSG can be a fixed polariser followed by a tuneable retarder28,63,79 (rotating quarter waveplate or LC components); the PSA can be a tuneable retarder (rotating quarter waveplate or LC components) followed by a fixed polariser28,63,79. b Snap-shot techniques: (i) Partial Stokes vector polarimetry; the PSA can be a fixed PBS assembly or circular polariser. (ii) Partial MM polarimetry; the PSA can be a fixed polariser array32. (iii) Full Stokes vector polarimetry; the types consist of division of the amplitude97,98,99, division of the wavefront71,100,101,102,103,104,105, division of the focus plane106,107,108,109,110,111,112,113,114,115,116,117,118,119,120,121,122, and so on. (iv) Full MM polarimetry; the types consist of division of wavefront72, division of the wavelength67,96, and so on. (LP: linear polariser; QWP: quarter waveplate; FDP: four-detector photopolarimeter; for more information about such Fresnel’s law based polarimetry refer to ref. 288; BS: beam splitter; PBS: polarisation beam splitter; CP: circular polariser; FQD: four quadrant detector; FQWP: four quadrant waveplate; UFPG: universal full Poincaré generator)

Ramella-Roman, J. C., Saytashev, I. & Piccini, M. A review of polarization-based imaging technologies for clinical and preclinical applications. J. Opt. 22, 123001 (2020).

Backman, V. et al. Polarized light scattering spectroscopy for quantitative measurement of epithelial cellular structures in situ. IEEE J. Sel. Top. Quantum Electron. 5, 1019–1026 (1999).

Sugita, M. et al. Retinal nerve fiber bundle tracing and analysis in human eye by polarization sensitive OCT. Biomed. Opt. Express 6, 1030–1054 (2015).

Marenko, V. & Molebnaya, T. Optimization of stokes polarimeters employing a measurement of 4intensities. Sov. J. Opt. Technol. 57, 452–455 (1990).

Li, P. C. et al. Characteristic Mueller matrices for direct assessment of the breaking of symmetries. Opt. Lett. 45, 706–709 (2020).

Ghosh, N. et al. Mueller matrix decomposition for polarized light assessment of biological tissues. J. Biophotonics 2, 145–156 (2009).

Wang, L. V., Coté, G. L. & Jacques, S. L. Special section guest editorial: tissue polarimetry. J. Biomed. Opt. 7, 278 (2002).

Haigh, J. A., Kinebas, Y. & Ramsay, A. J. Inverse conoscopy: a method to measure polarization using patterns generated by a single birefringent crystal. Appl. Opt. 53, 184–188 (2014).

Wang, X. D. & Wang, L. V. Propagation of polarized light in birefringent turbid media: a Monte Carlo study. J. Biomed. Opt. 7, 279–290 (2002).

Chen, D. S. et al. Study of optical clearing in polarization measurements by Monte Carlo simulations with anisotropic tissue-mimicking models. J. Biomed. Opt. 21, 081209 (2016).

Liu, Y. et al. Complementary fluorescence-polarization microscopy using division-of-focal-plane polarization imaging sensor. J. Biomed. Opt. 17, 116001 (2012).

Qi, J. & Elson, D. S. A high definition Mueller polarimetric endoscope for tissue characterisation. Sci. Rep. 6, 25953 (2016).

Abrahamsson, S. et al. MultiFocus polarization microscope (MF-PolScope) for 3D polarization imaging of up to 25 focal planes simultaneously. Opt. Express 23, 7734–7754 (2015).

Zhanghao, K. et al. Super-resolution dipole orientation mapping via polarization demodulation. Light.: Sci. Appl. 5, e16166 (2016).

Ossikovski, R. Analysis of depolarizing Mueller matrices through a symmetric decomposition. J. Optical Soc. Am. A 26, 1109–1118 (2009).

He, H. H. et al. A possible quantitative Mueller matrix transformation technique for anisotropic scattering media/Eine mögliche quantitative Müller-Matrix-Transformations-Technik für anisotrope streuende Medien. Photonics Lasers Med. 2, 129–137 (2013).

Numerous vectorial information measurement methods have been put forward in the past decades4,7,11,26,28,29,63. In this section, we categorise the polarisation measurement techniques into two types: time-sequenced and snap-shot approaches28,29,64,65,66,67 (see Fig. 3). For both cases, the preparation required before detection is similar and can be divided into three general steps: denoising, optimisation and calibration32,68,69,70,71,72,73,74,75,76,77,78 (see Fig. 4). The aim of those steps is to reduce the complex errors that would occur during the measurement process, hence obtaining imaging results with higher precision and accuracy32,69,70,71,72. The technical aspects of such advanced polarimetry are summarised in the review papers by Azzam28, Chipman63 and Tyo79.

Several parameters can be calculated from the Stokes vector directly (see Fig. 5a (i) and previous section): such as the degree of polarisation (DOP), degree of linear polarisation (DOLP) and degree of circular polarisation (DOCP) of light. For a single uniform light beam, the DOP is 1 for fully polarised, 0 for unpolarised or completely depolarised, and between 0 and 1 for partially polarised. The DOP cannot be larger than 1. Despite containing four elements, a Stokes vector contains fewer than four degrees of freedom due to physical constraints. The Stokes vector can also be considered as an incoherent superposition of a completely polarised part and an unpolarised part3. Those parameters have been adopted in different polarimetric applications16,20,22,23,32,80,81. The polarisation angle (PA) and intensity of the linear SOP also can be defined, with respect to dipole orientation applications153,154,155 (see Fig. 5a (iv)). For a beam generated via an incoherent light source (such as a LED), the Stokes vectors can be directly added by scalar calculation. Therefore, partially polarised light can be divided into two parts—fully polarised/depolarised components3, i.e., Stotal = Su + Sp; where Su and Sp represent fully depolarised and polarised components, respectively.

He, C. et al. Linear polarization optimized Stokes polarimeter based on four-quadrant detector. Appl. Opt. 54, 4458–4463 (2015).

Khaliq, A. et al. Comparative study of 3 × 3 Mueller matrix transformation and polar decomposition. Opt. Commun. 485, 126756 (2021).

Twietmeyer, K. M. & Chipman, R. A. Optimization of Mueller matrix polarimeters in the presence of error sources. Opt. Express 16, 11589–11603 (2008).

Ossikovski, R., De Martino, A. & Guyot, S. Forward and reverse product decompositions of depolarizing Mueller matrices. Opt. Lett. 32, 689–691 (2007).

a Stokes vector based approaches: (i) DOP, DOLP and DOCP2,3,4,5; (ii) Stokes vector projection approach158; (iii) Stokes vector location approach159; (iv) Dipole orientation differentiation approach. Polarisation angle (PA) and intensity are used as main parameters154,155. b MM based approaches: (i) MMPD method52: diattenuation (D), retardance (R) and depolarisation (Δ) (all of them maintain linear/circular components); and (ii) MMT method57: depolarisation (1-b; associated with small molecule scattering), level of linear anisotropy (t1), azimuth orientation of the anisotropy (α1), and more179; (iii) PFP method95; (iv) Property of symmetry and asymmetry168,169. (b(i)) Adapted with permission from ref. 167 © The Optical Society. (b(iii)) © [2021] IEEE. Reprinted, with permission, from ref. 95

Hsu, W. L. et al. Polarization microscope using a near infrared full-Stokes imaging polarimeter. Opt. Express 23, 4357–4368 (2015).

Several factors may contribute to depolarisation in experimental scenarios. We describe three main reasons here. (a) The first reason relates to the time domain. In general, the Stokes vector polarimeter is based on intensity measurement26, so in practice the intensity recorded at the detector includes a time-integration process. If the SOP changes rapidly, possibly due to multi-scattering induced by complex bio-media, then depolarisation would be measured. (b) This reason relates to the spatial domain. When imaging processes are involved, every point on the beam section is created through the integration of various sub-beams that could have different polarisation states. The superposition of these states leads to depolarisation. (c) The final reason is given in the spectral domain. Many processes that affect polarisation, such as birefringence and scattering, are also dependent on wavelength. Hence for different wavelengths, variations in amplitude and phase may also lead to depolarisation.

Bueno, J. M. Polarimetry using liquid-crystal variable retarders: theory and calibration. J. Opt. A: Pure Appl. Opt. 2, 216–222 (2000).

Chang, J. T. et al. Division of focal plane polarimeter-based 3×4 Mueller matrix microscope: a potential tool for quick diagnosis of human carcinoma tissues. J. Biomed. Opt. 21, 056002 (2016).

Ndagano, B. et al. Characterizing quantum channels with non-separable states of classical light. Nat. Phys. 13, 397–402 (2017).

Jacques, S. L., Ramella-Roman, J. C. & Lee, K. Imaging skin pathology with polarized light. J. Biomed. Opt. 7, 329–340 (2002).

Pierangelo, A. et al. Ex-vivo characterization of human colon cancer by Mueller polarimetric imaging. Opt. Express 19, 1582–1593 (2011).

Qi, J. & Elson, D. S. Mueller polarimetric imaging for surgical and diagnostic applications: a review. J. Biophotonics 10, 950–982 (2017).

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

The MM decomposition methods all require different assumptions (strong or weak) such as matrix reciprocity, the order that polarisation effects happen in the media, or homogeneity for the tissue analysis52,58,176,177,178. Therefore, their decomposed values are not strictly physically determined, if the assumptions do not hold in reality, which may well be the case, as biological tissue has high spatial complexity58. However, extraction through the MM polarisation parameters that have less assumptions and clearer physical meaning is always something to strive for. Several works pointed in such a direction: (1) Gil et al. and Li et al. proposed different polarisation parameters with physical determination via the asymmetric properties of the MM elements168,169, by considering assumptions about layer constructions or the presence of absence of specific vector properties such as polarisance or diattenuation; (2) Dong et al. employed a data-driven machine learning technique to fit several polarimetry feature parameters (PFPs) for characterising determined pathological applications, such as detection of the abnormal areas of breast carcinoma and cervical cancerous tissue slices95 (see Fig. 5b (iii)); 3) Breaking or restoring the symmetry (see Fig. 5b (iv)), based on analysis of different sub-regions of the MM, to extract determined information of the system is recently gaining interest179; The information extraction process is gradually developing from an analytical mathematics approach (equation-based, forward problem), to fitting or observing vectorial semantics/metrics (data-based, or shape/form-based inverse problem).

Sornsin, E. A. & Chipman, R. A. Mueller matrix polarimetry of electro-optic PLZT spatial light modulators. Proc. SPIE 2873, 196–201 (1996).

a General backscattering-mode polarimetric imaging: (i) MMs and original sample images for several ex vivo bulk specimens167. (ii) Polarisation properties of ex vivo bladder tissue256. (iii) Polarisation properties using combined wavelength information for in vivo characterisation of rat abdomen tissue92. b Certain specific backscattering-mode phantoms: (i) Full-depth MM confocal imaging of an unstained rat cornea265. (ii) PS-OCT imaging for a nerve head, which is conducted under the condition of in vivo eye imaging34. (iii) A birefringent phantom designed for bulk tissue research270. (a(ii)) Adapted with permission ref. 256. CC BY 4.0. (a(i), a(iii), b(i), b(ii) and b (iii)) Adapted with permission from refs. 34,92,167,265,270 © The Optical Society

Macdonald, C. M., Jacques, S. L. & Meglinski, I. V. Circular polarization memory in polydisperse scattering media. Phys. Rev. E 91, 033204 (2015).

Polarisation techniques can help improve the image contrast of the superficial layers of tissues by eliminating multiply scattered photons from the deep layers20,21,22,23,24. The previous literature shows that more than 85% of cancers originate from the superficial epithelium, which means that polarisation imaging methods have great potential in screening and identifying cancer at an early stage255. This would be specifically useful for in vivo clinical diagnosis, such as for minimally invasive surgery (MIS)24. Measurements in ex vivo thin tissue can use a transmissive geometry, whereas ex/in vivo bulk tissue detection would need backscattering configurations. Figure 7 gives a brief demonstration of certain current research topics related to bulk tissue polarimetry.

Gonzalez, M. et al. Design and implementation of a portable colposcope Mueller matrix polarimeter. J. Biomed. Opt. 25, 116006 (2020).

Lee, R. C. et al. Automated assessment of the remineralization of artificial enamel lesions with polarization-sensitive optical coherence tomography. Biomed. Opt. Express 5, 2950–2962 (2014).

De Boer, J. F., Hitzenberger, C. K. & Yasuno, Y. Polarization sensitive optical coherence tomography—a review. Biomed. Opt. Express 8, 1838–1873 (2017).

Azzam, R. M. A., Elminyawi, I. M. & El-Saba, A. M. General analysis and optimization of the four-detector photopolarimeter. J. Opt. Soc. Am. A 5, 681–689 (1988).

Arteaga, O. et al. Mueller matrix polarimetry with four photoelastic modulators: theory and calibration. Appl. Opt. 51, 6805–6817 (2012).

Golaraei, A. et al. Characterization of collagen in non-small cell lung carcinoma with second harmonic polarization microscopy. Biomed. Opt. Express 5, 3562–3567 (2014).

Alali, S. et al. Assessment of local structural disorders of the bladder wall in partial bladder outlet obstruction using polarized light imaging. Biomed. Opt. Express 5, 621–629 (2014).

Bargo, P. R. & Kollias, N. Measurement of skin texture through polarization imaging. Br. J. Dermatol. 162, 724–731 (2010).

Sun, M. H. et al. Characterizing the microstructures of biological tissues using Mueller matrix and transformed polarization parameters. Biomed. Opt. Express 5, 4223–4234 (2014).

The project was supported by the European Research Council (AdOMiS, no. 695140); Shenzhen Fundamental Research and Discipline Layout Project (JCYJ20170412170814624); The first author C.H. (also the figure maker) would like to thank Prof. Jessica C Ramella-Roman at Florida International University, Prof. Daniel S Elson at Imperial College London, Prof. Christoph K Hitzenberger at Medical University of Vienna, Prof. Johannes Fitzgerald de Boer at Vrije Universiteit Amsterdam, Prof. Yoshiaki Yasuno at University of Tsukuba, Prof. Martin Villiger at Harvard Medical School, Prof. Nan Zeng and Prof. Yonghong He at Tsinghua University, Dr. Long Chen at University of Oxford for useful discussions and/or support through the provision of the original graphics of their works.

Dai, H. & Yan, C. X. Measurement errors resulted from misalignment errors of the retarder in a rotating-retarder complete Stokes polarimeter. Opt. Express 22, 11869–11883 (2014).

Wood, M. F. G. et al. Polarization birefringence measurements for characterizing the myocardium, including healthy, infarcted, and stem-cell-regenerated tissues. J. Biomed. Opt. 15, 047009 (2010).

Tyo, J. S. Relation between system optimization and systematic errors in Stokes vector polarimeters. Proc. SPIE 4481, 22–30 (2002).

Scattering, especially through multiple-scattering processes, alters the degree of polarisation and SOP of the incident light beam22. While it is an insightful procedure for evaluating structural information of biomedical samples including tissues and cells16, it also introduces uncertainty in expected photon properties22. This characteristic largely hinders the development of modern tissue polarimetric techniques and related information analysis20,22,23. The turbidity of many tissue structures imposes randomness on the photons’ interaction processes, which complicates the detection and analysis of vectorial information20. Such phenomena also distinguish tissue polarimetry from the traditional polarisation measurement technique of ellipsometry22,23,24,25,26. As summarised in Fig. 1b, their comparison shows several commonalities and differences. The Jones formalism is used for clear and non-depolarising media such as thin films; it consists of the Jones vector (describing the polarisation property of the light) and Jones matrix (describing the polarisation transformation properties of the object). They have been widely used in ellipsometry techniques25,26 (see Fig. 1b; and summary in ref. 26). Another polarisation formalism is Stokes–Mueller, in which the Stokes vector and the Mueller matrix are used to describe the light beam and the object, respectively. Neither the Stokes vector nor the Mueller matrix maintain absolute phase information, but have the advantage of being able to represent depolarisation27,28. This is often essential in biomedical polarimetry, whose applications normally involve scattering induced light depolarisation20,21,22,23. There exists an increasing trend in both modern ellipsometry and polarimetry to deal with increasingly complex media, moving from isotropic and homogeneous media towards anisotropic and inhomogeneous ones20,21,22,23,24,25,26,27,28,29,30,31. While modern ellipsometry is developing towards full polarisation measurement using the Stokes–Mueller formalism, advanced polarimetry is gradually changing from full vectorial measurement to partial detection, as some key features of biomedical specimens could possibly be revealed through partial, rather than complete, measurements of vectorial information16,32,33.

Ghosh, N., Patel, H. S. & Gupta, P. K. Depolarization of light in tissue phantoms—effect of a distribution in the size of scatterers. Opt. Express 11, 2198–2205 (2003).

Azzam, R. M. A. Photopolarimetric measurement of the Mueller matrix by Fourier analysis of a single detected signal. Opt. Lett. 2, 148–150 (1978).

From the above expressions, we see that the Stokes vector can be calculated via intensity measurements that can be readily performed in an experiment47,48. The Jones vector, on the other hand, is defined by amplitude and phase that cannot be directly measured, which is another reason why the Jones approach is less well suited to biomedical polarimetry20,21,22,23. The intrinsic reason for the existence of depolarisation is due to temporal or spatial averaging16,20,21,22,23. If an extremely fast and small detector could monitor the vector properties of the light, then it would only detect polarised light. Such averaging properties can also be found in the definition of the Stokes vector47,48. Note the definitions of right-handed circular polarised light (clockwise rotation) and left-handed circular polarised light (anticlockwise rotation) are different in optics books and academic communities. It depends on whether the observer ‘sees’ the light from the source (Convention I), or from the detector (Convention II). Institute of Electrical and Electronics Engineers (IEEE) uses Convention I, so it is also widely used in engineering fields; Quantum physicists also use Convention I, to be consistent with the conventions for representing particle spin states49,50. However, for numerous optics books such as Principles of Optics (Born and Wolf48) and Handbook of Optics51, Convention II is used. In this review, we use Convention II in order to correspond to such scientific references.

The Jones vector has a graphical representation known as the polarisation ellipse47,48 (if we add the parameter DOP, the polarisation ellipse can also represent the Stokes vector (see Fig. 2b (i))). While for Stokes vector visualisation, the Poincaré sphere (PS) is commonly used47,48 (see Fig. 2b (ii)). SOPs are represented via the PS, which is defined in a three-dimensional coordinate system, whose coordinates correspond to the eigenbasis formed by \(Q\), \(U\) and \(V\) (each normalised by \(I\)). The PS is a unitary sphere that represents complete polarisation states on its surface and depolarised states inside the sphere. Any transformation of the SOP through a specimen is equivalent to manipulation of the original Stokes vector between different points on or inside the PS. Figure 2b (ii) gives a schematic demonstration of the PS. The length of the vector from the origin point to the SOP location denotes the DOP47,48. The letters \(H\), \(V\), \(M\), \(P\) are specific polarisation states: horizontally polarised (\(H\)), vertically polarised (\(V\)), 45° polarised (\(M\)) and −45° polarised (\(P\)). The polarisation ellipse parameters (χ and ψ) can be interpreted from the azimuth angle (and the polar angle) of the derived vector inside the PS.

Lu, S. Y. & Chipman, R. A. Interpretation of Mueller matrices based on polar decomposition. J. Optical Soc. Am. A 13, 1106–1113 (1996).

For polarimetric bulk tissue research, ex vivo detection plays an important role22,23,24. For example, collagen fibres, which widely exist in tissues and organs such as tendons, skin and bladder (from porcine, swine, lobster, calf or other animals156,256,257,258,259), skeletal and myocardial muscle fibres260,261,262, and elastin fibres are widely used due to their linear birefringence properties20,21,22,23,164. The alignment directions of all such fibrous structures are also linked with the fast axis orientation of the generated linear birefringence16,165,171. Furthermore, the scattering of bulk media is also studied via the extracted depolarisation16,167,263. The retardance and depolarisation related properties are the dominant parts of the vectorial properties of bulk tissues, as the magnitude of diattenuation for majority of tissue is typically very small158, with several exceptions like skeletal and myocardial muscles167 (see Fig. 7a (i)). Previous research analysing muscle tissue264 showed lower retardance compared with tendon tissue, owing to the cellularity of these tissues. Sections of the bulk myocardial fibre tissues showed two circularly aligned ring-shaped fibrous structures (see Fig. 7a (i)), revealing their anisotropic properties167. The different anisotropic vectorial information obtained from polarimetric measurements can be very helpful for the discrimination and identification of different fibrous structures in tissues164,165,167.

Liu, T. et al. Distinguishing structural features between Crohn’s disease and gastrointestinal luminal tuberculosis using Mueller matrix derived parameters. J. Biophotonics 12, e201900151 (2019).

Rivet, S., Bradu, A. & Podoleanu, A. 70 kHz full 4 x 4 Mueller polarimeter and simultaneous fiber calibration for endoscopic applications. Opt. Express 23, 23768–23786 (2015).

Note again that A is the instrument matrix for polarimetric measurement specifically, which is determined by the systematic configurations of polarisation optics and determines the error propagation amplification20,21,22,23,24,69,70,71,72; while I refers to the recorded intensity information. In Fig. 4, we take Stokes vector measurement equation (S = inv(A)·I)20,21,22,23,24,69,70,71,72 as a main illustration, to show the relationships between three steps in a picture, for simplicity. A similar structure (using the generalised equation: M = inv(A′)·I) can be derived for the MM measurement, which is based fundamentally on the Stokes vector measurement process.

Jacques, S. L., Roman, J. R. & Lee, K. Imaging superficial tissues with polarized light. Lasers Surg. Med. 26, 119–129 (2000).

Vizet, J. & Ossikovski, R. Symmetric decomposition of experimental depolarizing Mueller matrices in the degenerate case. Appl. Opt. 57, 1159–1167 (2018).

For the matrix-based calculation of Stokes polarimetry (within the scope of the above explanations), there exist two problems: first, the mathematical aspect of minimal error amplification through the matrix calculation; second, the practical aspect that the above-mentioned three separate procedures contribute to error accumulation separately, as they require different evaluation criteria and are normally based upon different assumptions. In fact, there exists the possibility to jump out of the domain of matrix calculation for Stokes polarimetry, circumventing those drawbacks. An interesting direction is the adoption of a full Poincaré beam, taking advantage of its feature that maps all SOPs in a single beam144. Vella, Zimmerman, He and others have proposed different measurement approaches harnessing such beams based on different phantoms such as stress engineered optics145,146,147,148,149,150,151 and graded index optics68. The full Poincaré beam Stokes vector technique has recently made it possible to have a clear information-based learning approach (such as the task of searching for the brightest points), combining the ‘end to end’ solution (a combination of above three processes—denoising, optimisation and calibration) together for an enhanced polarimetric measurement precision and accuracy68. In essence, this approach means that the Stokes vector retrieval process changes from matrix-based calculation to information-based image processing.

He, C. et al. Full Poincare mapping for ultra-sensitive polarimetry. Preprint at https://arxiv.org/abs/2101.09372 (2021).

Born, M. & Wolf, E. Principles of Optics: Electromagnetic Theory of Propagation, Interference and Diffraction of Light (Elsevier, 2013).

Many polarisation techniques have been harnessed for decades in biological and clinical research, each based upon measurement of the vectorial properties of light or the vectorial transformations imposed on light by objects. Various advanced vector measurement/sensing techniques, physical interpretation methods, and approaches to analyse biomedically relevant information have been developed and harnessed. In this review, we focus mainly on summarising methodologies and applications related to tissue polarimetry, with an emphasis on the adoption of the Stokes–Mueller formalism. Several recent breakthroughs, development trends, and potential multimodal uses in conjunction with other techniques are also presented. The primary goal of the review is to give the reader a general overview in the use of vectorial information that can be obtained by polarisation optics for applications in biomedical and clinical research.

Gladkova, N. et al. Evaluation of oral mucosa collagen condition with cross-polarization optical coherence tomography. J. Biophotonics 6, 321–329 (2013).

Iqbal, M. et al. Comparative study of Mueller matrix transformation and polar decomposition for optical characterization of turbid media. Optik 224, 165508 (2020).

Swami, M. K. et al. Polar decomposition of 3 × 3 Mueller matrix: a tool for quantitative tissue polarimetry. Opt. Express 14, 9324–9337 (2006).

Tyo, J. S. et al. Review of passive imaging polarimetry for remote sensing applications. Appl. Opt. 45, 5453–5469 (2006).

Du, E. et al. Two-dimensional backscattering Mueller matrix of sphere-cylinder birefringence media. J. Biomed. Opt. 17, 126016 (2012).

Donner, C. & Jensen, H. W. Light diffusion in multi-layered translucent materials. ACM Trans. Graph. 24, 1032–1039 (2005).

Both time-sequenced and snap-shot polarimetry techniques can be classified in two general ways: firstly, as either Stokes vector (light property) or MM (material property) measurement; and secondly, as partial or full vectorial measurement (Fig. 3). We will classify different techniques using the second criteria in later sections of this review.

a Light properties: intensity, wavelength, phase and polarisation. b The comparisons of the development of modern ellipsometry and polarimetry. The inset blue arrows in different colour boxes represent either changing directions or developing trends. c The structure of this review

a Several scattering types and media manipulation, showing schematically the way in which light is scattered in the differing scenarios16,17,18,19,20,21,22,23,24. b The Stokes vector representations. Top row: the polarisation ellipse and the DOP can be used to represent a given Stokes vector. Bottom row: the Poincaré sphere visualises all states of polarisation with linear states on the equator, circular states on the north/south poles and elliptical states in between2,3,4,5. c Vectorial properties that can be encoded in the Mueller matrix. Linear/circular diattenuator, linear/circular retarder etc. are fundamental polarisation elements; the arrows on the lines or circles represent the eigenbasis of modulated light beam passing through such media, where ‘before’ and ‘after’ illustrate the amplitude and/or relative phase change of the chosen eigenbasis. For further details refer to refs. 2,3,4,5

Azzam, R. M. A. Arrangement of four photodetectors for measuring the state of polarization of light. Opt. Lett. 10, 309–311 (1985).

a (i) Model of the dipole and target molecule with linker-related fundamental mechanism can be found in refs. 154,155. (ii) Fluorescence dipole orientation imaging181,182,233,289. The polarisation orientation is given as a demonstration. b (i) Model of the main properties in thin specimen for label-free polarimetric imaging. (ii) Unstained/polarimetric imaging of human cervical and liver carcinoma tissue samples32. (iii) Polarimetric parameters imaging results to distinguish between Crohn’s disease and gastrointestinal luminal tuberculosis tissues164. (iv) Polarimetric parameters images of human liver cirrhosis samples in different stages163. (a(ii) besides upper-right) Adapted with permission from refs. 181,182. CC BY 4.0. (a(ii); upper-right) Adapted from refs. 233,289. (b(ii) and b(iv)) Adapted from refs. 32,163. CC BY 4.0. (b(iii)) Copyright Wiley-VCH GmbH. Reproduced with permission ref. 164

Li, P. C. et al. Analysis of tissue microstructure with Mueller microscopy: logarithmic decomposition and Monte Carlo modeling. J. Biomed. Opt. 25, 015002 (2020).

Chan, K. H. et al. Use of 2D images of depth and integrated reflectivity to represent the severity of demineralization in cross-polarization optical coherence tomography. J. Biophotonics 8, 36–45 (2015).

Gurjar, R. S. et al. Imaging human epithelial properties with polarized light-scattering spectroscopy. Nat. Med. 7, 1245–1248 (2001).

Rapidly changing or dynamic objects need snap-shot detection, in order to correctly extract vectorial information that would be complicated by time-sequenced measurement. Snap-shot approaches are configured to take different measurements in parallel, as opposed to the serial measurement of sequential techniques. In general, snap-shot techniques must, to some degree, sacrifice alternative dimensions to enable simultaneous vector measurement67,96. Those methods include (see Fig. 3b): Stokes vector polarimeters with division-of-amplitude97,98,99, division-of-wavefront71,100,101,102,103,104,105 or division-of-focus-plane106,107,108,109,110,111,112,113,114,115,116,117,118,119,120,121,122— these fit in the category of spatial modulation with respect to different analysis channels (see Fig. 3b (iii)). Savart-plate-based polarimeters (Oka et al.) are in the category of Fourier frequency domain segmentation, which are interferometric systems where the polarisation information is encoded in the spatial carrier fringes123. If combined with the property of birefringence dispersion, spectroscopic polarimetry with channelled spectrum can also be presented67,96.

In order to understand the interactions between polarised photons and biological specimens, and link the parameters obtained via the Stokes vector or MM with the biomedical microstructural information, a software phantom—Monte Carlo (MC) simulation—was proposed to give plausible explanations for the originality of the observed physical phenomena45,46. While biomedical samples are considered as turbid media with complex structures, different fundamental units to mimic the microstructural architecture have been employed: spherical scatterers35,46; cylindrical scatterers41,46; birefringent intermedia37,38,39, multi-layered geometry45 and so on46. MC simulations have successfully reproduced most of the important polarimetric characteristic features for biomedical samples16,230,231.

Pham, H. T. T. et al. Optical parameters of human blood plasma, collagen, and calfskin based on the Stokes–Mueller technique. Appl. Opt. 57, 4353–4359 (2018).

Sosa, H. et al. ADP-induced rocking of the kinesin motor domain revealed by single-molecule fluorescence polarization microscopy. Nat. Struct. Biol. 8, 540–544 (2001).

Gil, J. J. Invariant quantities of a Mueller matrix under rotation and retarder transformations. J. Optical Soc. Am. A 33, 52–58 (2016).

Firstly, the fast development of machine learning (ML) is clearly going to have an impact on this field95,273,274. Such data-driven techniques may pave new directions for biomedical polarimetry, either through improving the quality of polarimetry (such as overcoming the numerous sources of error) or through enhanced information extraction68,95. One possibility is to use low-resolution information to reconstruct high-resolution patterns (following the spirit of works such as refs. 275,276). Secondly, while ML is geared towards improving the information processing aspects of polarimetry, new adaptive optics techniques can be used to extend the capabilities of polarimetry through full vectorial beam control. This could enable enhanced polarisation imaging resolution physically via beam shaping and compensation of polarisation errors277,278,279,280. Thirdly, the emerging techniques based on metasurfaces—subwavelength arrays of nano-scatterers that can modify polarisation—have been adopted for polarimetry281, as well as for 3D polarisation control282. Such developments may bring new opportunities for advanced biomedical polarimetry, such as forming compact vectorial sensors24,283 for deep tissue information extraction. Fourthly, second harmonic generation (SHG) and third harmonic generation (THG) based 3D MM techniques have been proposed284,285,286,287. These are described by extended MMs that are more complicated than 4 × 4 MMs used for linear scattering (4 × 9 and 4 × 16 elements, respectively, for SHG and THG)284,285. For these methods, further advanced information extraction and analysis approaches are of course intriguing. Finally, the intensity and wavelength have been utilised together with polarisation in polarimetry for a long time. However, the absolute phase information—especially geometric phase-related techniques8,10—may again open windows for new biomedical polarimetry approaches with multi-modal performance.

Hafi, N. et al. Fluorescence nanoscopy by polarization modulation and polarization angle narrowing. Nat. Methods 11, 579–584 (2014).

Zhanghao, K. et al. Super-resolution imaging of fluorescent dipoles via polarized structured illumination microscopy. Nat. Commun. 10, 4694 (2019).

For label-free biomedical polarimetric research, especially in clinical/pathological related topics, cancerous tissues detection is an important application22,23,24. In the past decades, such polarimetric techniques have assisted the diagnosis of various cancerous tissues, such as human skin cancer242, cervical cancer243,244,245,246, colon cancer166,247,248,249,250, liver cancer163,251, breast cancer and gastrointestinal cancer93,94,95,252. A typical bio-information analysis of polarimetric data is for quantitative evaluation of the fibrosis process among different stages of cancer development94,163. Beside the degree of fibrosis that can be quantified via biomedical polarimetry, the distribution of features in the fibrous regions also can serve as another characteristic parameter to assist the pathological diagnosis; this distribution can be readily extracted via polarisation information164,165,171. Intuitively, such structures contribute intrinsic birefringence mainly affecting the fourth row and fourth column of the target MM16. A good demonstration in ref. 164 shows how polarimetric textural mapping of retardance properties can distinguish between Crohn’s disease and gastrointestinal luminal tuberculosis tissues (see Fig. 6b). Some thin specimen phantoms, as found in ref. 152, target the fundamental understanding of the constitution of certain biomedical specimens, such as using nanoparticles or microspheres. Moreover, polarimetry has recently been applied to other diseases detections including Alzheimer’s disease and bladder outlet obstruction24,253,254.

Spandana, K. U., Mahato, K. K. & Mazumder, N. Polarization-resolved Stokes–Mueller imaging: a review of technology and applications. Lasers Med. Sci. 34, 1283–1293 (2019).

Compain, E. & Drevillon, B. Broadband division-of-amplitude polarimeter based on uncoated prisms. Appl. Opt. 37, 5938–5944 (1998).

Daly, S. M. & Leahy, M. J. ‘Go with the flow’: a review of methods and advancements in blood flow imaging. J. Biophotonics 6, 217–255 (2013).

Both Stokes vectors and MMs can represent the effects of time-averaged induced depolarisation16. An object may introduce two different classes of depolarisation: homogenous depolarisation and inhomogeneous depolarisation. The former one can lead to a similar DOP change for any SOP; such properties can be observed in media such as a polystyrene sphere solution. The later one can lead to different DOP change for different SOP; typical examples are found in complex biomedical tissue.

Yamanari, M. et al. Scleral birefringence as measured by polarization-sensitive optical coherence tomography and ocular biometric parameters of human eyes in vivo. Biomed. Opt. Express 5, 1391–1402 (2014).

Popescu, D. P. et al. Assessment of early demineralization in teeth using the signal attenuation in optical coherence tomography images. J. Biomed. Opt. 13, 054053 (2008).

Sase, I. et al. Axial rotation of sliding actin filaments revealed by single-fluorophore imaging. Proc. Natl Acad. Sci. USA 94, 5646–5650 (1997).

Zhanghao, K. et al. Super-resolution fluorescence polarization microscopy. J. Innovative Opt. Health Sci. 11, 1730002 (2018).

Dong, Y. et al. Probing variations of fibrous structures during the development of breast ductal carcinoma tissues via Mueller matrix imaging. Biomed. Opt. Express 11, 4960–4975 (2020).

Swami, M. K. et al. Effect of gold nanoparticles on depolarization characteristics of Intralipid tissue phantom. Opt. Lett. 38, 2855–2857 (2013).

Collett, E. Measurement of the four Stokes polarization parameters with a single circular polarizer. Opt. Commun. 52, 77–80 (1984).

Sridhar, S. & Da Silva, A. Enhanced contrast and depth resolution in polarization imaging using elliptically polarized light. J. Biomed. Opt. 21, 071107 (2016).

Kuranov, R. V. et al. Complementary use of cross-polarization and standard OCT for differential diagnosis of pathological tissues. Opt. Express 10, 707–713 (2002).

Kampmann, M. et al. Mapping the orientation of nuclear pore proteins in living cells with polarized fluorescence microscopy. Nat. Struct. Mol. Biol. 18, 643–649 (2011).

FiberOpticLightsfor Ceiling

Cruz, C. A. V. et al. Quantitative nanoscale imaging of orientational order in biological filaments by polarized superresolution microscopy. Proc. Natl Acad. Sci. USA 113, E820–E828 (2016).

Wang, W. F. et al. Roles of linear and circular polarization properties and effect of wavelength choice on differentiation between ex vivo normal and cancerous gastric samples. J. Biomed. Opt. 19, 046020 (2014).

Pierangelo, A. et al. Ex vivo photometric and polarimetric multilayer characterization of human healthy colon by multispectral Mueller imaging. J. Biomed. Opt. 17, 066009 (2012).

Qi, J. et al. Assessment of tissue polarimetric properties using Stokes polarimetric imaging with circularly polarized illumination. J. Biophotonics 11, e201700139 (2018).

Jagtap, J. et al. Quantitative Mueller matrix fluorescence spectroscopy for precancer detection. Opt. Lett. 39, 243–246 (2014).

Foreman, M. R. & Goudail, F. On the equivalence of optimization metrics in Stokes polarimetry. Opt. Eng. 58, 082410 (2019).

Ugryumova, N. et al. The collagen structure of equine articular cartilage, characterized using polarization-sensitive optical coherence tomography. J. Phys. D: Appl. Phys. 38, 2612–2619 (2005).

Zhanghao, K. et al. High-dimensional super-resolution imaging reveals heterogeneity and dynamics of subcellular lipid membranes. Nat. Commun. 11, 5890 (2020).

Before full Stokes vector/MM measurements became widely adopted, there was successful work using fixed input SOP and fixed analysers to perform partial vectorial detection for biomedical applications. Jacques et al. showed crossed-polarised light imaging to enhance surface contrast, detect skin cancer and other lesion margins80,81; Demos et al. added the dimension of wavelength based on crossed SOPs82,83; Groner et al. noted such techniques can enhance superficial vascular contrast, and hence adopt it into brain perfusion, pancreatic and further clinical diagnoses84; Bargo et al. took angle-dependency into consideration when measuring skin tissue85. Sridhar et al. also studied multiply scattered photons to enhance information extraction from biological specimens via elliptically polarised light86.

Zhao, X. J. et al. Patterned dual-layer achromatic micro-quarter-wave-retarder array for active polarization imaging. Opt. Express 22, 8024–8034 (2014).

York, T. & Gruev, V. Calibration method for division of focal plane polarimeters in the optical and near-infrared regime. Proc. SPIE 8012, 80120H (2011).

Kapsokalyvas, D. et al. In-vivo imaging of psoriatic lesions with polarization multispectral dermoscopy and multiphoton microscopy. Biomed. Opt. Express 5, 2405–2419 (2014).

Clancy, N. T. et al. Polarised stereo endoscope and narrowband detection for minimal access surgery. Biomed. Opt. Express 5, 4108–4117 (2014).

Tyo, J. S. Noise equalization in Stokes parameter images obtained by use of variable-retardance polarimeters. Opt. Lett. 25, 1198–1200 (2000).

Qiu, L. et al. Multispectral scanning during endoscopy guides biopsy of dysplasia in Barrett’s esophagus. Nat. Med. 16, 603–606 (2010).

York, T. et al. Bioinspired polarization imaging sensors: from circuits and optics to signal processing algorithms and biomedical applications. Proc. IEEE 102, 1450–1469 (2014).

For label-based direction, polarimetry has found use in scientific applications, such as biomedical microscopy16. The vectorial information of the dipole emitters is encoded in the SOP of the detected light153,154. The dipole orientation (and the fluorescence intensity) polarimetric detection technique plays an important role in thin biomedical sample analysis: e.g., in fluorescence polarisation microscopy (FPM)194,195,232,233; FPM can be used to study the nuclear pore complex subcomplexes and the relative orientations234, or be used to study different types of cytoskeleton such as actin, myosin, kinesin, microtubule and septin—those closely related with the performance of the dipole behaviours235,236,237,238,239—enabling research such as ATP and ADP binding237. Advanced research has been adopted in super-resolution imaging harnessing fluorescent dipoles via polarised illumination, with applications such as revealing heterogeneity and dynamics of subcellular lipid membranes181,240,241. These fluorescence anisotropy properties also belong to the fundamental polarisation properties that are encoded in the MM.

Forkey, J. N. et al. Three-dimensional structural dynamics of myosin V by single-molecule fluorescence polarization. Nature 422, 399–404 (2003).

Saytashev, I. et al. Self validating Mueller matrix Micro–Mesoscope (SAMMM) for the characterization of biological media. Opt. Lett. 45, 2168–2171 (2020).

Hitzenberger, C. K. et al. Measurement and imaging of birefringence and optic axis orientation by phase resolved polarization sensitive optical coherence tomography. Opt. Express 9, 780–790 (2001).

Vizet, J. et al. Optical fiber-based full Mueller polarimeter for endoscopic imaging using a two-wavelength simultaneous measurement method. J. Biomed. Opt. 21, 071106 (2016).

The structure of this review is given in Fig. 1c; it consists of introducing the basic polarisation optical tools, summarising the current vectorial information detection, extraction, and analysis approaches, and pointing out the possibilities for future multi-modal synergy with other cutting-edge technologies. Although biomedical polarimetry is still developing towards various research fields and applications, largely unexplored spaces still exist. We also hope this review could stimulate new explorations or breakthroughs in such prospective fields.

In the presence of depolarisation, Jones calculus, which represents only transitions between pure polarisation states, is of limited use as it cannot comprehensively describe the light properties, especially the degree of polarisation for partially polarised light2,3,4,5. Intrinsically, Jones calculus is based on the assumption that electric field vector holds a particular stationary state. For partially polarised light (or fully depolarised light), the variation of the electrical vector as the light propagates is semi-disordered (or completely disordered) so that more degrees of freedom are required to describe the light field47. In the scope of linear optics, the Stokes vector, which is a 4 × 1 vector, is used to characterise the SOP of the light beam47,48; while the Mueller matrix, which is a 4 × 4 matrix, describes the transformation properties of the object that affect the Stokes vector47,48. Hence, considering that the scope of this review focuses on tissue polarimetry, we place an emphasis on the Stokes–Mueller formalism.

Jiang, H. et al. Characterization of volume gratings based on distributed dielectric constant model using Mueller matrix ellipsometry. Appl. Sci. 9, 698 (2019).

It is worth noting that the use of biomedical polarimetry is expanding and has also been summarised in several recent reviews by Tuchin22, Ghosh & Vitkin20, Ramella-Roman et al.23, Qi & Elson24, De Boer et al.34, He et al.16. They have demonstrated the fast progress of this technique in the biomedical and clinical fields.

Measurement of the full vector properties of biomedical targets requires illumination with multiple SOPs in combination with multiple analysing SOPs20,21,22,23,24. As we have mentioned above, the individual MM elements lack clear physical meanings, or explicit associations with microstructures20,21,22,23,24. That is to say, vectorial characteristics of the object, like diattenuation, retardance, and depolarisation are encoded within the MM elements. For a complex optical system (like tissue), each MM element is always associated with more than one polarisation property. Hence, numerous MM decomposition methods were proposed to quantitatively characterise the optical and structural properties of the object52,53,54,55,56,57,58,59. One prevalent method is the Mueller matrix polar-decomposition (MMPD) proposed by Lu & Chipman52, which has been used and validated in lots of applications for characterisation of biomedical or material samples58,162,163,164,165,166,167 (see Fig. 5b (i)); He et al. put forward the Mueller matrix transformation (MMT) with validations using phantom experiments and Monte Carlo simulations57. Based on the MMT concept, more rotation invariant parameters were extracted from the MM and applied to biomedical sample characterisation167,168,169 (see Fig. 5b (ii)); Arteaga and colleagues derived Mueller matrix anisotropy coefficients (MMAC) to describe the level of different kinds of anisotropy for different polarisation systems56. Furthermore, other decomposition schemes were also developed, such as MM differential decomposition55, symmetric decomposition60,61 and Cloude decomposition62. Among the various decomposition approaches, different mathematical assumptions need to be made for different applications52,53,54,55,56,57,58,59,60,61,62, such as assuming a determined layer sequence of different fundamental polarisation components for a complex object, which in effect simplifies the matrix reciprocity problem52. Recently, those methods and related parameters have also been compared quantitatively with each other for the purpose of structural characterisation170,171,172. We can summarise the parameters derived via the above methods: MMPD: diattenuation (D), retardance (R) and depolarisation (Δ) (all of them maintain linear/circular components)52; MMT: depolarisation (1-b), level of linear anisotropy (t1), diattenuation property (t2), level of birefringence (t3) and fast axis orientation (x3) and more16,57; MMAC: horizontal linear anisotropy (α), 45° linear anisotropy (β) and circular anisotropy (γ) respectively with respect to the global anisotropy of the MM56.

Dubreuil, M. et al. Snapshot Mueller matrix polarimeter by wavelength polarization coding. Opt. Express 15, 13660–13668 (2007).

The Mueller matrix (MM) describes the vectorial transformation properties of an object16,20,21,22,23. As illustrated in Fig. 2c, the MM describes the transformation of one Stokes vector into another. The MM represents the full vector properties of an object through its 16 elements (\(m_{kl};k,l = 1,2,3,4\)). Among these, \(m_{11}\) represents the transformation of scalar intensity (absorption or other loss); the other 15 elements encode the vectorial properties of the object47,48. Direct physical meanings of these 15 elements taken individually are normally ambiguous16,20,21,22,23. As illustrated in Fig. 2c, several fundamental polarisation properties are encoded in (and can be extracted from) the MM. They are linear/circular diattenuation, linear/circular retardance, linear/circular polarisance, linear/circular depolarisation and so on2,3,4,5,52,53,54,55,56,57,58,59,60,61,62. The effects of each of these fundamental optical mechanisms on the light vectors along the propagation direction z are shown in Fig. 2c (i) to (vi), where ‘before’ and ‘after’ illustrate the amplitude and/or relative phase change of the chosen eigenbasis. Amongst these mechanisms, the diattenuator possesses two different absorption ratios for two polarisation directions; it in effect reduces the intensity of one polarisation compared to the other. The retarder exhibits different refractive indices for two polarised eigenvectors, in effect leading to an additional relative phase difference between the two vectors. The depolariser can modify the DOP of the light beams. For more detailed descriptions of the mechanisms and further examples see refs. 2,3,4,5.

Vella, A. J. Description and Applications of Space-variant Polarization States and Elements. PhD thesis. (University of Rochester, New York, 2018).

Schulz, M. et al. Giant intrinsic circular dichroism of prolinol-derived squaraine thin films. Nat. Commun. 9, 2413 (2018).

Rehbinder, J. et al. Ex vivo Mueller polarimetric imaging of the uterine cervix: a first statistical evaluation. J. Biomed. Opt. 21, 071113 (2016).

He, H. H. et al. Mueller matrix polarimetry—an emerging new tool for characterizing the microstructural feature of complex biological specimen. J. Lightwave Technol. 37, 2534–2548 (2019).

Lu, R. W. et al. A polarization-sensitive light field imager for multi-channel angular spectroscopy of light scattering in biological tissues. Quant. Imaging Med. Surg. 5, 1–8 (2015).

Li, X. P. et al. Polarimetric learning: a Siamese approach to learning distance metrics of algal Mueller matrix images. Appl. Opt. 57, 3829–3837 (2018).

LEDFiber Lights

Vella, A. & Alonso, M. A. Chapter Seven—maximum likelihood estimation in the context of an optical measurement. Prog. Opt. 65, 231–311 (2020).

Tariq, A. et al. Physically realizable space for the purity-depolarization plane for polarized light scattering media. Phys. Rev. Lett. 119, 033202 (2017).

Chen, Z. H. et al. A collinear reflection Mueller matrix microscope for backscattering Mueller matrix imaging. Opt. Lasers Eng. 129, 106055 (2020).

Abubakar, A. et al. A hybrid denoising algorithm of BM3D and KSVD for Gaussian noise in DoFP polarization images. IEEE Access 8, 57451–57459 (2020).

Yun, T. L. et al. Monte Carlo simulation of polarized photon scattering in anisotropic media. Opt. Express 17, 16590–16602 (2009).

Zhang, Z. G. et al. Nano-fabricated pixelated micropolarizer array for visible imaging polarimetry. Rev. Sci. Instrum. 85, 105002 (2014).

Fiberoptic lamp Vintage

Forward, S. et al. Flexible polarimetric probe for 3 × 3 Mueller matrix measurements of biological tissue. Sci. Rep. 7, 11958 (2017).

Information about the vectorial properties of a biological specimen can be derived partially from the polarisation properties of the light beam or, in a more complete fashion, from the polarisation properties of the tissue itself20,21,22,23,152. To extract information from the measured Stokes vector or MM (or part of them), different decomposition methods and parameters were proposed to represent meaningful physical processes, to extract information that could be used in subsequent analysis52,53,54,55,56,57,58,59,60,61,62.

Du, E. et al. Mueller matrix polarimetry for differentiating characteristic features of cancerous tissues. J. Biomed. Opt. 19, 076013 (2014).

De Martino, A. et al. General methods for optimized design and calibration of Mueller polarimeters. Thin Solid Films 455-456, 112–119 (2004).

Okoro, C. Second-harmonic Generation-based Mueller Matrix Polarization Analysis Of Collagen-rich Tissues. PhD thesis (University of Illinois, Urbana-Champaign, Urbana, 2018).

Ahlawat, S. et al. Polarized Raman spectroscopic investigations on hemoglobin ordering in red blood cells. J. Biomed. Opt. 19, 087002 (2014).

Soni, J. et al. Quantitative fluorescence and elastic scattering tissue polarimetry using an Eigenvalue calibrated spectroscopic Mueller matrix system. Opt. Express 21, 15475–15489 (2013).

Hu, Q. et al. Arbitrary vectorial state conversion using liquid crystal spatial light modulators. Opt. Commun. 459, 125028 (2020).

Banerjee, P. et al. Probing the fractal pattern and organization of Bacillus thuringiensis bacteria colonies growing under different conditions using quantitative spectral light scattering polarimetry. J. Biomed. Opt. 18, 035003 (2013).

Ghosh, N., Wood, M. F. G. & Vitkin, I. A. Mueller matrix decomposition for extraction of individual polarization parameters from complex turbid media exhibiting multiple scattering, optical activity, and linear birefringence. J. Biomed. Opt. 13, 044036 (2008).

DeWalt, E. L. et al. Polarization-modulated second harmonic generation ellipsometric microscopy at video rate. Anal. Chem. 86, 8448–8456 (2014).

Rosales-Guzmán, C., Ndagano, B. & Forbes, A. A review of complex vector light fields and their applications. J. Opt. 20, 123001 (2018).

De Boer, J. F. & Milner, T. E. Review of polarization sensitive optical coherence tomography and Stokes vector determination. J. Biomed. Opt. 7, 359–371 (2002).

For biomedical and clinical applications, characterising the vectorial properties of the outgoing light with a fixed incident SOP also showed great potential for structure identification22,23,24,32. Wu & Walsh reported that Stokes vector analysis with circular polarised illumination can reveal structural information about tissue156. Macdonald & Meglinski showed that turbid tissue can be quantitatively analysed via Stokes vector measurement with an optical clearing technique157. Qi et al. proposed a method of Stokes vector analysis for ex vivo porcine tongue, stomach, kidney and other tissues based on circular polarised illumination158; the most useful information was provided by circular depolarisation and linear retardance, which can normally be provided via MM decomposition52,58 (as Stokes vector projections shown in Fig. 5a (ii)). Kunnen et al. employed Stokes vector detection with circular and elliptical incident SOPs for differentiation between healthy and cancerous lung tissues specifically using a Poincaré sphere illustration159 (as Stokes vector locations shown in Fig. 5a (iii)). Note that the circular SOP illumination is especially useful for biomedical analysis, as its effects are independent of the orientation of the anisotropic components that widely exist in biomedical specimens20,21,22,23,24,32. What is more, its strong polarisation memory effect with respect to tissue-induced Mie-scattering has also gained attention160 (here the memory effect161 means that circular polarisation can survive many more scattering events than linear polarisation due to excessive forward scattering, hence it has higher probability to maintain the original information when passing through turbid tissue consisting of Mie-scattering particles that are comparable in size to the wavelength).

Guo, Y. H. et al. Study on retardance due to well-ordered birefringent cylinders in anisotropic scattering media. J. Biomed. Opt. 19, 065001 (2014).

Arteaga, O. & Canillas, A. Pseudopolar decomposition of the Jones and Mueller-Jones exponential polarization matrices. J. Optical Soc. Am. A 26, 783–793 (2009).

Rubin, N. A. et al. Matrix Fourier optics enables a compact full-Stokes polarization camera. Science 365, eaax1839 (2019).

Morio, J. & Goudail, F. Influence of the order of diattenuator, retarder, and polarizer in polar decomposition of Mueller matrices. Opt. Lett. 29, 2234–2236 (2004).

Dong, Y. et al. Quantitatively characterizing the microstructural features of breast ductal carcinoma tissues in different progression stages by Mueller matrix microscope. Biomed. Opt. Express 8, 3643–3655 (2017).

Kunnen, B. et al. Application of circularly polarized light for non-invasive diagnosis of cancerous tissues and turbid tissue-like scattering media. J. Biophotonics 8, 317–323 (2015).

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Song, B. K. et al. Broadband optical properties of graphene and HOPG investigated by spectroscopic Mueller matrix ellipsometry. Appl. Surf. Sci. 439, 1079–1087 (2018).

Marrucci, L., Manzo, C. & Paparo, D. Optical spin-to-orbital angular momentum conversion in inhomogeneous anisotropic media. Phys. Rev. Lett. 96, 163905 (2006).

Li, P. C. et al. Separating azimuthal orientation dependence in polarization measurements of anisotropic media. Opt. Express 26, 3791–3800 (2018).

Wang, X. D. & Wang, L. V. Propagation of polarized light in birefringent turbid media: time-resolved simulations. Opt. Express 9, 254–259 (2001).

Shukla, P. & Pradhan, A. Mueller decomposition images for cervical tissue: potential for discriminating normal and dysplastic states. Opt. Express 17, 1600–1609 (2009).

a Combination with machine learning techniques and big data. b Combination with vectorial adaptive optics. c Combination with metasurface based techniques. d Combination with nonlinear techniques, such as SHG and THG, in which a high order MM is required. e Combination with absolute phase information, such as geometric phase, with potential applications related with orbital angular momentum (OAM) and vector vortex beam (VVB) manipulations

Macias-Romero, C. & Török, P. Eigenvalue calibration methods for polarimetry. J. Eur. Opt. Soc. Rapid Publ. 7, 12004 (2012).

Institute of Biopharmaceutical and Health Engineering, Tsinghua Shenzhen International Graduate School, Tsinghua University, 518055, Shenzhen, China

Fu, Y. F. et al. Flexible 3×3 Mueller matrix endoscope prototype for cancer detection. IEEE Trans. Instrum. Meas. 67, 1700–1712 (2018).

Wang, L. H., Jacques, S. L. & Zheng, L. Q. MCML—Monte Carlo modeling of light transport in multi-layered tissues. Comput. Methods Prog. Biomed. 47, 131–146 (1995).

Ghosh, N. & Vitkin, I. A. Tissue polarimetry: concepts, challenges, applications, and outlook. J. Biomed. Opt. 16, 110801 (2011).

In order to reduce δA and δI32,68,69,70,71,72,73,74,75,76,77,78, a ‘denoising process’ is adopted. Figure 4 shows the approaches in time or spatial domain including time average and interpretation methods. To deal with the ΔA and ΔI, a ‘calibration process’ is required. Numerous polarimetric calibration methods have been proposed23,29,63,71,129; these can be divided into global and local calibration approaches. Note that the calibration process itself also suffers from the error transfer process. Hence, determining the SOPs for calibration, choosing the standard calibration samples, as well as designing specific calibration methods for different systems should be taken into consideration130,131. Figure 4 also shows the process of ‘optimisation’, which can deal with both types of errors, through global and local optimisation approaches. For this process, different evaluation standards have been put forward to estimate the systematic performance. Marenko et al. considered the condition number (CN) in polarimetric optimisation132, Ambirajan, Tyo and others have analysed CN-based optimisation on different phantoms70,133,134,135,136,137; and Sabatke et al. introduced equally weighted variance (EWV)69 into the polarimetric area; Azzam et al. and following researchers explained the usage of geometry optimisation based on Poincaré sphere internal volume (PSIV)138,139,140. Other useful criteria have also been proposed141,142. Such optimisation parameters can be used for evaluating the intrinsic error amplification of a polarimetry, which affect the accuracy and precision of the measurement23,29,63,71,129,130,131. If we consider the CN, the minimum CN value for a matrix-based Stokes polarimetry is \(\sqrt 3\), which is the theoretical limit for systematic error amplification68,70, as opposed to the minimum possible CN value (CN = 1) for matrix inversion. A similar error amplification also exists in MM polarimetry143. The three above-mentioned processes (denoising, optimisation, calibration) are vital for any biomedical polarimetry, as they determine the credibility of the information extraction and further analysis.

Gruev, V., Perkins, R. & York, T. CCD polarization imaging sensor with aluminum nanowire optical filters. Opt. Express 18, 19087–19094 (2010).

The MM contains fundamental physical characters like polarisance, diattenuation, retardance and depolarisation (shown in Fig. 2c); however, some concepts like anisotropy can be a combination of several fundamental polarisation processes16,57. It is worth mentioning that the depolarisation property—which is used for evaluating a SOP’s disorder, randomness, or uncertainty3,4,5—is also linked with the concept of entropy in polarimetric research11,173. While the above parameters are derived from a full MM (4 × 4); Ghosh et al. and Wang et al. also reported works on 3 × 3 MM decomposition methods, related simulations and experiments, with an emphasis on biomedical applications174,175. In summary, the decomposed linear depolarisation and linear retardance from a 3 × 3 MM display similar qualitative relationships to the changes with respect to the microstructure of the sample, such as the density, molecule size, and orientation distributions of the scatterers as well as birefringence level of the interstitial medium174,175.

Polarimetric techniques maintain unique advantages compared with other optical techniques: they can provide extra vectorial information through methods that are compatible with many existing optical systems, such as microscopes and endoscopes16,24,32,33,92,180. Much existing biomedical polarimetry research concerns sensing of bio-information in a label-free way without extraneous dyes16,22,24. In other areas, polarimetry can be used to characterise the vectorial information of fluorescence dyes, as the dipole orientation of the fluorophore is encoded in the polarisation state of the emitted light154,155. The SOP of such emission is always in a linear state; hence the polarisation angle (PA) and intensity of the linear SOP are quantities that can be harnessed, such as in biomedical applications in super-resolution microscopy153,181,182. Here we briefly summarise common phantoms used for biomedical polarimetric techniques. These techniques include polarised wide-field microscopy16,24,183, polarised light spatial frequency imaging184, polarimetric endoscopy185,186,187,188,189,190, spectral light scattering polarimetry18,82,191,192,193, polarised fluorescence spectroscopy194,195,196, polarised confocal microscopy197, polarised Raman-spectroscopy198,199, polarised super-resolution microscopy154,155, polarisation sensitive optical coherence tomography200,201,202,203,204,205,206,207,208,209,210,211,212,213,214,215,216,217,218, non-diffraction beam polarimetry (such as Bessel beam based)219, polarisation-resolved nonlinear microscopy (including second/third harmonic generation)220,221,222,223,224,225,226, and polarised speckle imaging213,227 (several techniques will be mentioned again in the Discussion). The relationship between incoherence and depolarisation of the light should be kept in mind when considering coherence based polarimetric techniques: they are different but related optical concepts. If a polarised coherent beam passing through a scattering medium becomes incoherent, it can result in either polarised light or depolarised light. If after such a medium a polarised coherent beam changes into depolarised, the coherence property may still be maintained. For more details see ref. 228. Several of the above techniques have also been adopted in three-dimensional (3D) imaging with signal integrations or sample segmentations229. However, numerous existing polarimetry techniques (within the scope of this review) fall into two-dimensional (2D) analysis23,24,25,26,27,28,29. With the completion of the cutting-edge mathematical interpretations and methodologies (see “Discussion”) there exists of course intriguing scope for further explorations.

Gao, S. K. & Gruev, V. Gradient-based interpolation method for division-of-focal-plane polarimeters. Opt. Express 21, 1137–1151 (2013).

Fan, C. M. & Yao, G. Imaging myocardial fiber orientation using polarization sensitive optical coherence tomography. Biomed. Opt. Express 4, 460–465 (2013).

Peinado, A. et al. Optimization and performance criteria of a Stokes polarimeter based on two variable retarders. Opt. Express 18, 9815–9830 (2010).

Wood, T. C. & Elson, D. S. Polarization response measurement and simulation of rigid endoscopes. Biomed. Opt. Express 1, 463–470 (2010).

Bhandari, P., Voss, K. J. & Logan, L. An instrument to measure the downwelling polarized radiance distribution in the ocean. Opt. Express 19, 17609–17620 (2011).

Ramkhalawon, R. D., Brown, T. G. & Alonso, M. A. Imaging the polarization of a light field. Opt. Express 21, 4106–4115 (2013).

Wang, J. W., Castellucci, F. & Franke-Arnold, S. Vectorial light–matter interaction: exploring spatially structured complex light fields. AVS Quantum Sci. 2, 031702 (2020).

Chen, Z. Y., Wang, X. & Liang, R. G. Calibration method of microgrid polarimeters with image interpolation. Appl. Opt. 54, 995–1001 (2015).

Peinado, A., Lizana, A. & Campos, J. Optimization and tolerance analysis of a polarimeter with ferroelectric liquid crystals. Appl. Opt. 52, 5748–5757 (2013).

Milione, G. et al. Measuring the self-healing of the spatially inhomogeneous states of polarization of vector Bessel beams. J. Opt. 17, 035617 (2015).

Alali, S., Gribble, A. & Vitkin, I. A. Rapid wide-field Mueller matrix polarimetry imaging based on four photoelastic modulators with no moving parts. Opt. Lett. 41, 1038–1041 (2016).

Dorrah, A. H. et al. Metasurface optics for on-demand polarization transformations along the optical path. Nat. Photonics 15, 287–296 (2021).

Vella, A. & Alonso, M. A. Optimal birefringence distributions for imaging polarimetry. Opt. Express 27, 36799–36814 (2019).

Macdonald, C. & Meglinski, I. Backscattering of circular polarized light from a disperse random medium influenced by optical clearing. Laser Phys. Lett. 8, 324–328 (2011).

Demos, S. G., Radousky, H. B. & Alfano, R. R. Deep subsurface imaging in tissues using spectral and polarization filtering. Opt. Express 7, 23–28 (2000).

Guangdong Engineering Center of Polarisation Imaging and Sensing Technology, Tsinghua Shenzhen International Graduate School, Tsinghua University, 518055, Shenzhen, China

\(I_0\), \(I_{90}\), \(I_{45}\), \(I_{ - 45}\) are the projection intensities (different linear components in directions of 0°, 90°, 45°, −45° with respect to the local coordinate system) of a light beam, \(I_L\) and \(I_R\) are components of left/right-handed circular polarised light, respectively. Note some other parameters can be defined with components of Stokes vector: degree of polarisation

The measurement precision and sensitivity are vital in polarimetric techniques, hence the errors need to be properly controlled28,32,68,69,70,71,72,73,74,75,76,77,78,125. However, as Stokes vector or MM measurements belong to high dimensional information detection with multiple components64,66,126, the error sources and error transfer process (such as accumulated amplifications through matrix calculations) are very complicated. Several previous analyses can be found in refs. 28,29,32,68,69,70,71,72,73,74,75,76,77,78,125,126,127,128. In Fig. 4, we summarise in one diagram an overview of the structure of the ‘denoising, optimisation, and calibration’ processes in polarimetric techniques with respect to random errors (δA, δI) and systematic errors (ΔA, ΔI). It also can be seen in the figure that three directions towards obtaining the correct vectorial measurements are still developing.

Ratliff, B. M., LaCasse, C. F. & Tyo, J. S. Interpolation strategies for reducing IFOV artifacts in microgrid polarimeter imagery. Opt. Express 17, 9112–9125 (2009).

DeMay, B. S. et al. Rapid and quantitative imaging of excitation polarized fluorescence reveals ordered septin dynamics in live yeast. Biophys. J. 101, 985–994 (2011).

He, H. H. et al. Application of sphere-cylinder scattering model to skeletal muscle. Opt. Express 18, 15104–15112 (2010).

Such a graphical representation excludes the absolute phase information, which is sometimes not addressed in typical vectorial beam analysis such as pure polarisation measurement or tissue polarimetry20,21,22,23,47,48. However, one type of the absolute phase variation that is referred to as geometric phase is related to the pathway of the SOP movement on the surface of the PS47,48, which features scope for the extension of the tissue polarimetric technique (also see “Discussion”).

Sample induced scattering is prevalent in biomedical imaging, particularly in tissues16,20,22,23. This introduces additional SOP modulations that affect diattenuation and retardance as well as depolarisation16,22,23. A complex scattering medium can often be modelled by several basic components, like spherical scatterers of different sizes35; cylindrical rod-like scatterers with different orientational distributions36; and birefringence for an interstitial medium37,38,39; combinations of these features can all be adjusted to mimic the real object36,40,41,42. Other physical conditions such as layer complexity (single-layered or multi-layered scattering43), or scattering type (elastic domain for Mie and Rayleigh scattering; or inelastic domain, like Raman scattering) are also described in the literature, e.g., see ref. 44. Modelling of the scattering assumptions can be conducted via Monte Carlo simulation45,46. This is a widely used statistical method for quantitative analysis of the interactions between polarised photons and complex biomedical media40,41,42, especially bulk media with multiple-scattering properties, for which the analytical solutions to describe the interactions cannot be obtained. In this review, we focus on the occurrence of elastic scattering in conjunction with other polarisation characteristics (see Fig. 2a) for biomedical polarimetry.

Mansfield, J. C. et al. Collagen fiber arrangement in normal and diseased cartilage studied by polarization sensitive nonlinear microscopy. J. Biomed. Opt. 13, 044020 (2008).

Tyo, J. S. Optimum linear combination strategy for an N-channel polarization-sensitive imaging or vision system. J. Opt. Soc. Am. A 15, 359–366 (1998).

Stokes polarimetry is clearly the basis for more advanced MM polarimetry. Both of their intrinsic mechanisms can be interpreted with respect to the instrument matrix (A)68,69,70,71,72 (see Fig. 4). This matrix represents the settings of the polarisation state generator (PSG) and polarisation state analyser (PSA) in the various measurement steps: for MM measurement it represents the PSG and PSA, for Stokes vector measurement it represents one the PSA. Combinations of the rotating waveplate and/or polariser are widely adopted in such approaches64,65,66. The original proposal for a Stokes vector measurement scheme (that using SOPs of \(H\), \(V\), \(M\), \(P\), \(L\) and \(R\)) was from Collett87 in 1984. Later it was adopted for biomedical information extraction or phantom analysis with ability of the full depolarisation information characterisation16,20,21,22,23.

Sugita, S. & Matsumoto, T. Quantitative measurement of the distribution and alignment of collagen fibers in unfixed aortic tissues. J. Biomech. 46, 1403–1407 (2013).

Wang, Y. et al. Mueller matrix microscope: a quantitative tool to facilitate detections and fibrosis scorings of liver cirrhosis and cancer tissues. J. Biomed. Opt. 21, 071112 (2016).

Specimens and their mimicking phantoms can be thin or bulky, which also, in general, determines the configurations of the biomedical polarimetry. A transmissive geometry is used for the thin cases (see Fig. 6) which are less scattering, thus most of the incident photons would be transmitted. A backscattering geometry (see Fig. 7) is preferred for the bulk cases (ex vivo and in vivo) which are highly scattering and depolarising, thus most of the incident photons would be backscattered. There is no clear boundary between what constitutes thin or bulk tissues. Indeed, intermediate or mixed states can exist, for which both the transmission and backscattering photons can be detected simultaneously16,22,23,152. In response to the beginning of this section, biomedical polarimetry can be used in labelled or label-free measurement; Fig. 6 gives a summary for two types of the use of thin tissue polarimetry.

Axelrod, D. Carbocyanine dye orientation in red cell membrane studied by microscopic fluorescence polarization. Biophys. J. 26, 557–573 (1979).

Dong, Y. et al. Deriving polarimetry feature parameters to characterize microstructural features in histological sections of breast tissues. IEEE Trans. Biomed. Eng. 68, 881–892 (2021).

Dai, Y. Y. et al. Active compensation of extrinsic polarization errors using adaptive optics. Opt. Express 27, 35797–35810 (2019).

He, H. H. et al. Monitoring microstructural variations of fresh skeletal muscle tissues by Mueller matrix imaging. J. Biophotonics 10, 664–673 (2017).

Mendoza-Galván, A. et al. Mueller matrix spectroscopic ellipsometry study of chiral nanocrystalline cellulose films. J. Opt. 20, 024001 (2018).

Light, as an electromagnetic wave, possesses several fundamental properties, which include intensity, wavelength, phase and polarisation1,2 (see Fig. 1a). While the former three are scalar quantities, polarisation has vectorial properties; its use has therefore required more advanced optical and computational approaches. Hence, studies of either the vector properties of light, described via the state of polarisation (SOP) or the full vectorial transformation properties of an object, have a shorter history in biomedical analysis compared with their scalar counterparts, and the extent of their application is still being explored3,4,5. So far, numerous intriguing areas of research have been enhanced through harnessing vectorial information acquired via polarisation optics; these range from fundamental research6,7,8,9,10, such as quantified polarisation entropy11, across quantum physics12, such as spin-orbital interaction of light13,14, to material characterisation (e.g. chiral characteristics15) or for biomedical studies and clinical applications (e.g. characterisation of structural features in tissue16,17,18,19,20,21).

Hu, Q., He, C. & Booth, M. J. Arbitrary complex retarders using a sequence of spatial light modulators as the basis for adaptive polarisation compensation. J. Opt. 23, 065602 (2021).

While thin samples can feature multiple-scattering process, such processes are of course more significant in bulk samples20,21,22. Considering the intriguing scope of the polarimetric technique for in vivo clinical diagnosis, beside the simulations, there is a need for complex phantoms, such as those exhibiting birefringence (see Fig. 7b (iii)) or depolarisation, to establish reliable processes for investigation of complex scattering mechanisms270,271. A recent review has summarised various phantoms for both thin and bulk samples152. Microspheres, silicon-based phantoms, nanoparticles, cylindrical scatterers, and birefringent/dichroism films35,158,270,272 have all been employed in various validations. In vivo biomedical polarimetry and its related applications clearly offer a large space for future exploration.

Ahmad, I. et al. Ex vivo characterization of normal and adenocarcinoma colon samples by Mueller matrix polarimetry. J. Biomed. Opt. 20, 056012 (2015).

Ahmad, J. E. & Takakura, Y. Error analysis for rotating active Stokes–Mueller imaging polarimeters. Opt. Lett. 31, 2858–2860 (2006).

Foreman, M. R., Favaro, A. & Aiello, A. Optimal frames for polarization state reconstruction. Phys. Rev. Lett. 115, 263901 (2015).

The use of rotating components has disadvantages, such as increasing measurement time and introducing unexpected errors from mechanical movements. However, such systems are easy to construct. Hence, numerous commercialised polarimeters still use this approach. In order to make improvements, researchers have tried to reduce the number of the rotating components (such as the dual-rotating waveplate MM polarimeter with fixed polarisers that was proposed by Azzam64 in 1978, which is widely used in tissue analysis16,22) or use fast modulation devices (such as Stokes or MM polarimeters enabled via liquid crystal variable retarders (LCVR)88, spatial light modulators (SLM)89, ferroelectric liquid crystals (FLC)90, or photoelastic modulators (PEM)91). Besides full MM detection, partial MM measurement, such as 3 × 3 MM imaging of linear polarisation states, also gained wide attention. Qi et al. used related methods in analysing linear depolarisation and retardance of rat tissue92.

Tyo, J. S., LaCasse, C. F. & Ratliff, B. M. Total elimination of sampling errors in polarization imagery obtained with integrated microgrid polarimeters. Opt. Lett. 34, 3187–3189 (2009).

Wang, S. & Larin, K. V. Optical coherence elastography for tissue characterization: a review. J. Biophotonics 8, 279–302 (2015).

Beckley, A. M. Polarimetry and Beam Apodization Using Stress-engineered Optical Elements. PhD thesis (University of Rochester, New York, 2012).

Manhas, S. et al. Demonstration of full 4 × 4 Mueller polarimetry through an optical fiber for endoscopic applications. Opt. Express 23, 3047–3054 (2015).

Mu, T. K. et al. Error analysis of single-snapshot full-Stokes division-of-aperture imaging polarimeters. Opt. Express 23, 10822–10835 (2015).

Similar to Stokes vector polarimeters, there exist concepts for snap-shot MM polarimeters, in which certain dimensions are sacrificed to enable simultaneous MM estimation (see Fig. 3b (iv)). Dubreuil et al.67 and Hagen et al. 96 utilised different wavelength-dependent birefringent media to resolve the MM in a single shot, within the limitation of the sample being achromatic. Piquero et al. 124 utilised full Poincaré beams as a PSG, enabling MM polarimetry with the division of the wavefront. He et al. applied a spatially segmented method with defocusing to measure statistically averaged properties of biomedical samples72. As complete snap-shot MM techniques are rather complex, their usage for extractions of biomedical information is less common than the use of single-shot partial MM or Stokes vector polarimetry. For instance, 3 × 4 MM imaging is also gaining attention using circularly polarised illumination32; Chang et al. brought such a technique into the human liver and cervical carcinoma tissue analysis32 (see Fig. 3b (ii)).

Dong, Y. et al. A quantitative and non-contact technique to characterise microstructural variations of skin tissues during photo-damaging process based on Mueller matrix polarimetry. Sci. Rep. 7, 14702 (2017).

Tanaka, Y. et al. Motion-artifact-robust, polarization-resolved second-harmonic-generation microscopy based on rapid polarization switching with electro-optic Pockells cell and its application to in vivo visualization of collagen fiber orientation in human facial skin. Biomed. Opt. Express 5, 1099–1113 (2014).

Hong, Y. J. et al. Optically buffered Jones-matrix-based multifunctional optical coherence tomography with polarization mode dispersion correction. Biomed. Opt. Express 6, 225–243 (2015).

Arteaga, O., Garcia-Caurel, E. & Ossikovski, R. Anisotropy coefficients of a Mueller matrix. J. Optical Soc. Am. A 28, 548–553 (2011).

Gan, Y. & Fleming, C. P. Extracting three-dimensional orientation and tractography of myofibers using optical coherence tomography. Biomed. Opt. Express 4, 2150–2165 (2013).

Although there are some applications that require high-speed operation, such as detection in dynamic situations like in vivo sensing for clinical diagnosis24, time-sequenced polarimeters still play an important role in modern polarimetric research, due to their mature state of development and simple configuration. Such applications include characterisation of complex vector fields6,7,9, or providing ground truth validation in tissue research (e.g., differentiating human breast cancer93,94,95).

Sun, T. et al. Distinguishing anisotropy orientations originated from scattering and birefringence of turbid media using Mueller matrix derived parameters. Opt. Lett. 43, 4092–4095 (2018).

Pierangelo, A. et al. Multispectral Mueller polarimetric imaging detecting residual cancer and cancer regression after neoadjuvant treatment for colorectal carcinomas. J. Biomed. Opt. 18, 046014 (2013).

He, C., He, H., Chang, J. et al. Polarisation optics for biomedical and clinical applications: a review. Light Sci Appl 10, 194 (2021). https://doi.org/10.1038/s41377-021-00639-x

Novikova, T. et al. The origins of polarimetric image contrast between healthy and cancerous human colon tissue. Appl. Phys. Lett. 102, 241103 (2013).

He, C. et al. Vectorial adaptive optics: correction of polarization and phase. in Adaptive Optics and Wavefront Control for Biological Systems VI. Vol. 11248 (eds Bifano, T. G., Gigan, S., & Ji, N.) (International Society for Optics and Photonics, 2020). https://www.spiedigitallibrary.org/conference-proceedings-of-spie/11248/1124808/Vectorial-adaptive-opticscorrection-of-polarization-and-phase/10.1117/12.2547715.full?SSO=1.

Azzam, R. M. A. Division-of-amplitude photopolarimeter (DOAP) for the simultaneous measurement of all four Stokes parameters of light. Opt. Acta.: Int. J. Opt. 29, 685–689 (1982).

Gao, S. K. & Gruev, V. Bilinear and bicubic interpolation methods for division of focal plane polarimeters. Opt. Express 19, 26161–26173 (2011).

Heinrich, C. et al. Mueller polarimetric imaging of biological tissues: classification in a decision-theoretic framework. J. Opt. Soc. Am. A 35, 2046–2057 (2018).

Shen, Y. X. et al. Comparative study of the influence of imaging resolution on linear retardance parameters derived from the Mueller matrix. Biomed. Opt. Express 12, 211–225 (2021).

FiberOptic Light Price

While ex vivo studies are mainly oriented towards fundamental research24,152,171 (e.g., understanding the vectorial properties characterisation; see Fig. 7a (i), (ii), b (i)), in vivo bulk tissue polarimetry is geared towards applications24,256,265. Typical backscattering-mode polarimetry includes polarisation endoscopy24, reflection MM microscopy16,266, MM colposcopy180,267, wide-field handheld polarimetry16, and PS-OCT268 (see Fig. 7a (iii), b (ii)), targeted to clinical diagnosis in vivo. As a promising in vivo, label-free diagnostic tool, polarisation endoscopic imaging has been implemented inside rat abdomen, revealing the small bowel, stomach, liver and fat with different polarisation characters92. Recent work also includes development of several different types of MM endoscope188,189,190,256, and extension into the spectral domain (with certain fixed wavelengths)92 (see Fig. 7a (iii)). PS-OCT200,201,202,203,204,205,206,207,208,209,210,211,212,213,214,215,216,217,218 is specifically used for in vivo ophthalmic imaging, where polarimetric data accompanied with clinical analysis has been demonstrated, for retinal imaging268 (see Fig. 7b (ii)). Other types of bulk tissue analysis, such as human lung cancerous tissue159 and skin tissue269, show good prospects for future clinical diagnosis16,22,23,24.

Brasselet, S. Polarization-resolved nonlinear microscopy: application to structural molecular and biological imaging. Adv. Opt. Photonics 3, 205–271 (2011).

Biomedical applications of polarimetry have attracted substantial attention. We hope this short review paper gives readers a general overview from fundamental polarisation concepts, through polarimetric techniques, to recent biomedical and clinical applications7,16,20,21,22,23,24,29,34,63,79,152. In addition to the summaries of recent research trends explained above, we provide here some further perspective on prospects in this application area, considering the use of polarimetry in a multimodal combination with other advanced technologies (see Fig. 8 for a summary).

Ghosh, N., Wood, M. F. G. & Vitkin, I. A. Influence of the order of the constituent basis matrices on the Mueller matrix decomposition-derived polarization parameters in complex turbid media such as biological tissues. Opt. Commun. 283, 1200–1208 (2010).

FiberOpticLightsfor Crafts

Chen, D. S. et al. Mueller matrix polarimetry for characterizing microstructural variation of nude mouse skin during tissue optical clearing. Biomed. Opt. Express 8, 3559–3570 (2017).

Chang, J. T. et al. Single-shot spatially modulated Stokes polarimeter based on a GRIN lens. Opt. Lett. 39, 2656–2659 (2014).

Mehta, S. B. et al. Dissection of molecular assembly dynamics by tracking orientation and position of single molecules in live cells. Proc. Natl Acad. Sci. USA 113, E6352–E6361 (2016).

Andreou, A. G. & Kalayjian, Z. K. Polarization imaging: principles and integrated polarimeters. IEEE Sens. J. 2, 566–576 (2002).

Yao, G. & Wang, L. V. Two-dimensional depth-resolved Mueller matrix characterization of biological tissue by optical coherence tomography. Opt. Lett. 24, 537–539 (1999).

Wood, M. F. et al. Proof-of-principle demonstration of a Mueller matrix decomposition method for polarized light tissue characterization in vivo. J. Biomed. Opt. 14, 014029 (2009).

Wu, P. J. & Walsh, J. T. Jr. Stokes polarimetry imaging of rat-tail tissue in a turbid medium using incident circularly polarized light. Lasers Surg. Med. 37, 396–406 (2005).

Novikova, T. et al. Special section guest editorial: polarized light for biomedical applications. J. Biomed. Opt. 21, 071001 (2016).

Borovkova, M. et al. Evaluating β-amyloidosis progression in Alzheimer’s disease with Mueller polarimetry. Biomed. Opt. Express 11, 4509–4519 (2020).

Yang, B. et al. Polarized light spatial frequency domain imaging for non-destructive quantification of soft tissue fibrous structures. Biomed. Opt. Express 6, 1520–1533 (2015).

Chue-Sang, J. et al. Use of combined polarization-sensitive optical coherence tomography and Mueller matrix imaging for the polarimetric characterization of excised biological tissue. J. Biomed. Opt. 21, 071109 (2016).

Zotter, S. et al. Measuring retinal nerve fiber layer birefringence, retardation, and thickness using wide-field, high-speed polarization sensitive spectral domain OCT. Investigative Ophthalmol. Vis. Sci. 54, 72–84 (2013).

Satapathi, S., Soni, J. & Ghosh, N. Fluorescent Mueller matrix analysis of a highly scattering turbid media. Appl. Phys. Lett. 104, 131902 (2014).

Chen, L. et al. Advances of super-resolution fluorescence polarization microscopy and its applications in life sciences. Comput. Struct. Biotechnol. J. 18, 2209–2216 (2020).

Wang, Y. et al. Differentiating characteristic microstructural features of cancerous tissues using Mueller matrix microscope. Micron 79, 8–15 (2015).

DeMay, B. S. et al. Septin filaments exhibit a dynamic, paired organization that is conserved from yeast to mammals. J. Cell Biol. 193, 1065–1081 (2011).

Groner, W. et al. Orthogonal polarization spectral imaging: a new method for study of the microcirculation. Nat. Med. 5, 1209–1212 (1999).