Determination of visible near-IR absorption coefficients of mammalian fat using time- and spatially resolved diffuse reflectance and transmission spectroscopy

In-vivo optical spectroscopy and the determination of tissue absorption and scattering properties have a central role in the development of novel optical diagnostic and therapeutic modalities in medicine. A number of techniques are available for the optical characterization of tissue in the visible near-IR region of the spectrum. An important consideration for many of these techniques is the reliability of the absorption spectrum of the various constituents of tissue. The availability of accurate absorption spectra in the range 600 to 1100 nm may allow for the determination of the concentration of key tissue constituents such as oxy- and deoxy-hemoglobin, water, and lipids. The objective of the current study is the determination of a reliable absorption spectrum of lipid(s) that can be used for component analysis of in-vivo spectra. We report the absorption spectrum of a clear purified oil obtained from pig lard. In the liquid phase above 36 degrees C, the oil is transparent and thus suitable for collimated transmission measurements. At room temperature, the oil is a solid grease that is highly scattering. The absorption and scattering properties in this solid phase are measured using time- and spatially resolved diffuse reflectance spectroscopy. Using these three independent measurement techniques, we have determined an accurate estimate for the absorption spectrum of mammalian fat.     

In vitro determination of normal and neoplastic human brain tissue optical properties using inverse adding-doubling

To complement a project towards the development of real-time optical biopsy for brain tissue discrimination and surgical resection guidance, the optical properties of various brain tissues were measured in vitro and correlated to features within clinical diffuse reflectance tissue spectra measured in vivo. Reflectance and transmission spectra of in vitro brain tissue samples were measured with a single-integrating-sphere spectrometer for wavelengths 400-1300 nm and converted to absorption and reduced scattering spectra using an inverse adding-doubling technique. Optical property spectra were classified as deriving from white matter, grey matter or glioma tissue according to histopathologic diagnosis, and mean absorption and reduced scattering spectra were calculated for the three tissue categories. Absolute reduced scattering and absorption values and their relative differences between histopathological groups agreed with previously reported results with the exception that absorption coefficients were often overestimated, most likely due to biologic variability or unaccounted light loss during reflectance/transmission measurement. Absorption spectra for the three tissue classes were dominated by haemoglobin absorption below 600 nm and water absorption above 900 nm and generally determined the shape of corresponding clinical diffuse reflectance spectra. Reduced scattering spectral shapes followed the power curve predicted by the Rayleigh limit of Mie scattering theory. While tissue absorption governed the shape of clinical diffuse reflectance spectra, reduced scattering determined their relative emission intensities between the three tissue categories.     

The determination of in vivo human tissue optical properties and absolute chromophore concentrations using spatially resolved steady-state diffuse reflectance spectroscopy

A method is described for measuring optical properties and deriving chromophore concentrations from diffuse reflection measurements at the surface of a turbid medium. The method uses a diffusion approximation model for the diffuse reflectance, in combination with models for the absorption and scattering coefficients. An optical fibre-based set-up, capable of measuring nine spectra from 400 to 1050 nm simultaneously, is used to test the method experimentally. Results of the analyses of phantom and in vivo measurements are presented. These demonstrate that in the wavelength range from 600 to 900 nm, tissue scattering can be described as a simple power dependence of the wavelength and that the tissue absorption can be accurately described by the addition of water, oxy- and deoxyhaemoglobin absorption.     

Spatial frequency domain spectroscopy of two layer media

Monitoring of tissue blood volume and oxygen saturation using biomedical optics techniques has the potential to inform the assessment of tissue health, healing, and dysfunction. These quantities are typically estimated from the contribution of oxyhemoglobin and deoxyhemoglobin to the absorption spectrum of the dermis. However, estimation of blood related absorption in superficial tissue such as the skin can be confounded by the strong absorption of melanin in the epidermis. Furthermore, epidermal thickness and pigmentation varies with anatomic location, race, gender, and degree of disease progression. This study describes a technique for decoupling the effect of melanin absorption in the epidermis from blood absorption in the dermis for a large range of skin types and thicknesses. An artificial neural network was used to map input optical properties to spatial frequency domain diffuse reflectance of two layer media. Then, iterative fitting was used to determine the optical properties from simulated spatial frequency domain diffuse reflectance. Additionally, an artificial neural network was trained to directly map spatial frequency domain reflectance to sets of optical properties of a two layer medium, thus bypassing the need for iteration. In both cases, the optical thickness of the epidermis and absorption and reduced scattering coefficients of the dermis were determined independently. The accuracy and efficiency of the iterative fitting approach was compared with the direct neural network inversion.     

