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

Evaluation of temporal moments and Fourier transformed data in time-domain diffuse optical tomography

Time-domain diffuse optical tomography (TD-DOT) uses near-infrared pulsed lasers as light sources to measure time-varying exitance on the boundary of the target. These are used to estimate optical properties of the imaged target. Several integral-transform-based moments of the time-resolved data have been utilized in TD-DOT, the most common being the mean time of flight and variance. Recently, it has been shown that Fourier transforming the time-domain data to frequency domain enables utilization of these data at one or several frequencies, producing equally as good estimates as the whole time-domain data. In this work, we present a systematic comparison of the usage of the temporal moments and Fourier transformed data in TD-DOT. Both absolute and difference imaging are evaluated using numerical simulations. The simulations show that utilizing temporal moments and Fourier transformed data in TD-DOT provides good quality reconstructions with a good estimation accuracy. These estimates are improved if more than one data type is used. Furthermore, the simulations show that the frequency-domain computations enable computationally cheaper and straightforward implementation of the inverse solver when compared to the temporal moments.

Time-domain diffuse optical tomography utilizing truncated Fourier series approximation

Diffuse optical tomography (DOT) uses near infrared light for in vivo imaging of spatially varying optical parameters in biological tissues. It is known that time-resolved measurements provide the richest information on soft tissues, among other measurement types in DOT such as steady-state and intensity-modulated measurements. Therefore, several integral-transform-based moments of the time-resolved DOT measurements have been considered to estimate spatially distributed optical parameters. However, the use of such moments can result in low-contrast images and cross-talks between the reconstructed optical parameters, limiting their accuracy. In this work, we propose to utilize a truncated Fourier series approximation in time-resolved DOT. Using this approximation, we obtained optical parameter estimates with accuracy comparable to using whole time-resolved data that uses low computational time and resources. The truncated Fourier series approximation based estimates also displayed good contrast and minimal parameter cross-talk, and the estimates further improved in accuracy when multiple Fourier frequencies were used.

Incorporating reflection boundary conditions in the Neumann series radiative transport equation: application to photon propagation and reconstruction in diffuse optical imaging

We propose a formalism to incorporate boundary conditions in a Neumann-series-based radiative transport equation. The formalism accurately models the reflection of photons at the tissue-external medium interface using Fresnel's equations. The formalism was used to develop a gradient descent-based image reconstruction technique. The proposed methods were implemented for 3D diffuse optical imaging. In computational studies, it was observed that the average root-mean-square error (RMSE) for the output images and the estimated absorption coefficients reduced by 38% and 84%, respectively, when the reflection boundary conditions were incorporated. These results demonstrate the importance of incorporating boundary conditions that model the reflection of photons at the tissue-external medium interface.

A Highly Similar Mathematical Model for Cerebral Blood Flow Velocity in Geriatric Patients with Suspected Cerebrovascular Disease

Cerebral blood flow velocity(CBFV) is an important parameter for study of cerebral hemodynamics. However, a simple and highly similar mathematical model has not yet been established for analyzing CBFV. To alleviate this issue, through TCD examination in 100 geriatric patients with suspected cerebrovascular disease (46 males and 54 females), we established a representative eighth-order Fourier function V_x(t) that simulates the CBFV. The measured TCD waveforms were compared to those derived from V_x(t), an illustrative Kolmogorov-Smirnov test was employed to determine the validity. The results showed that the TCD waves could been reconstructed for patients with different CBFVs by implementing their variable heart rates and the formulated maximum/minimum of V_x(t). Comparisons between derived and measured TCD waveforms suggest that the two waveforms are very similar. The results confirm that CBFV can be well-modeled through an eighth-order Fourier function. This function V_x(t) can be used extensively for a prospective study of cerebral hemodynamics in geriatric patients with suspected cerebrovascular disease.

Full time-resolved diffuse fluorescence tomography accelerated with parallelized Fourier-series truncated diffusion approximation

Of the three measurement schemes established for diffuse fluorescence tomography (DFT), the time-domain scheme is well known to provide the richest information about the distribution of the targeting fluorophore in living tissues. However, the explicit use of the full time-resolved data usually leads to a considerably lengthy time for image reconstruction, limiting its applications to three-dimensional or small-volume imaging. To cope with the adversity, we propose herein a computationally efficient scheme for DFT image reconstruction where the time-dependent photon density is expanded to a Fourier-series and calculated by solving the independent frequency-domain diffusion equations at multiple sampling frequencies with the support of a combined multicore CPU-based coarse-grain and multithread GPU-based fine-grain parallelization strategy. With such a parallelized Fourier-series truncated diffusion approximation, both the time- and frequency-domain inversion procedures are developed and validated for their effectiveness and accuracy using simulative and phantom experiments. The results show that the proposed method can generate reconstructions comparable to the explicit time-domain scheme, with significantly reduced computational time.

Instrumentation in Diffuse Optical Imaging

Diffuse optical imaging is highly versatile and has a very broad range of applications in biology and medicine. It covers diffuse optical tomography, fluorescence diffuse optical tomography, bioluminescence, and a number of other new imaging methods. These methods of diffuse optical imaging have diversified instrument configurations but share the same core physical principle – light propagation in highly diffusive media, i.e., the biological tissue. In this review, the author summarizes the latest development in instrumentation and methodology available to diffuse optical imaging in terms of system architecture, light source, photo-detection, spectral separation, signal modulation, and lastly imaging contrast.