Sub-diffusive scattering parameter maps recovered using wide-field high-frequency structured light imaging

Impact of scattering phase function and polarization on the accuracy of diffuse and sub-diffuse spatial frequency domain imaging

Significance

Although spatial frequency domain imaging (SFDI) has been well characterized under diffuse optical conditions, tissue measurements made outside the diffuse regime can provide new diagnostic information. Before such measurements can become clinically relevant, however, the behavior of sub-diffuse SFDI and its effect on the accuracy of derived tissue parameters must be assessed.

Aim

We aim to characterize the impact that both the assumed scattering phase function (SPF) and the polarization state of the illumination light source have on the accuracy of SFDI-derived optical properties when operating under diffuse or sub-diffuse conditions, respectively.

Approach

Through the use of a set of well-characterized optical phantoms, SFDI accuracy was assessed at four wavelengths (395, 545, 625, and 850 nm) and two different spatial frequencies (0.3 and 1.0    mm - 1 ), which provided a broad range of diffuse and sub-diffuse conditions, using three different SPFs. To determine the effects of polarization, the SFDI accuracy was assessed using both unpolarized and cross-polarized illumination.

Results

It was found that the assumed SPF has a direct and significant impact on the accuracy of the SFDI-derived optical properties, with the best choice of SPF being dictated by the polarization state. As unpolarized SFDI retains the sub-diffuse portion of the signal, optical properties were found to be more accurate when using the full SPF that includes forward and backscattering components. By contrast, cross-polarized SFDI yielded accurate optical properties when using a forward-scattering SPF, matching the behavior of cross-polarization to attenuate the immediate backscattering of sub-diffuse reflectance. Using the correct pairings of SPF and polarization enabled using a reflectance standard, instead of a more subjective phantom, as the reference measurement.

Conclusions

