Investigation of the source‐detector separation in near infrared spectroscopy for healthy and clinical applications

Understanding near infrared light propagation in tissue is vital for designing next generation optical brain imaging devices. Monte Carlo (MC) simulations provide a controlled mechanism to characterize and evaluate contributions of diverse near infrared spectroscopy (NIRS) sensor configurations and parameters. In this study, we developed a multilayer adult digital head model under both healthy and clinical settings and assessed light‐tissue interaction through MC simulations in terms of partial differential pathlength, mean total optical pathlength, diffuse reflectance, detector light intensity and spatial sensitivity profile of optical measurements. The model incorporated four layers: scalp, skull, cerebrospinal‐fluid and cerebral cortex with and without a customizable lesion for modeling hematoma of different sizes and depths. The effect of source‐detector separation (SDS) on optical measurements' sensitivity to brain tissue was investigated. Results from 1330 separate simulations [(4 lesion volumes × 4 lesion depths for clinical +3 healthy settings) × 7 SDS × 10 simulation = 1330)] each with 100 million photons indicated that selection of SDS is critical to acquire optimal measurements from the brain and recommended SDS to be 25 to 35 mm depending on the wavelengths to obtain optical monitoring of the adult brain function. The findings here can guide the design of future NIRS probes for functional neuroimaging and clinical diagnostic systems.      

A multipixel diffuse correlation spectroscopy system based on a single photon avalanche diode array

The autocorrelation of laser speckles from coherent near infrared light is used for noninvasive estimates of relative changes in blood perfusion in techniques such as laser Doppler flowmetry (LDF) and diffuse correlation spectroscopy (DCS). In this study, a 2D array of single photon avalanche diodes (SPADs) was used to combine the strengths of multiple detectors in LDF with high light sensitivity in DCS. The system was tested on milk phantoms with varying detector fiber diameter (200 and 600 μm), source‐detector fiber separation (4.6‐10.2 mm), fiber‐SPAD distance (2.5‐36.5 mm), contiguous measurement time per repetition for the autocorrelation (1‐33 ms) and temperature (15.6‐46.7°C). An in vivo blood occlusion test was also performed. The multipixel approach improved signal‐to‐noise ratio (SNR) and, in our setup, the use of a multimode detector fiber was beneficial for SNR. In conclusion, the multipixel system works, but improvements and further studies regarding, for example, the data acquisition and optimal settings are still needed.      

Light distribution and thermal effects in the rat brain under optogenetic stimulation

Optical brain stimulation gained a lot of attention in neuroscience due to its superior cell‐type specificity. In the design of illumination strategies, predicting the light propagation in a specific tissue is essential and requires knowledge of the optical properties of that tissue. We present the estimated absorption and reduced scattering in rodent brain tissue using non‐destructive contact spatially resolved spectroscopy (cSRS). The obtained absorption and scattering in the cortex, hippocampus and striatum are similar, but lower than in the thalamus, leading to a less deep but broader light penetration profile in the thalamus. Next, the light distribution was investigated for different stimulation protocols relevant for fiber‐optic based optogenetic experiments, using Monte Carlo simulation. A protocol specific analysis is proposed to evaluate the potential of thermally induced side effects.

     

One‐to‐one registration of en‐face optical coherence tomography attenuation coefficients with histology of a prostatectomy specimen

Optical coherence tomography (OCT), enables high‐resolution 3D imaging of the morphology of light scattering tissues. From the OCT signal, parameters can be extracted and related to tissue structures. One of the quantitative parameters is the attenuation coefficient; the rate at which the intensity of detected light decays in depth. To couple the quantitative parameters with the histology one‐to‐one registration is needed. The primary aim of this study is to validate a registration method of quantitative OCT parameters to histological tissue outcome through one‐to‐one registration of OCT with histology. We matched OCT images of unstained fixated prostate tissue slices with corresponding histology slides, wherein different histologic types were demarcated. Attenuation coefficients were determined by a supervised automated exponential fit (corrected for point spread function and sensitivity roll‐off related signal losses) over a depth of 0.32 mm starting from 0.10 mm below the automatically detected tissue edge. Finally, the attenuation coefficients corresponding to the different tissue types of the prostate were compared. From the attenuation coefficients, we produced the squared relative residue and goodness‐of‐fit metric R². This article explains the method to perform supervised automated quantitative analysis of OCT data, and the one‐to‐one registration of OCT extracted quantitative data with histopathological outcomes.      

