Cramer-Rao analysis of steady-state and time-domain fluorescence diffuse optical imaging

System configuration optimization for mesoscopic fluorescence molecular tomography

Tissue engineering applications demand 3D, non-invasive, and longitudinal assessment of bioprinted constructs. Current emphasis is on developing tissue constructs mimicking in vivo conditions; however, these are increasingly challenging to image as they are typically a few millimeters thick and turbid, limiting the usefulness of classical fluorescence microscopic techniques. For such applications, we developed a Mesoscopic Fluorescence Molecular Tomography methodology that collects high information content data to enable high-resolution tomographic reconstruction of fluorescence biomarkers at millimeters depths. This imaging approach is based on an inverse problem; hence, its imaging performances are dependent on critical technical considerations including optode sampling, forward model design and inverse solver parameters. Herein, we investigate the impact of the optical system configuration parameters, including detector layout, number of detectors, combination of detector and source numbers, and scanning mode with uncoupled or coupled source and detector array, on the 3D imaging performances. Our results establish that an MFMT system with a 2D detection chain implemented in a de-scanned mode provides the optimal imaging reconstruction performances.

Cramér-Rao analysis of lifetime estimations in time-resolved fluorescence microscopy

Measuring the lifetime of fluorescent emitters by time-correlated single photon counting (TCSPC) is a routine procedure in many research areas spanning from nanophotonics to biology. The precision of such measurement depends on the number of detected photons but also on the various sources of noise arising from the measurement process. Using Fisher information theory, we calculate the lower bound on the precision of lifetime estimations for mono-exponential and bi-exponential distributions. We analyse the dependence of the lifetime estimation precision on experimentally relevant parameters, including the contribution of a non-uniform background noise and the instrument response function (IRF) of the setup. We also provide an open-source code to determine the lower bound on the estimation precision for any experimental conditions. Two practical examples illustrate how this tool can be used to reach optimal precision in time-resolved fluorescence microscopy.

On the use of the Cramér-Rao lower bound for diffuse optical imaging system design

We evaluated the potential of the Cramér-Rao lower bound (CRLB) to serve as a design metric for diffuse optical imaging systems. The CRLB defines the best achievable precision of any estimator for a given data model; it is often used in the statistical signal processing community for feasibility studies and system design. Computing the CRLB requires inverting the Fisher information matrix (FIM), however, which is usually ill-conditioned (and often underdetermined) in the case of diffuse optical tomography (DOT). We regularized the FIM by assuming that the inhomogeneity to be imaged was a point target and assessed the ability of point-target CRLBs to predict system performance in a typical DOT setting in silico. Our reconstructions, obtained with a common iterative algebraic technique, revealed that these bounds are not good predictors of imaging performance across different system configurations, even in a relative sense. This study demonstrates that agreement between the trends predicted by the CRLBs and imaging performance obtained with reconstruction algorithms that rely on a different regularization approach cannot be assumed a priori. Moreover, it underscores the importance of taking into account the intended regularization method when attempting to optimize source-detector configurations.

Truncated Fourier-series approximation of the time-domain radiative transfer equation using finite elements

The radiative transfer equation (RTE) is widely accepted to accurately describe light transport in a medium with scattering particles, and it has been successfully applied as a light-transport model, for example, in diffuse optical tomography. Due to the computationally expensive nature of the RTE, most of these applications have been in the frequency domain. In this paper, an efficient solution method for the time-domain RTE is proposed. The method is based on solving the frequency-domain RTE at multiple modulation frequencies and using the Fourier-series representation of the radiance to obtain approximation of the time-domain solution. The approach is tested with simulations. The results show that the method can be used to obtain the solution of the time-domain RTE with good accuracy and with significantly fewer computational resources than are needed in the direct time-domain solution.