Blind deconvolution estimation of fluorescence measurements through quadratic programming

Computational tools for inversion and uncertainty estimation in respirometry

In many physiological systems, real-time endogeneous and exogenous signals in living organisms provide critical information and interpretations of physiological functions; however, these signals or variables of interest are not directly accessible and must be estimated from noisy, measured signals. In this paper, we study an inverse problem of recovering gas exchange signals of animals placed in a flow-through respirometry chamber from measured gas concentrations. For large-scale experiments (e.g., long scans with high sampling rate) that have many uncertainties (e.g., noise in the observations or an unknown impulse response function), this is a computationally challenging inverse problem. We first describe various computational tools that can be used for respirometry reconstruction and uncertainty quantification when the impulse response function is known. Then, we address the more challenging problem where the impulse response function is not known or only partially known. We describe nonlinear optimization methods for reconstruction, where both the unknown model parameters and the unknown signal are reconstructed simultaneously. Numerical experiments show the benefits and potential impacts of these methods in respirometry.

Blind deconvolution estimation by multi-exponential models and alternated least squares approximations: Free-form and sparse approach

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.

Extended Blind End-member and Abundance Extraction for Biomedical Imaging Applications

In some applications of biomedical imaging, a linear mixture model can represent the constitutive elements (end-members) and their contributions (abundances) per pixel of the image. In this work, the extended blind end-member and abundance extraction (EBEAE) methodology is mathematically formulated to address the blind linear unmixing (BLU) problem subject to positivity constraints in optical measurements. The EBEAE algorithm is based on a constrained quadratic optimization and an alternated least-squares strategy to jointly estimate end-members and their abundances. In our proposal, a local approach is used to estimate the abundances of each end-member by maximizing their entropy, and a global technique is adopted to iteratively identify the end-members by reducing the similarity among them. All the cost functions are normalized, and four initialization approaches are suggested for the end-members matrix. Synthetic datasets are used first for the EBEAE validation at different noise types and levels, and its performance is compared to state-of-the-art algorithms in BLU. In a second stage, three experimental biomedical imaging applications are addressed with EBEAE: m-FLIM for chemometric analysis in oral cavity samples, OCT for macrophages identification in post-mortem artery samples, and hyper-spectral images for in-vivo brain tissue classification and tumor identification. In our evaluations, EBEAE was able to provide a quantitative analysis of the samples with none or minimal a priori information.

Optical fiber-based dispersion for spectral discrimination in fluorescence lifetime imaging systems

The excited state lifetime of a fluorophore together with its fluorescence emission spectrum provide information that can yield valuable insights into the nature of a fluorophore and its microenvironment. However, it is difficult to obtain both channels of information in a conventional scheme as detectors are typically configured either for spectral or lifetime detection. We present a fiber-based method to obtain spectral information from a multiphoton fluorescence lifetime imaging (FLIM) system. This is made possible using the time delay introduced in the fluorescence emission path by a dispersive optical fiber coupled to a detector operating in time-correlated single-photon counting mode. This add-on spectral implementation requires only a few simple modifications to any existing FLIM system and is considerably more cost-efficient compared to currently available spectral detectors.

Temporal binning of time-correlated single photon counting data improves exponential decay fits and imaging speed

Time-correlated single photon counting (TCSPC) enables acquisition of fluorescence lifetime decays with high temporal resolution within the fluorescence decay. However, many thousands of photons per pixel are required for accurate lifetime decay curve representation, instrument response deconvolution, and lifetime estimation, particularly for two-component lifetimes. TCSPC imaging speed is inherently limited due to the single photon per laser pulse nature and low fluorescence event efficiencies (<10%) required to reduce bias towards short lifetimes. Here, simulated fluorescence lifetime decays are analyzed by SPCImage and SLIM Curve software to determine the limiting lifetime parameters and photon requirements of fluorescence lifetime decays that can be accurately fit. Data analysis techniques to improve fitting accuracy for low photon count data were evaluated. Temporal binning of the decays from 256 time bins to 42 time bins significantly (p<0.0001) improved fit accuracy in SPCImage and enabled accurate fits with low photon counts (as low as 700 photons/decay), a 6-fold reduction in required photons and therefore improvement in imaging speed. Additionally, reducing the number of free parameters in the fitting algorithm by fixing the lifetimes to known values significantly reduced the lifetime component error from 27.3% to 3.2% in SPCImage (p<0.0001) and from 50.6% to 4.2% in SLIM Curve (p<0.0001). Analysis of nicotinamide adenine dinucleotide-lactate dehydrogenase (NADH-LDH) solutions confirmed temporal binning of TCSPC data and a reduced number of free parameters improves exponential decay fit accuracy in SPCImage. Altogether, temporal binning (in SPCImage) and reduced free parameters are data analysis techniques that enable accurate lifetime estimation from low photon count data and enable TCSPC imaging speeds up to 6x and 300x faster, respectively, than traditional TCSPC analysis.