Work

Approaches for modelling autocorrelation function and data processing in time-domain diffuse correlation spectroscopy

Year: 2025

Type: article

Abstract: Time-domain diffuse correlation spectroscopy (TD-DCS) is a non-invasive optical technique for measuring tissue blood flow. Recovering the blood flow index ( αD b ) requires accurate modelling of the normalised electric field autocorrelation function ( g 1 ), and an optimised data processing approach to minimise noise. We quantitatively compared four modelling approaches for g 1 : (i) using momentum transfer ( Y ) and pathlengths ( L ) from Monte Carlo (MC) simulations, (ii) using L only, (iii) applying an analytical solution of the photon diffusion equation (DE) in time domain, and (iv) applying an analytical solution of the correlation diffusion equation (CDE) in steady state. The second and third approaches use solutions in near-infrared spectroscopy (NIRS) for modelling g 1 in DCS by assuming Y = μ′ s L . We computed g 1 curves using the first approach, considered the gold standard, and recovered αD b using the other three approaches for various source-detector distances ( ρ ) and scattering coefficients ( μ′ s ). Also, we investigated how the correlator time bin width ( T bin ), which is an adjustable parameter in data processing, affects the standard deviation of g 1 (or the normalised intensity autocorrelation function g 2 ). We used a more convenient version of the noise equation expressed as a function of g 1 (or g 2 ), which removes the need to know the decay rate. When using photons detected after ∼0.5 ns, all four approaches produced nearly identical g 1 curves. Using all detected photons, the DE solution produced negligible errors (up to ∼2%) in the recovered αD b across various ρ and μ′ s , while using L from MC simulations resulted in larger errors (up to ∼9% at ρ = 5 mm and ∼1.5% at ρ = 30 mm). The analysis of the probability distributions P ( Y ) and P ( μ′ s L ) explained these differences. As expected, the standard deviation of g 1 (or g 2 ) can be reduced during data processing by increasing T bin . To achieve the lowest standard deviation, T bin should be longer than the inverse of the photon count rate, indicating that the optimal T bin may vary across different time gates. The results provide quantitative insights into modelling g 1 (or g 2 ), and provide a direct guideline for minimising the standard deviation of g 1 (or g 2 ) in data processing.

Source: Biomedical Optics Express

DOI: 10.1364/boe.577785

Authors: Aleh Sudakou, Ilias Tachtsidis, Michal Kacprzak, Adam Liebert, Stanislaw Wojtkiewicz

Institutions: Institute of Biocybernetics and Biomedical Engineering, Polish Academy of Sciences, University College London

Language: English

Topic: Optical Imaging and Spectroscopy Techniques

Subfield: Radiology, Nuclear Medicine and Imaging

Field: Medicine

Domain: Health Sciences

Open Access status: gold

Estimated APC paid: $1,913

Funders: Narodowe Centrum Nauki, Narodowe Centrum Badań i Rozwoju, Innovate UK