Multiparameter wide‐field integrated optical imaging system‐based spatially modulated illumination and laser speckles in model of tissue injuries

In this report, an integrated optical platform based on spatial illumination together with laser speckle contrast technique was utilized to measure multiple parameters in live tissue including absorption, scattering, saturation, composition, metabolism, and blood flow. Measurements in three models of tissue injury including drug toxicity, artery occlusion, and acute hyperglycemia were used to test the efficacy of this system. With this hybrid apparatus, a series of structured light patterns at low and high spatial frequencies are projected onto the tissue surface and diffuse reflected light is captured by a CCD camera. A six position filter wheel, equipped with four bandpass filters centered at wavelengths of 650, 690, 800 and 880 nm is placed in front of the camera. Then, light patterns are blocked and a laser source at 650 nm illuminates the tissue while the diffusely reflected light is captured by the camera through the two remaining open holes in the wheel. In this manner, near‐infrared (NIR) and laser speckle images are captured and stored together in the computer for off‐line processing to reconstruct the tissue's properties. Spatial patterns are used to differentiate the effects of tissue scattering from those of absorption, allowing accurate quantification of tissue hemodynamics and morphology, while a coherent light source is used to study blood flow changes, a feature which cannot be measured with the NIR structured light. This combined configuration utilizes the strengths of each system in a complementary way, thus collecting a larger range of sample properties. In addition, once the flow and hemodynamics are measured, tissue oxygen metabolism can be calculated, a property which cannot be measured independently. Therefore, this merged platform can be considered a multiparameter wide‐field imaging and spectroscopy modality. Overall, experiments demonstrate the capability of this spatially coregistered imaging setup to provide complementary, useful information of various tissue metrics in a simple and noncontact manner, making it attractive for use in a variety of biomedical applications.     

Monte Carlo modeling of photon migration in realistic human thoracic tissues for noninvasive monitoring of cardiac hemodynamics

Noninvasive monitoring of cardiac hemodynamics remains challenging in cardiovascular medicine. The possibility of noninvasive optical monitoring of cardiac hemodynamics was theoretically investigated in this study. By utilizing the Monte Carlo simulation method for voxelized media (MCVM) and Visible Chinese Human dataset, we quantified and visualized the photon migration in human thoracic region. The light fluence distribution was showed to reach heart tissue (～3 cm depth underbody surface) and 12% of the total fluence was absorbed by the myocardium. The proportion of spatial sensitivity distribution (SSD) in cardiac tissue to the total SSD reached 0.0195%. The portion of SSD increased following with cardiac diastole and diffuse reflectance deceased linearly with increasing cardiac volume. The optimal separation between the light source and detector was provided to be 3.5 to 4.0 cm for future development of noninvasive cardiac hemodynamics monitoring. A pilot experimental study was conducted to measure the diffuse reflectance light and fingertip photoplethysmography. These data suggest that the fluctuation period of near‐infrared (NIR) diffuse reflectance was consistent with the cardiac cycle, while the fluctuation features of the NIR signal was not consistent with that of photoplethysmography. All results indicate the great potential of noninvasive optical monitoring of myocardial hemodynamics.     

Depth‐resolved multimodal imaging: Wavelength modulated spatially offset Raman spectroscopy with optical coherence tomography

A major challenge in biophotonics is multimodal imaging to obtain both morphological and molecular information at depth. We demonstrate a hybrid approach integrating optical coherence tomography (OCT) with wavelength modulated spatially offset Raman spectroscopy (WM‐SORS). With depth colocalization obtained from the OCT, we can penetrate 1.2‐mm deep into strong scattering media (lard) to acquire up to a 14‐fold enhancement of a Raman signal from a hidden target (polystyrene) with a spatial offset. Our approach is capable of detecting both Raman and OCT signals for pharmaceutical particles embedded in turbid media and revealing the white matter at depth within a 0.6‐mm thick brain tissue layer. This depth resolved label‐free multimodal approach is a powerful route to analyze complex biomedical samples.      

Functional near infrared spectroscopy using spatially resolved data to account for tissue scattering: A numerical study and arm‐cuff experiment