These results provide the foundation for a more thorough understanding of SFDI and enable new applications of this technology in which sub-diffuse conditions dominate (e.g.,μ a ≮ μ s ' ) or high spatial frequencies are required.

Model for the diffuse reflectance in spatial frequency domain imaging

Significance

In spatial frequency domain imaging (SDFI), tissue is illuminated with sinusoidal intensity patterns at different spatial frequencies. For low spatial frequencies, the reflectance is diffuse and a model derived by Cuccia et al. (doi 10.1117/1.3088140) is commonly used to extract optical properties. An improved model resulting in more accurate optical property extraction could lead to improved diagnostic algorithms.

Aim

To develop a model that improves optical property extraction for the diffuse reflectance in SFDI compared to the model of Cuccia et al.

Approach

We derive two analytical models for the diffuse reflectance, starting from the theoretical radial reflectance R ( ρ ) for a pencil-beam illumination under the partial current boundary condition (PCBC) and the extended boundary condition (EBC). We compare both models and the model of Cuccia et al. to Monte Carlo simulations.

Results

The model based on the PCBC resulted in the lowest errors, improving median relative errors compared to the model of Cuccia et al. by 45% for the reflectance, 10% for the reduced scattering coefficient and 64% for the absorption coefficient.

Conclusions

For the diffuse reflectance in SFDI, the model based on the PCBC provides more accurate results than the currently used model by Cuccia et al.

Radiance backscattered by a strongly scattering medium in the high spatial frequency limit

We study the radiative transfer of a spatially modulated plane wave incident on a half-space composed of a uniformly scattering and absorbing medium. For spatial frequencies that are large compared to the scattering coefficient, we find that first-order scattering governs the leading behavior of the radiance backscattered by the medium. The first-order scattering approximation reveals a specific curve on the backscattered hemisphere where the radiance is concentrated. Along this curve, the radiance assumes a particularly simple expression that is directly proportional to the phase function. These results are inherent to the radiative transfer equation at large spatial frequency and do not have a strong dependence on any particular optical property. Consequently, these results provide the means by which spatial frequency domain imaging technologies can directly measure the phase function of a sample. Numerical simulations using the discrete ordinate method along with the source integration interpolation method validate these theoretical findings.

Simple demodulation method for optical property extraction in spatial frequency domain imaging

Different demodulation methods affect the efficiency and accuracy of spatial frequency domain imaging (SFDI). A simple and effective method of sum-to-product identities (STPI) demodulation was proposed in this study. STPI requires one fewer image than conventional three-phase demodulation (TPD) at a spatial frequency. Numerical simulation and phantom experiments were performed. The result proved the feasibility of STPI and showed that STPI combined with subtraction can achieve high-precision demodulation in the low spatial frequency domain. Through extraction of phantom optical properties, STPI had similar accuracy compared with other demodulation methods in extracting optical properties in phantoms. STPI was also used to extract the optical properties of milk, and it had highly consistent results with TPD, which can distinguish milk with different fat content. The demodulation effect of this method in the low spatial frequencies is better than other fast demodulation methods.

Imaging sub-diffuse optical properties of cancerous and normal skin tissue using machine learning-aided spatial frequency domain imaging

SIGNIFICANCE

Sub-diffuse optical properties may serve as useful cancer biomarkers, and wide-field heatmaps of these properties could aid physicians in identifying cancerous tissue. Sub-diffuse spatial frequency domain imaging (sd-SFDI) can reveal such wide-field maps, but the current time cost of experimentally validated methods for rendering these heatmaps precludes this technology from potential real-time applications.

AIM

Our study renders heatmaps of sub-diffuse optical properties from experimental sd-SFDI images in real time and reports these properties for cancerous and normal skin tissue subtypes.

APPROACH

A phase function sampling method was used to simulate sd-SFDI spectra over a wide range of optical properties. A machine learning model trained on these simulations and tested on tissue phantoms was used to render sub-diffuse optical property heatmaps from sd-SFDI images of cancerous and normal skin tissue.

RESULTS

The model accurately rendered heatmaps from experimental sd-SFDI images in real time. In addition, heatmaps of a small number of tissue samples are presented to inform hypotheses on sub-diffuse optical property differences across skin tissue subtypes.

CONCLUSION

These results bring the overall process of sd-SFDI a fundamental step closer to real-time speeds and set a foundation for future real-time medical applications of sd-SFDI such as image guided surgery.

Modeling and Synthesis of Breast Cancer Optical Property Signatures With Generative Models

Is it possible to find deterministic relationships between optical measurements and pathophysiology in an unsupervised manner and based on data alone? Optical property quantification is a rapidly growing biomedical imaging technique for characterizing biological tissues that shows promise in a range of clinical applications, such as intraoperative breast-conserving surgery margin assessment. However, translating tissue optical properties to clinical pathology information is still a cumbersome problem due to, amongst other things, inter- and intrapatient variability, calibration, and ultimately the nonlinear behavior of light in turbid media. These challenges limit the ability of standard statistical methods to generate a simple model of pathology, requiring more advanced algorithms. We present a data-driven, nonlinear model of breast cancer pathology for real-time margin assessment of resected samples using optical properties derived from spatial frequency domain imaging data. A series of deep neural network models are employed to obtain sets of latent embeddings that relate optical data signatures to the underlying tissue pathology in a tractable manner. These self-explanatory models can translate absorption and scattering properties measured from pathology, while also being able to synthesize new data. The method was tested on a total of 70 resected breast tissue samples containing 137 regions of interest, achieving rapid optical property modeling with errors only limited by current semi-empirical models, allowing for mass sample synthesis and providing a systematic understanding of dataset properties, paving the way for deep automated margin assessment algorithms using structured light imaging or, in principle, any other optical imaging technique seeking modeling. Code is available.

Speckle illumination SFDI for projector-free optical property mapping

Spatial frequency domain imaging can map tissue scattering and absorption properties over a wide field of view, making it useful for clinical applications such as wound assessment and surgical guidance. This technique has previously required the projection of fully characterized illumination patterns. Here, we show that random and unknown speckle illumination can be used to sample the modulation transfer function of tissues at known spatial frequencies, allowing the quantitative mapping of optical properties with simple laser diode illumination. We compute low- and high-spatial frequency response parameters from the local power spectral density for each pixel and use a lookup table to accurately estimate absorption and scattering coefficients in tissue phantoms, in vivo human hand, and ex vivo swine esophagus. Because speckle patterns can be generated over a large depth of field and field of view with simple coherent illumination, this approach may enable optical property mapping in new form-factors and applications, including endoscopy.

Active line scan with spatial gating for sub-diffuse reflectance imaging of scatter microtexture

We examine the value of an active line scan with spatial gating for imaging sub-diffuse, wide-field reflectance microtexture. Line scanning combined with spatial gating and linear translation can be used for localized detection of features in the surface layer of a turbid target. The line scan provides broadband spatial frequency modulation, and the spatial gating effectively high-pass filters the reflectance. The major benefit of this approach is that of high dynamic range (70%-90%) signal preservation and high contrast to noise when imaging at high spatial frequencies. Alternative approaches, such as spatial frequency domain imaging, are degraded by low dynamic range in demodulated images, making it nearly impossible to image over a wide field of view at frequencies over 1.5mm-1 using commercial technology. As such, active line scanning with spatial gating presents as an inherently high sensitivity and high dynamic range method of imaging microscopic scattering features in only the surface layer of a turbid medium.

Optical properties of PlatSil SiliGlass tissue-mimicking phantoms

In this work, we revise the preparation procedure and conduct an in depth characterization of optical properties for the recently proposed silicone-based tissue-mimicking optical phantoms in the spectral range from 475 to 925 nm. The optical properties are characterized in terms of refractive index and its temperature dependence, absorption and reduced scattering coefficients and scattering phase function related quantifiers. The scattering phase function and related quantifiers of the optical phantoms are first assessed within the framework of the Mie theory by using the measured refractive index of SiliGlass and size distribution of the hollow silica spherical particles that serve as scatterers. A set of purely absorbing optical phantoms in cuvettes is used to evaluate the linearity of the absorption coefficient with respect to the concentration of black pigment that serves as the absorber. Finally, the optical properties in terms of the absorption and reduced scattering coefficients and the subdiffusive scattering phase function quantifier γ are estimated for a subset of phantoms from spatially resolved reflectance using deep learning aided inverse models. To this end, an optical fiber probe with six linearly arranged optical fibers is used to collect the backscattered light at small and large distances from the source fiber. The underlying light propagation modeling is based on the stochastic Monte Carlo method that accounts for all the details of the optical fiber probe.

Wide-field optical property mapping and structured light imaging of the esophagus with spatial frequency domain imaging

As the incidence of esophageal adenocarcinoma continues to rise, there is a need for improved imaging technologies with contrast to abnormal esophageal tissues. To inform the design of optical technologies that meet this need, we characterize the spatial distribution of the scattering and absorption properties from 471 to 851 nm of eight resected human esophagi tissues using Spatial Frequency Domain Imaging. Histopathology was used to categorize tissue types, including normal, inflammation, fibrotic, ulceration, Barrett's Esophagus and squamous cell carcinoma. Average absorption and reduced scattering coefficients of normal tissues were 0.211 ± 0.051 and 1.20 ± 0.18 mm-1 , respectively at 471 nm, and both values decreased monotonically with increasing wavelength. Fibrotic tissue exhibited at least 68% larger scattering signal across all wavelengths, while squamous cell carcinoma exhibited a 36% decrease in scattering at 471 nm. We additionally image the esophagus with high spatial frequencies up to 0.5 mm-1 and show strong reflectance contrast to tissue treated with radiation. Lastly, we observe that esophageal absorption and scattering values change by an average of 9.4% and 2.7% respectively over a 30 minute duration post-resection. These results may guide system design for the diagnosis, prevention and monitoring of esophageal pathologies.

Structured light imaging for breast-conserving surgery, part II: texture analysis and classification

Subdiffuse spatial frequency domain imaging (sd-SFDI) data of 42 freshly excised, bread-loafed tumor resections from breast-conserving surgery (BCS) were evaluated using texture analysis and a machine learning framework for tissue classification. Resections contained 56 regions of interest (RoIs) determined by expert histopathological analysis. RoIs were coregistered with sd-SFDI data and sampled into ∼4  ×  4  mm2 subimage samples of confirmed and homogeneous histological categories. Sd-SFDI reflectance textures were analyzed using gray-level co-occurrence matrix pixel statistics, image primitives, and power spectral density curve parameters. Texture metrics exhibited statistical significance (p-value  <  0.05) between three benign and three malignant tissue subtypes. Pairs of benign and malignant subtypes underwent texture-based, binary classification with correlation-based feature selection. Classification performance was evaluated using fivefold cross-validation and feature grid searching. Classification using subdiffuse, monochromatic reflectance (illumination spatial frequency of fx  =  1.37  mm  −  1, optical wavelength of λ  =  490  nm) achieved accuracies ranging from 0.55 (95% CI: 0.41 to 0.69) to 0.95 (95% CI: 0.90 to 1.00) depending on the benign–malignant diagnosis pair. Texture analysis of sd-SFDI data maintains the spatial context within images, is free of light transport model assumptions, and may provide an alternative, computationally efficient approach for wide field-of-view (cm2) BCS tumor margin assessment relative to pixel-based optical scatter or color properties alone.

Structured light imaging for breast-conserving surgery, part I: optical scatter and color analysis

Structured light imaging (SLI) with high spatial frequency (HSF) illumination provides a method to amplify native tissue scatter contrast and better differentiate superficial tissues. This was investigated for margin analysis in breast-conserving surgery (BCS) and imaging gross clinical tissues from 70 BCS patients, and the SLI distinguishability was examined for six malignancy subtypes relative to three benign/normal breast tissue subtypes. Optical scattering images recovered were analyzed with five different color space representations of multispectral demodulated reflectance. Excluding rare combinations of invasive lobular carcinoma and fibrocystic disease, SLI was able to classify all subtypes of breast malignancy from surrounding benign tissues (p-value  <  0.05) based on scatter and color parameters. For color analysis, HSF illumination of the sample generated more statistically significant discrimination than regular uniform illumination. Pathological information about lesion subtype from a presurgical biopsy can inform the search for malignancy on the surfaces of specimens during BCS, motivating the focus on pairwise classification analysis. This SLI modality is of particular interest for its potential to differentiate tissue classes across a wide field-of-view (∼100  cm2) and for its ability to acquire images of macroscopic tissues rapidly but with microscopic-level sensitivity to structural and morphological tissue constituents.

Optical sampling depth in the spatial frequency domain

We present a Monte Carlo (MC) method to determine depth-dependent probability distributions of photon visitation and detection for optical reflectance measurements performed in the spatial frequency domain (SFD). These distributions are formed using an MC simulation for radiative transport that utilizes a photon packet weighting procedure consistent with the two-dimensional spatial Fourier transform of the radiative transport equation. This method enables the development of quantitative metrics for SFD optical sampling depth in layered tissue and its dependence on both tissue optical properties and spatial frequency. We validate the computed depth-dependent probability distributions using SFD measurements in a layered phantom system with a highly scattering top layer of variable thickness supported by a highly absorbing base layer. We utilize our method to establish the spatial frequency-dependent optical sampling depth for a number of tissue types and also provide a general tool to determine such depths for tissues of arbitrary optical properties.

Spatial frequency domain imaging in 2019: principles, applications, and perspectives

Spatial frequency domain imaging (SFDI) has witnessed very rapid growth over the last decade, owing to its unique capabilities for imaging optical properties and chromophores over a large field-of-view and in a rapid manner. We provide a comprehensive review of the principles of this imaging method as of 2019, review the modeling of light propagation in this domain, describe acquisition methods, provide an understanding of the various implementations and their practical limitations, and finally review applications that have been published in the literature. Importantly, we also introduce a group effort by several key actors in the field for the dissemination of SFDI, including publications, advice in hardware and implementations, and processing code, all freely available online.

Spatial frequency domain imaging in 2019: principles, applications, and perspectives

Spatial frequency domain imaging (SFDI) has witnessed very rapid growth over the last decade, owing to its unique capabilities for imaging optical properties and chromophores over a large field-of-view and in a rapid manner. We provide a comprehensive review of the principles of this imaging method as of 2019, review the modeling of light propagation in this domain, describe acquisition methods, provide an understanding of the various implementations and their practical limitations, and finally review applications that have been published in the literature. Importantly, we also introduce a group effort by several key actors in the field for the dissemination of SFDI, including publications, advice in hardware and implementations, and processing code, all freely available online.

Light scattering methods for tissue diagnosis

Light scattering has become a common biomedical research tool, enabling diagnostic sensitivity to myriad tissue alterations associated with disease. Light-tissue interactions are particularly attractive for diagnostics due to the variety of contrast mechanisms that can be used, including spectral, angle-resolved, and Fourier-domain detection. Photonic diagnostic tools offer further benefit in that they are non-ionizing, non-invasive, and give real-time feedback. In this review, we summarize recent innovations in light scattering technologies, with a focus on clinical achievements over the previous ten years.

Light scattering measured with spatial frequency domain imaging can predict stromal versus epithelial proportions in surgically resected breast tissue

This study aims to determine if light scatter parameters measured with spatial frequency domain imaging (SFDI) can accurately predict stromal, epithelial, and adipose fractions in freshly resected, unstained human breast specimens. An explicit model was developed to predict stromal, epithelial, and adipose fractions as a function of light scattering parameters, which was validated against a quantitative analysis of digitized histology slides for N  =  31 specimens using leave-one-out cross-fold validation. Specimen mean stromal, epithelial, and adipose volume fractions predicted from light scattering parameters strongly correlated with those calculated from digitized histology slides (r  =  0.90, 0.77, and 0.91, respectively, p-value <1  ×  10  -  6). Additionally, the ratio of predicted epithelium to stroma classified malignant specimens with a sensitivity and specificity of 90% and 81%, respectively, and also classified all pixels in malignant lesions with 63% and 79%, at a threshold of 1. All specimens and pixels were classified as malignant, benign, or fat with 84% and 75% accuracy, respectively. These findings demonstrate how light scattering parameters acquired with SFDI can be used to accurately predict and spatially map stromal, epithelial, and adipose proportions in fresh unstained, human breast tissue, and suggest that these estimations could provide diagnostic value.

Review of structured light in diffuse optical imaging

Diffuse optical imaging probes deep living tissue enabling structural, functional, metabolic, and molecular imaging. Recently, due to the availability of spatial light modulators, wide-field quantitative diffuse optical techniques have been implemented, which benefit greatly from structured light methodologies. Such implementations facilitate the quantification and characterization of depth-resolved optical and physiological properties of thick and deep tissue at fast acquisition speeds. We summarize the current state of work and applications in the three main techniques leveraging structured light: spatial frequency-domain imaging, optical tomography, and single-pixel imaging. The theory, measurement, and analysis of spatial frequency-domain imaging are described. Then, advanced theories, processing, and imaging systems are summarized. Preclinical and clinical applications on physiological measurements for guidance and diagnosis are summarized. General theory and method development of tomographic approaches as well as applications including fluorescence molecular tomography are introduced. Lastly, recent developments of single-pixel imaging methodologies and applications are reviewed.

Spatial frequency domain imaging using an analytical model for separation of surface and volume scattering

A method to correct for surface scattering in spatial frequency domain imaging (SFDI) is presented. The use of a modified analytical solution of the radiative transfer equation allows calculation of the reflectance and the phase of a rough semi-infinite geometry so that both spatial frequency domain reflectance and phase can be applied for precise retrieval of the bulk optical properties and the surface scattering. For validation of the method, phantoms with different surface roughness were produced. Contrarily, with the modified theory, it was possible to dramatically reduce systematic errors due to surface scattering. The evaluation of these measurements with the state-of-the-art theory and measuring modality, i.e., using crossed linear polarizers, reveals large errors in the determined optical properties, depending on the surface roughness, of up to ≈100  %  . These results were confirmed with SFDI measurements on a phantom that has a structured rough surface.

Sampling depth of a diffuse reflectance spectroscopy probe for in-vivo physiological quantification of murine subcutaneous tumor allografts

Diffuse reflectance spectroscopy (DRS) is a probe-based spectral biopsy technique used in cancer studies to quantify tissue reduced scattering (μs') and absorption (μa) coefficients and vary in source-detector separation (SDS) to fine-tune sampling depth. In subcutaneous murine tumor allografts or xenografts, a key design requirement is ensuring that the source light interrogates past the skin layer into the tumor without significantly sacrificing signal-to-noise ratio (target of ≥15  dB). To resolve this requirement, a DRS probe was designed with four SDSs (0.75, 2.00, 3.00, and 4.00 mm) to interrogate increasing tissue volumes between 450 and 900 nm. The goal was to quantify percent errors in extracting μa and μs', and to quantify sampling depth into subcutaneous Balb/c-CT26 colon tumor allografts. Using an optical phantom-based experimental method, lookup-tables were constructed relating μa,μs', diffuse reflectance, and sampling depth. Percent errors were <10  %   and 5% for extracting μa and μs', respectively, for all SDSs. Sampling depth reached up to 1.6 mm at the first Q-band of hemoglobin at 542 nm, the key spectral region for quantifying tissue oxyhemoglobin concentration. This work shows that the DRS probe can accurately extract optical properties and the resultant physiological parameters such as total hemoglobin concentration and tissue oxygen saturation, from sufficient depth within subcutaneous Balb/c-CT26 colon tumor allografts. Methods described here can be generalized for other murine tumor models. Future work will explore the feasibility of the DRS in quantifying volumetric tumor perfusion in response to anticancer therapies.

Quantitative diagnosis of tissue microstructure with wide-field high spatial frequency domain imaging

Non-contact and minimally invasive endoscopic optical imaging is an invaluable diagnostic tool for tissue examination and cancer screening. The point sampling techniques with high sensitivity to the tissue microenvironment are time consuming and often not affordable in clinics. There is a major clinical need for a large field-of-view (FOV) rapid screening method to highlight subtle tissue microstructural alterations. To address this unmet need, we have developed High Spatial Frequency Domain Imaging (HSFDI)-a non-contact imaging modality that spatially maps the tissue microscopic scattering structures over a large field of view (>1cm2). Based on an analytical reflectance model of sub-diffusive light from forward-peaked highly scattering media, HSFDI quantifies the spatially-resolved parameters of the light scattering phase function (i.e., the backscattering probability and the light spreading length) from the reflectance of structured light modulated at high spatial frequencies. Enhanced signal to noise ratio (SNR) is achieved at even ultra-high modulation frequencies with single snapshot multiple frequency demodulation (SSMD). The variations in tissue microstructures, including the strength of the background (pudding) refractive index fluctuation and the prominent scattering structure (plum) morphology, can then be inferred. After validation with optical phantoms, measurements of fresh ex vivo tissue samples revealed significant contrast and differentiation of the phase function parameters between different types and disease states (normal, inflammatory, and cancerous) of tissue whereas tissue absorption and reduced scattering coefficients only show modest changes. HSFDI may provide wide-field images of microscopic structural biomarkers unobtainable with either diffuse light imaging or point-based optical sampling. Potential clinical applications include the rapid screening of excised tissue and the noninvasive examination of suspicious lesions during operation.

Efficient estimation of subdiffusive optical parameters in real time from spatially resolved reflectance by artificial neural networks

Subdiffusive reflectance captured at short source-detector separations provides increased sensitivity to the scattering phase function and hence allows superficial probing of the tissue ultrastructure. Consequently, estimation of subdiffusive optical parameters has been the subject of many recent studies focusing on lookup-table-based (LUT) inverse models. Since an adequate description of the subdiffusive reflectance requires additional scattering phase function related optical parameters, the LUT inverse models, which grow exponentially with the number of estimated parameters, become excessively large and computationally inefficient. Herein, we propose, to the best of our knowledge, the first artificial-neural-network-based inverse Monte Carlo model that overcomes the limitations of the LUT inverse models and thus allows efficient real-time estimation of optical parameters from subdiffusive spatially resolved reflectance. The proposed inverse model retains the accuracy, is about four orders of magnitude faster than the LUT inverse models, grows only linearly with the number of estimated optical parameters, and can be easily extended to estimate additional optical parameters.

Spectral refractive index assessment of turbid samples by combining spatial frequency near-infrared spectroscopy with Kramers-Kronig analysis

A practical algorithm for estimating the wavelength-dependent refractive index (RI) of a turbid sample in the spatial frequency domain with the aid of Kramers-Kronig (KK) relations is presented. In it, phase-shifted sinusoidal patterns (structured illumination) are serially projected at a high spatial frequency onto the sample surface (mouse scalp) at different near-infrared wavelengths while a camera mounted normally to the sample surface captures the reflected diffuse light. In the offline analysis pipeline, recorded images at each wavelength are converted to spatial absorption maps by logarithmic function, and once the absorption coefficient information is obtained, the imaginary part (k) of the complex RI (CRI), based on Maxell's equations, can be calculated. Using the data represented by k, the real part of the CRI (n) is then resolved by KK analysis. The wavelength dependence of n  (  λ  )   is then fitted separately using four standard dispersion models: Cornu, Cauchy, Conrady, and Sellmeier. In addition, three-dimensional surface-profile distribution of n is provided based on phase profilometry principles and a phase-unwrapping-based phase-derivative-variance algorithm. Experimental results demonstrate the capability of the proposed idea for sample's determination of a biological sample's RI value.

Wide-field color imaging of scatter-based tissue contrast using both high spatial frequency illumination and cross-polarization gating

This study characterizes the scatter-specific tissue contrast that can be obtained by high spatial frequency (HSF) domain imaging and cross-polarization (CP) imaging, using a standard color imaging system, and how combining them may be beneficial. Both HSF and CP approaches are known to modulate the sensitivity of epi-illumination reflectance images between diffuse multiply scattered and superficially backscattered photons, providing enhanced contrast from microstructure and composition than what is achieved by standard wide-field imaging. Measurements in tissue-simulating optical phantoms show that CP imaging returns localized assessments of both scattering and absorption effects, while HSF has uniquely specific sensitivity to scatter-only contrast, with a strong suppression of visible contrast from blood. The combination of CP and HSF imaging provided an expanded sensitivity to scatter compared with CP imaging, while rejecting specular reflections detected by HSF imaging. ex vivo imaging of an atlas of dissected rodent organs/tissues demonstrated the scatter-based contrast achieved with HSF, CP and HSF-CP imaging, with the white light spectral signal returned by each approach translated to a color image for intuitive encoding of scatter-based contrast within images of tissue. The results suggest that visible CP-HSF imaging could have the potential to aid diagnostic imaging of lesions in skin or mucosal tissues and organs, where just CP is currently the standard practice imaging modality.

Calibration and analysis of a multimodal micro-CT and structured light imaging system for the evaluation of excised breast tissue

A multimodal micro-computed tomography (CT) and multi-spectral structured light imaging (SLI) system is introduced and systematically analyzed to test its feasibility to aid in margin delineation during breast conserving surgery (BCS). Phantom analysis of the micro-CT yielded a signal-to-noise ratio of 34, a contrast of 1.64, and a minimum detectable resolution of 240 μm for a 1.2 min scan. The SLI system, spanning wavelengths 490 nm to 800 nm and spatial frequencies up to 1.37 [Formula: see text], was evaluated with aqueous tissue simulating phantoms having variations in particle size distribution, scatter density, and blood volume fraction. The reduced scattering coefficient, [Formula: see text] and phase function parameter, γ, were accurately recovered over all wavelengths independent of blood volume fractions from 0% to 4%, assuming a flat sample geometry perpendicular to the imaging plane. The resolution of the optical system was tested with a step phantom, from which the modulation transfer function was calculated yielding a maximum resolution of 3.78 cycles per mm. The three dimensional spatial co-registration between the CT and optical imaging space was tested and shown to be accurate within 0.7 mm. A freshly resected breast specimen, with lobular carcinoma, fibrocystic disease, and adipose, was imaged with the system. The micro-CT provided visualization of the tumor mass and its spiculations, and SLI yielded superficial quantification of light scattering parameters for the malignant and benign tissue types. These results appear to be the first demonstration of SLI combined with standard medical tomography for imaging excised tumor specimens. While further investigations are needed to determine and test the spectral, spatial, and CT features required to classify tissue, this study demonstrates the ability of multimodal CT/SLI to quantify, visualize, and spatially navigate breast tumor specimens, which could potentially aid in the assessment of tumor margin status during BCS.

Wide-field fluorescence tomography with composited epi-illumination of multi-frequency sinusoidal patterns

We present a spatial-frequency domain (SFD) fluorescence tomography (FT) for acquiring three-dimensional fluorophore distribution in turbid media. The approach uses a composited epi-illumination of multi-frequency sinusoidal patterns on a sample of semi-infinite geometry and demodulates the measured data with a generalized phase shifting scheme to calculate the modulation transfer function (MTF) at each frequency. This method results in a significantly reduced number of the optical field measurements, as compared to those with separate illumination of single-frequency sinusoidal patterns, and, thereby, achieves a fast data acquisition that is desired for a dynamic imaging application. Fluorescence yield images are recovered with the normalized Born formulated inversion of the diffusion model by simultaneously using the multi-frequency MTFs. Simulative and experimental reconstructions are performed in comparison with the single-frequency scheme to validate the proposed algorithm. The results suggest that adoption of the multi-frequency strategy to the SFD-FT can substantially improve the reconstruction quality, as well as its imaging resolution and quantitative accuracy.

Determination of the complex refractive index segments of turbid sample with multispectral spatially modulated structured light and models approximation

Spectral data enabling the derivation of a biological tissue sample's complex refractive index (CRI) can provide a range of valuable information in the clinical and research contexts. Specifically, changes in the CRI reflect alterations in tissue morphology and chemical composition, enabling its use as an optical marker during diagnosis and treatment. In the present work, we report a method for estimating the real and imaginary parts of the CRI of a biological sample using Kramers-Kronig (KK) relations in the spatial frequency domain. In this method, phase-shifted sinusoidal patterns at single high spatial frequency are serially projected onto the sample surface at different near-infrared wavelengths while a camera mounted normal to the sample surface acquires the reflected diffuse light. In the offline analysis pipeline, recorded images at each wavelength are converted to spatial phase maps using KK analysis and are then calibrated against phase-models derived from diffusion approximation. The amplitude of the reflected light, together with phase data, is then introduced into Fresnel equations to resolve both real and imaginary segments of the CRI at each wavelength. The technique was validated in tissue-mimicking phantoms with known optical parameters and in mouse models of ischemic injury and heat stress. Experimental data obtained indicate variations in the CRI among brain tissue suffering from injury. CRI fluctuations correlated with alterations in the scattering and absorption coefficients of the injured tissue are demonstrated. This technique for deriving dynamic changes in the CRI of tissue may be further developed as a clinical diagnostic tool and for biomedical research applications. To the best of our knowledge, this is the first report of the estimation of the spectral CRI of a mouse head following injury obtained in the spatial frequency domain.

Reference-free determination of tissue absorption coefficient by modulation transfer function characterization in spatial frequency domain

Background

Spatial frequency domain (SFD) measurement allows rapid and non-contact wide-field imaging of the tissue optical properties, thus has become a potential tool for assessing physiological parameters and therapeutic responses during photodynamic therapy of skin diseases. The conventional SFD measurement requires a reference measurement within the same experimental scenario as that for a test one to calibrate mismatch between the real measurements and the model predictions. Due to the individual physical and geometrical differences among different tissues, organs and patients, an ideal reference measurement might be unavailable in clinical trials. To address this problem, we present a reference-free SFD determination of absorption coefficient that is based on the modulation transfer function (MTF) characterization.

Methods

Instead of the absolute amplitude that is used in the conventional SFD approaches, we herein employ the MTF to characterize the propagation of the modulated lights in tissues. With such a dimensionless relative quantity, the measurements can be naturally corresponded to the model predictions without calibrating the illumination intensity. By constructing a three-dimensional database that portrays the MTF as a function of the optical properties (both the absorption coefficient μ_a and the reduced scattering coefficient \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu^{\prime}_{s}$$\end{document}μs′) and the spatial frequency, a look-up table approach or a least-square curve-fitting method is readily applied to recover the absorption coefficient from a single frequency or multiple frequencies, respectively.

Results

Simulation studies have verified the feasibility of the proposed reference-free method and evaluated its accuracy in the absorption recovery. Experimental validations have been performed on homogeneous tissue-mimicking phantoms with μ_a ranging from 0.01 to 0.07 mm⁻¹ and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu^{\prime}_{s}$$\end{document}μs′ = 1.0 or 2.0 mm⁻¹. The results have shown maximum errors of 4.86 and 7% for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu^{\prime}_{s}$$\end{document}μs′ = 1.0 mm⁻¹ and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu^{\prime}_{s}$$\end{document}μs′ = 2.0 mm⁻¹, respectively. We have also presented quantitative ex vivo imaging of human lung cancer in a subcutaneous xenograft mouse model for further validation, and observed high absorption contrast in the tumor region.

Conclusions

The proposed method can be applied to the rapid and accurate determination of the absorption coefficient, and better yet, in a reference-free way. We believe this reference-free strategy will facilitate the clinical translation of the SFD measurement to achieve enhanced intraoperative hemodynamic monitoring and personalized treatment planning in photodynamic therapy.

Microscope objective based 4π spectroscopic tissue scattering goniometry

The measurement of optical scattering as a function of angle, goniometry, can provide a wealth of information about tissue. The goniometry technique described here measures the intensity profile at the pupil planes of two microscope objectives with a scattering sample between them. The maximum observable scattering angle is extended by employing off-axis illumination. This configuration permits several advantages including: i) rapid measurement of scattering into 4π sr to characterize the entire scattering phase function in isotropic tissue, ii) sensitivity to axially asymmetric scattering from anisotropic fibrous tissue, iii) selective interrogation of small regions within spatially inhomogenous tissue, iv) concurrent measurement of scattering coefficient μ_s , and v) measurement of wavelength dependent scattering properties via spectrally tunable source. The instrument is validated by comparing measurements of microsphere suspensions to the Mie scattering solution. Instrument capabilities are demonstrated with samples of rat brain and mouse eye tissues.

