Blind instrument response function identification from fluorescence decays.

Time-resolved fluorescence spectroscopy plays a crucial role when studying dynamic properties of complex photochemical systems. Nevertheless, the analysis of measured time decays and the extraction of exponential lifetimes often requires either the experimental assessment or the modeling of the instrument response function (IRF). However, the intrinsic nature of the IRF in the measurement process, which may vary across measurements due to chemical and instrumental factors, jeopardizes the results obtained by reconvolution approaches. In this paper, we introduce a novel methodology, called blind instrument response function identification (BIRFI), which enables the direct estimation of the IRF from the collected data. It capitalizes on the properties of single exponential signals to transform a deconvolution problem into a well-posed system identification problem. To delve into the specifics, we provide a step-by-step description of the BIRFI method and a protocol for its application to fluorescence decays. The performance of BIRFI is evaluated using simulated and time-correlated single-photon counting data. Our results demonstrate that the BIRFI methodology allows an accurate recovery of the IRF, yielding comparable or even superior results compared with those obtained with experimental IRFs when they are used for reconvolution by parametric model fitting.

Kernelizing: A way to increase accuracy in trilinear decomposition analysis of multiexponential signals.

The unmixing of multiexponential decay signals into monoexponential components using soft modelling approaches is a challenging task due to the strong correlation and complete window overlap of the profiles. To solve this problem, slicing methodologies, such as PowerSlicing, tensorize the original data matrix into a three-way data array that can be decomposed based on trilinear models providing unique solutions. Satisfactory results have been reported for different types of data, e.g., nuclear magnetic resonance or time-resolved fluorescence spectra. However, when decay signals are described by only a few sampling (time) points, a significant degradation of the results can be observed in terms of accuracy and precision of the recovered profiles. In this work, we propose a methodology called Kernelizing that provides a more efficient way to tensorize data matrices of multiexponential decays. Kernelizing relies on the invariance of exponential decays, i.e., when convolving a monoexponential decaying function with any positive function of finite width (hereafter called "kernel"), the shape of the decay (determined by the characteristic decay constant) remains unchanged and only the preexponential factor varies. The way preexponential factors are affected across the sample and time modes is linear, and it only depends on the kernel used. Thus, using kernels of different shapes, a set of convolved curves can be obtained for every sample, and a three-way data array generated, for which the modes are sample, time and kernelizing effect. This three-way array can be afterwards analyzed by a trilinear decomposition method, such as PARAFAC-ALS, to resolve the underlying monoexponential profiles. To validate this new approach and assess its performance, we applied Kernelizing to simulated datasets, real time-resolved fluorescence spectra collected on mixtures of fluorophores and fluorescence-lifetime imaging microscopy data. When the measured multiexponential decays feature few sampling points (down to fifteen), more accurate trilinear model estimates are obtained than when using slicing methodologies.

Linear unmixing protocol for hyperspectral image fusion analysis applied to a case study of vegetal tissues.

Hyperspectral imaging (HSI) is a useful non-invasive technique that offers spatial and chemical information of samples. Often, different HSI techniques are used to obtain complementary information from the sample by combining different image modalities (Image Fusion). However, issues related to the different spatial resolution, sample orientation or area scanned among platforms need to be properly addressed. Unmixing methods are helpful to analyze and interpret the information of HSI related to each of the components contributing to the signal. Among those, Multivariate Curve Resolution-Alternating Least Squares (MCR-ALS) offers very suitable features for image fusion, since it can easily cope with multiset structures formed by blocks of images coming from different samples and platforms and allows the use of optional and diverse constraints to adapt to the specific features of each HSI employed. In this work, a case study based on the investigation of cross-sections from rice leaves by Raman, synchrotron infrared and fluorescence imaging techniques is presented. HSI of these three different techniques are fused for the first time in a single data structure and analyzed by MCR-ALS. This example is challenging in nature and is particularly suitable to describe clearly the necessary steps required to perform unmixing in an image fusion context. Although this protocol is presented and applied to a study of vegetal tissues, it can be generally used in many other samples and combinations of imaging platforms.

3D and 4D Image Fusion: Coping with Differences in Spectroscopic Modes among Hyperspectral Images.

Image fusion is often oriented to solve differences in spatial scale and orientation among different spectroscopic platforms. However, an additional problem arises when the nature of the spectroscopic information differs in dimensionality as well. Indeed, most imaging systems, e.g., Raman, IR, MS, etc., allow acquisition of 3D images, with a linear spectrum per pixel, but new platforms have emerged, such as the recent excitation-emission fluorescence imaging platforms that provide 4D images, with a 2D spectral landscape per pixel. A proper 3D/4D image fusion needs to take into account the difference in the dimension of the spectral information and in the underlying models of both measurements (bilinear for 3D images and trilinear for 4D images). This work solves this image fusion problem through a new dedicated variant of the multivariate curve resolution-alternating least squares (MCR-ALS) algorithm for multiset analysis based on the incorporation of a hybrid bilinear/trilinear model that can handle the image fused structure preserving the natural behavior of the 3D and 4D imaging techniques coupled. The example is illustrated on the fusion of real 3D Raman and 4D fluorescence images recorded on cross sections of rice leaf samples.

Understanding the Formation of Heartwood in Larch Using Synchrotron Infrared Imaging Combined With Multivariate Analysis and Atomic Force Microscope Infrared Spectroscopy.

Formation of extractive-rich heartwood is a process in live trees that make them and the wood obtained from them more resistant to fungal degradation. Despite the importance of this natural mechanism, little is known about the deposition pathways and cellular level distribution of extractives. Here we follow heartwood formation in Larix gmelinii var. Japonica by use of synchrotron infrared images analyzed by the unmixing method Multivariate Curve Resolution - Alternating Least Squares (MCR-ALS). A subset of the specimens was also analyzed using atomic force microscopy infrared spectroscopy. The main spectral changes observed in the transition zone when going from sapwood to heartwood was a decrease in the intensity of a peak at approximately 1660 cm-1 and an increase in a peak at approximately 1640 cm-1. There are several possible interpretations of this observation. One possibility that is supported by the MCR-ALS unmixing is that heartwood formation in larch is a type II or Juglans-type of heartwood formation, where phenolic precursors to extractives accumulate in the sapwood rays. They are then oxidized and/or condensed in the transition zone and spread to the neighboring cells in the heartwood.