Search tips
Search criteria 


Logo of boeaboutauthor infoeditorial boardsearchOptics InfoBaseboeThis article
Biomed Opt Express. 2010 September 1; 1(2): 463–470.
Published online 2010 August 3. doi:  10.1364/BOE.1.000463
PMCID: PMC3018013

Polarization response measurement and simulation of rigid endoscopes


Polarized light can reveal diagnostic information about tissue morphology. To promote easy adoption of polarization imaging techniques in the clinic it would be beneficial if they can be used with standard medical imaging instruments such as rigid endoscopes. We have characterized the polarization properties of two commercial laparoscopes and observed birefringence effects that complicate polarization imaging. Possible solutions are discussed that may be of interest to other tissue polarization imaging researchers.

OCIS codes: (170.0110) Imaging systems, (170.2150) Endoscopic imaging

1. Introduction

Several studies have shown that the interaction of polarized light with tissue can reveal additional information compared to standard unpolarized illumination. For instance simply viewing tissue under crossed polarizers will reduce specular highlights and increase the contrast between tissues [1], whereas fully characterizing the interaction of polarized light with tissue can reveal changes in its bulk structure caused by cancer or necrosis [2–5]. Light-scattering spectroscopy uses Mie-theory calculations to recover the average size of cell nuclei, which changes as tumours progress [6,7]. Polarization techniques provide morphological information and hence complement fluorescence, Raman or absorption sensing which reveal functional or metabolic information. Imaging these properties over a wide area will give clinicians a tool for assessing the boundaries between different tissue states. This is currently a complicated clinical task that requires a thorough visual examination, use of pre-operative medical imaging data or intra-operative biopsies depending on the procedure.

As with many biophotonics technologies, translating polarization imaging systems from the lab to the clinic raises many practical issues, partly because many current systems use a transmission geometry where the sample is located directly between the illumination source and the camera. Whilst useful for ex-vivo samples or small animal imaging this is not feasible for imaging organs in-vivo. Simple reflection geometries [3] are suitable for imaging skin but cannot image inside the body except during open surgery. Minimally invasive endoscopic techniques have been widely adopted in modern medicine due to better patient outcomes. The simplest way to bring polarization techniques into the clinic is to incorporate them into current endoscopic surgery instruments. We have previously used a rigid laparoscope to perform polarization spectroscopy, but after noticing interference fringes we have carried out a full polarization characterization of two commercial laparoscopes, one of which is shown in Fig. 1.

Fig. 1.
The Olympus 0° forward-viewing laparoscope characterized in this experiment (top), and a schematic showing a cross section (bottom). The length is 500 mm and diameter 10 mm. The proximal end is on the left and distal tip the right. Illumination ...

2. Methods

2.1. Rigid Endoscopes

We briefly outline the construction of rigid endoscopes in order to describe the polarization properties observed. Rigid endoscopes have superior image quality compared to flexible endoscopes as the rod lens relay system gives better light transmission than a fibre bundle. They are commonly used for procedures inside the abdominal cavity where they are called laparoscopes. We obtained two laparoscopes from prominent manufacturers, one from Karl Storz GmbH (10 mm 0° viewing angle, part no. 26003 AA) and one from Olympus Ltd (10 mm 0°, part no. A5254A). Similar considerations described in this paper apply to other rigid endoscopes, e.g. cystoscopes and arthroscopes as the fundamental design is the same.

The optics of a laparoscope can be divided into an objective lens, the rod lens system and an eyepiece. A schematic representation is shown in Fig. 1. The objective lens at the distal end demagnifies the incident image so it can be passed into the relay system. The relay system comprises several sets of rod lenses which allow parallel rays from image points to travel with shallow angles. This increases the distance a ray travels inside the lens without intersection with the laparoscope tube and reduces the number of lenses required. Finally the proximal eyepiece creates approximately parallel rays out of the laparoscope so the eye or another lens can bring them to focus. Both the distal and proximal lenses are covered by windows that can withstand the autoclaving process required to sterilize surgical instruments.