Finite element simulation of light transfer in turbid media under structured illumination

The spatial-frequency domain (SFD) imaging technique allows us to estimate the optical properties of biological tissues in a wide field of view. The technique is, however, prone to error in measurement because the two crucial assumptions used for deriving the analytical solution to the diffusion approximation cannot be met perfectly in practical applications. This research mainly focused on modeling light transfer in turbid media under the normal incidence of structured illumination using the finite element method (FEM). Finite element simulations were performed for 50 simulation samples with different combinations of optical absorption and scattering coefficients for varying spatial frequencies, and the results were then compared with the analytical method and Monte Carlo simulation. Relationships between diffuse reflectance and dimensionless absorption and dimensionless scattering coefficients were investigated. The results indicated that the FEM provided reasonable results for diffuse reflectance, compared with the analytical method. Both the FEM and the analytical method overestimated the reflectance for μ_tr/f_x values of greater than 2 and underestimated the reflectance for μ_tr/f_x values of smaller than 2. Larger values of μs'/μ_a yielded better diffuse reflectance estimations than did those of smaller than 10. The reflectance increased nonlinearly with the dimensionless scattering, whereas the reflectance decreased linearly with the dimensionless absorption. It was also observed that diffuse reflectance was relatively stable and insensitive to μs' when the dimensionless scattering was larger than 50. Overall results demonstrate that the FEM is effective for modeling light transfer in turbid media and can be used to explore the effects of crucial parameters for the SFD imaging technique.

