Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Integrated Mueller-matrix near-infrared imaging and point-wise spectroscopy improves colonic cancer detection

Open Access Open Access

Abstract

We report the development and implementation of a unique integrated Mueller-matrix (MM) near-infrared (NIR) imaging and Mueller-matrix point-wise diffuse reflectance (DR) spectroscopy technique for improving colonic cancer detection and diagnosis. Point-wise MM DR spectra can be acquired from any suspicious tissue areas indicated by MM imaging. A total of 30 paired colonic tissue specimens (normal vs. cancer) were measured using the integrated MM imaging and point-wise MM DR spectroscopy system. Polar decomposition algorithms are employed on the acquired images and spectra to derive three polarization metrics including depolarization, diattentuation and retardance for colonic tissue characterization. The decomposition results show that tissue depolarization and retardance are significantly decreased (p<0.001, paired 2-sided Student’s t-test, n = 30); while the tissue diattentuation is significantly increased (p<0.001, paired 2-sided Student’s t-test, n = 30) associated with colonic cancer. Further partial least squares discriminant analysis (PLS-DA) and leave-one tissue site-out, cross validation (LOSCV) show that the combination of the three polarization metrics provide the best diagnostic accuracy of 95.0% (sensitivity: 93.3%, and specificity: 96.7%) compared to either of the three polarization metrics (sensitivities of 93.3%, 83.3%, and 80.0%; and specificities of 90.0%, 96.7%, and 80.0%, respectively, for the depolarization, diattentuation and retardance metrics) for colonic cancer detection. This work suggests that the integrated MM NIR imaging and point-wise MM NIR diffuse reflectance spectroscopy has the potential to improve the early detection and diagnosis of malignant lesions in the colon.

© 2016 Optical Society of America

1. Introduction

Colorectal cancer (CRC) is the third most commonly diagnosed malignancy in males and the second in females worldwide, with an estimated 1.4 million cases and 693,900 deaths occurring in 2012 [1]. In Singapore, CRC has become the most frequent cancer with a total of 9,324 new cases diagnosed from 2010 to 2014 [2]. Current routine screening of CRC uses white light reflectance (WLR) colonoscopy which may reduce CRC incidence and mortality [3]. However, some individuals are still diagnosed with CRC despite recent colonoscopy [3]. This is probably because conventional WLR colonoscopy heavily relies on the visualization of gross mucosal features associated with neoplastic transformation [4]. Subtle tissue changes may not be apparent, limiting its diagnostic accuracy. Consequently, existing diagnostic guidelines recommend extensive but random biopsy samplings during colonoscopic inspections of patients [5], followed by the microscopic examination which is highly subjective and depends heavily on the experiences of the pathologists. Overall, the current approach for colonic tissue diagnosis is clinically labor intensive and a burden to the patients. There is a need to develop advanced optical diagnostic techniques for objective diagnosis and characterization of colonic tissue.

Polarized light imaging/spectroscopy has been comprehensively investigated for tissue diagnosis [6–17]. Polarized light implementation offers several compelling advantages as follows: (i) surface and beneath-the-surface detection of biological tissue taken from the tissue depolarization properties [11, 15]; (ii) tissue anisotropy analysis through the tissue diattentuation and retardance [13, 14]; (3) enhanced tissue diagnosis through the combination of complementary depolarization, diattentuation and retardance of the tissues [9]. Among the various polarized light imaging/spectroscopy techniques developed [6–14, 18], Mueller-Matrix polarimetry is capable of measuring the complete polarimetric transfer function [6–9], known as Mueller-matrix, of the bulk biological tissues which are optically inhomogeneous, birefringent, and absorbing media [19]. Currently, biomedical Mueller-Matrix polarimetry is mostly centered on the use of short visible wavelengths of illumination light that has a limited penetration depth and cannot detect lesions in deeper areas [8,9]. The near-infrared (NIR) light, on the other hand, penetrates much deeper into the tissue, and is well suited for deep tissue diagnosis [12, 20–22]. Further, the reported Mueller-Matrix polarimetries only acquired either the images [7,8] or the optical spectra [9] of the biological tissues alone. In this work, we report on the development of a unique integrated Mueller-Matrix NIR imaging and point-wise Mueller-Matrix spectroscopy system for colonic tissue diagnosis and characterization. Point-wise Mueller-Matrix spectra can be acquired from any suspicious tissue areas indicated by Mueller-Matrix imaging. Polar decomposition algorithms are employed on the acquired Mueller-Matrix images/spectra to derive three polarization metrics including depolarization, diattentuation and retardance. Partial least squares discriminant analysis (PLS-DA) and leave-one tissue site-out, cross-validation (LOSCV) are implemented on the derived spectroscopic polarization metrics (i.e., depolarization, diattentuation and retardance) to develop robust spectral diagnostic models for the differentiation between cancerous and normal colonic tissues.

2. Materials and methods

2.1 Integrated Mueller-matrix NIR imaging and point-wise diffuse reflectance spectroscopy system

