Tissue phantom for optical diagnostics based on a suspension of microspheres with a fractal size distribution

Investigations on correlations between changes of optical tissue properties and NMR relaxation times

BACKGROUND

Accurate light dosimetry is a complex remaining challenge in interstitial photodynamic therapy (iPDT) for malignant gliomas. The light dosimetry should ideally be based on the tissue morphology and the individual optical tissue properties of each tissue type in the target region. First investigations are reported on using NMR information to estimate changes of individual optical tissue properties.

METHODS

Porcine brain tissue and optical tissue phantoms were investigated. To the porcine brain, supplements were added to simulate an edema or high blood content. The tissue phantoms were based on agar, Lipoveneous, ink, blood and gadobutrol (Gd-based MRI contrast agent). The concentrations of phantom ingredients and tissue additives are varied to compare concentration-dependent effects on optical and NMR properties. A 3-tesla whole-body MRI system was used to determine T1 and T2 relaxation times. Optical tissue properties, i.e., the spectrally resolved absorption and reduced scattering coefficient, were obtained using a single integrating sphere setup. The observed changes of NMR and optical properties were compared to each other.

RESULTS

By adjusting the NMR relaxation times and optical tissue properties of the tissue phantoms to literature values, recipes for human brain tumor, white matter and grey matter tissue phantoms were obtained that mimic these brain tissues simultaneously in both properties. For porcine brain tissue, it was observed that with increasing water concentration in the tissue, both NMR-relaxation times increased, while µa decreased and µs' increased at 635 nm. The addition of blood to porcine brain samples showed a constant T1, while T2 shortened and the absorption coefficient at 635 nm increased.

CONCLUSIONS

In this investigation, by changing sample contents, notable changes of both NMR relaxation times and optical tissue properties have been observed and their relations examined. The developed dual NMR/optical tissue phantoms can be used in iPDT research, clinical training and demonstrations.

Development of a platform for broadband, spectra-fitted, tissue optical phantoms

SIGNIFICANCE

Current methods of producing optical phantoms are incapable of accurately capturing the wavelength-dependent properties of tissue critical for many optical modalities.

AIM

We aim to introduce a method of producing solid, inorganic phantoms whose wavelength-dependent optical properties can be matched to those of tissue over the wavelength range of 370 to 950 nm.

APPROACH

The concentration-dependent optical properties of 20 pigments were characterized and used to determine combinations that result in optimal fits compared to the target properties over the full spectrum. Phantoms matching the optical properties of muscle and nerve, the diffuse reflectance of pale and melanistic skin, and the chromophore concentrations of a computational skin model with varying oxygen saturation ( StO 2 ) were made with this method.

RESULTS

Both optical property phantoms were found to accurately mimic their respective tissues' absorption and scattering properties across the entire spectrum. The diffuse reflectance phantoms were able to closely approximate skin reflectance regardless of skin type. All three computational skin phantoms were found to have emulated chromophore concentrations close to the model, with an average percent error for the StO 2 of 4.31%.

CONCLUSIONS

This multipigment phantom platform represents a powerful tool for creating spectrally accurate tissue phantoms, which should increase the availability of standards for many optical techniques.

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.

Spectrophotometric determination of turbid optical parameters without using an integrating sphere

Spectrophotometric quantification of turbidity by multiple optical parameters has wide-ranging applications in material analysis and life sciences. A robust system design needs to combine hardware for precise measurement of light signals with software to accurately model measurement configuration and rapidly solve a sequence of challenging inverse problems. We have developed and validated a design approach and performed system validation based on radiative transfer theory for determination of absorption coefficient, scattering coefficient, and anisotropy factor without using an integrating sphere. Accurate and rapid determination of parameters and spectra is achieved for microsphere suspension samples by combining photodiode-based measurement of four signals with the Monte Carlo simulation and perturbation-based inverse calculations. The three parameters of microsphere suspension samples have been determined from the measured signals as functions of wavelength from 400 to 800 nm and agree with calculated results based on the Mie theory. It has been shown that the inverse problems in the cases of microsphere suspension samples are well posed with convex cost functions to yield unique solutions, and it takes about 1 min to obtain the three parameters per wavelength.

Modeling the propagation of light in realistic tissue structures with MMC-fpf: a meshed Monte Carlo method with free phase function