Illumination light is transmitted to the distal via a built-in fibre bundle with a connector at the proximal end as shown in Fig. 1. These fibres are not polarization preserving, hence to provide polarized illumination either these fibres must be replaced or a small polarizer must be mounted at the distal tip.

2.2. Measuring Polarization States

A Stokes’ vector S describes a polarization state using four components which are generally labelled I, Q, U, and V. These describe the illumination intensity, and then the difference between horizontal and vertical (H & V), +45° and -45° (P & M), and right- and left-hand circular (R & L) polarization respectively, see Eq. (1). A Stokes’ vector can be calculated from just six measurements using a standard analyser, as the total illumination intensity is just the sum of any two orthogonal states (conventionally the horizontal and vertical states are used) [8].

The polarization effects of an optical component can then be described by a four-by-four Mueller matrix M (1) which exhaustively details how each output polarization state depends on the input polarization state. To calculate the Stokes’ vector of a ray exiting a component the input vector is simply pre-multiplied by the Mueller matrix as in Eq. (2) . Hence an identity Mueller matrix implies that the component does not interact with polarized light at all, and the left-most column gives the output state for input unpolarized light. In order to record a complete Mueller matrix 36 separate measurements are required [9]. This can be reduced to 18 by using some algebraic manipulation but we have chosen to record the full 36 measurements in order to average out noise [3].


A matrix was recorded for each pixel in the field of view using a CCD camera and the set-up shown in Fig. 2. This used rotatable sheets of linear polarizer (TechSpec, Edmund Optics) and removable λ/4 waveplates (CVI-Melles Griot) to create the relevant polarizer and analyser orientations. Malus’ law was used to establish a polarization reference frame consistent with the laboratory and hence the CCD camera. Laser light was reflected obliquely from a slide onto a power meter through the rotatable analyser that was rotated until the incident illumination was at a minimum. This was labelled as the ‘Vertical’ polarization direction with ‘Horizontal’ perpendicular. The rotatable polarizer was then aligned to the analyser. Once the reference frame had been established the polarizer and analyser were placed at either end of the laparoscope with sufficient space to insert and remove the waveplates as required. The CCD and lens were placed behind the analyser so that the whole field of view was brought into focus. A diffuser was placed behind the polarizer to create a more even illumination field.

Fig. 2.
Experimental schematic. Laser light at 600 nm is incident on a diffuser to ensure even illumination. A rotatable polarizer and removable λ/4 waveplate are used to create linear and circular polarization states. This uniform polarization is then ...

The laser light was provided by a spectrally filtered supercontinuum laser source (Fianium SC400) [10]. The white light output of the laser was dispersed by a prism and focused onto a Digital Multi-Mirror Device (DMD, Texas Instruments Ltd.). Portions of the spectrum were selected by switching the DMD mirrors to reflect specific wavelengths through a second focusing lens into a liquid light guide (Karl Storz Ltd.). Liquid light guides are used widely to connect lamp sources to endoscopes instead of fibre-optic light guides as their saline cores are more flexible and durable.

For all but the most trivial optical interactions Mueller matrices can be hard to interpret. This effect is magnified if a single matrix contains the compound action of several components. Lu and Chipman previously showed that is possible to decompose a non-degenerate Mueller matrix into a set of matrices each representing a single optical effect [11]. From these Mueller matrices a single value can be calculated representing the depolarization (Δ), retardance (ψ) and diattenuation (D) effects [12–15] as in Eq. (3). Depolarization is the decrease in the degree of polarization defined as V=Q2+U2+V2/I due to effects such as scattering. The retardance is the additional phase difference introduced between two orthogonal polarization states. This effect is used in λ/4 waveplates to convert between linear and circular polarization. Diattenuation describes the difference in transmission between two orthogonal states. A linear polarizer is a perfect diattenuator as it completely passes one state while blocking the orthogonal state.


3. Results and Discussion

3.1. Measured Results