Preclinical evaluation of spatial frequency domain-enabled wide-field quantitative imaging for enhanced glioma resection

5-Aminolevelunic acid-induced protoporphyrin IX (PpIX) fluorescence-guided resection (FGR) enables maximum safe resection of glioma by providing real-time tumor contrast. However, the subjective visual assessment and the variable intrinsic optical attenuation of tissue limit this technique to reliably delineating only high-grade tumors that display strong fluorescence. We have previously shown, using a fiber-optic probe, that quantitative assessment using noninvasive point spectroscopic measurements of the absolute PpIX concentration in tissue further improves the accuracy of FGR, extending it to surgically curable low-grade glioma. More recently, we have shown that implementing spatial frequency domain imaging with a fluorescent-light transport model enables recovery of two-dimensional images of [PpIX], alleviating the need for time-consuming point sampling of the brain surface. We present first results of this technique modified for <italic<in vivo</italic< imaging on an RG2 rat brain tumor model. Despite the moderate errors in retrieving the absorption and reduced scattering coefficients in the subdiffusive regime of 14% and 19%, respectively, the recovered [PpIX] maps agree within 10% of the point [PpIX] values measured by the fiber-optic probe, validating its potential as an extension or an alternative to point sampling during glioma resection.

Determination of three optical properties from subdiffusive spatially resolved reflectance calculations

