Experimental spectro-angular mapping of light distribution in turbid media

High speed two-photon laser scanning stereomicroscopy for three-dimension tracking multiple particles simultaneously in three-dimension

In this paper, we will describe a video rate two-photon laser scanning stereomicroscopy for imaging-based three-dimensional particle tracking. Using a resonant galvanometer, we have now achieved 30 volumes per second (frame size 512 × 512) in volumetric imaging. Owing to the pulse multiplexing and demultiplexing techniques, the system does not suffer the speed loss for taking two parallax views of a volume. The switching time between left and right views is reduced to several nanoseconds. The extremely fast view switching and high volumetric imaging speed allow us to track fast transport processes of nanoparticles in deep light-scattering media. For instance, in 1%-intralipid solution and fibrillar scaffolds, the tracking penetration depth can be around 400 μm.

An experimental and numerical modelling investigation of the optical properties of Intralipid using deep Raman spectroscopy

In this study, Monte Carlo simulations were created to investigate the distribution of Raman signals in tissue phantoms and to validate the arctk code that was used. The aim was to show our code is capable of replicating experimental results in order to use it to advise similar future studies and to predict the outcomes. The experiment performed to benchmark our code used large volume liquid tissue phantoms to simulate the scattering properties of human tissue. The scattering agent used was Intralipid (IL), of various concentrations, filling a small quartz tank. A thin sample of PTFE was made to act as a distinct layer in the tank; this was our Raman signal source. We studied experimentally, and then reproduced via simulations, the variation in Raman signal strength in a transmission geometry as a function of the optical properties of the scattering agent and the location of the Raman material in the volume. We have also found that a direct linear extrapolation of scattering coefficients between concentrations of Intralipid is an incorrect assumption at lower concentrations when determining the optical properties. By combining experimental and simulation results, we have calculated different estimates of these scattering coefficients. The results of this study give insight into light propagation and Raman transport in scattering media and show how the location of maximum Raman signal varies as the optical properties change. The success of arctk in reproducing observed experimental signal behaviour will allow us in future to inform the development of noninvasive cancer screening applications (such as breast and prostate cancers) in vivo.

Computationally Efficient Forward Operator for Photoacoustic Tomography Based on Coordinate Transformations

Photoacoustic tomography (PAT) is an imaging modality that utilizes the photoacoustic effect. In PAT, a photoacoustic image is computed from measured data by modeling ultrasound propagation in the imaged domain and solving an inverse problem utilizing a discrete forward operator. However, in realistic measurement geometries with several ultrasound transducers and relatively large imaging volume, an explicit formation and use of the forward operator can be computationally prohibitively expensive. In this work, we propose a transformation-based approach for efficient modeling of photoacoustic signals and reconstruction of photoacoustic images. In the approach, the forward operator is constructed for a reference ultrasound transducer and expanded into a general measurement geometry using transformations that map the formulated forward operator in local coordinates to the global coordinates of the measurement geometry. The inverse problem is solved using a Bayesian framework. The approach is evaluated with numerical simulations and experimental data. The results show that the proposed approach produces accurate 3-D photoacoustic images with a significantly reduced computational cost both in memory requirements and time. In the studied cases, depending on the computational factors, such as discretization, over the 30-fold reduction in memory consumption was achieved without a reduction in image quality compared to a conventional approach.

Evaluation of hyperspectral NIRS for quantitative measurements of tissue oxygen saturation by comparison to time-resolved NIRS

Near-infrared spectroscopy (NIRS) is considered ideal for brain monitoring during preterm infancy because it is non-invasive and provides a continuous measure of tissue oxygen saturation (StO2). Hyperspectral NIRS (HS NIRS) is an inexpensive, quantitative modality that can measure tissue optical properties and oxygen saturation (StO2) by differential spectroscopy. In this study, experiments were conducted using newborn piglets to measure StO2 across a range of oxygenation levels from hyperoxia to hypoxia by HS and time-resolved (TR) NIRS for validation. A strong correlation between StO2 measurements from the two techniques was observed (R2 = 0.98, average slope of 1.02 ± 0.28); however, the HS-NIRS estimates were significantly higher than the corresponding TR-NIRS values. These regression results indicate that HS NIRS could become a clinically feasible method for monitoring StO2 in preterm infants.

Determination of optical properties of turbid medium from relative interstitial CW radiance measurements using the incomplete P3 approximation