Functional Near‐Infrared Spectroscopy (fNIRS) aims to recover changes in tissue optical parameters relating to tissue hemodynamics, to infer functional information in biological tissue. A widely‐used application of fNIRS relies on continuous wave (CW) methodology that utilizes multiple distance measurements on human head for study of brain health. The typical method used is spatially resolved spectroscopy (SRS), which is shown to recover tissue oxygenation index (TOI) based on gradient of light intensity measured between two detectors. However, this methodology does not account for tissue scattering which is often assumed. A new parameter recovery algorithm is developed, which directly recovers both the scattering parameter and scaled chromophore concentrations and hence TOI from the measured gradient of light‐attenuation at multiple wavelengths. It is shown through simulations that in comparison to conventional SRS which estimates cerebral TOI values with an error of ±12.3%, the proposed method provides more accurate estimate of TOI exhibiting an error of ±5.7% without any prior assumptions of tissue scatter, and can be easily implemented within CW fNIRS systems. Using an arm‐cuff experiment, the obtained TOI using the proposed method is shown to provide a higher and more realistic value as compared to utilizing any prior assumptions of tissue scatter.     

Polarization‐multiplexed,dual‐beam swept source optical coherence tomography angiography

A polarization‐multiplexed, dual‐beam setup is proposed to expand the field of view (FOV) for a swept source optical coherence tomography angiography (OCTA) system. This method used a Wollaston prism to split sample path light into 2 orthogonal‐polarized beams. This allowed 2 beams to shine on the cornea at an angle separation of ~14°, which led to a separation of ~4.2 mm on the retina. A 3‐mm glass plate was inserted into one of the beam paths to set a constant path length difference between the 2 polarized beams so the interferogram from the 2 beams are coded at different frequency bands. The resulting OCTA images from the 2 beams were coded with a depth separation of ~2 mm. A total of 5 × 5 mm² angiograms from the 2 beams were obtained simultaneously in 4 seconds. The 2 angiograms then were montaged to get a wider FOV of ~5 × 9.2 mm².      

Imaging depth variations in hyperspectral imaging: Development of a method to detect tumor up to the required tumor‐free margin width

Hyperspectral imaging is a promising technique for resection margin assessment during cancer surgery. Thereby, only a specific amount of the tissue below the resection surface, the clinically defined margin width, should be assessed. Since the imaging depth of hyperspectral imaging varies with wavelength and tissue composition, this can have consequences for the clinical use of hyperspectral imaging as margin assessment technique. In this study, a method was developed that allows for hyperspectral analysis of resection margins in breast cancer. This method uses the spectral slope of the diffuse reflectance spectrum at wavelength regions where the imaging depth in tumor and healthy tissue is equal. Thereby, tumor can be discriminated from healthy breast tissue while imaging up to a similar depth as the required tumor‐free margin width of 2 mm. Applying this method to hyperspectral images acquired during surgery would allow for robust margin assessment of resected specimens. In this paper, we focused on breast cancer, but the same approach can be applied to develop a method for other types of cancer.     

生物组织漫射光模拟计算的改进杂合模型

结合近源处Ｍｏｎｔｅ Ｃａｒｌｏ法的精确性和漫射理论计算的快速性,在用漫射理论计算漫反射系数时,采用了辐射对背半球的积分形式,即考虑光强度和通量两项对它的贡献,建立了基于Ｍｏｎｔｅ Ｃａｒｌｏ法－漫理论的更为精确、有效的改进的杂合模型,模拟计算了板状生物组织漫反射光的分布。结果表明,在保证高计算速度的前提下,进一步提高了计算精度;同时,此模型也适用于强吸收介质漫反射光分布的模拟计算。     

Quantitative subsurface spatial frequency‐domain fluorescence imaging for enhanced glioma resection

The rate of complete resection of glioma has improved with the introduction of 5‐aminolevulinic acid‐induced protoporphyrin IX (PpIX) fluorescence image guidance. Surgical outcomes are further enhanced when the fluorescence signal is decoupled from the intrinsic tissue optical absorption and scattering obtained from diffuse reflectance measurements, yielding the absolute PpIX concentration, [PpIX]. Spatial frequency domain imaging was used previously to measure [PpIX] in near‐surface tumors under blue fluorescence excitation. Here, we extend this to subsurface [PpIX] fluorescence under red‐light excitation. The decay rate of the modulation amplitude of the fluorescence signal was used to calculate the PpIX depth, which was then applied in a forward diffusion model to estimate [PpIX] at depth. For brain‐like optical properties in phantoms with PpIX fluorescent inclusions, the depth can be recovered up to depths of 9.5 mm ± 0.4 mm, with [PpIX] ranging from 5 to 15 μg/mL within an average deviation of 15% from the true [PpIX] value.      