We report a theoretical study on the determination of three optical properties from spatially resolved reflectance calculations. In particular, the reduced scattering coefficient μs′, the absorption coefficient μa, and the recently defined phase function parameter σ are identified. The solution of the inverse problem is based on the principal component analysis of a large set of reflectance profiles that were calculated using an analytical solution of the radiative transfer equation. Different phase function types were studied to test the method in the range of 0.63  mm−1≤μs′≤4.2  mm−1 and 0.002  mm−1≤μa≤0.1  mm−1. For curves impaired with noise, we were able to reconstruct μs′ and μa with relative median errors of 2.5% and 12%, respectively, and σ with an absolute median error of 0.04.

Modeling subdiffusive light scattering by incorporating the tissue phase function and detector numerical aperture

To detect small-scale changes in tissue with optical techniques, small sampling volumes and, therefore, short source–detector separations are required. In this case, reflectance measurements are not adequately described by the diffusion approximation. Previous studies related subdiffusive reflectance to ? or ? , which parameterize the phase function. Recently, it was demonstrated that ? predicts subdiffusive reflectance better than ? , and that ? becomes less predictive for lower numerical apertures (NAs). We derive and evaluate the parameter R p NA , which incorporates the NA of the detector and the integral of the phase function over the NA in the backward and forward directions. Monte Carlo simulations are performed for overlapping source/detector geometries for a range of phase functions, reduced scattering coefficients, NAs, and source/detector diameters. R p NA improves prediction of the measured reflectance compared to ? and ? . It is, therefore, expected that R p NA will improve derivation of optical properties from subdiffusive measurements.

