RESUMO

The deconvolution process is a key step for quantitative evaluation of fluorescence lifetime imaging microscopy (FLIM) samples. By this process, the fluorescence impulse responses (FluoIRs) of the sample are decoupled from the instrument response (InstR). In blind deconvolution estimation (BDE), the FluoIRs and InstR are jointly extracted from a dataset with minimal a priori information. In this work, two BDE algorithms are introduced based on linear combinations of multi-exponential functions to model each FluoIR in the sample. For both schemes, the InstR is assumed with a free-form and a sparse structure. The local perspective of the BDE methodology assumes that the characteristic parameters of the exponential functions (time constants and scaling coefficients) are estimated based on a single spatial point of the dataset. On the other hand, the same exponential functions are used in the whole dataset in the global perspective, and just the scaling coefficients are updated for each spatial point. A least squares formulation is considered for both BDE algorithms. To overcome the nonlinear interaction in the decision variables, an alternating least squares (ALS) methodology iteratively solves both estimation problems based on non-negative and constrained optimizations. The validation stage considered first synthetic datasets at different noise types and levels, and a comparison with the standard deconvolution techniques with a multi-exponential model for FLIM measurements, as well as, with two BDE methodologies in the state of the art: Laguerre basis, and exponentials library. For the experimental evaluation, fluorescent dyes and oral tissue samples were considered. Our results show that local and global perspectives are consistent with the standard deconvolution techniques, and they reached the fastest convergence responses among the BDE algorithms with the best compromise in FluoIRs and InstR estimation errors.

Assuntos

RESUMO

Spectral unmixing is the process of breaking down data from a sample into its basic components and their abundances. Previous work has been focused on blind unmixing of multi-spectral fluorescence lifetime imaging microscopy (m-FLIM) datasets under a linear mixture model and quadratic approximations. This method provides a fast linear decomposition and can work without a limitation in the maximum number of components or end-members. Hence this work presents an interactive software which implements our blind end-member and abundance extraction (BEAE) and quadratic blind linear unmixing (QBLU) algorithms in Matlab. The options and capabilities of our proposed software are described in detail. When the number of components is known, our software can estimate the constitutive end-members and their abundances. When no prior knowledge is available, the software can provide a completely blind solution to estimate the number of components, the end-members and their abundances. The characterization of three case studies validates the performance of the new software: ex-vivo human coronary arteries, human breast cancer cell samples, and in-vivo hamster oral mucosa. The software is freely available in a hosted webpage by one of the developing institutions, and allows the user a quick, easy-to-use and efficient tool for multi/hyper-spectral data decomposition.

Assuntos

RESUMO

Time-deconvolution of the instrument response from fluorescence lifetime imaging microscopy (FLIM) data is usually necessary for accurate fluorescence lifetime estimation. In many applications, however, the instrument response is not available. In such cases, a blind deconvolution approach is required. An iterative methodology is proposed to address the blind deconvolution problem departing from a dataset of FLIM measurements. A linear combination of a base conformed by Laguerre functions models the fluorescence impulse response of the sample at each spatial point in our formulation. Our blind deconvolution estimation (BDE) algorithm is formulated as a quadratic approximation problem, where the decision variables are the samples of the instrument response and the scaling coefficients of the basis functions. In the approximation cost function, there is a bilinear dependence on the decision variables. Hence, due to the nonlinear nature of the estimation process, an alternating least-squares scheme iteratively solves the approximation problem. Our proposal searches for the samples of the instrument response with a global perspective, and the scaling coefficients of the basis functions locally at each spatial point. First, the iterative methodology relies on a least-squares solution for the instrument response, and quadratic programming for the scaling coefficients applied just to a subset of the measured fluorescence decays to initially estimate the instrument response to speed up the convergence. After convergence, the final stage computes the fluorescence impulse response at all spatial points. A comprehensive validation stage considers synthetic and experimental FLIM datasets of ex vivo atherosclerotic plaques and human breast cancer cell samples that highlight the advantages of the proposed BDE algorithm under different noise and initial conditions in the iterative scheme and parameters of the proposal.

Assuntos

RESUMO

In this paper, we investigate novel low-dimensional and model-free representations for multi-spectral fluorescence lifetime imaging microscopy (m-FLIM) data. We depart from the classical definition of the phasor in the complex plane to propose the extended output phasor (EOP) and extended phasor (EP) for multi-spectral information. The frequency domain properties of the EOP and EP are analytically studied based on a multiexponential model for the impulse response of the imaged tissue. For practical implementations, the EOP is more appealing since there is no need to perform deconvolution of the instrument response from the measured m-FLIM data, as in the case of EP. Our synthetic and experimental evaluations with m-FLIM datasets of human coronary atherosclerotic plaques show that low frequency indexes have to be employed for a distinctive representation of the EOP and EP, and to reduce noise distortion. The tissue classification of the m-FLIM datasets by EOP and EP also improves with low frequency indexes, and does not present significant differences by using either phasor.

RESUMO

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.

Assuntos

RESUMO

Multi-Spectral Fluorescent Lifetime Imaging Microscopy (m-FLIM) is a technique that aims to perform noninvasive in situ clinical diagnosis of several diseases. It measures the endogenous fluorescence of molecules, recording their lifetime decay in different wavelength bands. This signal is a mixed response of multiple fluorescent components present in a tissue sample. The goal is to decompose the mixture and estimate the proportional contributions of its constituents. Estimation of such quantitative description will help to characterize the molecular constitution of a given sample. This paper presents a new method to estimate the abundances of multiple components present in a mixture measured using m-FLIM data. It provides a closed-form solution under the fully constrained linear unmixing model and assuming the number of components as well as their ideal lifetime decays are known. Its performance is tested using synthetic samples with three components, where performance can be measured accurately and the percentage error is around 6%. The algorithm was also validated performing unmixing of ex vivo data samples from atherosclerotic human tissue containing collagen, elastin and low-density lipoproteins. These experiments were validated against ground-truth maps, which only give a quantitative description, and the estimated accuracy was around 88%.

Assuntos