The Resolution Matrix in Tomographic Multiplexing: Optimization of Inter-parameter Cross-talk, Relative Quantitation and Localization

Objective: We use a resolution matrix-based Bayesian framework to compare inversion methods for tomographic fluorescence lifetime multiplexing in a diffuse medium, such as biological tissue. Methods: We consider three inversion methods; an asymptotic time domain (ATD) approach, based on a multiexponential analysis of time domain data, a direct time domain (DTD) approach, which is a minimum error solution, and a cross-talk constrained time domain (CCTD) inversion, which is a solution to an optimization problem that minimizes both error and cross-talk. We compare these methods using Monte Carlo simulations and time domain fluorescence measurements with tissue-mimicking phantoms. Results: The ATD approach provides high accuracy of relative quantitation and spatial localization of two fluorophores embedded in a 18-mm thick turbid medium, with concentration ratios of up to 1:4.25. DTD leads to significant errors in relative quantitation and localization. CCTD provides improved quantitation accuracy over DTD, and better spatial resolution compared to ATD. We present a rigorous theoretical basis for these results and provide a complete derivation of the CCTD estimator. The Bayesian analysis also leads to a formula for rapid computation of the DTD inverse operator for large-scale tomography measurements. Conclusion: The ATD and CCTD inversion methods provide significant advantages over DTD for accurately estimating multiple overlapping fluorophores. Significance: Time domain fluorescence tomography, using zero cross-talk estimators, can serve as a powerful tool for quantifying multiple fluorescently labeled biological processes. The Bayesian framework presented here can be applied to general multiparameter inverse problems for the quantitative estimation of multiple overlapping parameters.