Sex,but not skin tone affects penetration of red‐light (660 nm) through sites susceptible to sports injury in lean live and cadaveric tissues

Red‐light treatment is emerging as a novel therapy for promoting tissue recovery but data on red‐light penetration through human tissues are lacking. We aimed to: (1) determine the effect of light irradiance, tissue thickness, skin tone, sex and bone/muscle content on 660 nm light penetration through common sites of sports injuries, and (2) establish if cadaver tissues serve as a useful model for predicting red‐light penetration in live tissues. Live and cadaver human tissues were exposed to 660 nm light at locations across the skull, spinal cord and upper and lower limbs. Red‐light was produced by a light emitting diode array of various irradiances (15‐500 mW/cm²) and measured by a light‐probe positioned on the tissue surface opposite to the light emitting diodes. 100 mW/cm² successfully penetrated tissue <50 mm thick; a disproportionate irradiance increase was required to achieve deeper penetration. Penetration was unaffected by skin tone, increased with irradiance and relative bone/muscle composition, and decreased with greater tissue thickness and in males. Live and cadaveric tissue penetration did not differ statistically for tissues <50 mm but cadavers required more red‐light to penetrate >50 mm. These results assist clinicians and researchers in determining red‐light treatment intensities for penetrating human tissues.      

The influence of the blood vessel diameter on the full scattering profile from cylindrical tissues: experimental evidence for the shielding effect

Optical methods for detecting physiological state based on light–tissue interaction are noninvasive, inexpensive, simplistic, and thus very useful. The blood vessels in human tissue are the main cause of light absorbing and scattering. Therefore, the effect of blood vessels on light–tissue interactions is essential for optically detecting physiological tissue state, such as oxygen saturation, blood perfusion and blood pressure. We have previously suggested a new theoretical and experimental method for measuring the full scattering profile, which is the angular distribution of light intensity, of cylindrical tissues. In this work we will present experimental measurements of the full scattering profile of heterogenic cylindrical phantoms that include blood vessels. We show, for the first time that the vessel diameter influences the full scattering profile, and found higher reflection intensity for larger vessel diameters accordance to the shielding effect. For an increase of 60% in the vessel diameter the light intensity in the full scattering profile above 90° is between 9% to 40% higher, depending on the angle. By these results we claim that during respiration, when the blood‐vessel diameter changes, it is essential to consider the blood‐vessel diameter distribution in order to determine the optical path in tissues.

A CT scan of the measured silicon‐based phantoms. The phantoms contain the same blood volume in different blood‐vessel diameters.     

Cone–shell Raman spectroscopy (CSRS) for depth‐sensitive measurements in layered tissue

We report the development of a depth‐sensitive Raman spectroscopy system using the configuration of cone–shell excitation and cone detection. The system uses a 785 nm diode laser and three identical axicons for Raman excitation of the target sample in the form of a hollow conic section. The Raman scattered light from the sample, passed through the same (but solid) conic section, is collected for detection. Apart from its ability of probing larger depths (? few mm), an important attraction of the system is that the probing depths can be varied by simply varying the separation between axicons in the excitation arm. Furthermore, no adjustment is required in the sample arm, which is a significant advantage for noncontact, depth‐sensitive measurement. Evaluation of the performance of the developed setup on nonbiological phantom and biological tissue sample demonstrated its ability to recover Raman spectra of layers located at depths of ?2–3 mm beneath the surface.

     

Inside Cover: Real‐time optical spectroscopic monitoring of nonirrigated lesion progression within atrial and ventricular tissues (J. Biophotonics 4/2019)

There is a clinical need in the treatment of atrial fibrillation by radiofrequency ablation to provide lesion set validation. The integration of optical fibers into radiofrequency ablation catheters can address this issue. Thermally‐induced, dynamic changes in tissue diffuse reflectance are measured through transcatheter near‐infrared spectroscopy. A regression model based on spectral features is proposed and utilized to validate the presence of thermal injury, measure lesion depth, and assess gaps. Further details can be found in the article by Rajinder P. Singh‐Moon, Xinwen Yao, Vivek Iyer, et al. ( e201800144 ).

     

Variability in odor-modulated flight by moths

