Efficient quantitative hyperspectral image unmixing method for large-scale Raman micro-spectroscopy data analysis

Vibrational micro-spectroscopy is a powerful optical tool, providing a non-invasive label-free chemically specific imaging for many chemical and biomedical applications. However, hyperspectral image produced by Raman micro-spectroscopy typically consists of thousands discrete pixel points, each having individual Raman spectrum at thousand wavenumbers, and therefore requires appropriate image unmixing computational methods to retrieve non-negative spatial concentration and corresponding non-negative spectra of the image biochemical constituents. Here, we present a new efficient Quantitative Hyperspectral Image Unmixing (Q-HIU) method for large-scale Raman micro-spectroscopy data analysis. This method enables to simultaneously analyse multi-set Raman hyperspectral images in three steps: (i) Singular Value Decomposition with innovative Automatic Divisive Correlation which autonomously filters spatially and spectrally uncorrelated noise from data; (ii) a robust subtraction of fluorescent background from the data using a newly developed algorithm called Bottom Gaussian Fitting; (iii) an efficient Quantitative Unsupervised/Partially Supervised Non-negative Matrix Factorization method, which rigorously retrieves non-negative spatial concentration maps and spectral profiles of the samples' biochemical constituents with no a priori information or when one or several samples' constituents are known. As compared with state-of-the-art methods, our approach allows to achieve significantly more accurate results and efficient quantification with several orders of magnitude shorter computational time as verified on both artificial and real experimental data. We apply Q-HIU to the analysis of large-scale Raman hyperspectral images of human atherosclerotic aortic tissues and our results show a proof-of-principle for the proposed method to retrieve and quantify the biochemical composition of the tissues, consisting of both high and low concentrated compounds. Along with the established hallmarks of atherosclerosis including cholesterol/cholesterol ester, triglyceride and calcium hydroxyapatite crystals, our Q-HIU allowed to identify the significant accumulations of oxidatively modified lipids co-localizing with the atherosclerotic plaque lesions in the aortic tissues, possibly reflecting the persistent presence of inflammation and oxidative damage in these regions, which are in turn able to promote the disease pathology. For minor chemical components in the diseased tissues, our Q-HIU was able to detect the signatures of calcium hydroxyapatite and β-carotene with relative mean Raman concentrations as low as 0.09% and 0.04% from the original Raman intensity matrix with noise and fluorescent background contributions of 3% and 94%, respectively.