1 May 2010 Rotating linear polarization imaging technique for anisotropic tissues
Author Affiliations +
Abstract
A novel rotating linear polarization imaging technique is developed to characterize the anisotropic properties of tissues. Differences of orthogonal linear polarization with different incident and detection polarization angles are fitted to an analytical function to retrieve a set of parameters. Experiments with different tissues and Monte Carlo simulations indicate that two of the parameters, G and φ3/2, are correlated to the anisotropic property and the orientation angle of the fibrous structure in the media. The technique can be used for clinical diagnosis.

1.

Introduction

There has been growing interest in studies of polarization imaging in turbid media for achieving higher imaging quality or better characterization of the samples.1, 2, 3, 4 The primary advantage for clinical application of polarization imaging is that it suppresses backscattered photons from deep tissues and improves the image contrast and resolution in the superficial layers where many cancers start to develop. Degree of polarization (DOP) imaging with linearly polarized light has been successfully used in clinical diagnosis of many cancers.2 Other polarization measurements, such as Mueller matrix polarimetry5, 6 and Stokes polarimetry,7 have also been used to retrieve structural and optical information of turbid media.

Most biomedical tissues are anisotropic, and many contain prominent fibrous structures.8 It has been found in our work and by other authors that linear polarization images with single incident polarization can be very different when the incident polarization or sample orientation is changed.7 It could be uncertain whether the observed differences represent physiological differences or the optical artifacts due to changes in incident polarization or sample orientations for different tests.9, 10, 11 Such uncertainty may affect objective diagnostics in clinical applications. Objective and quantitative characterization of anisotropic tissues remains a challenge.

In this work, we present a novel technique, rotating linear polarization imaging (RLPI), which provides a set of parameters to objectively and quantitatively characterize the anisotropic properties of the scattering media. It is based on a reflective linear polarization imaging system but with independently rotating incident and detection linear polarization angles. Linear differential polarization (LDP) images are fitted at each pixel to an analytical function of the incident and detection polarization angles to retrieve a set of new images of the fitting parameters. Both experiments on different tissues and Monte Carlo simulations show that at least two of the parameters have clear physics meanings, representing the structure anisotropy and the orientation angle of the fibrous structure in the samples. The images are all independent of the incident and detection polarization angles.

2.

Experiment Setup and Data Processing

The experimental setup (Fig. 1 ) is similar to a typical reflective linear polarization imaging system.2 650-nm light from a 1-W LED is collimated by a lens L1, propagates though a linear polarizer P1, and illuminates the sample at 25deg to the normal of sample surface, which is defined as the Z axis. Backscattered photons from the sample pass through the analyzer P2, are collected by lens L2, and recorded by a 12-bit charge-coupled device (CCD) (Retiga EXi, QImaging, Surrey, British Columbia, Canada). The image size is 520×696pixels and each pixel is 13×13μm2 . Driven by two-step motors, both P1 and P2 can rotate around their optical axis to vary the polarization angles for illumination θi and for detection θs . The samples are different types of tissues, such as muscles, connective tissues, and fats from different animals, including bovine, chicken, and porcine. A glass plate is placed in close contact to the sample to direct surface glare away from the detector.

Fig. 1

(a) Schematics of experimental setup and (b) reference frames of experiment. θi and θs represent the incidence and detection polarization angles.

036014_1_035003jbo1.jpg

In each experiment, we record a series of images I(θi,θs) . For each incident polarization angle θi , we take a number of image pairs corresponding to orthogonal detection polarization angles θs and θs+π2 . However, at least once we choose θs=θi to obtain the DOP image. Then we calculate LDP at each pixel (x,y) :

1

LDP(θi,θs)=I(θi,θs)I(θi,θs+π2).

For a fixed θs , the LDPx,y(θi) at each pixel (x,y) is cosinoidal [Fig. 2 ]:

2

LDPx,y(Iin2)=h1cos(2θih2)+h3,
where Iin is the incident intensity.

Fig. 2

Cosinoidal relationship of LDP (a) versus θi when θs=0 , and (b) h12Iin24 and (c) h3Iin2 versus θs . Three curves correspond to three different pixels on the same image series of bovine skeletal muscle. Solid curves are the fitting results of the experimental data to Eqs. 2, 3, 5.

036014_1_035003jbo2.jpg

It was found in the experiments that Eq. 2 stays valid for all the different samples we have tested. In fact, a similar cosinoidal expression was also proposed for polarization measurements in transmission mode.12 The authors speculated a connection between the parameters in the expression and tissue characteristics, but did not give experimental evidences or detailed analysis.