Monte Carlo methods commonly used in tissue optics are limited to a layered tissue geometry and thus provide only a very rough approximation for many complex media such as biological structures. To overcome these limitations, a Meshed Monte Carlo method with flexible phase function choice (fpf-MC) has been developed to function in a mesh. This algorithm can model the light propagation in any complexly shaped structure, by attributing optical properties to the different mesh elements. Furthermore, this code allows the use of different discretized phase functions for each tissue type, which can be simulated from the microstructural properties of the tissue, in combination with a tool for simulating the bulk optical properties of polydisperse suspensions. As a result, the scattering properties of tissues can be estimated from information on the microstructural properties of the tissue. This is important for the estimation of the bulk optical properties that can be used for the light propagation model, since many types of tissue have never been characterized in literature. The combination of these contributions, made it possible to use the MMC-fpf for modeling the light porapagation in plant tissue. The developed Meshed Monte Carlo code with flexible phase function choice (MMC-fpf) was successfully validated in simulation through comparison with the Monte Carlo code in Multi-Layered tissues (R2 > 0.9999) and experimentally by comparing the measured and simulated reflectance (RMSE = 0.015%) and transmittance (RMSE = 0.0815%) values for tomato leaves.

Characterization of thin poly(dimethylsiloxane)-based tissue-simulating phantoms with tunable reduced scattering and absorption coefficients at visible and near-infrared wavelengths

Optical phantoms are used in the development of various imaging systems. For certain applications, the development of thin phantoms that simulate the physical size and optical properties of tissue is important. Here, we demonstrate a method for producing thin phantom layers with tunable optical properties using poly(dimethylsiloxane) (PDMS) as a substrate material. The thickness of each layer (between 115 and 880 μm) was controlled using a spin coater. The reduced scattering and absorption coefficients were controlled using titanium dioxide and alcohol-soluble nigrosin, respectively. These optical coefficients were quantified at six discrete wavelengths (591, 631, 659, 691, 731, and 851 nm) at varying concentrations of titanium dioxide and nigrosin using spatial frequency domain imaging. From the presented data, we provide lookup tables to determine the appropriate concentrations of scattering and absorbing agents to be used in the design of PDMS-based phantoms with specific optical coefficients. In addition, heterogeneous phantoms mimicking the layered features of certain tissue types may be fabricated from multiple stacked layers, each with custom optical properties. These thin, tunable PDMS optical phantoms can simulate many tissue types and have broad imaging calibration applications in endoscopy, diffuse optical spectroscopic imaging, and optical coherence tomography, etc.

Detectability of absorption and reduced scattering coefficients in frequency-domain measurements using a realistic head phantom

Detection limits of the changes in absorption and reduced scattering coefficients were investigated using a frequency-domain near-infrared system in a realistic head phantom. The results were quantified in terms of the maximum detectable depth for different activation volumes in the range of 0.8-20 microliters. The non-linear relation between the maximum detectable depth and the magnitude of changes in the absorption coefficient conform well with the Born approximation to the diffusion equation. The minimal detectable changes in the reduced scattering coefficient measured in terms of the phase signal were found to be approximately twice as large as that of the absorption coefficient using the AC signal for the same volume and at the same depth. The phase delay, which can be used to quantify the fast neuronal optical response in the human brain, showed a linear dependence on the reciprocal of the reduced scattering coefficient, as predicted by the Rytov approximation.

Biomedical tissue phantoms with controlled geometric and optical properties for Raman spectroscopy and tomography

To support the translation of Raman spectroscopy into clinical applications, synthetic models are needed to accurately test, optimize and validate prototype fiber optic instrumentation. Synthetic models (also called tissue phantoms) are widely used for developing and testing optical instrumentation for diffuse reflectance, fluorescence, and Raman spectroscopies. While existing tissue phantoms accurately model tissue optical scattering and absorption, they do not typically model the anatomic shapes and chemical composition of tissue. Because Raman spectroscopy is sensitive to molecular composition, Raman tissue phantoms should also approximate the bulk tissue composition. We describe the fabrication and characterization of tissue phantoms for Raman tomography and spectroscopy. These phantoms have controlled chemical and optical properties, and also multilayer morphologies which approximate the appropriate anatomic shapes. Tissue phantoms were fabricated to support on-going Raman studies by simulating the human wrist and rat leg. Surface meshes (triangle patch models) were generated from computed tomography (CT) images of a human arm and rat leg. Rapid prototyping was used to print mold templates with complex geometric patterns. Plastic casting techniques used for movie special effects were adapted to fabricate molds from the rapid prototypes, and finally to cast multilayer gelatin tissue phantoms. The gelatin base was enriched with additives to model the approximate chemistry and optical properties of individual tissue layers. Additional studies were performed to determine optimal casting conditions, phantom stability, layer delamination and chemical diffusion between layers. Recovery of diffuse reflectance and Raman spectra in tissue phantoms varied with probe placement. These phantoms enable optimization of probe placement for human or rat studies. These multilayer tissue phantoms with complex geometries are shown to be stable, with minimal layer delamination and chemical diffusion.