Multimodal optical setup based on spectrometer and cameras combination for biological tissue characterization with spatially modulated illumination

The application of optical techniques as tools for biomedical research has generated substantial interest for the ability of such methodologies to simultaneously measure biochemical and morphological parameters of tissue. Ongoing optimization of optical techniques may introduce such tools as alternative or complementary to conventional methodologies. The common approach shared by current optical techniques lies in the independent acquisition of tissue’s optical properties (i.e., absorption and reduced scattering coefficients) from reflected or transmitted light. Such optical parameters, in turn, provide detailed information regarding both the concentrations of clinically relevant chromophores and macroscopic structural variations in tissue. We couple a noncontact optical setup with a simple analysis algorithm to obtain absorption and scattering coefficients of biological samples under test. Technically, a portable picoprojector projects serial sinusoidal patterns at low and high spatial frequencies, while a spectrometer and two independent CCD cameras simultaneously acquire the reflected diffuse light through a single spectrometer and two separate CCD cameras having different bandpass filters at nonisosbestic and isosbestic wavelengths in front of each. This configuration fills the gaps in each other’s capabilities for acquiring optical properties of tissue at high spectral and spatial resolution. Experiments were performed on both tissue-mimicking phantoms as well as hands of healthy human volunteers to quantify their optical properties as proof of concept for the present technique. In a separate experiment, we derived the optical properties of the hand skin from the measured diffuse reflectance, based on a recently developed camera model. Additionally, oxygen saturation levels of tissue measured by the system were found to agree well with reference values. Taken together, the present results demonstrate the potential of this integrated setup for diagnostic and research applications.

Lookup table-based sampling of the phase function for Monte Carlo simulations of light propagation in turbid media

Analytical expressions for sampling the scattering angle from a phase function in Monte Carlo simulations of light propagation are available only for a limited number of phase functions. Consequently, numerical sampling methods based on tabulated values are often required instead. By using Monte Carlo simulated reflectance, we compare two existing and propose an improved numerical sampling method and show that both the number of the tabulated values and the numerical sampling method significantly influence the accuracy of the simulated reflectance. The provided results and guidelines should serve as a good starting point for conducting computationally efficient Monte Carlo simulations with numerical phase function sampling.

Monochromatic subdiffusive spatial frequency domain imaging provides in-situ sensitivity to intratumoral morphological heterogeneity in a murine model