Further examinations reveal that h12 and h3 are also cosinoidal functions of the detection polarization angle θs [Figs. 2 and 2]:

3

h12=Acos(4θsφ1)+B,

4

h2=φ2(θs),

5

h3=Ccos(2θsφ3).
Experimental data do not indicate a simple relation between h2 and θs .

Combining Eqs. 2, 3, 4, 5, we obtain an analytical expression of LDP:

6

LDP(θi,θs)=12Iin[Acos(4θsφ1)+B]12cos[2θiφ2(θs)]+12IinCcos[2θsφ3].
For the LDP image series, fitting the data at each pixel to Eq. 6 using nonlinear least squares method results in a set of new parameters A , B , C , φ1 , and φ3 . These parameters are independent of incident and detection polarization angles and are related to different structural or optical properties of the sample. The new parameters form their own images as shown in Fig. 3 . Figure 3 corresponds to a new parameter defined as:

7

GAB.
More discussions about these parameters are given in following sections.

Fig. 3

Different images of the chicken heart cross section. (a): intensity image, and (b) and (c) are DOP images corresponding to incident polarization angle θi and θi+π4 . (d) through (h) are images of A , B , C , G , and φ32 . (a) through (g) are normalized to the maximum value of the image. The color bar for (h) is π2 to π2 . Background is set to zero in (a) through (h). Mean values of φ32 are calculated within the equally spaced black squares along the arc in (h) and are shown in (i).

036014_1_035003jbo3.jpg

3.

Interpretation of the Fitted Parameters

In this section, we examine in detail the parameters obtained in RLPI to understand their physical origins, particularly their relations with the structure of the biological samples.

The sample is a chicken heart (Fig. 4 ) whose lower part is removed to reveal a 9-mm -diam cross sectional area of the left ventricle. Typically, 36 images are captured for each experiment, corresponding to six incident and detection polarization angles, respectively. Each experiment takes about 8s limited only by the rotation mechanics of the polarizers.

Fig. 4

(a) Picture and (b) schematic anatomical diagram of a chicken heart. The dashed line indicates the cut position.

036014_1_035003jbo4.jpg

Figure 3 shows a CCD image of the cross section. It is divided into three circular sections marked 1, 2, and 3. It is known from anatomy that the highly fibrous heart muscles are aligned to different directions at different sections of the heart. The muscle fibers are aligned perpendicular to the imaging plane at zone 1, and parallel to the imaging plane at zone 2. Close to the center at zone 3, alignment of the heart muscles varies.

In the RLPI experiment, at each incident polarization angle θi we record a series of image pairs corresponding to detection polarization angles θs and θs+π2 . θs can take any value, but has to include the special case θs=θi . Therefore, during the RLPI experiment, we can obtain the DOP and a series of corresponding LDP different incident polarization states. Figures 3 and 3 show two DOP images of the same sample but with different incident polarization angles θi and θi+π4 . It is clearly demonstrated that DOP images with different incident polarization can be very different.

Figures 3, 3, 3, 3, 3 show the images of A , B , C , G , and φ32 . Figure 3 shows clearly the three zones: G is small at zone 1 and high at zone 2. It shows more complex structures at the central area, zone 3. Considering the anatomical structure of chicken heart, one may suspect that G is higher when the fibrous structure of the sample is well aligned in the image plane. Figure 3 shows that in zone 2 where G stays high, φ32 varies continuously along the circle. We choose a series of equally spaced 7×7pixel square areas ( 0.1×0.1mm2 at the sample) along the circle, as shown in Fig. 3, and calculate the mean value of φ32 . Figure 3 shows that φ32 varies almost linearly with the angle along the arc and may represent the orientation of a fibrous structure.

RLPI is independent on the incident polarization angles. Since the experimental setup is not axially symmetric, the fitted parameters may vary slightly as the sample rotates. We examine this effect by taking RLPI measurements as the sample rotates around the Z axis. For quantitative examination, we choose a 7×7pixel rectangle region on the images ( 0.1×0.1mm2 at the sample) as the region of interest (ROI) and calculate the mean values of DOP, A , B , C , φ32 , and G for different sample orientations. The results are shown in Fig. 5 .

Fig. 5

Dependence of fitting parameters to sample orientation. The sample rotates around the Z axis. In (a), curves for A and B are both normalized to the maximum of B , and C is normalized to its own maximum.

036014_1_035003jbo5.jpg