Interstitial determination of the tissue optical properties is important in biomedicine, especially for interstitial laser therapies. Continuous wave (CW) radiance techniques which examine light from multiple directions have been proposed as minimally invasive methods for determining the optical properties under an interstitial probe arrangement. However, both the fitting algorithm based on the P3 approximation and the analytical method based on the diffusion approximation (DA), which are currently used recovery algorithms, cannot extract the optical properties of tissue with low transport albedos accurately from radiance measurements. In this paper, we proposed an incomplete P3 approximation for the radiance, the P3_in for short, which is the asymptotic part of the solution for the P3 approximation. The relative differences between the P3_in and the P3 were within 0.48% over a wide range of clinically relevant optical properties for measurements at source detector separations (SDS) from 5 mm to 10 mm and angles from 0° to 160°. Based on the P3_in, we developed an analytical method for extracting the optical properties directly using simple expressions constructed from the radiance measurements at only two SDSs and four angles. The developed recovery algorithm was verified by simulated and experimental radiance data. The results show that both the absorption and reduced scattering coefficients were recovered accurately with relative errors within 5.28% and 3.86%, respectively, from the simulated data and with relative errors within 10.82% and 10.67%, respectively, from the experimental data over a wide range of albedos from 0.5 to 0.99. Since the developed P3_in-based radiance technique can obtain the optical properties rapidly from the measurements at only two SDSs and four angles, it is expected to be used for in vivo and in situ determination of the optical properties in online treatment planning during laser therapies.

Optical absorption and scattering properties of bulk porcine muscle phantoms from interstitial radiance measurements in 650-900 nm range

We demonstrated the application of relative radiance-based continuous wave (cw) measurements for recovering absorption and scattering properties (the effective attenuation coefficient, the diffusion coefficient, the absorption coefficient and the reduced scattering coefficient) of bulk porcine muscle phantoms in the 650-900 nm spectral range. Both the side-firing fiber (the detector) and the fiber with a spherical diffuser at the end (the source) were inserted interstitially at predetermined locations in the phantom. The porcine phantoms were prostate-shaped with ∼4 cm in diameter and ∼3 cm thickness and made from porcine loin or tenderloin muscles. The described method was previously validated using the diffusion approximation on simulated and experimental radiance data obtained for homogenous Intralipid-1% liquid phantom. The approach required performing measurements in two locations in the tissue with different distances to the source. Measurements were performed on 21 porcine phantoms. Spectral dependences of the effective attenuation and absorption coefficients for the loin phantom deviated from corresponding dependences for the tenderloin phantom for wavelengths <750 nm. The diffusion constant and the reduced scattering coefficient were very close for both phantom types. To quantify chromophore presence, the plot for the absorption coefficient was matched with a synthetic absorption spectrum constructed from deoxyhemoglobin, oxyhemoglobin and water. The closest match for the porcine loin spectrum was obtained with the following concentrations: 15.5 µM (±30% s.d.) Hb, 21 µM (±30% s.d.) HbO2 and 0.3 (±30% s.d.) fractional volume of water. The tenderloin absorption spectrum was best described by 30 µM Hb (±30% s.d), 19 µM (±30% s.d.) HbO2 and 0.3 (±30% s.d.) fractional volume of water. The higher concentration of Hb in tenderloin was consistent with a dark-red appearance of the tenderloin phantom. The method can be applied to a number of biological tissues and organs for interstitial optical interrogation.

Tagging photons with gold nanoparticles as localized absorbers in optical measurements in turbid media

We analyze a role of a localized inclusion as a probe for spatial distributions of migrating photons in turbid media. We present new experimental data and two-dimensional analysis of radiance detection of a localized absorptive inclusion formed by gold nanoparticles in Intralipid-1% when the target is translated along the line connecting the light source and detector. Data are analyzed using the novel analytical expression for the relative angular photon distribution function for radiance developed by extending the perturbation approach for fluence. Obtained photon maps allow predicting conditions for detectability of inclusions for which proximity to the detector is essential.

Radiance detection of non-scattering inclusions in turbid media

Detection of non-scattering domains (voids) is an area of active research in biomedical optics. To avoid complexities of image reconstruction algorithms and requirements of a priori knowledge of void locations inherent to diffuse optical tomography (DOT), it would be useful to establish specific experimental signatures of voids that would help identify and detect them by other means. To address this, we present a radiance-based spectro-angular mapping approach that identifies void locations in the angular domain and establishes their spectral features. Using water-filled capillaries in scattering Intralipid as a test platform, we demonstrate perturbations in the directional photon density distribution produced by individual voids.

Separation of absorption and scattering properties of turbid media using relative spectrally resolved cw radiance measurements

We present a new method for extracting the effective attenuation coefficient and the diffusion coefficient from relative spectrally resolved cw radiance measurements using the diffusion approximation. The method is validated on both simulated and experimental radiance data sets using Intralipid-1% as a test platform. The effective attenuation coefficient is determined from a simple algebraic expression constructed from a ratio of two radiance measurements at two different source-detector separations and the same 90° angle. The diffusion coefficient is determined from another ratio constructed from two radiance measurements at two angles (0° and 180°) and the same source-detector separation. The conditions of the validity of the method as well as possible practical applications are discussed.