The measured Mueller matrices for the laparoscopes are shown in Fig. 3. Each sub-image shows one element of the Mueller matrix across the circular image field of the endoscope. To correct for the illumination intensity the elements are normalized to element m11. Instead of showing the resulting unity image for m11, we have normalized to its own maximum value in the figure to illustrate the illumination field. This shows a radial fall-off, as usually observed in endoscopes due to the uneven illumination and image apodization.

Fig. 3.
Measured Mueller matrices for the Olympus and Storz laparoscope (a)–(c) and a simulation for a sheet of sapphire (d). Each sub-image shows one element of the matrix across the whole field. Matrices were measured with the system illustrated in ...

Uniform polarization effects across the field of view would result in elements of uniform value and hence flat images displaying no structure. This is clearly not the case as highly structured patterns are visible across the field of view, meaning the polarization properties change with image co-ordinate. The Olympus endoscope showed circular arcs in elements m22 to m43, which were observed to rotate and change value as the endoscope was rotated about its long axis [compare Figs. 3(a) and 3(b)]. The Karl Storz endoscope showed patterns with high amounts of circular symmetry which did not change if it is rotated (hence only one orientation is shown). Both patterns exhibited a simple wavelength dependence, expanding outward as the wavelength increased (Not shown).

Both laparoscopes shared the property that elements m12 to m14 (top-row) and m21 to m41 (left-column) were zero. From this we can deduce that they do not contain polarizing elements and only interact with light that is already polarized. Therefore these effects are not apparent when illuminated using standard sources such as halogen lamps.

3.2. Simulated Results

In order to better understand where these effects may originate we considered the arrangement of rays passing through the entrance and exit windows of the laparoscopes as shown in Fig. 2. Rays originating from an image point are aligned approximately parallel to each other when they pass through the exit window with points at the edge of the image passing at greater angles. This is equivalent to the conoscopic geometry used in crystallography [17] and illustrated in Fig. 4. The important feature of this arrangement is that it directly interchanges image coordinates for angles passing through a sample. This allows the birefringence of a sample to be easily assessed at all angles, and any symmetry properties of the resulting conoscopic images are related to the orientation and symmetry properties of the crystal lattice. In our case the exit window is in an equivalent plane to the crystal sample.

Fig. 4.
The standard conoscopic geometry. Points in the image correspond to angles through the sample, usually a crystal. Any symmetry in the captured images will correspond to symmetries in the crystal lattice of the sample. These will change with the orientation ...

The windows of the Karl Storz laparoscope are made from sapphire which is a birefringent crystal. To confirm that the observed patterns could be explained by this we simulated the Mueller matrix that would be expected in a conoscopic geometry if a thin sheet of sapphire had its slow axis aligned parallel to the optic axis. The half-angle of the field of view for the laparoscope was calculated as approximately 38.6°. A grid of angles matching this was created and used to find the effective refractive indices for the ordinary and extra-ordinary rays, which were subtracted to find the effective birefringence at a wavelength of 600 nm. The resulting Jones matrix was calculated using an assumed thickness of 0.5 mm for the sapphire sheet at all incidence angles. These were pre- and post-multiplied by the Jones vectors corresponding to the measurement polarizations to produce simulated images of a conoscopy experiment. The same technique as for the real experiment was then used to convert these into a Mueller matrix. The Jones calculus was used as it simplifies calculations where only fully polarized light is used. The simulation is presented in Fig. 3(d) which shows closely matching patterns for the top left 3×3 sub-matrix with Fig. 3(c), however they do not match in the right-most column and bottom row (the elements of a Mueller matrix relating to circular polarizations).

This indicated that the λ/4 waveplates might cause the discrepancy between measurement and simulation as they are only used to measure these elements. Achromatic waveplates have an angular dependence on their retardance and so are not suited to imaging applications. The angular dependence means that at points away from the centre we are no longer dealing with purely circular polarization but some elliptical state. The Karl Storz matrix was decomposed using Lu and Chipman’s technique, shown in Fig. 5. The diattenuation and depolarization were approximately zero across the field of view, as expected. The retardance exhibited angular symmetry with a sinusoidal radial profile. However all of the parameters had X-shaped patterns where they were discontinuous. This effect is strongest at the distal end of the laparoscope as the rays from the object crossing the waveplate are at steeper angles than those at the proximal end.