Figure 5 shows that DOP varies with the sample orientation as a sinusoidal curve with π2 period. A and B also show periodical but nonsinusoidal variation with a period of π . The magnitude of the variation is comparable with DOP. The variation of C with sample rotation is more prominent than A and B , and the period varies at different parts of the image. G stays almost constant during sample rotation, which is consistent with the assumption that G is related to the order of alignment of the fibrous structure. φ32 varies almost linearly with the sample’s rotation angle, proving that it indeed represents the orientation angle. It is also expected that the φ32 curve should be subject to “smoothing” as shown in Fig. 5, since the orientation of the fibrous structure may not be the same within the ROI.

More evidences on the relation between G and anisotropy of the sample are provided by RLPI experiments on tissues of different fibrous structure and alignment. Mean values and relative standard deviations (RSD=SDmean*100%) of G are calculated over 7×7pixel ROI and listed in Table 1 .

Table 1

G for different samples.

GroupSample G (mean±RSD)
1Bovine skeletal muscle withregular fiber alignmentparallel to sample surface 0.91±8.8%
2The same as group 1, but thefibers are aligned perpendicularto sample surface 0.19±68%
3Porcine skeletal muscle withregular fiber alignmentto sample surface 0.95±7.4%
4Porcine skeletal muscle withirregular fiber alignment tosample surface 0.63±33%
5Gallinaceous skeletal musclewith regular fiber alignment tosample surface 0.96±4.2%
6Gallinaceous skeletal musclewith irregular fiber alignment tosample surface 0.65±38%
7Porcine fatty tissue 0.17±53%
8Porcine liver 0.07±57%

It is clearly seen from Table 1 that there is a strong correlation between G and the order of alignment, or anisotropy, of the fibrous structure in the imaging plane ( x-y plane). All the different skeletal muscles in groups 1, 3, and 5 contain prominent and well-aligned fibrous structures in the imaging plane. They all show very high G and small RSD. Samples in group 2 are the same bovine skeletal muscle as that in group 1, but are rotated to make the well-aligned muscular fibers run perpendicular to the imaging plane. This results in a sharp reduction of anisotropy in the x-y plane. Therefore, G is much smaller and not well defined as the RSD increases sharply. Samples in groups 4 and 6 are porcine and gallinaceous skeletal muscles, the same as those in group 3 and 5, but their fibrous structures look less regular. The structure is less ordered and the anisotropy should be smaller. Table 1 does show significant reduction in the measured G and an increase in the RSD. Porcine liver and fat tissues in groups 7 and 8 are known to be isotropic.13 The measured Gs are small and the RSDs are large.

Therefore, we can conclude that G and φ32 correlate quantitatively to the anisotropy of the fibrous structure and the fiber orientation angle in the sample. Images of these two parameters provide more objective information on the microstructure and optical property of tissues.

4.

Discussion

Although the analytical expression for LDP, Eq. 6, was first obtained empirically, it can also be derived from the text book scattering theory for polarized photons.14 For simplicity, we consider normal incidence and detection, which do not affect the validity of the conclusion.

Stokes vector of the detected photons Sout can be expressed as the products of the Mueller matrixes of the analyzer Ms , sample M , polarizer Mi , and the Stokes vector of the incident light Sin :

8

Sout=MsMMiSin,
where,

9

Ms=12(1cos2θssin2θs0cos2θscos22θssin2θscos2θs0sin2θssin2θscos2θssin22θs00000),

10

M=(M11M12M13M14M21M22M23M24M31M32M33M34M41M42M43M44),

11

Mi=12(1cos2θisin2θi0cos2θicos22θisin2θicos2θi0sin2θisin2θicos2θisin22θi00000),

12

Sout=[Iout,Qout,Uout,Vout]T,

13

Sin=[Iin,Qin,Uin,Vin]T.
If the incident light is nonpolarized, Qin=Uin=Vin=0 , we get

14

LDP(θi,θs)=Iout(θi,θs)Iout(θi,θs+π2)=12Iin[Acos(4θsφ1)+B]12cos[2θiφ2(θs)]+12IinCcos(2θsφ3),
where,

15

B=(M222+M232+M322+M332)2;

16

A=[B2(M22M33M23M32)2]12;

17

C=M212+M312;

18

tanφ1=2M22M32+2M23M33M222+M232(M322+M332);

19

tanφ2=M23cos(2θs)+M33sin(2θs)M22cos(2θs)+M32sin(2θs);

20

tanφ3=M31M21.
From Eqs 7, 15, 16, 17, 20, we find that: A , B , and G relate to the four central elements of the Mueller matrix, which couple the Stokes components Q and U . Since both Q and U are defined as polarization differences, A , B , and G should be sensitive, as in DOP, only to the structure and optical properties of the superficial layer. For φ32 and C , they are linked to Mueller matrix elements in the first column, which couple I to Q and U . Since I contains contributions from all the backscattered photons, φ32 and C should be sensitive to both superficial and deep layer properties of the sample.