For the first time, spatially resolved quantitative metrics of light scattering recovered with sub-diffusive spatial frequency domain imaging (sd-SFDI) are shown to be sensitive to changes in intratumoral morphology and viability by direct comparison to histopathological analysis. Two freshly excised subcutaneous murine tumor cross-sections were measured with sd-SFDI, and recovered optical scatter parameter maps were co-registered to whole mount histology. Unique clustering of the optical scatter parameters μs' vs. γ (i.e. diffuse scattering vs. relative backscattering) evaluated at a single wavelength showed complete separation between regions of viable tumor, aggresive tumor with stromal growth, varying levels of necrotic tumor, and also peritumor muscle. The results suggest that with further technical development, sd-SFDI may represent a non-destructive screening tool for analysis of excised tissue or a non-invasive approach to investigate suspicious lesions without the need for exogenous labels or spectrally resolved imaging.

Miniature spectral imaging device for wide-field quantitative functional imaging of the morphological landscape of breast tumor margins

We have developed a portable, breast margin assessment probe leveraging diffuse optical spectroscopy to quantify the morphological landscape of breast tumor margins during breast conserving surgery. The approach presented here leverages a custom-made 16-channel annular photodiode imaging array (arranged in a 4 × 4 grid), a raster-scanning imaging platform with precision pressure control, and compressive sensing with an optimized set of eight wavelengths in the visible spectral range. A scalable Monte-Carlo-based inverse model is used to generate optical property [ ? s ? ( ? ) and ? a ( ? ) ] measures for each of the 16 simultaneously captured diffuse reflectance spectra. Subpixel sampling (0.75 mm) is achieved through incremental x , y raster scanning of the imaging probe, providing detailed optical parameter maps of breast margins over a 2 × 2 ?? cm 2 area in ? 9 ?? min . The morphological landscape of a tumor margin is characterized using optical surrogates for the fat to fibroglandular content ratio, which has demonstrated diagnostic utility in delineating tissue subtypes in the breast.

Deep imaging of absorption and scattering features by multispectral multiple scattering low coherence interferometry

We have developed frequency domain multispectral multiple scattering low coherence interferometry (ms2/LCI) for deep imaging of absorption and scattering contrast. Using tissue-mimicking phantoms that match the full scattering phase function of human dermal tissue, we demonstrate that ms2/LCI can provide a signal/noise ratio (SNR) improvement of 15.4 dB over conventional OCT at an imaging depth of 1 mm. The enhanced SNR and penetration depth provided by ms2/LCI could be leveraged for a variety of clinical applications including the assessment of burn injuries where current clinical classification of severity only provides limited accuracy. The utility of the approach was demonstrated by imaging a tissue phantom simulating a partial-thickness burn revealing good spectroscopic contrast between healthy and injured tissue regions deep below the sample surface. Finally, healthy rat skin was imaged in vivo with both a commercial OCT instrument and our custom ms2/LCI system. The results demonstrate that ms2/LCI is capable of obtaining spectroscopic information far beyond the penetration depth provided by conventional OCT.

Estimation of optical properties by spatially resolved reflectance spectroscopy in the subdiffusive regime

We propose and objectively evaluate an inverse Monte Carlo model for estimation of absorption and reduced scattering coefficients and similarity parameter ? from spatially resolved reflectance (SRR) profiles in the subdiffusive regime. The similarity parameter ? carries additional information on the phase function that governs the angular properties of scattering in turbid media. The SRR profiles at five source-detector separations were acquired with an optical fiber probe. The inverse Monte Carlo model was based on a cost function that enabled robust estimation of optical properties from a few SRR measurements without a priori knowledge about spectral dependencies of the optical properties. Validation of the inverse Monte Carlo model was performed on synthetic datasets and measured SRR profiles of turbid phantoms comprising molecular dye and polystyrene microspheres. We observed that the additional similarity parameter ? substantially reduced the reflectance variability arising from the phase function properties and significantly improved the accuracy of the inverse Monte Carlo model. However, the observed improvement of the extended inverse Monte Carlo model was limited to reduced scattering coefficients exceeding ?15??cm?1, where the relative root-mean-square errors of the estimated optical properties were well within 10%.

Wide-field quantitative imaging of tissue microstructure using sub-diffuse spatial frequency domain imaging

Localized measurements of scattering in biological tissue provide sensitivity to microstructural morphology but have limited utility to wide-field applications, such as surgical guidance. This study introduces sub-diffusive spatial frequency domain imaging (sd-SFDI), which uses high spatial frequency illumination to achieve wide-field sampling of localized reflectances. Model-based inversion recovers macroscopic variations in the reduced scattering coefficient [Formula: see text] and the phase function backscatter parameter (γ). Measurements in optical phantoms show quantitative imaging of user-tuned phase-function-based contrast with accurate decoupling of parameters that define both the density and the size-scale distribution of scatterers. Measurements of fresh ex vivo breast tissue samples revealed, for the first time, unique clustering of sub-diffusive scattering properties for different tissue types. The results support that sd-SFDI provides maps of microscopic structural biomarkers that cannot be obtained with diffuse wide-field imaging and characterizes spatial variations not resolved by point-based optical sampling.

Low-cost tissue simulating phantoms with adjustable wavelength-dependent scattering properties in the visible and infrared ranges

We present a method for low-cost fabrication of polydimethylsiloxane (PDMS) tissue simulating phantoms with tunable scattering spectra, spanning visible, and near-infrared regimes. These phantoms use optical polishing agents (aluminum oxide powders) at various grit sizes to approximate in vivo tissue scattering particles across multiple size distributions (range: 17 to 3  μm). This class of tunable scattering phantoms is used to mimic distinct changes in wavelength-dependent scattering properties observed in tissue pathologies such as partial thickness burns. Described by a power-law dependence on wavelength, the scattering magnitude of these phantoms scale linearly with particle concentration over a physiologic range [μs′=(0.5 to 2.0  mm−1)] whereas the scattering spectra, specific to each particle size distribution, correlate to distinct exponential coefficients (range: 0.007 to 0.32). Aluminum oxide powders used in this investigation did not detectably contribute to the absorption properties of these phantoms. The optical properties of these phantoms are verified through inverse adding-doubling methods and the tolerances of this fabrication method are discussed.

Color structured light imaging of skin

We illustrate wide-field imaging of skin using a structured light (SL) approach that highlights the contrast from superficial tissue scattering. Setting the spatial frequency of the SL in a regime that limits the penetration depth effectively gates the image for photons that originate from the skin surface. Further, rendering the SL images in a color format provides an intuitive format for viewing skin pathologies. We demonstrate this approach in skin pathologies using a custom-built handheld SL imaging system.

Quantifying phase function influence in subdiffusively backscattered light

Light backscattering at short source-detector separations is considerably influenced by the scattering phase function of a turbid medium. We seek to more precisely relate a medium's subdiffusive backscattering to the angular scattering characteristics of its microstructure. First, we demonstrate the inability of the scattering asymmetry g1 = < cos θ > to predict phase function influence on backscattering and reveal ambiguities related to the established phase function parameter γ. Through the use of high-order similarity relations, we introduce a new parameter that more accurately relates a scattering phase function to its subdiffusive backscattering intensity. Using extensive analytical forward calculations based on solutions to the radiative transfer equation in the spatial domain and spatial frequency domain, we demonstrate the superiority of our empirically derived quantifier σ over the established parameter γ.