Based on previous studies of odor-modulated flight where track parameter data was lumped and averaged, the speed and orientation of the moths' movement along their flight tracks have been said to be controlled to maintain certain “preferred” values. The results from our fine-scaled analysis of this behavior show that none of the track parameters typically measured are held constant. The moths' speed along the flight track is modulated substantially and predictably: fastest along the straight legs and slowest around the turns. In addition, about half of the individuals studied progressively reduced the peak speed along the straight legs as they approached the pheromone source. While most of the track legs between the turns were directed upwind, their orientations were widely distributed, indicating no preferred direction. Small fluctuations of orientation along some straight legs suggest corrective maneuvers to stabilize flight direction about an internal set point. The visual inputs hypothesized to control steering and speed, transverse and longitudinal image flow, changed continuously during upwind flight in pheromone, but no regular relationship between them was observed. We found that the orientation of the longitudinal body axis and the direction of thrust (course angle) were only rarely coincident during upwind flight to the odor source, suggesting that moths receive sensory input which differs quantitatively from that calculated by conventional methods. Our results strongly suggest that the long-accepted hypothetical mechanisms of control for this behavior do not operate in the manner in which they have been proposed. Accepted: 11 July 1997     

Interplay between Bimolecular Recombination and Carrier Transport Distances in Bulk Heterojunction Organic Solar Cells

In this work, it is demonstrated that bimolecular recombination depends on the distance that free carriers are required to travel in transit to the electrodes in bulk heterojunction organic solar cells. By employing semi‐transparent devices, the carrier transport distance can be controlled via the local light absorption profile with an appropriate choice of the illumination side and incident wavelength. Using a series of light intensity‐dependent measurements, bimolecular recombination is shown to depend on the distance electrons or holes are required to transit the active layer. This effect is demonstrated for three different bulk heterojunction blends, where the restrictive carrier that causes the onset of recombination is identified. The mobility‐lifetime products of the limiting carriers are also estimated using a simple model for carrier extraction, where similar values are obtained regardless of the absorption profile. Implications for 1‐sun performance are also discussed, which provide guidelines for fabricating devices with thicker active layers capable of maximizing light absorption without succumbing to recombination losses.     

Measurement of light gradients and spectral regime in plant tissue with a fiber optic probe

A method is described in which light gradients and spectral regime can be measured within plant tissue using fiber optics. A fiber optic probe was made by modifying a single optical fiber (200 μm diameter) so that it had a light harvesting end that was a truncated tip 20–70 μm in diameter. The probe was a directional sensor with a half-band acceptance angle of 17–20°. Light measurements were made as the fiber optic probe was driven through plant tissue by a motorized micromanipulator, and the light that entered the fiber tip was piped to a spectroradiometer. By irradiating green leaf tissue of the succulent Crassula falcata L. with collimated light and inserting the probe from different directions, it was possible to measure light quality and quantity at different depths. Collimated light was scattered completely by the initial 1.0 mm of leaf tissue, which also greatly attenuated all light except the green and far-red. Light scatter contributed significantly to light quantity and had a pronounced spectral structure. Immediately beneath the irradiated surface the amount of light at 550 nm was 1.2 times that of the incident light. The light gradient declined rapidly to 0.5 times incident light at 1.4 mm depth. In contrast, the amount of light at 750 nm increased during the initial 0.5 mm to 2.9 times incident light and then declined linearly to 0.5 times incident light at the dark side of the leaf (4.5 mm). The implications of the magnitude of the contribution of light scatter to the light gradient is also discussed.     

Interaction Mechanisms of Cavitation Bubbles Induced by Spatially and Temporally Separated fs-Laser Pulses

The emerging use of femtosecond lasers with high repetition rates in the MHz regime together with limited scan speed implies possible mutual optical and dynamical interaction effects of the individual cutting spots. In order to get more insight into the dynamics a time-resolved photographic analysis of the interaction of cavitation bubbles is presented. Particularly, we investigated the influence of fs-laser pulses and their resulting bubble dynamics with various spatial as well as temporal separations. Different time courses of characteristic interaction effects between the cavitation bubbles were observed depending on pulse energy and spatio-temporal pulse separation. These ranged from merely no interaction to the phenomena of strong water jet formation. Afterwards, the mechanisms are discussed regarding their impact on the medical application of effective tissue cutting lateral to the laser beam direction with best possible axial precision: the mechanical forces of photodisruption as well as the occurring water jet should have low axial extend and a preferably lateral priority. Furthermore, the overall efficiency of energy conversion into controlled mechanical impact should be maximized compared to the transmitted pulse energy and unwanted long range mechanical side effects, e.g. shock waves, axial jet components. In conclusion, these experimental results are of great importance for the prospective optimization of the ophthalmic surgical process with high-repetition rate fs-lasers.