Optofluidic phantom mimicking optical properties of porcine livers

One strategy for assessing efficacy of a liver transplant is to monitor perfusion and oxygenation after transplantation. An implantable optical sensor is being developed to overcome inadequacies of current monitoring approaches. To facilitate sensor design while minimizing animal use, a polydimethylsiloxane (PDMS)-based liver phantom was developed to mimic the optical properties of porcine liver in the 630-1000 nm wavelength range and the anatomical geometry of liver parenchyma. Using soft lithography to construct microfluidic channels in pigmented elastomer enabled the 2D approximation of hexagonal liver lobules with 15mm sinusoidal channels, which will allow perfusion with blood-mimicking fluids to facilitate the development of the liver perfusion and oxygenation monitoring system.

Physical interpretation of the phase function related parameter γ studied with a fractal distribution of spherical scatterers

The optical properties within limited volumes of diffusive media can be probed by carrying spatially-resolved measurements of diffused light at short source-detector separation (typically one scattering mean free path). At such distance, analytical models only relying on the absorption and reduced scattering coefficients fail at correctly predicting reflectance and it was demonstrated that adding a third optical coefficient γ improves the description of light propagation conditions near the source. In an attempt to relate the γ coefficient to physical properties of turbid media, this paper uses a fractal distribution law for modeling scatterers' sizes distributions and investigates numerically and experimentally how γ is related to the fractal power α. The results indicate that within the range of γ typically encountered in biological samples, this coefficient is approximately linearly correlated with α.

Use of Monte Carlo simulations for propagation of light in biomedical tissues

In problems relating to light propagation in biomedical tissues, the tissue is generally modeled as a turbid medium and Monte Carlo (MC) simulation is employed to compute quantities such as diffuse reflectance, fluence, and transmittance. Two prescriptions are available in the literature for MC simulations. The first prescription considers all input quantities, including phase function, as an average over the particle size distribution, and the second prescription considers the phase function of each scatterer individually. The two prescriptions have been compared and contrasted in this paper for a given soft tissue model. It is demonstrated that, in general, the two recipes do not yield identical results. The source of this disagreement has been traced.

Monte Carlo model of spatially offset Raman spectroscopy for breast tumor margin analysis

We have previously demonstrated the discrimination of two layers of soft tissue, specifically normal breast tissue overlying breast tumor, using spatially offset Raman spectroscopy (SORS). In this report, a Monte Carlo code for evaluating SORS in soft tissues has been developed and compared to experimental results. The model was employed to investigate the effects of tissue and probe geometry on SORS measurements and therefore to develop the design strategies of applying SORS for breast tumor surgical margin evaluation. The model was used to predict SORS signals for different tissue geometries difficult to precisely control experimentally, such as varying normal and tumor layer sizes and the addition of a third layer. The results from the model suggest that, using source-detector separations of up to 3.75 mm, SORS can detect sub-millimeter-thick tumors under a 1 mm normal layer, and tumors at least 1 mm thick can be detected under a 2 mm normal layer.

Monte-Carlo-based model for the extraction of intrinsic fluorescence from turbid media

A Monte-Carlo-based model of fluorescence is developed that is capable of extracting the intrinsic fluorescence properties of tissue, which are independent of the absorption and scattering properties of tissue. This model is flexible in its applicability to different illumination-collection geometries and is also valid for a wide range of optical properties, representative of tissue in the UV-visible spectrum. This is potentially useful in a variety of biomedical applications, including cancer diagnostics and monitoring the physiological response to therapy. The model is validated using phantoms composed of hemoglobin (absorber), polystyrene spheres (scatterer), and furan-2 (fluorophore). It is found that this model is able to retrieve the intrinsic fluorescence spectra of the tissue phantoms and recover the intrinsic fluorescence intensity of furan within the phantoms to within a mean error of less than 10%.

A soft deformable tissue-equivalent phantom for diffuse optical tomography

A recipe is presented for the manufacture of highly compressible phantoms for diffuse optical tomography. The recipe is based on polyvinyl alcohol (PVA) slime, a viscoelastic fluid which readily deforms under moderate pressure. Scattering particles and absorbing compounds can be added to provide a uniform material with stable and reproducible optical properties. A linear relationship between the concentration of scattering particles (either titanium dioxide or microspheres) and the transport scatter coefficient is demonstrated. Phantoms of an arbitrary size and shape may be produced by containing the slime within a thin latex shell, and a stability over a period of at least 3 months has been established. The deformable phantoms may be used to test and calibrate optical tomography systems designed for use on patients with irregular or variable geometries.