Fig. 5.
The decomposed parameters of the Storz laparoscope (a)–(c) and the calculated retardance of a magnesium fluoride plate (d) for similar viewing angles, showing the same X-shaped patterns in both. Adapted from [16].

We then simulated the retardance through a sheet of magnesium fluoride, a material used to make achromatic waveplates using the same technique as for sapphire above. The result is shown in Fig. 5(d) and demonstrates the same X-shaped zones. This implies that the waveplates do affect the measurement of the Mueller matrix and hence the mismatch with simulation. Using zero-order waveplates would reduce the angular dependence although they would introduce a wavelength dependence, meaning the results would only be correct close to the central wavelength of the plate.

The close match with the simulation of the sapphire sheet shows that the polarization effects of the Karl Storz laparoscope come primarily from the birefrigent window material. We are currently unable to provide a similarly convincing explanation of the behaviour of the Olympus laparoscope. We are unaware of the exact window material used but the lack of any symmetry means that the windows cannot be made from a birefrigent crystal oriented as a λ/4 or λ/2 waveplate, although the orientation dependence and interaction with circularly polarized light does display waveplate-like properties.

3.3. Possible Solutions

In order to use commercial endoscopes with polarized imaging any polarization effects must be calibrated and corrected. In theory this can be done simply by measuring the Mueller matrix of a particular laparoscope before use in theatre, and provided this matrix is not degenerate its inverse could be post-multiplied by any Stokes’ measurements made through the laparoscope to extract the polarization state at the distal tip. When no circular polarizations are present in a system it is possible to use just a three-by-three Mueller matrix and consider only the linear states. However the birefringent nature of the crystals used means that most rays exiting the endoscope will have a circular component, and this cannot just be ignored. Hence the full Mueller matrix would have to be measured for every pixel in the image and at all wavelengths used in the system. But as described above achromatic waveplates have a strong angular dependency, whereas zero-order waveplates have a strong wavelength dependency. This means that the circular components of the Mueller matrix cannot be measured accurately for all pixels and wavelengths.

A simpler way of removing the polariztion effects would be to remove the birefrigent crystals and replace them with non-birefrigent alternatives. Unfortunately simple fused silica is unsuitable since the bonding agents used with it cannot withstand the autoclave sterilization process. A good alternative would be diamond, as this uses the same bonding materials as sapphire but is not birefringent.

4. Conclusions

Current commercial laparoscopes are difficult to use with polarized light due to the birefringent materials used for the entrance and exit windows. Limitations in the characteristics of the waveplates necessary to measure circularly polarized light complicates the calibration of the polarization response using a Mueller matrix. This will necessitate the use of new, non-birefrigent rigid endoscopes in future work, allowing the simple measurement of polarization data in-vivo using an instrument already familiar to surgeons.


We would like to acknowledge Karl Storz GmbH for the loan of an endoscope system and the UK EPSRC and Technology Strategy Board grants DT/F003064/1 and EP/E06342X/1.

References and links