Figure 1 shows the schematic of the integrated Mueller-Matrix NIR imaging and point-wise diffuse reflectance (DR) spectroscopy system developed for tissue measurements. The light from a tungsten halogen lamp (HL-2000, Ocean Optics Inc., Dunedin, FL) is coupled into a 200 fiber and passes through a beam expander, a long-pass filter (FEL0700, Thorlabs, Newton, NJ), a polarizer (LPNIR100-MP2, Thorlabs, Newton, NJ), and a quarter waveplate (AQWP10M-980, Thorlabs, Newton, NJ) for tissue illumination. The NIR diffuse reflectance photons backscattered from the tissue pass through a quarter waveplate (AQWP10M-980, Thorlabs, Newton, NJ), a polarizer (LPNIR100-MP2, Thorlabs, Newton, NJ), a collection lens, and a specially designed point-spectrum optical adaptor [23] before detected by a CCD camera (Pixis 1024, Princeton Instruments, Trenton, NJ). The customized optical adaptor comprises three lenses (f = 50 mm), a thin quartz glass plate (25 × 25 × 1 mm3) coated with a gold mirror (diameter of 100 µm, reflection of ~99% in 700-1100 nm) and a 2-D motorized translational stage. During tissue measurement, a small portion of the backscattered light was reflected by the optical adaptor and collected by a spectrometer (QE65000, Ocean Optics Inc., Dunedin, FL) for tissue spectroscopic analysis [23].

 figure: Fig. 1

Fig. 1 Schematic of the integrated Mueller-Matrix NIR imaging and point-wise Mueller-Matrix diffuse reflectance spectroscopy system developed for colonic tissue measurements: collimator (C); long pass filter (LP); lens (L); polarizer (P); quarter waveplate (QWP).

Download Full Size | PDF

To acquire the 4 by 4 Mueller-Matrix DR images/spectra, the fast axis of the polarizers (P1, P2) is fixed while the quarter waveplates (QWP1, QWP2) were rotated automatically by using two pairs of gears and two sets of step motors in the excitation and collection paths, respectively. The rotation speed ratio of the two quarter waveplates is fixed at 1:5. The detected intensity is Fourier modulated as follows [24, 25]:

I=a0+n=112(ancosnwt+bnsinnwt)
where ω is the rotation speed of QWP1 (Fig. 1), t is the exposure time of the camera, and a0, an, bn are the Fourier coefficients which can be derived through the detected intensity I. With the integrated Mueller-Matrix NIR imaging and point-wise spectroscopy system developed, a set of 25 Mueller-Matrix images/spectra can be acquired from colonic tissues in tandem within 5 s; hence the 4 by 4 Mueller-Matrix imaging/point-wise spectroscopy can be generated for tissue diagnosis and characterization. Further automatic motorization of the small gold mirror coated on the quartz plate together with the point-wise spectral measurement module enables a rapid movement of the dark spot (of 0.2 mm in diameter due to the reflection of gold mirror in the point spectrum optical adapter) to any spot of the imaged tissue of interest on the Mueller-Matrix image, and the subsequent 4 by 4 Mueller-Matrix point-wise spectroscopy can be realized within 1 s. One notes that the Eq. (1) applies only when the retarders used are quarter waveplates, and their angular speed ratio is fixed at 1:5. More general formalism for Mueller-Matrix measurements with dual rotating retarders rotating with different speed ratios can be found in [26].

To derive the colonic tissue polarization metrics (i.e., diattentuation D, depolarization Δ, and retardance R), polar decomposition [27] is implemented on the 4 × 4 Mueller-Matrix images/spectra acquired with the system developed (Fig. 1). Briefly, the tissue Mueller-Matrix M is expressed as the product of three 4 by 4 matrices: the diattentuation matrix (MD), the depolarization matrix (MΔ), and the retardance matrix (MR) [27]:

M=[m11m12m13m14m21m22m23m24m31m32m33m34m41m42m43m44]
M=MΔMRMD

The diattentuation vector D¯ is determined by the elements on the first row of the Mueller-Matrix [27]:

D¯=1m11(m12,m13,m14)

The diattentuation, D, and the diattentuation matrix, MD, can thus be determined as [27]:

D=1m11m122+m132+m142
MD=m11[1DTDmD]
where mD=1D2I+(11D2)D^D^T, and I is the 3 by 3 identity matrix, D^=D¯/|D^| denotes the unit vector along D¯.

From Eqs. (3) and (6), a Mueller-Matrix can be defined based on M [27]:

M'=MMD1

Note that M' has no diattentuation, and it can be further decomposed as a retardance matrix (MR) followed by the depolarization matrix (MΔ) [27]:

