This paper proposes a new blind end-member and abundance extraction (BEAE) method for multispectral fluorescence lifetime imaging microscopy (m-FLIM) data. The chemometrical analysis relies on an iterative estimation of the fluorescence decay end-members and their abundances. The proposed method is based on a linear mixture model with positivity and sum-to-one restrictions on the abundances and end-members to compensate for signature variability. The synthesis procedure depends on a quadratic optimization problem, which is solved by an alternating least-squares structure over convex sets. The BEAE strategy only assumes that the number of components in the analyzed sample is known a spriori. The proposed method is first validated by using synthetic m-FLIM datasets at 15, 20, and 25 dB signal-to-noise ratios. The samples simulate the mixed response of tissue containing multiple fluorescent intensity decays. Furthermore, the results were also validated with six m-FLIM datasets from fresh postmortem human coronary atherosclerotic plaques. A quantitative evaluation of the BEAE was made against two popular techniques: minimum volume constrained nonnegative matrix factorization (MVC-NMF) and multivariate curve resolution-alternating least-squares (MCR-ALS). Our proposed method (BEAE) was able to provide more accurate estimations of the end-members: 0.32% minimum relative error and 13.82% worst-case scenario, despite different initial conditions in the iterative optimization procedure and noise effect. Meanwhile, MVC-NMF and MCR-ALS presented more variability in estimating the end-members: 0.35% and 0.34% for minimum errors and 15.31% and 13.25% in the worst-case scenarios, respectively. This tendency was also maintained for the abundances, where BEAE obtained 0.05 as the minimum absolute error and 0.12 in the worst-case scenario; MCR-ALS and MVC-NMF achieved 0.04 and 0.06 for the minimum absolute errors, and 0.15 and 0.17 under the worst-case conditions, respectively. In addition, the average computation time was evaluated for the synthetic datasets, where MVC-NMF achieved the fastest time, followed by BEAE and finally MCR-ALS. Consequently, BEAE improved MVC-NMF in convergence to a local optimal solution and robustness against signal variability, and it is roughly 3.6 time faster than MCR-ALS.

This paper presents a new unmixing methodology of multispectral fluorescence lifetime imaging microscopy (m-FLIM) data, in which the spectrum is defined as the combination of time-domain fluorescence decays at multiple emission wavelengths. The method is based on a quadratic constrained optimization (CO) algorithm that provides a closed-form solution under equality and inequality restrictions. In this paper, it is assumed that the time-resolved fluorescence spectrum profiles of the constituent components are linearly independent and known a priori. For comparison purposes, the standard least squares (LS) solution and two constrained versions nonnegativity constrained least squares (NCLS) and fully constrained least squares (FCLS) (Heinz and Chang, 2001) are also tested. Their performance was evaluated by using synthetic simulations, as well as imaged samples from fluorescent dyes and ex vivo tissue. In all the synthetic evaluations, the CO obtained the best accuracy in the estimations of the proportional contributions. CO could achieve an improvement ranging between 41% and 59% in the relative error compared to LS, NCLS, and FCLS at different signal-to-noise ratios. A liquid mixture of fluorescent dyes was also prepared and imaged in order to provide a controlled scenario with real data, where CO and FCLS obtained the best performance. The CO and FCLS were also tested with 20 ex vivo samples of human coronary arteries, where the expected concentrations are qualitatively known. A certainty measure was employed to assess the confidence in the estimations made by each algorithm. The experiments confirmed a better performance of CO, since this method is optimal with respect to equality and inequality restrictions in the linear unmixing formulation. Thus, the evaluation showed that CO achieves an accurate characterization of the samples. Furthermore, CO is a computational efficient alternative to estimate the abundance of components in m-FLIM data, since a global optimal solution is always guaranteed in a closed form.

This study proposes an automatic method for the sleep-wake staging in normal and pathologic sleep based only on respiratory effort acquired from a Pressure Bed Sensor (PBS). Motion and respiratory movements were obtained through a PBS and sleep-wake staging was evaluated from those time series. 20 all night polysomnographies, with annotations, used as gold standard and the time series coming from the PBS were used to develop and to evaluate the automatic wake-sleep staging. The database was built up by: 10 healthy subjects and 10 patients with severe sleep apnea. The agreement of the statistical measures between the automatic classification and the human scoring were: 83.59 ± 6.79 of sensitivity, 83.60 ± 15.13 of specificity and 81.91 ± 6.36 of accuracy. These results suggest that some important indexes, such as sleep efficiency, could be computed through a contactless technique.

Multispectral fluorescence lifetime imaging (m-FLIM) can potentially allow identifying the endogenous fluorophores present in biological tissue. Quantitative description of such data requires estimating the number of components in the sample, their characteristic fluorescent decays, and their relative contributions or abundances. Unfortunately, this inverse problem usually requires prior knowledge about the data, which is seldom available in biomedical applications. This work presents a new methodology to estimate the number of potential endogenous fluorophores present in biological tissue samples from time-domain m-FLIM data. Furthermore, a completely blind linear unmixing algorithm is proposed. The method was validated using both synthetic and experimental m-FLIM data. The experimental m-FLIM data include in-vivo measurements from healthy and cancerous hamster cheek-pouch epithelial tissue, and ex-vivo measurements from human coronary atherosclerotic plaques. The analysis of m-FLIM data from in-vivo hamster oral mucosa identified healthy from precancerous lesions, based on the relative concentration of their characteristic fluorophores. The algorithm also provided a better description of atherosclerotic plaques in term of their endogenous fluorophores. These results demonstrate the potential of this methodology to provide quantitative description of tissue biochemical composition.

The increasing prevalence of vascular diseases encourages the development of minimally invasive approaches to assess tissue perfusion. A significant challenge facing current state-of-the-art methods is their validation against clinical data. In this study, we introduce an open-source database designed to evaluate tissue perfusion during the application of an occlusion protocol. The database comprises sequences of multi-spectral images (visible and near-infrared region) from the subjects’ predominant hand and their photoplethysmography data for validation. Our study recruited 45 healthy participants, including 21 females, with an age range between 18-24 years old (standard deviation equal to 1.73). The database was evaluated using two methods for estimating skin perfusion parameters based on multi-spectral images: a Kubelka-Munk model, and a linear regression. Meanwhile, for validation purposes, the changes in oxygenated and deoxygenated hemoglobin were evaluated by photoplethysmography data as baseline perfusion parameters. The Pearson correlation between plethysmography-based perfusion parameters and those extracted from multi-spectral images was evaluated in all cases as a validation metric. Our findings demonstrated a strong Pearson correlation (ρ > 0.7) between changes in oxygenated and deoxygenated hemoglobin and multi-spectral based perfusion parameters, suggesting that the database is useful for further research related to in-vivo perfusion assessment. The primary objective of this database is to provide open-source data from a controlled occlusion protocol to evaluate new approaches based on multi-spectral images in the visible and near-infrared regions. In addition, the validation by photoplethysmography data facilitates the development and assessment of innovative tissue perfusion estimation techniques.