1. Anderson R. R., “Polarized light examination and photography of the skin,” Arch. Dermatol. 127(7), 1000–1005 (1991).10.1001/archderm.127.7.1000 [PubMed] [Cross Ref]
2. Sun C. W., Lu L. S., Yang C. C., Kiang Y. W., Su M. J., “Myocardial tissue characterization based on the time-resolved Stokes-Mueller formalism,” Opt. Express 10(23), 1347–1353 (2002). [PubMed]
3. Baba J. S., Chung J. R., DeLaughter A. H., Cameron B. D., Coté G. L., “Development and calibration of an automated Mueller matrix polarization imaging system,” J. Biomed. Opt. 7(3), 341–349 (2002).10.1117/1.1486248 [PubMed] [Cross Ref]
4. Swami M. K., Manhas S., Buddhiwant P., Ghosh N., Uppal A., Gupta P. K., “Polar decomposition of 3 x 3 Mueller matrix: a tool for quantitative tissue polarimetry,” Opt. Express 14(20), 9324–9337 (2006).10.1364/OE.14.009324 [PubMed] [Cross Ref]
5. Wood M. F., Ghosh N., Moriyama E. H., Wilson B. C., Vitkin I. A., “Proof-of-principle demonstration of a Mueller matrix decomposition method for polarized light tissue characterization in vivo,” J. Biomed. Opt. 14(1), 014029 (2009).10.1117/1.3065545 [PubMed] [Cross Ref]
6. Backman V., Gurjar R., Badizadegan K., Itzkan I., Dasari R., Perelman L., Feld M., “Polarized light scattering spectroscopy for quantitative measurement of epithelial cellular structures in situ,” IEEE J. Sel. Top. Quantum Electron. 5(4), 1019–1026 (1999).10.1109/2944.796325 [Cross Ref]
7. Gurjar R. S., Backman V., Perelman L. T., Georgakoudi I., Badizadegan K., Itzkan I., Dasari R. R., Feld M. S., “Imaging human epithelial properties with polarized light-scattering spectroscopy,” Nat. Med. 7(11), 1245–1248 (2001).10.1038/nm1101-1245 [PubMed] [Cross Ref]
8. Bickel W. S., Bailey W. M., “Stokes vectors, Mueller matrices, and polarized scattered light,” Am. J. Phys. 53(5), 468–478 (1985).10.1119/1.14202 [Cross Ref]
9. Anderson R., “Measurement of Mueller matrices,” Appl. Opt. 31(1), 11–13 (1992).10.1364/AO.31.000011 [PubMed] [Cross Ref]
10. T. Wood, S. Thiemjarus, K. Koh, D. Elson, and G. Yang, “Optimal feature selection applied to multispectral fluorescence imaging,” Medical Image Computing and Computer-Assisted Intervention MICCAI 2008 2, 222–229 (2008). [PubMed]
11. Lu S., Chipman R., “Interpretation of Mueller matrices based on polar decomposition,” J. Opt. Soc. Am. A 13(5), 1106–1113 (1996).10.1364/JOSAA.13.001106 [Cross Ref]
12. Oberemok E. A., Savenkov S. N., “Structure of deterministic Mueller matrices and their reconstruction in the method of three input polarizations,” J. Appl. Spectrosc. 70(2), 224–229 (2003).10.1023/A:1023813117458 [Cross Ref]
13. Morio J., Goudail F., “Influence of the order of diattenuator, retarder, and polarizer in polar decomposition of Mueller matrices,” Opt. Lett. 29(19), 2234–2236 (2004).10.1364/OL.29.002234 [PubMed] [Cross Ref]
14. Shukla P., Pradhan A., “Mueller decomposition images for cervical tissue: potential for discriminating normal and dysplastic states,” Opt. Express 17(3), 1600–1609 (2009).10.1364/OE.17.001600 [PubMed] [Cross Ref]
15. Ghosh N., Wood M. F., Li S. H., Weisel R. D., Wilson B. C., Li R. K., Vitkin I. A., “Mueller matrix decomposition for polarized light assessment of biological tissues,” J Biophotonics 2(3), 145–156 (2009).10.1002/jbio.200810040 [PubMed] [Cross Ref]
16. T. C. Wood, and D. S. Elson, “Polarization characterisation of laparoscope systems for polarization resolved tissue imaging,” in “Biomedical Optics,” (Optical Society of America, 2010), p. BTuD29.
17. MacGregor A., “Method for computing homogeneous liquid-crystal conoscopic figures,” J. Opt. Soc. Am. B 7(3), 337–347 (1990).10.1364/JOSAA.7.000337 [Cross Ref]

Articles from Biomedical Optics Express are provided here courtesy of Optical Society of America