MΔMR=MMD1=[10TPΔmΔ][10T0mR]=[10TPΔmΔmR]=[10TPΔmΔmR]=[10TPΔm']=M'
Here, m' is a 3 by 3 submatrix of M'. Letλ1,λ2,λ3 be the eigenvalues of (m')(m')T. The submatrix and of MΔ can be determined:

mΔ=±[m'(m')T+(λ1λ2+λ2λ3+λ3λ1)I]1×[(λ1+λ2+λ3)m'(m')T+λ1λ2λ3I]

If the determinant of m' is negative, the minus sign is applied. Otherwise, the plus sign is applied [27].

Using Eqs. (8) and (9), MΔ can be determined, and MR can be obtained by:

MR=MΔ1M'

Finally, the depolarization Δ, and retardance R can be calculated as follows:

Δ=1|trace(MΔ1)|3
R=cos1[trace(MR)21]

To validate the performances of the system developed, the Mueller-Matrix NIR spectra of a half waveplate and a quarter waveplate were measured and decomposed. The difference between the measured retardance and that provided by the manufacturer is less than 3%, confirming the robustness of the system developed.

2.2. Statistical analysis

The unpaired two-sided Student’s t-test was used to evaluate the decomposed Mueller-Matrix spectroscopic differences between cancer and normal colonic tissues [28]. Partial least squares (PLS) - discriminant analysis (DA) was applied on the derived spectroscopic polarization metrics for developing spectral diagnosis models [28]. Leave-one-tissue site out, cross-validation (LOSCV) was further used to assess and optimize the PLS-DA model complexity, while reducing the risk of over-fitting. The above multivariate statistical analysis was performed using in-house written scripts in the Matlab programming environment (Mathworks. Inc., Natick, MA).

2.3. Colonic tissue specimens

A total of 30 paired (i.e., normal vs cancer) colonic tissue specimens (average size of ~6 x 3 x 3 mm3) were collected from 30 patients (18 men and 12 women with a mean age of 56) who underwent partial colectomy or surgical resections with clinically suspicious lesions or histopathologically proven malignancies in the colon. All patients preoperatively signed an informed consent permitting the investigative use of the tissue, and this study was approved by the Institutional Review Board (IRB) of the National Healthcare Group (NHG) of Singapore. Immediately after surgical resections, tissue specimens were put into vials with physiological saline solution which were stored in a flask with ice (−4 °C) for sending to the Optical Bioimaging Laboratory within 10 minutes for Mueller-Matrix (MM) NIR imaging and point-wise spectroscopy measurements. The paired tissue specimens from each patient were placed on a quartz glass slide (26 × 76 × 1.2 mm3) (cancer tissue was placed at the lower part of the slide while the normal one was placed at upper part of the slide) for MM NIR imaging measurements. After the MM NIR imaging/spectra acquisitions, the tissue specimens were fixed in 10% formalin solution and then sent back to the hospital for histopathological examinations. The histopathological examinations confirmed that 30 tissue specimens were normal, and 30 tissue specimens were cancer (moderately differentiated adenocarcinoma). Figure 2 shows the NIR diffuse reflectance image of one typical paired colonic tissues (~6 x 3 x 3 mm3).

 figure: Fig. 2

Fig. 2 Example of an NIR diffuse reflectance image of a paired colon tissues (~6 x 3 x 3 mm3). Normal tissue is on the upper part, while cancer tissue is in the lower part of the image. One notes that the dark spot (~0.2 mm in diameter) is due to the refection of gold mirror coated in the point spectrum optical adapter (Fig. 1).

Download Full Size | PDF

3. Results and discussion

With the integrated Mueller-Matrix NIR imaging and point-wise spectroscopy technique developed, we are able to acquire 4 by 4 Mueller-Matrix NIR images of 30 paired colonic tissues. Figure 3 shows the representative normalized Muller matrix images of the paired (normal vs. cancer) colonic tissue sample as confirmed by histological examinations. All the Mueller-Matrix elements (except m11) are normalized by m11. It is observed from Fig. 3 that the diagonal elements of the Mueller-Matrix are much higher than the non-diagonal elements, reflecting a high depolarization power of the colon tissue. Interestingly, m22 and m33 element are identical not only for normal colon tissues, but also for the cancerous colon tissues. Besides, the magnitude of m22 and m33 are much higher than that of m44, indicating that the backscattered light is less depolarized when the incident light is linearly rather than circularly polarized. The magnitudes of m22 and m33 are higher for colon cancer than normal colon tissue, indicating a lower depolarization power of colon cancer tissues. The results we observed are consistent with previous reports in literature [9, 17]. We also found the existence of non-diagonal Mueller-Matrix images (i.e., m34, and m43), demonstrating the anisotropy of colonic tissues. The colonic anisotropy might originate from NIR light propagation in birefringent collagen located in the submucosa layer of the colonic tissue [29]. This is because the NIR light penetrates deeper (~1 mm) [12, 20–22] inside the colonic tissue, and the diffuse reflectance signal of the deeper collagen [29] can be detected and reflected in the non-diagonal Mueller-Matrix images.

 figure: Fig. 3