Thermal coagulation-induced changes of the optical properties of normal and adenomatous human colon tissues in vitro in the spectral range 400-1,100 nm

The absorption coefficients, the reduced scattering coefficients and the optical penetration depths for native and coagulated human normal and adenomatous colon tissues in vitro were determined over the range of 400-1,100 nm using a spectrophotometer with an internal integrating sphere system, and the inverse adding-doubling method was applied to calculate the tissue optical properties from diffuse reflectance and total transmittance measurements. The experimental results showed that in the range of 400-1,100 nm there were larger absorption coefficients (P < 0.01) and smaller reduced scattering coefficients (P < 0.01) for adenomatous colon tissues than for normal colon tissues, and there were smaller optical penetration depths for adenomatous colon tissues than for normal colon tissues, especially in the near-infrared wavelength. Thermal coagulation induced significant increase of the absorption coefficients and reduced scattering coefficients for the normal and adenomatous colon tissues, and significantly reduced decrease of the optical penetration depths for the normal and adenomatous colon tissues. The smaller optical penetration depth for coagulated adenomatous colon tissues is a disadvantage for laser-induced thermotherapy (LITT) and photodynamic therapy (PDT). It is necessary to adjust the application parameters of lasers to achieve optimal therapy.     

Fast imaging of high-resolution two-dimensional effective attenuation profile from diffuse reflectance

Biological tissue is non-homogeneous in nature. It is difficult to measure its optical properties due to non-uniformity throughout the tissue being tested. To obtain the spatial distribution of optical parameters, conventional approaches use an array of light sources and detectors to reconstruct the image, thus, there is very limited spatial resolution. In contrast, solutions that provide high resolution have a high computational complexity. In this paper, we propose a fast, simple scheme to resolve the effective attenuation profile from the spatial diffuse reflectance. Rather than giving one single value for the absorption and reduced scattering coefficients, a novel algorithm is proposed for the reconstruction of an effective attenuation profile in 2-dimension from a diffuse reflectance curve. This technique is applied to the reconstruction of a 2-D effective attenuation profile. By obtaining the diffuse reflectance image from a camera and using the algorithm developed here, fast imaging of the effective attenuation profile with relatively high spatial resolution can be achieved.     

Optical characteristics of cartilage at a wavelength of 1560 nm and their dynamic behavior under laser heating conditions

A double-integrating-sphere system was used to measure the diffuse transmittance, diffuse reflectance, and collimated transmittance of cartilage and polyacrylamide hydrogel samples as a function of temperature under 1560-nm laser heating conditions. The dynamic behavior of the absorption and scattering coefficients and scattering anisotropy of the biomaterials was calculated by the inverse Monte Carlo method. The absorption coefficient of the cartilage and hydrogel samples proved to be linear in temperature. Raising the temperature of the cartilage samples to 80°C caused their absorption coefficient to decrease by some 25%. The temperature-induced change of the absorption spectrum of the interstitial water was found to be responsible for the clarification of the cartilage tissue observed to occur under 1560-nm laser heating conditions. The temperature field produced in the tissue by the laser energy deposited therein was calculated using a bioheat transfer equation with temperature-dependent parameters. The calculation results demonstrated that the temperature-induced changes of the optical parameters of biological tissues should be taken into account to make their 1560-nm laser treatment effective and safe.     

Spectroscopic method for determination of the absorption coefficient in brain tissue

I use Monte Carlo simulations and phantom measurements to characterize a probe with adjacent optical fibres for diffuse reflectance spectroscopy during stereotactic surgery in the brain. Simulations and measurements have been fitted to a modified Beer-Lambert model for light transport in order to be able to quantify chromophore content based on clinically measured spectra in brain tissue. It was found that it is important to take the impact of the light absorption into account when calculating the apparent optical path length, lp, for the photons in order to get good estimates of the absorption coefficient, μa. The optical path length was found to be well fitted to the equation lp=a+b ln(Is)+c ln(μa)+d ln(Is)ln(μa), where Is is the reflected light intensity for scattering alone (i.e., zero absorption). Although coefficients a-d calculated in this study are specific to the probe used here, the general form of the equation should be applicable to similar probes.     

Radiative transport in the delta-P1 approximation for semi-infinite turbid media