Image-derived arterial input function for quantitative fluorescence imaging of receptor-drug binding in vivo

Receptor concentration imaging (RCI) with targeted-untargeted optical dye pairs has enabled in vivo immunohistochemistry analysis in preclinical subcutaneous tumors. Successful application of RCI to fluorescence guided resection (FGR), so that quantitative molecular imaging of tumor-specific receptors could be performed in situ, would have a high impact. However, assumptions of pharmacokinetics, permeability and retention, as well as the lack of a suitable reference region limit the potential for RCI in human neurosurgery. In this study, an arterial input graphic analysis (AIGA) method is presented which is enabled by independent component analysis (ICA). The percent difference in arterial concentration between the image-derived arterial input function (AIFICA ) and that obtained by an invasive method (ICACAR ) was 2.0 ± 2.7% during the first hour of circulation of a targeted-untargeted dye pair in mice. Estimates of distribution volume and receptor concentration in tumor bearing mice (n = 5) recovered using the AIGA technique did not differ significantly from values obtained using invasive AIF measurements (p = 0.12). The AIGA method, enabled by the subject-specific AIFICA , was also applied in a rat orthotopic model of U-251 glioblastoma to obtain the first reported receptor concentration and distribution volume maps during open craniotomy.

Fiber-bundle microendoscopy with sub-diffuse reflectance spectroscopy and intensity mapping for multimodal optical biopsy of stratified epithelium

Early detection of structural or functional changes in dysplastic epithelia may be crucial for improving long-term patient care. Recent work has explored myriad non-invasive or minimally invasive "optical biopsy" techniques for diagnosing early dysplasia, such as high-resolution microendoscopy, a method to resolve sub-cellular features of apical epithelia, as well as broadband sub-diffuse reflectance spectroscopy, a method that evaluates bulk health of a small volume of tissue. We present a multimodal fiber-based microendoscopy technique that combines high-resolution microendoscopy, broadband (450-750 nm) sub-diffuse reflectance spectroscopy (sDRS) at two discrete source-detector separations (374 and 730 μm), and sub-diffuse reflectance intensity mapping (sDRIM) using a 635 nm laser. Spatial resolution, magnification, field-of-view, and sampling frequency were determined. Additionally, the ability of the sDRS modality to extract optical properties over a range of depths is reported. Following this, proof-of-concept experiments were performed on tissue-simulating phantoms made with poly(dimethysiloxane) as a substrate material with cultured MDA-MB-468 cells. Then, all modalities were demonstrated on a human melanocytic nevus from a healthy volunteer and on resected colonic tissue from a murine model. Qualitative in vivo image data is correlated with reduced scattering and absorption coefficients.

Multifrequency synthesis and extraction using square wave projection patterns for quantitative tissue imaging

We present a method for spatial frequency domain data acquisition utilizing a multifrequency synthesis and extraction (MSE) method and binary square wave projection patterns. By illuminating a sample with square wave patterns, multiple spatial frequency components are simultaneously attenuated and can be extracted to determine optical property and depth information. Additionally, binary patterns are projected faster than sinusoids typically used in spatial frequency domain imaging (SFDI), allowing for short (millisecond or less) camera exposure times, and data acquisition speeds an order of magnitude or more greater than conventional SFDI. In cases where sensitivity to superficial layers or scattering is important, the fundamental component from higher frequency square wave patterns can be used. When probing deeper layers, the fundamental and harmonic components from lower frequency square wave patterns can be used. We compared optical property and depth penetration results extracted using square waves to those obtained using sinusoidal patterns on an in vivo human forearm and absorbing tube phantom, respectively. Absorption and reduced scattering coefficient values agree with conventional SFDI to within 1% using both high frequency (fundamental) and low frequency (fundamental and harmonic) spatial frequencies. Depth penetration reflectance values also agree to within 1% of conventional SFDI.

Detecting structural information of scatterers using spatial frequency domain imaging

We demonstrate optical phantom experiments on the phase function parameter γ using spatial frequency domain imaging. The incorporation of two different types of scattering particles allows for control of the optical phantoms’ microscopic scattering properties. By laterally structuring areas with either TiO2 or Al2O3 scattering particles, we were able to obtain almost pure subdiffusive scattering contrast in a single optical phantom. Optical parameter mapping was then achieved using an analytical radiative transfer model revealing the microscopic structural contrast on a macroscopic field of view. As part of our study, we explain several correction and referencing techniques for high spatial frequency analysis and experimentally study the sampling depth of the subdiffusive parameter γ.

Optical phantoms with adjustable subdiffusive scattering parameters

A new epoxy-resin-based optical phantom system with adjustable subdiffusive scattering parameters is presented along with measurements of the intrinsic absorption, scattering, fluorescence, and refractive index of the matrix material. Both an aluminium oxide powder and a titanium dioxide dispersion were used as scattering agents and we present measurements of their scattering and reduced scattering coefficients. A method is theoretically described for a mixture of both scattering agents to obtain continuously adjustable anisotropy values g between 0.65 and 0.9 and values of the phase function parameter γ in the range of 1.4 to 2.2. Furthermore, we show absorption spectra for a set of pigments that can be added to achieve particular absorption characteristics. By additional analysis of the aging, a fully characterized phantom system is obtained with the novelty of g and γ parameter adjustment.

Model-based analysis on the influence of spatial frequency selection in spatial frequency domain imaging

Frequency variation in spatial frequency domain imaging is a powerful tool for adjusting the penetration depth of the imaging signal and the parameter sensitivity toward absorption and diffusive and subdiffusive scattering. Through our computational analysis, using an analytical solution of the radiative transfer equation, we add quantitation to this tool by linking the different spatial frequency regimes to their relative information content and to their absolute depth sensitivity. Special focus is placed on high spatial frequencies by analysis of the phase function parameter γ and its significance and ambiguity in describing subdiffusive scattering.

Polarization influence on reflectance measurements in the spatial frequency domain

In this work, we quantify the influence of crossed polarizers on reflectance measurements in the spatial frequency domain. The use of crossed polarizers is a very common approach for suppression of specular surface reflections. However, measurements are typically evaluated using a non-polarized scalar theory. The consequences of this discrepancy are the focus of our study, and we also quantify the related errors of the derived optical properties. We used polarized Monte Carlo simulations for forward calculation of the reflectance from different samples. The samples' scatterers are assumed to be spherical, allowing for the calculation of the scattering functions by Mie theory. From the forward calculations, the reduced scattering coefficient [Formula: see text] and the absorption coefficient μa were derived by means of a scalar theory, as commonly used. Here, we use the analytical solution of the scalar radiative transfer equation. With this evaluation approach, which does not consider polarization, we found large errors in [Formula: see text] and μa in the range of 25% and above. Furthermore, we investigated the applicability of the use of a reference measurement to reduce these errors as suggested in literature. We found that this method is not able to generally improve the accuracy of measurements in the spatial frequency domain. Our general recommendation is to apply a polarized theory when using crossed polarizers.