Fig. 3 Normalized Muller matrix NIR images of the paired (normal vs cancer) colonic tissues. All the Mueller-Matrix elements (except m11) are normalized to m11.

Download Full Size | PDF

Given the abundance of the diagnostic information contained in tissue Mueller-Matrix images acquired (Fig. 3), the quantitative biophysical polarization metrics (i.e., diattentuation, depolarization, and retardance) were derived using the polar decomposition algorithms [27] (Fig. 4). As shown in Fig. 4(a), the diattentuation of colonic cancer tissue is higher than that of normal colonic tissue. The diattentuation profile [Fig. 4(d)] confirms the significantly increased diattentuation for the colon cancer. One notes that biomolecules such as amino acids, proteins and nucleic acids exhibit diattentuation effects [30]. The higher magnitude observed for the diattentuation of cancerous tissue compared to normal tissue may be due to the enlarged nuclei and increased concentrations of chromatin (hence, nucleic acids) during colonic cancer development [31], which led to the increase in diattentuation effects in colonic cancer. Further, the decomposed depolarization image confirms that the cancerous colon clearly exhibits less depolarization effects [Figs. 4(b), 4(e)]. The decreased depolarization of cancer tissue can be attributed to the multiple scattering effects of polarized incident light in the bulk colonic tissue, originating from variations in the refractive indices of microstructures in cancer tissue [18]. Since an increase in cellular and nuclear sizes is accompanied with high cellular density and vascularization during cancer progression, an enhancement in anisotropic or Mie scattering (directionally dependent) of light in cancerous tissue could result in less depolarizing effects as compared to more isotropic or Rayleigh scattering in normal tissue [18, 31]. It is also observed that the decomposed retardance distributions are different between the normal and cancer colonic tissues. Overall, the very different diattentuation, depolarization and retardance images between normal and cancer colonic tissue observed (as shown in Fig. 4) demonstrate the potential of Mueller-matrix NIR imaging for early diagnosis and characterization of colonic cancer.

 figure: Fig. 4

Fig. 4 The processed (a) diattentuation image, (b) depolarization image, and (c) retardance image. Intensity profiles (d-f) along the black dotted lines as indicated on the decomposed Mueller-Matrix NIR images in (a-c).

Download Full Size | PDF

To determine the specific biochemical changes associated with colonic cancer, we have further acquired 60 sets (normal: n = 30; cancer: n = 30) of point-wise Mueller-matrix DR spectra from the suspicious tissue regions as indicated by the Mueller-Matrix images (Fig. 3). Figure 5 shows the typical 4 by 4 Mueller-Matrix DR spectra acquired from the histopathologically confirmed normal and cancerous colonic tissues. Obviously, the values of m22 and m33 are identical, which are consistent with Mueller-Matrix imaging (Fig. 3). Besides, the magnitude of spectroscopic m22 and m33 are higher in the colon cancer tissues, reconfirming the lower depolarizing power of cancer tissue. We also found non-diagonal spectroscopic Mueller-matrix elements (i.e., m34) to be non-zero, substantiating the anisotropy of colonic tissue. One notes that m11 generally represents the overall diffuse reflectance spectra of colonic tissue when unpolarized light is used [32]. Further analysis conducted on m11 (Fig. 5) casts light on the biochemical changes associated with colonic cancer. Overall, the NIR (700-1100 nm) DR spectra of both the normal and cancerous colonic tissues are dominated by the characteristic of water (970 nm) [33] and hemoglobin (940 nm) absorption bands [34]. Specifically, prominent water absorption valley is observed at 970 nm for both tissue types (Fig. 5). The water absorption band at 970 nm is due to the combination of the first harmonic of the O-H symmetric stretch vibration and the fundamental anti-symmetric stretch vibration from hydrogen bound O-H [33], making it a sensitive indicator of the local environment of water molecules [35]. Further, we found the water absorption valley is more obvious on the cancer tissue than that on the normal one, indicating increased water content for the cancerous colonic tissues. The enhanced metabolic rate in colonic cancer [36] contributes to the increased water content as water provides the conversion of mechanical energy developed by contractile proteins into the chemical energy useful in cell process [37]. The increased water for the cancerous colonic tissue has also been observed in other cancer tissues (e.g., esophagus [28], stomach [38], cervix [39–41] and brain [42]) by using Raman spectroscopy [4, 28, 38–42]. We have also found the decreased DR of hemoglobin band (940 nm) for the colonic cancer, signifying the increased hemoglobin content associated with colonic cancer tissue. The increased hemoglobin content for colonic cancer could be attributed to the increased microvasculatures in malignant tumors [36, 43]. The above observation is also consistent with previous NIR autofluorescence study in colonic cancer [44]. Further, the changes of tissue microstructural scattering properties (e.g. nucleus size, refractive index, microvasculature, etc) may also attribute to the significant differences of DR spectra between normal and colonic cancer tissues [45]. However, how the complex nature of the scattering process in pathologic tissue contributing to tissue DR spectroscopy still warrants further investigations.

 figure: Fig. 5