More evidence on the physics origin of the RLPI parameters is also given by Monte Carlo simulations. In a sphere-cylinder scattering model, general anisotropic tissues are approximated to a mixture of solid spherical and infinitely long cylindrical scatterers suspended in liquid.15, 16 Orientation angle of the cylindrical scatterers fluctuates around a fixed direction following a Gaussian probability distribution. The fraction of spherical scatterers and the width of the Gaussian distribution vary independently to change the anisotropy of the scattering media. Using a Monte Carlo program,15 we simulate the propagation of polarized photons in such samples and obtain the RLPI parameters following the same experimental condition and data analysis procedures as described in previous sections. For the results shown in Fig. 6 , the diameters of the spherical and cylindrical scatterers are 1.5 and 2.0μm , respectively, refractive indices of the scatterers and the liquid are 1.4 and 1.35, and the total scattering coefficient is 100cm1 .

Fig. 6

(a) Dependence of G on the fraction of spherical scatterers and the width of the Gaussian distribution. (b) Dependence of φ32 on orientation angle of the cylindrical scatterers.

036014_1_035003jbo6.jpg

As the fraction of spherical scatterers in the mixture or the width of the Gaussian distribution increases, the sample becomes less anisotropic. As shown in Fig. 6, G derived from the simulated data also decreases monotonously, which confirms that G is related to the structural anisotropy of the sample.

Then we examine the relation between φ32 and the cylinder orientation. The cylindrical scatterers take 80% of the total number of scatterers and their orientation fluctuation within a Gaussian distribution of 25-deg full-width at half maxium (FWHM). As shown in Fig. 6, φ32 derived from the simulated data varies linearly with the most probable orientation of the cylinders, which confirm that φ32 represents the fiber orientation in the sample.

Monte Carlo simulation is also used to study the dependence of the RLPI parameters on the backscattered photons from different depths. Figure 7 shows the relation between A , B , C , and the thickness of the sample. As the thickness of the sample increases, A and B are quickly saturated, but C keeps increasing. It is evident that backscattered photons from the superficial layer make the most significant contribution to A and B , while photons from the much deeper layer still contribute to C . Therefore, both A and B represent the information from the superficial layer, and C carries the information from deeper layers.

Fig. 7

Dependences of A , B , and C on the sample’s thickness. Curves for A and B are both normalized to the maximum of B , and C is normalized to its own maximum.

036014_1_035003jbo7.jpg

5.

Conclusion

We develop a rotating linear polarization imaging (RLPI) technique. It is shown that linear polarization difference (LDP) follows an analytical relation to incident and detection polarization angles. By fitting LDP images pixel by pixel to an analytical expression, we obtain a set of new images of the fitting parameters that are independent on incident polarization. Both experiments on tissues and Monte Carlo simulations prove that two of the parameters, G and φ32 , correlate to the anisotropy and orientation angle of the fibrous structure in tissues. The new technique can be used for objective and quantitative characterization of anisotropic tissues, and for clinical diagnosis of skin diseases.

Acknowledgment

This work has been supported by the National Natural Science Foundation of China (grants 60778044 and 10974114), the Ministry of Science and Technology (grant 2006CB70570), and National High Technology Research and Development Program of China (863 program number 2006AA06Z402).

References

1.  S. G. Demos and R. R. Alfano, “Optical polarization imaging,” Appl. Opt.0003-6935 36(1), 150–155 (1997). 10.1364/AO.36.000150 Google Scholar

2.  S. L. Jacques, J. R. Roman, and K. Lee, “Imaging skin pathology with polarized light,” J. Biomed. Opt.1083-3668 7(3), 329–340 (2002). 10.1117/1.1484498 Google Scholar

3.  S. P. Morgan and I. M. Stockford, “Surface-reflection elimination in polarization imaging of superficial tissue,” Opt. Lett.0146-9592 28(2), 114–116 (2003). 10.1364/OL.28.000114 Google Scholar

4.  H. R. Shao, Y. H. He, W. Li, and H. Ma, “Polarization-degree imaging contrast in turbid media: a quantitative study,” Appl. Opt.0003-6935 45(15), 4491–4496 (2006). 10.1364/AO.45.004491 Google Scholar