We have developed an analytic solution for spatially resolved diffuse reflectance within the deltaP1 approximation to the radiative transport equation for a semi-infinite homogeneous turbid medium. We evaluate the performance of this solution by comparing its predictions with those provided by Monte Carlo simulations and the standard diffusion approximation. We demonstrate that the delta-P1 approximation provides accurate estimates for spatially resolved diffuse reflectance in both low and high scattering media. We also develop a multi-stage nonlinear optimization algorithm in which the radiative transport estimates provided by the delta-P1 approximation are used to recover the optical absorption (microa), reduced scattering (micros'), and single-scattering asymmetry coefficients (g1) of liquid and solid phantoms from experimental measurements of spatially resolved diffuse reflectance. Specifically, the delta-P1 approximation can be used to recover microa, micros', and g1 with errors within +/- 22%, +/- 18%, and +/- 17%, respectively, for both intralipid-based and siloxane-based tissue phantoms. These phantoms span the optical property range 4 < (micros' /microa) < 117. Using these same measurements, application of the standard diffusion approximation resulted in the recovery of microa and micros' with errors o f +/- 29% and +/- 25%, respectively. Collectively, these results demonstrate that the delta-P1 approximation provides accurate radiative transport estimates that can be used to determine accurately the optical properties of biological tissues, particularly in spectral regions where tissue may display moderate/low ratios of reduced scattering to absorption (micros'/microa).     

In vivo water state measurements in breast cancer using broadband diffuse optical spectroscopy

Structural changes in water molecules are related to physiological, anatomical and pathological properties of tissues. Near infrared (NIR) optical absorption methods are sensitive to water; however, detailed characterization of water in thick tissues is difficult to achieve because subtle spectral shifts can be obscured by multiple light scattering. In the NIR, a water absorption peak is observed around 975 nm. The precise NIR peak's shape and position are highly sensitive to water molecular disposition. We introduce a bound water index (BWI) that quantifies shifts observed in tissue water absorption spectra measured by broadband diffuse optical spectroscopy (DOS). DOS quantitatively measures light absorption and scattering spectra and therefore reveals bound water spectral shifts. BWI as a water state index was validated by comparing broadband DOS to magnetic resonance spectroscopy, diffusion-weighted MRI and conductivity in bound water tissue phantoms. Non-invasive DOS measurements of malignant and normal breast tissues performed in 18 subjects showed a significantly higher fraction of free water in malignant tissues (p < 0.0001) compared to normal tissues. BWI of breast cancer tissues inversely correlated with Nottingham-Bloom-Richardson histopathology scores. These results highlight broadband DOS sensitivity to molecular disposition of water and demonstrate the potential of BWI as a non-invasive in vivo index that correlates with tissue pathology.     

Spectra from 2.5-15 microm of tissue phantom materials,optical clearing agents and ex vivo human skin: implications for depth profiling of human skin

Infrared measurements have been used to profile or image biological tissue, including human skin. Usually, analysis of such measurements has assumed that infrared absorption is due to water and collagen. Such an assumption may be reasonable for soft tissue, but introduction of exogenous agents into skin or the measurement of tissue phantoms has raised the question of their infrared absorption spectrum. We used Fourier transform infrared spectroscopy in attenuated total reflection mode to measure the infrared absorption spectra, in the range of 2-15 microm, of water, polyacrylamide, Intralipid, collagen gels, four hyperosmotic clearing agents (glycerol, 1,3-butylene glycol, trimethylolpropane, Topicare), and ex vivo human stratum corneum and dermis. The absorption spectra of the phantom materials were similar to that of water, although additional structure was noted in the range of 6-10 microm. The absorption spectra of the clearing agents were more complex, with molecular absorption bands dominating between 6 and 12 microm. Dermis was similar to water, with collagen structure evident in the 6-10 microm range. Stratum corneum had a significantly lower absorption than dermis due to a lower content of water. These results suggest that the assumption of water-dominated absorption in the 2.5-6 microm range is valid. At longer wavelengths, clearing agent absorption spectra differ significantly from the water spectrum. This spectral information can be used in pulsed photothermal radiometry or utilized in the interpretation of reconstructions in which a constant mu(ir) is used. In such cases, overestimating mu(ir) will underestimate chromophore depth and vice versa, although the effect is dependent on actual chromophore depth.     

A spectroscopic approach to imaging and quantification of cartilage lesions in human knee joints

We have previously described a technology based on diffuse reflectance of broadband light for measuring joint articular cartilage thickness, utilizing that optical absorption is different in cartilage and subchondral bone. This study is the first evaluation of the technology in human material. We also investigated the prospects of cartilage lesion imaging, with the specific aim of arthroscopic integration. Cartilage thickness was studied ex vivo in a number of sites (n = 87) on human knee joint condyles, removed from nine patients during total knee replacement surgery. A reflectance spectrum was taken at each site and the cartilage thickness was estimated using the blue, green, red and near-infrared regions of the spectrum, respectively. Estimated values were compared with reference cartilage thickness values (taken after sample slicing) using an exponential model. Two-dimensional Monte Carlo simulations were performed in a theoretical analysis of the experimental results. The reference cartilage thickness of the investigated sites was 1.60 ± 1.30 mm (mean ± SD) in the range 0-4.2 mm. Highest correlation coefficients were seen for the calculations based on the near-infrared region after normalization to the red region (r = 0.86) and for the green region (r = 0.80).     

Experimental validation of Monte Carlo modeling of fluorescence in tissues in the UV-visible spectrum

The goal of the work is to experimentally verify Monte Carlo modeling of fluorescence and diffuse reflectance measurements in turbid, tissue phantom models. In particular, two series of simulations and experiments, in which one optical parameter (absorption or scattering coefficient) is varied while the other is fixed, are carried out to assess the effect of the absorption coefficient (mu(a)) and scattering coefficient (mu(s)) on the fluorescence and diffuse reflectance measured from a turbid medium. Moreover, simulations and experiments are carried out for several fiber optic probe geometries that are designed to sample small tissue volumes. Additionally, a group of conversion expressions are derived to convert the optical properties and fluorescence quantum yield measured from tissue phantoms for use in Monte Carlo simulations. The conversions account for the differences between the definitions of the absorption coefficient and fluorescence quantum yield of fluorophores in a tissue phantom model and those in a Monte Carlo simulation. The results indicate that there is good agreement between the simulated and experimentally measured results in most cases. This dataset can serve as a systematic validation of Monte Carlo modeling of fluorescent light propagation in tissues. The simulations are carried out for a wide range of absorption and scattering coefficients as well as ratios of scattering coefficient to absorption coefficient, and thus would be applicable to tissue optical properties over a wide wavelength range (UV-visible/near infrared). The fiber optic probe geometries that are modeled in this study include those commonly used for measuring fluorescence from tissues in practice.     

Optical properties of normal and diseased human breast tissues in the visible and near infrared

The optical absorption and scattering coefficients have been determined for specimens of normal and diseased human breast tissues over the range of wavelengths from 500 to 1100 nm. Total attenuation coefficients were measured for thin slices of tissue cut on a microtome. The diffuse reflectance and transmittance were measured for 1.0 mm thick samples of these tissues, using standard integrating sphere techniques. Monte Carlo simulations were performed to derive the scattering and absorption coefficients, as well as the mean cosine of the scattering angle. The results indicate that scatter exceeds absorption by at least two orders of magnitude. Absorption is most significant at wavelengths below 600 nm. The scattering coefficients lie in the range 30-90 mm-1 at 500 nm, and fall smoothly with increasing wavelength to between 10 and 50 mm-1 at 1100 nm. The scattering coefficient for adipose tissue differs, in that it is invariant with wavelength over this spectral range. For all tissues examined, the scattered light is highly forward peaked, with the mean cosine of the scattering angle in the range 0.945-0.985. Systematic differences between the optical properties of some tissue types are demonstrated.     

Broadband absorption spectroscopy of turbid media using a dual step steady-state method

We present a method for the determination of the absorption coefficient of turbid media in a broad wavelength range with high spectral resolution using a dual step method. First, the reduced scattering coefficient is determined for a few wavelengths with spatially resolved reflectance measurements. The reduced scattering coefficient for the intermediate wavelengths is interpolated by fitting a power law. Second, the absorption coefficient is obtained from measurements of the total reflectance using the a priori knowledge of the reduced scattering coefficient. By applying a white light source and a spectrometer to measure the total reflectance, the absorption coefficient is determined with a high spectral resolution. The methodology is verified by comparing the absorption coefficients determined by the spatially resolved reflectance measurements with those obtained by the dual step method. The influence of an unknown refractive index and phase function on the determination of the optical properties is investigated. In addition, the optical properties of Intralipid/ink phantoms and the fat layer of porcine rind were determined. The absorption coefficient of the investigated phantoms varying by four orders of magnitude could be determined with an average error of less than 10%.     

Imaging of in vitro chicken leg using time-resolved near-infrared optical tomography

Near-infrared optical imaging gains much attention because of its noninvasiveness and deep penetration depths into tissue. Here, we report near-infrared optical tomographic imaging of an in vitro chicken leg from time-resolved measurements. The in vitro chicken leg, dipped in a cylindrical container filled with diluted Intralipid-10% solution, was imaged with a multichannel time-resolved imaging system. A two-dimensional reconstruction algorithm based on a modified generalized pulse spectrum technique has been developed to reconstruct the images of both the absorption and reduced scattering coefficients simultaneously and quickly. The results demonstrate that a simultaneous reconstruction of absorption and reduced scattering coefficients from time-resolved measurement has a potential to reveal the changes in the optical properties associated with not only the physiological information but also the anatomical structure of the organ.     

Visualizing depth and thickness of a local blood region in skin tissue using diffuse reflectance images

A method is proposed for visualizing the depth and thickness distribution of a local blood region in skin tissue using diffuse reflectance images at three isosbestic wavelengths of hemoglobin: 420, 585, and 800 nm. Monte Carlo simulation of light transport specifies a relation among optical densities, depth, and thickness of the region under given concentrations of melanin in epidermis and blood in dermis. Experiments with tissue-like agar gel phantoms indicate that a simple circular blood region embedded in scattering media can be visualized with errors of 6% for the depth and 22% for the thickness to the given values. In-vivo measurements on human veins demonstrate that results from the proposed method agree within errors of 30 and 19% for the depth and thickness, respectively, with values obtained from the same veins by the conventional ultrasound technique. Numerical investigation with the Monte Carlo simulation of light transport in the skin tissue is also performed to discuss effects of deviation in scattering coefficients of skin tissue and absorption coefficients of the local blood region from the typical values of the results. The depth of the local blood region is over- or underestimated as the scattering coefficients of epidermis and dermis decrease or increase, respectively, while the thickness of the region agrees well with the given values below 1.2 mm. Decreases or increases of hematocrit value give over- or underestimation of the thickness, but they have almost no influence on the depth.     

Polarized angular dependent spectroscopy of epithelial cells and epithelial cell nuclei to determine the size scale of scattering structures

An understanding of the relationship between tissue structures and light scattering from tissue will help facilitate the development and acceptance of noninvasive optical diagnostics including elastic scattering spectroscopy, diffuse reflectance, and optical coherence tomography. For example, a quantitative model of the structures that scatter light in epithelial cells would allow determination of what structures control the characteristics of in vivo light transport measurements and subsequently could provide a detailed relationship between cellular structures and optical measurements. We have determined the size distribution of refractive index structure variations in epithelial cells as well as in nuclei isolated from epithelial cells from measurements of the angular dependence of polarized light scattering. The quantitative size distributions we obtained for both whole cells and isolated nuclei include particles with effective radii of 2 microm to 10 nm or less and contain orders of magnitude more small particles than large particles. These results demonstrate that not only are biological cells very heterogeneous, but so are the nuclei within them. Light scattering is likely sensitive to structures smaller than those commonly investigated by standard pathology methods.     

Optical properties of porcine skin dermis between 900 nm and 1500 nm

The weak absorption of shortwave infrared light by skin tissues between 700 and 1500 nm offers an important window for diagnosis by optical means. The strong scattering of shortwave infrared light by the skin, however, presents a challenge to the modelling of light propagation through the skin and the understanding of skin optics. We have measured the collimated and diffuse transmittance and diffuse reflectance of porcine skin dermis samples within 30 h post-mortem. Monte Carlo simulations have been performed to inversely determine the absorption coefficient, scattering coefficient and anisotropy factor of the dermis samples in the spectral range from 900 to 1500 nm. We further analyse the sensitivity of the values of the parameters to the experimental errors and inverse calculation procedures. The state of the cellular integrity of the skin samples following optical measurements was verified using transmission electron microscopy. These results were correlated to study post-mortem effects on the in vitro optical properties of porcine dermis. We concluded that for samples stored within crushed ice for up to 30 h post-mortem the wavelength dependence of optical properties of the dermis remains unchanged while the values of the parameters vary moderately due to modification of the water content of the tissue.