Fig. 5 Representative 4 by 4 Muller-Matrix NIR DR spectra recorded of normal and cancer colonic tissue. (Normal: __________ ; Cancer: _________).

Download Full Size | PDF

To develop robust multivariate spectral diagnostic algorithms for optical diagnosis of colonic cancer, we fully utilize all the diagnostic biochemical information contained in the Mueller-matrix DR spectra acquired (Fig. 5) by implementing PLS-DA and LOSCV technique on the quantitative Mueller-matrix metrics (i.e., diattentuation, depolarization and retardance) derived using polar decomposition algorithms [27]. Figure 6 shows the decomposed spectroscopic diattentuation ± 1 standard error (SE) (shaded area), depolarization ± 1 SE, and retardance ± 1 SE metrics of 30 paired colonic tissues. As consistent with the decomposed Mueller-Matrix images [Figs. 4(a), 4(b), 4(d), 4(e))], a significantly increased diattentuation (p<0.01) while a much reduced depolarization is observed in colonic cancer, demonstrating the potential of the Mueller-Matrix spectroscopy for colon cancer diagnosis. Remarkably, the decomposed retardance spectra [Fig. 5(b)] show a clear decrease for the colonic cancer. The changes in retardance are likely caused by the decreased collagen content in the colon cancer tissues [46] as the retardance effects are mainly attributed by the anisotropic orientation of collagen fibers in the concentric lamina propria and submucosa layers of the cross-section of a colon wall [9]. The prominent differences in the decomposed retardance between normal and cancer colon tissues [Fig. 6(c)] reconfirm the capability of Mueller-Matrix NIR diffuse reflectance spectroscopy for colonic cancer detection. Further PLS-DA and LOSCV analysis implemented on the 3 derived spectroscopic polarimetric metrics (Fig. 6) shows that colon cancer was identified with an accuracy of 90.0%, 91.7%, and 80.0%, respectively, by using diattentuation, depolarization, and retardance metrics (Table 1). Remarkably, the combination of the three polarization metrics with majority voting [38] provides an enhanced colonic cancer detection with an accuracy of 95.0% (sensitivity of 93.3%, and specificity of 96.7%), superior to using either of the three polarization metrics alone (Table 1).

 figure: Fig. 6

Fig. 6 (a) Mean diattentuation ± 1 standard error (SE) (shaded area); (b) Mean depolarization ± 1 SE, and (c) Mean retardance ± 1 SE, for the paired (normal (n = 30) vs cancer (n = 30)) colonic tissue.

Download Full Size | PDF

Tables Icon

Table 1. Diagnostic results of colonic cancer by using Mueller-matrix DR spectroscopy together with PLS-DA and LOSCV

4. Conclusions

We have developed an integrated Mueller-matrix NIR imaging and point-wise Mueller-matrix spectroscopy system for optical diagnosis and characterization of colonic cancer. Point-wise Mueller-Matrix spectra can be acquired immediately under Mueller-matrix imaging guidance. Significantly increased diattentuation while significantly reduced depolarization and retardance effects were observed associated with colonic cancer. Using the combined decomposed spectroscopic polarimetric metrics (i.e., diattentuation, depolarization, and retardance), colonic cancer can be identified with a high diagnostic accuracy (~95%). This work shows that Mueller-matrix NIR imaging and point-wise Mueller-Matrix NIR diffuse reflectance spectroscopy technique may open a new avenue for enhancing optical detection and diagnosis of colonic cancer in gastrointestinal tracts.

Acknowledgments

This work was supported by the National Medical Research Council (NMRC) (Grant Numbers: CIRG/1331/2012; BnB13dec037; BnB/0012c/2014), and the Academic Research Fund (AcRF)-Tier 2 from Ministry of Education (MOE) (Grant Number: MOE2014-T2-1-010), Singapore.

References and links

1. L. A. Torre, F. Bray, R. L. Siegel, J. Ferlay, J. Lortet-Tieulent, and A. Jemal, “Global cancer statistics, 2012,” CA Cancer J. Clin. 65(2), 87–108 (2015). [CrossRef]   [PubMed]  

2. N. R. O. D. Office, “Trends in Cancer Incidence in Singapore, 2010-2014,” Singapore Cancer Registry Interim Annual Report 1–56 (2015).

3. D. J. Robertson, D. A. Lieberman, S. J. Winawer, D. J. Ahnen, J. A. Baron, A. Schatzkin, A. J. Cross, A. G. Zauber, T. R. Church, P. Lance, E. R. Greenberg, and M. E. Martínez, “Colorectal cancers soon after colonoscopy: a pooled multicohort analysis,” Gut 63(6), 949–956 (2014). [CrossRef]   [PubMed]  