5.  B. D. Cameron, M. J. Rakovic, M. Mehrubeoglu, G. W. Kattawar, S. Rastegar, L. V. Wang, and G. L. Coté, “Measurement and calculation of the two-dimensional backscattering Mueller matrix of a turbid medium,” Opt. Lett.0146-9592 23(7), 485–487 (1998); 10.1364/OL.23.000485 Google Scholar B. D. Cameron, M. J. Rakovic, M. Mehrubeoglu, G. W. Kattawar, S. Rastegar, L. V. Wang, and G. L. Coté, “Measurement and calculation of the two-dimensional backscattering Mueller matrix of a turbid medium: errata,” Opt. Lett.0146-9592 23(20), 1630–487 (1998). 10.1364/OL.23.001630 Google Scholar

6.  S. Manhas, M. K. Swami, P. Buddhiwant, N. Ghosh, P. K. Gupta, and K. Singh, “Mueller matrix approach for determination of optical rotation in chiral turbid media in backscattering geometry,” Opt. Express1094-4087 14(1), 190–202 (2006); 10.1364/OPEX.14.000190 Google Scholarsee  http://www.opticsinfobase.org/oe/abstract.cfm?URI=OPEX-14-1-190Google Scholar

7.  P. J. Wu and J. T. Walsh, “Stokes polarimetry imaging of rat-tail tissue in a turbid medium: Degree of linear polarization image-maps using incident linearly polarized light,” J. Biomed. Opt.1083-3668 11(1), 014031 (2006). 10.1117/1.2162851 Google Scholar

8.  T. Binzoni, C. Courvoisier, R. Giust, G. Tribillon, T. Gharbi, J. C. Hebden, T. S. Leung, J. Roux, and D. T. Delpy, “Anisotropic photon migration in human skeletal muscle,” Phys. Med. Biol.0031-9155 51(5), N79–N90 (2006). 10.1088/0031-9155/51/5/N01 Google Scholar

9.  W. Cheong, S. A. Prahl, and A. J. Welch, “A review of the optical properties of biological tissues,” IEEE J. Quantum Electron.0018-9197 26(12), 2166–2185 (1990). 10.1109/3.64354 Google Scholar

10.  A. Kienle, F. K. Forster, R. Diebolder, and R. Hibst, “Light propagation in dentin: influence of microstructure on anisotropy,” Phys. Med. Biol.0031-9155 48(2), N7–N14 (2003). 10.1088/0031-9155/48/2/401 Google Scholar

11.  A. Kienle and R. Hibst, “Light guiding in biological tissue due to scattering,” Phys. Rev. Lett.0031-9007 97(1), 018104 (2006). 10.1103/PhysRevLett.97.018104 Google Scholar

12.  S. K. Gayen, M. E. Zevallos, M. Alrubaiee, J. M. Evans, and R. R. Alfano, “Two-dimensional near-infrared transillumination imaging of biomedical media with a chromium-doped forsterite laser,” Appl. Opt.0003-6935 37(22), 5327–5336 (1998). 10.1364/AO.37.005327 Google Scholar

13.  S. L. Jacques, J. R. Roman, and K. Lee, “Imaging superficial tissues with polarized light,” Lasers Surg. Med.0196-8092 26(2), 119–129 (2000). 10.1002/(SICI)1096-9101(2000)26:2<119::AID-LSM3>3.0.CO;2-Y Google Scholar

14.  C. F. Bohden and D. Huffman, “Electromagnetic theory,” Chap. 2 in Absorption and Scattering of Light by Small Particles, pp 44–56, Wiley, New York (1983). Google Scholar

15.  T. L. Yun, N. Zeng, W. Li, D. Z. Li, X. Y. Jiang, and H. Ma, “Monte Carlo simulation of polarized photon scattering in anisotropic media,” Opt. Express1094-4087 17(19), 16590–16602 (2009); 10.1364/OE.17.016590 Google Scholarsee  http://www.opticsinfobase.org/oe/abstract.cfm?uri=oe-17-19-16590Google Scholar

16.  N. Zeng, X. Y. Jiang, Q. Gao, Y. H. He, and H. Ma, “Linear polarization difference imaging and its potential applications,” Appl. Opt.0003-6935 48(35), 6734–6739 (2009). 10.1364/AO.48.006734 Google Scholar

© (2010) Society of Photo-Optical Instrumentation Engineers (SPIE)
Ran Liao, Ran Liao, Nan Zeng, Nan Zeng, Xiaoyu Jiang, Xiaoyu Jiang, Dongzhi Li, Dongzhi Li, Tianliang Yun, Tianliang Yun, Yonghong He, Yonghong He, Hui Ma, Hui Ma, } "Rotating linear polarization imaging technique for anisotropic tissues," Journal of Biomedical Optics 15(3), 036014 (1 May 2010). https://doi.org/10.1117/1.3442730 . Submission:
JOURNAL ARTICLE
6 PAGES


SHARE
Back to Top