4. M. S. Bergholt, K. Lin, J. Wang, W. Zheng, H. Xu, Q. Huang, J. Ren, K. Y. Ho, M. Teh, S. Srivastava, B. Wong, K. G. Yeoh, and Z. Huang, “Simultaneous fingerprint and high-wavenumber fiber-optic Raman spectroscopy enhances real-time in vivo diagnosis of adenomatous polyps during colonoscopy,” J. Biophotonics 201, 1–10 (2015).

5. S. Winawer, R. Fletcher, D. Rex, J. Bond, R. Burt, J. Ferrucci, T. Ganiats, T. Levin, S. Woolf, D. Johnson, L. Kirk, S. Litin, C. Simmang, and Gastrointestinal Consortium Panel, “Colorectal cancer screening and surveillance: clinical guidelines and rationale-update based on new evidence,” Gastroenterology 124(2), 544–560 (2003). [CrossRef]   [PubMed]  

6. S. Alali and A. Vitkin, “Polarized light imaging in biomedicine: emerging Mueller matrix methodologies for bulk tissue assessment,” J. Biomed. Opt. 20(6), 061104 (2015). [CrossRef]   [PubMed]  

7. J. Qi, M. Ye, M. Singh, N. T. Clancy, and D. S. Elson, “Narrow band 3 × 3 Mueller polarimetric endoscopy,” Biomed. Opt. Express 4(11), 2433–2449 (2013). [CrossRef]   [PubMed]  

8. M. Sun, H. He, N. Zeng, E. Du, Y. Guo, S. Liu, J. Wu, Y. He, and H. Ma, “Characterizing the microstructures of biological tissues using Mueller matrix and transformed polarization parameters,” Biomed. Opt. Express 5(12), 4223–4234 (2014). [CrossRef]   [PubMed]  

9. I. Ahmad, M. Ahmad, K. Khan, S. Ashraf, S. Ahmad, and M. Ikram, “Ex vivo characterization of normal and adenocarcinoma colon samples by Mueller matrix polarimetry,” J. Biomed. Opt. 20(5), 056012 (2015). [CrossRef]   [PubMed]  

10. R. S. Gurjar, V. Backman, L. T. Perelman, I. Georgakoudi, K. Badizadegan, I. Itzkan, R. R. Dasari, and M. S. Feld, “Imaging human epithelial properties with polarized light-scattering spectroscopy,” Nat. Med. 7(11), 1245–1248 (2001). [CrossRef]   [PubMed]  

11. S. G. Demos and R. R. Alfano, “Optical polarization imaging,” Appl. Opt. 36(1), 150–155 (1997). [CrossRef]   [PubMed]  

12. X. Shao, W. Zheng, and Z. Huang, “Polarized near-infrared autofluorescence imaging combined with near-infrared diffuse reflectance imaging for improving colonic cancer detection,” Opt. Express 18(23), 24293–24300 (2010). [CrossRef]   [PubMed]  

13. R. D. Allen, J. Brault, and R. D. Moore, “A new method of polarization microscopic analysis I. Scanning with a birefringence detection system,” J. Cell Biol. 18(2), 223–235 (1963). [CrossRef]   [PubMed]  

14. S. B. Mehta, M. Shribak, and R. Oldenbourg, “Polarized light imaging of birefringence and diattenuation at high resolution and high sensitivity,” J. Opt. 15(9), 094007 (2013). [CrossRef]   [PubMed]  

15. K. Sokolov, R. Drezek, K. Gossage, and R. Richards-Kortum, “Reflectance spectroscopy with polarized light: is it sensitive to cellular and nuclear morphology,” Opt. Express 5(13), 302–317 (1999). [CrossRef]   [PubMed]  

16. D. S. Kliger and J. W. Lewis, Polarized light in optics and spectroscopy (Elsevier, 2012).

17. A. Pierangelo, A. Benali, M.-R. Antonelli, T. Novikova, P. Validire, B. Gayet, and A. De Martino, “Ex-vivo characterization of human colon cancer by Mueller polarimetric imaging,” Opt. Express 19(2), 1582–1593 (2011). [CrossRef]   [PubMed]  

18. M.-R. Antonelli, A. Pierangelo, T. Novikova, P. Validire, A. Benali, B. Gayet, and A. De Martino, “Mueller matrix imaging of human colon tissue for cancer diagnostics: how Monte Carlo modeling can help in the interpretation of experimental data,” Opt. Express 18(10), 10200–10208 (2010). [CrossRef]   [PubMed]  

19. S. L. Jacques, J. R. Roman, and K. Lee, “Imaging superficial tissues with polarized light,” Lasers Surg. Med. 26(2), 119–129 (2000). [CrossRef]   [PubMed]  

20. X. Han, H. Lui, D. I. McLean, and H. Zeng, “Near-infrared autofluorescence imaging of cutaneous melanins and human skin in vivo,” J. Biomed. Opt. 14(2), 024017 (2009). [CrossRef]   [PubMed]  

21. Z. Huang, H. Zeng, I. Hamzavi, A. Alajlan, E. Tan, D. I. McLean, and H. Lui, “Cutaneous melanin exhibiting fluorescence emission under near-infrared light excitation,” J. Biomed. Opt. 11(3), 034010 (2006). [CrossRef]   [PubMed]  

22. J. Wang, M. S. Bergholt, W. Zheng, and Z. Huang, “Development of a beveled fiber-optic confocal Raman probe for enhancing in vivo epithelial tissue Raman measurements at endoscopy,” Opt. Lett. 38(13), 2321–2323 (2013). [CrossRef]   [PubMed]  

23. K. Lin, W. Zheng, and Z. Huang, “Integrated autofluorescence endoscopic imaging and point-wise spectroscopy for real-time in vivo tissue measurements,” J. Biomed. Opt. 15(4), 040507 (2010). [CrossRef]   [PubMed]  

24. R. M. A. Azzam, “Photopolarimetric measurement of the Mueller matrix by Fourier analysis of a single detected signal,” Opt. Lett. 2(6), 148–150 (1978). [CrossRef]   [PubMed]  

25. D. H. Goldstein, “Mueller matrix dual-rotating retarder polarimeter,” Appl. Opt. 31(31), 6676–6683 (1992). [CrossRef]   [PubMed]  

26. K. Ichimoto, K. Shinoda, T. Yamamoto, and J. Kiyohara, “Photopolarimetric measurement system of Mueller matrix with dual rotating waveplates,” Publ. Natl. Astron. Obs. Jpn. 9, 11–19 (2006).

27. S.-Y. Lu and R. A. Chipman, “Interpretation of Mueller matrices based on polar decomposition,” J. Opt. Soc. Am. A 13(5), 1106–1113 (1996). [CrossRef]  

28. J. Wang, K. Lin, W. Zheng, K. Y. Ho, M. Teh, K. G. Yeoh, and Z. Huang, “Simultaneous fingerprint and high-wavenumber fiber-optic Raman spectroscopy improves in vivo diagnosis of esophageal squamous cell carcinoma at endoscopy,” Sci. Rep. 5, 12957 (2015).

29. H. J. Thomson, A. Busuttil, M. A. Eastwood, A. N. Smith, and R. A. Elton, “The submucosa of the human colon,” J. Ultrastruct. Mol. Struct. Res. 96(1-3), 22–30 (1986). [CrossRef]   [PubMed]  

30. N. Ghosh and I. A. Vitkin, “Tissue polarimetry: concepts, challenges, applications, and outlook,” J. Biomed. Opt. 16(11), 110801 (2011). [CrossRef]   [PubMed]  

31. M. Fleming, S. Ravula, S. F. Tatishchev, and H. L. Wang, “Colorectal carcinoma: Pathologic aspects,” J. Gastrointest. Oncol. 3(3), 153–173 (2012). [PubMed]  

32. J. M. Bueno, “Measurement of parameters of polarization in the living human eye using imaging polarimetry,” Vision Res. 40(28), 3791–3799 (2000). [CrossRef]   [PubMed]  

33. J. G. Bayly, V. B. Kartha, and W. H. Stevens, “The absorption spectra of liquid phase H 2 O, HDO and D 2 O from 0.7 μm to 10 μm,” Infrared Phys. 3(4), 211–222 (1963). [CrossRef]  

34. K. Sardana, and V. K. Garg, Lasers in Dermatological Practice (JP Medical Ltd, 2014).

35. S. H. Chung, A. E. Cerussi, S. I. Merritt, J. Ruth, and B. J. Tromberg, “Non-invasive tissue temperature measurements based on quantitative diffuse optical spectroscopy (DOS) of water,” Phys. Med. Biol. 55(13), 3753–3765 (2010). [CrossRef]   [PubMed]  

36. T. J. Saclarides, “Angiogenesis in colorectal cancer,” Surg. Clin. North Am. 77(1), 253–260 (1997). [CrossRef]   [PubMed]  

37. J. H. Ali, W. B. Wang, M. Zevallos, and R. R. Alfano, “Near infrared spectroscopy and imaging to probe differences in water content in normal and cancer human prostate tissues,” Technol. Cancer Res. Treat. 3(5), 491–497 (2004). [CrossRef]   [PubMed]  

38. J. Wang, K. Lin, W. Zheng, K. Y. Ho, M. Teh, K. G. Yeoh, and Z. Huang, “Fiber-optic Raman spectroscopy for in vivo diagnosis of gastric dysplasia,” Faraday Discuss. 1039, 151 (2015). [CrossRef]  

39. J. Mo, W. Zheng, J. J. Low, J. Ng, A. Ilancheran, and Z. Huang, “High wavenumber Raman spectroscopy for in vivo detection of cervical dysplasia,” Anal. Chem. 81(21), 8908–8915 (2009). [CrossRef]   [PubMed]  

40. S. Duraipandian, W. Zheng, J. Ng, J. J. Low, A. Ilancheran, and Z. Huang, “Simultaneous fingerprint and high-wavenumber confocal Raman spectroscopy enhances early detection of cervical precancer in vivo,” Anal. Chem. 84(14), 5913–5919 (2012). [CrossRef]   [PubMed]  

41. S. Duraipandian, W. Zheng, J. J. Low, J. Ng, A. Ilancheran, and Z. Huang, “Noninvasive analysis of hormonal variations and vagifem treatment effect on postmenopausal women using in vivo high wavenumber confocal Raman spectroscopy,” Analyst 138(14), 4120–4128 (2013). [CrossRef]  

42. R. Wolthuis, M. van Aken, K. Fountas, J. S. Robinson Jr, H. A. Bruining, and G. J. Puppels, “Determination of water concentration in brain tissue by Raman spectroscopy,” Anal. Chem. 73(16), 3915–3920 (2001). [CrossRef]   [PubMed]  

43. Y. S. Fawzy, M. Petek, M. Tercelj, and H. Zeng, “In vivo assessment and evaluation of lung tissue morphologic and physiological changes from non-contact endoscopic reflectance spectroscopy for improving lung cancer detection,” J. Biomed. Opt. 11(4), 044003 (2006). [CrossRef]   [PubMed]  

44. X. Shao, W. Zheng, and Z. Huang, “In vivo diagnosis of colonic precancer and cancer using near-infrared autofluorescence spectroscopy and biochemical modeling,” J. Biomed. Opt. 16(6), 067005 (2011). [CrossRef]   [PubMed]  

45. Y. S. Fawzy and H. Zeng, “Determination of scattering volume fraction and particle size distribution in the superficial layer of a turbid medium by using diffuse reflectance spectroscopy,” Appl. Opt. 45(16), 3902–3912 (2006). [CrossRef]   [PubMed]  

46. J. Turnay, N. Olmo, J. G. Gavilanes, and M. A. Lizarbe, “Collagen metabolism in human colon adenocarcinoma,” Connect. Tissue Res. 23(4), 251–260 (1989). [CrossRef]   [PubMed]  

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (6)

Fig. 1
Fig. 1 Schematic of the integrated Mueller-Matrix NIR imaging and point-wise Mueller-Matrix diffuse reflectance spectroscopy system developed for colonic tissue measurements: collimator (C); long pass filter (LP); lens (L); polarizer (P); quarter waveplate (QWP).
Fig. 2
Fig. 2 Example of an NIR diffuse reflectance image of a paired colon tissues (~6 x 3 x 3 mm3). Normal tissue is on the upper part, while cancer tissue is in the lower part of the image. One notes that the dark spot (~0.2 mm in diameter) is due to the refection of gold mirror coated in the point spectrum optical adapter (Fig. 1).
Fig. 3
Fig. 3 Normalized Muller matrix NIR images of the paired (normal vs cancer) colonic tissues. All the Mueller-Matrix elements (except m11) are normalized to m11.
Fig. 4
Fig. 4 The processed (a) diattentuation image, (b) depolarization image, and (c) retardance image. Intensity profiles (d-f) along the black dotted lines as indicated on the decomposed Mueller-Matrix NIR images in (a-c).
Fig. 5
Fig. 5 Representative 4 by 4 Muller-Matrix NIR DR spectra recorded of normal and cancer colonic tissue. (Normal: __________ ; Cancer: _________).
Fig. 6
Fig. 6 (a) Mean diattentuation ± 1 standard error (SE) (shaded area); (b) Mean depolarization ± 1 SE, and (c) Mean retardance ± 1 SE, for the paired (normal (n = 30) vs cancer (n = 30)) colonic tissue.

Tables (1)

Tables Icon

Table 1 Diagnostic results of colonic cancer by using Mueller-matrix DR spectroscopy together with PLS-DA and LOSCV

Equations (12)

Equations on this page are rendered with MathJax. Learn more.

I= a 0 + n=1 12 ( a n cosnwt+ b n sinnwt)
M = [ m 11 m 12 m 13 m 14 m 21 m 22 m 23 m 24 m 31 m 32 m 33 m 34 m 41 m 42 m 43 m 44 ]
M = M Δ M R M D
D ¯ = 1 m 11 ( m 12 , m 13 , m 14 )
D= 1 m 11 m 12 2 + m 13 2 + m 14 2
M D = m 11 [ 1 D T D m D ]
M'=M M D 1
M Δ M R =M M D 1 =[ 1 0 T P Δ m Δ ][ 1 0 T 0 m R ] =[ 1 0 T P Δ m Δ m R ] =[ 1 0 T P Δ m Δ m R ] =[ 1 0 T P Δ m' ]=M'
m Δ =± [ m' (m') T +( λ 1 λ 2 + λ 2 λ 3 + λ 3 λ 1 )I ] 1 ×[ ( λ 1 + λ 2 + λ 3 )m' (m') T + λ 1 λ 2 λ 3 I ]
M R = M Δ 1 M'
Δ=1 | trace( M Δ 1) | 3
R= cos 1 [ trace( M R ) 2 1]
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.